Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 235 records · Page 13

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics↗

Assessment of Pioneer Venus Entry Heating with Coupled Radiation and Ablation

The objective of this work is to reconstruct thermocouple temperature measurements taken at the stagnation point of the Pioneer Venus large probe using uncoupled and coupled ablation calculations. The time history of the large probe forebody entry heating was obtained from a trajectory-based analysis that employed a material response solver to model equilibrium and nonequilibrium ablation chemistry at the surface. In the uncoupled approach, equilibrium ablation chemistry was assumed in the material response calculations, but the material response outputs were not returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer due to ablation (e.g., convective blockage and radiative energy absorption) did not influence the flowfield solutions. In the coupled approach, nonequilibrium ablation chemistry was performed and the material response outputs were returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer were accounted for in the flowfield solutions. Both uncoupled and coupled approaches produced reconstructed temperatures that were within 10% of the flight measurements through the entry heat pulse, with the coupled results comparing more favorably with the measurements. Overall, the reconstructed thermocouple temperatures in this work significantly improved on the results of previous reconstructions for the large probe.

Jannuel V V Cabrera↗

Bernal model - A simple equilibrium theory of close-packed liquids.

The Bernal model of a hard-sphere liquid is used in conjunction with an extension of the cell method to predict the behavior of liquid argon near melting. The entropy of disorder associated with the Bernal state is deduced. The model is found to be applicable to liquid metals.

Caron, L. G.↗

Computer program to solve two-dimensional shock-wave interference problems with an equilibrium chemically reacting air model

The computer program EASI, an acronym for Equilibrium Air Shock Interference, was developed to calculate the inviscid flowfield, the maximum surface pressure, and the maximum heat flux produced by six shock wave interference patterns on a 2-D, cylindrical configuration. Thermodynamic properties of the inviscid flowfield are determined using either an 11-specie, 7-reaction equilibrium chemically reacting air model or a calorically perfect air model. The inviscid flowfield is solved using the integral form of the conservation equations. Surface heating calculations at the impingement point for the equilibrium chemically reacting air model use variable transport properties and specific heat. However, for the calorically perfect air model, heating rate calculations use a constant Prandtl number. Sample calculations of the six shock wave interference patterns, a listing of the computer program, and flowcharts of the programming logic are included.

Glass, Christopher E.↗

Modeling the Spin Equilibrium of Neutron Stars in LMXBs Without Gravitational Radiation

In this paper we discuss the spin-equilibrium of accreting neutron stars in LMXBs. We demonstrate that, when combined with a naive spin-up torque, the observed data leads to inferred magnetic fields which are at variance with those of galactic millisecond radiopulsars. This indicates the need for either additional spin-down torques (eg. gravitational radiation) or an improved accretion model. We show that a simple consistent accretion model can be arrived at by accounting for radiation pressure in rapidly accreting systems (above a few percent of the Eddington accretion rate). In our model the inner disk region is thick and significantly sub-Keplerian, and the estimated equilibrium periods are such that the LMXB neutron stars have properties that accord well with the galactic millisecond radiopulsar sample. The implications for future gravitational-wave observations are also discussed briefly.

Andersson, N.↗

More accurate predictions with transonic Navier-Stokes methods through improved turbulence modeling

Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Johnson, Dennis A.↗

Application of space-marching methods to hypersonic forebody flow fields

The heavy reliance on CFD in the design of hypersonic flight vehicles emphasizes the importance of choosing computationally efficient techniques. As a result, the development and application of space-marching flow solvers for the parabolized Navier-Stokes equations have received considerable attention in the past decade. This paper summarizes a space-marching approach used at NASA Ames Research Center for the analysis of hypersonic vehicle forebodies. Aspects of the simulation process are addressed that contain features unique to space-marching flow solvers, such as grid generation, initial condition specification, and finite-rate chemistry models. For completeness, sections are also included on turbulence modeling and equilibrium air chemistry modeling. The discussion includes an overview of the computational methods utilized as well as a number of relevant results.

Lawrence, Scott L.↗

Evidence for a temperature rise in the outer layers of alpha Lyrae, from Copernicus observations of Lyman-alpha

Copernicus satellite observations of the Ly-alpha profiles in alpha Lyrae (Vega) are used to determine whether classical radiative-equilibrium LTE model atmospheres can fit the thermal structure in the outer layers of that star. Two plane-parallel LTE model photospheres of alpha Lyrae are considered: a line-blanketed radiative-equilibrium model with an effective temperature of 9650 K and log g of 4.05, and the same model with a temperature of 9500 K and log g of 4.0. The profiles of the Ly-alpha wings are computed, and it is found that classical LTE models are unable to predict either the observed violet wing or the red wing longwards of 1239 A, regardless of the line source function. It is concluded that the electron temperature must increase outwards over the surface value reached in radiative equilibrium.

Praderie, F.↗

Modeling of gamma/gamma-prime phase equilibrium in the nickel-aluminum system

