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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A study of the physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

During this period of performance (1 Jan. 1992 - 31 Dec. 1992), the investigation of the following tasks was accomplished: (1) application of the Non-linear Force-free Field (NLFF) Model to study the active region evolution; (2) coronal dynamical responses due to emerging flux including the transition region; (3) loss of MHD equilibrium due to footpoint motions; and (4) two-dimensional MHD global coronal model: steady-state streamers. The detailed description of these studies are included.

Wu, S. T.↗
NASA NTRS2017

Validation of a Thermo-Ablative Model of Elastomeric Internal Insulation Materials

In thermo-ablative material modeling, as in many fields of analysis, the quality of the existing models significantly exceeds that of the experimental data required for their validation. In an effort to narrow this gap, a laboratory-scale internal insulation test bed was developed that exposes insulation samples to realistic solid rocket motor (SRM) internal environments while being instrumented to record real-time rates of both model inputs (i.e., chamber pressure, total surface heat flux, and radiative heat flux) as well as model outputs (i.e., material decomposition depths (MDDs) and in-depth material temperatures). In this work, the measured SRM internal environment parameters were used in conjunction with equilibrium thermochemistry codes as inputs to one-dimensional thermo-ablative models of the PBINBR and CFEPDM insulation samples used in the lab-scale test firings. The computed MDD histories were then compared with those deduced from real-time X-ray radiography of the insulation samples, and the calculated in-depth temperatures were compared with those measured by embedded thermocouples. The results of this exercise emphasize the challenges of modeling and testing elastomeric materials in SRM environments while illuminating the path forward to improved fidelity.

Martin, Heath T.↗
DOE OSTI2022

New Empirical Formulation for the Sublimational Breakup of Graupel and Dendritic Snow

Abstract Ice fragments are generated by sublimation of ice particles in subsaturated conditions in natural clouds. Conceivably, such sublimational breakup would be expected to cause ice multiplication in natural clouds. Any fragment that survives will grow to become ice precipitation that may sublimate and fragment further. As a first step toward assessing this overlooked process, a formulation is proposed for the number of ice fragments from sublimation of ice particles for an atmospheric model. This is done by amalgamating laboratory observations from previously published studies. The concept of a “sublimated mass activity spectrum” for the breakup is applied to the dataset. The number of ice fragments is determined by the relative humidity over ice and the initial size of the parent ice particles. The new formulation applies to dendritic crystals and heavily rimed particles only. Finally, a thought experiment is performed for an idealized scenario of subsaturation with in-cloud descent. Scaling analysis yields an estimate of an ice enhancement ratio of about 5 (10) within a weak deep convective downdraft of about 2 m s −1 , for an initial monodisperse population of dendritic snow (graupel) particles of 3 L −1 and 2 mm. During descent, there is a dynamic equilibrium between continual emission of fragments and their depletion by sublimation. A simplified bin microphysics parcel model exhibits this dynamical quasi equilibrium, consistent with the thought experiment. The fragments have average lifetimes of around 90 and 70 s for dendrites and graupel, respectively. Sublimational breakup is predicted to cause significant secondary ice production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Replica symmetry breaking in spin glasses in the replica-free Keldysh formalism

We show that the algebra of Parisi ultrametric matrices is recovered by the real-time, replica-free, Dyson-Keldysh equations of infinite-range quantum spin glasses in the late time glassy limit. This connects to earlier results on classical and quantum systems showing how ultrametricity emerges from the persistent slow aging dynamics of the glass phase. The stationary spin glass state thereby spontaneously breaks thermal symmetry, or the Kubo-Martin-Schwinger relation of a state in global thermal equilibrium. We describe the Keldysh path integral of the infinite-range Ising model in transverse and longitudinal fields, and in the context of the Landau expansion of the action functional, show how the long-time limit connects to the full replica symmetry breaking obtained in the equilibrium formalism. We also illustrate our formalism by applying it to the spherical quantum p p -spin model, which only exhibits one-step replica symmetry breaking.

