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

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A comparison of machine learning surrogate models of street-scale flooding in Norfolk, Virginia

Low-lying coastal cities, exemplified by Norfolk, Virginia, face the challenge of street flooding caused by rainfall and tides, which strain transportation and sewer systems and can lead to personal and property damage. While high-fidelity, physics-based simulations provide accurate predictions of urban pluvial flooding, their computational complexity renders them unsuitable for real-time applications. Using data from Norfolk rainfall events between 2016 and 2018, this study compares the performance of a previous surrogate model based on a random forest algorithm with two deep learning models: Long Short-Term Memory (LSTM) and Gated Recurrent Unit (GRU). The comparison of deep learning to the random forest algorithm is motivated by the desire to utilize a machine learning architecture that allows for the future inclusion of common uncertainty quantification techniques and the effective integration of relevant, multi-modal features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling and implementation of vertical excursion FFA in the Zgoubi ray-tracing code

Vertical Fixed Field Accelerators (vFFAs) feature complex and highly non-linear magnetic fields, which require simulation codes allowing step-wise particle tracking. Methods to model the 3D magnetic field of scaling vFFAs have been developed in the ray-tracing code Zgoubi. The field modeling and particle tracking methods include the field non-linearities, the fringe fields, and the field superposition of neighboring magnets. The procedure implements the vFFA analytical field expressions, allowing design studies and parameter optimizations using the Zgoubi built-in fit method. The vFFA procedure has been applied to a ten-fold symmetry ring with a triplet focusing structure designed to accelerate protons from 3 MeV to 12 MeV and studied under the ISIS-II proton driver prototype project. Results from particle tracking in externally generated 3D semi-analytical field maps and the developed vFFA analytical model are shown to be in excellent agreement.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

The dispersion and propagation of topological Langmuir-cyclotron waves in cold magnetized plasmas

The topological Langmuir-cyclotron wave (TLCW) is a recently identified topological surface excitation in magnetized plasmas. We show that TLCW originates from the topological phase transition at the Langmuir wave-cyclotron wave resonance. By isofrequency surface analysis and two- and three-dimensional time-dependent simulations, we demonstrate that the TLCW can propagate robustly along complex phase transition interfaces in a unidirectional manner and without scattering. Because of these desirable features, the TLCW could be explored as an effective mechanism to drive current and flow in magnetized plasmas. As a result, the analysis also establishes a close connection between the newly instituted topological phase classification of plasmas and the classical Clemmow-Mullaly-Allis (CMA) diagram of plasma waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Designing accurate emulators for scientific processes using calibration-driven deep models

Abstract Predictive models that accurately emulate complex scientific processes can achieve speed-ups over numerical simulators or experiments and at the same time provide surrogates for improving the subsequent analysis. Consequently, there is a recent surge in utilizing modern machine learning methods to build data-driven emulators. In this work, we study an often overlooked, yet important, problem of choosing loss functions while designing such emulators. Popular choices such as the mean squared error or the mean absolute error are based on a symmetric noise assumption and can be unsuitable for heterogeneous data or asymmetric noise distributions. We propose Learn-by-Calibrating, a novel deep learning approach based on interval calibration for designing emulators that can effectively recover the inherent noise structure without any explicit priors. Using a large suite of use-cases, we demonstrate the efficacy of our approach in providing high-quality emulators, when compared to widely-adopted loss function choices, even in small-data regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Formation of secondary organic aerosol from wildfire emissions enhanced by long-time ageing

Wildfire smoke, consisting primarily of organic aerosols, has profound impacts on air quality, climate and human health. Wildfire organic aerosol evolves over long-time photochemical oxidation due to the formation and ageing of secondary organic aerosol, which substantially changes its magnitude and properties. However, there are large uncertainties in the long-time ageing of wildfire organic aerosol because of the distinct ageing behaviours of the complex organic emissions. Here we developed an oxidation model that simulates the ageing of wildfire organic emissions in the full volatility range on a precursor level and integrated insights from single-species ageing and wildfire emissions ageing experiments and field plume observations to constrain the long-time ageing of wildfire organic aerosol. The model captured the enhancement of organic aerosol mass (2–8 times) and oxygen-to-carbon ratio (1–4 times) in the wildfire ageing experiments. It also reconciled a long-standing discrepancy between field and laboratory observations of the magnitude of secondary organic aerosol formation. The model indicated large emissions-driven variations in precursor contributions to secondary organic aerosol, which further evolve with long-time ageing. In conclusion, the estimated global wildfire secondary organic aerosol production (139 ± 34 Tg per year) was much higher than previous studies omitting or under-constraining long-time ageing.

