Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble 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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High yield polar direct drive fusion neutron sources at the National Ignition Facility

Polar direct drive neutron source experiments were performed at the National Ignition Facility showing substantial improvement in total neutron yield and efficiency of conversion of laser energy to fusion output. Plastic capsules 3–4 mm in diameter were filled with 1.5 mg/cc of deuterium–tritium (DT) fuel and imploded with laser beam pointing and defocus designed to compensate for polar asymmetry introduced by the facility beam entrance angles. Radiation-hydrodynamics simulations were employed to optimize the multi-dimensional laser and target parameter space, within facility and target fabrication constraints. Ensembles of 1D simulations tuned to match the outputs of early shots in the series were used to design subsequent shots in the series. This allowed the later shots to be designed based on empirically motivated sensitivities to laser and target input parameters, while eliminating the need to explicitly model phenomena such as hydrodynamic instabilities and nonlinear laser–plasma interactions. Additionally, one experiment with a 3.0 mm diameter CH capsule produced 13.6 kJ (4.81 × 10 15 DT neutrons) from a laser input below the NIF optics damage threshold at 585 kJ, 328 TW. Two experiments with 4.0 mm capsules produced 31.3 and 33.6 kJ of fusion output (1.11 × 10 16 and 1.19 × 10 16 DT neutrons) with 1.10 MJ, 390 TW and 1.26 MJ, 425 TW of laser input, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2022

Are Soybean Models Ready for Climate Change Food Impact Assessments?

An accurate estimation of crop yield under climate change scenarios is essential to quantify our ability to feed a growing population and develop agronomic adaptations to meet future food demand. A coordinated evaluation of yield simulations from process-based eco-physiological models for climate change impact assessment is still missing for soybean, the most widely grown grain legume and the main source of protein in our food chain. In this first soybean multi-model study, we used ten prominent models capable of simulating soybean yield under varying temperature and atmospheric CO 2 concentration [CO 2 ] to quantify the uncertainty in soybean yield simulations in response to these factors. Models were first parametrized with high quality measured data from five contrasting environments. We found considerable variability among models in simulated yield responses to increasing temperature and [CO 2 ]. For example, under a + 3 °C temperature rise in our coolest location in Argentina, some models simulated that yield would reduce as much as 24%, while others simulated yield increases up to 29%. In our warmest location in Brazil, the models simulated a yield reduction ranging from a 38% decrease under + 3 °C temperature rise to no effect on yield. Similarly, when increasing [CO 2 ] from 360 to 540 ppm, the models simulated a yield increase that ranged from 6% to 31%. Model calibration did not reduce variability across models but had an unexpected effect on modifying yield responses to temperature for some of the models. The high uncertainty in model responses indicates the limited applicability of individual models for climate change food projections. However, the ensemble mean of simulations across models was an effective tool to reduce the high uncertainty in soybean yield simulations associated with individual models and their parametrization. Ensemble, ensemble mean yield responses to temperature and [CO 2 ] were similar to those reported from the literature. Our study is the first demonstration of the benefits achieved from using an ensemble of grain legume models for climate change food projections, and highlights that further soybean model development with experiments under elevated [CO 2 ] and temperature is needed to reduce the uncertainty from the individual models.

Agricultural Model Inter-comparison and Improvemen↗
DOE OSTI2024

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗
DOE OSTI2021

Characterizing the acceleration time of laser-driven ion acceleration with data-informed neural networks

Peak ion energy is an important figure-of-merit in short-pulse, laser-driven ion acceleration and is dependent on an associated acceleration time. Standard metrics for these quantities depend on analytical results such as the self-similar fluid model or empirical models based on relatively small experimental and simulation datasets. In this work we attempt to use a data-informed neural network (NN) as a surrogate model for a large ensemble of PIC simulations to investigate an effective acceleration time. We explore the application of a stacked convolutional and recurrent NN architecture for improved regression by incorporating the time dependencies of the data into the training process. Of particular note is how pretraining a network on lower fidelity data, e.g. 1D analytical results, greatly improves the network's ability to learn more complex, higher fidelity data. Finally, the dependency of the acceleration time on various laser and plasma parameters is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

The Effects of Lightning NO(x) Production during the July 21 EULINOX Storm studied with a 3-D Cloud-scale Chemical Transport Model

The July 21,1998 thunderstonn observed during the European Lightning Nitrogen Oxides Project (EULINOX) project was simulated using the three-dimensional Goddard Cumulus Ensemble (GCE) model. The simulation successfully reproduced a number of observed storm features including the splitting of the original cell into a southern cell which developed supercell characteristics, and a northern cell which became multicellular. Output from the GCE simulation was used to drive an offline cloud-scale chemical transport model which calculates tracer transport and includes a parameterization of lightning NO(x) production which uses observed flash rates as input. Estimates of lightning NO(x) production were deduced by assuming various values of production per intracloud and production per cloud-to-ground flash and comparing the results with in-cloud aircraft observations. The assumption that both types of flashes produce 360 moles of NO per flash on average compared most favorably with column mass and probability distribution functions calculated from observations. This assumed production per flash corresponds to a global annual lightning NOx source of 7 Tg N per yr. Chemical reactions were included in the model to evaluate the impact of lightning NO(x), on ozone. During the storm, the inclusion of lightning NOx in the model results in a small loss of ozone (on average less than 4 ppbv) at all model levels. Simulations of the chemical environment in the 24 hours following the storm show on average a small increase in the net production of ozone at most levels resulting from lightning NO(x), maximizing at approximately 5 ppbv per day at 5.5 km. Between 8 and 10.5 km, lightning NO(x) causes decreased net ozone production.

