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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Verification of a quasi-linear model for gyrokinetic turbulent transport

The verification and calibration of a new quasi-linear transport model with a large database of gyrokinetic turbulence simulations is presented in this paper. In a previous paper (Staebler et al 2020 Plasma Phys. Control. Fusion 63 015013), a model for the saturated spectrum of electric potential fluctuations was developed based on the properties of the non-linear 3D spectrum. In this paper, a modification to the overall multiplicative factor in this model is found to be necessary to improve the fit to scans of the temperature and density gradients. The error in the fit of the quasi-linear fluxes of electron and ion energy fluxes is significantly better than for previous saturation models. The spectral shift model for the impact of equilibrium E × B velocity shear (Staebler et al 2013 Phys. Rev. Lett. 110 055003) and the zonal flow mixing model for electron-scale turbulence (Staebler et al 2016 Phys. Plasmas 23 062518) are both revised to be compatible with this new model. Furthermore, the models for the loss of bounce averaging and electron collisions in the TGLF reduced linear equations (Staebler et al 2005 Phys. Plasmas 12 102508) are also changed to improve the linear eigenmodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Modeling the Chromosphere and Transition Region of Planet-hosting Star GJ 436

Ahead of upcoming space missions intending to conduct observations of low-mass stars in the ultraviolet (UV) spectral region it becomes imperative to simultaneously conduct atmospheric modeling from the UV to the visible (VIS) and near-infrared (NIR). Investigations on extended spectral regions will help to improve the overall understanding of the diversity of spectral lines arising from very different atmospheric temperature regions. Here we investigate atmosphere models with a chromosphere and transition region for the M2.5V star GJ 436, which hosts a close-in Hot Neptune. The atmosphere models are guided by observed spectral features from the UV to the VIS/NIR originating in the chromosphere and transition region of GJ 436. High-resolution observations from the Hubble Space Telescope and Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs (CARMENES) are used to obtain an appropriate model spectrum for the investigated M dwarf. We use a large set of atomic species considered in nonlocal thermodynamic equilibrium conditions within our PHOENIX model computations to approximate the physics within the low-density atmospheric regions. In order to obtain an overall match for the nonsimultaneous observations, it is necessary to apply a linear combination of two model spectra, where one of them better reproduces the UV lines while the other better represents the lines from the VIS/NIR range. This is needed to adequately handle different activity states across the observations.

stars↗
DOE OSTI2023

Imaging Field–Driven Melting of a Molecular Solid at the Atomic Scale

Solid–liquid phase transitions are basic physical processes, but atomically resolved microscopy has yet to capture their full dynamics. A new technique is developed for controlling the melting and freezing of self–assembled molecular structures on a graphene field–effect transistor (FET) that allows phase–transition behavior to be imaged using atomically resolved scanning tunneling microscopy. This is achieved by applying electric fields to 2,3,5,6–tetrafluoro–7,7,8,8–tetracyanoquinodimethane–decorated FETs to induce reversible transitions between molecular solid and liquid phases at the FET surface. Nonequilibrium melting dynamics are visualized by rapidly heating the graphene substrate with an electrical current and imaging the resulting evolution toward new 2D equilibrium states. An analytical model is developed that explains observed mixed–state phases based on spectroscopic measurement of solid and liquid molecular energy levels. The observed nonequilibrium melting dynamics are consistent with Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A modified thermodynamic sea ice model and its application

This paper presents a modified thermodynamic sea ice model based on the Winton’s three-layer model framework. Several improvements focused on the vertical thermodynamics have been made in the model. Results from a series of one-dimensional experiments show that equilibrium ice thickness in our modified model is increased by about 45 cm when compared with the original Winton’s model. Sensitivity tests indicate that all modifications mentioned above contribute to this change of ice thickness, among which the increase of ice layer has the most significant effect followed by the sea ice conductivity parameterization. The modified thermodynamic sea ice model is then coupled into the sea ice component SIS in GFDL MOM4 to exert long-term integration experiments forced by CORE2.0 data. Results indicate that the modified model shows significant improvement in simulating the Arctic sea ice, including an increase in both the sea ice thickness and volume over the whole Arctic region. This confirms the above founding from 1-D simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Geometric transformation and three-dimensional hopping of Hopf solitons

Arising in many branches of physics, Hopf solitons are three-dimensional particle-like field distortions with nontrivial topology described by the Hopf map. Despite their recent discovery in colloids and liquid crystals, the requirement of applied fields or confinement for stability impedes their utility in technological applications. Here we demonstrate stable Hopf solitons in a liquid crystal material without these requirements as a result of enhanced stability by tuning anisotropy of parameters that describe energetic costs of different gradient components in the molecular alignment field. Nevertheless, electric fields allow for inter-transformation of Hopf solitons between different geometric embodiments, as well as for their three-dimensional hopping-like dynamics in response to electric pulses. Numerical modelling reproduces both the equilibrium structure and topology-preserving out-of-equilibrium evolution of the soliton during switching and motions. Our findings may enable myriads of solitonic condensed matter phases and active matter systems, as well as their technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Signatures of hot carriers and hot phonons in the re-entrant metallic and semiconducting states of Moiré-gapped graphene

