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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Impurity Extraction by Droplets

The goals are to model and to measure the phase equilibrium properties of a finely divided fluid containing a large number of chemically similar species. The objective is to develop an accurate, usable model for such phenomena as pollutant extraction of rain clouds, industrial separation in spray towers, and separation in emulsions. The project was designed as a hierarchy of complementary theoretical and experimental steps. A theory was developed to describe the segregation of complex impurities at the interface of a solvent. This phenomenon is important in phase behavior when a large fraction of molecules in a material are near an interface, the situation in a finely divided material. The theory will be modified to account for the effect of surface curvature on the surface tension. The study of mixtures differs from pure fluids not only because of the surface effects but also because composition differences between the droplet and the surrounding vapor can stabilize a droplet with respect to a bulk phase.

Morrison, G.↗
NASA NTRS1986

Para hydrogen equilibration in the atmospheres of the outer planets

The thermodynamic behavior of the atmospheres of the Jovian planets is strongly dependent on the extent to which local thermal equilibration of the ortho and para states of molecular hydrogen is achieved. Voyager IRIS data from Jupiter imply substantial departures of the para hydrogen fraction from equilibrium in the upper troposphere at low latitudes, but with values approaching equilibrium at higher latitudes. Data from Saturn are less sensitive to the orth-para ratio, but suggest para hydrogen fractions near the equilibrium value. Above approximately the 200 K temperature level, para hydrogen conversion can enhance the efficiency of convection, resulting in a substantial increase in overturning times on all of the outer planets. Currently available data cannot definitively establish the ortho-para ratios in the atmospheres of Uranus and Neptune, but suggest values closer to local equilibrium than to the 3.1 normal ratio. Modeling of sub-millimeter wavelength measurements of these planets suggest thermal structures with frozen equilibrium lapse rates in their convective regions.

Conrath, Barney J.↗
DOE OSTI2021

Lithium-ion battery physics and statistics-based state of health model

A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.

Crawford, Aladsair J.↗
NASA NTRS1977

Generalized computer-aided discrete time domain modeling and analysis of dc-dc converters

A generalized discrete time domain modeling and analysis technique is presented for all types of switching regulators using any type of duty-cycle controller, and operating in both continuous and discontinuous inductor current. State space techniques are employed to derive an equivalent nonlinear discrete time model that describes the converter exactly. The system is linearized about its equilibrium state to obtain a linear discrete time model for small signal performance evaluations, such as stability, audiosusceptibility and transient response. The analysis makes extensive use of the digital computer as an analytical tool. It is universal, exact and easy to use.

Lee, F. C.↗
NASA NTRS1983

Ethane ocean on Titan

Voyager I radio occultation data is employed to develop a qualitative model of an ethane ocean on Titan. It is suggested that the ocean contains 25 percent CH4 and that the ocean is in dynamic equilibrium with an N2 atmosphere. Previous models of a CH4 ocean are discounted due to photolysis rates of CH4 gas. Tidal damping of Titan's orbital eccentricity is taken as evidence for an ocean layer approximately 1 km deep, with the ocean floor being covered with a solid C2H2 layer 100 to 200 m thick. The photolytic process disrupting the CH4, if the estimates of the oceanic content of CH4 are correct, could continue for at least one billion years. Verification of the model is dependent on detecting CH4 clouds in the lower atmosphere, finding C2H6 saturation in the lower troposphere, or obtaining evidence of a global ocean.

Lunine, J. I.↗
DOE OSTI2025

Measurement of directed flow in Au + Au collisions at $\sqrt{s{NN}}$ = 19.6 and 27 GeV with the STAR event plane detector

In heavy-ion collision experiments, the global collectivity of final-state particles can be quantified by anisotropic flow coefficients (𝑣 𝑛 ) . The first-order flow coefficient, also referred to as the directed flow (𝑣 1 ), describes the collective sideward motion of produced particles and nuclear fragments in heavy-ion collisions. It carries information on the very early stage of the collision, especially at large pseudorapidity (𝜂), where it is believed to be generated during the nuclear passage time. Directed flow therefore probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to models of the pre-equilibrium stage. In 2018, the Event Plane Detector (EPD) was installed in STAR and used for the Beam Energy Scan phase-II (BES-II) data taking. The combination of EPD (2.1 < |𝜂| < 5.1) and high-statistics BES-II data enables us to extend the 𝑣 1 measurement to the forward and backward 𝜂 regions. In this paper, we present the measurement of 𝑣 1 over a wide 𝜂 range in Au + Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 27 GeV using the STAR EPD. The results of the analysis at $\sqrt{s_{NN}}$ = 19.6 GeV exhibit excellent consistency with the previous PHOBOS measurement, while elevating the precision of the overall measurement. The increased precision of the measurement also revealed finer structures in heavy-ion collisions, including a potential observation of the first-order event-plane decorrelation. Multiple physics models were compared to the experimental results. Only a transport model and a three-fluid hybrid model can reproduce a sizable 𝑣 1 at large 𝜂 as was observed experimentally. The model comparison also indicates 𝑣 1 at large 𝜂 might be sensitive to the QGP phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Model for photothermal ionization and molecular recombination during pulsed ablation of polyethylene

