Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Radiative Forcing of Nitrate Aerosols From 1975 to 2010 as Simulated by MOSAIC Module in CESM2-MAM4

In this work, we incorporate the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) module in the Community Earth System Model version 2 (CESM2) Community Atmosphere Model Version 6 with interactive chemistry (CAM6-chem), and couple it with the four mode version of the Modal Aerosol Module (MAM4). The MOSAIC module is used to simulate the thermodynamics of the gas-aerosol mass exchange, with a special focus on simulating nitrate aerosol. By comparing against ground and satellite observations, we found that the MOSAIC/MAM4 scheme performs reasonably well in simulating spatiotemporal distributions of aerosols, including nitrate aerosol. We conducted a series of model experiments with and without nitrate aerosols, and examined the radiative effect (RE) associated with nitrate aerosols in 1975, 2000, and 2010, and accessed the radiative forcing (RF) of nitrate aerosols between the present day and pre-industrial periods. Comparing with the nitrate aerosol RE, we predicted relatively small RF of anthropogenic nitrate aerosol from aerosol-radiation interactions (RF ari : -0.014 W m -2 ) and large RF from aerosol-cloud interactions (RF aci : -0.219 W m -2 ). Regional signatures of nitrate RE/RF are noticeable and important: for instance, very small changes in RE ari in Europe and USA, but 2.8–3 times increases in RE ari in India and China from 1975 to 2010, while RE aci /RF aci in China is a warming effect due to the competing effect between sulfate and nitrate aerosols as cloud condensation nuclei.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global Dust Cycle and Direct Radiative Effect in E3SM Version 1: Impact of Increasing Model Resolution

Quantification of dust aerosols in Earth System Models (ESMs) has important implications for water cycle and biogeochemistry studies. This study examines the global life cycle and direct radiative effects (DREs) of dust in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1), and the impact of increasing model resolution both horizontally and vertically. The default 1° E3SMv1 captures the spatial and temporal variability in the observed dust aerosol optical depth (DAOD) reasonably well, but overpredicts dust absorption in the shortwave (SW). Simulations underestimate the dust vertical and long-range transport, compared with the satellite dust extinction profiles. After updating dust refractive indices and correcting for a bias in partitioning size-segregated emissions, both SW cooling and longwave (LW) warming of dust simulated by E3SMv1 are increased and agree better with other recent studies. The estimated net dust DRE of –0.42 Wm –2 represents a stronger cooling effect than the observationally based estimate –0.2 Wm –2 (–0.48 to +0.2), due to a smaller LW warming. Constrained by a global mean DAOD, model sensitivity studies of increasing horizontal and vertical resolution show strong influences on the simulated global dust burden and lifetime primarily through the change of dust dry deposition rate; there are also remarkable differences in simulated spatial distributions of DAOD, DRE, and deposition fluxes. Thus, constraining the global DAOD is insufficient for accurate representation of dust climate effects, especially in transitioning to higher- or variable-resolution ESMs. Better observational constraints of dust vertical profiles, dry deposition, size, and LW properties are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Stable Simulation of the Community Atmosphere Model Using Machine‐Learning Physical Parameterization Trained With Experience Replay

In recent years, machine learning (ML) models have been used to improve physical parameterizations of general circulation models (GCMs). A significant challenge of integrating ML models into GCMs is the online instability when they are coupled for long‐term simulation. We present a new strategy that demonstrates robust online stability when the physical parameterization package of an atmospheric GCM is replaced by a deep ML model. The method uses experience replay with a multistep training scheme of the ML model in which the model's own output at the previous time step is used in the training. Predicted physics tendencies in the replay buffer with the most recent errors in the training iterations are reused, making the ML model learn from its own errors. The training method reduces the gap between the offline and online environments of the ML model. The method is used to train the ML model as the physical parameterization of the Community Atmosphere Model (CAM5) with training data from the Multi‐scale Modeling Framework high resolution simulations. Three 6‐year online simulations of the CAM5 are carried out by using the ML physics package. The simulated spatial distributions of precipitation, surface temperature and zonally averaged atmospheric fields demonstrate overall better accuracy than that of the standard CAM5 and benchmark model even without the use of additional physical constraints or tuning. This work is the first to demonstrate a solution to address the online instability problem in climate modeling with ML physics by using experience replay.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Laboratory measurements of the physics of auroral electron acceleration by Alfvén waves

