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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dynamical dark energy from lattice quantum gravity

We study the behavior of the vacuum in Euclidean dynamical triangulations (EDT). Algorithmic improvements and better lattice spacing determinations allow us to test the properties of the emergent de Sitter geometries of our simulations to higher precision than previously possible. Although the agreement with de Sitter is good, the improved precision reveals deviations that can be interpreted as nontrivial vacuum dynamics, well described by a cosmological constant that runs with scale. The simulations show that the dominant running is quadratic and that the scale can be identified with the Hubble rate. Several key cross-checks support this picture, including consistent results across multiple lattice spacings and the fact that the null energy condition is not violated. The parameters of the running are fully determined by simulations, enabling predictions when extrapolated to the scales relevant for our Universe. This leads to a model for dark energy that is compatible with current observations, but which predicts deviations from the Standard Model of cosmology at the 𝒪⁡(10 −3 ) level in observables that could be tested with future improvements in precision measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2020

Lone Pair Rotational Dynamics in Solids

Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom — here lone pairs — exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Finally, such fluctuations in the charge density should be observable, in principle via synchrotron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Electric-Field-Induced Spatially Dynamic Heterogeneity of Solvent Motion and Cation Transference in Electrolytes

While electric fields primarily result in migration of charged species in electrolytic solutions, the solutions are dynamically heterogeneous. Solvent molecules within the solvation shells of the cation will be dragged by the field while free solvent molecules will not. Here we combine electrophoretic NMR measurements of ion and solvent velocities under applied electric fields with molecular dynamics simulations to interrogate different solvation motifs in a model liquid electrolyte. Measured values of the cation transference number $(t^0_+)$ agree quantitatively with simulation-based predictions over a range of electrolyte concentrations. Solvent-cation interactions strongly influence the concentration-dependent behavior of $t^0_+$. We identify a critical concentration at which most of the solvent molecules lie within solvation shells of the cations. The dynamic heterogeneity of solvent molecules is minimized at this concentration where $t^0_+$ is approximately equal to zero.

42 ENGINEERING↗

DOE OSTI2022

Direct Measurement of Ice-Ablator Interface Motion for Instability Mitigation in Indirect Drive ICF Implosions

In indirect drive inertial confinement fusion (ICF) implosions hydrodynamic instability growth at the imploding capsule ablator-DT fuel interface can reduce fuel compressibility and inject ablator into the hot spot hence reducing hot spot pressure and temperature. As a mitigation strategy, a gentle acceleration of this interface is predicted by simulations and theory to significantly reduce this instability growth in the early stage of the implosion. We have performed high-contrast, time-resolved x-ray refraction enhanced radiography (RER) to accurately measure the level of acceleration as a function of the initial laser drive time history for indirect-drive implosions on the National Ignition Facility. Here, we demonstrate a transition from no acceleration to 20 ± 1.8μm ns -2 acceleration by tweaking the drive that should reduce the initial instabilities by an order of magnitude at high modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Polaritonic Chern Insulators in Monolayer Semiconductors

Systems with strong light-matter interaction open up new avenues for studying topological phases of matter. Examples include exciton polaritons, mixed light-matter quasiparticles, where the topology of the polaritonic band structure arises from the collective coupling between matter wave and optical fields strongly confined in periodic dielectric structures. Distinct from light-matter interaction in a uniform environment, the spatially varying nature of the optical fields leads to a fundamental modification of the well-known optical selection rules, which were derived under the plane wave approximation. Here, in this study, we identify polaritonic Chern insulators by coupling valley excitons in transition metal dichalcogenides to photonic Bloch modes in a dielectric photonic crystal slab. We show that polaritonic Dirac points, which are markers for topological phase transition points, can be constructed from the collective coupling between valley excitons and photonic Dirac cones in the presence of both time-reversal and inversion symmetries. Lifting exciton valley degeneracy by breaking time-reversal symmetry leads to gapped polaritonic bands with nonzero Chern numbers. Through numerical simulations, we predict polaritonic chiral edge states residing inside the topological gaps. Our Letter paves the way for the further study of strong exciton-photon interaction in nanophotonic structures and for exploring polaritonic topological phases and their practical applications in polaritonic devices.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Kilonova light-curve interpolation with neural networks

Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22 248 high-fidelity simulations (composed of 412 unique ejecta parameter combinations evaluated at 54 viewing angles), we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multiband inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

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↗