The surface erosion of polyethylene is studied under conditions characteristic of extremely high-rate heating expected in laser-driven and high energy density experiments. A coupled model for photoionization and volumetric photothermal degradation of polyethylene is developed that includes finite rate chemistry and molecular recombination of radicals. First, the model is used to explore the evolution of radicals and hydrocarbons during isochoric heating, for which it is found that polyethylene is not in chemical equilibrium for rates >108K/s. Then, the model is used to explore the cooperativity between photoionization and photothermal ionization during one-dimensional ablation from a pulsed heating source, accounting for coupled energy deposition, thermokinetics, thermochemistry, hydrodynamics, mechanics, and thermal conduction in the ALEGRA multi-physics code. The ablation depth per energy pulse is found to agree with an analytical model for instantaneous energy deposition and absent molecular recombination or photoionization. Otherwise, molecular recombination of radicals reduces the depth of ablation per pulse, owing to thermal stabilization on forming branched molecules. Radical formation through photoionization can be compensated for by recombination, but it accelerates ablation at high enough photoionization rates. Finally, the influence of molecular recombination on the ablation depth is found to decrease for long pulses, owing to decreased recombination reaction rates at lower vapor densities near the ablated surface. These and other findings highlight competition between the processes of scission, recombination, and photoionization during pulsed ablation of polyethylene.

Moore, N. W. (ORCID:0000000272259674)↗
NASA NTRS1991

Modern CFD applications for the design of a reacting shear layer facility

The RPLUS2D code, capable of calculating high speed reacting flows, has been adopted to design a compressible shear layer facility. In order to create reacting shear layers at high convective Mach numbers, hot air streams at supersonic speeds, render by converging-diverging nozzles, must be provided. A finite rate chemistry model is used to simulate the nozzle flows. Results are compared with one-dimensional, chemically equilibrium solutions. Additionally, a two equation turbulence model with compressibility effects has been successfully incorporated with the RPLUS code. The model has been applied to simulate a supersonic shear layer. Preliminary results show favorable comparisons with the experimental data.

Yu, S. T.↗
DOE OSTI2021

The Hubble PanCET Program: Transit and Eclipse Spectroscopy of the Strongly Irradiated Giant Exoplanet WASP-76b

Ultra-hot Jupiters with equilibrium temperatures greater than 2000 K are uniquely interesting targets as they provide us crucial insights into how atmospheres behave under extreme conditions. This class of giant planets receives intense radiation from their host star and usually has strongly irradiated and highly inflated atmospheres. At such a high temperature, cloud formation is expected to be suppressed and thermal dissociation of water vapor could occur. We observed the ultra-hot Jupiter WASP-76b with seven transits and five eclipses using the Hubble Space Telescope and the Spitzer Space Telescope (Spitzer) for a comprehensive study of its atmospheric chemical and physical processes. We detected TiO and H{sub 2}O absorption in the optical and near-infrared transit spectrum. Additional absorption by a number of neutral and ionized heavy metals like Fe, Ni, Ti, and SiO help explain the short-wavelength transit spectrum. The secondary eclipse spectrum shows muted water feature but a strong CO emission feature in Spitzer’s 4.5 μm band indicating an inverted temperature pressure profile. We analyzed both the transit and eclipse spectra with a combination of self-consistent PHOENIX models and atmospheric retrieval. Both spectra were well fitted by the self-consistent PHOENIX forward atmosphere model in chemical and radiative equilibrium at solar metallicity, adding to the growing evidence that both TiO/VO and NUV heavy metals opacity are prominent NUV-optical opacity sources in the stratospheres of ultra-hot Jupiters.

47 OTHER INSTRUMENTATION↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

A Comparative Study Using CFD to Predict Iced Airfoil Aerodynamics

