Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

Simulating sticky particles: A Monte Carlo method to sample a stratification

Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.

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

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Holistic Numerical Calibration of the iMESA Electrostatic Analyzer

A drifted Maxwellian velocity distribution is the most common model used to interpret the data from low-energy charged-particle instruments onboard spacecraft that are used to investigate the ambient plasma environment in the low Earth orbit (LEO). An original method is presented for determining the flow parameters (density, temperature, and flow energy) of such a distribution from the output of the integrated miniaturized electrostatic analyzer, which has been successfully flown on several LEO missions. Rather than attempting to deconvolve from the on-orbit data the analyzer’s response to an ideal, monoenergetic input, numerical simulation is used to predict and parameterize the response of the device to an input distribution that includes an isotropic, non-zero temperature, yielding a straightforward method for extracting the flow parameters from the spacecraft data. The method is computationally simple enough to be incorporated into a robust algorithm suitable for rapid batch processing or real-time analysis of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Absolute calibration of the Lyman-α measurement apparatus at DIII-D

The LLAMA (Lyman-Alpha Measurement Apparatus) diagnostic was recently installed on the DIII-D tokamak [Rosenthal et al., Rev. Sci. Instrum. (submitted) (2020)]. LLAMA is a pinhole camera system with a narrow band Bragg mirror, a bandpass interference filter, and an absolute extreme ultraviolet photodiode detector array, which measures the Ly-α brightness in the toroidal direction on the inboard, high field side (HFS) and outboard, low field side (LFS). This contribution presents a setup and a procedure for an absolute calibration near the Ly-α line at 121.6 nm. The LLAMA in-vacuum components are designed as a compact, transferable setup that can be mounted in an ex situ vacuum enclosure that is equipped with an absolutely calibrated Ly-α source. The spectral purity and stability of the Ly-α source are characterized using a vacuum ultraviolet spectrometer, while the Ly-α source brightness is measured by a NIST-calibrated photodiode. The non-uniform nature of the Ly-α source emission was overcome by performing a calibration procedure that scans the Ly-α source position and employs a numerical optimization to determine the emission pattern. Nominal and measured calibration factors are determined and compared, showing agreement within their uncertainties. A first conversion of the measured signal obtained from DIII-D indicates that the Ly-α brightness on the HFS and LFS is on the order of 1020 Ph sr-1 m-2 s-1. The established calibration setup and procedure will be regularly used to re-calibrate the LLAMA during DIII-D vents to monitor possible degradation of optical components and detectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A polar direct drive liquid deuterium–tritium wetted foam target concept for inertial confinement fusion

We propose a new approach to inertial confinement fusion (ICF) that could potentially lead to ignition and propagating thermonuclear burn at the National Ignition Facility (NIF). The proposal is based upon a combination of two concepts, referred to as polar direct drive and liquid deuterium–tritium wetted foam capsules. With this new concept, 2D radiation hydrodynamic simulations indicate that ICF ignition and propagating thermonuclear burn are possible with the laser power and energy capabilities available today on the NIF.

3D printing↗
DOE OSTI2022

Hybrid gyrokinetic ion/fluid electron simulation of toroidal tearing modes

The effects of toroidicity and kinetic ions on the resistive tearing mode are systematically studied with the gyrokinetic particle-in-cell simulation code GEM and compared with analytic theory. A new field solver in toroidal geometry has been developed for the simulation of low-n (n = 1, 2) modes in tokamaks. It is found that the toroidal effect significantly reduces the growth rate of the tearing mode. The toroidal effect can also increase the radial width of the tearing mode and change the scaling between the radial mode width and resistivity due to the toroidal pressure term in the electron continuity equation. The kinetic effects of ions can decrease the growth rate of the tearing mode. The plasma flux-surface shaping is found to have significant effect on the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamical electron correlation and the chemical bond. II. Recoupled pair bonds in the $a$ 4 Σ - states of CH and CF

