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 505 records · Page 28

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗

Nearest neighbour distributions: New statistical measures for cosmological clustering

ABSTRACT The use of summary statistics beyond the two-point correlation function to analyse the non-Gaussian clustering on small scales, and thereby, increasing the sensitivity to the underlying cosmological parameters, is an active field of research in cosmology. In this paper, we explore a set of new summary statistics – the k-Nearest Neighbour Cumulative Distribution Functions (kNN-CDF). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the k-nearest data points, and is sensitive to all connected N-point correlations in the data. The kNN-CDF can be used to measure counts in cell, void probability distributions, and higher N-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets – both discrete points, and the generalization for continuous fields. We use data from a large suite of N-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of kNN-CDF improves the constraints on the cosmological parameters by more than a factor of 2 when applied to the clustering of dark matter in the range of scales between 10 and $40\, h^{-1}\, {\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of haloes in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the kNN-CDF are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data, provided the higher order correlation functions are sensitive to cosmology on the scales of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity of the Atmospheric Response to Warm Pool El Nino Events to Modeled SSTs and Future Climate Forcings

Warm pool El Nino (WPEN) events are characterized by positive sea surface temperature (SST) anomalies in the central equatorial Pacific. Under present-day climate conditions, WPEN events generate poleward propagating wavetrains and enhance midlatitude planetary wave activity, weakening the stratospheric polar vortices. The late 21st century extratropical atmospheric response to WPEN events is investigated using the Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM), version 2. GEOSCCM simulations are forced by projected late 21st century concentrations of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) and by SSTs and sea ice concentrations from an existing ocean-atmosphere simulation. Despite known ocean-atmosphere model biases, the prescribed SST fields represent a best estimate of the structure of late 21st century WPEN events. The future Arctic vortex response is qualitatively similar to that observed in recent decades but is weaker in late winter. This response reflects the weaker SST forcing in the Nino 3.4 region and subsequently weaker Northern Hemisphere tropospheric teleconnections. The Antarctic stratosphere does not respond to WPEN events in a future climate, reflecting a change in tropospheric teleconnections: The meridional wavetrain weakens while a more zonal wavetrain originates near Australia. Sensitivity simulations show that a strong poleward wavetrain response to WPEN requires a strengthening and southeastward extension of the South Pacific Convergence Zone; this feature is not captured by the late 21st century modeled SSTs. Expected future increases in GHGs and decreases in ODSs do not affect the polar stratospheric responses to WPEN.

stratosphere↗

Integrating Growth Stage Deficit Irrigation into a Process Based Crop Model

Current rates of agricultural water use are unsustainable in many regions, creating an urgent need to identify improved irrigation strategies for water limited areas. Crop models can be used to quantify plant water requirements, predict the impact of water shortages on yield, and calculate water productivity (WP) to link water availability and crop yields for economic analyses. Many simulations of crop growth and development, especially in regional and global assessments, rely on automatic irrigation algorithms to estimate irrigation dates and amounts. However, these algorithms are not well suited for water limited regions because they have simplistic irrigation rules, such as a single soil-moisture based threshold, and assume unlimited water. To address this constraint, a new modeling framework to simulate agricultural production in water limited areas was developed. The framework consists of a new automatic irrigation algorithm for the simulation of growth stage based deficit irrigation under limited seasonal water availability; and optimization of growth stage specific parameters. The new automatic irrigation algorithm was used to simulate maize and soybean in Gainesville, Florida, and first used to evaluate the sensitivity of maize and soybean simulations to irrigation at different growth stages and then to test the hypothesis that water productivity calculated using simplistic irrigation rules underestimates WP. In the first experiment, the effect of irrigating at specific growth stages on yield and irrigation water use efficiency (IWUE) in maize and soybean was evaluated. In the reproductive stages, IWUE tended to be higher than in the vegetative stages (e.g. IWUE was 18% higher than the well watered treatment when irrigating only during R3 in soybean), and when rainfall events were less frequent. In the second experiment, water productivity (WP) was significantly greater with optimized irrigation schedules compared to non-optimized irrigation schedules in water restricted scenarios. For example, the mean WP across 38 years of maize production was 1.1 kg/cu m for non-optimized irrigation schedules with 50 mm of seasonal available water and 2.1 kg/cu m optimized ion schedules, a 91% improvement in WP with optimized irrigation schedules. The framework described in this work could be used to estimate WP for regional to global assessments, as well as derive location specific irrigation guidance.

crop model↗

Elimination LArTPC Simulation Uncertainty With Modifcations to TPC Wire Waveforms

