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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Twin sterile neutrino dark matter

We propose that the dark matter of our Universe could be sterile neutrinos which reside within the twin sector of a mirror twin Higgs model. In our scenario, these particles are produced through a version of the Dodelson-Widrow mechanism that takes place entirely within the twin sector, yielding a dark matter candidate that is consistent with x-ray and gamma-ray line constraints. Furthermore, this scenario can naturally avoid the cosmological problems that are typically encountered in mirror twin Higgs models. In particular, if the sterile neutrinos in the Standard Model sector decay out of equilibrium, they can heat the Standard Model bath and reduce the contributions of the twin particles to N eff . Such decays also reduce the effective temperature of the dark matter, thereby relaxing constraints from large-scale structure. The sterile neutrinos included in this model are compatible with the seesaw mechanism for generating Standard Model neutrino masses. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1967

An analysis of CH in the solar atmosphere.

CH molecule solar spectral observations interpreted using several upper photosphere models and assuming local thermodynamic equilibrium approximation

THERMODYNAMIC EQUILIBRIUM↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS1991

Addition of equilibrium air to an upwind Navier-Stokes code and other first steps toward a more generalized flow solver

An upwind three-dimensional volume Navier-Stokes code is modified to facilitate modeling of complex geometries and flow fields represented by proposed National Aerospace Plane concepts. Code enhancements include an equilibrium air model, a generalized equilibrium gas model and several schemes to simplify treatment of complex geometric configurations. The code is also restructured for inclusion of an arbitrary number of independent and dependent variables. This latter capability is intended for eventual use to incorporate nonequilibrium/chemistry gas models, more sophisticated turbulence and transition models, or other physical phenomena which will require inclusion of additional variables and/or governing equations. Comparisons of computed results with experimental data and results obtained using other methods are presented for code validation purposes. Good correlation is obtained for all of the test cases considered, indicating the success of the current effort.

Rosen, Bruce S.↗
DOE OSTI2024

Core Model Proposal #399: Updating the SSP Database (v3.0) (Population, GDP, and Labor Force) and Labor Productivity (KLEM)

This Core Model Proposal (CMP) updates the Shared Socioeconomic Pathway (SSP) database to a recent version (v3.0.1; released in 2024) within GCAM. Currently, GCAM relies on socioeconomic drivers, including population, GDP, and labor force projections, from the original SSP database version released in 2013. These projections, provided by independent socioeconomic dynamic models (e.g., multi-dimensional demographic models and macroeconomic models of convergence growth), may need regular updates when (1) near-term observations become available and (2) there are updates and advancements in the socioeconomic modeling. Timely updates of socioeconomic drivers in global economic equilibrium and multisector dynamic modeling will ensure (1) alignment of historical years and near-term projections with observations, enhancing base year calibrations, including calibration parameters and labor productivity, and (2) improvement of long-term projections with updated socioeconomic drivers, which set the scale of the economy. This CMP updates the SSP data (from v2013 to v2024) and also fixes/reconciles historical GDP data sources in GCAM. We investigate the impact of these updates on GCAM projections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Time-Varying Model for Predicting Formaldehyde Emission Rates in Homes

Recent studies have succeeded in relating emissions of various volatile organic compounds to material mass diffusion transfer using detailed empirical characteristics of each of the individual emitting materials. While significant, the resulting models are often scenario specific and/or require a host of individual component parameters to estimate emission rates. This study developed an approach to estimate aggregated emissions rates based on a wide number of field measurements. We used a multi-parameter regression model based on previous mass transfer models to predict formaldehyde emission rate for a whole dwelling using field-measured, time-resolved formaldehyde concentrations, air exchange rates, and indoor environmental parameters in 63 California single-family houses built between 2011 and 2017. The resulting model provides time-varying formaldehyde emission rates, normalized by floor area, for each study home, assuming a well-mixed mass balance transport model of the home, and a well-mixed layer transport model of indoor surfaces. The surface layer model asserts an equilibrium concentration within the surface layer of the emitted materials that is a function of temperature and RH; the dwelling ventilation rate serves as a surrogate for indoor concentration. We also developed a more generic emission model that is suitable for broad prediction of emission for a population of buildings. This model is also based on measurements aggregated from 27 homes from the same study. We showed that errors in predicting household formaldehyde concentrations using this approach were substantially less than those using a traditional constant emission rate model, despite requiring less unique building information.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