While the aurora has attracted attention for millennia, important questions remain unanswered. Foremost is how auroral electrons are accelerated before colliding with the ionosphere and producing auroral light. Powerful Alfvén waves are often found traveling Earthward above auroras with sufficient energy to generate auroras, but there has been no direct measurement of the processes by which Alfvén waves transfer their energy to auroral electrons. Here, we show laboratory measurements of the resonant transfer of energy from Alfvén waves to electrons under conditions relevant to the auroral zone. Experiments are performed by launching Alfvén waves and simultaneously recording the electron velocity distribution. Numerical simulations and analytical theory support that the measured energy transfer process produces accelerated electrons capable of reaching auroral energies. The experiments, theory, and simulations demonstrate a clear causal relationship between Alfvén waves and accelerated electrons that directly cause auroras.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Upper troposphere dust belt formation processes vary seasonally and spatially in the Northern Hemisphere

Dust aerosols impact global energy balance substantially by acting as efficient ice nuclei to alter cold cloud properties. However, the estimate of dust indirect effect remains uncertain due to simulating dust distributions poorly and lacking reliable dust observations, especially in the upper-troposphere. Here, we characterize and understand upper-troposphere dust sources and transport with an improved dust dataset derived from A-train satellite lidar and radar measurements and an air parcel trajectory model. The distinct upper-troposphere dust belt over the northern hemisphere has seasonally varying base and top heights of 3.65 ± 2.84 and 8.35 ± 1.50 km above mean sea level and its column loading is strongest during spring (March-April-May). The out-of-phase annual cycles of mid-level dust concentration and westerly wind over source regions control the seasonal upper-tropospheric dust loading variations. African deserts contribute the most (46.3%) to the upper-troposphere dust belt in spring and the synoptic trough is the leading (49%) dust lifting mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Deep‐learning‐based power distribution network switch action identification leveraging dynamic features of distributed energy resources

This study proposes a data‐driven approach for identifying switch actions in power distribution networks. Simulated micro‐phasor measurement unit data is utilised to train a convolutional neural network (CNN) model. The trained CNN model can identify multi‐phase multi‐switch actions. Instead of working as a blackbox, the proposed approach extracts the features from the hidden layers of the trained CNN for engineering interpretation and error check. In addition, a random‐forest‐based feature ranking algorithm is proposed to identify the most important features. The proposed approach is validated on the IEEE 123‐node feeder modelled in GridLAB‐D. The CNN model is built and trained using TensorFlow. The proposed approach achieves 96.57% identification accuracy.

42 ENGINEERING↗
DOE OSTI2014

Modeling delay in genetic networks: From delay birth-death processes to delay stochastic differential equations

Delay is an important and ubiquitous aspect of many biochemical processes. For example, delay plays a central role in the dynamics of genetic regulatory networks as it stems from the sequential assembly of first mRNA and then protein. Genetic regulatory networks are therefore frequently modeled as stochastic birth-death processes with delay. Here, we examine the relationship between delay birth-death processes and their appropriate approximating delay chemical Langevin equations. We prove a quantitative bound on the error between the pathwise realizations of these two processes. Our results hold for both fixed delay and distributed delay. Simulations demonstrate that the delay chemical Langevin approximation is accurate even at moderate system sizes. It captures dynamical features such as the oscillatory behavior in negative feedback circuits, cross-correlations between nodes in a network, and spatial and temporal information in two commonly studied motifs of metastability in biochemical systems. Overall, these results provide a foundation for using delay stochastic differential equations to approximate the dynamics of birth-death processes with delay.

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