This paper introduces a novel approach to reducing systematic uncertainties in LArTPC detector simulations, which arise from known detector effects that are challenging and computationally intensive to model. The proposed method involves modifying simulated waveforms based on observed discrepancies in ionization signals from the TPC between actual data and simulations. This approach reduces the sensitivity to detailed detector modeling and decreases the computational resources required for accurate simulations.

Mktchyan, Karen↗

Comparison of Commonly-Used Microwave Radiative Transfer Models for Snow Remote Sensing

This paper reviews four commonly-used microwave radiative transfer models that take different electromagnetic approaches to simulate snow brightness temperature (T(sub B)): the Dense Media Radiative Transfer - Multi-Layer model (DMRT-ML), the Dense Media Radiative Transfer - Quasi-Crystalline Approximation Mie scattering of Sticky spheres (DMRT-QMS), the Helsinki University of Technology n-Layers model (HUT-nlayers) and the Microwave Emission Model of Layered Snowpacks (MEMLS). Using the same extensively measured physical snowpack properties, we compared the simulated T(sub B) at 11, 19 and 37 GHz from these four models. The analysis focuses on the impact of using different types of measured snow microstructure metrics in the simulations. In addition to density, snow microstructure is defined for each snow layer by grain optical diameter (Do) and stickiness for DMRT-ML and DMRT-QMS, mean grain geometrical maximum extent (D(sub max)) for HUT n-layers and the exponential correlation length for MEMLS. These metrics were derived from either in-situ measurements of snow specific surface area (SSA) or macrophotos of grain sizes (D(sub max)), assuming non-sticky spheres for the DMRT models. Simulated T(sub B) sensitivity analysis using the same inputs shows relatively consistent T(sub B) behavior as a function of Do and density variations for the vertical polarization (maximum deviation of 18 K and 27 K, respectively), while some divergences appear in simulated variations for the polarization ratio (PR). Comparisons with ground based radiometric measurements show that the simulations based on snow SSA measurements have to be scaled with a model-specific factor of Do in order to minimize the root mean square error (RMSE) between measured and simulated T(sub B). Results using in-situ grain size measurements (SSA or D(sub max), depending on the model) give a mean T(sub B) RMSE (19 and 37 GHz) of the order of 16-26 K, which is similar for all models when the snow microstructure metrics are scaled. However, the MEMLS model converges to better results when driven by the correlation length estimated from in-situ SSA measurements rather than D(sub max) measurements. On a practical level, this paper shows that the SSA parameter, a snow property that is easy to retrieve in-situ, appears to be the most relevant parameter for characterizing snow microstructure, despite the need for a scaling factor.

Snow↗

The Ozone Budget in the Upper Troposphere from Global Modeling Initiative (GMI)Simulations

Ozone concentrations in the upper troposphere are influenced by in-situ production, long-range tropospheric transport, and influx of stratospheric ozone, as well as by photochemical removal. Since ozone is an important greenhouse gas in this region, it is particularly important to understand how it will respond to changes in anthropogenic emissions and changes in stratospheric ozone fluxes.. This response will be determined by the relative balance of the different production, loss and transport processes. Ozone concentrations calculated by models will differ depending on the adopted meteorological fields, their chemical scheme, anthropogenic emissions, and treatment of the stratospheric influx. We performed simulations using the chemical-transport model from the Global Modeling Initiative (GMI) with meteorological fields from (It)h e NASA Goddard Institute for Space Studies (GISS) general circulation model (GCM), (2) the atmospheric GCM from NASA's Global Modeling and Assimilation Office(GMAO), and (3) assimilated winds from GMAO . These simulations adopt the same chemical mechanism and emissions, and adopt the Synthetic Ozone (SYNOZ) approach for treating the influx of stratospheric ozone -. In addition, we also performed simulations for a coupled troposphere-stratosphere model with a subset of the same winds. Simulations were done for both 4degx5deg and 2degx2.5deg resolution. Model results are being tested through comparison with a suite of atmospheric observations. In this presentation, we diagnose the ozone budget in the upper troposphere utilizing the suite of GMI simulations, to address the sensitivity of this budget to: a) the different meteorological fields used; b) the adoption of the SYNOZ boundary condition versus inclusion of a full stratosphere; c) model horizontal resolution. Model results are compared to observations to determine biases in particular simulations; by examining these comparisons in conjunction with the derived budgets, we may pinpoint deficiencies in the representation of chemical/dynamical processes.

Rodriquez, J.↗

Sensitivity Analysis of the Human Research Program’s Impact 1.0 Model

