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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Confusion-Driven Machine Learning of Structural Phases of a Flexible, Magnetic Stockmayer Polymer

We use a semisupervised, neural-network-based machine learning technique, the confusion method, to investigate structural transitions in magnetic polymers, which we model as chains of magnetic colloidal nanoparticles characterized by dipole–dipole and Lennard-Jones interactions. As input for the neural network, we use the particle positions and magnetic dipole moments of equilibrium polymer configurations, which we generate via replica-exchange Wang–Landau simulations. We demonstrate that by measuring the classification accuracy of neural networks, we can effectively identify transition points between multiple structural phases without any prior knowledge of their existence or location. We corroborate our findings by investigating relevant conventional order parameters. Our study furthermore examines previously unexplored low-temperature regions of the phase diagram, where we find new structural transitions between highly ordered helicoidal polymer configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of a Systematic and Extensible Force Field for Peptoids (STEPs)

Peptoids (N-substituted glycines) are a class of biomimetic polymers that have attracted significant attention due to their accessible synthesis and enzymatic and thermal stability relative to their naturally occurring counterparts (polypeptides). While these polymers provide the promise of more robust functional materials via hierarchical approaches, they present a new challenge for computational structure prediction for material design. The reliability of calculations hinges on the accuracy of interactions represented in the force field used to model peptoids. For proteins, structure prediction based on sequence and de novo design has made dramatic progress in recent years; however, these models are not readily transferable for peptoids. Current efforts to develop and implement peptoid-specific force fields are spread out, leading to replicated efforts and a fragmented collection of parameterized sidechains. Here, we developed a peptoid-specific force field containing 70 different side chains, using GAFF2 as starting point. The new model is validated based on the generation of Ramachandran-like plots from DFT optimization compared against force field reproduced potential energy and free energy surfaces as well as the reproduction of equilibrium cis/trans values for some residues experimentally known to form helical structures. In conclusion, equilibrium cis/trans distributions (Kct) are estimated for all parameterized residues to identify which residues have an intrinsic propensity for cis or trans states in the monomeric state.

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

Comment on “Microhydration of Biomolecules: Revealing the Native Structures by Cold Ion IR Spectroscopy”

This paper presents a re-examination of the conclusions of a study reported in The Journal of Physical Chemistry Letters (Saparbaev et al., 2021, 12, 907) that compared the structure of microsolvated ions formed by electrospray ionization to those formed in the gas-phase via a previously published cryogenic ion trap approach. We conducted additional experiments that clearly show that most of the observed differences in the IR spectra can be accounted for by considering the different spectroscopic action schemes used to obtain them. In particular, the presence of the D 2 -tag induces shifts in some of the N-H and O-H peaks which need to be carefully considered before comparing spectra. Furthermore, once these spectral effects are taken into account, we show that both clustering approaches yields similar clusters structures for the small GlyH + (H 2 O)n species. Using unimolecular reaction rate theory, we also show that for the small complexes considered here, only the gas-phase equilibrium distribution of conformers should be expected in both experimental approaches. In addition, the barrier heights necessary to kinetically trap high-energy conformers at 298K is explored using a series of model polyalanine chains.

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

Kibble–Zurek mechanism of Ising domains

The formation of topological defects after a symmetry-breaking phase transition is an overarching phenomenon that encodes the underlying dynamics. The Kibble-Zurek mechanism (KZM) describes these non-equilibrium dynamics of second-order phase transitions and predicts a power-law relationship between the cooling rates and the density of topological defects. It has been verified as a successful model in a wide variety of physical systems, including structure formation in the early Universe and condensed-matter materials. However, it is uncertain if the KZM mechanism is also valid for topologically trivial Ising domains, one of the most common and fundamental types of domain in condensed-matter systems. Here we show that the cooling rate dependence of Ising domain density follows the KZM power law in two different three-dimensional structural Ising domains: ferro-rotation domains in NiTiO 3 and polar domains in BiTeI. However, although the KZM slope of NiTiO 3 agrees with the prediction of the 3D Ising model, the KZM slope of BiTeI exceeds the theoretical limit, providing an example of steepening KZM slope with long-range dipolar interactions. Finally, our results demonstrate the validity of KZM for Ising domains and reveal an enhancement of the power-law exponent for transitions of non-topological quantities with long-range interactions. The Kibble-Zurek mechanism is shown to apply to structural Ising domains in three-dimensional materials. Long-range interactions modify the critical exponents away from theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Air Carbon for Algae Production (AirCAP) – Expanding algae resource potential via direct (in-pond) air-CO2 capture.

