Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive simulation”

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 145 records · Page 8

Foundation models for atomistic simulation of chemistry and materials

Conventional computational methods for modeling chemical and materials systems are limited by system size and timescale, forcing a trade-off between quantum-mechanical accuracy and the sampling needed for realistic observables. Large language and vision foundation models — pre-trained on massive datasets using transformer architectures — have revolutionized many fields. It is thus interesting to ask whether a foundation model — subject to suitable data, parameter scaling and training — could enable learned simulations of chemistry and materials. Here, in this study, we review the field of machine-learned interatomic potentials (MLIPs) and posit that scaling up large and diverse chemical and materials datasets and highly expressive architectures using advanced training strategies should result in models that are: more efficient, transferable, robust to out-of-distribution scenarios, and easier to fine-tune to a variety of downstream physical observables than models trained from scratch on small datasets corresponding to specific, targeted atomistic simulation tasks. We provide specific criteria for creating such large-scale MLIP foundation models, coordinated strategies for their development, evaluation and deployment, and highlight potential emergent capabilities that could transform predictive simulations in chemistry and materials science and accelerate discovery across multiple technological domains.

Yuan, Eric C.-Y. [University of California, Berkel↗

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure and Stability of Iron Fluoride at High Pressure–Temperature and Implication for a New Reservoir of Fluorine in the Deep Earth

Fluorine (F) is the most abundant halogen in the bulk silicate Earth. F plays an important role in geochemical and biological systems, but its abundance and distribution in the terrestrial mantle are still unclear. Recent studies suggested that F reservoirs in the deep mantle are potentially hosted in terrestrial oxide minerals, especially in aluminous bridgmanite. However, the knowledge about the formation and stability field of fluoride in the Earth’s interior is rare. In this study, we combine in situ laser-heated diamond anvil cell, synchrotron X-ray diffraction, and first-principles structure search to show that a new tetragonal structure of FeF 3 is stable at pressures of 78–130 GPa and temperatures up to ~1900 K. Simulation predicted the tetragonal phase takes a much denser structure due to the rotation of FeF 6 octahedral units. The equations of states of tetragonal FeF3 are determined by experiment and verified by simulation. Our results indicate that FeF 3 can be a potential key phase for storing F in the Earth’s lower mantle and may explain some mantle-derived magma with high F concentration.

58 GEOSCIENCES↗

Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions

The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.

43 PARTICLE ACCELERATORS↗

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs’ advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM + ][TFSI – ]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ~1 solvent molecule. Finally, the scarcity of solvent–ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicted performance of a tangential viewing hard x-ray camera for the DIII-D high field side lower hybrid current drive experiment

High field side launch of lower hybrid current drive (LHCD) has improved accessibility and penetration over low field side launch on DIII-D. Simulations predict single pass absorption under a wide range of plasma conditions. Hard x-ray (HXR) measurement of LHCD generated fast electron bremsstrahlung (50–250 keV) will validate wave propagation and absorption. Emissivity profiles are recovered from one-dimensional inversion of HXR brightness to determine LH damping location, fast electron slowing down time, and some indication of the fast electron energy. The camera will be implemented by populating 32 tangential sightlines of the existing Gamma Ray Imager with Kromek SPEARTM Cadmium Zinc Telluride (CZT) detectors sensitive to 10–1000 keV photons with 10 keV energy resolution. Expected count rates allow for <0.5 ms time resolution. Pulses are processed using 50 ns shaping time Cremat CR-200 Gaussian shaping modules and are digitized by 25 MHz D-TACQ ACQ216 digitizers. The performance of the HXR camera is evaluated by comparing predicted fast electron density profiles and inverted synthetic brightnesses obtained from the ray-tracing/Fokker–Planck codes GENRAY/CQL3D. Inversions closely matched predicted fast electron profiles for a range of experimental parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗

Theoretically Informed Kinetics (ThInK): Establishing a modern C 0 -C 3 mechanism for combustion modeling

In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.

Energy transfer↗

Numerical assessment of the impact of the guiding-centre approximation on fast ion simulations in NSTX

Guiding-centre (GC) and full-orbit (FO) simulations of the beam-injected fast ion distribution and the corresponding neutron emissivity have been carried out for magnetohydrodynamics-quiescent National Spherical Torus eXperiment (NSTX) plasmas, using a combination of ASCOT5 and DRESS, to assess the suitability of the GC approximation for fast ions in NSTX. It was found that GC and FO simulations predicted substantially different steady-state distributions in both position and velocity space and different neutron emissivity profiles, leading to a 15% reduction in the predicted global neutron rate for FO relative to GC. These changes accompany a higher magnetic moment in FO, and correspond to a change in particle orbits from co-passing to trapped and stagnation orbits. ASCOT5 was also benchmarked against TRANSP/NUBEAM with input loaded entirely from TRANSP/NUBEAM output files, with agreement found between the GC simulations when finite Larmor radius (FLR) corrections were omitted. ASCOT5 FO and TRANSP/NUBEAM with FLR produced fast ion distributions which differed in localised regions, but predicted global neutron rates which agree within 3%.