Resolution of the vibrational energy distribution function using a direct simulation Monte Carlo-master equation approach

The direct simulation Monte Carlo (DSMC) method is the primary numerical technique for analysis of rarefied gas flows. While recent progress in computational chemistry is beginning to provide vibrationally resolved transition and reaction cross sections that can be employed in DSMC calculations, the particle nature of the standard DSMC method makes it difficult to use this information in a statistically significant way. The current study introduces a new technique that makes it possible to resolve all of the vibrational energy levels by using a master equation approach along with temperature-dependent transition rates. The new method is compared to the standard DSMC technique for several heat bath and shock wave conditions and demonstrates the ability to resolve the full vibrational manifold at the expected overall rates of relaxation. The ability of the new master equation approach to the DSMC method for resolving, in particular, the high-energy states addresses a well-known, longstanding deficiency of the standard DSMC method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Transducer selection and application in magnetoacoustic tomography with magnetic induction

As an acoustic receiver, transducer plays a vital role in signal acquisition and image reconstruction for magnetoacoustic tomography with magnetic induction (MAT-MI). In order to optimize signal acquisition, the expressions of acoustic pressure detection and waveform collection are theoretically studied based on the radiation theory of acoustic dipole and the reception pattern of transducer. Pressure distributions are simulated for a cylindrical phantom model using a planar piston transducer with different radii and bandwidths. The proposed theory is also verified by the experimental measurements of acoustic waveform detection for an aluminum foil cylinder. It is proved that acoustic pressure with sharp and clear boundary peaks can be detected by the large-radius transducer with wide bandwidth, reflecting the differential of the induced Lorentz force accurately, which is helpful for precise conductivity reconstruction. To detect acoustic pressure with acceptable pressure amplitude, peak pressure ratio, amplitude ratio, and improved signal to noise ratio, the scanning radius of 5–10 times the radius of the object should be selected to improve the accuracy of image reconstruction. This study provides a theoretical and experimental basis for transducer selection and application in MAT-MI to obtain reconstructed images with improved resolution and definition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Formation of snowflake domains during fast cooling of lithium tantalate crystals

Formation of the original dendrite snowflake-shape domains during fast cooling after heating above phase transition temperature by pulse laser irradiation was revealed in congruent lithium tantalate crystals. The effect was attributed to polarization reversal under the action of spatially nonuniform pyroelectric field. Two stages of the domain shape evolution at the surface were separated: (1) growth of circular domains by sideways motion of the domain walls and (2) backswitching leading to formation of the snowflake domains. The simulated spatial distribution of the pyroelectric field in regular two-dimensional structure was used for an explanation of the obtained results. The backswitching process in the surface layer has been attributed to change of the sign of the pyroelectric field at the domain wall. The snowflake domain shape is caused by the formation of isolated nanodomain fingers and hampering of their merging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

High-speed nanoscale characterization of dewetting via dynamic transmission electron microscopy

The dewetting of thin films can occur in either the solid or the liquid state for which different mass transport mechanisms are expected to control morphological changes. Traditionally, dewetting dynamics have been examined on time scales between several seconds to hours, and length scales ranging between nanometers and millimeters. The determination of mass transport mechanisms on the nanoscale, however, requires nanoscale spatial resolution and much shorter time scales. This study reports the high-speed observation of dewetting phenomena for kinetically constrained Ni thin films on crystalline SrTiO{sub 3} substrates. Movie-mode Dynamic Transmission Electron Microscopy (DTEM) was used for high-speed image acquisition during thin film dewetting at different temperatures. DTEM imaging confirmed that the initial stages of film agglomeration include edge retraction, hole formation, and growth. Finite element modeling was used to simulate temperature distributions within the DTEM samples after laser irradiation with different energies. For pulsed laser irradiation at 18 μJ, experimentally observed hole growth suggests that Marangoni flow dominates hole formation in the liquid nickel film. After irradiation with 13.8 μJ, however, the observations suggest that dewetting was initiated by nucleation of voids followed by hole growth through solid-state surface diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A co-axial velocity map imaging spectrometer for electrons

