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 631 records · Page 35

Spectra from pair-equilibrium plasmas

A numerical model of relativistic nonmagnetized plasma with uniform temperature and electron density distributions is considered, and spectra from plasma in pair equilibrium are studied. A range of dimensionless temperature (T) greater than about 0.2 is considered. The spectra from low pair density plasmas in pair equilibrium vary from un-Comptonized bremsstrahlung spectra at Thomson cross section tau(N) much less than one to Comptonized bremsstrahlung spectra with tau(N) over one. For high pair density plasmas the spectra are flat for T greater than about one, and have broad intensity peaks at energy roughly equal to 3T for T less than one. In the latter region the total luminosity is approximately twice the annihilation luminosity. All spectra are flat in the X-ray region, in contradiction to observed AGN spectra. For dimensionless luminosity greater than about 100, the cooling time becomes shorter than the Thomson time.

Zdziarski, A. A.↗

238 U PFNS Evaluation Update [Slides]

Experimental 238 U PFNS were reviewed in detail including Chi-Nu PFNS. 7/10 data sets have been accepted. We have performed a detailed UQ using LANL code ARIADNE and templates of expected PFNS measurement uncertainties. we have included Chi-Nu PFNS in experimental data base. The extended Los Alamos model in CoH was used for prior calculation. Multiple chance fission and pre-equilibrium (with exciton model) contributions were modeled.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiative forcing of the Venus mesosphere. II - Thermal fluxes, cooling rates, and radiative equilibrium temperatures

A radiative heat-transfer model is presently used to ascertain the way in which radiative forcing contributes to the up to 20 K higher temperature of the Venus polar regions, by comparison with the tropics, in the 60-100 km mesospheric levels. Model global-mean radiative equilibrium temperatures for 55-100 km are compared with observations to show how each opacity source contributes to the thermal structure. The results obtained from latitude-dependent radiative equilibrium experiments indicate that meridional variations in radiative forcing obliterate observed mesospheric temperature gradients and yield polar temperatures up to 40 K cooler than the tropics.

Crisp, David↗

A separations and purification process for improving yields and meeting fuel contaminant specifications for high-octane gasoline produced from dimethyl-ether over a Cu/BEA catalyst

In this work, we have been developing a three-step conversion of biomass-derived syngas to methanol to dimethyl-ether (DME) to non-aromatic hydrocarbons for use as high-octane gasoline and sustainable aviation fuel. This process produces branched alkanes from DME using a Cu/BEA catalyst and is a promising alternative to other syngas conversion processes such as Fischer-Tropsch to linear alkanes and traditional ZSM-5 catalyzed methanol to aromatic gasoline. In this short article we describe some advances in our understanding related to separations and purification via the use of more detailed experimental speciation in an updated process model involving multiple phase equilibrium-based separation steps. Primary modeled reactor outlet constituents (and weight %) are: C3 and lighter hydrocarbon gases (11.1%), C4s (54.5%), H 2 (1.2%), CO 2 (2.9%), water (5.0%), unreacted DME (16.5%), methanol (2.3%), and C5+ hydrocarbons (6.4%). DME (the primary reactant) and H 2 recycle and reuse are important for the overall process efficiency, and the recycle of C4s is important to increase the C5+ yield via reactivation and homologation. Thus H 2 , C4s, and DME are targeted for recycle, while methanol and water need to be removed from the product to conform with fuel specifications. Model predictions from Aspen Plus using the NRTL-RK property method indicate a fuel composition with C5+ content of 97.1 wt%, with minor constituents: 2.4 wt% C4s, 0.3 wt% methanol, 0.1 wt% DME, 0.03 wt% water, and 0.01 wt% C3s. These ranges of minor components conform with fuel quality requirements, and the modeled product is amenable for unconstrained blending to boost gasoline octane ratings.

09 BIOMASS FUELS↗

Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

histories↗

Modeling The DLOFC Accident Scenario of HTR-PM Equilibrium Core Using NEAMS Tools

High-Temperature Gas-cooled Reactors (HTGRs) have excellent characteristics in terms of safety and high thermal efficiency, and they are gaining a large interest from the industry as a candidate of Gen-IV reactors for a wide range of applications. The High Temperature gas-cooled Reactor Pebble-bed Module project (HTR-PM) is one of these designs and where helium gas is used to cool the pebble-bed region that consists of spherical fuel elements moderated with graphite. The HTR-PM design is based on the combined experience from the German pebble-bed reactor program from the 1960s through the 1990s and the HTR-10 experience in China during the 2000s. Idaho National Laboratory has a long experience in modeling of HTGRs working in developing neutronics and thermal hydraulics tools for the proper modeling of these reactors. The neutronics code Griffin has the capability to model pebble depletion . While the thermal hydraulics code Pronghorn was developed mainly to model the pebble bed reactors with the porous media assumption. In this work, an equilibrium core Multiphysics model was developed for the HTR-PM reactor to analyze the depressurized loss of forced cooling accident scenario (DLOFC). This paper is organized as follows: First, a brief description of the reactor and model specifications are provided. Then, the developed Multiphysics model is discussed. Finally, verification results of the steady-state equilibrium core and DLOFC accident are presented followed by a summary of the conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Idealized model for stratocumulus cloud layer thickness

