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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗
DOE OSTI2022

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

BEPAM-E Model Code and CABBI Simulation Results for "Repeal of the Clean Power Plan: Social Cost and Distributional Implications"

The dataset consists of results and various input data that are used in the GAMS model for the publication "Repeal of the Clean Power Plan: Social Cost and Distributional Implications". All the data are either excel files or in the .inc format which can be read within GAMS or Notepad. Main data sources include: agriculture, transportation and electricity data. Model details can be found in the paper and the GAMS model package.

carbon abatement↗
DOE OSTI2022

New equations of state for saturated tuff based on the CHEST model

The CHemical Equilibrium for Saturated Tuff (CHEST) model was originally presented in Ref. 1, with numerical results depicted in Figure 1 and tabulated in Appendix B of that work. The results consisted of specific internal energy as a function of pressure and specific volume, E(P, V), spanning the domain P = 10 –4 – 10 4 kbar and V = 0.17 – 752 cm 3 /g (ρ = 1.33 × 10 –3 – 5.9 g/cm 3 ).

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

Dynamic Linking of Upstream Energy and Freight Demands for Bio and Fossil Energy Pathways in the Global Change Analysis Model

Comprehensive study of the environmental impacts associated with demand for an energy resource or carrier in any one sector requires a full consideration of the direct and indirect impacts on the rest of the regional and global energy system. It is important to consider the energy that is consumed in producing primary energy resources and energy consumed in transporting these energy resources internationally and domestically in order to produce a complete picture of these impacts. Biofuels are especially complex since they have feedbacks not just to the energy system but also to regional and global crop markets. Different modeling and analysis strategies have been applied to consider these impacts, which occur “upstream” of final energy consumption, with some success. Traditional life cycle models allow for rich technological detail in linking upstream impacts of energy demand with downstream final consumption; however, they cannot by themselves represent dynamic economic feedbacks across multiple sectors and regions over time. Computable general equilibrium (CGE) modeling does account for economic feedback among all sectors and regions; however, CGE modeling does not lend itself to physical representation of technological detail. Partial equilibrium (PE) modeling is a heterogeneous category describing economic models that focus on a subset of the economy, and they differ in coverage and ease in incorporating more sectors and economic links. In this study, we present a strategy for dynamically including the direct and indirect impacts of energy demand with physical and technology detail by explicitly adding these upstream energy and transportation links to GCAM, a PE model of global energy, land use, and emissions. We incorporate the following inter-sectoral linkages: energy inputs to crop production, energy inputs to fossil resource production, and freight transport requirements of energy and agricultural commodities. We assess the implications of explicitly including these links by measuring the global impacts of increased corn ethanol demand in the United States with and without these links included. Although the net global impact of the upstream links on energy and emissions are relatively modest in the scenarios we studied, the inclusion of these links illustrates interesting trade-offs in energy and transportation demand among fossil fuel and agriculture sectors. These sectoral interactions suggest that this level of modeling detail could be important in evaluating future analytical questions.

Sampedro Martinez de Estivariz, Jon↗
NASA NTRS2001

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.↗
DOE OSTI2025

Temperature sensitivity of the equilibrium neutronics and accident analysis of the HTR-10

