Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global sensitivity analysis”

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 163 records · Page 9

A Neural Network Parametrization of Volumetric Cloud Fraction Profiles Using Satellite Observations and MERRA-2 Reanalysis Meteorological Data

Clouds play a crucial role in regulating the hydrologic cycle and Earth's radiative energy budget, yet they are often poorly represented in global climate models (GCMs). This study applies deep machine learning techniques to develop a physical parameterization of volumetric cloud fraction (VCF), the fraction of a 3-D grid volume occupied by clouds using satellite lidar-radar measurements. The neural network (NN) captures the complicated relationships between observed VCF profiles and collocated meteorological variables from MERRA-2 reanalysis data. Our results show that the NN model, particularly a sequence-to-sequence long short-term memory (LSTM) network with a sixfactor loss function, effectively learns the underlying cloud physical processes. The NN model outperforms MERRA-2 reanalysis in representing low-level clouds in tropical and subtropical regions and low- and middle-level clouds over midlatitude storm-track regions, and also improves VCF histograms. These improvements are reflected in the vertical distributions of zonally, meridionally, and globally averaged VCFs, geographic distributions of low-, middle-, and high-level clouds, and seasonal variations in monthly-mean VCF. Furthermore, the NN predictions effectively capture the El Niño-Southern Oscillation (ENSO) effects and other interannual variations. The NN parameterization is further evaluated through a sensitivity analysis, in which a single predictor is perturbed at a time. This reveals that relative humidity (RH) is the dominant factor influencing variations in globally averaged VCF at low and middle altitudes, followed by temperature. At higher altitudes, temperature becomes the primary driver of VCF through its effect on RH. Changes in wind components had minimal impact on globally averaged VCF.

Shan Zeng↗

Linearized ray-trace analysis

A new, coordinate-free version of the exact ray-trace equations for optical systems consisting of conic reflecting, refracting and reference surfaces is presented. These equations are differentiated to obtain closed-form optical sensitivity dyadics. For computation, the sensitivities are evaluated in a single global coordinate frame and combined in linearized ray-trace matrix difference equations that propagate the rays and the sensitivities from element to element. One purpose of this analysis is to create optical models that can be directly integrated with models of the instrument structure and control systems for dynamic simulation.

Redding, David C.↗

A Type II Hamiltonian Variational Principle and Adjoint Systems for Lie Groups

We present a novel Type II variational principle on the cotangent bundle of a Lie group which enforces Type II boundary conditions, i.e., fixed initial position and final momentum. In general, such Type II variational principles are only globally defined on vector spaces or locally defined on general manifolds; however, by left translation, we are able to define this variational principle globally on cotangent bundles of Lie groups. Type II boundary conditions are particularly important for adjoint sensitivity analysis, which is our motivating application. As such, we additionally discuss adjoint systems on Lie groups, their properties, and how they can be used to solve optimization problems subject to dynamics on Lie groups.

97 MATHEMATICS AND COMPUTING↗

An Integrated System for Water Budget Closure Over the Pan-Tropical Land Mass: A TRMM Validation Effort

A concerted effort was undertaken this past year to analyze the National Centers for Environmental Prediction (NCEP) and Department of Energy (DOE) Reanalysis II tropical hydrologic cycle. The topics included in this progress report are: (1) analysis of the surface hydrologic cycle; (2) long term CO2 sensitivity simulations with NCAR's CCM3 and global recycling rates; (3) tropical atmospheric hydrologic cycle; (4) GCIP (Global Continental Scale International Project) WEBS (Water and Energy Budget Synthesis); (5) the impact of precipitation uncertainties on global water balance calculations.

Roads, J.↗

Assimilation of Satellite Soil Moisture for Improved Atmospheric Reanalyses

A newly developed, weakly coupled land and atmosphere data assimilation system for NASA’s Global Earth Observing System model is presented, and used to demonstrate the benefit of assimilating satellite soil moisture into an atmospheric reanalysis. Specifically, Advanced Scatterometer and Soil Moisture Ocean Salinity soil moisture retrievals are assimilated into a system that uses the same model, atmospheric assimilation system, and atmospheric observations as the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2). The atmosphere is sensitive to soil moisture only under certain conditions. Hence, while the globally averaged model improvements were small, regionally, the soil moisture assimilation induced some substantial improvements. For example, in a large region spanning from western Europe across southern Russia, the soil moisture assimilation decreased the RMSE against independent station observations of daily maximum 2-m temperature (T(sup 2m, sub max)) by up to 0.4 K, and of 2-m specific humidity (q(exp 2m)) by up to 0.5 g/kg. Over all available stations, the mean T(sup 2m, sub max) RMSE was reduced from 2.82 to 2.79 K, while the mean q(exp 2m) RMSE was reduced from 1.25 to 1.20 g/kg. The soil moisture assimilation also reduced the mean RMSE across 29 flux tower sites from 34.2 to 32.6 Wm(exp −2) for latent heating, and from 37.7 to 36.5 Wm(exp −2) for sensible heating. For all variables evaluated, the soil moisture assimilation improved the model at monthly to seasonal, rather than daily, time scales. Based on the above experiments, it is recommended that satellite soil moisture be assimilated into future reanalyses, including the follow-on to MERRA-2.

