Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic 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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

End-to-End Differentiable Modeling and Management of the Environment

Focal Area: (2) Data acquisition and assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization. We emphasize the importance of leveraging optimization techniques from AI/machine learning (ML) to solve challenging problems in Earth system modeling. Science Challenge: Automatic differentiation has had a transformative effect on ML by allowing the calculation of gradients of arbitrary functions in an incredibly large class of models. We can potentially realize similar improvements in parameter estimation and control for Earth system models (ESMs) by reimplementing them in computational frameworks from ML. Practitioners working with large (>10 7 parameters) models in ML and AI can obtain good predictive performance in a range of spatiotemporal tasks by making use of optimization via stochastic gradient descent and incorporating prior knowledge at multiple levels. We propose writing ESMs in open-source computational frameworks such as Torch, Tensorflow, and JAX to greatly expand the scope of environmental forecasting and management challenges, which can be addressed by leveraging automatic differentiation and gradient descent-like algorithms. We do not call for a wholesale replacement of physical models with data-driven surrogates, but rather advocate for interleaving physical and empirical equations in a manner that is most faithful to the extent of our scientific knowledge and observational data. Central to this topic is the merging of differentiable physical simulations with differentiable optimization layers, which are now both beginning to come to the forefront.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗
NASA NTRS1989

Computer simulation of dust grain evolution

The latest results are reported from a Monte Carlo code that is being developed at NASA Ames. The goal of this program, is to derive from the observed and presumed properties of the interstellar medium (ISM) the following information: (1) the size spectrum of interstellar dust; (2) the chemical structure of interstellar dust; (3) interstellar abundances; and (4) the lifetime of a dust grain in the ISM. Presently this study is restricted to refractory interstellar material, i.e., the formation and destruction of ices are not included in the program. The program is embedded in an analytic solution for the bulk chemical evolution of a two-phase interstellar medium in which stars are born in molecular clouds, but new nucleosynthesis products and stellar return are entered into a complementary intercloud medium. The well-mixed matter of each interstellar phase is repeatedly cycled stochastically through the complementary phase and back. Refractory dust is created by thermal condensation as stellar matter flows away from sites of nucleosynthesis such as novae and supernovae and/or from the matter returned from evolved intermediate stars. The history of each particle is traced by standard Monte Carlo techniques as it is sputtered and fragmented by supernova shock waves in the intercloud medium. It also accretes an amorphous mantle of gaseous refractory atoms when its local medium joins with the molecular cloud medium. Finally it encounters the possibility of astration (destruction by star formation) within the molecular clouds.

Liffman, K.↗
DOE OSTI2024

Heat Transport Hysteresis Generated Through Frequency Switching of a Time-Dependent Temperature Gradient

A stochastic energetics framework is applied to examine how periodically shifting the frequency of a time-dependent oscillating temperature gradient affects heat transport in a nanoscale molecular model. We specifically examine the effects that frequency switching, i.e., instantaneously changing the oscillation frequency of the temperature gradient, has on the shape of the heat transport hysteresis curves generated by a particle connected to two thermal baths, each with a temperature that is oscillating in time. Analytical expressions are derived for the energy fluxes in/out of the system and the baths, with excellent agreement observed between the analytical expressions and the results from nonequilibrium molecular dynamics simulations. We find that the shape of the heat transport hysteresis curves can be significantly altered by shifting the frequency between fast and slow oscillation regimes. We also observe the emergence of features in the hysteresis curves such as pinched loops and complex multi-loop patterns due to the frequency shifting. The presented results have implications in the design of thermal neuromorphic devices such as thermal memristors and thermal memcapacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Surrogate modeling for efficiently, accurately and conservatively estimating measures of risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. Risk measures combine objective probabilities with the subjective values of a decision maker to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical approximations obtained from the training data. These surrogate models limit over-confidence in reliability and safety assessments and produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to a stakeholder’s risk preferences and then present an approach, based on stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. Our surrogate models include biases that permit them to conservatively estimate the target risk measures. We provide theoretical results that show that these biases decay at the same rate as the L 2 error in the surrogate model. Numerical demonstrations confirm that risk-adapted surrogate models do indeed overestimate the target risk measures while converging at the expected rate.

