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 865 records · Page 48

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

79 ASTRONOMY AND ASTROPHYSICS↗

Status of Multiphysics Modeling of Critical DPCs in a Geological Repository

A summary of previous work focused on the first two scenarios is presented as an established workflow. These canisters could reach an equilibrium at a low quasi–steady-state power level that could be sustained for thousands of years. The majority of this report describes an effort to expand the initial methodology to include the scenario in which a DPC has enough excess reactivity to boil the water inside the canister. This scenario involves a variety of tools and modeling approaches to capture the two-phase behavior, including RELAP5-3D, STAR-CCM+, and TRACE. Efforts using the high-fidelity computational fluid dynamics results to verify the accuracy of RELAP5-3D and TRACE for boiling DPCs are also presented. A new modeling strategy using TRACE appears to be the most promising approach for expanding the analysis of critical DPCs beyond the subcooled scenarios. Finally, this report describes ongoing work investigating the fourth scenario, in which a critical DPC is sealed in saturated bentonite. The primary goal in this research is to understand the power output and associated pressure increase to inform geomechanics simulations that could predict the impact on a geological repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Energy Migration Processes in Re(I) MLCT Complexes Featuring a Chromophoric Ancillary Ligand

We present herein the synthesis, structural characterization, electronic structure calculations, and ultrafast and supra-nanosecond photophysical properties of a series of five Re(I) bichromophores exhibiting metal to ligand charge transfer (MLCT) excited states based on the general formula fac -[Re(N ^ N)(CO) 3 (PNI-py)]PF 6, where PNI-py is 4-piperidinyl-1,8-naphthalimidepyridine and N ^ N is a diimine ligand ( Re1 – 5 ), along with their corresponding model chromophores where 4-ethylpyridine was substituted for PNI-py ( Mod1 – 5 ). The diimine ligands used include 1,10-phenanthroline (phen, 1 ), 2,9-dimethyl-4,7-diphenyl-1,10-phenanthroline (bcp, 2 ), 4,4'-di- tert -butyl-2,2'-bipyridine (dtbb, 3 ), 4,4'-diethyl ester-2,2'-bipyridine (deeb, 4 ), and 2,2'-biquinoline (biq, 5 ). In these metal–organic bichromophores, structural modification of the diimine ligand resulted in substantial changes to the observed energy transfer efficiencies between the two chromophores as a result of the variation in 3 MLCT excited-state energies. The photophysical properties and energetic pathways of the model chromophores were explored in parallel to accurately track the changes that arose from introduction of the organic chromophore pendant on the ancillary ligand. All relevant photophysical and energy transfer processes were probed and characterized using time-resolved photoluminescence spectroscopy, ultrafast and nanosecond transient absorption spectroscopy, and time-dependent density functional theory calculations. Of the five bichromophores in this study, four ( Re1 – 4 ) exhibited a thermal equilibrium between the 3 PNI-py and the 3 MLCT excited state, drastically extending the lifetimes of the parent model chromophores.

14 SOLAR ENERGY↗

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13 C and 2 H Tracers

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with 2 H and 13 C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli . The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellummthan in T. saccharolyticum. These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway. Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae , the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling of plasma beta effects on the island divertor transport in the standard configuration of W7-X

The influence of plasma beta effects on the edge plasma transport in the Wendelstein 7-X standard configuration is studied systematically by using EMC3-EIRENE combined with a 3D equilibrium code named HINT. The magnetic topology changes induced by plasma beta effects are significantly reflected in plasma transport behaviors and heat flux patterns on divertor targets. After validating the modeling strategy by comparisons with experimental data, the extended simulations for high performance plasmas show that the threshold separatrix density for accessing the power detachment is reduced in higher beta plasmas. Compared with the vacuum field case, the impurity radiation distributions with finite beta effects are modified in the magnetic island region. The divertor heat flux is distributed more evenly along the toroidal direction on the strike line at the vertical target. The strike line on the horizontal target moves towards the pumping gap with an increase in the plasma beta. In addition, the different pressure profiles with the same central beta also result in a modified heat flux pattern on the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

EFIT‐AI: Machine Learning and Artificial Intelligence Assisted Equilibrium Reconstruction for Tokamak Experiments and Burning Plasmas (Final Report)