Microalgae biofuel production potential in the United States is constrained by a lack of suitable cultivation locations within a close proximity to CO2 point-sources. Initial estimates suggest that removing the requirement of a concentrated CO2 source expands resource potential ten-fold, allowing the projection of production levels that approach levels of current U.S. diesel consumption. Decoupling cultivation locations from concentrated CO2 point-sources requires direct capture of air-CO2 into the pond bulk at a rate sufficient to meet maximum expected carbon assimilation rates (10 g C/m2-day and higher). Enhancement of the ingassing rate via chemical reaction, the direct reaction of dissolved CO2 with hydroxide ions, has been proposed as a mechanism to increase the ingassing rate of air-CO2. The central objectives of this project are to validate whether chemically enhanced air-CO2 ingassing rates in large raceway ponds can meet or exceed the level required to support high levels of algae productivity, then subsequently identify alkaliphilic strains able to withstand such conditions while maintaining a biomass productivity that enables process economics. Full scale ponds (~1 acre, unlined) will be used to validate ingassing rates as a function of pH expected at scale in a series of abiotic experiments. To date, mass-transfer dynamics have been characterized in ponds ranging from 1 to 1,000 m2. Initial results will be presented describing the mass transfer coefficient in 1-acre (~4,000 m2), unlined ponds. The maximum expected carbon ingassing rate will be measured during ingassing trials, in which pond pH is displaced from equilibrium via the addition of a strong base, with the subsequent decrease in pH (and increase in dissolved inorganic carbon) corresponding to the air-CO2 absorption rate. A mass-transfer model, validated with experimental results, will be used to project ingassing rates expected at pond scales (10 acres) required for production of commodity chemicals.

09 BIOMASS FUELS↗
DOE OSTI2022

Measurements of Electron Beam-Induced Spin-Relaxation in Frozen-Spin Hydrogen-Deuteride (HD)

Targets of solid hydrogen-deuteride (HD) can be prepared in a frozen-spin state with spin relaxation times (T1) in excess of a year. The present work has studied their potential use in experiments with electron beams. Polarization relaxation rates have been measured during exposure to sub-nanoAmp currents of 9.7 MeV electron beams from the recently commissioned Upgraded Injector Test Facility (UITF) at Jefferson Lab (JLab). These UITF measurements can be used to anticipate the expected performance at the GeV energies typical of JLab experiments, since the energy deposition in a target is almost independent of beam energy. An in-beam dilution refrigerator equipped with superconducting solenoids has been used to maintain solid HD samples at ?0.1 K and ?1 tesla, and internal NMR coils have been used to monitor hydrogen polarization. Spin-relaxation in two frozen-spin targets, with initial H-polarizations of 40% and 34%, have been tracked while exposed to beams under varying conditions of current, dose, beam duty factor and temperature. At a fixed accumulated dose, the spin-relaxation rates drop with current, suggesting depolarization by the charge cloud of the beam. Imposing a duty factor on the otherwise continuous UITF beam, with forced millisecond scale gaps, has shown no obvious correlation to polarization loss rates. After an accumulated dose of ?2 ?C/cm2, beam-off spin-relaxation rates drop from their immeasurably long pre-irradiation values to the order of weeks (with beam-on T1 values typically an order of magnitude shorter), reflecting a buildup of paramagnetic charge centers within the HD lattice. The accumulated polarization loss was approximately proportional to dose in both targets, dropping to 1/e of their initial values after ?6 ?C/cm2. Thermal equilibrium polarizations of targets not in the frozen-spin state (with intentionally short T1) have been used to deduce the in situ temperature of solid HD while under electron bombardment. A model for depolarization by beam-associated paramagnetic impurities largely accounts for the data, and suggests that improvements in heat removal could lead to significant increases in the in-beam T1.

