Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation framework”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Black carbon mixing state impacts on cloud microphysical properties: Effects of aerosol plume and environmental conditions

Black carbon (BC) is usually mixed with other aerosol species within individual aerosol particles. This mixture, along with the particles' size and morphology, determines the particles' optical and cloud condensation nuclei properties, and hence black carbon's climate impacts. In this study the particle-resolved aerosol model PartMC-MOSAIC was used to quantify the importance of black carbon mixing state for predicting cloud microphysical quantities. Based on a set of about 100 cloud parcel simulations a process level analysis framework was developed to attribute the response in cloud microphysical properties to changes in the underlying aerosol population ("plume effect") and the cloud parcel cooling rate ("parcel effect"). It shows that the response of cloud droplet number concentration to changes in BC emissions depends on the BC mixing state. When the aerosol population contains mainly aged BC particles an increase in BC emission results in increasing cloud droplet number concentrations ("additive effect"). In contrast, when the aerosol population contains mainly fresh BC particles they act as sinks for condensable gaseous species, resulting in a decrease in cloud droplet number concentration as BC emissions are increased ("competition effect"). Additionally, we quantified the error in cloud microphysical quantities when neglecting the information on BC mixing state, which is often done in aerosol models. Here, the errors ranged from -12% to +45% for the cloud droplet number fraction, from 0% to +1022% for the nucleation-scavenged black carbon (BC) mass fraction, from -12% to +4% for the effective radius, and from -30% to +60% for the relative dispersion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2023

How Climate and Data Quality Impact Photovoltaic Performance Loss Rate Estimations

Different data pipelines and statistical methods are applied to photovoltaic (PV) performance datasets to quantify the performance loss rate (PLR). Since the real values of PLR are unknown, a variety of unvalidated values are reported. As such, the PV industry commonly assumes PLR based on statistically extracted ranges from the literature. However, the accuracy and uncertainty of PLR depend on several parameters including seasonality, local climatic conditions, and the response of a particular PV technology. In addition, the specific data pipeline and statistical method used affect the accuracy and uncertainty. To provide insights, a framework of (≈200 million) synthetic simulations of PV performance datasets using data from different climates is developed. Time series with known PLR and data quality are synthesized, and large parametric studies are conducted to examine the accuracy and uncertainty of different statistical approaches over the contiguous US, with an emphasis on the publicly available and “standardized” library, RdTools . In the results, it is confirmed that PLRs from RdTools are unbiased on average, but the accuracy and uncertainty of individual PLR estimates vary with climate zone, data quality, PV technology, and choice of analysis workflow. Best practices and improvement recommendations based on the findings of this study are provided.

14 SOLAR ENERGY↗
DOE OSTI2025

End-to-end microgrid protection using distributed data-driven methods

This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Informing Geologic CO 2 Storage Site Management Decisions under Uncertainty: Demonstration of NRAP's Integrated Assessment Model (NRAP-IAM-CS) Application

Managing risks at geologic CO 2 storage sites requires making decisions based on the quantitative predictions of long-term performance that not only take into account various subsurface hydrologic, geochemical processes and interactions but also the related uncertainties. The US DOE (United States Department of Energy) funded National Risk Assessment Partnership (NRAP) has developed a quantitative approach and a model to predict the long-term behavior of GCS sites. The approach uses an integrated assessment model (NRAP-IAM-CS) that couples reduced order models (ROMs) for various GCS components (such as storage reservoir, wells, seals, and groundwater aquifer) in a stochastic modeling framework that allows for rapid simulations of entire system behavior over thousands of years while probing uncertainties. We demonstrate application of NRAP-IAM-CS using two examples that focus on CO 2 leakage and resulting impacts. The examples are used to demonstrate how different uncertain parameters affect overall performance objectives.

03 NATURAL GAS↗
DOE OSTI2018

Mass, momentum, and energy transfer in supersonic aerosol deposition processes

