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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A complementarity-based vapor-liquid equilibrium formulation for equation-oriented simulation and optimization

Vapor-Liquid Equilibrium (VLE) is a cornerstone of computer aided process engineering (CAPE). Embedded within process system models, VLE calculations are inherently procedural with non-smooth behavior that frequently require discrete decisions. Traditionally, these features resist the incorporation of VLE within efficient, large-scale equation-oriented (EO) process simulation and optimization strategies. On the other hand, recent reformulation of VLE models through the incorporation of complementarity constraints has broadened its scope to deal seamlessly with phase transitions and even supercritical excursions in process simulation and optimization. In this study, we extend these VLE complementarity models to EO frameworks where procedural thermodynamic property libraries are still required. Here, we develop an efficient, non-intrusive, and intuitive “square-flash” equation system that has been implemented within the IDAES Integrated Platform (IDAES-IP). Further, the effectiveness of this modular approach is demonstrated on case studies for non-ideal flash calculations and distillation optimization with disappearing phases and supercritical transitions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Assessing mean climate change signals in the global CORDEX-CORE ensemble

The new Coordinated Output for Regional Evaluations (CORDEX-CORE) ensemble provides high-resolution, consistent regional climate change projections for the major inhabited areas of the world. It serves as a solid scientific basis for further research related to vulnerability, impact, adaptation and climate services in addition to existing CORDEX simulations. The aim of this study is to investigate and document the climate change information provided by the CORDEX-CORE simulation ensemble, as a part of the World Climate Research Programme (WCRP) CORDEX community. An overview of the annual and monthly mean climate change information in selected regions in different CORDEX domains is presented for temperature and precipitation, providing the foundation for detailed follow-up studies and applications. Initially, two regional climate models (RCMs), REMO and RegCM were used to downscale global climate model output. The driving simulations by AR5 global climate models (AR5-GCMs) were selected to cover the spread of high, medium, and low equilibrium climate sensitivity at a global scale. The CORDEX-CORE ensemble has doubled the spatial resolution compared to the previously existing CORDEX simulations in most of the regions (25 $\,\mathrm {km}$ (0.22 $^{\circ }$ /jats:inline-formula>) versus 50 $\,\mathrm {km}$ (0.44 $^{\circ }$ )) leading to a potentially improved representation of, e.g., physical processes in the RCMs. The analysis focuses on changes in the IPCC physical climate reference regions. The results show a general reasonable representation of the spread of the temperature and precipitation climate change signals of the AR5-GCMs by the CORDEX-CORE simulations in the investigated regions in all CORDEX domains by mostly covering the AR5 interquartile range of climate change signals. The simulated CORDEX-CORE monthly climate change signals mostly follow the AR5-GCMs, although for specific regions they show a different change in the course of the year compared to the AR5-GCMs, especially for RCP8.5, which needs to be investigated further in region specific process studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Structural and thermodynamic effects of hydration in Na-zeolite A (LTA) from low-temperature heat capacity

Zeolite A (Linde Type A; LTA) is an industrially important porous mineral that has recently been shown to exhibit framework flexibility upon changes in hydration level. Here, to investigate the flexibility transition from a thermodynamic perspective, we have performed heat capacity measurements on sodium zeolite A at seven incremental hydration levels ranging from zero to equilibrium with ambient air. Excess low-frequency vibrations beyond the predictions of the Debye model are found in all samples, and the frequency of these vibrations increases as a function of hydration level. This suggests that an increase in hydration causes a decrease in at least one type of framework flexibility for sodium zeolite A. In addition, a subtle excess heat capacity contribution from 150 to 280 K is observed only for low and intermediate hydration levels, which may arise from a transformation tied to framework flexibility previously observed in zeolite A via gas absorption calorimetry. Values of the standard thermodynamic functions C p,m °, Δ 0 T S m °, Δ 0 T H m °, and Φ m ° are also reported.

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

Defect-induced formation and frustration-driven multiple magnetic transitions in Gd 2 Co 0.90 Si 2.90