We extended our studies of the effect of dynamical electron correlation on the covalent bonds in the AH and AF series (A = B–F) to the recoupled pair bonds in the excited $a$ 4 Σ - states of the CH and CF molecules. Dynamical correlation is energetically less important in the $a$ 4 Σ - states than in the corresponding X 2 Π states for both molecules, which is reflected in smaller changes in bond energies ($D$ e ). Changes in the equilibrium bond distance ($R$ e ) and vibrational frequency ($ω$ e ), on the other hand, are influenced by the changes in the slope and curvature of the dynamical electron correlation energy as a function of the internuclear distance $R$, $E$ DEC ($R$). In the CH($a$ 4 Σ - ) state, these changes are much smaller than in the CH(X 2 Π) state, but in the CF($a$ 4 Σ - ) state, they are larger, reflecting a significant difference in the shapes of $E$ DEC ($R$) curves. At large $R$, the shape of $E$ DEC ($R$) curves for covalent and recoupled pair bonds is similar although different in magnitude. For the CH($a$ 4 Σ - ) state, $E$ DEC ($R$) has a minimum at $R$ = $R$ e + 0.72 Å as the orbitals associated with the formation of the recoupled pair bond switch places. $E$ DEC (R) for the CF($a$ 4 Σ - ) state decreases continuously throughout the bound region of the potential energy curve because the dynamical correlation energy associated with the electrons in the lone pair orbitals is increasing. In conclusion, these results support our earlier conclusion that we still have much to learn about the nature of dynamical electron correlation in molecules.

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

Flatfielding of hybrid pixel detectors in tender x-ray scattering

The ability of the soft matter interfaces beamline at National Synchrotron Light Source II to access x-ray energy in the tender x-ray regime, i.e., from 2.1 to 5 keV, enables new resonant x-ray scattering studies at the sulfur K-edge and others. We present a new approach to correct data acquired in the tender x-ray regime with a Pilatus3 detector in order to improve the data quality and to correct the various artifacts inherent to hybrid pixel detectors, such as variations in modules’ efficiency or noisy detector module junctions. Further, this new flatfielding significantly enhances the data quality and enables detection of weak scattering signals.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Development of the high energy engineering X-ray (HEX) superconducting wiggler, magnetic measurement, installation, and commissioning

The High energy Engineering X-ray (HEX) diffraction beamline at the National Synchrotron Light Source II (NSLS-II) at Brookhaven National Lab (BNL) is the first high-energy beamline capable of reaching 200 keV for a monochromatic beam. With the 3 GeV electron beam energy for the NSLS-II ring, only the superconducting wiggler (SCW) producing greater than 4 T peak field can cover these ranges with a sufficient number of photons. The 1.2 m-long HEX-SCW has a period length of 70 mm and a field strength on-axis of 4.3 T. It utilizes no liquid helium, and the vertical aperture size of the electron beam vacuum chamber is 8 mm. Unlike regular undulators/wigglers, there is no standard configuration for the magnetic measurement system for superconducting insertion devices. The NSLS-II Insertion Devices group has developed, in collaboration with the vacuum group, a novel in-vacuum Hall mapper with a 1.75 m in-vacuum linear motor and an in-vacuum flip coil system utilizing many commercial-off-the-shelf products. The measurements were conducted at the BNL, and the device was installed in the ring and commissioned. Here, this paper provides a description of the SCW and its magnetic measurement systems, as well as a brief account of the installation and commissioning efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Universal and interpretable classification of atomistic structural transitions via unsupervised graph learning

Materials processing often occurs under extreme dynamic conditions leading to a multitude of unique structural environments. These structural environments generally occur at high temperatures and/or high pressures, often under non-equilibrium conditions, which results in drastic changes in the material's structure over time. Computational techniques, such as molecular dynamics simulations, can probe the atomic regime under these extreme conditions. However, characterizing the resulting diverse atomistic structures as a material undergoes extreme changes in its structure has proved challenging due to the inherently non-linear relationship between structures as large-scale changes occur. Here, we introduce SODAS++, a universal graph neural network framework, that can accurately and intuitively quantify the atomistic structural evolution corresponding to the transition between any two arbitrary phases. We showcase SODAS++ for both solid–solid and solid–liquid transitions for systems of increasing geometric and chemical complexity, such as colloidal systems, elemental Al, rutile and amorphous TiO 2 , and the non-stoichiometric ternary alloy Ag 26 Au 5 Cu 19 . Finally, we show that SODAS++ can accurately quantify all transitions in a physically interpretable manner, showcasing the power of unsupervised graph neural network encodings for capturing the complex and non-linear pathway, a material's structure takes as it evolves.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GradDFT. A software library for machine learning enhanced density functional theory