Aerosol deposition (AD) can lead to the formation of dense coatings via deposition of particles from a gas flow; in AD the aerosol is passed through a converging-diverging nozzle, facilitating inertial particle impaction on a desired substrate at supersonic particle velocities. Unique from thermal spray methods, AD can be applied near room temperature and unique from cold spray, in AD the aerosol is typically at atmospheric pressure upstream of the nozzle. Though AD has been successfully demonstrated previously, a number of aspects related to particle motion in AD systems remain poorly understood. In this work, we simulated compressible flow field profiles and particle trajectories for typical AD working conditions for a slit type converging-diverging nozzle with a planar substrate. In examining the fluid flow profile, we show that the velocity and pressure profiles, as well as the shock structure are sensitive to the upstream and downstream operating pressures of the nozzle. These ultimately affect particle impaction speed. Importantly, in AD, the particle drag regime is dynamic; both particle Knudsen numbers and Mach numbers can vary by orders of magnitude. To aid particle trajectory simulations, we trained a neural network to predict the drag force on the particles based on existing experimental data, theoretical limits, and new direct simulation Monte Carlo (DMSC) results. The neural network based drag law, which depends upon both Mach and Knudsen numbers, shows better agreement with the DSMC simulation data than pre-existing correlations. With it, particle trajectory simulation results reveal that for a given particle density, there exists an optimal particle diameter to maximize particle impaction speed. We also find that in AD particles undergo size dependent inertial focusing, i.e. there is a particular particle diameter where the particle deposition linewidth is minimized. Particles smaller than this diameter are underfocused, and particles larger than this are overfocused, and hence have larger deposition linewidths in both cases. Using trajectory simulations, we additionally develop a framework that can be used to evaluate the position-dependent mass, momentum and kinetic energy fluxes to the deposition substrate for any aerosol size distribution function upstream the nozzle. It is shown for typical aerosol concentrations achievable in the laboratory, the kinetic energy flux can approach a magnitude normally observed in convective heat transfer with phase change, hence translational kinetic energy to thermal energy transfer in AD is likely the key contributor to the formation of dense coatings.

42 ENGINEERING↗
DOE OSTI2026

Enabling event-by-event precision in γ-ray cascades for neutron-induced reactions

Neutron-induced γ-ray spectra provide key inputs for modern active interrogation applications. A precise modeling of the nuclear reaction and subsequent emission of γ rays is challenging and often impossible due to limitations on evaluated data file formats and nuclear transport simulation codes. We present a framework that addresses these challenges by combining experimental data and reaction-model calculation outputs into an extended candidate version of the Generalized Nuclear Data Structure (GNDS) file, the successor format for the legacy Evaluated Nuclear Data File (ENDF-6). This proposed GNDS hierarchical format contains all the necessary ingredients for inline γ-ray cascade reproduction with event-by-event precision, including continuum–continuum and continuum–discrete transitions following neutron-capture and inelastic neutron scattering reactions. Cascade-event generation based on our approach demonstrates improved energy conservation on an event-by-event basis and permits the use of γ-γ coincidences in applications. This work offers, for the first time, a method to generate neutron-capture and inelastic neutron-scattering γ-ray cascades where energy conservation, correlations, and experimental primaries are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Plasma-assisted chemical-looping combustion: Mechanistic insights into low temperature methane oxidation with CuO

The low-temperature oxidation of CH 4 by CuO in a coaxial, fixed bed, double dielectric barrier discharge (DBD) reactor was investigated with time-dependent species measurements by an electron-ionization molecular beam mass spectrometer (EI-MBMS). In the experiment, 10% methane carried by noble gasses was flown at 50 sccm through 1 g CuO dispersed in quartz wool both under plasma and non-plasma conditions, while time-dependent gas-phase species profiles were collected. Plasma conditions were explored from 300 to 600 °C while the non-plasma conditions were set from 300 to 900 °C. Mechanistic insights into the oxidation of CH 4 by CuO with plasma discharge at lower temperatures (≤ 600 °C) were obtained by quantifying the fuel oxidation, intermediate species, and CO 2 production in comparison to the non-plasma conditions. Here, we observed significant enhancement of fuel oxidation from the plasma discharge between 400 and 500 °C. The CO 2 production at 500 °C with plasma was greater than that at 700 °C without plasma, reducing fuel oxidation temperature by 200+ °C. During tests, three distinct reaction stages were observed: a gas-phase transport limited stage, a surface reaction limited kinetic stage, and an oxygen ion diffusion limited stage. It was observed that plasma greatly improved the reactivity of the second stage at low temperature. In addition, no carbon deposits were observed on the resultant particles, even under the presence of plasma. M. species such as C 4 H 2 and C 6 H 6 not previously observed or predicted in CuO/CH 4 chemical looping were observed, with some species such as CH 3 OH only becoming detectable as total flowrate was increased from 50 to 1500sccm. A non-plasma reaction pathway for CH 4 based the observed species from the MBMS spectrum and previous predictions from reactive molecular dynamics simulations was created, providing a framework from which more complex plasma CuO mechanisms can be crafted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluating the Role of Metastable Surfaces in Mechanochemical Reduction of Molybdenum Oxide

Mechanochemistry and mechanocatalysis are gaining increasing attention as environmentally friendly chemical processes because of their solvent-free nature and scalability. Significant effort has been devoted for studying continuum-scale phenomena in mechanochemistry, such as temperature and pressure gradients, but the atomic-scale mechanisms remain relatively unexplored. In this work, we focus on the mechanochemical reduction of MoO3 as a case study. We use experimental techniques to determine the mechanochemical reduction conditions and density functional theory (DFT) simulations to establish an atomistic framework for identifying the metastable surfaces that are most likely to enable this process. Our results show that metastable surfaces can significantly lower or remove thermodynamic barriers for surface reduction and that kinetic energy from milling can facilitate the formation of metastable surfaces that have high surface fracture energies and are not thermally accessible. These findings indicate that metastable surfaces are an important aspect of mechanochemistry along with hot spots and other continuum-scale phenomena.

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