O'Connell, Thomas↗
DOE OSTI2022

Optimization Modeling for Advanced Syngas to Olefin Reactive Systems

Reactor designs with mixed catalysts play an important role in transforming a multiple reactor system to single-shot reactors. In addition to savings in capital and ease of implementation, single-shot reactors are useful to break equilibrium limitations, thereby increasing the yield and selectivity of desired product as shown in previous studies. However, the nonlinear and highly exothermic nature of mixed-catalyst systems makes it difficult for commercial process simulation and optimization tools to optimize these systems. This study describes the development and application of optimization strategies for mixed-catalyst, single-shot reactors for syngas to olefin (STO) processes. Finding the optimal catalyst distribution is challenging and requires advanced solution strategies for singular optimal control problems, which are poorly conditioned and often lead to flat response surfaces. The graded bed and partial-moving finite-element approaches are used to find the optimal catalyst distribution that maximizes the olefins yield. A 1.3% increase in the yield is observed from one zone to three zones. The yield further improves from three zones to the exact infinite dimensional solution by 0.2%. This improvement can be realized in practice by changing only the catalyst distribution, without any extra investment. Lastly, the results suggest that a mixed-catalyst single shot reactor bed can be applied to other reaction mechanisms to increase reactor performance.

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

Nanoscale mapping of point defect concentrations with $\mathrm{4D-STEM}$

Vacancies are missing atoms in a crystalline material, and occur both at equilibrium (varying with temperature) and out of equilibrium such as when crystalline materials are damaged with radiation or corrosion. While we know of their importance, particularly regarding diffusive mechanisms, it is not straightforward to experimentally measure their concentration or to visualize them directly. Traditionally, measurements such as positron annihilation spectroscopy and to some degree X-ray diffraction can measure average concentrations, but generally lack the ability to visualize or quantify a heterogenous concentration of vacancies that can occur at the level of individual defects and microstructural features. Here, we present a method to map vacancy concentrations and their distribution using local lattice parameter measurements with a high-resolution electron microscope. Our method utilizes a Au thin film as a model to demonstrate the method via four-dimensional scanning transmission electron microscopy (4D-STEM) by correlating the differences between changes in lattice parameter and the volumetric thermal expansion during in situ heating experiments. The vacancy mapping methodology is also applied to non-equilibrium defects accumulated in pure Al via knock-on electron beam irradiation. Our method demonstrates the ability to map point defect concentrations in heterogeneous systems in situ with nanometer spatial resolution. The result is a technique that can provide direct measurements of vacancy concentrations at the level of individual defects in studies of materials in and out of equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Ion–ion interactions enhance aluminum solubility in alkaline suspensions of nano-gibbsite (α-Al(OH) 3 ) with sodium nitrite/nitrate