The EFIT-AI project is creating a modern advanced equilibrium reconstruction code suitable for tokamak experiments of burning plasmas. EFIT [1,2] was the first and is the most extensively used equilibrium reconstruction code in the world. This project builds on the production-level experience and adds key elements as follows. 1. A Model Order Reduction (MOR) version of the two-dimensional (2D) Grad-Shafranov equation solver (EFIT-MORNN) using physics-informed neural networks. 2. Improved optimization and data analysis capabilities using a Bayesian framework enhanced with machine learning. 3. A MOR version of the three-dimensional (3D) perturbed equilibrium reconstruction tool.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of plasma squareness on RMP response in MAST-U

The effect of the squareness of the plasma boundary shape on the plasma response to the applied n = 2 (n is the toroidal mode number) resonant magnetic perturbation, and consequently on controlling the edge-localized mode (ELM) in MAST-U in the double-null configuration, is numerically systematically investigated. The modeling is performed utilizing three different plasma response models: the single-fluid model, the two-fluid model and the magnetohydrodynamic-kinetic hybrid model. The effect of the plasma squareness on ELM control is found to depend on the equilbrium kinetic profiles. Fixing the equilibrium radial profiles to that from a recent MAST-U discharge, the experimental plasma shape is found to be close to optimal for ELM control in terms of squareness. The optimal coil phasing does not significantly vary with varying plasma squareness. These findings robustly hold with (i) two figures of merit (FoMs) defined for predicting the ELM control feasibility, and (ii) three plasma response models employed for evaluating these FoMs.

MAST-U↗

Thermodynamic Model and Initial Experimental Investigation of Air Dehumidification through Electrically Charged Vapor Capturing Electrostatic Droplets

The removal of water vapor from the air to reduce relative humidity is a well known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation with the goal to reduce the energy consumption associated with dehumidification. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. In this thesis, a broad theoretical investigation to study the electrostatically enhanced condensation processes was carried out. These processes consist of the use of highly charged particles, preferably highly charged water droplets to attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. The electrical charge promotes the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation on curved surfaces and, consequently, the equilibrium between evaporation and condensation is shifted towards condensation. This investigation resulted in the development of a thermodynamics model, which was able to predict an effective size range of the charged droplets for optimal dehumidification, under ideal conditions. The range resulted in about 2 to 4 m in diameter, while the electrical charge was kept to the maximum limit predicted by the Rayleigh model. A sensitivity analysis on the variation of the size and the charge of the charged droplets was also performed. In terms of dehumidification rates, when six electrospray heads, i.e. the maximum number of heads tested in the experimental campaign, were considered, the model predicted very limited rates, way lower than the target of this research, i.e. a dehumidification rate in terms of relative humidity of 5 % with an air flow rate of 5 cfm. While ways to increase this rate existed, their implementation proved difficult for this initial investigation. According to the requirements predicted by the model, the use of electrosprays appeared the most suitable solution for the production of small but highly charged droplets. The electrospray features and operational modes were studied in detail and an electrospray assembly was designed to be tested in this initial experimental investigation. For this initial investigation, an experimental system was designed and built to validate the thermodynamics model and with the goal to achieve a dehumidification rate of 5 % in terms of relative humidity with an air flow rate of 5 cfm. The experimental system consisted of a wind tunnel test apparatus, dew point sensors, thermocouples for temperature measurements, pressure sensors and other devices. The experimental system was equipped with a computerized data acquisition and storage system. The system was suitable to evaluate the air water content differential before and after the test section where the electrospray heads were installed. The air water content differential was evaluated in terms of the ∆ω/ω_1 ratio, to eliminate the dependence on the dry bulb temperature along the test apparatus. This work presented initial experimental data for electrostatically enhanced dehumidification processes. Working parameters, such as air and water flow rates, high voltage potentials and polarity, deionized water types with different electrical conductivity, were varied to find the best combination for an improved dehumidification. With the set up used and the conditions considered, it was never possible to achieve a 5 % dehumidification rate with 5 cfm. In general, the dehumidification rate was always limited and lower than 1 % for all the air flow rates considered and mainly within the uncertainty of the dew point sensors. These results confirmed the already limited predictions of the thermodynamics model, for ideal conditions and the interesting and promising results obtained in some experimental investigations available in the open literature were not achievable in this thesis. The more probable reasons were the higher air flow rates considered and the complete absence of the use of cooling power, which was mainly implemented to facilitate the vapor condensation.