Clara Draper↗

The use of NASA GEOS Global Analysis in MM5/WRF Initialization: Current Studies and Future Applications

An effort has been made at NASA/GSFC to use the Goddard Earth Observing system (GEOS) global analysis in generating the initial and boundary conditions for MM5/WRF simulation. This linkage between GEOS global analysis and MM5/WRF models has made possible for a few useful applications. As one of the sample studies, a series of MM5 simulations were conducted to test the sensitivity of initial and boundary conditions to MM5 simulated precipitation over the eastern; USA. Global analyses horn different operational centers (e.g., NCEP, ECMWF, I U ASA/GSFCj were used to provide first guess field and boundary conditions for MM5. Numerical simulations were performed for one- week period over the eastern coast areas of USA. the distribution and quantities of MM5 simulated precipitation were compared. Results will be presented in the workshop. In addition,other applications from recent and future studies will also be addressed.

Pu, Zhao-Xia↗

Sensitivity analysis for characterizing the impact of HNGD model on the prediction of hydrogen redistribution in Zircaloy cladding using BISON code

Hydrogen in zirconium cladding is able to precipitate into zirconium hydrides which impacts cladding integrity. The Hydride Nucleation-Growth-Dissolution (HNGD) model in the BISON code accounts for the precipitation and dissolution kinetics of hydride in Zircaloy material. This paper presents global sensitivity analyses of the HNGD model aiming to enhance our understanding of the hydride precipitation phenomena by quantifying the variance that key parameters have on the prediction of hydrogen behavior under various environmental conditions. Here, model predictions are compared to experimental data obtained under two different conditions: 1) with uniformly precharged specimens subjected to a linear thermal gradient, and 2) specimens precharged with a cathodically applied hydride rim at one end of the sample and subjected to an asymmetric thermal gradient. The Sobol sensitivity analysis identifies the key parameters in the HNGD model for both types of specimens. For linear temperature cases, the heat of transport dominates the accuracy of predictions when no precipitation occurs at the cold end, while Terminal Solid Solubility for Dissolution (TSSD) is the most important parameter when precipitation occurs. A large variation in the predicted hydrogen concentration profiles is found in the range of high TSSD due to the occurrence of precipitation. For asymmetric temperature cases, the solubility coefficient gives the largest impact on the predicted hydrogen distribution, as it determines the amount of solute hydrogen dissolved from the initially applied hydride rim. A large discrepancy in hydrogen distribution between simulations and experiments exists with the asymmetric specimens because BISON simulations fail to predict the precipitation of hydride at the cooler end. Comparative studies using former and updated models verifies the significant impact of the hydride growth mechanism on predicted hydrogen concentration profiles. In particular, when hydride initially exists, changes in TSSD generate a large variation in the predicted amount of precipitation by hydride growth, giving large uncertainty in predicting the hydrogen distribution over the sample length. The outputs characterize the significant impact of the hydride growth mechanism in the HNGD model on predicting hydrogen behavior, and improve the understanding of the precipitation of hydride in Zircaloy cladding within a range of expected environmental conditions. The analyses indicate work is still needed to improve the hydride solvus models in the BISON code to accurately predict experimentally observed hydride concentrations and distributions.

36 MATERIALS SCIENCE↗

The First Global 883 GHz Cloud Ice Survey: IceCube Level 1 Data Calibration, Processing and Analysis

Sub-millimeter (sub-mm, 200 - 1000 GHz) wavelengths contribute a unique capability to fill-in the sensitivity gap between operational visible/infrared (VIS/IR) and microwave (MW) remote sensing for atmosphere cloud ice and snow. Being able of penetrating cloud to measure cloud ice mass and microphysical properties in the middle to upper troposphere, this is a critical spectrum range for us to understand the connection between cloud ice and precipitation processes.5As the first space-borne 883 GHz radiometer, IceCube mission was NASA’s latest effort in spaceflight demonstration of a commercial sub-mm radiometer technology. Successfully launched from the International Space Station, IceCube is essentially a free-running radiometer and collected valuable 15-month measurements of atmosphere and cloud ice. This paper describes the detailed procedures for Level 1 data calibration, processing and validation. The scientific quality and values of IceCube data are then discussed, including radiative transfer model validation and evaluation, as well as the unique spatial distribution10and diurnal cycle of cloud ice that are revealed for the first time on a quasi-global scale at this frequency.