WIND, Fluent, and PowerFLOW were used to predict the lift, drag, and moment coefficients of a business-jet airfoil with a rime ice (rough and jagged, but no protruding horns) and with a glaze ice (rough and jagged end has two or more protruding horns) for angles of attack from zero to and after stall. The performance of the following turbulence models were examined by comparing predictions with available experimental data. Spalart-Allmaras (S-A), RNG k-epsilon, shear-stress transport, v(sup 2)-f, and a differential Reynolds stress model with and without non-equilibrium wall functions. For steady RANS simulations, WIND and FLUENT were found to give nearly identical results if the grid about the iced airfoil, the turbulence model, and the order of accuracy of the numerical schemes used are the same. The use of wall functions was found to be acceptable for the rime ice configuration and the flow conditions examined. For rime ice, the S-A model was found to predict accurately until near the stall angle. For glaze ice, the CFD predictions were much less satisfactory for all turbulence models and codes investigated because of the large separated region produced by the horns. For unsteady RANS, WIND and FLUENT did not provide better results. PowerFLOW, based on the Lattice Boltzmann method, gave excellent results for the lift coefficient at and near stall for the rime ice, where the flow is inherently unsteady.

Chi, x.↗
DOE OSTI2020

Analytic investigation of combined effects of anisotropic thermal transport and energetic particles on stability of resistive plasma resistive wall mode

The combined effects of anisotropic thermal transport and trapped energetic particles (EPs) on the stability of the resistive plasma resistive wall mode (RPRWM) are investigated by an energy-principle based analytical model. The results qualitatively confirm that of a recent toroidal modeling study [Bai et al., Phys. Plasmas 27, 072502 (2020)], in which the thermal transport can stabilize the RPPWM depending on the parameters of both the plasma equilibrium and energetic particles. The analytical model predicts a complete stabilization of the RPRWM in highly resistive plasmas, at sufficiently high EPs’ pressure, and finite plasma flow. The stabilizing effect of thermal transport originates from its enhancement of energy dissipations associated with both the resistive layer and the trapped energetic particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

RCEMIP‐ACI: Aerosol‐Cloud Interactions in a Multimodel Ensemble of Radiative‐Convective Equilibrium Simulations

Aerosol‐cloud interactions are a persistent source of uncertainty in climate research. This study presents findings from a model intercomparison project examining the impact of aerosols on clouds and climate in convection‐permitting radiative‐convective equilibrium (RCE) simulations. Specifically, 11 different modeling teams conducted RCE simulations under varying aerosol concentrations, domain configurations, and sea surface temperatures (SSTs). We analyze the response of domain‐mean cloud and radiative properties to imposed aerosol concentrations across different SSTs. Additionally, we explore the potential impact of aerosols on convective aggregation and large‐scale circulation in large‐domain simulations. The results reveal that the cloud and radiative responses to aerosols vary substantially across models. However, a common trend across models, SSTs, and domain configurations is that increased aerosol loading tends to suppress warm rain formation, enhance cloud water content in the mid‐troposphere, and consequently increase mid‐tropospheric humidity and upper‐tropospheric temperature, thereby impacting static stability. The warming of the upper troposphere can be attributed to reduced lateral entrainment effects due to the higher environmental humidity in the mid‐troposphere. However, models do not agree on aerosol impacts on convective updraft velocity based on the preliminary examination of high‐percentiles of vertical velocity at a single mid‐troposheric layer (500 hPa). In large‐domain simulations, where convection tends to self‐organize, aerosol loading does not consistently influence self‐organization but tends to reduce the intensity of large‐scale circulation forming between convective clusters and dry regions. This reduction in circulation intensity can be explained by the increase in static stability due to the upper tropospheric warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation of Li-pellet triggered ELMs in EAST with an impurity model implemented under BOUT++ framework

A simple impurity model has been developed under the BOUT++ framework to investigate the Li-pellet triggered edge localized mode (ELM) in EAST configuration. The present impurity model decouples the ion pressure enhancement effect (IPEE) (which affects radial E X B flow shear and gyro-viscosity of deuterium ions and is implicitly included in our previous work) to make the simulation of Li pellet injection more physically reasonable. In addition, the present impurity model also includes the impurity equilibrium effect (IEE) (which induces modifications on vorticity and gyro-viscosity). The simulation results show that without IPEE a turbulent ELM induced by Li pellet is triggered by multiple peeling-ballooning modes (PBMs) rather than by a single dominant mode, and the simultaneous growth of multiple modes is conducive to reducing the time for pedestal entering the energy loss state. When a turbulent ELM occurs, the nonlinearly dominant modes undergo a secondary fast growth during fast crash phase. The finding explains the secondary increase of D α emission observed in the DIII-D Li injection experiments. During the evolution of a turbulent ELM, the continuous pedestal collapse dominates the pedestal energy loss in turbulent transport phase; however, in saturation phase, the stabilizing effects on PBMs (n ≠ 0) by n = 0 mode and radial E X B flow shear within the steepest gradient region are more prominent. It is also found that with IEE and without IPEE, the size threshold of ELM triggering and the magnitude of ELM sizes both show a good agreement with EAST and DIII-D experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