Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATOR”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Partially coherent wavefront propagation simulations: Mirror and monochromator crystal quality assessment

We have applied fully-and partially-coherent synchrotron radiation wavefront propagation simulations, implemented in the “Synchrotron Radiation Workshop” (SRW) computer code, to analyse the effects of imperfect mirrors and monochromator at the Coherent Hard X-ray beamline. This beamline is designed for X-ray Photon Correlation Spectroscopy, a technique that heavily relies on the partial coherence of the X-ray beam and benefits from a careful preservation of the X-ray wavefront. We present simulations and a comparison with the measured beam profile at the sample position, which show the impact of imperfect optics on the wavefront.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Phonocatalysis. An ab initio simulation experiment

Using simulations, we postulate and show that heterocatalysis on large-bandgap semiconductors can be controlled by substrate phonons, i.e., phonocatalysis. With ab initio calculations, including molecular dynamic simulations, the chemisorbed dissociation of XeF{sub 6} on h-BN surface leads to formation of XeF{sub 4} and two surface F/h-BN bonds. The reaction pathway and energies are evaluated, and the sorption and reaction emitted/absorbed phonons are identified through spectral analysis of the surface atomic motion. Due to large bandgap, the atomic vibration (phonon) energy transfer channels dominate and among them is the match between the F/h-BN covalent bond stretching and the optical phonons. We show that the chemisorbed dissociation (the pathway activation ascent) requires absorption of large-energy optical phonons. Then using progressively heavier isotopes of B and N atoms, we show that limiting these high-energy optical phonons inhibits the chemisorbed dissociation, i.e., controllable phonocatalysis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Density-functional theory computer simulations of CZTS{sub 0.25}Se{sub 0.75} alloy phase diagrams

Density-functional theory simulations of CZTS, CZTSe, and CZTS{sub 0.25}Se{sub 0.75} photovoltaic compounds have been performed to investigate the stability of the CZTS{sub 0.25}Se{sub 0.75} alloy vs. decomposition into CZTS, CZTSe, and other secondary compounds. The Gibbs energy for vibrational contributions was estimated by calculating phonon spectra and thermodynamic properties at finite temperatures. It was demonstrated that the CZTS{sub 0.25}Se{sub 0.75} alloy is stabilized not by enthalpy of formation but primarily by the mixing contributions to the Gibbs energy. The Gibbs energy gains/losses for several decomposition reactions were calculated as a function of temperature with/without intermixing and vibration contributions to the Gibbs energy. A set of phase diagrams was built in the multidimensional space of chemical potentials at 300 K and 900 K temperatures to demonstrate alloy stability and boundary compounds at various chemical conditions. It demonstrated for CZTS{sub 0.25}Se{sub 0.75} that the chemical potentials for stability differ between typical processing temperature (∼900 K) and operating temperature (300 K). This implies that as cooling progresses, the flux/concentration of S should be increased in MBE growth to maintain the CZTS{sub 0.25}Se{sub 0.75} in a thermodynamically stable state to minimize phase decomposition.

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

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Case Study on Pathogen Transport, Deposition, Evaporation and Transmission: Linking High-Fidelity Computational Fluid Dynamics Simulations to Probability of Infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose–response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Formation of compact clusters from high resolution hybrid cosmological simulations

The early universe hosted a large population of small dark matter 'minihalos' that were too small to cool and form stars on their own. These existed as static objects around larger galaxies until acted upon by some outside influence. Outflows, which have been observed around a variety of galaxies, can provide this influence in such a way as to collapse, rather than disperse, the minihalo gas. Gray and Scannapieco performed an investigation in which idealized spherically symmetric minihalos were struck by enriched outflows. Here we perform high-resolution cosmological simulations that form realistic minihalos, which we then extract to perform a large suite of simulations of outflow-minihalo interactions including non-equilibrium chemical reactions. In all models, the shocked minihalo forms molecules through non-equilibrium reaction, and then cools to form dense, chemically homogenous clumps of star-forming gas. The formation of these high-redshift clusters may be observable with the next generation of telescopes and the largest of them should survive to the present-day, having properties similar to halo globular clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Numerical simulations of optically thick accretion onto a black hole. II. Rotating flow

