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,207 records · Page 67

Arctic Sea Ice in CMIP6

Abstract We examine CMIP6 simulations of Arctic sea‐ice area and volume. We find that CMIP6 models produce a wide spread of mean Arctic sea‐ice area, capturing the observational estimate within the multimodel ensemble spread. The CMIP6 multimodel ensemble mean provides a more realistic estimate of the sensitivity of September Arctic sea‐ice area to a given amount of anthropogenic CO 2 emissions and to a given amount of global warming, compared with earlier CMIP experiments. Still, most CMIP6 models fail to simulate at the same time a plausible evolution of sea‐ice area and of global mean surface temperature. In the vast majority of the available CMIP6 simulations, the Arctic Ocean becomes practically sea‐ice free (sea‐ice area <1 × 10 6 km 2 ) in September for the first time before the Year 2050 in each of the four emission scenarios SSP1‐1.9, SSP1‐2.6, SSP2‐4.5, and SSP5‐8.5 examined here.

Notz, Dirk↗

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗

Evaluating Climate Models with the CLIVAR 2020 ENSO Metrics Package

El Niño–Southern Oscillation (ENSO) is the dominant mode of interannual climate variability on the planet, with far-reaching global impacts. It is therefore key to evaluate ENSO simulations in state-of-the-art numerical models used to study past, present, and future climate. Recently, the Pacific Region Panel of the International Climate and Ocean: Variability, Predictability and Change (CLIVAR) Project, as a part of the World Climate Research Programme (WCRP), led a community-wide effort to evaluate the simulation of ENSO variability, teleconnections, and processes in climate models. The new CLIVAR 2020 ENSO metrics package enables model diagnosis, comparison, and evaluation to 1) highlight aspects that need improvement; 2) monitor progress across model generations; 3) help in selecting models that are well suited for particular analyses; 4) reveal links between various model biases, illuminating the impacts of those biases on ENSO and its sensitivity to climate change; and to 5) advance ENSO literacy. By interfacing with existing model evaluation tools, the ENSO metrics package enables rapid analysis of multipetabyte databases of simulations, such as those generated by the Coupled Model Intercomparison Project phases 5 (CMIP5) and 6 (CMIP6). The CMIP6 models are found to significantly outperform those from CMIP5 for 8 out of 24 ENSO-relevant metrics, with most CMIP6 models showing improved tropical Pacific seasonality and ENSO teleconnections. Only one ENSO metric is significantly degraded in CMIP6, namely, the coupling between the ocean surface and subsurface temperature anomalies, while the majority of metrics remain unchanged.

54 ENVIRONMENTAL SCIENCES↗

Microphysical Sensitivity of Superparameterized Precipitation Extremes in the Contiguous United States Due to Feedbacks on Large-Scale Circulation

Superparameterized (SP) global climate models have been shown to better simulate various features of precipitation relative to conventional models, including its diurnal cycle as well as its extremes. While various studies have focused on the effect of differing microphysics parameterizations on precipitation within limited-area cloud-resolving models, we examine here the effect on contiguous U.S. (CONUS) extremes in a global SP model. We vary the number of predicted moments for hydrometeor distributions, the character of the rimed ice species, and the representation of raindrop self-collection and breakup. Using a likelihood ratio test and accounting for the effects of multiple hypothesis testing, we find that there are some regional differences, particularly during spring and summer in the Southwest and the Midwest, in both the current climate and a warmer climate with uniformly increased sea surface temperatures. These differences are most statistically significant and widespread when the number of moments is changed. To determine whether these results are due to (fast) local effects of the different microphysics or the (slower) ensuing feedback on the large-scale atmospheric circulation, we run a series of short, 5-day simulations initialized from reanalysis data. We find that the differences largely disappear in these runs and therefore infer that the different parameterizations impact precipitation extremes indirectly via the large-scale circulation. Finally, we compare the present-day results with hourly rain gauge data and find that SP underestimates extremes relative to observations regardless of which microphysics scheme is used given a fixed model configuration and resolution.

54 ENVIRONMENTAL SCIENCES↗

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State