ASCOT↗

Root Cause Analysis of Degradation in Protonic Ceramic Electrochemical Cell with Interfacial Electrical Sensors Using Data‐Driven Machine Learning

Abstract Protonic ceramic electrochemical cells (PCECs) offer promising paths for energy storage and conversion. Despite considerable achievements made, PCECs still face challenges such as physiochemical compatibility between componenets and suboptimal solid–solid contact at the interfaces between the electrolytes and electrodes. In this study, a novel approach is proposed that combines in situ electrochemical characterization of interfacial electrical sensor embedded PCECs and machine learning to quantify the contributions of different cell components to total degradation, as well as to predict the remaining useful life. The experimental results suggest that the overpotential induced by the oxygen electrode is 48% less than that of oxygen electrode/electrolyte interfacial contact for up to 1171 h. The data‐driven machine learning simulation predicts the RUL of up to 2132 h. The root cause of degradation is overpotential increase induced by oxygen electrode, which accounts for 82.9% of total cell degradation. The success of the failure diagnostic model is demonstrated by its consistency with degradation modes that do not manifest in electrolysis fade during early real operations. This synergistic approach provides valuable insights into practical failure diagnosis of PCECs and has the potential to revolutionize their development by enabling improved performance prediction and material selection for enhanced durability and efficiency.

25 ENERGY STORAGE↗

Structure of the Bastnäsite (001) Surface by Crystal Truncation Rod X-ray Diffraction and Ab Initio Molecular Dynamics: Implications for Separations of a Rare Earth Ore Mineral

Bastnäsite ((Ce,La)FCO 3 ) is the primary mineral source of light rare earth elements, but its surface structure is not well understood. This presents a major challenge in improving beneficiation strategies. In this work, a synergistic combination of X-ray scattering and ab initio molecular dynamics (AIMD) was used to gain atomistic insight into the interfacial structure of bastnäsite. Surface X-ray scattering was used to measure crystal truncation rods (CTRs) of the bastnäsite (001) surface, a significant crystal face with a previously unknown termination. The best-fit atomic-scale model of the CTR data features a carbonate layer at the surface, which is stabilized by the relaxation of carbonate groups from their bulk structural positions. AIMD simulations predict similar surface relaxations, which are shown to be influenced by the protonation of oxygen atoms at the surface. Evidence of ordered water at the interface is also observed in the best-fit model and AIMD simulations. Further, the presence of a carbonate layer at this dominant crystal surface is significant for improving separation technologies because most commonly used ligands utilize anionic functional groups to chelate metal cations at particle surfaces. Without modification, anionic ligands are expected to have poor affinity for the carbonate-terminated (001) surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Super-X and conventional divertor configurations in MAST-U ohmic L-mode; a comparison facilitated by interpretative modelling

Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,\textrm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16 $\,\pm\,7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e = 5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream n e and T e profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,\textrm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

Cloud radiative effects significantly increase wintertime atmospheric blocking in the Euro-Atlantic sector

Reliable simulation, prediction, and complete theoretical understanding of atmospheric blocking remain challenging despite its significant socio-economic impacts. Generations of climate models have notoriously underestimated blocking frequency, particularly over the Euro-Atlantic sector. Identifying factors controlling blocking frequency and dynamics is therefore essential for improving its simulation. Here, using a cloud-locking experiment, we show that cloud radiative effects (CREs) significantly increase the frequency of Euro-Atlantic blocking. CREs enhance upstream diabatic source of wave activity, both directly through longwave heating and indirectly through their feedback on latent heating, with the latter playing the dominant role. The resulting increase in the upstream diabatic source feeds into local wave activity downstream and promotes blocking formation. Qualitatively similar results are shown by multi-model experiments with radiatively inactive clouds to longwave radiation, albeit with a larger impact from mean-state changes. The results underscore the necessity of accurately representing cloud-radiation interactions in weather and climate models for improved prediction of blocking events.

Atmospheric dynamics↗

Role of self-generated magnetic fields in the inertial fusion ignition threshold

Magnetic fields spontaneously grow at unstable interfaces around hot-spot asymmetries during inertial confinement fusion implosions. Although difficult to measure, theoretical considerations and numerical simulations predict field strengths exceeding 5 kT in current National Ignition Facility experiments. Magnetic confinement of electrons then reduces the rate of hot-spot heat loss by >5%. We demonstrate this via magnetic post-processing of two-dimensional xRAGE hydrodynamic simulation data at bang time. We then derive a model for the self-magnetization, finding that it varies with the square of the hot-spot temperature and inversely with the areal density. The self-magnetized Lawson analysis then gives a slightly reduced ignition threshold. Time-dependent hot-spot energy balance models corroborate this finding, with the magnetic field quadrupling the fusion yield for near-threshold parameters. The inclusion of magnetized multi-dimensional fluid instabilities could further alter the ignition threshold and will be the subject of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