A Theory for Rapid Charging Events on the International Space Station

The Floating Potential Measurement Unit (FPMU) has detected high negative amplitude rapid charging events (RCEs) on the International Space Station (ISS) at the morning terminator. These events are larger and more rapid than the ISS morning charging events first seen by the Floating Potential Probe (FPP) on ISS in 2001. In this paper, we describe a theory for the RCEs that further elucidates the nature of spacecraft charging in low Earth orbit (LEO) in a non-equilibrium situation. The model accounts for all essential aspects of the newly discovered phenomenon, and is amenable to testing on-orbit. Predictions of the model for the amplitude of the ISS RCEs for the full set of ISS solar arrays and for the coming solar cycle are given, and the results of modeling by the Environments WorkBench (EWB) are compared to the observed events to show that the phenomenon can be explained by solar array driven charging. The situation is unique because the coverglasses have not yet reached equilibrium with the surrounding plasma during the RCEs. Finally, a prescription for further use of the ISS for investigating fundamental plasma physics in LEO is given. Already, plasma and charging monitoring instruments on ISS have taught us much about spacecraft interactions with the dense LEO plasma, and we expect they will continue to yield more valuable science when the Japanese Experiment Module (JEM) is in place.

Ferguson, Dale C.↗
NASA NTRS1997

