SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A multiple-scale eddy viscosity model is described. This model splits the energy spectrum into a high wave number regime and a low wave number regime. Dividing the energy spectrum into multiple regimes simplistically emulates the cascade of energy through the turbulence spectrum. The constraints on the model coefficients are determined by examining decaying turbulence and homogeneous turbulence. A direct link between the partitioned energies and the energy transfer process is established through the coefficients. This new model was calibrated and tested for boundary-free turbulent shear flows. Calculations of mean and turbulent properties show good agreement with experimental data for two mixing layers, a plane jet and a round jet.
This paper reports the progress being made towards complete unsteady turbopump simulation capability by using overset grid systems. A computational model of a turbo-pump impeller is used as a test case for the performance evaluation of the MPI, hybrid MPI/Open-MP, and MLP versions of the INS3D code. Relative motion of the grid system for rotor-stator interaction was obtained by employing overset grid techniques. Unsteady computations for a turbo-pump, which contains 114 zones with 34.3 Million grid points, are performed on Origin 2000 systems at NASA Ames Research Center. The approach taken for these simulations, and the performance of the parallel versions of the code are presented.
This viewgraph presentation provides information on mathematical modelling of the SSME (space shuttle main engine). The unsteady SSME-rig1 start-up procedure from the pump at rest has been initiated by using 34.3 million grid points. The computational model for the SSME-rig1 has been completed. Moving boundary capability is obtained by using DCF module in OVERFLOW-D. MPI (Message Passing Interface)/OpenMP hybrid parallel code has been benchmarked.
Performance of the two commonly used numerical procedures, one based on artificial compressibility method and the other pressure projection method, are compared. These formulations are selected primarily because they are designed for three-dimensional applications. The computational procedures are compared by obtaining steady state solutions of a wake vortex and unsteady solutions of a curved duct flow. For steady computations, artificial compressibility was very efficient in terms of computing time and robustness. For an unsteady flow which requires small physical time step, pressure projection method was found to be computationally more efficient than an artificial compressibility method. This comparison is intended to give some basis for selecting a method or a flow solution code for large three-dimensional applications where computing resources become a critical issue.
During the grant period from January 1,2002 to December 31,2002 work was carried out on three projects to extend the range of applicability of advanced turbulence models. First, a new transition-sensitized turbulence model was tested and refined. Second, the influence of compressibility on the pressure-strain rate correlation was studied. Third, the relationship between time-filtered large eddy simulation (TLES) and Reynolds-averaged Navier Stokes (RANS) modeling was investigated leading to submission of the article. The transition-sensitized turbulence model encompasses the early-stage transition and turbulent flow regimes describing the evolution of the ensemble mean disturbance energy and dissipation rate. It is founded on a consistent mathematical description of the laminar regime with its linear disturbances and the fully turbulent regime with its stochastic fluctuations. The unified description is provided by the ensemble viewpoint.
The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.
The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.
A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.
Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.
A conference was held on Computational Fluid Dynamics (CFD) and produced related papers. Topics included CFD algorithms, transition and turbulent flow, hypersonic reacting flow, incompressible flow, two phase flow and combustion, internal flow, compressible flow, grid generation and adaption, boundary layers, environmental and industrial applications, and non-Newtonian flow.
Theory suggests that three distinct types of turbulence can occur in the low Mach number limit of polytropic flow: nearly incompressible flows dominated by vorticity, nearly pure acoustic turbulence dominated by compression, and flows characterized by near statistical equipartition of vorticity and compressions. Distinctions between these kinds of turbulence are investigated here by direct numerical simulation of two-dimensional compressible hydrodynamic turbulence. Dynamical scalings of density fluctuations, examination of the ratio of transverse to longitudinal velocity fluctuations, and spectral decomposition of the fluctuations are employed to distinguish the nature of these low Mach number solutions. A strong dependence on the initial data is observed, as well as a tendency for enhanced effects of compressibility at later times and at higher wave numbers, as suggested by theories of nearly incompressible flows.
For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.
Turbomachines for propulsion applications operate with many different working fluids and flow conditions. The flow may be incompressible, such as in the liquid hydrogen pump in a rocket engine, or supersonic, such as in the turbine which may drive the hydrogen pump. Separate codes have traditionally been used for incompressible and compressible flow solvers. The General Equation Set (GES) method can be used to solve both incompressible and compressible flows, and it is not restricted to perfect gases, as are many compressible-flow turbomachinery solvers. An unsteady GES turbomachinery flow solver has been developed and applied to both air and water flows through turbines. It has been shown to be an excellent alternative to maintaining two separate codes.
The full Navier-Stokes equations for incompressible turbulent flow must be solved to accurately represent all flow phenomena which occur in a high Reynolds number incompressible flow. A two layer algebraic eddy viscosity turbulence model is used to represent the Reynolds stress in the primitive variable formulation. The development of the boundary-fitted coordinate systems makes the numerical solution of these equations feasible for arbitrarily shaped bodies. The nondimensional time averaged Navier-Stokes equations, including the turbulence mode, are represented by finite difference approximations in the transformed plane. The resulting coupled system of nonlinear algebraic equations is solved using a point successive over relaxation iteration. The test case considered was a NACA 64A010 airfoil section at an angle of attack of two degrees and a Reynolds number of 2,000,000.
Incompressible laminar boundary layer development over two-dimensional symmetric body in presence of shear flow