Dynamical analysis of non-equilibrium forced convective boiling flows
Dynamical analysis of nonequilibrium, forced convective, boiling liquid metal flows
Pan, C. H. T.↗
SEARCH · Search NASA
Results for “non-equilibrium dynamics”
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.
Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo
The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.
Rarefied Gas Dynamics↗
The fragmentation of proto-globular clusters. I - Thermal instabilities
The metal abundances among the stars within a typical globular cluster are remarkably homogeneous. This indicates that star formation in these systems was a globally coordinated event which occurred over a time span less than or comparable to the collapse time scale of the cluster. This issue is addressed by assuming that the fragmentation of a proto-globular cluster cloud proceeded in two steps. In the first step, thermal instability led to the rapid growth of initial fluctuations. This led to a large contrast in the dynamical time scales between the perturbations and the parent cloud, and the perturbations then underwent gravitational instabilities on short time scales. This process is modeled using one-dimensional hydrodynamic simulations of clouds both with and without external heat sources and self-gravity. The models include the effects of a non-equilibrium H2 abundance. The results indicate that fragmentation can occur on time scales significantly less than the dynamical time scale of the parent cloud.
Murray, Stephen D.↗
Non-equilibrium chemistry in the solar nebula and early solar system: Implications for the chemistry of comets
Theoretical models of solar nebula and early solar system chemistry which take into account the interplay between chemical, physical, and dynamical processes have great utility for deciphering the origin and evolution of the abundant chemically reactive volatiles (H, O, C, N, S) observed in comets. In particular, such models are essential for attempting to distinguish between presolar and solar nebula products and for quantifying the nature and duration of nebular and early solar system processing to which the volatile constituents of comets have been subjected. The diverse processes and energy sources responsible for chemical processing in the solar nebula and early solar system are discussed. The processes considered include homogeneous and heterogeneous thermochemical and photochemical reactions, and disequilibration resulting from fluid transport, condensation, and cooling whenever they occur on timescales shorter than those for chemical reactions.
Fegley, Bruce, Jr.↗
Aperture synthesis observations of the molecular ring in the galactic center
Reported are 88 GHz aperture synthesis observations of HCN J=1 yields 0 emission and absorption in the central 5 pc of the Galaxy. The data, taken by the Hat Creek mm-interferometer at 5" to 10" spatial and 4 km/s spectral resolution, show a complete, clumpy ring of molecular gas surrounding the ionized central 2 pc of the Galaxy. The ring is the inner edge of a larger disk extending to about 5 pc. Comparison with sub-mm line data suggests that the HCN 1-0 line is slightly optically thick and originates in subthermally populated gas. The clumpy line emission distribution reflects a combination of hydrogen volume and column density variations. The new data clearly show a close physical relation between the molecular and the ionized gas in the central cavity. The western arc appears to be the ionized inner surface of the molecular ring, and the northern arm and bar may be streamers of ionized gas falling from the ring toward the center. The dominant large scale velocity pattern of the majority of the molecular gas in the inner 5 pc is rotation. No overall radial motion of the ring greater than about 20 km/s is apparent. The rotation is perturbed in several ways; (1) there is a very large local velocity dispersion, (2) the ring shows changes in position angle and inclination (warps), (3) there is a bright, redshifted cloud which appears to be located in the western part of the ring but does not participate in the rotation. These characteristics and the high degree of clumpiness indicate a non-equilibrium configuration of short (less than or approx. 10 to the 4th power to 10 to the 5th power y) dynamical lifetime. The warping and tilting of the structure and the short dynamical lifetime make an accurate determination of equilibrium rotation velocity uncertain.
Guesten, R.↗
Thermal instabilities in protogalactic clouds
The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.
Murray, S. D.↗
Study of non-equilibrium transport phenomena
Nonequilibrium phenomena due to real gas effects are very important features of low density hypersonic flows. The shock shape and emitted nonequilibrium radiation are identified as the bulk flow behavior parameters which are very sensitive to the nonequilibrium phenomena. These parameters can be measured in shock tubes, shock tunnels, and ballistic ranges and used to test the accuracy of computational fluid dynamic (CFD) codes. Since the CDF codes, by necessity, are based on multi-temperature models, it is also desirable to measure various temperatures, most importantly, the vibrational temperature. The CFD codes would require high temperature rate constants, which are not available at present. Experiments conducted at the NASA Electric Arc-driven Shock Tube (EAST) facility reveal that radiation from steel contaminants overwhelm the radiation from the test gas. For the measurement of radiation and the chemical parameters, further investigation and then appropriate modifications of the EAST facility are required.
Sharma, Surendra P.↗
The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions
A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.
Grossman, B.↗
Infrared spectroscopy of Jupiter and Saturn
High resolution infrared spectoscopy provides unique insights into the chemistry and dynamics of the atmospheres of Jupiter and Saturn. The 5 micrometer spectral region, which is transparent to deep levels, is particularly useful for the identification of molecules that are present at very low (parts per billion) concentrations. These are tracers of convective and strongly non-equilibrium processes in the atmosphere. High resolution ground-based spectroscopy complements Voyager and Galileo measurements. Spectroscopy is sensitive to lower mixing levels for selected molecules, while the on-board mass spectrometers probe molecules that are spectroscopically inaccessible. Analysis and modeling of the 4.7 micrometer carbon monoxide in Jupiter was completed. CO is present at a mole fraction of 1.6 plus or minus 0.3 x 10 to the 9th power and concentrated in the troposphere. At this abundance, it must be convected upward from much deeper levels in Jupiter where the temperature is near 1100 K. Thus CO is a tracer of the deep atmosphere which is otherwise unobservable. The oxygen abundance in Jupiter (as measured by the CO abundance) is near solar. Chemical or physical process must deplete the major oxygen carrier, water. Germane, GeH4, was discovered on Saturn at amole fraction of 4 plus or minus 2 x 10 to the 10th power.
Knacke, Roger↗
Rarefied-continuum gas dynamics transition for SUMS project
This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.
Cheng, Sin-I↗
Theoretical investigation of non-equilibrium chemistry and optical radiation in hypersonic flow fields
Future space vehicles returning from distant missions or high earth orbits may enter the upper regions of the atmosphere and use aerodynamic drag to reduce their velocity before they skip out of the atmosphere and enter low earth orbit. The Aeroassist Flight Experiment (AFE) is designed to explore the special problems encountered in such entries. A computer code was developed to calculate the radiative transport along line-or-sight in the general 3-D flow field about an arbitrary entry vehicle, if the temperatures and species concentrations along the line-of-sight are known. The radiative heating calculation at the stagnation point of the AFE vehicle along the entry trajectory was performed, including a detailed line-by-line accounting of the radiative transport in the vacuum ultraviolet (below 200 nm) by the atomic N and O lines. A method was developed for making measurements of the haze particles in the Titan atmosphere above 200 km altitude. Several other tasks of a continuing nature, to improve the technical ability to calculate the nonequilibrium gas dynamic flow field and radiative heating of entry vehicles, were completed or advanced.
Whiting, Ellis E.↗
Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows
A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.
stagnation heat flux↗
Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows
A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.
inviscid theory↗
Fluid dynamic modeling and numerical simulation of low-density hypersonic flow
The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.
Cheng, H. K.↗
Stellar kinematics of elliptical galaxies in pairs
In both galaxy pairs Arp 166 and 3C 278 the authors find radially increasing velocity dispersions indicating a perturbed, non-equilibrium state of the galaxies after the tidal interaction. In all galaxies, the increase is most pronounced in the regions which correspond to the centers of the outer isophotes. The authors suggest a scenario in which the galaxies are strongly decelerated on their orbits during the encounter. The deceleration depends on the radial position in the perturbed galaxy and vanishes in the center of the perturbed galaxy (Spitzer, 1958). In addition, the crossing time of the stars near the center is very short, implying that the tidal perturbations can be averaged over several orbital periods (e.g., Binney and Tremaine, 1987). In consequence, the central parts are not affected by the tidal interaction while the outer parts are strongly decelerated. This leads to a displacement of the central parts of the galaxies with respect to their envelopes in an anti-symmetrical way for the two components of each galaxy pair. The motions of the central parts subsequently are opposed by dynamical friction with the surrounding envelopes. Due to dynamical friction, the density of the stars increases in the wakes of the moving central parts (Mulder, 1983). The overdensity of stars in the wakes of the moving central parts efficiently decelerates the motions of the central parts. The reaction of the stars in the overdensity regions leads to an increase of the velocity dispersion mainly along the orbits of the moving central parts. The presented observations, especially the asymmetrical luminosity profiles and the radially increasing velocity dispersions support consistently the above scenario of tidal interaction between galaxies. Further spectroscopic observations are necessary in order to investigate the degree of anisotropy in the kinematically perturbed regions.
Madejsky, Rainer↗
Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)
Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.
Bohdan O Wesely↗
A new formulation of the conservation equations of fluid dynamics
The computation of time-dependent flows has inspired a new, higher-dimensional formulation of the conservation equations of fluid dynamics in which time is treated as a fourth coordinate. The formulation is derived for a constant-density flow, and then extended to a variable-density flow by introducing a fifth, fictitious coordinate. This new coordinate can also act as a source coordinate, so that external source terms can be included. The analysis is carried out for both incompressible, stratified flow, and compressible equilibrium flow. The results are then extended to non-equilibrium and magnetohydrodynamic flows. Several applications of the new formulation to the computation of time-dependent flows are discussed.
Vinokur, M.↗
A dynamic model of filament eruptions and two ribbon flares
Two basically different models for the filament equilibrium by Kippenhahn and Schluter (1957) and Kuperus and Raadu (1974) have appeared in the literature. A further analyses by van Tend and Kuperus (1978) added the force due to the horizontal component of the background field to the Kuperus and Raadu model. In order to obtain a better model which actually describes these phenomena, the evolution of the filament has to be considered in detail. A first attempt was recently presented by Kaastra. Kaastra did not formulate the precise energy balance equations for the problem, as is done in the present work. In the present model not only the force balance, but also the energy balance of the filament is taken into account. Thus a fully closed system of equations is obtained, that describes the evolution of the filament, first in force equilibrium during the current build-up phase, then in the non-equilibrium phase before the eruption, and the eruption itself. A neutral point appears above the photospheric surface in the non-equilibrium phase, but long before the eruption. It was found that although the filament itself may be in non-equilibrium, the evolution may still be slow up to the height where the eruption takes place. The eruption of the filament itself causes a large induced electric field at the neutral point which leads to the observed flare phenomena.
Kuin, N. Paul M.↗