A new ternary compound, Gd 2 Co 0.90 Si 2.90 , has been synthesized in a chemically single phase by deliberately introducing lattice vacancies in the Co and Si sites. The system is characterized by DC magnetization, heat capacity, resistivity and density functional theory (DFT) calculations. Here, a detailed experimental study reveals that the system exhibits an antiferromagnetic transition below T N = 11.3 K followed by spin freezing behaviour below T B ~ 4.8 K. An additional high-temperature magnetic transition could also be detected at T H ~ 150 K, which is short range in character and is associated with defect-induced polarization of the conduction electrons. The spin-glass-like state formation in the system is additionally bolstered by the observation of magnetic relaxation and associated aging phenomena as well as the magnetic memory effect. On the basis of non-equilibrium dynamical behaviour, we argue that the glassy state in this compound favours the hierarchical model over the droplet model. DFT calculations and experimental outcomes establish that the ground state is magnetically frustrated due to the small energy difference between distinct spin configurations having closely spaced ground-state energies.

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

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Topological band structure transitions and goniopolar transport in honeycomb antimonene as a function of buckling

The electronic band topology of monolayer β-Sb (antimonene) is studied from the flat honeycomb to the equilibrium buckled structure using first-principles calculations and analyzed using a tight-binding model and low-energy Hamiltonians. In flat monolayer Sb, the Fermi level occurs near the intersection of two warped Dirac cones, one associated with the pz orbitals, and one with the {p x , p y } orbitals. The differently oriented threefold warping of these two cones leads to an unusually shaped nodal line, which leads to anisotropic in-plane transport properties and goniopolarity. A slight buckling opens a gap along the nodal line except at six remaining Dirac points, protected by symmetry. Under increasing buckling, pairs of Dirac points of opposite winding number annihilate at a critical buckling angle. At a second critical angle, the remaining Dirac points disappear when the band structure opens a gap. Spin-orbit coupling and edge states are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonequilibrium DMFT + CPA for correlated disordered systems

We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. Finally, the real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Floquet insulators and lattice fermions beyond naive time discretization

Periodically driven quantum systems known as Floquet insulators can host topologically protected bound states known as “ π modes” that exhibit response at half the frequency of the drive. Such states can also appear in undriven lattice field theories when time is discretized as a result of fermion doubling, raising the question of whether these two phenomena could be connected. Recently we demonstrated such a connection at the level of an explicit mapping between the spectra of a continuous-time Floquet model and a discrete-time undriven lattice fermion model. However, this mapping relied on a symmetry of the single-particle spectrum that is not present for generic drive parameters. Inspired by the example of the temporal Wilson term in lattice field theory, in this paper we extend this mapping to the full drive parameter space by allowing the parameters of the discrete-time model to be frequency-dependent. The spectra of the resulting lattice fermion models exactly match the quasienergy spectrum of the Floquet model in the thermodynamic limit. Our results demonstrate that spectral features characteristic of beyond-equilibrium physics in Floquet systems can be replicated in static systems with appropriate time discretization. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Velocity Measurements in Channel Gas Flows in the Slip Regime by means of Molecular Tagging Velocimetry

Direct measurements of the slip velocity in rarefied gas flows produced by local thermodynamic non-equilibrium at the wall represent crucial information for the validation of existing theoretical and numerical models. In this work, molecular tagging velocimetry (MTV) by direct phosphorescence is applied to argon and helium flows at low pressures in a 1-mm deep channel. MTV has provided accurate measurements of the molecular displacement of the gas at average pressures of the order of 1 kPa. To the best of our knowledge, this work reports the very first flow visualizations of a gas in a confined domain and in the slip flow regime, with Knudsen numbers up to 0.014. MTV is cross-validated with mass flowrate measurements by the constant volume technique. The two diagnostic methods are applied simultaneously, and the measurements in terms of average velocity at the test section are in good agreement. Moreover, preliminary results of the slip velocity at the wall are computed from the MTV data by means of a reconstruction method

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
DOE OSTI2022

Molecular basis of force-pCa relation in MYL2 cardiomyopathy mice: Role of the super-relaxed state of myosin

