Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex simulations”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

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

ezAlign

The ezAlign is aimed at clustering coarse grain simulation to find common or uncommon occurrences and convert coarse (CG) grained coordinate and topology files to atomistic formats. We use a PointNet based approach to map individual frames of simulation to points in a latent space. These points are then clustered using a variety of clustering methods. Clusters are analyzed to associate them with states in the simulation. Frames can be chosen from these clusters based on proximity to cluster centers. ezAlign takes CG coordinate and topology files and converts and outputs their corresponding atomistic formats using an alignment and relaxation procedure. ezAlign is designed to convert complex, solvated biological systems including lipid membranes with drug-like molecules using GROMACS. A GROMACS checkpoint (.cpt) file is also outputted to enable continuation simulations that retain the equilibrated atomic velocities. Independent atomistic coordinates and topologies for every molecule must already be included in ezAlign/files. A number of commonly simulated biological molecules are currently provided.

Bennett, WilliamF.↗
DOE OSTI2022

Extracting forces from noisy dynamics in dusty plasmas

Extracting environmental forces from noisy data is a common yet challenging task in complex physical systems. Machine learning (ML) represents a robust approach to this problem, yet is mostly tested on simulated data with known parameters. Here we use supervised ML to extract the electrostatic, dissipative, and stochastic forces acting on micron-sized charged particles levitated in an argon plasma (dusty plasma). By tracking the sub-pixel motion of particles in subsequent images, we successfully estimated these forces from their random motion. The experiments contained important sources of non-Gaussian noise, such as drift and pixel-locking, representing a data mismatch from methods used to analyze simulated data with purely Gaussian noise. Our model was trained on simulated particle trajectories that included all of these artifacts, and used more than 100 dynamical and statistical features, resulting in a prediction with 50\% better accuracy than conventional methods. Lastly, in systems with two interacting particles, the model provided non-contact measurements of the particle charge and Debye length in the plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗
DOE OSTI2021

Simulations of applications using diaboloid mirrors

The diaboloid is a reflecting surface that converts a spherical wave to a cylindrical wave. This complex surface may find application in new Advanced Light Source bending-magnet beamlines or in other beamlines that now use toroidal optics for astigmatic focusing. Here, the numerical implementation of diaboloid mirrors is described, and the benefit of this mirror in beamlines exploiting diffraction-limited storage rings is studied by ray tracing. The use of diaboloids becomes especially interesting for the new low-emittance storage rings because the reduction of aberration becomes essential for such small sources. Here, the validity of the toroidal and other mirror surfaces approximating the diaboloid, and the effect of the mirror magnification, are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Short- to Intermediate-Range Structure, Transport, and Thermophysical Properties of LiF–NaF–ZrF4 Molten Salts

LiF–NaF–ZrF 4 multicomponent molten salts are identified as promising candidates for coolant salts in molten salt reactors and advanced high-temperature reactors. This study focused on low-melting point salt compositions of interest: 38LiF–51NaF–11ZrF 4 , 42LiF–29NaF–29ZrF 4 , and 26LiF–37NaF–37ZrF 4 . Ab-initio molecular dynamics (AIMD) calculations were performed and compared with available experimental data to assess the ability of rigid ion models (RIM) to reproduce short- to intermediate-range structure, transport, and thermophysical properties of the LiF–NaF–ZrF 4 salt mixtures. It is found that as ZrF 4 mol% increases, the average cation–anion coordination number (CN) of monovalent cations (Li + , Na + ) obtained from RIM calculations decreases, while multivalent Zr 4+ CN varied from 15% to 19% in comparison to corresponding AIMD values. In addition, RIM is found to predict the existence of 7, 8, and 9 coordinated fluorozirconate complexes, while AIMD and the available experimental data showed an occurrence of 6, 7, and 8 coordinated complexes in the melt. The intermediate-range structure analysis revealed that while the RIM parameters are able to reproduce a local structure for lower ZrF 4 mol% salts such as in 38LiF–51NaF–11ZrF 4 , an extensive fluorozirconate network formation is observed in RIM simulations for higher ZrF 4 mol% compositions. The network generated by RIM parameters is found to be mainly connected by “corner-sharing” fluorozirconate complexes as opposed to both “edge-sharing” and “corner-sharing” connectively portrayed by AIMD. It is found that a close agreement between AIMD and the RIM salt structure for the 11-mol% ZrF 4 salt resulted in good agreement in the calculated Zr diffusivities and the viscosity values. However, due to the inaccurate short- to intermediate-range structure prediction by RIM for higher ZrF 4 mol% compositions, thermophysical properties such as densities and heat capacity differ by up to 26% and 27%, respectively, upon comparison with AIMD and experimental values. Also, the network-dominated properties such as diffusion coefficients and viscosities differed by up to two and three orders of magnitude, respectively. This study signifies the importance of accurate salt structure generation for an accurate prediction of transport and thermophysical properties of multicomponent molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Temperature-dependent kinetic pathways featuring distinctive thermal-activation mechanisms in structural evolution of ice VII

