Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUILIBRIUM FLOW”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Predicted aerodynamics for a proposed personnel launch vehicle

The Langley Aerothermodynamic Upwind Relaxation Algorithm has been modified to compute equilibrium flows over shapes at hypersonic velocities. As an initial test of the modified code, solutions were obtained for a paraboloid at Mach 6 wind-tunnel test conditions in CF4 gas and at Mach 16.85 flight conditions in equilibrium air that approximated the normal shock density ratio for the experimental test. Solutions for a lifting-body concept that is proposed for a Personnel Launch System were compared with test results at Mach 6 in CF4 over an angle-of-attack range of 15 to 35 deg and with results for nominal entry flight conditions at Mach 15.

Greene, Francis A.↗
NASA NTRS1991

Nonadiabatic and three-dimensional effects in compressible turbulent boundary layers

A defect stream function formulation for nonadiabatic flow with small crossflow is developed. The first-integral property of this formulation provides for two removal of the streamline curvature term in the governing equation so that the form of the reduced equation for small crossflow is the same as that for two dimensional flow. The combined law of the wall and wake is used in place of the no-slip boundary condition. The tangential velocity equation for law-of-the-wall flow is shown to be the same for three-dimensions as for two when the Boussinesq approximation applies, and a closed form solution for the crossflow angle in the inner region is obtained. Analytic solutions for nonadiabatic, compressible, equilibrium flow with a Clauser outer-region eddy-viscosity model are obtained, and excellent agreement with experimental skin friction and velocity profile data for nonadiabatic, compressible flat-plate flow is achieved. An analytic solution for a linear inner-region eddy-viscosity model is also obtained; the wake function part of this solution is found to be inconsistent with the empirically established law of the wake.

Barnwell, Richard W.↗
NASA NTRS1993

The addition of algebraic turbulence modeling to program LAURA

The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) is modified to allow the calculation of turbulent flows. This is accomplished using the Cebeci-Smith and Baldwin-Lomax eddy-viscosity models in conjunction with the thin-layer Navier-Stokes options of the program. Turbulent calculations can be performed for both perfect-gas and equilibrium flows. However, a requirement of the models is that the flow be attached. It is seen that for slender bodies, adequate resolution of the boundary-layer gradients may require more cells in the normal direction than a laminar solution, even when grid stretching is employed. Results for axisymmetric and three-dimensional flows are presented. Comparison with experimental data and other numerical results reveal generally good agreement, except in the regions of detached flow.

Cheatwood, F. Mcneil↗
NASA NTRS1979

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗
NASA NTRS1987

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗
NASA NTRSDate not supplied

Shock-Layer Radiation Insights Available Through Flowfield Property Binning

This paper introduces a method for assessing similarity between available shock-layer radiation measurements and the radiative heating to a vehicle of interest. This method will be referred to as the flowfield property binning (FPB) approach. The challenge in developing this approach is that radiometer measurements from a surface normal line of sight, observed radiation measurements of the entire shock-layer from a distance, or measured shock-tube nonequilibrium profiles are notably different geometrically than the radiative flux to a vehicle surface. Therefore, the relative contributions of various regions of the shock-layer to the measurement, and how these relative contributions compare with the radiative heating to a vehicle of interest, are unclear. By tracking the radiation contribution from each flowfield location, and then grouping this contribution based on two local flowfield properties, these geometric differences are removed. This provide a two dimensional group, or bin, distribution that represents the radiation environment in the same format for various types of radiation simulations, such as observed radiation or radiative flux to a vehicle surface. In the limits of optically-thin radiation transport and thermochemical equilibrium flow, the choice of pressure and temperature as the binning flowfield properties is shown to provide a bin distribution that uniquely defines the radiation simulation. This concept is found to be sufficient for establishing similarity between two cases even when the optically-thin radiation transport and thermochemical equilibrium assumptions are not satisfied. A simple metric of comparison between bin distributions from different cases is provided by the histogram intersection algorithm. This algorithm results in a histogram intersection factor (H) between 0 and 1, which represents the fraction of the radiation that is similar between the two cases. Therefore, if H evaluated between bin distributions for two completely different radiation problems is near 1, then they represent essentially the same radiation problem. This means that they depend on the same weighting of various radiation modeling parameters evaluated at the same pressure and temperature. If one of these cases agrees well with a radiation measurement, then the matching bin distributions provide justification for applying this experimental validation to the other case. Even for H values as low as 0.2, the measurement case could be considered relevant if there are no other useful measurements. This approach allows observed radiation measurements to directly inform radiative heating simulations, which allows NASA’s past and future investments in observed radiation measurements to directly inform the radiative heating uncertainty or margins for future vehicle design. This FPB approach is also used to determine shock-tube conditions for informing backshell radiative heating simulations, where the normal shock velocity and Lagrangian time from the shock is used for the binning variables instead of pressure and temperature.