IceCube↗

The First Global 883 GHz Cloud Ice Survey: IceCube Level 1 Data Calibration, Processing and Analysis

Sub-millimeter (sub-mm, 200 - 1000 GHz) wavelengths contribute a unique capability to fill-in the sensitivity gap between operational visible/infrared (VIS/IR) and microwave (MW) remote sensing for atmosphere cloud ice and snow. Being able of penetrating cloud to measure cloud ice mass and microphysical properties in the middle to upper troposphere, this is a critical spectrum range for us to understand the connection between cloud ice and precipitation processes.5As the first space-borne 883 GHz radiometer, IceCube mission was NASA’s latest effort in spaceflight demonstration of a commercial sub-mm radiometer technology. Successfully launched from the International Space Station, IceCube is essentially a free-running radiometer and collected valuable 15-month measurements of atmosphere and cloud ice. This paper describes the detailed procedures for Level 1 data calibration, processing and validation. The scientific quality and values of IceCube data are then discussed, including radiative transfer model validation and evaluation, as well as the unique spatial distribution10and diurnal cycle of cloud ice that are revealed for the first time on a quasi-global scale at this frequency

Jie Gong↗

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↗

Jet charge with global event shapes: Probing quark flavor dynamics

We propose measuring the jet electric charge of jet regions, defined within the framework of global event shapes, as a probe of quark flavor dynamics within the nucleon and the hadronization process. In particular, we consider a measurement of the jet region charge while simultaneously keeping track of the energy flow throughout the event, as characterized by the global event shape. As a concrete example, we focus on the measurement of the (Q), the jet charge of the jet region (𝐽) defined within the framework of the 1-Jettiness global event shape (𝜏 1 ) for the Deep Inelastic Scattering (DIS) process, 𝑒 − +𝑝→𝑒 − +𝐽+𝑋, with unpolarized or longitudinally polarized protons. The 1-Jettiness distribution, binned according to jet charge, allows for enhanced quark flavor separation of the initial state unpolarized or polarized PDFs. On the other hand, the jet charge distribution binned by 1-Jettiness can serve as a probe of quark flavor dynamics in the final state hadronization process. We derive a factorization theorem for simultaneous measurements of 𝜏 1 and Q in the resummation region, 𝜏 1 ≪𝑃 𝐽 𝑇 , where 𝑃 𝐽 𝑇 denotes the transverse momentum of the jet region. The correlation between the jet region charge and the charge of the struck quark is expected to be the strongest in the resummation region. The factorization theorem contains a new universal , generalizing the standard jet function to include a jet charge measurement. Therefore these universal functions can be extracted from a global analysis of N-jettiness and thrust at 𝑒 + ⁢𝑒 − colliders. We provide simulation studies to demonstrate the sensitivity of the 1-Jettiness jet charge observable to quark flavor dynamics in nucleon structure and explore the possibility of probing the final state hadronization process. This observable is well-suited for applications with existing HERA data and the future Electron-Ion Collider (EIC).

Mantry, Sonny [University of North Georgia, Dahlon↗

Engineering well-expressed, V2-immunofocusing HIV-1 envelope glycoprotein membrane trimers for use in heterologous prime-boost vaccine regimens

HIV-1 vaccine immunofocusing strategies may be able to induce broadly-reactive neutralizing antibodies (NAbs). Here, we engineered a panel of diverse, membrane-resident native HIV-1 trimers vulnerable to two broad targets—the V2 apex and fusion peptide (FP). Selection criteria included i) high expression and ii) infectious function, so that trimer neutralization sensitivity can be profiled in pseudovirus (PV) assays. Initially, we boosted the expression of 17 candidate trimers by truncating gp41 and introducing a gp120-gp41 SOS disulfide to prevent gp120 shedding. "Repairs" were made to fill glycan holes and eliminate other strain-specific aberrations. A new neutralization assay allowed PV infection when our standard assay was insufficient. Trimers with exposed V3 loops, a target of non-NAbs, were discarded. To try to increase V2-sensitivity, we removed clashing glycans and modified the C-strand. Notably, a D167N mutation improved V2-sensitivity in several cases. Glycopeptide analysis of JR-FL trimers revealed near complete sequon occupation and that filling the N197 glycan hole was well-tolerated. In contrast, sequon optimization and inserting/removing glycans at other positions frequently had global "ripple" effects on glycan maturation and sequon occupation throughout the gp120 outer domain and gp41. V2 MAb CH01 selectively bound to trimers with small high mannose glycans near the base of the V1 loop, thereby avoiding clashes. Knocking in a rare N49 glycan was found to perturb gp41 glycans, increasing FP NAb sensitivity—and sometimes improving expression. Finally, a biophysical analysis of VLPs revealed that i) ~25% of particles bear Env spikes, ii) spontaneous particle budding is high and only increases 4-fold upon Gag transfection, and iii) Env+ particles express ~30–40 spikes. Taken together, we identified 7 diverse trimers with a range of sensitivities to two targets to allow rigorous testing of immunofocusing vaccine concepts.