Pebble-Bed High-Temperature Gas-cooled Reactors (PB-HTGR) are moderated by the graphite in the fuel pebbles and the graphite reflector surrounding the pebble-bed. Because graphite is by far the most abundant material in PB-HTGRs and the primary moderator, accurate modeling of the graphite material, including density, impurities, and temperatures, is crucial for accurate computational modeling and simulation of these reactors. While main characteristics of the graphite components are often known, the local temperature is less well known and often averaged over all components. Here, this work studies the impact of considering accurate temperature profiles in the graphite material on the generation of a small PB-HTGR model at the state of equilibrium operation and on short-term accident progression. The fuel compositions for the PB-HTGR were determined using a jump-in equilibrium modeling method, the Axial Radial Zone Equilibrium Modeling (AR-ZEM) method. In contrast to previous work, the AR-ZEM method was used considering thermal-hydraulic feedback from the MELCOR code to determine temperatures of the fuel pebbles and the surrounding graphite reflector. The consideration of an axial and radial temperature profile in the core and reflector, as opposed to uniform material temperatures, had an impact of almost 1,300 pcm on the equilibrium core eigenvalue and caused significant differences in the discharged plutonium fuel inventory with up to 4.9% and 11.0% for Pu-239 and Pu-242, respectively. To assess the impact on short-term accident progression, two Anticipated Transient Without SCRAM (ATWS) events, a Pressurized Loss of Forced Coolant (PLOFC) and a Control Rod Withdrawal (CRW) with loss of flow, were simulated with MELCOR. The use of temperature profiles in the equilibrium core models did not reveal a significant impact on the temperature, power, or reactivity responses during the transients. In conclusion, a need for consideration of accurate temperature profiles, in particular for the graphite reflector, was found for the generation of equilibrium PB-HTGRs core models using jump-in methods, but detailed temperature profiles may not be necessary when performing conservative transient analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Self-consistent and detailed opacities from a non-equilibrium average-atom model

Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
NASA NTRS1985

Oxidation states of iron in the terrestrial planets: Evidence and implications for accretion models

Analyses of Mars spectra reveal that primary minerals may be more iron rich than average basalts on Earth, and the oxidation state of iron in the pyroxines suggests equally high or higher oxygen fugacities on Mars than Earth. Analysis of Mercury spectra reveal that silicate iron contents are substantially lower than those on Mars. Mercury, however, probably has a substantial iron core. This is consistent with a lower oxidation state of iron on Mercury than on Mars. These findings are consistent with predictions of models of equilibrium condensation and homogeneous accretion. Those models predict that Mercury would have negligible Fe(2+) in silicates (Fe/Fe+Mg approx. 0.5). For Mercury iron should be virtually all metallic, while for Mars the iron would be all oxidized as FeS and silicate. Surface measurements of Fe/Fe+Mg on Venus, Earth, and Mars further support the predicted relative increase in Fe/Fe+Mg ratio (increased oxidation state) with distance from the Sun. The implied high oxygen fugacities on Mars (Fe(3+) in silicates) provide additional consistence with the model predictions.

Huguenin, R. L.↗
NASA NTRS1972

Equilibrium configurations of Vlasov plasmas carrying a current component along an external magnetic field.

A model of equilibrium configurations of Vlasov plasmas is considered which represents a combination of the models of Harris (1962) and Nicholson (1963). These plasma configurations carry a current component along an external magnetic field. The considered slab model contains a diamagnetic current and a field-aligned current for an arbitrary ratio of particle pressure to magnetic pressure of the applied constant field. For a fixed pressure ratio and field-aligned current, the model admits a family of equilibrium solutions in which the diamagnetic currents range from zero to a maximum value. The amount of diamagnetic current flowing in a machine depends on the width of the machine, the field-aligned current and other plasma parameters.

Kan, J. R.↗
NASA NTRS1989

Model formulation of non-equilibrium gas radiation for hypersonic flight vehicles

Several radiation models for low density nonequilibrium hypersonic flow are studied. It is proposed that these models should be tested by the 3-D VRFL code developed at NASA/JSC. A modified and optimized radiation model may be obtained from the testing. Then, the current VRFL code could be expanded to solve hypersonic flow problems with nonequilibrium thermal radiation.

Chang, Ing↗
NASA NTRS2019

Computational Modeling for Non-Equilibrium Shock Tube Flows

Time accurate simulation of non-equilibrium flows inside shock tube facilities presents several challenges from both physical and mathematical aspects. Furthermore, the large computational cost makes it impractical to support a real-time experimental test campaign. In this work, we explore other methods for modeling the shock tube problem with the main focus on the post-shock region and the absolute radiation emanating from it. The proposed alternative approach is several orders of magnitude less computationally expensive while still accurate enough with regards to the quantities of interest. Excellent agreement is found with the established stagnation-line approach. Comparison with time-accurate simulations shows good agreement close to the peak values and disagreement of the temperatures relaxation and radiance profiles toward equilibrium.

