Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION METHODS”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

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

Assessing reservoir performance for geologic carbon sequestration in offshore saline reservoirs

Carbon capture and storage (CCS) is a technological strategy to reduce CO 2 emissions from hydrocarbon-based energy production and non-energy industrial emitters, like fertilizer and cement production. To expand the geographic opportunities for CCS, offshore geologic settings are increasingly considered for saline CO 2 storage, particularly where legacy oil and gas infrastructure may be repurposed for CO 2 transport and injection. In this context, the Gulf of Mexico (GoM), United States (U.S.), may offer tremendous opportunity for CCS because the surrounding states are among the largest CO 2 producers in the U.S. and there has been extensive oil and gas development over the last half century. Nevertheless, there remains significant uncertainty in how GoM reservoir/seal systems will respond to industrial-scale CO 2 injections, particularly in the context of meso-scale (1–10s of meters) permeability variations that are difficult to identify prior to injection. Such permeability variations can have significant impacts on CO 2 migration and fluid pressure propagation, both of which govern the overall effectiveness of the CCS project. This research uses ensemble simulation methods to quantify the influence of spatially variable and a priori unknown permeability fields on CO 2 plume development and fluid pressure propagation when CO 2 is injected at 1 MMt per year for ten years. Results show a number of characteristic patterns that reveal the substantial influence of both near- and far-field permeability on reservoir performance. By focusing on basin characteristics that are common among offshore basins worldwide, results from this study suggest that offshore CCS may be a feasible carbon management strategy in regions where offshore basins are proximal to industrial CO 2 sources.

42 ENGINEERING↗
DOE OSTI2024

Modeling Electric Vehicle Charging Load Using Origin-Destination Data

The accelerating adoption of electric vehicles (EVs) poses challenges to the power grid, necessitating precise representation of mobility patterns for effective infrastructure upgrades. Traditional simulation-based charging demand estimation faces limitations in generating trip chains reflective of actual travel patterns without complex network modeling. Hence, an innovative agent-based trip chain generation model is introduced to overcome these challenges. Drawing from the National Household Travel Survey (NHTS) and the NextGen NHTS origin-destination add-on data for Clarke County, Georgia, this study proposes a simulation method capturing both temporal and spatial mobility patterns without relying on extensive network topology data. The resulting trip chains predict EV charging load at the Census Block Group level, validated with a 1.03 correlation to actual trip counts, affirming their reflective accuracy. Two charging scenarios, residential-only and charging-everywhere, reveal distinct demand profiles. The charging-everywhere scenario aligns closely with the trip profile, while the residential-only scenario exhibits an afternoon peak slightly surpassing the former. This study contributes a data-driven charging demand estimation methodology, offering critical insights for grid resiliency planning amid the evolving landscape of EV adoption.

Pan, Melrose↗
DOE OSTI2013

Molecular Dynamics Study of Freezing Point and Solid-Liquid Interfacial Free Energy of Stockmayer Fluids

Freezing temperatures of Stockmayer fluids with different dipolar strength at zero pressure are estimated and computed using three independent molecular-dynamics (MD) simulation methods, namely, the superheating-undercooling method, the constant-pressure and constant-temperature (NPT) two phase coexistence method, and the constant-pressure and constant-enthalpy (NPH) coexistence method. The best estimate of the freezing temperature (in reduced unit) for the Stockmayer (SM) fluid with a reduced dipole moment is 0.656±0.001, 0.726±0.002 and 0.835±0.005, respectively. The freezing temperature increases with the dipolar strength. The solid-liquid interfacial free energies of the (111), (110) and (100) interface are calculated for the first time using two independent methods, namely, the cleaving-wall method and the interfacial fluctuation method. Both methods predict that the interfacial free energy increases with the dipole moment. Although the interfacial fluctuation method suggests a weaker interfacial anisotropy, particularly for strongly dipolar SM fluids, both methods predicted the same trend of interfacial anisotropy, that is, .

Wang, J.↗
DOE OSTI2014

Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways

Solid-liquid interfaces are at the heart of many modern-day technologies and provide a challenge to many materials simulation methods. A realistic first-principles computational study of such systems entails the inclusion of solvent effects. In this work, we implement an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code Vienna ab initio Software Package. The implicit solvation model follows the framework of joint density functional theory. We describe the framework, our algorithm and implementation, and benchmarks for small molecular systems. We apply the solvation model to study the surface energies of different facets of semiconducting and metallic nanocrystals and the S N 2 reaction pathway. We find that solvation reduces the surface energies of the nanocrystals, especially for the semiconducting ones and increases the energy barrier of the S N 2 reaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Design of a low-emittance lattice for 3.5 GeV synchrotron light source

