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 1,135 records · Page 63

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of Electro-static Assisted Air Dehumidification

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. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. Also, the findings from the literature focused mainly on small flow rates on the order of one to two cfm. Electrically-enhanced condensation consists of the use of highly charged particles, preferably highly charged water droplets, that attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. An effect of the electric charge is the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation. Consequently, the equilibrium between evaporation and condensation is shifted towards condensation. Following the application of the modified Kelvin-Thomson theory, we developed a preliminary physics-based model to predict an effective size range of the charged droplets for optimal dehumidification. The range resulted in about 2 to 4 μm in diameter, under few simplifying assumptions. The effect of the size and the charge of the electrosprayed droplets on the overall dehumidification rate was briefly discussed. The use of electrosprays to produce small but highly charged droplets was broadly discussed in this paper. The cone-jet mode was identified as the most suitable electrospray operational mode, and it generated droplets of small size and high electrical charge. The cone-jet stability was also analyzed in detail. The preliminary data of the present work and the model results indicated that several electrospray heads were required to achieve a 5% dehumidification rate for airflow rates of about 5 cfm.

Morcelli, Stefano↗

Analysis of New Data of Electro-static Assisted Air Dehumidification Processes

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. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. Also, the findings from the literature focused mainly on small flow rates on the order of one to two cfm. Electrically-enhanced condensation consists of the use of highly charged particles, preferably highly charged water droplets, that attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. An effect of the electric charge is the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation. Consequently, the equilibrium between evaporation and condensation is shifted towards condensation. Following the application of the modified Kelvin-Thomson theory, we developed a preliminary physics-based model to predict an effective size range of the charged droplets for optimal dehumidification. The range resulted in about 2 to 4 m in diameter, under few simplifying assumptions. The effect of the size and the charge of the electrosprayed droplets on the overall dehumidification rate was briefly discussed. The use of electrosprays to produce small but highly charged droplets was broadly discussed in this paper. The cone-jet mode was identified as the most suitable electrospray operational mode, and it generated droplets of small size and high electrical charge. The cone-jet stability was also analyzed in detail. The preliminary data of the present work and the model results indicated that several electrospray heads were required to achieve a 5% dehumidification rate for airflow rates of about 5 cfm.

Morcelli, Stefano↗

Development of a Representative Molten Chloride Fast Reactor Model to Assess the Impact of Nuclear Data

The SCALE code system was employed to conduct a preliminary investigation of nuclear data impacts for a fast spectrum molten chloride salt reactor. A computationally effective depletion model that is representative of the reactor system was successfully developed and used to conduct fuel depletion simulations. Development of this model draws from the SLICE method that was developed at Oak Ridge National Laboratory to enable generation of fuel compositions for an advanced reactor core at equilibrium operation. Eigenvalue uncertainty calculations using the ENDF/B-VII.1 nuclear data library were performed for the reactor in the fresh fuel state and an irradiated fuel state. It was determined that the primary driver of eigenvalue uncertainty was the uncertainty in the 235U (n, 𝛾) cross section. Uncertainty calculation results from this study were compared to results available for a different fast system, a sodium-cooled fast reactor, to confirm similarities and identify differences with respect to nuclear data impacts between the two fast advanced reactor systems.

Hirji, Rakim [Georgia Institute of Technology]↗

Projections of H-mode access and edge pedestal in the SPARC tokamak

In order to inform core performance projections and divertor design, the baseline SPARC tokamak plasma discharge is evaluated for its expected H-mode access, pedestal pressure and edge-localized mode (ELM) characteristics. A clear window for H-mode access is predicted for full field DT plasmas, with the available 25 MW of design auxiliary power. Additional alpha heating is likely needed for H-mode sustainment. Pressure pedestal predictions in the developed H-mode are surveyed using the EPED model. The projected SPARC pedestal would be limited dominantly by peeling modes and may achieve pressures in excess of 0.3 MPa at a density of approximately 3 × 10 20 m -3 . High pedestal pressure is partially enabled by strong equilibrium shaping, which has been increased as part of recent design iterations. Edge-localized modes (ELMs) with >1 MJ of energy are projected, and approaches for reducing the ELM size, and thus the peak energy fluence to divertor surfaces, are under consideration. The high pedestal predicted for SPARC provides ample margin to satisfy its high fusion gain (Q) mission, so that even if ELM mitigation techniques result in a 2x reduction of the pedestal pressure, Q > 2 is still predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory

Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assembly of polyelectrolyte star block copolymers at the oil–water interface