Molecular Hydrogen in Diffuse Interstellar Clouds of Arbitrary Three-Dimensional Geometry

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen in diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photodissociation rate of H2 in such clouds was computed with a 26 ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H2 line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 micron continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 micron continuum intensity to H I column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be the largest. Using a standard chi square analysis to assess the goodness of fit of our models, we find (at the 60 a level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be much smaller than its extent in the plane of the sky. If D is the distance to G236 + 39, and given standard assumptions about the rate of grain-catalyzed H2 formation, we find that the cloud has an extent along the line of sight that is 0.9 +/- 0.1 times its mean extent projected onto the plane of the sky and a gas density of 53 +/- 8 (100 pc/D) H nuclei/cc and is illuminated by a radiation field of 1.1 +/- 0.2 (100 pc/D) times the mean interstellar radiation field. The derived 100 micron emissivity per nucleon is 1.13 +/- 0.06 x 10(exp -20) MJy/sr sq cm.

Spaans, Marco↗
NASA NTRS1974

Descartes Mountains and Cayley Plains - Composition and provenance

Trace element compositions of petrographically characterized 2-4 mm lithic fragments from Apollo 16 soil samples are used to calculate initial REE concentrations in liquids in equilibrium with lunar anorthosites and to discuss the provenance of the Cayley Formation. Lithic fragments may be subdivided into four groups: (1) ANT rocks, (2) K- and SiO2-rich mesostasis-bearing rocks, (3) poikiloblastic rocks, and (4) (spinel) troctolites. Model liquids in equilibrium with essentially monominerallic anorthosites have initial REE concentrations 5-8 times those of chondrites. The REE contents of K- and SiO2-rich mesostasis-bearing rocks and poikiloblastic rocks are dominated by the mesostasis phases. ANT rocks appear to be more abundant in the Descartes Mountains, while poikiloblastic rocks appear to be more abundant in the Cayley Plains. Poikiloblastic rocks have intermediate to high LIL-element concentrations yet the low gamma-ray activity of Mare Orientale implies low LIL-element concentrations. Consequently, it is unlikely that the Cayley Formation is Orientale ejecta. A local origin as ejecta from smaller impacts is a more plausible model for the deposition of the Cayley Formation.

Drake, M. J.↗
DOE OSTI2025

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
NASA NTRS1979

Relaxation to photochemical equilibrium of charged species within displaced air parcels in the D-region

The effects of the relaxation to photochemical equilibrium are investigated for the electrons and NO(+) ions in the upper D-region and for the electrons and several species of positive and negative ions present in the lower D-region. In the upper D-region the photochemical model includes electrons, nitric oxide, and NO(+). A self-consisting model of the equilibrium concentrations is developed for the lower D-region, which includes electrons, nitric oxide, O2(-), heavy negative and positive ions, and the oxygen species O, O3, and O2(1 delta g). At altitudes of 75 and 80 km the time dependence of the oxygen species has a significant effect on the equilibration of the electron and O2(-) concentrations. The relaxation to photochemical equilibrium proceeds by five e-folding time scales in the lower D-region, where the time required for final equilibration due to recombination is very long.

Hill, R. J.↗
NASA NTRS2013

Ablator Response Model Development: From Flight Data Back to Fundamental Experiments

The successful Mars atmospheric entry by the Mars Science Laboratory (MSL-Curiosity) combined with the success of the Earth atmospheric entry by the Stardust capsule have established PICA as a major Thermal Protection Systems (TPS) material. We expect that this class of materials will be on the short list selected by NASA for any atmospheric entry missions and that it will be the lead of that list of materials in any planning, feasibility studies or flight readiness studies. In addition to NASAs successes, the Dragon capsule, the successful commercial space vehicle built by SpaceX, uses PICA-X, while the European Space Agency is considering ASTERM for its exploration missions that involve atmospheric entries, both of these materials are of the same family as PICA. In the talk, a high-fidelity model will be detailed and discussed. The model tracks the chemical composition of the gases produced during pyrolysis. As in the conventional models, it uses equilibrium chemistry to determine the recession rate at high temperatures but switches to in-volume finite-rate ablation for lower temperatures. It also tracks the time evolution of the porosity of the material. Progress in implementing this high-fidelity model in a code will be presented. In addition, a set of basic experimental data being supported for model validation will be summarized. The validation process for the model development will be discussed. Preliminary results will be presented for a case where detailed pyrolysis product chemistry is computed. Finally, a wish list for a set of validation experiments will be outlined and discussed.

Mansour, Nagi N.↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗
NASA NTRSDate not supplied

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗
NASA NTRS1991

How much hydrogen is there in a white dwarf?

Stratified hydrogen/helium envelope models in diffusive equilibrium are calculated for a 0.6-solar-mass white dwarf for effective temperatures between 10,000 and 80,000 K in order to investigate the observational constraints placed on the total hydrogen mass. Convective mixing is included ab initio in the calculations, and synthetic spectra are used for comparing these models with observational materials. It is shown that evolutionary changes in the surface composition of white dwarfs cannot be explained by a model in which a small amount of hydrogen floats to the surface from initially being mixed in the outer parts of a helium envelope. It is pointed out that the shape of the hydrogen lines can be used for constraining theories of convective overshoot.

Macdonald, James↗
NASA NTRS2003

Multibody Parachute Flight Simulations for Planetary Entry Trajectories Using "Equilibrium Points"

A method has been developed to reduce numerical stiffness and computer CPU requirements of high fidelity multibody flight simulations involving parachutes for planetary entry trajectories. Typical parachute entry configurations consist of entry bodies suspended from a parachute, connected by flexible lines. To accurately calculate line forces and moments, the simulations need to keep track of the point where the flexible lines meet (confluence point). In previous multibody parachute flight simulations, the confluence point has been modeled as a point mass. Using a point mass for the confluence point tends to make the simulation numerically stiff, because its mass is typically much less that than the main rigid body masses. One solution for stiff differential equations is to use a very small integration time step. However, this results in large computer CPU requirements. In the method described in the paper, the need for using a mass as the confluence point has been eliminated. Instead, the confluence point is modeled using an "equilibrium point". This point is calculated at every integration step as the point at which sum of all line forces is zero (static equilibrium). The use of this "equilibrium point" has the advantage of both reducing the numerical stiffness of the simulations, and eliminating the dynamical equations associated with vibration of a lumped mass on a high-tension string.

Raiszadeh, Ben↗