The COVID-19 pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the “D-form” to the “G-form”) that carried an amino acid substitution D614G in its “Spike” protein. The G-form is more infectious in vitro and associated with increased viral loads. To gain insight into the molecular-level underpinnings of these characteristics, we employed microsecond all-atom simulations. Here we show that changes in the protein energetics favor a higher population of infection-capable (open) states through release of hydrogen bonds of an asymmetry present in the D-form but not the G-form. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are significant for vaccine design. The Molecular Dynamics Simulations datasets generated in this study, namely that of the soluble form SARS-CoV-2 Spike protein are made available here. Details of the trajectory and supporting files are given in the README.txt file.

59 BASIC BIOLOGICAL SCIENCES↗

Spatial Decorrelation of Young Stars and Dense Gas as a Probe of the Star Formation–Feedback Cycle in Galaxies

The spatial decorrelation of dense molecular gas and young stars observed on ≲1 kpc scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. In this work, we explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC 300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC 300, and we show that this agreement is due to different aspects of feedback, including H 2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affects the correlation on ≲100 pc scales, while supernovae play a significant role on ≳ 100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per free fall time, ϵ ff , which provides a strong observational constraint on ϵ ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗

Nucleosynthesis of an 11.8 M ⊙ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-lived Radionuclides in the Early Solar System

Based on a 3D supernova simulation of an 11.8 M ⊙ progenitor model with an initial solar composition, we study the nucleosynthesis using tracers covering the innermost 0.1 M ⊙ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of 45 Sc and 64 Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an 11.8 M ⊙ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on 41 Ca, 53 Mn, and 60 Fe without exceeding those on the other SLRs.

Astronomy & Astrophysics↗

Modeling Interactive Effects of Manganese Bioavailability, Nitrogen Deposition, and Warming on Soil Carbon Storage

Manganese (Mn) is a redox–active micronutrient that has been shown to accelerate plant litter decomposition; however, the effect of Mn–promoted decomposition on soil C storage is unclear. We present a novel biogeochemical model simulating how Mn bioavailability influences soil organic C (SOC) stocks in a soil profile (<50 cm) within a temperate forest. In our model, foliar Mn increased in response to increasing soluble Mn released through Mn–oxide (birnessite) dissolution in mineral soil layers. The ensuing Mn enrichment in leaf litter redistributed Mn to the surface forest floor layer, promoted enzymatic oxidation of lignin, and decreased SOC stocks. Total SOC loss was partially mitigated by accumulation of lignin–oxidation products as mineral–associated organic C. We also explored how Mn–driven changes to C storage interacted with effects of N deposition and warming. Nitrogen enrichment inhibited Mn–dependent lignin degradation, increasing SOC stocks and weakening their dependence on Mn bioavailability. Warming stimulated decomposition and reduced C stocks but was less effective at low Mn bioavailability. Our model results suggest that SOC stocks are sensitive to Mn bioavailability because increased plant uptake redistributes Mn to surface soils where it can enhance litter decomposition. Finally, based on our simulations, we predict that Mn becomes limiting to litter decomposition where Mn is poorly soluble. Depletion of bioavailable Mn or other cofactors that are critical to decomposition could limit the response of organic C stocks to warming over time, but quantitative projections of the role of Mn bioavailability in regulating decomposition requires additional measurements to constrain model uncertainties.

58 GEOSCIENCES↗

Using CESM-RESFire to understand climate–fire–ecosystem interactions and the implications for decadal climate variability