Lang, Johannes (ORCID:0000000283533392)↗
DOE OSTI2022

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Modeling and Validation of a Propellant Mixer for Controller Design

A mixing chamber used in rocket engine testing at the NASA Stennis Space Center is modelled by a system of two nonlinear ordinary differential equations. The mixer is used to condition the thermodynamic properties of cryogenic liquid propellant by controlled injection of the same substance in the gaseous phase. The three inputs of the mixer are the positions of the valves regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the flow of conditioned propellant. Mixer operation during a test requires the regulation of its internal pressure, exit mass flow, and exit temperature. A mathematical model is developed to facilitate subsequent controller designs. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as uniqueness of the equilibrium point, feedback linearizability and local stability are shown to hold under conditions having direct physical interpretation. The existence of fixed valve positions that attain a desired operating condition is also shown. Validation of the model against real data is likewise provided.

Richter, Hanz↗
NASA NTRS1982

Planetary geology: Impact processes on asteroids

The fundamental geological and geophysical properties of asteroids were studied by theoretical and simulation studies of their collisional evolution. Numerical simulations incorporating realistic physical models were developed to study the collisional evolution of hypothetical asteroid populations over the age of the solar system. Ideas and models are constrained by the observed distributions of sizes, shapes, and spin rates in the asteroid belt, by properties of Hirayama families, and by experimental studies of cratering and collisional phenomena. It is suggested that many asteroids are gravitationally-bound "rubble piles.' Those that rotate rapidly may have nonspherical quasi-equilibrium shapes, such as ellipsoids or binaries. Through comparison of models with astronomical data, physical properties of these asteroids (including bulk density) are determined, and physical processes that have operated in the solar system in primordial and subsequent epochs are studied.

Chapman, C. R.↗
NASA NTRS2011

The Lack of Chemical Equilibrium does not Preclude the Use of the Classical Nucleation Theory in Circumstellar Outflows

Classical nucleation theory has been used in models of dust nucleation in circumstellar outflows around oxygen-rich asymptotic giant branch stars. One objection to the application of classical nucleation theory (CNT) to astrophysical systems of this sort is that an equilibrium distribution of clusters (assumed by CNT) is unlikely to exist in such conditions due to a low collision rate of condensable species. A model of silicate grain nucleation and growth was modified to evaluate the effect of a nucleation flux orders of magnitUde below the equilibrium value. The results show that a lack of chemical equilibrium has only a small effect on the ultimate grain distribution.

Paquette, John A.↗
NASA NTRS1989

Collisional-radiative non-equilibrium in high-temperature helium