Bensassi, Khalil↗
NASA NTRS2019

Computational Modeling for Non-equilibrium Shock Tube Flows

Time accurate simulation of non-equilibrium flows inside shock tube facilities presents several challenges from both physical and mathematical aspects. Furthermore, the drastic computational cost makes it non-practical to support real-time experimental test campaign. In this work, we explore other methods for modeling the shock tube prob- lem with the main focus on the post-shock region and the absolute radiation emanating from it. The proposed alternative approach is several orders of magnitude less computa- tionally expansive while still accurate enough with regards to the quantities of interest. Excellent agreement is found with the well-established stagnation-line approach. Comparison with the time-accurate simulation shows good agreement close to the peak values and disagreement of the temperatures relaxation and radiance profiles toward equilibrium, due to shock speed unsteadiness.

Bensassi, Khalil↗
DOE OSTI2020

Modified Isotherm Modeling to Predict Cs Exchange with Crystalline Silicotitanate in Tank Waste Simulants

The U.S. Department of Energy is working to expedite processing of Hanford tank waste supernate at the Hanford Waste Treatment and Immobilization Plant. To support this goal, Washington River Protection Solutions is designing a Tank Side Cesium Removal (TSCR) system for suspended solids and cesium (Cs/ 137 Cs) removal from Hanford tank waste supernate. The ion exchange media selected for Cs removal at TSCR is crystalline silicotitanate (CST) that is manufactured in a nearly spherical form by Honeywell UOP LLC (UOP; Des Plaines, IL) as product IONSIV® R9140-B (Na form). The Zheng Anthony Miller (ZAM) isotherm model (Zheng et al. 1997) is a multicomponent ion exchange model used to predict the exchange of Group I metals onto CST. The ZAM isotherm has historically been used to predict Cs distribution values from Hanford and Savannah River Site (SRS) tank waste simulants. Figure S.1 summarizes model predictions from the ZAM isotherm that indicate poor prediction of Cs distribution values for simple and complex simulants with the engineered form of CST, IONSIV® R9140-B, and IONSIV® R9120-B where the solid line indicates a perfect fit by the model. The dotted lines indicate ±20% error. Batch contact testing with Hanford tank waste complex and simple simulants was used in conjunction with SRS simulants to experimentally determine Cs distribution values using a modification to the original isotherm model. The experimentally determined maximum Cs capacity for IONSIV® R9140-B CST in both the simple and complex matrices was found to be 0.53±0.3 mmoles Cs/g of CST. This value is not drastically different from the maximum Cs capacity of 0.58 mmoles Cs/g TAM-5 reported by Zheng et al. (1997). However, it is important to note that TAM-5 (commercially IONSIV® IE-910) is a powder. Hamm et al. (2002) determined that a dilution factor was needed to account for the Zr(OH)2 binder in the engineered form of CST. Hamm et al. determined that a dilution factor of 0.68 was appropriate to account for binder contribution and correct overprediction of Cs exchange on the engineered form of CST in ZAM calculations. This reduced the total capacity from 0.58 mmol/g with TAM-5 to 0.39 mmol/g for the engineered form of CST (Hamm et al. 2002). Despite substituting the experimentally determined maximum Cs capacity of 0.55 mmoles Cs/g for the literature-reported capacity of 0.39 mmoles Cs/g, it was determined that additional modifications to the model’s equilibrium rate constants were necessary in refining the isotherm model. The modified model overpredicted K+ uptake by the CST when compared to digested CST results described by Campbell et al. (2019). The modified model was further revised to omit three of the five K+ exchange equilibrium reactions described by ZAM to reduce the additional K+ loading seen by the model. Figure S.2 summarizes the revised model isotherm predictions plotted against measured Kd values.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