Sensitivity analysis estimates the relative contribution of the uncertainty in input values to the uncertainty of model outputs. Partial Rank Correlation Coefficient (PRCC) and Leave One Out (LOO) are methods of conducting sensitivity analysis on non-linear simulation models like the IMPACT Model. The PRCC method estimates the sensitivity using partial correlation of the ranks of the generated input values to each generated output value. The “partial” part is due to the adjustments made for the linear effects of all the other input values in the calculation of correlation between a particular input and each output. LOO removes a medical condition from test suite and calculates the change in the outcome. This is used to identify the influence of the condition input in the model The inputs to the sensitivity procedures include the number of occurrences of each of the one hundred plus IMPACT medical conditions generated over the simulations, and the IMPACT outputs from the MEDPRAT mathematical model. The outputs total task time lost (TTL), number of return to definitive care (RTDC), and number of loss of crew lives (LOCL). The IMPACT team will report the results of using PRCC and LOO on IMPACT 1.0. Tornado plots will assist in the visualization of the condition-related input sensitivities to each of the main outcomes. The outcomes of this sensitivity analysis will drive review focus by identifying conditions where changes in uncertainty, and input values could drive changes in overall model output uncertainty. These efforts are an integral part of the overall verification, validation, and credibility review of IMPACT 1.0.

Sensitivity Anaylsis↗

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗

Satellite Hyperspectral Imaging Simulation

