Search NASA⌕ Search

SEARCH · Search NASA

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 757 records · Page 42

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↗

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↗

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↗

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↗

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↗

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↗

A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state: A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state

The development of numerical approaches to perform direct numerical simulations of compressible multiphase flows has been an active field of research for several years. Proper treatment of fluid interfaces is crucial as important physics occur in this infinitesimally small region. Furthermore, the compressibility of the fluid requires proper treatment of discontinuities. Artificial diffusivity is among a number of methods widely used for compressible flows. This study develops a general form of consistent artificial diffusion fluxes and extends the localized artificial diffusivity method for high-order central schemes to solve multiphase flows with an interface-capturing method. These fluxes ensure an oscillation-free interface for pressure, velocity, and temperature without employing a sharpening technique. Moreover, the high-order representation of all scales in the flow helps capture the wide range of instabilities inherent in these flows. The goal is to develop an approach capable of performing high-fidelity simulations supported by physics-driven validation. lastly, this is achieved by solving the five-equation model with the stiffened-gas equation of state using the proposed method for multicomponent and multiphase flows on a variety of 1D and 2D problems.

42 ENGINEERING↗

FAST.Farm Development and Validation of Structural Load Prediction Against Large Eddy Simulations

FAST.Farm is a mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. FAST.Farm is based on several principles of the dynamic wake meandering (DWM) model, but also addresses limitations of previous DWM implementations. Previous FAST.Farm studies have shown the similarities and differences between FAST.Farm and large eddy simulations for rigid turbine cases. The objective of this work is to quantify the ability of FAST.Farm to accurately predict turbine structural response in a small wind farm. This is done by comparing FAST.Farm structural response to SOWFA-OpenFAST results for three laterally-aligned turbines. The purpose of this study is to characterize the similarities and differences between FAST.Farm and a higher fidelity model for predicting turbine structural response in a wind farm for differing atmospheric inflows. Strong statistical agreement was found between FAST.Farm and SOWFA-OpenFAST structural response for the non-waked upstream turbine, and good agreement was found for the downstream turbines for most structural quantities. Higher differences were seen for downstream turbines with low ambient turbulence intensity or yawed turbines, suggesting areas for FAST.Farm wake dynamics modeling improvements. For all cases and turbines, small statistical differences were seen between blade deflections and bending moments, with larger differences for tower-top and tower-base bending moments. Overall, the results establish confidence for applying FAST.Farm to wind farm power and loads analyses and identify areas where further model validations and model improvements should be targeted.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Development of a Low-Level 39Ar Calibration Standard – Analysis by Absolute Gas Counting Measurements Augmented with Simulation

This paper describes the generation of 39Ar, via reactor irritation of potassium carbonate, followed by quantitative analysis (length-compensated proportional counting) to yield calibration standards that are approximately 50 and 3 times atmospheric background levels. Measurements were performed in Pacific Northwest National Laboratory’s shallow underground counting laboratory studying the effect of gas density on beta-transport; these results are compared with simulation. The total expanded uncertainty of the specific activity for the ~50x 39Ar in P10 standard is 3.6% (k=2).

Absolute radiation measurement↗

Molecular dynamics simulation of sodium aluminosilicate glass structures and glass surface-water reactions using the reactive force field (ReaxFF)

Reactive potentials are increasingly used to study the properties of glasses and glass water reactions in a reactive molecular dynamics (MD) framework. In this study, we have simulated a ternary sodium aluminosilicate glass and investigated the initial stages of the glass surface-water reactions at 300 K using reactive force field (ReaxFF). On comparison of the simulated glass structures generated using ReaxFF and classical Buckingham potentials, our results show that the atomic density profiles calculated for the surface glass structures indicate a bond-angle distribution dependency. The atomic density profiles also show higher concentrations of non-bridging oxygens (NBOs) and sodium ions at the glass surface. Additionally, we present our results of formation of silanol species and the diffusion of water molecules at the glass surface using ReaxFF.

Chemistry↗

Evaluation of a multi-scale WRF-CAM5 simulation during the 2010 East Asian Summer Monsoon

The Weather Research and Forecasting model with Chemistry (WRF-Chem) with the physics package of the Community Atmosphere Model Version 5 (CAM5) has been applied at multiple scales over Eastern China (EC) and the Yangtze River Delta (YRD) to evaluate how increased horizontal resolution with physics designed for a coarser resolution climate model impacts aerosols and clouds, and the resulting precipitation characteristics and performance during the 2010 East Asian Summer Monsoon (EASM). Despite large underpredictions in surface aerosol concentrations and aerosol optical depth, there is good spatial agreement with surface observations of chemical predictions, and increasing spatial resolution tends to improve performance. Model bias and normalized root mean square values for precipitation predictions are relatively small, but there are significant differences when comparing modeled and observed probability density functions for precipitation in EC and YRD. Increasing model horizontal resolution tends to reduce model bias and error for precipitation predictions. The surface and column aerosol loading is maximized between about 32N and 42N in early to mid-May during the 2010 EASM, and then shifts north while decreasing in magnitude during July and August. Changing model resolution moderately changes the spatiotemporal relationships between aerosols, cloud properties, and precipitation during the EASM, thus demonstrating the importance of model grid resolution in simulating EASM circulation and rainfall patterns over EC and the YRD. Results from this work demonstrate the capability and limitations in the aerosol, cloud, and precipitation representation of WRF-CAM5 for regional-scale applications down to relatively fine horizontal resolutions. Further WRF-CAM5 model development and application in this area is needed.

Campbell, Patrick↗