42 ENGINEERING↗
DOE OSTI2022

Measuring characteristic differences between high- and low-performing discharges on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF

A dense plasma focus (DPF) is a compact coaxial plasma gun, which completes its discharge as a Z-pinch, producing short (<100 ns) pulses of ions, x rays, and/or neutrons. Lawrence Livermore National Laboratory recently constructed and began operating a new device, the MJOLNIR (MegaJOuLe Neutron Imaging Radiography) DPF, which is designed for single-pulse flash neutron radiography. This device has achieved neutron yields of up to 4.1×10 11 neutrons/pulse at 3.3 MA peak current, and higher-current commissioning is under way. Like most DPFs, MJOLNIR exhibits variable yields in some configurations. Here, we present evidence of the role of parasitic current paths within the gun in stochastically influencing the yield. First through “conditioning shots,” where new hardware has been introduced, we show that increased run-down and run-in speeds correlate with higher yields. These observations are consistent with current being delivered to the electrodes but not to the main plasma sheath, degrading the implosion-driving force. Once nominal conditions are established, we correlate low-performing discharges with smaller current dip and associated voltage spike for a fixed machine configuration. A snow-plow model is able to recreate small-magnitude current dips through the introduction of a parasitic current path, and particle-in-cell simulations establish how parasitic current paths lower the ion beam energy available to produce neutrons. Finally, we observe an increased likelihood of shots with low yield and smaller current dip with increasing fill pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Allowing for Horizontally Heterogeneous Clouds and Generalized Overlap in an Atmospheric GCM

While fully accounting for 3D effects in Global Climate Models (GCMs) appears not realistic at the present time for a variety of reasons such as computational cost and unavailability of 3D cloud structure in the models, incorporation in radiation schemes of subgrid cloud variability described by one-point statistics is now considered feasible and is being actively pursued. This development has gained momentum once it was demonstrated that CPU-intensive spectrally explicit Independent Column Approximation (lCA) can be substituted by stochastic Monte Carlo ICA (McICA) calculations where spectral integration is accomplished in a manner that produces relatively benign random noise. The McICA approach has been implemented in Goddard's GEOS-5 atmospheric GCM as part of the implementation of the RRTMG radiation package. GEOS-5 with McICA and RRTMG can handle horizontally variable clouds which can be set via a cloud generator to arbitrarily overlap within the full spectrum of maximum and random both in terms of cloud fraction and layer condensate distributions. In our presentation we will show radiative and other impacts of the combined horizontal and vertical cloud variability on multi-year simulations of an otherwise untuned GEOS-5 with fixed SSTs. Introducing cloud horizontal heterogeneity without changing the mean amounts of condensate reduces reflected solar and increases thermal radiation to space, but disproportionate changes may increase the radiative imbalance at TOA. The net radiation at TOA can be modulated by allowing the parameters of the generalized overlap and heterogeneity scheme to vary, a dependence whose behavior we will discuss. The sensitivity of the cloud radiative forcing to the parameters of cloud horizontal heterogeneity and comparisons of CERES-derived forcing will be shown.

Lee, D.↗
DOE OSTI2021

Electrodeposition of ligand-free copper nanoparticles from aqueous nanodroplets