Here, we present the design of an electron velocity map imaging (VMI) spectrometer where the ionizing laser source propagates along the symmetry axis of the spectrometer. The co-axial geometry is useful in a variety of experiments, because it provides a unique 2-dimensional projection of the 3-dimensional electron momentum distribution. Initial simulations show that this co-axial VMI can work with both high energy (more than 100 eV) and low energy (tens of eV) electrons. We demonstrate the performance of this co-axial VMI spectrometer at the Linac Coherent Light Source X-ray Free Electron Laser facility.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Reconnection from a turbulence perspective

The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power law behavior, with a slope near −5/3 at wave numbers smaller than the inverse ion inertial length and a slope steeper than −8/3 for larger wave numbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is the fact that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70° to about 48° (nearly isotropic) by the end of the simulation. The distribution of energy over length scales is further analyzed by dividing the domain into spatial subregions and employing structure functions.

Adhikari, S. (ORCID:0000000321607066)↗
DOE OSTI2023

Modeling of the spray-induced wall stress acting on the ignition assistance device

This research introduces a novel wall-stress model called the Spray-Induced Wall Stress (SIWS) model, which considers the effects of spray-wall impingement and the resulting formation of wall stress within the Lagrangian spray modeling framework. The primary objective of this paper is to provide a mathematical description of the fundamental physics underlying the model. Subsequently, the proposed model is validated using existing experimental data. Here, the remainder of the study focuses on the practical application of the model to an ignition assistance device. Specifically, this device is installed in a compression ignition engine and designed to enhance ignition in aviation-fueled high-altitude aircraft propulsion systems. The research sheds light on the mechanical impulse caused by the high-speed impact of the spray jet, leading to the accumulation of mechanical stress on the rigid body of the ignition assistance device. Previous studies on fluid-structure interaction have only considered the interaction between the gas phase and the solid wall. However, the SIWS model incorporates the additional impact of the impinging liquid spray jet. Consequently, the simulated stress distribution on the ignition assistance device can be estimated by considering both the gas-phase-induced term and the spray-induced term simultaneously.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

THE CLOSE BINARY FRACTION OF DWARF M STARS

We describe a search for close spectroscopic dwarf M star binaries using data from the Sloan Digital Sky Survey to address the question of the rate of occurrence of multiplicity in M dwarfs. We use a template-fitting technique to measure radial velocities from 145,888 individual spectra obtained for a magnitude-limited sample of 39,543 M dwarfs. Typically, the three or four spectra observed for each star are separated in time by less than four hours, but for {approx}17% of the stars, the individual observations span more than two days. In these cases we are sensitive to large-amplitude radial velocity variations on timescales comparable to the separation between the observations. We use a control sample of objects having observations taken within a four-hour period to make an empirical estimate of the underlying radial velocity error distribution and simulate our detection efficiency for a wide range of binary star systems. We find the frequency of binaries among the dwarf M stars with a < 0.4 AU to be 3%-4%. Comparison with other samples of binary stars demonstrates that the close binary fraction, like the total binary fraction, is an increasing function of primary mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ON THE EFFECTS OF THE EVOLUTION OF MICROBIAL MATS AND LAND PLANTS ON THE EARTH AS A PLANET. PHOTOMETRIC AND SPECTROSCOPIC LIGHT CURVES OF PALEO-EARTHS