In this paper, we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general M {sub 1} closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work, we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

WHISTLER TURBULENCE FORWARD CASCADE VERSUS INVERSE CASCADE: THREE-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS

We present the results of the first fully three-dimensional particle-in-cell simulations of decaying whistler turbulence in a magnetized, homogeneous, collisionless plasma in which both forward cascades to shorter wavelengths, and inverse cascades to longer wavelengths are allowed to proceed. For the electron beta β {sub e} = 0.10 initial value considered here, the early-time rate of inverse cascade is very much smaller than the rate of forward cascade, so that at late times the fluctuation energy in the regime of the inverse cascade is much weaker than that in the forward cascade regime. Similarly, the wavevector anisotropy in the inverse cascade regime is much weaker than that in the forward cascade regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SciDAC Center for Simulation of Wave-Plasma Interactions - Iterated Finite-Orbit Monte Carlo Simulations with Full-Wave Fields for Modeling Tokamak ICRF Wave Heating Experiments - Final Report

This final report describes the work performed under U.S. Department of Energy Cooperative Agreement DE-FC02-08ER54954 for the period April 1, 2011 through March 31, 2013. The goal of this project was to perform iterated finite-orbit Monte Carlo simulations with full-wall fields for modeling tokamak ICRF wave heating experiments. In year 1, the finite-orbit Monte-Carlo code ORBIT-RF and its iteration algorithms with the full-wave code AORSA were improved to enable systematical study of the factors responsible for the discrepancy in the simulated and the measured fast-ion FIDA signals in the DIII-D and NSTX ICRF fast-wave (FW) experiments. In year 2, ORBIT-RF was coupled to the TORIC full-wave code for a comparative study of ORBIT-RF/TORIC and ORBIT-RF/AORSA results in FW experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Star cluster formation in cosmological simulations. I. Properties of young clusters

We present a new implementation of star formation in cosmological simulations by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha \approx 1.8\mbox{–}2$, while the cutoff at high mass scales with the star formation rate (SFR). A related trend is a positive correlation between the surface density of the SFR and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major-merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. As a result, comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Emulating Simulations of Cosmic Dawn for 21 cm Power Spectrum Constraints on Cosmology, Reionization, and X-Ray Heating

Current and upcoming radio interferometric experiments are aiming to make a statistical characterization of the high-redshift 21 cm fluctuation signal spanning the hydrogen reionization and X-ray heating epochs of the universe. However, connecting 21 cm statistics to the underlying physical parameters is complicated by the theoretical challenge of modeling the relevant physics at computational speeds quick enough to enable exploration of the high-dimensional and weakly constrained parameter space. In this work, we use machine learning algorithms to build a fast emulator that can accurately mimic an expensive simulation of the 21 cm signal across a wide parameter space. We embed our emulator within a Markov Chain Monte Carlo framework in order to perform Bayesian parameter constraints over a large number of model parameters, including those that govern the Epoch of Reionization, the Epoch of X-ray Heating, and cosmology. As a worked example, we use our emulator to present an updated parameter constraint forecast for the Hydrogen Epoch of Reionization Array experiment, showing that its characterization of a fiducial 21 cm power spectrum will considerably narrow the allowed parameter space of reionization and heating parameters, and could help strengthen Planck's constraints on σ{sub 8}. We provide both our generalized emulator code and its implementation specifically for 21 cm parameter constraints as publicly available software.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A review of the remote sensing of lower tropospheric thermodynamic profiles and its indispensable role for the understanding and the simulation of water and energy cycles: REMOTE SENSING OF THERMODYNAMIC PROFILES