Equilibrium solutions are constructed using an idealized model for mixed and partially mixed boundary layers to show the dependence of the equilibrium cloud layer thickness on different physical processes. The condition when a mixed layer becomes cloud-free is reduced to a comparison between three pressure scales, associated with the surface wind and divergence, the radiative cooling and the subsaturation of the entrained air. The evaporation of layer cloud by cloud-top mixing is also formulated in terms of simple pressure scales.

Betts, Alan K.↗

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗

IDAES-PSE 2.7.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.7.0 Release Highlights New features: AutoScaler and CustomScalerBase classes: Such tools are the core of the new scaling framework being implemented in IDAES. Wider adoption of scaling tools among users will result in quicker and more robust model solutions. Scaler for equilibrium reactor and saponification properties: These scaler models are examples to follow for how to use the new scaling tools. ONNX Surrogate support from Optimization & Machine Learning Toolkit (OMLT): ONNX is an open standard format to save and load ML/AI models that is widely supported by all major frameworks. This capability makes it easier for IDAES users to create surrogate models and use them without having to support each framework individually. 1D Membrane Model for CO2 Capture and Utilization: Supports ongoing efforts for modeling and optimizing polymer membrane processes for CO2 capture and conversion into formic acid. StreamScaler unit model: Unrelated to the CustomScalerBase, this unit model allows a stream’s extensive variables to be scaled by a fixed factor. This allows streams being processed by multiple units in parallel to be scaled down to unit scale and scaled back up to process scale. Bug fixes or improvements: Scaling, EoS, Diagnostics tool, Modular Properties, tests & documentation Deprecations: Old Cubic EoS

AS↗

Statistical equilibrium calculations for silicon in early-type model stellar atmospheres

Line profiles of 36 multiplets of silicon (Si) II, III, and IV were computed for a grid of model atmospheres covering the range from 15,000 to 35,000 K in effective temperature and 2.5 to 4.5 in log (gravity). The computations involved simultaneous solution of the steady-state statistical equilibrium equations for the populations and of the equation of radiative transfer in the lines. The variables were linearized, and successive corrections were computed until a minimal accuracy of 1/1000 in the line intensities was reached. The common assumption of local thermodynamic equilibrium (LTE) was dropped. The model atmospheres used also were computed by non-LTE methods. Some effects that were incorporated into the calculations were the depression of the continuum by free electrons, hydrogen and ionized helium line blocking, and auto-ionization and dielectronic recombination, which later were found to be insignificant. Use of radiation damping and detailed electron (quadratic Stark) damping constants had small but significant effects on the strong resonance lines of Si III and IV. For weak and intermediate-strength lines, large differences with respect to LTE computations, the results of which are also presented, were found in line shapes and strengths. For the strong lines the differences are generally small, except for the models at the hot, low-gravity extreme of our range. These computations should be useful in the interpretation of the spectra of stars in the spectral range B0-B5, luminosity classes III, IV, and V.

Kamp, L. W.↗

A class of analytic solutions for the thermally balanced magnetostatic prominence sheet

A theoretical study is presented for the nonlinear interplay between magnetostatic equilibrium and energy balance in a Kippenhahn-Schlueter type solar prominence sheet. A class of theoretical models is presented, expressed in closed analytic forms, thus facilitating the direct illustration of the nonlinear physical properties. The model couples the equilibrium between magnetic field, plasma pressure, and weight on the one hand, with the balance between a rho-squared T radiative loss, a rho wave heating (where rho equals plasma density, and T equals plasma temperature), and thermal conduction channeled along magnetic field lines on the other. The steady solutions are divided into three classes, and are characterized by the total wave heating in the prominence sheet which is greater than, equal to, or less than the total radiative loss. The compaction of the plasma along the field lines, under its own weight, and the energy transport effects determine which of the three basic behaviors obtains in a particular situation. A discussion is presented of the implications of the steady solutions for the formation of prominences.

Low, B. C.↗

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗

Phase transition to blob-hole coherent structure in the Hasegawa–Mima model for plasmas

An equilibrium statistical mechanics theory for the Hasegawa–Mima equations of toroidal plasmas, with canonical constraint on energy and microcanonical constraint on potential enstrophy, is solved exactly as a spherical model. The use of a canonical energy constraint instead of a fixed-energy microcanonical approach is justified by the preference for viewing real plasmas as an open system. A significant consequence of the results obtained from the partition function, free energy and critical temperature, is the condensation into a ground state exhibiting a blob-hole-like structure observed in real plasmas.

Physics↗

Computer codes for the evaluation of thermodynamic properties, transport properties, and equilibrium constants of an 11-species air model

The computer codes developed provide data to 30000 K for the thermodynamic and transport properties of individual species and reaction rates for the prominent reactions occurring in an 11-species nonequilibrium air model. These properties and the reaction-rate data are computed through the use of curve-fit relations which are functions of temperature (and number density for the equilibrium constant). The curve fits were made using the most accurate data believed available. A detailed review and discussion of the sources and accuracy of the curve-fitted data used herein are given in NASA RP 1232.

Thompson, Richard A.↗