We investigate a lattice design for a low-emittance, high-brilliance intermediate energy synchrotron light source that is being intensively performed in the world. We present the design results in detail for a non-achromatic double bend achromat lattice with an emittance of 1.1 nm and 597.6 m circumference. Each superperiod has 4 straight sections, 7 m, 5.4 m, and 2 × 2.8 m in length, giving a total of 60 straight sections in the ring. The lattice is designed to be flexible and also provides an achromatic optics with an emittance of 3.2 nm. The lattice provides high brilliance at a photon energy of a few tens of keV, which meets the requirements of hard X-ray synchrotron radiation users. We investigated the dynamic aperture in the lattice with machine errors using a simulation method that achieves the optimal tune for the ring. The low-emittance lattice provides sufficient dynamic aperture after closed orbit distortion correction. We show an estimate of variations of emittance, energy spread, and dynamic aperture due to a wiggler in the straight sections. We present the design of an injection scheme for full energy injection, the space allocation in injection section, and the particle dynamics of injected beam. Particle tracking after beam injection was performed to examine the validity of the injection scheme. Our design lattice shows good optimization in terms of the emittance, number of straight sections, brilliance, and circumference as a light source for 3.5 GeV beam energy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Parametric study of broadband terahertz radiation generation based on interaction of two-color ultra-short laser pulses

In this work, using a two-dimensional kinetic model based on particle in cell-Monte Carlo collision simulation method, the influence of different parameters on the broadband intense Terahertz (THz) radiation generation via application of two-color laser fields, i.e., the fundamental and second harmonic modes, is studied. These two modes are focused into the molecular oxygen (O{sub 2}) with uniform density background gaseous media and the plasma channels are created. Thus, a broadband THz pulse that is around the plasma frequency is emitted from the formed plasma channel and co-propagates with the laser pulse. For different laser pulse shapes, the THz electric field and its spectrum are both calculated. The effects of laser pulse and medium parameters, i.e., positive and negative chirp pulse, number of laser cycles in the pulse, laser pulse shape, background gas pressure, and exerted DC electric field on THz spectrum are verified. Application of a negatively chirped femtosecond (40 fs) laser pulse results in four times enhancement of the THz pulse energy (2 times in THz electric field). The emission of THz radiation is mostly observed in the forward direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Freely cooling granular gases with short-ranged attractive potentials

We treat the case of an undriven gas of inelastic hard-spheres with short-ranged attractive potentials via an extension of the pseudo-Liouville operator formalism. New evolution equations for the granular temperature and coordination number are obtained. The granular temperature exhibits deviation from both Haff’s law and the case of long-ranged potentials. We verify this departure using soft-sphere discrete element method simulations. Excellent agreement is found for the duration of the simulation even beyond where exclusively binary collisions are expected. Simulations show the emergence of strong spatial-velocity correlations on the length scale of the last peak in the pair-correlation function but do not show strong correlations beyond this length scale. We argue that molecular chaos may remain an adequate approximation if the system is modelled as a Smoluchowski type equation with aggregation and break-up processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Frequency spirals

We study the dynamics of coupled phase oscillators on a two-dimensional Kuramoto lattice with periodic boundary conditions. For coupling strengths just below the transition to global phase-locking, we find localized spatiotemporal patterns that we call “frequency spirals.” These patterns cannot be seen under time averaging; they become visible only when we examine the spatial variation of the oscillators' instantaneous frequencies, where they manifest themselves as two-armed rotating spirals. In the more familiar phase representation, they appear as wobbly periodic patterns surrounding a phase vortex. Unlike the stationary phase vortices seen in magnetic spin systems, or the rotating spiral waves seen in reaction-diffusion systems, frequency spirals librate: the phases of the oscillators surrounding the central vortex move forward and then backward, executing a periodic motion with zero winding number. We construct the simplest frequency spiral and characterize its properties using analytical and numerical methods. Simulations show that frequency spirals in large lattices behave much like this simple prototype.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Simulation of stellarator divertors

The nested magnetic surfaces that confine a fusion plasma can be developed to be bounded by a limiter or a divertor. For a limiter, confining surfaces extend until they intercept a part of the surrounding structure. For a divertor, an outermost confining magnetic surface exists, which is well separated from the surrounding structures. The only designs that are thought to be fusion relevant have divertors that direct field lines from the plasma edge into chambers where the particle exhaust can be pumped and the residual heat exhaust can be handled. The topological properties of magnetic field lines just outside the outermost confining surface determine much of the physics of divertors. Axisymmetric tokamak divertors are well-known, and the outermost confining surface is defined by a sharp separatrix. The topology of the magnetic field lines associated with a stellarator divertor is far more subtle. Related subtleties arise in tokamak divertors when subjected to sufficiently strong non-axisymmetric perturbations.Here, an efficient simulation method for carrying out topological studies of non-axisymmetric divertors is created and illustrative results are given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Dissipative particle dynamics for directed self-assembly of block copolymers