A theoretical model is proposed for the determination of phase equilibrium in alloys, taking into consideration dissimilar lattice parameters. Volume-dependent pair interactions are introduced by means of phenomenological Lennard-Jones potentials and the configurational entropy of the system is treated in the tetrahedron approximation of the cluster variation method. The model is applied to the superalloy-relevant, nickel-rich, gamma/gamma-prime phase region of the Ni-Al phase diagram. The model predicts reasonable values for the lattice parameters and the enthalpy of formation as a function of composition, and the calculated phase diagram closely approximates the experimental diagram.

Sanchez, J. M.↗

Ionization and excitation in cool giant stars. I - Hydrogen and helium

The influence that non-LTE radiative transfer has on the electron density, ionization equilibrium, and excitation equilibrium in model atmospheres representative of both oxygen-rich and carbon-rich red giant stars is demonstrated. The radiative transfer and statistical equilibrium equations are solved self-consistently for H, H(-), H2, He I, C I, C II, Na I, Mg I, Mg II, Ca I, and Ca II in a plane-parallel static medium. Calculations are made for both radiative-equilibrium model photospheres alone and model photospheres with attached chromospheric models as determined semiempirically with IUE spectra of g Her (M6 III) and TX Psc (C6, 2). The excitation and ionization results for hydrogen and helium are reported.

Luttermoser, Donald G.↗

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗

Vapor-liquid equilibrium thermodynamics of N2 + CH4 - Model and Titan applications

A thermodynamic model is presented for vapor-liquid equilibrium in the N2 + CH4 system, which is implicated in calculations of the Titan tropospheric clouds' vapor-liquid equilibrium thermodynamics. This model imposes constraints on the consistency of experimental equilibrium data, and embodies temperature effects by encompassing enthalpy data; it readily calculates the saturation criteria, condensate composition, and latent heat for a given pressure-temperature profile of the Titan atmosphere. The N2 content of condensate is about half of that computed from Raoult's law, and about 30 percent greater than that computed from Henry's law.

Thompson, W. R.↗

Re-Entry Aeroheating Analysis of Tile-Repair Augers for the Shuttle Orbiter

Computational re-entry aerothermodynamic analysis of the Space Shuttle Orbiter s tile overlay repair (TOR) sub-assembly is presented. Entry aeroheating analyses are conducted to characterize the aerothermodynamic environment of the TOR and to provide necessary inputs for future TOR thermal and structural analyses. The TOR sub-assembly consists of a thin plate and several augers and spacers that serve as the TOR fasteners. For the computational analysis, the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) is used. A 5-species non-equilibrium chemistry model with a finite rate catalytic recombination model and a radiation equilibrium wall condition are used. It is assumed that wall properties are the same as reaction cured glass (RCG) properties with a surface emissivity of epsilon = 0.89. Surface heat transfer rates for the TOR and tile repair augers (TRA) are computed at a STS-107 trajectory point corresponding to Mach 18 free stream conditions. Computational results show that the average heating bump factor (BF), which is a ratio of local heat transfer rate to a design reference point located at the damage site, for the auger head alone is about 1.9. It is also shown that the average BF for the combined auger and washer heads is about 2.0.

Mazaheri, Ali R.↗

Internal structure and shape of Mimas

The observed triaxial ellipsoid shape of Mimas, in conjunction with a mean density of 1137 kg/cu cm, is presently used as the basis of hydrostatic equilibrium and nonhydrostatic models of internal structure whose outer layers are sufficiently strong to support part of the tidal load. Such plausible internal structures are consistent with the shape and density constraints as (1) a two-layer model in hydrostatic equilibrium, (2) a hydrostatic model with four layers of equal depth, and (3) a nonhydrostatic constant-density model with an elastic lithosphere; these models admit either homogeneous or heterogeneous accretion, as well as either a cold thermal history or melting and differentiation.

Ross, M. N.↗

Equilibrium Climate Sensitivity Estimated by Equilibrating Climate Models

The methods to quantify equilibrium climate sensitivity are still debated. We collect millennial length simulations of coupled climate models and show that the global mean equilibrium warming is higher than those obtained using extrapolation methods from shorter simulations. Specifically, 27 simulations with 15 climate models forced with a range of CO2 concentrations show a median 17% larger equilibrium warming than estimated from the first 150 years of the simulations. The spatial patterns of radiative feedbacks change continuously, in most regions reducing their tendency to stabilizing the climate. In the equatorial Pacific, however, feedbacks become more stabilizing with time. The global feedback evolution is initially dominated by the tropics, with eventual substantial contributions from the midlatitudes. Time dependent feedbacks underscore the need of a measure of climate sensitivity that accounts for the degree of equilibration, so that models, observations, and paleo proxies can be adequately compared and aggregated to estimate future warming.

Maria Rugenstein↗

Turbulence modeling for non-equilibrium flows

Two projects are reported. The first is the development and testing of an eddy viscosity transport model. This project also is a starting point for our work on developing computational tools for solving turbulence models in complex geometries. The second project is a stochastic analysis of the realizability of Reynolds stress transport models.

Durbin, Paul A.↗