Unified theoretical framework for black carbon mixing state allows greater accuracy of climate effect estimation

Black carbon (BC) plays an important role in the climate system because of its strong warming effect, yet the magnitude of this effect is highly uncertain owing to the complex mixing state of aerosols. Here we build a unified theoretical framework to describe BC’s mixing states, linking dynamic processes to BC coating thickness distribution, and show its self-similarity for sites in diverse environments. The size distribution of BC-containing particles is found to follow a universal law and is independent of BC core size. A new mixing state module is established based on this finding and successfully applied in global and regional models, which increases the accuracy of aerosol climate effect estimations. Our theoretical framework links observations with model simulations in both mixing state description and light absorption quantification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unravelling crystal growth of nanoparticles

Crystal growth from nanoscale constituents is a ubiquitous phenomenon in biology, geology and materials science. Numerous studies have focused on understanding the onset of nucleation and on producing high-quality crystals by empirically sampling constituents with different attributes and varying the growth conditions. However, the kinetics of post-nucleation growth processes, an important determinant of crystal morphology and properties, have remained underexplored due to experimental challenges associated with real-space imaging at the nanoscale. Here we report the imaging of the crystal growth of nanoparticles of different shapes using liquid-phase transmission electron microscopy, resolving both lateral and perpendicular growth of crystal layers by tracking individual nanoparticles. We observe that these nanoscale systems exhibit layer-by-layer growth, typical of atomic crystallization, as well as rough growth prevalent in colloidal systems. Surprisingly, the lateral and perpendicular growth modes can be independently controlled, resulting in two mixed crystallization modes that, until now, have received only scant attention. Combining analytical considerations with molecular dynamics and kinetic Monte Carlo simulations, we develop a comprehensive framework for our observations, which are fundamentally determined by the size and shape of the building blocks. Furthermore, these insights unify the understanding of crystal growth across four orders of magnitude in particle size and suggest novel pathways to crystal engineering.

Imaging techniques↗
DOE OSTI2014

Dynamic mesoscale model of dipolar fluids via fluctuating hydrodynamics

Fluctuating hydrodynamics (FHD) is a general framework of mesoscopic modeling and simulation based on conservational laws and constitutive equations of linear and nonlinear responses. However, explicit representation of electrical forces in FHD has yet to appear. In this work, we devised an Ansatz for the dynamics of dipole moment densities that is linked with the Poisson equation of the electrical potential ϕ in coupling to the other equations of FHD. The resulting ϕ-FHD equations then serve as a platform for integrating the essential forces, including electrostatics in addition to hydrodynamics, pressure-volume equation of state, surface tension, and solvent-particle interactions that govern the emergent behaviors of molecular systems at an intermediate scale. This unique merit of ϕ-FHD is illustrated by showing that the water dielectric function and ion hydration free energies in homogeneous and heterogenous systems can be captured accurately via the mesoscopic simulation. Furthermore, we show that the field variables of ϕ-FHD can be mapped from the trajectory of an all-atom molecular dynamics simulation such that model development and parametrization can be based on the information obtained at a finer-grained scale. With the aforementioned multiscale capabilities and a spatial resolution as high as 5 Å, the ϕ-FHD equations represent a useful semi-explicit solvent model for the modeling and simulation of complex systems, such as biomolecular machines and nanofluidics.

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

Contact stiffness and damping of liquid films in dynamic atomic force microscope

The mechanical properties and dissipation behaviors of nanometers confined liquid films have been long-standing interests in surface force measurements. The correlation between the contact stiffness and damping of the nanoconfined film is still not well understood. We establish a novel computational framework through molecular dynamics (MD) simulation for the first time to study small-amplitude dynamic atomic force microscopy (dynamic AFM) in a simple nonpolar liquid. Through introducing a tip driven dynamics to mimic the mechanical oscillations of the dynamic AFM tip-cantilever assembly, we find that the contact stiffness and damping of the confined film exhibit distinct oscillations within 6-7 monolayer distances, and they are generally out-of-phase. For the solid-like film with integer monolayer thickness, further compression of the film before layering transition leads to higher stiffness and lower damping, while much lower stiffness and higher damping occur at non-integer monolayer distances. These two alternating mechanisms dominate the mechanical properties and dissipation behaviors of simple liquid films under cyclic elastic compression and inelastic squeeze-out. Our MD simulations provide a direct picture of correlations between the structural property, mechanical stiffness, and dissipation behavior of the nanoconfined film.

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

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
DOE OSTI2025

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium↗