Significance The kinetics of the structural evolution from crystalline ice (ice VII or ice VIII) toward thermodynamically stable ice I at ambient pressure have been characterized using time-resolved in situ X-ray diffraction. The time resolution captures events that were not observed previously and reveals a complex kinetic pathway with three distinctive transitions, i.e., amorphization, amorphous–amorphous, and crystallization of amorphous ice into ice I. The first two transitions exhibit different thermally activated mechanism above and below 110 K. Molecular-dynamics simulations show the transition involves a large change in density and substantial displacements of water at the nanoscale. This study presents a perspective on the metastability and complexity of the energy landscape in temperature/time-dependent structural evolution in amorphous ices.

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

About the Complex Natural and Engineered Systems Pillar

Los Alamos National Laboratory is a world leader in applying multi-disciplinary science to complex systems within the Complex Natural and Engineered Systems (CNES) challenge areas. Our search for solutions requires science and technology innovation, as well as an integrated experiment, theory, and modeling and simulation approach. Our research and development spans from improving engineered systems such as nuclear weapons and the power grid, to understanding the interface of human and engineered systems from the subsurface to space, to studying how complex natural systems such as disease and climate impact humanity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

The design of a photoionization front experiment using the Z-Machine as a driving source and estimated measurements

Radiation-driven heat fronts are present in the early universe during reionization, the circumstellar medium of supernovae, and in high-energy-density physics experiments. Dedicated experiments to observe and diagnose the behavior of these types of heat fronts can improve our understanding of these phenomena. A simulation study of photoionization fronts using the HELIOS-CR radiation hydrodynamics code provides an experimental design for the Z-Machine at Sandia National Laboratory using a measurement-calibrated input radiation flux to drive the photoionization front. The simulations use detailed atomic physics and non-diffusive radiation transport in 1D to determine an optimal gas pressure of 0.75 atm for an experiment in N gas as well as the effects of increasing the thickness of the window that seals the gas cell. Post-processing of these simulations demonstrates that ratios of atomic rate coefficients place the heat front in a physics regime where photoionization dominates the energy deposition. Here, to see the sensitivity of the simulations to changes in the model and spatial grid, this analysis performed resolution, atomic model detail, and radiation transport angular grid studies showing less than 10% deviation from the nominal model for increased complexity, when possible. An effort to emulate 3D geometric effects on the radiation flux using an artificial attenuation scheme has shown that, even for conservative estimates of the flux, simulations still produce a photoionization front. Estimations of a streaked, visible spectroscopy measurement using SPECT3D showed that line emission measurements are present early in time and that later in time thermal emission should become dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Partitioning Quantum Chemistry Simulations with Clifford Circuits