To understand and resolve adsorption, reconfiguration, and equilibrium conformations of charged star copolymers, we carried out an integrated experimental and coarse-grained molecular dynamics simulation study of the assembly process at the oil–water interface. This is important to guide development of novel surfactants or amphiphiles for chemical transformations and separations. The star block copolymer consisted of arms that are comprised of hydrophilic–hydrophobic block copolymers that are covalently tethered via the hydrophobic blocks to one point. The hydrophobic core represents polystyrene (PS) chains, while the hydrophilic corona represents quaternized poly(2-vinylpyridine) (P2VP) chains. The P2VP is modeled to become protonated when in contact with an acidic aqueous phase, thereby massively increasing the hydrophilicity of this block, and changing the nature of the star at the oil–water interface. This results in a configurational change whereby the chains comprising the hydrophilic corona are significantly stretched into the aqueous phase, while the hydrophobic core remains solubilized in the oil phase. In the simulations, we followed the kinetics of the anchoring and assembly of the star block copolymer at the interface, monitoring the lateral assembly, and the subsequent reconfiguration of the star via changes in the interfacial tension that varies as the degree-of-protonation increases. At low fractions of protonation, the arm cannot fully partition into the aqueous side of the interface and instead interacts with other arms in the oil phase forming a network near the interface. These insights were used to interpret the non-monotonic dependence of pH with the asymptotic interfacial tension from pendant drop tensiometry experiments and spectral signatures of aromatic stretches seen in vibrational sum frequency generation (SFG) spectroscopy. We describe the relationship of interfacial tension to the star assembly via the Frumkin isotherm, which phenomenologically describes anti-cooperativity in adsorbing stars to the interface due to crowding. Although our model explicitly considers long-range electrostatics, the contribution of electrostatics to interfacial tension is small and brought about by strong counterion condensation at the interface. Finally, these results provide key insights into resolving the adsorption, reconfiguration, and equilibrium conformations of charged star block copolymers as surfactants.

36 MATERIALS SCIENCE↗

Phase-field simulations of fission gas bubbles in high burnup UO2 during steady-state and LOCA transient conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantifying metal-binding specificity of Cc NikZ-II from Clostridium carboxidivorans in the presence of competing metal ions

Many proteins bind transition metal ions as cofactors to carry out their biological functions. Despite binding affinities for divalent transition metal ions being predominantly dictated by the Irving-Williams series for wild-type proteins, in vivo metal ion binding specificity is ensured by intracellular mechanisms that regulate free metal ion concentrations. However, a growing area of biotechnology research considers the use of metal-binding proteins in vitro to purify specific metal ions from wastewater, where specificity is dictated by the protein's metal binding affinities. A goal of metalloprotein engineering is to modulate these affinities to improve a protein's specificity towards a particular metal; however, the quantitative relationship between the affinities and the equilibrium metal-bound protein fractions depends on the underlying binding mechanisms. Here we demonstrate a high-throughput intrinsic tryptophan fluorescence quenching method to validate binding models in multi-metal solutions for CcNikZ-II, a nickel-binding protein from Clostridium carboxidivorans. Using our validated models, we quantify the relationship between binding affinity and specificity in different classes of metal-binding models for CcNikZ-II. In conclusion, we further illustrate the potential relevance of data-informed models to predicting engineering targets for improved specificity.

59 BASIC BIOLOGICAL SCIENCES↗

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Harnessing dimethyl ether with ultra-low-grade heat for scaling-resistant brine concentration and fractional crystallization