58 GEOSCIENCES↗
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Feasibility study of a compact neutron resonance transmission analysis instrument

Neutron Resonance Transmission Analysis (NRTA) uses resonant absorption of neutrons to infer the absolute isotopic composition of a target object, enabling applications in a broad range of fields such as archaeology, enrichment analysis of nuclear fuel, and arms control treaty verification. In the past, NRTA involved large user facilities and complex detector systems. However, recent advances in the intensity of compact neutron sources have made compact neutron imaging designs increasingly feasible. This work describes the Monte Carlo (MC) based design of a compact epithermal NRTA radiographic instrument which uses a moderated, compact deuterium-tritium (DT) neutron source and an epithermal neutron detector. Such an instrument would have a wide range of applications, and would be especially impactful for such scenarios as nuclear inspection and arms control verification exercises, where system complexity and mobility may be of critical importance. The MC simulations presented in this work demonstrate accurate time-of-flight (TOF) reconstructions for transmitted neutron energies, capable of differentiating isotopic compositions of nuclear material with high levels of accuracy. A new generation of miniaturized and increasingly more intense neutron sources will allow this technique to achieve measurements with greater precision and speed, with significant impact on a variety of engineering and societal problems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Competitive formation of intercalated versus supported metal nanoclusters during deposition on layered materials with surface point defects

Intercalated metal nanoclusters (NCs) can be formed under the surface of graphite after sputtering to generate surface “portal” defects that allow deposited atoms to reach the subsurface gallery. Although, there is a competition between formation of supported NCs on top of the surface and intercalated NCs under the surface, the latter only dominating at sufficiently high temperature. A stochastic model incorporating appropriate system thermodynamics and kinetics is developed to capture this complex and competitive nucleation and growth process. Kinetic Monte Carlo simulation shows that the model captures experimental trends observed for Cu and other metals and reveals that higher temperatures are needed to facilitate detachment of atoms from supported NCs enabling them to reach the gallery.

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

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

A generalized approach to x-ray data modeling for high-energy-density plasma experiments

Accurate understanding of x-ray diagnostics is crucial for both interpreting high-energy-density experiments and testing simulations through quantitative comparisons. X-ray diagnostic models are complex. Past treatments of individual x-ray diagnostics on a case-by-case basis have hindered universal diagnostic understanding. Here, in this study, we derive a general formula for modeling the absolute response of non-focusing x-ray diagnostics, such as x-ray imagers, one-dimensional space-resolved spectrometers, and x-ray power diagnostics. The present model is useful for both data modeling and data processing. It naturally accounts for the x-ray crystal broadening. The new model verifies that standard approaches for a crystal response can be good approximations, but they can underestimate the total reflectivity and overestimate spectral resolving power by more than a factor of 2 in some cases near reflectivity edge features. We also find that a frequently used, simplified-crystal-response approximation for processing spectral data can introduce an absolute error of more than an order of magnitude and the relative spectral radiance error of a factor of 3. The present model is derived with straightforward geometric arguments. It is more general and is recommended for developing a unified picture and providing consistent treatment over multiple x-ray diagnostics. Such consistency is crucial for reliable multi-objective data analyses.

Nagayama, Taisuke↗
DOE OSTI2024

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. Here, we illustrate PigmentHunter’s capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna–Matthews–Olson protein, and ring deformation in photosystems I and II.

14 SOLAR ENERGY↗
DOE OSTI2022

Opinion dynamics in financial markets via random networks