Large wildfires exert strong disturbance on regional and global climate systems and ecosystems by perturbing radiative forcing as well as the carbon and water balance between the atmosphere and land surface, while short- and long-term variations in fire weather, terrestrial ecosystems, and human activity modulate fire intensity and reshape fire regimes. The complex climate–fire–ecosystem interactions were not fully integrated in previous climate model studies, and the resulting effects on the projections of future climate change are not well understood. Here we use the fully interactive REgion-Specific ecosystem feedback Fire model (RESFire) that was developed in the Community Earth System Model (CESM) to investigate these interactions and their impacts on climate systems and fire activity. We designed two sets of decadal simulations using CESM-RESFire for present-day (2001–2010) and future (2051–2060) scenarios, respectively, and conducted a series of sensitivity experiments to assess the effects of individual feedback pathways among climate, fire, and ecosystems. Our implementation of RESFire, which includes online land–atmosphere coupling of fire emissions and fire-induced land cover change (LCC), reproduces the observed aerosol optical depth (AOD) from space-based Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products and ground-based AErosol RObotic NETwork (AERONET) data; it agrees well with carbon budget benchmarks from previous studies. We estimate the global averaged net radiative effect of both fire aerosols and fire-induced LCC at –0.59±0.52 W m –2 , which is dominated by fire aerosol–cloud interactions (–0.82±0.19 W m –2 ), in the present-day scenario under climatological conditions of the 2000s. The fire-related net cooling effect increases by ~170 % to –1.60±0.27 W m –2 in the 2050s under the conditions of the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Such considerably enhanced radiative effect is attributed to the largely increased global burned area (+19 %) and fire carbon emissions (+100 %) from the 2000s to the 2050s driven by climate change. The net ecosystem exchange (NEE) of carbon between the land and atmosphere components in the simulations increases by 33 % accordingly, implying that biomass burning is an increasing carbon source at short-term timescales in the future. High-latitude regions with prevalent peatlands would be more vulnerable to increased fire threats due to climate change, and the increase in fire aerosols could counter the projected decrease in anthropogenic aerosols due to air pollution control policies in many regions. We also evaluate two distinct feedback mechanisms that are associated with fire aerosols and fire-induced LCC, respectively. On a global scale, the first mechanism imposes positive feedbacks to fire activity through enhanced droughts with suppressed precipitation by fire aerosol–cloud interactions, while the second one manifests as negative feedbacks due to reduced fuel loads by fire consumption and post-fire tree mortality and recovery processes. These two feedback pathways with opposite effects compete at regional to global scales and increase the complexity of climate–fire–ecosystem interactions and their climatic impacts.

54 ENVIRONMENTAL SCIENCES↗

Quantitative modeling of gate-tunable cavity exciton-polaritons in transition-metal dichalcogenides

We performed a comprehensive real-space simulation study of propagative exciton-polaritons (EPs) of tungsten disulfide (WS 2 ) inside a Fabry–Pérot cavity. Through our simulations, we were able to construct the dispersion relation and determine the propagation length (L p ) of EPs, both of which show sensitive dependence on the applied gate voltage. We found that L p increases at higher gate voltages when excited at the exciton energy but shows opposite gate dependence at the trion energy. These observations are direct evidence of the coupling between the cavity mode with excitons and trions of WS 2 . Furthermore, we demonstrate a proof-of-concept device based on the tunable cavity EPs of WS 2 , which is promising for application as photonic transistors, modulators, or electrical-to-optical converters. Furthermore, our work unveils the real-space transport properties of gate-tunable cavity EPs and paves the way for future applications of EPs in nanophotonic devices and circuits.

42 ENGINEERING↗

Application of Neutron Multiplicity Counting Experiments to Optimal Cross-Section Adjustments

This paper presents the first application of model calibration to neutron multiplicity counting (NMC) experiments for cross-section optimization that is informed by adjoint-based sensitivity analysis (SA) and first-order uncertainty quantification (UQ). We summarize previous work on SA applied to NMC and describe notable modifications and additions. We give the procedure for first-order UQ and Bayesian-inference-based parameter estimation (PE). We then discuss model calibration applied to NMC of a 4.5-kg sphere of weapons-grade, alpha-phase plutonium metal (the BeRP ball) with the nPod neutron multiplicity counter. For the BeRP ball in bare and polyethylene-reflected configurations, we discuss the sensitivity of the first- and second-moment detector responses (i.e., first and second moments of the NMC distribution, respectively) to the cross sections. We describe the sources of uncertainty in the measured and simulated responses. Specifically, uncertainty in the measured responses is due to both random and systematic sources. Uncertainty in the simulated responses is due to the cross-section covariances. We describe in detail the adjustment to the cross sections and cross-section covariances due to the optimization. Due to the contribution of systematic uncertainties to the measured response uncertainties, the adjustment to the cross sections is similar in trend but larger in magnitude compared to that recommended by previous work. We compare the measured responses to responses simulated with nominal and optimized cross sections, demonstrating that the best-estimate cross sections produce simulations of NMC experiments that are more accurate with reduced uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗

Feature extraction for subtle anomaly detection using semi-supervised learning

