Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compressible flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form were iterated in time until a steady solution was reached. Evidence was given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at the domain boundaries was proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects were accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. The flow in a developing S-duct was then solved in the laminar regime in a Reynolds number (Re) of 790 and in the turbulent regime at Re equals 40,000 by using the Baldwin-Lomax model. The Stanitz elbow was then solved by using an invicid version of the same code at M sub inlet equals 0.4. Grid dependence and convergence rate were investigated, showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re equals 2.5 times 10(exp 6) was solved with the Baldwin-Lomax and the q-omega models. Both approaches show satisfactory agreement with experiments, although the q-omega model was slightly more accurate.

Michelassi, V.↗

Control of unsteady separated flow associated with the dynamic stall of airfoils

The two principal objectives of this research were to achieve an improved understanding of the mechanisms involved in the onset and development of dynamic stall under compressible flow conditions, and to investigate the feasibility of employing adaptive airfoil geometry as an active flow control device in the dynamic stall engine. Presented here are the results of a quantitative (PDI) study of the compressibility effects on dynamic stall over the transiently pitching airfoil, as well as a discussion of a preliminary technique developed to measure the deformation produced by the adaptive geometry control device, and bench test results obtained using an airfoil equipped with the device.

Wilder, Michael C.↗

Stability analysis for laminar flow control, part 2

Topics covered include: (1) optimization of the numerics of the SALLY stability analysis code; (2) relation between temporal and spatial stability theory; (3) compressible flow stability calculations; (4) spectral methods for the boundary layer equations; and (5) numerical study of nonlinear, nonparallel stability of incompressible flows.

Orszag, S. A.↗

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗

Application of a generalized minimal residual method to 2D unsteady flows

A generalized minimum residual scheme (GMRES), previously developed for solving nonlinear and linear systems of equations, has been applied to the numerical solution of 2D unsteady compressible flows. It is found that the use of GMRES significantly increases the time step that may be used, compared to noniterative implicit schemes. The feasibility of reducing the memory requirements of the GMRES scheme using a multigrid strategy has also been explored. Several sample steady and unsteady viscous flow applications are presented.

Hixon, Ray↗

A New Finite Element Approach for Prediction of Aerothermal Loads: Progress in Inviscid Flow Computations

Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Two dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Research progress is described for: (1) utilization of a commerically available program to construct flow solution domains and display computational results, (2) development of an explicit Taylor-Galerkin solution algorithm, (3) closed form evaluation of finite element matrices, (4) vector computer programming strategies, and (5) validation of solutions. Two test problems of interest to NASA Langley aerothermal research are studied. Comparisons of finite element solutions for Mach 6 flow with other solution methods and experimental data validate fundamental capabilities of the approach for analyzing high speed inviscid compressible flows.

Bey, K. S.↗

A new finite element approach for prediction of aerothermal loads - Progress in inviscid flow computations

Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Two dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Research progress is described for: (1) utilization of a commercially available program to construct flow solution domains and display computational results, (2) development of an explicit Taylor-Galerkin solution algorithm, (3) closed form evaluation of finite element matrices, (4) vector computer programming strategies, and (5) validation of solutions. Two test problems of interest to NASA Langley aerothermal research are studied. Comparisons of finite element solutions for Mach 6 flow with other solution methods and experimental data validate fundamental capabilities of the approach for analyzing high speed inviscid compressible flows.

Bey, K. S.↗

Solutions to Kuessner's integral equation in unsteady flow using local basis functions

The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.

Fromme, J. A.↗

Unified approach for incompressible flows

An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Chang, Tyne-Hsien↗

Development and calibration of buried wire gages for wall shear stress measurements in fluid flow

Special methods were developed to arrange 'Buried Wire Gage' inserts flush to the contoured flow surfaces of instrument plugs of a boundary-layer flow apparatus. The fabrication process was aimed at producing proper bonding of the sensor wire to the substrate surface, without causing excessive surface waviness. A large number of gages were built and first calibrated for the resistance-temperature characteristics. The gages were then installed in a flow calibration apparatus and operated from a constant temperature anemometer system for a series of flow settings to derive the calibration constants of each of the gages. The flow settings included a range of subsonic freestream Mach numbers in order to help establish the gage calibration characteristics for compressible flow fields. This paper provides a description of the buried wire gage technique, an explanation of the method evolved for making proper gages, the procedure for calibrating the gages and the results of measurements performed for determining the calibration constants.

