Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity 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 1,099 records · Page 61

Development of a High Resolution Liquid Xenon Imaging Telescope for Medium Energy Gamma Ray Astrophysics

In the third year of the research project, we have (1) tested a 3.5 liter prototype of the Liquid Xenon Time Projection Chamber, (2) used a prototype having a 4.4 cm drift gap to study the charge and energy resolution response of the 3.5 liter chamber, (3) obtained an energy resolution as good as that previously measured by us using chambers with drift gaps of the order of millimeters, (4) observed the induction signals produced by MeV gamma rays, (4) used the 20 hybrid charge sensitive preamplifiers for a nondestructive readout of the electron image on the induction wires, (5) performed extensive Monte Carlo simulations to obtain results on efficiency, background rejection capability, and source flux sensitivity, and (6) developed a reconstruction algorithm for events with multiple interaction points.

Aprile, Elena↗

Quench Simulation and Validation for EIC Direct Wind Magnets

The Electron-Ion Collider (EIC), a new facility to be built in the United States at the U.S. Department of Energy’s Brookhaven National Laboratory in collaboration with Thomas Jefferson National Accelerator Facility-will be the only new particle collider in the U.S. for the next 10∼20 years. The EIC project will investigate the origin of mass and building blocks of the Universe. It will accelerate and collide beams of electrons with protons and nuclei. Unlike any collider that has ever existed before, both EIC’s particle beams will be polarized, meaning their spins will be aligned in a controlled way to allow for precision measurements. Different types of new superconducting magnets will be required near the interaction region (IR) of EIC. A 6-around-1 NbTi cable will be used for some of those new superconducting magnets. Given the relatively smaller conductors used for those direct wind magnets, the overall engineering current density will be 2–3 times higher in comparison to traditional collared superconducting magnets in other accelerator magnets. This is technically challenging and potentially represents high risks for quench protection. Here, in this paper we use different tools and methods for quench simulation for several EIC direct wind magnets. We also present results of a sensitivity study, exploring different material properties and Jc curves. This paper also discusses and validates the quench simulation results with measurement results from a small-scale serpentine magnet tested at BNL in October 2024.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.

Computational Fluid Dynamics↗

Direct Impact Split-Hopkinson Bar Experiment for Compression Testing at Strain Rate of 50,000 S-1

A direct impact Hopkinson compression bar experiment is used for compression testing at strain rates in the order of 50,000 s-1. A short specimen mounted on the cross-section surface of a transmitter bar is impacted directly by a projectile (the incident bar in the classical split Hopkinson bar is eliminated).The quantities measured in the test are the projectile impact velocity, the force in the transmitter bar and the strain on the surface of the specimen, which is measured using the digital image correlation method with high-speed cameras. The deformation measurements show a nonuniform deformation in the specimen. To study the strain rate sensitivity properties of the specimen’s material, the experiment is simulated using an assumed material model properties which is validated by comparing the calculated and measured quantities. For aluminum 2024-T351 the results show a small strain rate sensitivity (increase of stress with strain rate) at strain rates above 10,000 s-1.

Nathan Spulak↗

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗

Subdiffusive High-Pressure Hydrogen Gas Dynamics in Elastomers

We report elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene–propylene–diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.

36 MATERIALS SCIENCE↗

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pliocene Model Intercomparison Project Phase 3 (PlioMIP3) – Science Plan and Experimental Design

The Pliocene Model Intercomparison Project (PlioMIP) was initiated in 2008. Over two phases PlioMIP has helped co-ordinate the experimental design and publication strategy of the community, which has included an increasing number of climate models and modelling groups from around the world. It has engaged with palaeoenvironmental scientists to foster new data synthesis supporting the construction of new model boundary conditions, as well as to facilitate new data-model comparisons. The work has advanced our understanding of Pliocene climates and environments, enhanced our knowledge regarding the ability of complex climate and Earth System models to accurately simulate climate change, and helped to refine our estimates of how sensitive the climate system is to forcing conditions. In this community protocol paper, we outline the scientific plan for PlioMIP Phase 3 (PlioMIP3). This plan provides the required guidance to participating modelling groups from around the world to successfully set up and perform PlioMIP3 climate model experiments. The project is open to new participants from the scientific community (both from the climate modelling and geosciences communities). In PlioMIP3, we retain the PlioMIP2 Core experiments (Eoi 400 , E 280 ) and extend the Core requirements to include either an experiment focussed on the Early Pliocene or an alternative Late Pliocene simulation (or both). These additions (a) allow a comparison of Early and Late Pliocene warm intervals and help build research connections and synergy with the MioMIP (Miocene Model Intercomparison Project - also known as DeepMIP-Miocene) and PlioMioVAR projects (Pliocene-Miocene Variability Working Group), and (b) create an alternative time slice simulation for 3.205 Ma (MIS KM5c) through removal of some of the largest palaeogeographic differences introduced between PlioMIP1 and 2 resulting in minimal land-sea mask variations from the modern. In addition, we present ten optional experiments designed to enhance our assessment of climate sensitivity and to explore the uncertainty in greenhouse gas-related forcing. For the first time, we introduce orbital sensitivity experiments into the science plan, as well as simulations incorporating dynamic vegetation-climate feedbacks and an experiment designed to examine the potential significance of East Antarctic Ice Sheet boundary condition uncertainty. These changes enhance palaeo-to-future scientific connections and enable an exploration of the significance of palaeogeographic uncertainties on climate simulations.