Ott, Lesley E.↗
NASA NTRS1994

Measurement uncertainty and feasibility study of a flush airdata system for a hypersonic flight experiment

Presented is a feasibility and error analysis for a hypersonic flush airdata system on a hypersonic flight experiment (HYFLITE). HYFLITE heating loads make intrusive airdata measurement impractical. Although this analysis is specifically for the HYFLITE vehicle and trajectory, the problems analyzed are generally applicable to hypersonic vehicles. A layout of the flush-port matrix is shown. Surface pressures are related airdata parameters using a simple aerodynamic model. The model is linearized using small perturbations and inverted using nonlinear least-squares. Effects of various error sources on the overall uncertainty are evaluated using an error simulation. Error sources modeled include boundarylayer/viscous interactions, pneumatic lag, thermal transpiration in the sensor pressure tubing, misalignment in the matrix layout, thermal warping of the vehicle nose, sampling resolution, and transducer error. Using simulated pressure data for input to the estimation algorithm, effects caused by various error sources are analyzed by comparing estimator outputs with the original trajectory. To obtain ensemble averages the simulation is run repeatedly and output statistics are compiled. Output errors resulting from the various error sources are presented as a function of Mach number. Final uncertainties with all modeled error sources included are presented as a function of Mach number.

Whitmore, Stephen A.↗
DOE OSTI2025

Multiscale Machine-Learned Modeling Infrastructure

The Multiscale Machine-Learned Modeling Infrastructure (MuMMI) is a multiscale workflow management infrastructure that can concurrently orchestrate thousands of molecular dynamics (MD) simulations operating at different time and/or length scales, spanning nanoseconds to seconds and nanometers to micrometers. MuMMI uses machine learning (backed by biology experiments) to guide a massive ensemble of MD simulations that capture biologically relevant time and length scales with unprecedented resolution. MuMMI supports multiple MD codes such as GROMACS and ddcMD and can be fully deployed using the HPC package manager Spack. MuMMI has been used in many publications to run hundreds of thousands simulations, leading to significant biology breakthroughs.