Copper nanoparticles have emerged as promising electrocatalysts for energy storage and conversion. Generally, homogenous nanoparticle synthesis requires a stabilizing ligand, which may influence the electrocatalysis. Ligands can be avoided by direct nanoparticle electrodeposition. Here, we extend nanodroplet-mediated electrodeposition to the electrodeposition of copper nanoparticles from aqueous nanodroplets suspended in 0.1 M tetrabutylammonium hexafluorophosphate ([TBA][PF 6 ]) and 1,2-dichloroethane. Stochastic electrochemistry on microelectrodes was used to elucidate nanoparticle growth kinetics for various solution conditions and substrate materials. A study on the effect of surfactants on nanoparticle morphology and growth kinetics demonstrated the surfactant's role as an avenue for morphological control. Nanoparticle morphology was studied by Scanning Electron Microscopy (SEM) and Transmission Electron Microscopy (TEM). Energy-dispersive X-ray Spectroscopy (EDX) confirmed the presence of copper nanoparticles. TEM and Selected Area Electron Diffraction (SAED) confirmed nanoparticle size, morphology, and polycrystallinity. We demonstrate tunable nanoparticle size using cyclic voltammetry on single aqueous nanodroplets by altering the voltammetric sweeping potentials of copper electroreduction. Finite element simulations validate the voltammetric response and ability to control the nanoparticle size with nanometer resolution during a voltammetric sweep. These results inform the electrodeposition of copper nanoparticles from aqueous nanodroplets for important applications, such as the conversion of CO 2 to hydrocarbon fuels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reaction dynamics and in-medium nucleon-nucleon cross section with 12 C+ 1 H at 95 MeV/nucleon

Reaction dynamics and effective in-medium nucleon-nucleon (NN) cross section are studied, using the antisymmetrized molecular dynamics (AMD) model with one of the simplest heavy ion collisions, 12 C+ ⁢1 H. The roles of different stochastic processes are studied, which are installed in the original AMD. One is the diffusion process as a part of the quantum fluctuation originated from the Fermi motion during the time evolution of the wave packets in the mean field (AMD/D). The second is a coalescence treatment for light clusters with A ≤ 8 (AMD/D-COALS). The third is the collision process, installed in AMD with different form of the in-medium NN cross sections. These different processes are examined, using the experimental angular distributions and energy spectra of light charged particles and intermediate mass fragments. The effective in-medium NN cross sections are evaluated from the simulated AMD events and compared with the experimentally extracted values. A good agreement is obtained between the AMD values from the 12 C + 1 H reaction at 95 MeV/nucleon and the experimental results for nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.

Zhang, Xin↗
DOE OSTI2021

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗
DOE OSTI2021

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗
DOE OSTI2022

A data science approach for analysis and reconstruction of spinodal-like composition fields in irradiated FeCrAl alloys

A statistical method for the analysis of continuously distributed data representative of composition fluctuations in irradiated FeCrAl alloys acquired using Energy Dispersive X-ray Spectroscopy (EDS) method is presented. Using probability distribution functions, direct and cross-covariances between the elemental compositions, the effects of alloy composition and irradiation dose were investigated on the spatial distribution and length scale of composition fluctuations at the nanoscale. We have observed that, for neutron-irradiated FeCrAl alloys, the distribution of Fe and Cr followed a left-skewed and right-skewed distribution, respectively for all (average) alloy compositions and irradiation doses. The analysis also revealed enhanced spatial gradients in the elemental compositions at higher irradiation dose. Direct and cross-covariance estimates of the experimental data were also utilized for reconstruction of composition data through fitting it to a parametric form of the covariance functions. Linear Model of Coregionalization was used to determine the parameters of the covariance functions. Subsequently, a spectral method was utilized for simulating a realization of the alloy compositions. Close correspondence was observed between the experimental and the reconstructed data which was analyzed using probability distribution functions and covariance functions. Composition space of the experimental and reconstructed data and dislocation velocities as a function of applied stress and line directions over the entire composition maps were also examined.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf galaxy evolution

ABSTRACT We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR), respectively, to study 78 radio active galactic nuclei (AGN) in nearby (z < 0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages, and accretion rates are ∼1035 W, ∼5 Myr, and ∼10−3.4 M⊙ yr−1, respectively. The median mechanical energy output of these sources is ∼100 times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability, and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf galaxy evolution, as is the case in massive galaxies.

Davis, F. (ORCID:0000000285224980)↗