Murthy, Sreedhara V.↗

CMOTT Turbulence Module for NPARC

This is a user's manual of the CMOTT turbulence module, version 2.0, developed for the NPARC code. The module is written in a self-contained manner so that the user can use any turbulence model in the module without concern as to how it is implemented and solved. Three two-equation turbulence models have been built into the module: Chien, Shih-Lumley and CMOTT models, and all of them have both the low Reynolds number and wall function options. Unlike Chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +) in the low Reynolds number approach, an advantage for separated flow calculations. The Van Driest transformation is used so that the wall functions can be applied to both incompressible and compressible flows. The manual gives the details of the turbulence models used and their numerical implementation. It also gives two application examples, one for subsonic and the other for transonic flow, for demonstration. The module can be easily linked to the NPARC code for practical applications.

Zhu, J.↗

Application of a Real-Fluid Turbomachinery Analysis to Rocket Turbopump Geometries

A three-dimensional flow solver has been developed for turbomachinery components utilizing real fluid properties. The code is applicable to both incompressible and compressible flow fields. In this study, the code has been applied to the analysis of inducer and impeller geometries representative of those used in rocket engine applications. The predicted results show good agreement with the available experimental data.

Dorney, Daniel J.↗

Simulation Of Flow In A Turbine Cascade

Report discusses application of lower-upper (LU) implicit numerical-integration scheme, together with interactive grid-generation method, to simulation of inviscid, compressible flow in two-dimensional model of turbine cascade. Study part of continuing effort to develop techniques for computation of complicated three-dimensional flows in turbomachinery.

Choo, Yung K.↗

Shock Tube Radiation Measurement in Expanding Air Flow

The EAST facility has been widely used in the past to provide experimental data aiming to validate radiative heating models. However, all the past studies involved measurements of compressing flows. In order to study recombining flows, an expansion cone has been installed in the EAST facility. Four spectrometers are mounted at different axial positions of the cone to measure the evolution of the radiation during expansion. Quantitative data from optical emission spectroscopy will be used to get the densities of radiative species. These data are intended to be compared with simulations in order to improve the current models used by the entry community.

Augustin Tibere-Inglesse↗

Shock Tube Radiation Measurement in Expanding Air Flows

The EAST facility has been widely used in the past to provide experimental data aiming to validate radiative heating models. However, all the past studies involved measurements of compressing flows. In order to study recombining flows, an expansion cone has been installed in the EAST facility. Four spectrometers are mounted at different axial positions of the cone to measure the evolution of the radiation during expansion. Quantitative data from optical emission spectroscopy will be used to get the densities of radiative species. These data are intended to be compared with simulations in order to improve the current models used by the entry community.

Augustin Tibère-Inglesse↗

Simultaneous measurements of velocity and pressure fields in subsonic and supersonic flows through image-intensified detection of laser-induced fluorescence

An optical technique is presented for combined, spatially resolved measurements of two-dimensional velocity and pressure fields in compressible flows. The single-mode frequency of an argon laser is fixed in the wing of an absorption line of iodine molecules, seeded in an underexpanded round jet of nitrogen gas. The emitted fluorescence, being proportional to the amount of absorbed radiation and hence the absorption line-shape function, is detected with an intensified 100 x 100 photodiode array camera. A single-microchannel-plate image intensifier is fiber-optically coupled to the array in order to improve time resolution and SNR. Three components of the velocity vector in a cross-sectional plane are sequentially probed with four laser sheets from three different directions. By shifting the laser frequency in one of the sheets with a piezo-tuned intra-cavity etalon, the slope of the absorption line can be measured in situ in order to provide the required scaling factor for the velocity measurement. With its short measurement times of less than 250 ms, this method is well suited for blow-down wind tunnel experiments.

Hiller, B.↗