Abstract Stacking of graphene with hexagonal boron nitride (h-BN) can dramatically modify its bands from their usual linear form, opening a series of narrow minigaps that are separated by wider minibands. While the resulting spectrum offers strong potential for use in functional (opto)electronic devices, a proper understanding of the dynamics of hot carriers in these bands is a prerequisite for such applications. In this work, we therefore apply a strategy of rapid electrical pulsing to drive carriers in graphene/h-BN heterostructures deep into the dissipative limit of strong electron-phonon coupling. By using electrical gating to move the chemical potential through the “Moiré bands”, we demonstrate a cyclical evolution between metallic and semiconducting states. This behavior is captured in a self-consistent model of non-equilibrium transport that considers the competition of electrically driven inter-band tunneling and hot-carrier scattering by strongly non-equilibrium phonons. Overall, our results demonstrate how a treatment of the dynamics of both hot carriers and hot phonons is essential to understanding the properties of functional graphene superlattices.

42 ENGINEERING↗
DOE OSTI2021

Adaptive four-level modeling of laser cooling of solids

Laser cooling of rare-earth doped solids has been demonstrated across a wide range of material platforms, inspiring the development of simple phenomenological models such as the four-level model to elucidate the universal properties of laser cooling under various operating conditions. However, these models usually require the input of full absorption spectra that must be provided experimentally or by additional complicated atomic modeling. In this Letter, we propose that a four-level model, when extended to admit effective energy levels adaptive to the pumping photon energy, can accurately predict the cooling efficiency as a function of temperature and pumping frequency using only few inputs such as the absorption coefficient measured at a single frequency and temperature. Our model exploits the quasi-equilibrium properties of the excitation of rare-earth ions for the determination of the effective four energy levels. The model is validated against published experimental results for a span of materials including ytterbium/thulium-doped glass and crystals. With the verified model, we derive explicit expressions for the optimal frequency and the operating bandwidth of pumping lasers. Our model significantly simplifies the modeling process of laser cooling and is expected to stimulate further development of optical refrigeration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Drag exerted by drift-acoustic radiation on a magnetic island

The Island Equilibrium and Transport (IslET) model is extended to allow for the effects of ion compressibility albeit for a constant electron temperature. The model describes fully relaxed, saturated, steady-state magnetic islands arising from the evolution of a tearing mode caused either by an instability or an externally imposed resonant perturbation. It assumes that the collision frequency lies in the semi-collisional regime. A salient new feature of the extended model is an equation describing the transport of parallel momentum. The model is applied to evaluate the drag exerted on a thin island by the emission of drift-acoustic waves.

Physics↗
DOE OSTI2025

A numerically exact description of ultrafast vibrational decoherence in vibration-coupled electron transfer

Broadband pump–probe spectroscopy has been widely used to measure vibrational decoherence associated with the reaction coordinate in photoinduced ultrafast vibration-coupled electron transfer (VCET) reactions. These experiments provide insight into the interplay of intramolecular coordinates along the reaction coordinate. However, a general theoretical foundation for analyzing, and even for explaining rigorously, these data is lacking. In this work, we study vibrational decoherence in a model VCET reaction using the nearly exact time-dependent density matrix renormalization group simulation method. We explore how analyzing the density matrix with quantum information measures can help elucidate the evolution of vibrational coherence in simulations of dynamics. We examine how vibrational coherence is affected by electron transfer on the timescale of approximately 100 femtoseconds. Our results suggest that electron transfer, in the nonadiabatic model, changes the vibrational equilibrium position abruptly—an example of a “quantum quench” event. This explains the concomitant vibrational decoherence. We find that abrupt vibrational decoherence can be mitigated by wavepacket motion occurring on the timescale of the electron transfer.

Science & Technology - Other Topics↗
DOE OSTI2022

Nonequilibrium plasma in water vapor

A kinetic model of a non-equilibrium water vapor plasma created by a nanosecond pulsed discharge at elevated temperatures and at low pressure is presented. The model includes equations for electron, vibrational and gas temperatures combined with equations for neutral, charged, electronically and vibrationally excited species. Simulation results are compared with the available experimental data and show good agreement in dynamics of the plasma decay and electron drift velocities. It is shown that the plasma decay in pure water vapor is controlled by the dissociative recombination. Under the conditions simulated, the concentration of electronically excited species is low, but we observe a formation of vibrationally excited species. Finally, a significant decrease in the electrical conductivity of water plasma with time, which depends on the level of initial ionization, is shown.

42 ENGINEERING↗
DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗
DOE OSTI2021

$\mathcal{g}$ modes of neutron stars with hadron-to-quark crossover transitions