DOE OSTI2021

Early Evaluation of Fugaku A64FX Architecture Using Climate Workloads

The Energy Exascale Earth System Model (E3SM) Project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project that targets efficient utilization of U.S. Department of Energy’s (DOE) supercomputers to meet the science needs of the nation and the mission needs of DOE. This work focuses on our early evaluation of the A64FX architecture on Fugaku supercomputer using E3SM benchmarks. We will present results that track hardware trends, facilitate architecture comparison and the specific impact on our workload using an atmospheric model benchmark. We have two variants of the code written in Fortran and C++/Kokkos respectively which were used to collect data on a variety of CPU and GPU platforms. Furthermore, we have conducted a comparative evaluation of the compilers on the A64FX architecture and found GNU to be the best performer for our workload. Our experience so far indicates that Fugaku/A64FX shows promising energy efficiency (performance/Watt) with further performance gains possible through architecture-aware optimization efforts.

Sreepathi, Sarat↗

DOE OSTI2023

Warming decreases desert ecosystem functioning by altering biocrusts in drylands

Warming and precipitation fluctuations are changing desert ecosystems in global drylands. However, the effects of climate change on keystone species such as cryptogamic biocrust in drylands remain relatively under-investigated, even though biocrusts play a vital role in desert ecosystems. Here we conducted a long-term experiment (14 years) to simulate the responses of two main types of biocrusts to warming coupled with reduction in precipitation that was achieved by open-top chambers (OTCs) to simulate the predicted warming and precipitation decreasing under climate change scenario. We also conduct a structural analysis to evaluate the resulting changes in desert ecosystem functioning (carbon and nitrogen cycling). Neither warming and corresponding rainfall reduction treatments had a negative effect on lichen species richness, but both treatments reduced lichen cover and biomass. The negative effects of warming on moss-dominated crusts were much greater than those on lichen-dominated crusts. Although mosses and lichens had varying degree responses to warming, the loss of mosses and decreased lichen cover and biomass, as well as the shortening of the wet time, resulted in a reduction in carbon and nitrogen fixation, soil enzyme activity and water-holding capacity of biocrusts and topsoil. These impacts collectively change the water balance of drylands and weaken the hydrological and biogeochemical function of biocrusts. Synthesis and applications: Results from this long-term experiment suggest that the ecosystem C and N cycling and water balance of global drylands may be highly impacted by climate change, in part because of the response of biocrusts, which contribute an important implication for both dryland restoration and earth system dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Bayesian, Multifidelity Operator Learning for Complex Engineering Systems–A Position Paper

Abstract Deep learning has significantly improved the state-of-the-art in computer vision and natural language processing, and holds great potential to design effective tools for predicting and simulating complex engineering systems. In particular, scientific machine learning seeks to apply the power of deep learning to scientific and engineering tasks, with operator learning (OL) emerging as a particularly effective tool. OL can approximate nonlinear operators arising in complex engineering systems, making it useful for simulating, designing, and controlling those systems. In this position paper, we provide a comprehensive overview of OL, including its potential applications to complex engineering domains. We cover three variations of OL approaches: deterministic OL for modeling nonautonomous systems, OL with uncertainty quantification (UQ) capabilities, and multifidelity OL. For each variation, we discuss drawbacks and potential applications to engineering, in addition to providing a detailed explanation. We also highlight how multifidelity OL approaches with UQ capabilities can be used to design, optimize, and control engineering systems. Finally, we outline some potential challenges for OL within the engineering domain.

Computer Science↗

DOE OSTI2022

RASER MRI: Magnetic resonance images formed spontaneously exploiting cooperative nonlinear interaction