Climate modelling↗

Files and scripts to support manuscript Needham et al. Canopy Gradients of Respiration

This dataset includes the parameter files, relevant output files, and scripts to perform analysis with Jupyter notebooks that support the manuscript Needham et al 2025 “Canopy Gradients of Respiration Drive Plant Carbon Budgets and Leaf Area Index.” We add functionality to the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) to allow flexible vertical gradients of leaf maintenance respiration (Rdark) and maximum carboxylation rate (Vcmax) through the canopy. We test the sensitivity of FATES to canopy gradients in Rdark, both in global simulations to assess broad scale impacts on leaf area index (LAI) and vegetation carbon, and in single site simulations where we assess impacts on plant functional type (PFT) competitive dynamics. Parameter files are netcdf files that can be converted to human readable .cdl files using NCO tools. Analysis scripts are Jupyter notebook files. These can be opened and run using the open source Jupyter notebook software. Model outputs are netcdf files.

54 ENVIRONMENTAL SCIENCES↗

Diurnal Variation of Precipitation During MC3E Campaign

The diurnal variation of precipitation processes in the United States (US) is well recognized but incompletely understood (Cabone et al. 2002). The diurnal cycle of precipitation has been studied using surface rainfall data, radar reflectivity data, and satellite-derived cloudiness and precipitation (Wallace 1975; Dai et al. 1999; Carbone et al. 2002; Carbone and Tuttle, 2008; Parker and Ahijevych, 2007; Matsui et al. 2010 and others). These observations indicate that the summer-time precipitation most of the North America and typically feature late-afternoon precipitation maxima. These diurnal variation of precipitation can also be generally categorized into three different types: 1) afternoon rainfall maxima due to mesoscale and local circulations over the south and east of the Mississippi and Ohio valleys, 2) nocturnal rainfall maxima from eastward-propagating mesoscale convective systems (MCSs) over the Lee side of Rocky Mountain regions and 3) afternoon rainfall maxima in the Appalachian Mountains, and then propagate eastward toward the coast. The main objective of this paper is to use a regional cloud-scale model with very high-resolution (i.e., WRF) to examine the WRF ability to simulate diurnal variation of precipitation. Specifically, the study will (1) identify the physical processes responsible to diurnal variation of precipitation, (2) examine the sensitivity of resolution (2, 6, 18, and 30 km) to model simulated diurnal variation of precipitation and (3) identify the relationships between microphysics and cumulus parameterization schemes.

Tao, Wei-Kuo↗

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗

Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

54 ENVIRONMENTAL SCIENCES↗

Natural variability contributes to model–satellite differences in tropical tropospheric warming

Significance Climate models have, on average, simulated substantially more tropical tropospheric warming than satellite data, with few simulations matching observations. It has been suggested that this discrepancy arises because climate models are overly sensitive to greenhouse gas increases. Tropical tropospheric temperature trends from a large ensemble of simulations performed with a single climate model span a wide range that is solely due to natural climate variability. A subset of these simulations have warming rates in accord with satellite observations. Simulations with diminished tropical tropospheric warming due to climate variability exhibit characteristic patterns of surface warming that are similar to the observed record. Our results indicate that multidecadal variability can explain current model–observational differences in the rate of tropical tropospheric warming.

58 GEOSCIENCES↗

New tree-ring data from Canadian boreal and hemi-boreal forests provide insight for improving the climate sensitivity of terrestrial biosphere models