Despite widespread industrial importance, predicting metal solubilities in highly concentrated, multicomponent aqueous solutions is difficult due to poorly understood ion–ion and ion–solvent interactions. Aluminum hydroxide solid phase solubility in concentrated sodium hydroxide (NaOH) solutions is one such case, with major implications for ore refining, as well as processing of radioactive waste stored at U.S. Department of Energy legacy sites, such as the Hanford Site, Washington State. The solubility of gibbsite (α-Al(OH) 3 ) is often not well predicted because other ions affect the activity of hydroxide (OH – ) and aluminate (Al(OH) 4 – ) anions. In the study reported herein, we systematically examined the influence of key anions, nitrite (NO 2 – ) and nitrate (NO 3 – ), as sodium salts on the solubility of α-Al(OH) 3 in NaOH solutions taking care to establish equilibrium from both under- and oversaturation. Rapid equilibration was enabled by use of a highly pure and crystalline synthetic nano-gibbsite of well-defined particle size and shape. Measured dissolved aluminum concentrations were compared with those predicted by an α-Al(OH) 3 solubility model derived for simple Al(OH) 4 – /OH – systems. Specific anion effects were expressed as an enhancement factor (Al enhc ) conveying the excess of dissolved aluminum. At 45 °C, NaNO 2 and NaNO 3 -containing systems exhibited Al enhc values of 2.70 and 1.88, respectively, indicating significant enhancement. The solutions were examined by Raman and high-field 27 Al NMR spectroscopy, indicating specific interactions including Al(OH) 4 – –Na + contact ion pairing and Al(OH) 4 – –NO 2 – /NO 3 – ion–ion interactions. Dynamic evolution of the α-Al(OH) 3 particles including growth and agglomeration was observed revealing the importance of dissolution/reprecipitation in establishing equilibrium. These studies indicate that incomplete ion hydration, as a result of the low water activity in these concentrated electrolytes, results in: (i) enhanced reactivity of the hydroxide ion with respect to α-Al(OH) 3 ; (ii) increased concentrations of Al(OH) 4 – in solution; and (iii) stronger ion–ion interactions that act to stabilize the supersaturated solutions. This information on the mechanisms by which α-Al(OH) 3 becomes supersaturated is essential for more energy-efficient aluminum processing technologies, including the treatment of millions of gallons of Al(OH) 4 – -rich high-level radioactive waste.

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

Positional ordering induced by dynamic steric interactions in superparamagnetic rods

The dynamic motion produced by precessing magnetic fields can drive matter into far-from-equilibrium states. We predict 1D periodic ordering in systems of precessing rods when magnetic interactions between rods remain insignificant. The precession angle of the rods is completely determined by the field's precession angle and the ratio of the field's precession frequency and the characteristic response frequency of the rods. Here, we develop a molecular dynamics model that explicitly calculates magnetic interactions between particles, and we also simulate rods in the limit of a strong and fast precessing magnetic field where inter-rod magnetic interactions are negligible, using a purely steric model. Our simulations show how steric interactions drive the rods from a positionally disordered phase (nematic) to a layered (smectic) phase. As the rod precession angle increases, the nematic–smectic transition density significantly decreases. The minimization of unfavorable steric interactions also induces phase separation in binary mixtures of rods of different lengths. This effect is general to any force that produces precession in elongated particles. This work will advance the understanding and control of out-of-equilibrium soft matter systems.

42 ENGINEERING↗
DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Consequences of Lattice Mismatch for Phase Equilibrium in Heterostructured Solids

Lattice mismatch - the difference in natural bond length between two species comprising a solid is a common occurrence in modern materials. Spatial patterns in certain nanocrystal heterostructures, for instance, are thought to be a result of elastic forces arising from lattice mismatch. However, the role of lattice mismatch in mediating the arrangement of atoms in such materials is incompletely understood. Here we consider a simple microscopic model for lattice mismatch, in which the difference in natural lengths between bonded atoms produces interactions mediated by elastic strain. Computer simulations reveal that the model exhibits rich phase behavior, including structures with periodic order and unusual coexistence scenarios. To explain this phase behavior, we derive an effective pairwise interaction potential between atoms, revealing preferred spatial arrangements of atoms driven by spatial variations in the interaction potential. We then develop a mean field theory based on this effective interaction which qualitatively captures observed transition to phases with modulated order. Finally, we explain the observed scenarios of coexistence between using a graphical construction, based on the realization that the free energy cost of phase separation in elastic systems grows with system size. Here, these results clarify the equilibrium effects of lattice mismatch in macroscopic solids and suggest a role for lattice mismatch in creating spatially heterogeneous compositions in nanoscale materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

DNA Reaction–Diffusion Attractor Patterns

Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction-diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and ‘hill’-shaped patterns formed and stabilized into shapes and at time scales consistent with reaction-diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady state profiles. Recovery also occurred after repeated perturbations. As a result, by designing the far-from-equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration

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

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