60 APPLIED LIFE SCIENCES↗

Venus global radar reflectivity and correlations with elevation

A geological analysis of the global radar reflectivity data for Venus is presented. The statistical and map pattern distribution of the radar reflectivity rho is first analyzed in the context of the porosity structure of the surface. Rho is then correlated with elevation in order to examine the degree of topographic control on material properties which could depend on PT-sensitive chemical weathering phenomena. Comparison of the rms slope alpha and rho correlations with elevation demonstrates that there are distinct topographic zones defined on the basis of the rho and alpha trends. Finally, hierarchical clustering analysis is applied to the global rho, alpha, and 100-km baseline regional slope data sets in order to define radar property-related subregions.

Garvin, J. B.↗

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Weather Sensitive High Spatio-Temporal Resolution Transportation Electric Load Profiles For Multiple Decarbonization Pathways

Electrification of transport compounded with climate change will transform hourly load profiles and their response to weather. We present a novel approach to generating hourly electric load profiles that considers charging strategies and evolving sensitivity to temperature. The approach consists of downscaling annual state-scale sectoral load projections from the multisectoral Global Change Analysis Model (GCAM) into hourly electric load profiles leveraging high resolution climate and population datasets. Profiles are developed and evaluated at the Balancing Authority scale, with a 5-year increment until 2050 over the Western U.S. Interconnect for multiple decarbonization pathways and climate scenarios. The datasets are readily available for production cost model analysis. Our open source approach is transferable to other regions.

Decarbonization, tranportation, Electric vehicle c↗

Water Ice Clouds in the Martian Atmosphere: A Comparison of Two Methods and Eras

Similar cloud features are seen in maps generated with each method with no obvious outliers. The temperature differencing method appears to possibly be somewhat more sensitive to weaker water ice signatures. We have also generated correlation plots comparing the two methods. At strong delta-T signals, the correlation between the two methods is quite good, and therefore extraction of opacities from earlier Viking data may be possible for these stronger detection levels. Weaker detections do not, however, show such a good correlation. We are currently analyzing why the correlation becomes poor at weak signal levels, though it may be due to the fact that the differencing method may be more sensitive to thin cloud hazes. Results of this ongoing analysis will be presented. A comparison of the Viking and Mars Global Surveyor (MGS) eras are also presented.

A S Hale↗

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

Anthromes and forest carbon responses to global change

Human effects on ecosystems date back thousands of years, and anthropogenic biomes—anthromes—broadly incorporate the effects of human population density and land use on ecosystems. Forests are integral to the global carbon cycle, containing large biomass carbon stocks, yet their responses to land use and climate change are uncertain but critical to informing climate change mitigation strategies, ecosystem management, and Earth system modeling. Using an anthromes perspective and the site locations from the Global Forest Carbon (ForC) Database, we compare intensively used, cultured, and wildland forest lands in tropical and extratropical regions. We summarize recent past (1900-present) patterns of land use intensification, and we use a feedback analysis of Earth system models from the Coupled Model Intercomparison Project Phase 6 to estimate the sensitivity of forest carbon stocks to CO 2 and temperature change for different anthromes among regions. Modeled global forest carbon stock responses are positive for CO 2 increase but neutral to negative for temperature increase. Across anthromes (intensively used, cultured, and wildland forest areas), modeled forest carbon stock responses of temperate and boreal forests are less variable than those of tropical forests. Tropical wildland forest areas appear especially sensitive to CO 2 and temperature change, with the negative temperature response highlighting the potential vulnerability of the globally significant carbon stock in tropical forests. The net effect of anthropogenic activities—including land-use intensification and environmental change and their interactions with natural forest dynamics—will shape future forest carbon stock changes. These interactive effects will likely be strongest in tropical wildlands.

54 ENVIRONMENTAL SCIENCES↗