The dissipative particle dynamics (DPD) simulation method has been shown to be a promising tool to study self-assembly of soft matter systems. In particular, it has been used to study block copolymer (BCP) self-assembly. However, previous parameterizations of this model are not able to capture most of the rich phase behaviors of BCPs in thin films nor in directed self-assembly (chemoepitaxy or graphoepitaxy). Here, we extend the applicability of the DPD method for BCPs to make it applicable to thin films and directed self-assembly. Our new reparameterization not only is able to reproduce the bulk phase behavior but also manages to predict thin film structures obtained experimentally from chemoepitaxy or graphoepitaxy. A number of different complex structures, such as bilayer nanomeshes, 90° bend structures, circular cylinders/lamellae and Frank-Kasper phases directed by trenches, and post arrays or chemically patterned substrates, have all been reproduced in this work. This reparameterized DPD model should serves as a powerful tool to predict BCP self-assembly, especially in some complex systems where it is difficult to implement self-consistent field theory.

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

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Phase-field modeling of aging-induced microstructure evolution in pentaerythritol tetranitrate thin films and ramifications for shock initiation

Aging of energetic materials may change performance and affect their safety and reliability, but the relationship between microstructure changes induced by aging and consequent performance changes has not been fully established. This work presents results of phase-field method simulations used to model microstructure evolution of vapor-deposited pentaerythritol tetranitrate (PETN) thin films. Simulated aging is shown to induce grain coarsening and substantial changes of the configuration of porosity in the film: Specifically, we show that porosity tends to concentrate in large pores to a greater degree in aged films, a state that is arrived at by closure or consolidation of small pores. To evaluate the performance of the as-deposited and aged films, we perform two-dimensional hydrocode flyer-film impact simulations that incorporate the phase-field output microstructures directly, permitting us to connect features therein to changes in reactivity, a key metric of energy output for shock initiation. The results demonstrate that declining sensitivity obtained for the simulated aged films can be correlated with the loss of fine-structured pores relatively early in the aging process, while long-term microstructure evolution that gradually alters the shape of large, branching pores is less impactful. Finally, we discuss commonalities and discrepancies between our simulation results and high-throughput initiation experiments on shock initiation of aged PETN thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Formation of HONO from the NH 3 -promoted hydrolysis of NO 2 dimers in the atmosphere

One challenging issue in atmospheric chemistry is identifying the source of nitrous acid (HONO), which is believed to be a primary source of atmospheric “detergent” OH radicals. Herein, we show a reaction route for the formation of HONO species from the NH 3 -promoted hydrolysis of a NO 2 dimer (ONONO 2 ), which entails a low free-energy barrier of 0.5 kcal/mol at room temperature. Our systematic study of HONO formation based on NH 3 + ONONO 2 + n H 2 O and water droplet systems with the metadynamics simulation method and a reaction pathway searching method reveals two distinct mechanisms: ( i ) In monohydrates ( n = 1), tetrahydrates ( n = 4), and water droplets, only one water molecule is directly involved in the reaction (denoted the single-water mechanism); and ( ii ) the splitting of two neighboring water molecules is seen in the dihydrates ( n = 2) and trihydrates ( n = 3) (denoted the dual-water mechanism). A comparison of the computed free-energy surface for NH 3 -free and NH 3 -containing systems indicates that gaseous NH 3 can markedly lower the free-energy barrier to HONO formation while stabilizing the product state, producing a more exergonic reaction, in contrast to the endergonic reaction for the NH 3 -free system. More importantly, the water droplet reduces the free-energy barrier for HONO formation to 0.5 kcal/mol, which is negligible at room temperature. We show that the entropic contribution is important in the mechanism by which NH 3 promotes HONO formation. Furthermore, this study provides insight into the importance of fundamental HONO chemistry and its broader implication to aerosol and cloud processing chemistry at the air–water interface.

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

Path Toward a Breeding, Proliferation-Resistant, Thermal-Spectrum Molten Salt Reactor

Liquid-fueled, thermal-spectrum molten salt breeder reactors (TS-MSBRs) offer the potential for affordable, safe, inexhaustible energy with minimal potential for nuclear material misuse and without significant actinide waste generation. Realizing the full set of TS-MSBR capabilities is only now becoming possible with the advent of advanced fuel-salt processing techniques, improved materials, and a more detailed understanding of fuel-salt properties. Additionally, modern higher-fidelity modeling and simulation methods enable a more detailed evaluation of TS-MSBR design options. TS-MSBRs, however, remain immature and will require substantial, sustained development resources.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