A review of remote sensing technology for lower tropospheric thermodynamic (TD) profiling is presented with focus on high accuracy and high temporal-vertical resolution. The contributions of these instruments to the understanding of the Earth system are assessed with respect to radiative transfer, land surface-atmosphere feedback, convection initiation, and data assimilation. We demonstrate that for progress in weather and climate research, TD profilers are essential. These observational systems must resolve gradients of humidity and temperature in the stable or unstable atmospheric surface layer close to the ground, in the mixed layer, in the interfacial layer—usually characterized by an inversion—and the lower troposphere. A thorough analysis of the current observing systems is performed revealing significant gaps that must be addressed to fulfill existing needs. We analyze whether current and future passive and active remote sensing systems can close these gaps. A methodological analysis and demonstration of measurement capabilities with respect to bias and precision is executed both for passive and active remote sensing including passive infrared and microwave spectroscopy, the global navigation satellite system, as well as water vapor and temperature Raman lidar and water vapor differential absorption lidar. Whereas passive remote sensing systems are already mature with respect to operational applications, active remote sensing systems require further engineering to become operational in networks. However, active remote sensing systems provide a smaller bias as well as higher temporal and vertical resolutions. For a suitable mesoscale network design, TD profiler system developments should be intensified and dedicated observing system simulation experiments should be performed.

58 GEOSCIENCES↗
DOE OSTI2015

How does increasing horizontal resolution in a global climate model improve the simulation of aerosol-cloud interactions?

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25 degrees, with the meteorology nudged towards the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the northern hemisphere mid-latitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m⁻² (-30%) compared to the 2° model, leading to a 0.26 W m⁻² reduction (-15%) in the global annual average AIF.

Ma, Po -Lun↗
DOE OSTI2015

Modifications to WRF's dynamical core to improve the treatment of moisture for large-eddy simulations: WRF DY-CORE MOISTURE TREATMENT

Yamaguchi and Feingold (2012) note that the cloud fields in their Weather Research and Forecasting (WRF) large-eddy simulations (LESs) of marine stratocumulus exhibit a strong sensitivity to time stepping choices. In this study, we reproduce and analyze this sensitivity issue using two stratocumulus cases, one marine and one continental. Results show that (1) the sensitivity is associated with spurious motions near the moisture jump between the boundary layer and the free atmosphere, and (2) these spurious motions appear to arise from neglecting small variations in water vapor mixing ratio (qv) in the pressure gradient calculation in the acoustic sub­stepping portion of the integration procedure. We show that this issue is remedied in the WRF dynamical core by replacing the prognostic equation for the potential temperature θ with one for the moist potential temperature θm=θ(1+1.61qv), which allows consistent treatment of moisture in the calculation of pressure during the acoustic sub­steps. With this modification, the spurious motions and the sensitivity to the time stepping settings (i.e., the dynamic time step length and number of acoustic sub­steps) are eliminated in both of the example stratocumulus cases. This modification improves the applicability of WRF for LES applications, and possibly other models using similar dynamical core formulations, and also permits the use of longer time steps than in the original code.

Xiao, Heng↗
DOE OSTI2017

Global three-dimensional simulation of Earth's dayside reconnection using a two-way coupled magnetohydrodynamics with embedded particle-in-cell model: initial results: 3D MHD-EPIC simulation of magnetosphere

Here, we perform a three-dimensional (3D) global simulation of Earth's magnetosphere with kinetic reconnection physics to study the flux transfer events (FTEs) and dayside magnetic reconnection with the recently developed magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC). During the one-hour long simulation, the FTEs are generated quasi-periodically near the subsolar point and move toward the poles. We also find the magnetic field signature of FTEs at their early formation stage is similar to a ‘crater FTE’, which is characterized by a magnetic field strength dip at the FTE center. After the FTE core field grows to a significant value, it becomes an FTE with typical flux rope structure. When an FTE moves across the cusp, reconnection between the FTE field lines and the cusp field lines can dissipate the FTE. The kinetic features are also captured by our model. A crescent electron phase space distribution is found near the reconnection site. A similar distribution is found for ions at the location where the Larmor electric field appears. The lower hybrid drift instability (LHDI) along the current sheet direction also arises at the interface of magnetosheath and magnetosphere plasma. Finally, the LHDI electric field is about 8 mV/m and its dominant wavelength relative to the electron gyroradius agrees reasonably with MMS observations.

58 GEOSCIENCES↗