Solvent-driven separations may enable scalable concentration of hypersaline brines, supporting a circular resource economy from the extraction of lithium and rare earth elements from spent battery and magnet leachates. This work analyses a novel solvent-driven water extraction (SDWE) system employing dimethyl ether (DME) and ultra-low-grade heat for brine concentration and fractional crystallization. SDWE exploits DME’s unique properties: (1) a low dielectric constant that promotes water solubility over charged solutes by a factor of 10 3 , and (2) a high volatility that facilitate efficient DME reconcentration with ultra-low-grade heat. The techno-economic viability of SDWE is assessed with a computational framework that encompasses a liquid–liquid separator and a solvent concentrator. We integrate the extended universal quasichemical model with the virial equation of state to predict the compositions of the complex three-phase DME-water mixture at vapor–liquid and liquid–liquid equilibrium. Subsequently, we optimize the thermodynamic and economic performance of SDWE, by controlling the interstage flash pressure, heat source temperature, and the number of concentrating stages. DME-based SDWE concentrates an input saline feed to 5.5 M and regenerates over 99 % of the DME using ultra-low-grade heat below 50 °C, with a DME/water selectivity ratio of 125. Here our calculations reveal that optimal performance is achieved at interstage flash pressures of 0.4 – 0.5 bar for heat source temperatures between 323–373 K, with improved exergetic efficiencies at lower temperatures. At a heat source temperature of 323 K and an interstage pressure of 0.489 bar, DME-driven SDWE achieves an optimal thermodynamic efficiency of 20.5 % and a projected specific cost of US$ 1.93 m -3 . These specific costs suggest that SDWE is competitive with commercialized thermal distillation technologies, while mitigating the traditional risks associated with scaling in heat and mass exchangers with hypersaline brines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural Heterogeneity and Hydrodynamics of an Intrinsically Disordered Protein Condensate

Biology demonstrates precise control over the free-energy landscape through the selective partitioning of biomacromolecules into membraneless organelles, enabling essential functions such as biochemical transformations, signaling cascades, and mechanical reinforcement. Although the function of these condensates depends on their underlying structure and hydrodynamics, molecular-scale information on these systems remains sparse. Here, in this study, neutron scattering is used to probe the organization and dynamics of the intrinsically disordered N-terminal domain of Galectin-3, an extracellular lectin responsible for facilitating liquid–liquid phase separation on the cellular surface, in both dilute and condensed phases. Dilute solutions contain isolated protein chains in equilibrium with mesoscopic clusters, whereas the condensed phase adopts a bicontinuous, microemulsion-like morphology. The dilute phase behavior is quantitatively described by coarse-grained polymer models from soft-matter physics, demonstrating their predictive power for complex biological proteins. At elevated concentrations, the proteins self-assemble akin to block copolymers, microphase separating through the aggregation of hydrophobic domains along the protein contour. The resulting condensate remains fluid-like despite a 25-fold increase in concentration; its internal hydrodynamics slow by only a factor of 3 relative to dilute protein chains. These results provide a molecular-level framework for how disordered proteins achieve both the structural complexity and dynamic fluidity of biomolecular condensates.

Carrick, Brian R. [Massachusetts Inst. of Technolo↗

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of V o complex

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a Vo proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo–electron microscopy (cryo-EM) structure of yeast V o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V o that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

59 BASIC BIOLOGICAL SCIENCES↗

Framework For Performing Time-dependent Pebble Bed Reactor Simulations

The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. The code kugelpy is within the GitHub repository `pyrates`.

Stewart, RyanH. [Idaho National Laboratory (INL), ↗

Theoretical Studies of Polar Systems near Ferroelectrics Quantum Critical Points (Final Report)

The PI has studied emergent quantum phases near polar quantum critical points. A first challenge was to develop a theoretical description of the observation that many polar materials undergo classical first-order transitions while displaying quantum criticality. The possibility of novel metallic states near polar quantum critical points was then studied; identification and characterization of a non-Fermi liquid phase was made in a multiband system. In dilute quantum critical polar metals, electronic coupling to polar energy fluctuations was shown to result in attractive electron-electron interactions and to superconductivity. A spin-phonon resonance measurement in applied magnetic field was proposed to determine the magnitude of spin-orbit mediated electron-phonon coupling; furthermore in the polar phase new field-induced phonon collective modes were identified with specific signatures for experiment. In order to keep current, the PI has also worked on other types of strongly correlated problems. With experimental groups, she has developed phenomenologies to identify and explain observed dynamical critical behaviors. The PI has studied the interplay of topology, fractionalization and deconfinement in a minimalist spin model. An emergent phase transition in a 1 + 1 frustrated spin nanotube has also been identified and characterized. The PI has also begun exploring criticality out of equilibrium, characterizing rich dynamical phases associated with photoinduced polar transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