Understanding boreal/hemi-boreal forest growth sensitivity to seasonal variations in temperature and water availability provides important basis for projecting the potential impacts of climate change on the productivity of these ecosystems. Our best available information currently comes from a limited number of field experiments and terrestrial biosphere model (TBM) simulations of varying predictive accuracy. Here, we assessed the sensitivity of annual boreal/hemi-boreal forest growth in Canada to yearly fluctuations in seasonal climate variables using a large tree-ring dataset and compared this to the climate sensitivity of annual net primary productivity (NPP) estimates obtained from fourteen TBMs. We found that boreal/hemi-boreal forest growth sensitivity to fluctuations in seasonal temperature and precipitation variables changed along a southwestern to northeastern gradient, with growth limited almost entirely by temperature in the northeast and west and by water availability in the southwest. We also found a lag in growth climate sensitivity, with growth largely determined by the climate during the summer prior to ring formation. Analyses of NPP sensitivity to the same climate variables produced a similar southwest to northeast gradient in growth climate sensitivity for NPP estimates from all but three TBMs. However, analyses of growth from tree-ring data and analyses of NPP from TBMs produced contrasting evidence concerning the key climate variables limiting growth. While analyses of NPP primarily indicated a positive relationship between growth and seasonal temperature, tree-ring analyses indicated negative growth relationships to temperature. Also, the positive effect of precipitation on NPP derived from most TBMs was weaker than the positive effect of precipitation on tree-ring based growth: temperature had a more important limiting effect on NPP than tree-ring data indicated. These mismatches regarding the key climate variables limiting growth suggested that characterization of tree growth in TBMs might need revision, particularly regarding the effects of stomatal conductance and carbohydrate reserve dynamics.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulation of Infrasound Resonance in Underground Tunnel Structures

Remote observation of infrasound resonant signals emanating from underground tunnel structures could potentially allow the remote quantification of the geometry of the underground structures in which the signals were generated. However, the sensitivity of these observations to tunnel geometry and changes in that geometry are unknown. In this report we outline a numerical simulation study with the following three objectives: 1. Can we model infrasound resonance using numerical simulations? 2. What is the sensitivity of the modeled observations to changes in tunnel geometry or boundary conditions? And 3. Can we accurately model resonant observations of explosions occurring in Redmond Salt Mine in Central Utah. In this report we outline affirmative answers to the first two objectives, but we were unable to accurately model the Redmond Explosions due to numerical instability in the model of the complex structure of the mine. Recommendations for future work emphasize the need to acquire additional datasets and to explore a more data-based approach in which changes in data signatures are detected as a first step towards developing a method which inverts resonant infrasound signals for tunnel geometry.

47 OTHER INSTRUMENTATION↗

The latitudinal gradient of the NO peak density

Results are presented from SME observations of the latitudinal gradients of peak NO densities at about 110-km altitude during the solstice and equinox periods from 1982 through 1985. It is shown that the response of the peak NO densities to the declining level of solar activity varies with latitude, with the polar regions exhibiting low sensitivity and the low-latitude regions responding strongly. The SME data also revealed marked asymmetries in the latitudinal structure of the two hemispheres for each season and considerable day-to-day variations in the NO densities. The solar cycle minimum data for June were simulated using a two-dimensional model; results of sensitivity studies performed with varied quenching rate and eddy diffusion coefficient are presented.

Fesen, C. G.↗

Radiative, Hydrologic, and Circulation Responses to Warming in Cess‐Potter Simulations Using the Global 3.25‐km SCREAM

Using the global 3.25-km Simple Cloud Resolving E3SM Atmosphere Model (SCREAM 3 km), a pair of 13-month Cess-Potter simulations are performed to quantify the radiative feedbacks and the hydrologic and circulation responses to warming. Large-scale aspects of SCREAM 3 km's top-of-atmosphere radiative fluxes, precipitation rates, and circulations are in good agreement with observations and reanalysis, with notable differences, including a drier lower free-troposphere in the Tropics, reduced precipitation and humidity over the Tropical West Pacific, and poleward shifted Southern Hemisphere midlatitude jet. In response to warming, SCREAM 3 km predicts a total radiative feedback within the top 15% of the CMIP5 and CMIP6 models, which puts it substantially higher than the feedback reported by other kilometer-scale models. SCREAM 3 km's high radiative feedback stems from a strongly positive shortwave cloud feedback, most prominent over the mid- and high-latitudes. SCREAM 3 km's high precipitation response also puts it among the highest of CMIP models, whereas its circulation response is within the spread of CMIP models. An ensemble of five perturbed initial condition Cess-Potter simulations with a 12 km version of SCREAM (SCREAM 12 km) is performed to characterize uncertainty and resolution sensitivity. It suggests that the uncertainty from analyzing a pair of 1-year simulations is small compared to the inter-model spread in feedbacks and precipitation response. SCREAM 12 km also produces a strong precipitation response to warming but a much lower cloud feedback and total radiative feedback. The results from these experiments suggest that the spread in climate feedbacks will likely persist in the next generation of kilometer-scale models.

54 ENVIRONMENTAL SCIENCES↗