A self-consistent model of a thermally balanced quiescent prominence in magnetostatic equilibrium in a uniform gravitational field
Explore the source record for details and available documents.
Lerche, I.↗
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.
Pluto's atmosphere - Models based on refraction, inversion, and vapor-pressure equilibrium
Viking spacecraft radio-occultation measurements indicate that, irrespective of substantial differences, the polar ice cap regions on Mars have inversions similar to those of Pluto, and may also share vapor pressure equilibrium characteristics at the surface. This temperature-inversion phenomenon occurs in a near-surface boundary layer; surface pressure-temperature may correspond to the vapor-pressure equilibrium with CH4 ice, or the temperature may be slightly higher to match the value derived from IRAS data.
Eshleman, Von R.↗
A family of dynamic models for large-eddy simulation
Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.
Carati, D.↗
Global biomass supply modeling for long-run management of the climate system
Bioenergy is projected to have a prominent, valuable, and maybe essential, role in climate management. However, there is significant variation in projected bioenergy deployment results, as well as concerns about the potential environmental and social implications of supplying biomass. Bioenergy deployment projections are market equilibrium solutions from integrated modeling, yet little is known about the underlying modeling of the supply of biomass as a feedstock for energy use in these modeling frameworks. We undertake a novel diagnostic analysis with ten global models to elucidate, compare, and assess how biomass is supplied within the models used to inform long-run climate management. With experiments that isolate and reveal biomass supply modeling behavior and characteristics (costs, emissions, land use, market effects), we learn about biomass supply tendencies and differences. The insights provide a new level of modeling transparency and understanding of estimated global biomass supplies that informs evaluation of the potential for bioenergy in managing the climate and interpretation of integrated modeling. For each model, we characterize the potential distributions of global biomass supply across regions and feedstock types for increasing levels of quantity supplied, as well as some of the potential societal externalities of supplying biomass. We also evaluate the biomass supply implications of managing these externalities. Finally, we interpret biomass market results from integrated modeling in terms of our new understanding of biomass supply. Overall, we find little consensus between models on where biomass could be cost-effectively produced and the implications. We also reveal model specific biomass supply narratives, with results providing new insights into integrated modeling bioenergy outcomes and differences. The analysis finds that many integrated models are considering and managing emissions and land use externalities of supplying biomass and estimating that environmental and societal trade-offs in the form of land emissions, land conversion, and higher agricultural prices are cost-effective, and to some degree a reality of using biomass, to address climate change.
09 BIOMASS FUELS↗
Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations
There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.
79 ASTRONOMY AND ASTROPHYSICS↗
Venus climate stability and volcanic resurfacing rates
The climate of Venus is to a large degree controlled by the radiative properties of its massive atmosphere. In addition, outgassing due to volcanic activity, exospheric escape processes, and surface/atmosphere interactions may all be important in moderating the abundances of atmospheric CO2 and other volatiles. We have developed an evolutionary climate model for Venus using a systems approach that emphasizes feedbacks between elements in the climate system. Modules for atmospheric radiative transfer, surface/atmosphere interactions, tropospheric chemistry, and exospheric escape processes have so far been developed. Climate feedback loops result from interconnections between modules, in the form of the environmental parameters pressure, temperature, and atmospheric mixing ratios. The radiative transfer module has been implemented by using Rosseland mean opacities in a one dimensional grey radiative-convective model. The model has been solved for the static (time independent) case to determine climate equilibrium points. The dynamics of the model have also been explored by employing reaction/diffusion kinetics for possible surface atmosphere heterogeneous reactions over geologic timescales. It was found that under current conditions, the model predicts that the climate of Venus is at or near an unstable equilibrium point. The effects of constant rate volcanism and corresponding exsolution of volatiles on the stability of the climate model were also explored.
Bullock, M. A.↗
On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow
In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.
Sweby, P. K.↗
Computational investigation of water glasses using machine-learning potentials
The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Validation of ablation model for polyethylene using pulsed x-ray and proton exposures
The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.
Moore, N. W. (ORCID:0000000272259674)↗
Rotating cylinder electrode in reactive CO 2 capture: Identifying active C species via transport, VLE models and kinetics
Here, this article explores technical challenges and potential methodologies for understanding electrochemical Reactive CO 2 Capture (RCC) mechanisms. RCC offers potential energy cost advantages by directly converting captured CO 2 into fuels and chemicals, unlike traditional carbon capture and utilization (CCU) processes that require sequential capture, concentration, and compression. However, direct conversion of captured CO 2 introduces complexity due to additional equilibrium buffer reactions, making it challenging to identify active species for reduction in electrochemical studies. This article discusses methods to integrate transport, thermodynamics, and kinetics concepts to identify active carbon sources in RCC. Vapor‐Liquid Equilibrium (VLE) and transport models are validated against experimental results obtained in a gastight rotating cylinder electrode reactor and are shown as useful tools for studying RCC in heterogeneous electrocatalysts across different capture agents, solvents, and temperatures. This article establishes an experimental framework for advancing research in electrochemical RCC.
Electrocatalysis↗
Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study
We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.
08 HYDROGEN↗
The Riemann disks. I - Equilibrium and secular evolution
The two-dimensional analogs to Riemann S-type ellipsoids are considered as simple fluid models for galactic bars. Their equilibrium is discussed and their evolution under the effects of viscosity, gravitational radiation, and external torques from other components of the galaxy is computed. As a bar loses angular momentum due to external torques, it generally becomes less elongated, and its figure speed may either increase or decrease.
Weinberg, M. D.↗
A model for the CI (609 micron) emission of Orion
A numerical model energy balance and chemical equilibrium in the photodissociation regions at the edge of molecular clouds is presented. The model is used to calculate the emergent intensities of the following fine-structure lines: OI (at 63, 145 microns); CI (at 609, 370 microns); C II (at 158 microns); and the low-lying rotational transitions of CO. It is shown that column densities in the range 2 x 10 to the 17th to 2 x 10 to the 18th per sq cm can be obtained for the C(+)/C/CO transition region at the edges of molecular clouds. The difference in the column densities is attributed to changes in the charge exchange reactions of C(+) with SiO and S, and to the process of carbon self-healing. It is found that the calculations are in good agreement with the observed conditions in the photodissociation regions behind Orion (1) C Ori, and near the surface of OMC 1.
Tielens, A. G. G. M.↗
The effects of nonequilibrium ionization on the radiative losses of the solar corona
The emissivity of the ions of carbon and oxygen has been recalculated for a set of solar coronal loop models with a steady state siphon flow. The ion densities were calculated from the plasma velocities, temperatures, and densities of the models, and large departures from equilibrium were found. For purposes of comparison, the emissivity was calculated with and without the approximation of ionization equilibrium. Considerable differences in the radiative loss function Lambda(T) curve between equilibrium and nonequilibrium conditions were found. The nonequilibrium Lambda(T) function was then used to solve again the steady state flow equations of the loop models. The differences in the structure of these models with respect to the models calculated adopting the Lambda(T) curve in equilibrium are discussed.
Spadaro, D.↗
Rotational state-to-state transition rate coefficients for H 2 O + H 2 O collisions at nonequilibrium conditions
Aims.The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H 2 O deviates from local thermodynamic equilibrium. Methods.A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures,T kin andT rot . Results.Upon introducing our new parameterization of the collisional rates, taking into account their dependence on bothT kin andT rot , we find a change of up to 20% in the H 2 O rotational level populations for both ortho and para-H 2 O for the part of the cometary coma where the nonequilibrium regime occurs.
Astronomy & Astrophysics↗
FLARE: field line analysis and reconstruction for 3D boundary plasma modeling
The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Sensitivity Analysis, Reduced-order Modeling, and Optimization of a Gas-Cooled Pebble Bed Reactor using Equilibrium-Core and DLOFC Performance
This work presents and applies a workflow for performing design optimization on gas-cooled pebble-bed reactors. Based on previous research, a representative equilibrium core of a pebble-bed reactor and a depressurized loss-of-forced-cooling model are created. These applications are built using the Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically utilizing Griffin, Pronghorn, and Bison. After defining design-related parameters and quantities of interest regarding reactor safety and efficiency, this multiphysics model is sampled using the MOOSE stochastic tools module. The result is a comprehensive dataset of configurations, enabling sensitivity analysis and the generation of reduced-order models. Subsequently, the dataset and reduced-order models are employed in an optimization study aimed at maximizing fuel utilization while adhering to safety and operational constraints. The optimization process leads to an improvement of fuel utilization by approximately 10\%, compared to engineering-judgment-based nominal conditions.
97 - MATHEMATICS AND COMPUTING↗
A new step forward in realistic cluster lens mass modelling: analysis of Hubble Frontier Field Cluster Abell S1063 from joint lensing, X-ray, and galaxy kinematics data
We present a new method to simultaneously and self-consistently model the mass distribution of galaxy clusters that combines constraints from strong lensing features, X-ray emission, and galaxy kinematics measurements. We are able to successfully decompose clusters into their collisionless and collisional mass components thanks to the X-ray surface brightness, as well as use the dynamics of cluster members, to obtain more accurate masses exploiting the fundamental plane of elliptical galaxies. Knowledge from all observables is included through a consistent Bayesian approach in the likelihood or in physically motivated priors. We apply this method to the galaxy cluster Abell S1063 and produce a mass model that we publicly release with this paper. The resulting mass distribution presents different ellipticities for the intra-cluster gas and the other large-scale mass components as well as deviation from elliptical symmetry in the main halo. We assess the ability of our method to recover the masses of the different elements of the cluster using a mock cluster based on a simplified version of our Abell S1063 model. Thanks to the wealth of mutliwavelength information provided by the mass model and the detected X-ray emission, we also found evidence for an ongoing merger event with gas sloshing from a smaller infalling structure into the main cluster. In agreement with previous findings, the total mass, gas profile, and gas mass fraction are all consistent with small deviations from the hydrostatic equilibrium. This new mass model for Abell S1063 is publicly available, as the lenstool extension used to construct it.
79 ASTRONOMY AND ASTROPHYSICS↗