Understanding the spectral and photometric variability of the Earth and the rest of the solar system planets has become of utmost importance for the future characterization of rocky exoplanets. As this is not only interesting at present times but also along the planetary evolution, we studied the effect that the evolution of microbial mats and plants over land has had on the way our planet looks from afar. As life evolved, continental surfaces changed gradually and non-uniformly from deserts through microbial mats to land plants, modifying the reflective properties of the ground and most likely the distribution of moisture and cloudiness. Here, we used a radiative transfer model of the Earth, together with geological paleo-records of the continental distribution and a reconstructed cloud distribution, to simulate the visible and near-IR radiation reflected by our planet as a function of Earth's rotation. We found that the evolution from deserts to microbial mats and to land plants produces detectable changes in the globally averaged Earth's reflectance. The variability of each surface type is located in different bands and can induce reflectance changes of up to 40% in period of hours. We conclude that by using photometric observations of an Earth-like planet at different photometric bands it would be possible to discriminate between different surface types. While recent literature proposes the red-edge feature of vegetation near 0.7 {mu}m as a signature for land plants, observations in near-IR bands can be equally or even better suited for this purpose.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Jet emission in young radio sources: A Fermi large area telescope gamma-ray view

We investigate the contribution of the beamed jet component to the high-energy emission in young and compact extragalactic radio sources, focusing for the first time on the γ-ray band. We derive predictions on the γ-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high-energy emission is produced via Compton scattering by the relativistic electrons in a spherical region at the considered scales (≲10 kpc). Simulations show a wide range of γ-ray luminosities, with intensities up to ∼10{sup 46}-10{sup 48} erg s{sup –1} depending on the assumed jet parameters. We find a highly linear relation between the simulated X-ray and γ-ray luminosities that can be used to select candidates for γ-ray detection. We compare the simulated luminosity distributions in the radio, X-ray, and γ-ray regimes with observations for the largest sample of X-ray-detected young radio quasars. Our analysis of ∼4-yr Fermi Large Area Telescope (LAT) data does not yield any statistically significant detections. However, the majority of the model-predicted γ-ray fluxes for the sample are near or below the current Fermi-LAT flux threshold and compatible with the derived upper limits. Our study gives constraints on the minimum jet power (L {sub jet,} {sub kin}/L {sub disk} > 0.01) of a potential jet contribution to the X-ray emission in the most compact sources (≲ 1 kpc) and on the particle-to-magnetic field energy density ratio that are in broad agreement with equipartition assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Gravitational wave source counts at high redshift and in models with extra dimensions

Gravitational wave (GW) source counts have been recently shown to be able to test how gravitational radiation propagates with the distance from the source. Here, we extend this formalism to cosmological scales, i.e. the high redshift regime, and we discuss the complications of applying this methodology to high redshift sources. We also allow for models with compactified extra dimensions like in the Kaluza-Klein model. Furthermore, we also consider the case of intermediate redshifts, i.e. 0 < z ∼< 1, where we show it is possible to find an analytical approximation for the source counts dN / d ( S / N ). This can be done in terms of cosmological parameters, such as the matter density Ω {sub m} {sub ,0} of the cosmological constant model or the cosmographic parameters for a general dark energy model. Our analysis is as general as possible, but it depends on two important factors: a source model for the black hole binary mergers and the GW source to galaxy bias. This methodology also allows us to obtain the higher order corrections of the source counts in terms of the signal-to-noise S / N . We then forecast the sensitivity of future observations in constraining GW physics but also the underlying cosmology by simulating sources distributed over a finite range of signal-to-noise with a number of sources ranging from 10 to 500 sources as expected from future detectors. We find that with 500 events it will be possible to provide constraints on the matter density parameter at present Ω {sub m} {sub ,0} on the order of a few percent and with the precision growing fast with the number of events. In the case of extra dimensions we find that depending on the degeneracies of the model, with 500 events it may be possible to provide stringent limits on the existence of the extra dimensions if the aforementioned degeneracies can be broken.

79 ASTRONOMY AND ASTROPHYSICS↗