We investigate financial market dynamics by introducing a heterogeneous agent-based opinion formation model. In this work, we organize individuals in a financial market according to their trading strategy, namely, whether they are noise traders or fundamentalists. The opinion of a local majority compels the market exchanging behavior of noise traders, whereas the global behavior of the market influences the decisions of fundamentalist agents. We introduce a noise parameter, q , to represent the level of anxiety and perceived uncertainty regarding market behavior, enabling the possibility of adrift financial action. We place individuals as nodes in an Erdös-Rényi random graph, where the links represent their social interactions. At any given time, individuals assume one of two possible opinion states ±1 regarding buying or selling an asset. The model exhibits fundamental qualitative and quantitative real-world market features such as the distribution of logarithmic returns with fat tails, clustered volatility, and the long-term correlation of returns. We use Student’s t distributions to fit the histograms of logarithmic returns, showing a gradual shift from a leptokurtic to a mesokurtic regime depending on the fraction of fundamentalist agents. Furthermore, we compare our results with those concerning the distribution of the logarithmic returns of several real-world financial indices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Understanding the electromagnetic topology during the ohmic breakdown in tokamaks considering self-generated electric fields

We report that the physical mechanisms of the ohmic breakdown in a tokamak have been understood based on the classical Townsend avalanche theory. However, a new systematic theory (Yoo et al 2018 Nat. Commun.9 3523) recently demonstrated that electron avalanches during the ohmic breakdown are completely different from the Townsend avalanche due to strong self-generated electric fields. In this study, we elucidate the multi-dimensional effects of the self-generated electric field on plasma dynamics during the ohmic breakdown. We also propose a novel electromagnetic topology analysis method that can easily predict the overall plasma behavior and where the main plasma is generated. The topology analysis method is validated by a state-of-art particle simulation for various magnetic configurations. New physical insights into the complex electromagnetic topology would facilitate designing more reliable and optimized ohmic breakdown scenarios in future tokamaks, such as ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Preliminary Introduction of a Free Vortex Wake Method Into OpenFAST

Increasing wind turbine size can have substantial economic benefits. However, increasing wind turbine turbine size requires the use of models that accurately capture the aerodynamic response to large deflections and nonstraight blade geometries. Additionally, the installation of offshore floating turbines requires turbine models that can handle large rotor motion. Free vortex wake (FVW) methods are capable of modeling this type of complex physics while remaining computationally tractable to perform the many simulations necessary for turbine design. In this work, a FVW model is added to the National Renewable Energy Laboratory engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. Code implementation is not complete and thus no results are presented. However, areas of potential modeling improvement are discussed, such as modifications to the circulation calculation and the need for different FVW input model specifications when high-thrust cases are being modeled.

17 WIND ENERGY↗
DOE OSTI2022

Copper migration and surface oxidation of CuxBi2Se3 in ambient pressure environments

Abstract Chemical modifications such as intercalation can be used to modify surface properties or to further functionalize the surface states of topological insulators (TIs). Using ambient pressure x-ray photoelectron spectroscopy, we report copper migration in C u x B i 2 S e 3 , which occurs on a timescale of hours to days after initial surface cleaving. The increase in near-surface copper proceeds along with the oxidation of the sample surface and large changes in the selenium content. These complex changes are further modeled with core-level spectroscopy simulations, which suggest a composition gradient near the surface which develops with oxygen exposure. Our results shed light on a new phenomenon that must be considered for intercalated TIs—and intercalated materials in general—that surface chemical composition can change when specimens are exposed to ambient conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Smooth periodic gauge satisfying crystal symmetry and periodicity to study high-harmonic generation in solids

Intense lasers can easily drive nonadiabatic transitions of excited electron wave packets across the Brillouin zones, thus transition dipole moments (TDM) between energy bands of solids should be continuous, satisfying crystal symmetry, and periodic at zone boundaries. While current ab initio algorithms are powerful in calculating band structures of solids, they all introduced random phases into the eigenfunctions at each crystal momentum k. In this work, we show how to choose a “smooth-periodic” gauge where TDMs can be smooth versus k, preserving crystal symmetry, as well as maintaining periodic at boundaries. The symmetry properties of TDMs with respect to k ensure the absence of even-order harmonics from MgO with inversion symmetry, while the TDM in the “smooth-periodic” gauge for broken-symmetry ZnO is responsible for even harmonics that were underestimated in previous simulations. These results reveal the importance of correctly treating the complex TDMs that satisfy crystal symmetry and continuous across zone boundaries in nonlinear laser-solid interactions, which has been elusive in most theories so far.

36 MATERIALS SCIENCE↗