Di Natale, Francesco [Lawrence Livermore National ↗
DOE OSTI2021

An assessment of tropical cyclones in North American $\mathrm{CORDEX}$ $\mathrm{WRF}$ simulations

This work presents an assessment of tropical cyclones (TCs) in the 25 km and 50 km resolution reanalysis-forced and baseline and future (RCP8.5) global climate model (GCM) forced simulations produced for the North American branch of the international Coordinated Regional climate Downscaling Experiment (NA-CORDEX) using the Weather Research and Forecasting (WRF) model. A set of complementary 12 km resolution simulations produced as a part of a different project is also included in this assessment. Before examining the projections from the GCM-driven simulations, the ability of the simulations to minimally produce a realistic spatial distribution of historical TC occurrence was assessed in simulations forced by reanalysis and the three different GCMs used herein. Then, projections for occurrence, TC related mean precipitation and precipitation intensity, storm duration, the intensity measured by minimum pressure and maximum wind speed, storm size, and translation speed were examined. Several of these characteristics show little to no change in the future in trend or in distribution across the ensemble. However, many simulations suggest a westward shift or increase in TC occurrence over the East Pacific basin and a decrease in occurrence over the Caribbean and Gulf of Mexico. Increases (decreases) in total storm-related precipitation are projected where TC occurrence increases (decreases). TC precipitation intensity is found to increase in all simulations over the East Pacific, but projections are mixed over the North Atlantic. Finally, the ensemble projects a distribution shift towards more intense TC over the East Pacific, and a shift toward faster translation speeds over the North Atlantic.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

A Modeling Study of the Spring 2011 Extreme US Weather Activity

The spring of 2011 was characterized by record-breaking tornadic activity with substantial loss of life and destruction of property. While a waning La Nina and other atmospheric teleconnections have been implicated in the development of these extreme weather events, a quantitative assessment of their causes is still lacking. This study uses high resolution (1/4 lat/lon) GEOS-5 AGCM experiments to quantify the role of SSTs and soil moisture in the development of the extreme weather activity with a focus on April - the month of peak tornadic activity. The simulations, consisting of 22-member ensembles of three-month long simulations (initialized March 1st) reproduce the main features of the observed large-scale changes including the below-normal temperature and above-normal precipitation in the Central US, and the hot and dry conditions to the south. Various sensitivity experiments are conducted to separate the roles of the SST, soil moisture and the initial atmospheric conditions in the development and predictability of the atmospheric conditions (wind shear, moisture, etc.) favoring the severe weather activity and flooding.

Schubert, S.↗
NASA NTRS2015

Multimodel Ensembles of Wheat Growth: More Models are Better than One

Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop models can give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 24-38% for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

wheat↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2015

Multimodel Ensembles of Wheat Growth: Many Models are Better than One

Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop model scan give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 2438 for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

model intercomparison↗
DOE OSTI2020

Changes in Precipitation From North Atlantic Tropical Cyclones Under RCP Scenarios in the Variable-Resolution Community Atmosphere Model

Decreasing climate models' grid spacing improves the representation of tropical cyclones at decadal time scales. In this study, a variable-resolution (VR) version of the Community Atmosphere Model 5 (CAM5-VR) is utilized to study North Atlantic tropical cyclone climatology in ensemble historical climate simulations and under two Representative Concentration Pathway (RCP) projections (RCP4.5 and RCP8.5). Basin-wide tropical cyclone counts decrease in the RCP simulations, although landfalling storm counts do not show as straightforward of a pattern, especially when focusing on regional changes. Lifetime maximum intensity metrics suggest that tropical cyclones increase in strength in the RCP ensembles. However, despite increases in tropical cyclone-related precipitation rates and the amount of precipitation produced per storm with warming, the annual average Rx5day from tropical cyclones over the eastern United States decreases due to less landfalling storms. This work is part of a continued effort to quantify how tropical cyclone-induced hazards may change in future climates.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗
NASA NTRS2001

A Regional Model Study of Synoptic Features Over West Africa

Synoptic weather features over West Africa were studied in simulations by the regional simulation model (RM) at the NASA/Goddard Institute for Space Studies. These pioneering simulations represent the beginning of an effort to adapt regional models for weather and climate prediction over West Africa. The RM uses a cartesian grid with 50 km horizontal resolution and fifteen vertical levels. An ensemble of four simulations was forced with lateral boundary conditions from ECMWF global analyses for the period 8-22 August 1988. The simulated mid-tropospheric circulation includes the skillful development and movement of several African wave disturbances. Wavelet analysis of mid-tropospheric winds detected a dominant periodicity of about 4 days and a secondary periodicity of 5-8 days. Spatial distributions of RM precipitation and precipitation time series were validated against daily rain gauge measurements and ISCCP satellite infrared cloud imagery. The time-space distribution of simulated precipitation was made more realistic by combining the ECMWR initial conditions with a 24-hr spin-up of the moisture field and also by damping high frequency gravity waves by dynamic initialization. Model precipitation "forecasts" over the Central Sahel were correlated with observations for about three days, but reinitializing with observed data on day 5 resulted in a dramatic improvement in the precipitation validation over the remaining 9 days. Results imply that information via the lateral boundary conditions is not always sufficient to minimize departures between simulated and actual precipitation patterns for more than several days. In addition, there was some evidence that the new initialization may increase the simulations' sensitivity to the quality of lateral boundary conditions.

Druyan, Leonard M.↗
DOE OSTI2023

Historical and future trends in South Asian monsoon low pressure systems in a high-resolution model ensemble

Historical trends in monsoon low pressure systems (LPS), the dominant rain-bearing weather system of South Asia, have been difficult to assess due to changes in the observing network. Future projections have also remained uncertain because prior studies concluded that many coarse-resolution climate models do not accurately simulate LPS. Here, we examine changes in South Asian monsoon LPS simulated by an ensemble of global models, including some with high spatial resolution, that we show skillfully represent LPS. In the ensemble mean, the number of strong LPS (monsoon depressions) decreased over the last 65 years (1950–2014) by about 15% while no trend was detected for weaker LPS (monsoon lows). The reduction in depression counts then moderated, yielding no trend in the periods 1980–2050 or 2015–2050. The ensemble mean projects a shift in genesis from ocean to land and an increase in LPS precipitation of at least 7% K −1 , which together contribute to a projected increase in seasonal mean and extreme precipitation over central India.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Timescales of Soil Moisture Anomalies: Results from Two GCMs

Soil moisture anomalies dissipate over timescales that may span weeks to months. Characterizing the geographical and seasonal variations in these timescales can have important practical benefit; significant soil moisture "memory" allows long-lead forecasts of soil moisture, which have been found in recent studies to be essential for useful Ion--lead forecasts of precipitation in many regions. In this talk, we will present and compare the soil moisture timescales derived in two separate general circulation model (GCM) studies. Both studies employ multiple ensembles of short-term climate simulations. Timescales at a given point are effectively estimated by determining how quickly the soil moisture distribution generated in one ensemble of simulations (characterized by a unique set of initial soil moisture conditions) approaches that produced by another ensemble (characterized by a different set of initial soil moisture conditions). The talk will include a discussion of why the timescales produced by the two GCMs differ in some regions, and it will also describe the impact of soil moisture memory on simulated precipitation.

Koster, Randal D.↗