Shock Layer Radiation↗
NASA NTRS1993

DSMC simulation of the interaction between rarefied free jets

This paper presents a direct simulation Monte Carlo (DSMC) calculation of two interacting free jets exhausting into vacuum. The computed flow field is compared against available experimental data and shows excellent agreement everywhere except in the very near field (less than one orifice diameter downstream of the jet exhaust plane). The lack of agreement in this region is attributed to having assumed an inviscid boundary condition for the orifice lip. The results serve both to validate the DSMC code for a very complex, three dimensional non-equilibrium flow field, and to provide some insight as to the complicated nature of this flow.

Dagum, Leonardo↗
NASA NTRS1985

A three-dimensional Navier-Stokes/Euler code for blunt-body flow computations

The formulation computation method of an improved version of the three-dimensional Navier-Stokes/Euler computation algorithm of Li (1981) for the numerical simulation of blunt-body reentry flows are discussed, and results for five sample problems are presented graphically. The vector notation of the coordinate systems is defined; the governing equations are presented in full; the Jacobian matrices and damping terms of the factorization technique (based on the alternating-direction implicit procedure of Beam and Warming, 1978) are explained; and the capabilities, limitations and proper use of the code are summarized. Examples presented include the equilibrium-flow problem for Shuttle-orbiter reentry at Mach 22 and angle of attack 40.8 deg and the perfect-gas problem of an aerobraking orbital-transfer vehicle with an ellipsoidal/60-deg cone and a toroidal sonic shoulder. The advantages of the improved code in terms of accuracy and computation time are indicated.

Li, C. P.↗
NASA NTRS1976

Accretion onto magnetized neutron stars - Structure and interchange instability of a model magnetosphere

A self-consistent model is analyzed for the spherical infall of weakly magnetized plasma into the magnetosphere of a slowly rotating, strongly magnetized neutron star. It is shown that spherical infall is probably a good approximation for X-ray sources which accrete from a stellar wind. The location of the standoff shock which halts the hypersonic infall is estimated along with the emission from the shocked layer. The location of the equilibrium magnetopause and the structure of the magnetic field within it are calculated; it is found that the magnetic poles are true cusps and that the entry of gas due to equilibrium flow across a cusp is almost certainly dominated by the interchange instability near the magnetic equator. The energy principle is applied to derive necessary conditions for the occurrence of this instability. The results indicate that the strong magnetic-pressure gradient stabilizes the gas unless moderately strong radiative cooling takes place and that the cooled plasma enters the magnetosphere as long filaments capable of moving between field lines. The rate at which the equilibrium magnetopause can 'absorb' mass and momentum is derived, the validity of the approximations employed is discussed, and the likely evolution of the sinking filaments is outlined to show that the spatial distribution of the plasma is determined mainly by the dynamics and thermodynamics of the filaments rather than the magnetic-field structure.

Arons, J.↗
NASA NTRS1994

A random distribution reacting mixing layer model

A methodology for simulation of molecular mixing and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and non-reacting shear layer present in the facility given basic assumptions about turbulence properties.

Jones, Richard A.↗
NASA NTRS1994

A random distribution reacting mixing layer model

A methodology for simulation of molecular mixing, and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Research Center Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and nonreacting shear layer present in the facility given basic assumptions about turbulence properties.

Jones, Richard A.↗
NASA NTRS1993

An extended supersonic combustion model for the dynamic analysis of hypersonic vehicles