A stationary, nonlinear collisional-radiative model for determination of production of electrons, excited atoms and atomic line intensities for high-temperature helium is presented. The populations of several atomic levels and the intensities of spectral and continuum lines are calculated in a wide range of conditions (Te = 12,000-22,000 K, Nt = 10 to the 12 - 10 to the 20th per cu cm. Transport of radiation is included by coupling the rate equations for production of the electrons and excited atoms with the radiation escape factors which are not constant but depend on plasma conditions.

Kunc, J. A.↗
NASA NTRS1990

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Nguyen, H. L.↗
NASA NTRS1973

Non-reacting and equilibrium chemically reacting turbulent boundary-layer flows.

Turbulent boundary-layer flows of nonreacting gases are predicted for both internal (nozzle) and external flows. Effects of favorable pressure gradients on two eddy viscosity models were studied in rocket and hypervelocity wind tunnel flows. Nozzle flows of equilibrium air with stagnation temperatures up to 10,000 K were computed. Predictions of equilibrium nitrogen flows through hypervelocity nozzles were compared with experimental data. A slender spherically blunted cone was studied at 70,000 ft altitude and 19,000 ft/sec in the earth's atmosphere. Comparisons with available experimental data showed good agreement. A computer program was developed and fully documented during this investigation for use by interested individuals.

Anderson, E. C.↗
DOE OSTI2023

TOUGH3/EOS7MG (EOS7MG) v1.1

EOS7MG is a new TOUGH3 module for modeling water, brine, and multiple gas mixtures in the subsurface multiphase flow system. The module allows mixture of as more as three gas components selecting from a list of predefined gases (CH4, C2H6, C3H8, H2S, CO2, N2, NO2, H2, and AIR). EOS7MG was developed based on EOS7C and inherited the general modeling capabilities of the TOUGH3 platform. It uses cubic equation of states for gas density, viscosity, and enthalpy calculation. The module provides the selection options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state. The subsurface flow and transport of gas and aqueous phase mixtures over a wide range of pressures and temperatures is modeled with a multiphase Darcy's Law and an accurate chemical equilibrium approach as originally implemented in EOS7C is used for calculation of the solubilities of different mass components. EOS7MG was designed to allow flexible number of components based on the number of gases included in the model. The potential application of this module may include subsurface gas storage, energy storage, CCS, enhanced gas recovery, and more.

Zhang, Keni↗
NASA NTRS2001

Modeling Feedbacks Between Water and Vegetation in the Climate System

Not only is water essential for life on earth, but life itself affects the global hydrologic cycle and consequently the climate of the planet. Whether the global feedbacks between life and the hydrologic cycle tend to stabilize the climate system about some equilibrium level is difficult to assess. We use a global climate model to examine how the presence of vegetation can affect the hydrologic cycle in a particular region. A control for the present climate is compared with a model experiment in which the Sahara Desert is replaced by vegetation in the form of trees and shrubs common to the Sahel region. A second model experiment is designed to identify the separate roles of two different effects of vegetation, namely the modified albedo and the presence of roots that can extract moisture from deeper soil layers. The results show that the presence of vegetation leads to increases in precipitation and soil moisture in western Sahara. In eastern Sahara, the changes are less clear. The increase in soil moisture is greater when the desert albedo is replaced by the vegetation albedo than when both the vegetation albedo and roots are added. The effect of roots is to withdraw water from deeper layers during the dry season. One implication of this study is that the insertion of vegetation into the Sahara modifies the hydrologic cycle so that the vegetation is more likely to persist than initially.

Miller, James R.↗
NASA NTRS1990

Dynamical evolution of twisted magnetic flux tubes. I - Equilibrium and linear stability

The three-dimensional dynamical evolution of twisted magnetic flux tubes is studied using a time-dependent magnetohydrodynamic (MHD) model. The flux tubes are intended to model solar coronal loops, and include the stabilizing effect of photospheric line tying. The model permits the complete evolution of flux tubes to be followed self-consistently, including the formation, equilibrium, linear instability, and nonlinear behavior. Starting from an initial uniform background magnetic field, a twisted flux tube is created by the application of slow, localized photospheric vortex flows. The flux tube evolves quasi-statically through sequences of equilibria with increasing twist, until it becomes linearly unstable to an ideal MHD kink mode. In this paper, the equilibrium properties and the linear stability behavior are discussed. The application of the method to the uniform-twist, Gold-Hoyle field confirms the previous stability threshold for kink instability and provides estimates of the resulting growth rate.

Mikic, Zoran↗
DOE OSTI2023

Water dissociation at the water–rutile TiO 2 (110) interface from ab initio-based deep neural network simulations

The interaction of water with TiO 2 surfaces is of crucial importance in various scientific fields and applications, from photocatalysis for hydrogen production and the photooxidation of organic pollutants to self-cleaning surfaces and bio-medical devices. In particular, the equilibrium fraction of water dissociation at the TiO 2 –water interface has a critical role in the surface chemistry of TiO 2 , but is difficult to determine both experimentally and computationally. Among TiO 2 surfaces, rutile TiO 2 (110) is of special interest as the most abundant surface of TiO 2 ’s stable rutile phase. While surface-science studies have provided detailed information on the interaction of rutile TiO 2 (110) with gas-phase water, much less is known about the TiO 2 (110)–water interface, which is more relevant to many applications. In this work, we characterize the structure of the aqueous TiO 2 (110) interface using nanosecond timescale molecular dynamics simulations with ab initio-based deep neural network potentials that accurately describe water/TiO 2 (110) interactions over a wide range of water coverages. Simulations on TiO 2 (110) slab models of increasing thickness provide insight into the dynamic equilibrium between molecular and dissociated adsorbed water at the interface and allow us to obtain a reliable estimate of the equilibrium fraction of water dissociation. We find a dissociation fraction of 22 ± 6% with an associated average hydroxyl lifetime of 7.6 ± 1.8 ns. These quantities are both much larger than corresponding estimates for the aqueous anatase TiO 2 (101) interface, consistent with the higher water photooxidation activity that is observed for rutile relative to anatase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Entanglement Hamiltonian of Many-Body Dynamics in Strongly Correlated Systems

A powerful perspective in understanding nonequilibrium quantum dynamics is through the time evolution of its entanglement content. Yet apart from a few guiding principles for the entanglement entropy,to date, much less is known about the refined characteristics of entanglement propagation. In this paper we unveil signatures of the entanglement evolving and information propagating out of equilibrium, from the view of the entanglement Hamiltonian. We investigate quantum quench dynamics of prototypical Bose-Hubbard model using state-of-the-art numerical technique combined with conformal field theory. Before reaching equilibrium, it is found that a current operator emerges in the entanglement Hamiltonian, implying that entanglement spreading is carried by particle flow. In the long-time limit the subsystem enters a steady phase, evidenced by the dynamic convergence of the entanglement Hamiltonian to the expectation of a thermal ensemble. Importantly, the entanglement temperature in steady state is spatially independent,which provides an intuitive trait of equilibrium. These findings not only provide crucial information on how equilibrium statistical mechanics emerges in many-body dynamics, but also add a tool to exploring quantum dynamics from the perspective of the entanglement Hamiltonian.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Areal Distribution of the Oxygen-Isotope Ratio in Greenland

Mean values of the oxygen-isotope ratio relative to standard mean ocean water reported for 46 sites on the Greenland ice sheet are compiled together with data on mean annual surface temperature, latitude, 6180 elevation, and mean annual shortest distance to the open ocean denoted by the 10% sea-ice concentration boundary. Stepwise regression analyses, with 6180 as the dependent variable, define two robust models. In the forward mode at the 99.9% confidence level, only temperature enters the model. In the backward mode at the 95% confidence level, only temperature, latitude, and distance to the open ocean remain in the model. Inversions of the models on the basis of 160 gridpoint locations 100 km apart in the area delimited by the surface equilibrium line produce four contoured distributions of 6"0. Two distributions are based on the bivariate model and two on the multivariate model. The second distribution for each model is obtained substituting mean annual surface-temperature values obtained from the Nimbus-7 Temperature Humidity Infrared Radiometer (THIR) database. All four distributions are considered valid, and differences between them are evaluated using contoured anomaly maps. It is suggested that the inversion of the multivariate model using THIR data provides the more reliable pattern for studies of atmospheric advection or for the derivation of ice-flow adjustments for 6180 series obtained from deep-core or ablation-zone sites.

Zwally, H. Jay↗
NASA NTRS2020

The Pliocene Model Intercomparison Project Phase 2: Large-scale Climate Features and Climate Sensitivity

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ∼400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.7 and 5.2 ∘C relative to the pre-industrial era with a multi-model mean value of 3.2 ∘C. Annual mean total precipitation rates increase by 7 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases by 4.3 ∘C over land and 2.8 ∘C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60∘ N and 60∘ S exceeding the global mean warming by a factor of 2.3. In the Atlantic and Pacific oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. There is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (equilibrium climate sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble Earth system response to a doubling of CO2 (including ice sheet feedbacks) is 67 % greater than ECS; this is larger than the increase of 47 % obtained from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea surface temperatures are used to assess model estimates of ECS and give an ECS range of 2.6–4.8 ∘C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports.

Pliocene Model Intercomparison Project↗