The demand for automated and effective monitoring techniques has soared with the increased digitization of industrial monitoring systems. State-of-the-art machine learning methods are effectively detecting abrupt changes in system states. However, these methods lack comparable maturity in detecting subtle changes that may be signs of incipient faults. This manuscript argues that the current anomaly detection methods can be enhanced by exploring weak patterns to enable subtle variation detection. Specifically, the concept of semi-supervised learning is employed, with labels representing knowledge about some anomalous conditions of a system. The basic idea is to extract a candidate set of weak patterns discarded by state-of-the-art baselining algorithms. With few labeled anomalous data, the algorithm selects the weak patterns and allows for their possible fusion using the highest sensitivity to the labeled anomalies. Here, the method’s applicability is demonstrated using a representative pressurized water reactor (PWR) model simulated by Dymola.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Squeezing Out Interfacial Solvation: The Role of Hydrogen-Bonding in the Structural and Orientational Freedom of Molecular Self-Assembly

We report bioinspired membrane molecules with improved physical properties and enhanced stability can serve as functional models for conventional lipid or amphiphilic species. Importantly, these molecules can also provide new insights into emergent phenomena that manifest during self-assembly at interfaces. Here, we elucidate the structural response and mechanistic steps underlying the self-assembly of the amphiphilic, charged oligodimethylsiloxane imidazolium cation (ODMS-MIM+) at the air–aqueous interface using Langmuir trough methods with coincident surface-specific vibrational sum-frequency generation (SFG) spectroscopy. We find evidence for a new compression-induced desolvation step that precedes commonly known disordered-to-ordered phase transitions to form nanoscopic assemblies. The experimental data was supported by atomistic molecular dynamics (MD) simulations to provide a detailed mechanistic picture underlying the assembly and the role of water in these phase transitions. The sensitivity of the hydrophobic ODMS tail conformations to compression-owing to distinct water–ODMS interactions and tail–tail solvation properties-offers new strategies for the design of interfaces that can be further used to develop soft-matter electronics and low-dimensional materials using physical and chemical controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating Climate Models’ Cloud Feedbacks Against Expert Judgment

Abstract The persistent and growing spread in effective climate sensitivity (ECS) across global climate models necessitates rigorous evaluation of their cloud feedbacks. Here we evaluate several cloud feedback components simulated in 19 climate models against benchmark values determined via an expert synthesis of observational, theoretical, and high‐resolution modeling studies. We find that models with smallest feedback errors relative to these benchmark values generally have moderate total cloud feedbacks (0.4–0.6 W m −2 K −1 ) and ECS (3–4 K). Those with largest errors generally have total cloud feedback and ECS values that are too large or too small. Models tend to achieve large positive total cloud feedbacks by having several cloud feedback components that are systematically biased high rather than by having a single anomalously large component, and vice versa. In general, better simulation of mean‐state cloud properties leads to stronger but not necessarily better cloud feedbacks. The Python code base provided herein could be applied to developmental versions of models to assess cloud feedbacks and cloud errors and place them in the context of other models and of expert judgment in real‐time during model development.

58 GEOSCIENCES↗

Stratospheric aerosol injection can weaken the carbon dioxide greenhouse effect

Abstract Stratospheric aerosol injection is a proposed method for offsetting greenhouse gas-induced warming by introducing scattering aerosols into the lower stratosphere to reflect sunlight. Here we explore a potentially more efficient alternative: weakening the Earth’s greenhouse effect by deploying absorptive aerosols in the upper stratosphere (~10 hPa). These aerosols warm the carbon dioxide emission level—where outgoing longwave radiation is most sensitive to temperature—thereby enhancing top-of-atmosphere infrared emission without altering atmospheric carbon dioxide concentrations. Idealized climate model simulations indicate that this approach can reduce global temperatures an order of magnitude more efficiently per unit aerosol mass than conventional scattering-based interventions. Although based on simplified model experiments lacking interactive aerosol processes and operational constraints, our results identify a distinct physical mechanism for climate intervention, arguing for further research into the impacts—especially potential unintended side effects—of injecting absorptive aerosols into the upper stratosphere as an alternative solar radiation management strategy.

Environmental Sciences & Ecology↗