Here, in this study, we investigated the role of the super-relaxed (SRX) state of myosin in the structure–function relationship of sarcomeres in the hearts of mouse models of cardiomyopathy-bearing mutations in the human ventricular regulatory light chain (RLC, MYL2 gene). Skinned papillary muscles from hypertrophic (HCM–D166V) and dilated (DCM–D94A) cardiomyopathy models were subjected to small-angle X-ray diffraction simultaneously with isometric force measurements to obtain the interfilament lattice spacing and equatorial intensity ratios (I 11 /I 10 ) together with the force-pCa relationship over a full range of [Ca 2+ ] and at a sarcomere length of 2.1 μm. In parallel, we studied the effect of mutations on the ATP-dependent myosin energetic states. Compared with wild-type (WT) and DCM–D94A mice, HCM–D166V significantly increased the Ca 2+ sensitivity of force and left shifted the I 11 /I 10 -pCa relationship, indicating an apparent movement of HCM–D166V cross-bridges closer to actin-containing thin filaments, thereby allowing for their premature Ca 2+ activation. The HCM–D166V model also disrupted the SRX state and promoted an SRX-to-DRX (super-relaxed to disordered relaxed) transition that correlated with an HCM-linked phenotype of hypercontractility. While this dysregulation of SRX ↔ DRX equilibrium was consistent with repositioning of myosin motors closer to the thin filaments and with increased force-pCa dependence for HCM–D166V, the DCM–D94A model favored the energy-conserving SRX state, but the structure/function–pCa data were similar to WT. Our results suggest that the mutation-induced redistribution of myosin energetic states is one of the key mechanisms contributing to the development of complex clinical phenotypes associated with human HCM–D166V and DCM–D94A mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Real-time chiral dynamics at finite temperature from quantum simulation

In this study, we explore the real-time dynamics of the chiral magnetic effect (CME) at a finite temperature in the (1+1)-dimensional QED, the massive Schwinger model. By introducing a chiral chemical potential μ 5 through a quench process, we drive the system out of equilibrium and analyze the induced vector currents and their evolution over time. The Hamiltonian is modified to include the time-dependent chiral chemical potential, thus allowing the investigation of the CME within a quantum computing framework. We employ the quantum imaginary time evolution (QITE) algorithm to study the thermal states, and utilize the Suzuki-Trotter decomposition for the real-time evolution. This study provides insights into the quantum simulation capabilities for modeling the CME and offers a pathway for studying chiral dynamics in low-dimensional quantum field theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Step in understanding the S 8 tension

Models of dark sectors with a mass threshold can have important cosmological signatures. If, in the era prior to recombination, a relativistic species becomes nonrelativistic and is then depopulated in equilibrium, there can be measurable impacts on the cosmic microwave background as the entropy is transferred to lighter relativistic particles. In particular, if this “step” occurs near z~20, 000, the model can naturally accommodate larger values of H 0 . If this stepped radiation is additionally coupled to dark matter, there can be a meaningful impact on the matter power spectrum as dark matter can be coupled via a species that becomes nonrelativistic and depleted. This can naturally lead to suppressed power at scales inside the sound horizon before the step, while leaving conventional cold dark matter signatures for power outside the sound horizon. We study these effects and show such models can naturally provide lower values of S 8 than scenarios without a step. Furthermore, this suggests these models may provide an interesting framework to address the S 8 tension, both in concert with the H 0 tension and without.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The Role of Oxygen Concentration in Enabling High Gradients in Niobium SRF Cavities

We studied the role of O concentration with depth in the performance of Nb SRF cavities. An ensemble of electropolished 1.3 GHz cavities, which initially showed high field Q-slope (HFQS), was subjected to sequential testing and treatment with in-situ low temperature baking at various temperatures. We find that increasing the bake duration causes (i) an increase in the onset of HFQS until it is absent up to quench (ii) a non-monotonic relationship with the quench field (iii) an evolution of the RBCS toward a non-equilibrium behavior that drives anti-Q slope. Our data is qualitatively explained by assuming an O diffusion model and suggests that the mitigation of HFQS that arises from 120 C in-situ LTB is mediated by the diffusion of O from the native oxide which prevents the precipitation of proximity-coupled Nb nano-hydrides, in turn enabling higher quench fields. The decrease in quench field for cavities in which O has been diffused >90 nm from the RF surface may be due to a reduction of the field limit in the SS bilayer structure. We also suggest that the evolution of the RBCS occurs due to the absence of proximity coupled inclusions, bringing about non-equilibrium effects.

Bafia, D.↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2025

Efficient Simulation of Cascading Outages Using an Energy Function-Embedded Quasi-Steady-State Model

Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.

Guo, Zhenping [Univ. of Tennessee, Knoxville, TN (↗