Current quantum computing hardware is restricted by the availability of only few, noisy qubits which limits the investigation of larger, more complex molecules in quantum chemistry calculations on quantum computers in the near term. Here, in this work, we investigate the limits of their classical and near-classical treatment while staying within the framework of quantum circuits and the variational quantum eigensolver. To this end, we consider naive and physically motivated, classically efficient product ansatz for the parametrized wavefunction adapting the separable-pair ansatz form. We combine it with post-treatment to account for interactions between subsystems originating from this ansatz. The classical treatment is given by another quantum circuit that has support between the enforced subsystems and is folded into the Hamiltonian. To avoid an exponential increase in the number of Hamiltonian terms, the entangling operations are constructed from purely Clifford or near-Clifford circuits. While Clifford circuits can be simulated efficiently classically, they are not universal. In order to account for missing expressibility, near-Clifford circuits with only few, selected non-Clifford gates are employed. The exact circuit structure to achieve this objective is molecule-dependent and is constructed using simulated annealing and genetic algorithms. We demonstrate our approach on a set of molecules of interest and investigate the extent of our methodology’s reach.

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

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Discrete-event Simulation Process Model for the Pyrochemical Processing of Plutonium at Los Alamos National Laboratory

The pyrochemical metal production operations that occur in the Plutonium Facility at Los Alamos National Laboratory perform plutonium purification with the aim to provide plutonium metal for a variety of defense- and non-defense missions within the National Nuclear Security Administration. The demands and constraints associated with the pyrochemical processing of plutonium are complex, making decision analyses challenging for program managers who require plutonium production for their mission applications. The construction of a discrete-event simulation process model is proposed to measure and report the process capacity, material throughput, equipment requirements, and dose accumulation for operators of the pyrochemical metal production operations. The process model, constructed in the ExtendSim™ software, will represent the cause-and-effect relationships between the pyrochemical processing environment and the process constraints, including criticality limitations, material control and accountability measures, chemical analysis requirements, and equipment availability. An accurate representation of the pyrochemical metal production process capacity through simulation modeling will be helpful to program managers in their efforts to forecast plutonium availability for mission applications. Furthermore, the proposed process model will be vital for future analyses that will measure the interactions between the pyrochemical metal production operations and the aqueous reprocessing operations and their ability to minimize transuranic waste disposal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The “Hole” Story in Ionized Water from the Perspective of Ehrenfest Dynamics

The radiolysis of liquid water and the radiation-matter interactions that happen in aqueous environments are important to the elds of chemistry, materials, and environmental sciences, as well as biological and physiological response to extreme conditions and medical treatments. The initial stage of radiolysis is the ultrafast response, or hole dynamics, that trigger chemical processes within complex energetic landscapes that may include reactivity. A fundamental understanding necessitates the use of theoretical methods that are capable of simulating both ultrafast coherence and non-adiabatic energy transfer pathways. In this work, we carry out an ab initio Ehrenfest dynamic study in order to provide a more complete description of the ultrafast dynamics and reactive events initiated by photoionization of water. Furthermore, after sudden ionization, the range of processes that include hole trapping and transfer, large OH oscillations, proton transfer and subsequent relay, formation of the metastable Zundel complex, and long-lived coherence, are identied and new insight into their driving forces are elucidated.

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

Quantum Computing and Simulations for Energy Applications

While quantum computing (QC) is considered as a paradigm shift in our basic understanding of physical computation, effective implementation of QC in energy applications also depends on progress and development in the dimensions of both QC hardware and algorithms. To fully address the status and future challenges of QC applied within the energy sector, in this presentation, we firstly summarize recent advancements on the applications of QC to energy infrastructure and materials, complex energy system processes, advanced manufacturing, and energy system security. Then, we will demonstrate the results of QC performed both on a simulator and a quantum device targeting on energy-related applications.

Paudel, Hari P.↗
DOE OSTI2021

Modeling of aqueous species interaction energies prior to nucleation in cement-based gel systems

Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.

36 MATERIALS SCIENCE↗