Morcelli, Stefano↗

Sunny.jl

Sunny is a Julia package designed to serve the needs of the quantum magnetism community. It supports the specification of a very broad class of spin models and a diverse suite of numerical solvers. These include powerful methods for simulating spin dynamics both in and out of equilibrium. Uniquely, it features a broad generalization of classical and semiclassical approaches to SU(N) coherent states, which is useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects. Sunny also offers a well-developed framework for calculating the dynamical spin structure factor, enabling direct comparison with scattering experiments. Ease of use is a priority, with tools for symmetry-guided modeling and interactive visualization.

Dahlbom, David↗

Sunny.jl: A Julia Package for Spin Dynamics

Sunny is a Julia package designed to serve the needs of the quantum magnetism community. It supports the specification of a very broad class of spin models and a diverse suite of numerical solvers. These include powerful methods for simulating spin dynamics both in and out of equilibrium. Uniquely, it features a broad generalization of classical and semiclassical approaches to SU(N) coherent states, which is useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects. Sunny also offers a well-developed framework for calculating the dynamical spin structure factor, enabling direct comparison with scattering experiments. Ease of use is a priority, with tools for symmetry-guided modeling and interactive visualization.

97 MATHEMATICS AND COMPUTING↗

A Machine Learning Assisted Development of a Model for the Populations of Convective and Stratiform Clouds

Abstract Traditional parameterizations of the interaction between convection and the environment have relied on an assumption that the slowly varying large‐scale environment is in statistical equilibrium with a large number of small and short‐lived convective clouds. They fail to capture nonequilibrium transitions such as the diurnal cycle and the formation of mesoscale convective systems as well as observed precipitation statistics and extremes. Informed by analysis of radar observations, cloud‐permitting model simulation, theory, and machine learning, this work presents a new stochastic cloud population dynamics model for characterizing the interactions between convective and stratiform clouds, with the goal of informing the representation of these interactions in global climate models. Fifteen wet seasons of precipitating cloud observations by a C‐band radar at Darwin, Australia are fed into a machine learning algorithm to obtain transition functions that close a set of coupled equations relating large‐scale forcing, mass flux, the convective cell size distribution, and the stratiform area. Under realistic large‐scale forcing, the derived transition functions show that, on the one hand, interactions with stratiform clouds act to dampen the variability in the size and number of convective cells and therefore in the convective mass flux. On the other, for a given convective area fraction, a larger number of smaller cells is more favorable for the growth of stratiform area than a smaller number of larger cells. The combination of these two factors gives rise to solutions with a few convective cells embedded in a large stratiform area, reminiscent of mesoscale convective systems.

54 ENVIRONMENTAL SCIENCES↗

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗

BBN photodisintegration constraints on gravitationally produced vector bosons

Gravitational production of massive particles due to cosmic expansion can be significant during the inflationary and reheating period of the Universe. If the particle also has non-gravitational interactions that do not significantly affect its production, numerous observational probes open up, including cosmological probes. In this work, we focus on the gravitational production of light vector bosons that couple feebly to the Standard Model (SM) particles. Due to the very feeble coupling, the light vector bosons never reach thermal equilibrium, and if the Hubble scale at the end of inflation is above 10 8 GeV, the gravitational production can overwhelm the thermal production via the freeze-in mechanism by many orders of magnitude. As a result, much stronger constraints from the Big Bang Nucleosynthesis (BBN) can be placed on the lifetime and mass of the vector bosons compared to the scenario where only thermal production is considered. As an example, we study the sub-GeV scale dark photons, which couple to the SM only through kinetic mixing, and derive constraints on the mass and kinetic mixing parameter of the dark photon from the photodisintegration effects on the light element abundances relevant at the end of the BBN when the cosmic age was around 10 4 s.

79 ASTRONOMY AND ASTROPHYSICS↗

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