Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAMMPS”

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

Particle Sensitivity Analysis

We propose to develop a computational sensitivity analysis capability for Monte Carlo sampling-based particle simulation relevant to Aleph, Cheetah-MC, Empire, Emphasis, ITS, SPARTA, and LAMMPS codes. These software tools model plasmas, radiation transport, low-density fluids, and molecular motion. Our report demonstrates how adjoint optimization methods can be combined with Monte Carlo sampling-based adjoint particle simulation. Our goal is to develop a sensitivity analysis to drive robust design-based optimization for Monte Carlo sampling-based particle simulation - a currently unavailable capability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Exploring the Role of Viscosity in Inertial Confinement Fusion Implosions [Slides]

For one-component materials, adding pickets to the pulse shape did not affect the viscosity substantially. While Be seems to be most effective in dampening instabilities while Cr is the least effective. In MD-simulation of TCP, LAMMPS TCP input deck was successfully verified against OCP simulation, and we successfully calculated the viscosity of a CH TCP for a range of Γ values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

CephFS experiments on stria.sandia.gov

This report is an institutional record of experiments conducted to explore performance of a vendor installation of CephFS on the SNL stria cluster. Comparisons between CephFS, the Lustre parallel file system, and NFS were done using the IOR and MDTEST benchmarking tools, a test program which uses the SEACAS/Trilinos IOSS library, and the checkpointing activity performed by the LAMMPS molecular dynamics simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal Neutron Scattering Cross Sections for Graphitic Amorphous Carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm 3 . From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical features and the optical graphitic peak at 0.25 eV for the ideal graphite DOS disappear for graphite-like amorphous carbon, which shows good agreement with the expected structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2025

Development and Validation of Dense Discrete Phase Flow Models for Concentrating Solar Power Particle Receivers

This report summarizes a recent project aimed at developing and validating the necessary tools to enable more accurate modeling of denser and more complex particle flows in next-generation particle receivers used in concentrating solar power towers. A newly developed CFD/DEM simulation capability was created by coupling existing the modeling and simulation tools Sierra and LAMMPS. This new capability permitted the inclusion of additional physics for particle drag and particle collisions to model

14 SOLAR ENERGY↗
DOE OSTI2022

Simulations of graphite boronization: A molecular dynamics study of amorphization resulting from bombardment

The molecular dynamics code LAMMPS was used to simulate the bombardment of a graphite structure by atomic boron with impact energies ranging from 50–250 eV. The transient structural evolution, penetration depth, and amorphous layer thickness were analyzed. Simulations show that larger impact energies lead to a greater volume of amorphization and penetration of boron, but that the growth rate of the amorphous layer decreases with increasing fluence. Furthermore, the change in surface chemistry of the amorphized structures was studied using the ReaxFF formalism, which found that the amorphization process introduces dangling bonds thus increasing reactivity in the amorphous region.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Effect of Grain Boundary Misorientation on Spall Strength in Ta via Shock-Free Simulations with Relatively Few Atoms

A suite of 37 molecular dynamics simulations is conducted at two system sizes to systematically characterize the role of grain boundary (GB) misorientation on spall strength in pure BCC tantalum (Ta). The systems studied consist of bicrystals with a single [110] symmetric tilt grain boundary. Two loading conditions are compared: (i) homogeneous extension under uniaxial strain simulated in this study and (ii) piston/flyer impact of sample, which induces heterogeneous deformation via shockwave propagation along the length of the sample. The piston/flyer impact is taken from the literature and run on the same set of GB misorientation angles using LAMMPS. The major finding here is that both methods result in similar spall strength predictions, but the homogeneous extension method generally requires two to three orders of magnitude fewer atoms and similar reductions in computational costs. Spall strength results systematically overpredict using this method, by about 10% for the dataset three orders of magnitude smaller than piston/flyer simulations, and 5% for the dataset two orders of magnitude smaller. Lastly, the effect of system size and pre-compression magnitude on spall strength is systematically characterized.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data from "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions"

Data and input files related to the paper "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions" (https://doi.org/10.1116/6.0004027). This includes the final DP model used in all simulations, training data set, example input files to run DeepMD (with LAMMPS), and data tables summarizing the results obtained from the simulations.

machine learning models↗
DOE OSTI2024

Data and Code for Atomic Scale Etching of Diamond: Insights from Molecular Dynamics Simulations

This work investigates the effects of argon ions, hydrogen atoms, and hydrogen ions on the diamond (100) surface using classical molecular dynamics simulations. The purpose of this investigation was to asses plasma processing techniques for applications in quantum device manufacturing. The simulations suggest that combining argon ion smoothing with selective, near threshold energy H removal of amorphous C could be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices. Results were found to differ significantly with interatomic potential, and an analysis of these differences was also carried out. Included in this repository are LAMMPS source files, input scripts, and plotting scripts required to reproduce the data. Also included are the output data required to make all the plots included in the associated publication.

Brenner↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2024

Modeling Dense Particle Flow in Multistage and Obstructed Flow Receivers Using High Fidelity Simulations

Particles are a leading contender for next-generation, concentrating solar power technologies, and the design of the particle receiver is critical to minimize the levelized cost of electricity. Falling particle receivers (FPRs) are a viable receiver concept, but many new designs feature complex particle obstructions that include dense discrete phase flows. This creates additional challenges for modeling as particle-to-particle interactions (i.e., collisions) and particle drag become more complex. To improve upon existing modeling strategies, a CFD-DEM simulation capability was created by coupling two independent codes: Sierra/Fuego and LAMMPS. A suitable receiver model was then defined using a traditional continuum-based model for the air and a granular model for the particle curtain. A sensitivity study was executed using this model to determine the relevance of different granular model inputs on important quantities of interest in obstructed flow FPRs: the particle velocity and curtain opacity. The study showed that the granular model inputs had little effect on the particle velocity magnitude and curtain opacity after an obstruction.

Mills, Brantley↗
DOE OSTI2023

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