Density functional theory (DFT) stands as a cornerstone method in computational quantum chemistry and materials science due to its remarkable versatility and scalability. Yet, it suffers from limitations in accuracy, particularly when dealing with strongly correlated systems. To address these shortcomings, recent work has begun to explore how machine learning can expand the capabilities of DFT: an endeavor with many open questions and technical challenges. In this work, we present GradDFT a fully differentiable JAX-based DFT library, enabling quick prototyping and experimentation with machine learning-enhanced exchange–correlation energy functionals. GradDFT employs a pioneering parametrization of exchange–correlation functionals constructed using a weighted sum of energy densities, where the weights are determined using neural networks. Moreover, GradDFT encompasses a comprehensive suite of auxiliary functions, notably featuring a just-in-time compilable and fully differentiable self-consistent iterative procedure. To support training and benchmarking efforts, we additionally compile a curated dataset of experimental dissociation energies of dimers, half of which contain transition metal atoms characterized by strong electronic correlations. The software library is tested against experimental results to study the generalization capabilities of a neural functional across potential energy surfaces and atomic species, as well as the effect of training data noise on the resulting model accuracy.

Chemistry↗
DOE OSTI2024

Enhanced sensitivity to target offset when using cross-beam energy transfer mitigation techniques in direct-drive inertial confinement fusion implosions

In direct-drive inertial confinement fusion, target offset from the target chamber center (or center of beam convergence) may lead to significant implosion asymmetry and fusion yield degradation. In addition, cross-beam energy transfer (CBET) has been shown to be a significant source of laser energy scattering and leads to a reduction in implosion velocity and yield. To improve energy coupling and implosion performance, several techniques for CBET mitigation have been proposed. Recent simulations, however, have shown that CBET also substantially mitigates the effect of target offset on implosion asymmetry and yield [Anderson et al., Phys. Plasmas 27, 112713 (2020)]. Furthermore, the inclusion of CBET models in radiation-hydrodynamics codes was shown to greatly improve agreement between simulations and experiments involving substantial target offset distances. This paper explores the intensity dependence of this CBET–offset effect. In addition, it is shown that enhanced sensitivity to target offset can be expected when CBET-mitigation techniques are used in direct-drive implosions. This is shown through simulations of two such CBET-mitigation techniques on the OMEGA laser: (1) decreased beam-to-target radius, and (2) beam-to-beam frequency detuning. For the typical target offset distances (<15 μm) observed in experiments on OMEGA, however, overall yield is still anticipated to be substantially higher when CBET-mitigation techniques are employed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Imaging nanoscale carrier, thermal, and structural dynamics with time-resolved and ultrafast electron energy-loss spectroscopy

Time-resolved and ultrafast electron energy-loss spectroscopy (EELS) is an emerging technique for measuring photoexcited carriers, lattice dynamics, and near-fields across femtosecond to microsecond timescales. When performed in either a specialized scanning transmission electron microscope or ultrafast electron microscope (UEM), time-resolved and ultrafast EELS can directly image charge carriers, lattice vibrations, and heat dissipation following photoexcitation or applied bias. Yet, recent advances in theoretical calculations and electron optics are often required to realize the full potential of ultrafast EEL spectrum imaging. Here, in this review, we present a comprehensive overview of the recent progress in the theory and instrumentation of time-resolved and ultrafast EELS. We begin with an introduction to the technique, followed by a physical description of the loss function. We outline approaches for calculating and interpreting ground-state and transient EEL spectra spanning low-loss plasmons to core-level excitations analogous to x-ray absorption. We then survey the current state of time-resolved and ultrafast EELS techniques beyond photon-induced near-field electron microscopy, highlighting abilities to image carrier and thermal dynamics. Finally, we examine future directions enabled by emerging technologies, including electron beam monochromation, in situ and operando cells, laser-free UEM, and high-speed direct electron detectors. These advances position time-resolved and ultrafast EELS as a critical tool for uncovering nanoscale dynamic processes in quantum materials and solar energy conversion devices.

Computational methods↗
DOE OSTI2020

Peculiarities of the formation of photonic nanojets by a matrix of dielectric microtoroids

We report the results of theoretical modelling of photonic nanojets (PNJs) formed by laser radiation scattering on a single-layer ordered assembly of dielectric microtoroids placed on the surface of a transparent matrix (silicone film). Using the method of computational electrodynamics (FDTD), the main PNJ parameters (length, width, peak intensity) are analysed under the conditions of mutual influence of the light fields of neighbouring microparticles. It is shown that the main factor affecting the PNJ characteristics under study is the spatial configuration of the radiation-scattering particle, namely, the internal diameter of its cross section. It is found that for certain configurations of toroid placement in a cluster, a PNJ ensemble is implemented with parameters significantly better than those for a single toroid. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