Here we perform the first study of the principal core g-mode oscillation in hybrid stars containing quark matter, utilizing a crossover model for the hadron-to-quark transition inspired by lattice QCD. The ensuing results are compared with our recent findings of g-mode frequencies in hybrid stars with a first-order phase transition using Gibbs constructions. We find that models using Gibbs construction yield g-mode amplitudes and the associated gravitational energy radiated that dominate over those of the chosen crossover model owing to the distinct behaviors of the equilibrium and adiabatic sound speeds in the various models. Based on our results, we conclude that were g modes to be detected in upgraded LIGO and Virgo detectors it would indicate a first-order phase transition akin to a Gibbs construction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accuracy Guarantees and Quantum Advantage in Analog Open Quantum Simulation with and without Noise

Many-body open quantum systems, described by Lindbladian master equations, are a rich class of physical models that display complex equilibrium and out-of-equilibrium phenomena which remain to be understood. In this paper, we theoretically analyze noisy analog quantum simulation of geometrically local open quantum systems and provide evidence that this problem both is hard to simulate on classical computers and could be approximately solved on near-term quantum devices. First, given a noiseless quantum simulator, we show that the dynamics of local observables and the fixed-point expectation values of rapidly mixing local observables in geometrically local Lindbladians can be obtained to a precision of ϵ in time that is poly ( ϵ − 1 ) and uniform in system size. Furthermore, we establish that the quantum simulator would provide a superpolynomial advantage, in run-time scaling with respect to the target precision and either the evolution time (when simulating dynamics) or the Lindbladian’s decay rate (when simulating fixed points), over any classical algorithm for these problems, assuming BQP ≠ BPP . We then consider the presence of noise in the quantum simulator in the form of additional geometrically local Lindbladian terms. We show that the simulation tasks considered in this paper are stable to errors; i.e., they can be solved to a noise-limited, but system-size independent, precision. Finally, we establish that, assuming BQP ≠ BPP , there are stable geometrically local Lindbladian simulation problems such that, as the noise rate on the simulator is reduced, classical algorithms must take time superpolynomially longer in the inverse noise rate to attain the same precision as the analog quantum simulator. Published by the American Physical Society 2025

Kashyap, Vikram (ORCID:0000000208195207)↗
DOE OSTI2022

Computational Fluid Dynamics Modeling of Low Temperature Ignition Processes From a Nanosecond Pulsed Discharge at Quiescent Conditions

Recent interest in nonequilibrium plasma discharges as sources of ignition for the automotive industry has not yet been accompanied by the availability of dedicated models to perform this task in computational fluid dynamics (CFD) engine simulations. The need for a low-temperature plasma (LTP) ignition model has motivated much work in simulating these discharges from first principles. Most ignition models assume that an equilibrium plasma comprises the bulk of discharge kernels. LTP discharges, however, exhibit highly nonequilibrium behavior. In this work, a method to determine a consistent initialization of LTP discharge kernels for use in engine CFD codes like CONVERGE is proposed. The method utilizes first principles discharge simulations. Such an LTP kernel is introduced in a flammable mixture of air and fuel, and the subsequent plasma expansion and ignition simulation is carried out using a reacting flow solver with detailed chemistry. Finally, the proposed numerical approach is shown to produce results that agree with experimental observations regarding the ignitability of methane-air and ethylene-air mixtures by LTP discharges.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Climate Sensitivity and the Direct Effect of Carbon Dioxide in a Limited-Area Cloud-Resolving Model

Even in a small domain, it can be prohibitively expensive to run cloud-resolving greenhouse gas warming experiments due to the long equilibration time. Here, a technique is introduced that reduces the computational cost of these experiments by an order of magnitude: instead of fixing the carbon dioxide concentration and equilibrating the sea surface temperature (SST), this technique fixes the SST and equilibrates the carbon dioxide concentration. Using this approach in a cloud-resolving model of radiative–convective equilibrium (RCE), the equilibrated SST is obtained as a continuous function of carbon dioxide concentrations spanning 1 ppmv to nearly 10 000 ppmv, revealing a dramatic increase in equilibrium climate sensitivity (ECS) at higher temperatures. This increase in ECS is due to both an increase in forcing and a decrease in the feedback parameter. In addition, the technique is used to obtain the direct effects of carbon dioxide (i.e., the rapid adjustments) over a wide range of SSTs. Overall, the direct effect of carbon dioxide offsets a quarter of the increase in precipitation from warming, reduces the shallow cloud fraction by a small amount, and has no impact on convective available potential energy (CAPE).

Romps, David M.↗
DOE OSTI2020

Quantitative defect analysis in MOCVD GaN-on-GaN using cathodoluminescence

Cathodoluminescence (CL) is used as a quantitative characterization technique to probe impurities at the metal-organic chemical vapor deposition (MOCVD) grown GaN-on-GaN homoepitaxial interfaces. CL intensity contrast shows a strong correlation with the interfacial impurity concentrations. Based on the analysis of recombination mechanisms of electron beam induced non-equilibrium carriers, an analytical model is proposed to quantitatively determine the impurity concentrations from CL intensity. The extracted interfacial impurity concentrations from the analytical model show a good agreement with the compensation levels obtained from capacitance-voltage measurement, signifying the potential of CL for probing the quantitative impurity levels in GaN-on-GaN structures. This approach can also be extended to be applied in other material systems.

Chen, Zhaoying↗