The development of an advanced dynamic model for aeroelastic hypersonic vehicles powered by air breathing engines requires an adequate engine model. This report provides a discussion of some of the more important features of supersonic combustion and their relevance to the analysis and design of supersonic ramjet engines. Of particular interest are those aspects of combustion that impact the control of the process. Furthermore, the report summarizes efforts to enhance the aeropropulsive/aeroelastic dynamic model developed at the Aerospace Research Center of Arizona State University by focusing on combustion and improved modeling of this flow. The expanded supersonic combustor model described here has the capability to model the effects of friction, area change, and mass addition, in addition to the heat addition process. A comparison is made of the results from four cases: (1) heat addition only; (2) heat addition plus friction; (3) heat addition, friction, and area reduction, and (4) heat addition, friction, area reduction, and mass addition. The relative impact of these effects on the Mach number, static temperature, and static pressure distributions within the combustor are then shown. Finally, the effects of frozen versus equilibrium flow conditions within the exhaust plume is discussed.

Bossard, J. A.↗
NASA NTRS1995

Effects of Sonic Line Transition on Aerothermodynamics of the Mars Pathfinder Probe

Flow field solutions over the Mars Pathfinder Probe spanning the trajectory through the Martian atmosphere at angles of attack from 0 to 11 degrees are obtained. Aerodynamic coefficients derived from these solutions reveal two regions where the derivative of pitching moment with respect to angle of attack is positive at small angles of attack. The behavior is associated with the transition of the sonic line location between the blunted nose and the windside shoulder of the 70 degree half-angle cone in a gas with a low effective ratio of specific heats. The transition first occurs as the shock layer gas chemistry evolves from highly nonequilibrium to near equilibrium, above approximately 6.5 km/s and 40 km altitude, causing the effective specific heat ratio to decrease. The transition next occurs in an equilibrium flow regime as velocities decrease through 3.5 km/s and the specific heat ratio increases again with decreasing enthalpy. The effects of the expansion over the shoulder into the wake are more strongly felt on the fustrum when the sonic line sits on the shoulder. The transition also produces a counter-intuitive trend in which windside heating levels decrease with increasing angle of attack resulting from an increase in the effective radius of curvature. Six-degree-of-freedom trajectory analyses utilizing the computed aerodynamic coefficients predict a moderate, 3 to 4 degree increase in total angle of attack as the probe, spinning at approximately 2 revolutions per minute, passes through these regions.

Gnoffo, Peter A.↗
NASA NTRS2017

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad/Separation Bolt Wedge Tests

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Panel test articles included a metallic separation bolt imbedded in the compression-pad and heat shield materials, resulting in a circular protuberance over a flat plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the non-equilibrium flow field in the facility nozzle, test box, and flow field over test articles, and comparisons with the measured calibration data.

TPS↗
NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC↗
NASA NTRS2017

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad Separation Bolt Wedge Tests

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Each panel test article included a metallic separation bolt imbedded in Orion compression-pad and heatshield materials, resulting in a circular protuberance over a flat plate. The protuberances produce complex model flowfields, containing shock-shock and shock-boundary layer interactions, and multiple augmented heating regions on the test plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the non-equilibrium flow field in the facility nozzle, test box, and flow field over test articles, and comparisons with the measured calibration data.

CFD↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of NBI-induced energetic ions on internal kink stability in HL-2M

A systematic investigation is carried out, studying the effect of the neutral beam injection induced energetic particles (EPs) on the n = 1 (n is the toroidal mode number) internal kink (IK) instability in the HL-2M tokamak, utilizing the MARS-F/K code [Liu et al., Phys. Plasmas 7, 3681 (2000) and 15, 112503 (2008)]. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code [Breslau et al. Computer Software (2018)], is chosen for this study. Compared to the fluid model, non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid computations with MARS-K show a generally stabilization effect on the IK, due to drift kinetic resonances associated with EPs. Further, the bounce resonance of trapped EPs has minor influence on the mode stability. In the absence of the plasma equilibrium flow and with the assumed particle pitch distribution, the transit resonance of co-current (countercurrent) passing EPs destabilizes (stabilizes) the IK. With plasma flow, both co- and countercurrent passing EPs tend to stabilize the mode, but the effect is stronger with the countercurrent particles. These modeling results provide useful guidance for interpreting MHD instabilities in the future high-performance experiments in HL-2M.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