The spatial resolution of magnetic resonance imaging (MRI) is limited by the width of Lorentzian point spread functions associated with the transverse relaxation rate 1/T 2 * . Here, we show a different contrast mechanism in MRI by establishing RASER (radio-frequency amplification by stimulated emission of radiation) in imaged media. RASER imaging bursts emerge out of noise and without applying radio-frequency pulses when placing spins with sufficient population inversion in a weak magnetic field gradient. Small local differences in initial population inversion density can create stronger image contrast than conventional MRI. This different contrast mechanism is based on the cooperative nonlinear interaction between all slices. On the other hand, the cooperative nonlinear interaction gives rise to imaging artifacts, such as amplitude distortions and side lobes outside of the imaging domain. Contrast mechanism and artifacts are explored experimentally and predicted by simulations on the basis of a proposed RASER MRI theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Functional soft materials from blue phase liquid crystals

Blue phase (BP) liquid crystals are chiral fluids wherein millions of molecules self-assemble into cubic lattices that are on the order of hundred nanometers. As the unit cell sizes of BPs are comparable to the wavelength of light, they exhibit selective Bragg reflections in the visible. The exploitation of the photonic properties of BPs for technological applications is made possible through photopolymerization, a process that renders mechanical robustness and thermal stability. We review here the preparation and characterization of stimuli-responsive, polymeric photonic crystals based on BPs. We highlight recent studies that demonstrate the promise that polymerized BP photonic crystals hold for colorimetric sensing and dynamic light control. We review using Landau–de Gennes simulations for predicting the self-assembly of BPs and the potential for using theory to guide experimental design. Finally, opportunities for using BPs to synthesize new soft materials, such as highly structured polymer meshes, are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Disentangling the sources of ionizing radiation in superconducting qubits

Abstract Radioactivity was recently discovered as a source of decoherence and correlated errors for the real-world implementation of superconducting quantum processors. In this work, we measure levels of radioactivity present in a typical laboratory environment (from muons, neutrons, and $$\gamma $$ γ -rays emitted by naturally occurring radioactive isotopes) and in the most commonly used materials for the assembly and operation of state-of-the-art superconducting qubits. We present a GEANT-4 based simulation to predict the rate of impacts and the amount of energy released in a qubit chip from each of the mentioned sources. We finally propose mitigation strategies for the operation of next-generation qubits in a radio-pure environment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Analytical and Numerical Analysis of Lithium Plating Onset in Single and Bilayer Graphite Electrodes during Fast Charging

Lithium plating is commonly observed in anodes charged at fast rates, and can lead to capacity loss and battery safety issues. The increased risk of plating has been attributed to transport limitations, and architectured electrodes may reduce plating risk. However, while theoretical studies have shown that reaction non-uniformity arises due to interplay of transport limitations, anode open circuit voltage behavior and reaction kinetics, its effect on lithium plating has not been studied. We use analytic and numerical simulations to predict onset of plating in graphite anode half-cells at high C-rates and demonstrate how anodes with layered porosities can delay plating. Simplified analytical models identify trends for plating onset and predictions are calibrated against numerical models. A calibrated numerical model of graphite demonstrates qualitative agreement with analytical model predictions. This reaction inhomogeneity mechanism occurs in the absence of lithium ion depletion, indicating that these mechanisms may contribute to capacity loss independently or simultaneously. A bilayer model of graphite exhibits delayed plating onset, and an optimization procedure is presented. This theoretical work presents quantitative and mechanistic insight on how reaction inhomogenity affects lithium metal plating onset and can be used as a guide to engineer anodes resistant to lithium plating.

Mijailovic, Aleksandar S.↗

DOE OSTI2026

FROMP-from-Monomer

This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.

Chua, Lauren [Massachusetts Inst. of Technology (M↗

DOE OSTI2023

MAM4xx

SAND2023-05377O The Modal Aerosol Model with 4 fixed modes (MAM4xx) is a performance-portable, C++ implementation for modeling the chemical reactions of aerosols in the atmosphere. MAM4xx is part of the Energy Exascale Earth System Model (E3SM). The E3SM project is an ongoing Earth system modeling, simulation, and prediction project that optimizes the use of DOE laboratory resources to meet the science needs of the nation and the mission needs of DOE. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Overfelt, James↗