Simulation of generic pushbroom satellite hyperspectral sensors have been performed to evaluate the potential performance and validation techniques for satellite systems such as COIS(NEMO), Warfighter-1(OrbView-4) and Hyperion(EO-1). The simulations start with a generation of synthetic scenes from material maps of studied terrain. Scene-reflected radiance is corrected for atmospheric effects and convolved with sensor spectral response using MODTRAN 4 radiance and transmissions calculations. Scene images are further convolved with point spread functions derived from Optical Transfer Functions (OTF's) of the sensor system. Photon noise and etectorr/electronics noise are added to the simulated images, which are also finally quantized to the sensor bit resolution. Studied scenes include bridges and straight roads used for evaluation of sensor spatial resolution, as well as fields of minerals, vegetation and manmade materials used for evaluation of sensor radiometric response and sensitivity. The scenes are simulated with various seasons and weather conditions. Signal-to-noise ratios and expected performance are estimated for typical satellite system specifications and are discussed for all the scenes.

Zanoni, Vicki↗

Satellite Hyperspectral Imaging Simulation

Simulations of generic pushbroom satellite hyperspetral sensors have been performed to evaluate the potential performance and validation techniques for satellite systems such as COIS (NEMO), Warfighter-1 (OrbView-4), and Hyperion (EO-1). The simulaitons start with a generation of synthetic scenes form material maps of studied terrain. Scene-reflected radiance is corrected for atmospheric effects and convolved with sensor spectral response uwing MODTRAN 4 radiance and transmission calculations. Scene images are further convolved with point spread functions derived from Optical Transfer Functions (OTF's) of the sensor system. Photon noise and detector/electronics noise are added to the simulated images, which are also finally quantized to the sensor bit resolution. Studied scenes include bridges and straight roads used for evaluation of sensor spatial resolution, as well as fields of minerals, vegetation, and manmade materials used for evaluation of sensor radiometric response and sensitivity. The scenes are simulated with various seasons and weather conditions. Signal-to-noise ratos and expected performance are estimated for typica satellite system specifications and are discussed for all the scenes.

Zanoni, Vicki↗

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗

Failure Bounding And Sensitivity Analysis Applied To Monte Carlo Entry, Descent, And Landing Simulations

In the study of entry, descent, and landing, Monte Carlo sampling methods are often employed to study the uncertainty in the designed trajectory. The large number of uncertain inputs and outputs, coupled with complicated non-linear models, can make interpretation of the results difficult. Three methods that provide statistical insights are applied to an entry, descent, and landing simulation. The advantages and disadvantages of each method are discussed in terms of the insights gained versus the computational cost. The first method investigated was failure domain bounding which aims to reduce the computational cost of assessing the failure probability. Next a variance-based sensitivity analysis was studied for the ability to identify which input variable uncertainty has the greatest impact on the uncertainty of an output. Finally, probabilistic sensitivity analysis is used to calculate certain sensitivities at a reduced computational cost. These methods produce valuable information that identifies critical mission parameters and needs for new technology, but generally at a significant computational cost.

Gaebler, John A.↗

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗

Experimental demonstration and modeling of a robotic neutron detector with spectral and directional sensitivity for treaty verification

International safeguards and arms control agreements often require labor-intensive, intrusive onsite inspections to perform verification tasks. The ability to localize a neutron source and/or characterize a neutron field may be imperative for identifying anomalies. We are interested in the role of autonomous mobile robots, which, if designed properly, may be more effective and efficient and less intrusive than their human counterparts. Toward developing such a capability, we, for this study, have constructed the N-SpecDir Bot, comprised of three boron-coated straw detectors azimuthally-distributed within a cylinder of high-density polyethylene, which is mounted on an omni-directional robotic platform. Our N-SpecDir Bot is specifically designed to provide spectral and directional sensitivity, in addition to gross counts, by utilizing the signals from the three detectors. The detection system has been extensively characterized by MCNP modeling, which has been benchmarked to experiments conducted at the Princeton Plasma Physics Laboratory. We demonstrate the spectral and directional sensitivity experimentally and in simulation, and provide a simple yet robust method for directional measurements.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 2: Kinetic energy and resolution analysis

A kinetic energy (KE) analysis of the forcing of a mesoscale upper-tropospheric jet streak by organized diabatic processes within the simulated convective system (SCS) that was discussed in Part 1 is presented in this study. The relative contributions of the ageostrophic components of motion to the generation of KE of the convectively generated jet streak are compared, along with the KE generation by the rotational (nondivergent) and irrotational (divergent) mass transport. The sensitivity of the numerical simulations of SCS development to resolution is also briefly examined. Analysis within isentropic coordinates provides for an explicit determination of the influence of the diabatic processes on the generation of KE. The upper-level production of specific KE is due predominatly to the inertial advective ageostrophic component (IAD), and as such represents the primary process through which the KE of the convectively generated jet streak is realized. A secondary contribution by the inertial diabatic (IDI) term is observed. Partitioning the KE generation into its rotational and irrotational components reveals that the latter, which is directly linked to the diabatic heating within the SCS through isentropic continuity requirements, is the ultimate source of KE generation as the global area integral of generation by the rotational component vanishes. Comparison with an identical dry simulation reveals that the net generation of KE must be attributed to latent heating. Both the IAD and IDI ageostrophic components play important roles in this regard. Examination of results from simulations conducted at several resolutions supports the previous findings in that the effects of diabatic processes and ageostrophic motion on KE generation remain consistent. Resolution does impact the location and timing of SCS development, a result that has important implications in forecasting the onset of convection that develops from evolution of the large-scale flow and moisture transport. Marked differences are observed in the momentum field aloft subsequent to the life cycle of the SCS in the 1 deg, 30-level base case (MP130) simulation discussed in Part 1 versus its 2 deg counterparts in that the MP130 simulation with higher spatial resolution contains from 14% to 30% more total KE.

Wolf, Bart J.↗

Sensitivity of the dynamic-shell target to laser drive nonuniformities

The dynamic-shell concept for inertial confinement fusion (ICF) uses an initially homogeneous target and a carefully shaped laser pulse to form a shell and implode it. The laser pulse consists of a series of pickets that drive shocks into the target. The first few shocks converge inwards and rebound from the center of the target, creating an expanding, low-density plasma. Subsequent shocks are launched into the expanding plasma and eventually coalesce to form a shell, which is then imploded with a traditional ICF laser pulse. This study describes radiation-hydrodynamic simulations that investigate the sensitivity of dynamic-shell targets to imperfections in the laser drive. A one-dimensional (1D) study looks at mistiming and power variations in the pickets and a two-dimensional (2D) study examines irradiation perturbations imposed by the laser-beam geometry. Simulations show that less than ~2% power imbalance or 200 ps timing variation in the pickets is sufficient to keep the yield above 90% of the maximum. Additionally, the 2D simulations show that 72 or more beams are required to keep irradiation nonuniformities low enough to obtain fusion yields close to that of 1D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elucidating Norrish type I reactive pathways by ultrafast X-ray absorption spectroscopy

Norrish type I reactions selectively cleave carbon–carbon bonds directly adjacent to carbonyl groups. Despite their broad use in combination with aromatic carbonyls for additive manufacturing and dental UV curing applications, the nature of the photochemically active state and its population mechanism remain insufficiently understood. Detailed mechanistic insight requires mapping of the photoexcited population flow involving internal conversion and intersystem crossing. We present a time-domain study of gas phase acetophenone as a prototypical aromatic carbonyl combining soft X-ray time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy at the oxygen K-edge with ab initio multiple spawning (AIMS) simulations. Exploiting the specific sensitivity of TR-NEXAFS spectroscopy to states with nπ* character, we observe population transfer from the initially excited 1 ππ* state to the 1 nπ* state with a time constant of (0.13 ± 0.02) ps after an initial induction period of (0.12 ± 0.02) ps without population transfer, in quantitative agreement with the AIMS simulations. The population in the 1 nπ* state subsequently decays via intersystem crossing, likely mediated by a 3 ππ* state, within (3.17 ± 0.66) ps to a long-lived 3 nπ* state, which is presumed to be active towards Norrish type I chemistry.

Graßl, Martin [SLAC National Accelerator Laborator↗