Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable importance”

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 487 records · Page 27

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗

Evidence-Based Approach to the Analysis of Serious Decompression Sickness with Application to EVA Astronauts

It is important to understand the risk of serious hypobaric decompression sickness (DCS) in order to develop procedures and treatment responses to mitigate the risk. Since it is not ethical to conduct prospective tests about serious DCS with humans, the necessary information was gathered from 73 published reports. We hypothesize that a 4-hr 100% oxygen (O2) prebreathe results in a very low risk of serious DCS, and test this through analysis. We evaluated 258 tests containing information from 79,366 exposures in attitude chambers. Serious DCS was documented in 918 men during the tests. Serious DCS are signs and symptoms broadly classified as Type II DCS. A risk function analysis with maximum likelihood optimization was performed to identify significant explanatory variables, and to create a predictive model for the probability of serious DCS [P(serious DCS)]. Useful variables were Tissue Ratio, the planned time spent at altitude (T(sub alt)), and whether or not repetitive exercise was performed at altitude. Tissue Ratio is P1N2/P2, where P1N2 is calculated nitrogen (N2) pressure in a compartment with a 180-min half-time for N2 pressure just before ascent, and P2 is ambient pressure after ascent. A prebreathe and decompression profile Shuttle astronauts use for extravehicular activity (EVA) includes a 4-hr prebreathe with 100% O2, an ascent to P2 = 4.3 lb per sq. in. absolute, and a T(sub alt) = 6 hr. The P(serious DCS) is: 0.0014 (0.00096 - 0.00196, 95% confidence interval) with exercise and 0.00025 (0.00016 - 0.00035) without exercise. Given 100 Shuttle EVAs to date and no report of serious DCS, the true risk is less than 0.03 with 95% confidence (Binomial Theorem). It is problematic to estimate the risk of serious DCS since it appears infrequently, even if the estimate is based on thousands of altitude chamber exposures. The true risk to astronauts may lie between the extremes of the confidence intervals (0.00016 - 0.00196) since the contribution of other factors, particularly exercise, to the risk of serious DCS during EVA is unknown. A simple model that only accounts for four important variables in retrospective data is still helpful to increase our understanding about the risk of serious DCS.

Conkin, Johnny↗

Mid-Pliocene El Niño/Southern Oscillation Suppressed By Pacific Intertropical Convergence Zone Shift

The El Niño/Southern Oscillation (ENSO), the dominant driver of year-to-year climate variability in the equatorial Pacific Ocean, impacts climate pattern across the globe. However, the response of the ENSO system to past and potential future temperature increases is not fully understood. Here we investigate ENSO variability in the warmer climate of the mid-Pliocene (~3.0–3.3 Ma), when surface temperatures were ~2–3 °C above modern values, in a large ensemble of climate models—the Pliocene Model Intercomparison Project. We show that the ensemble consistently suggests a weakening of ENSO variability, with a mean reduction of 25% (±16%). We further show that shifts in the equatorial Pacific mean state cannot fully explain these changes. Instead, ENSO was suppressed by a series of off-equatorial processes triggered by a northward displacement of the Pacific intertropical convergence zone: weakened convective feedback and intensified Southern Hemisphere circulation, which inhibit various processes that initiate ENSO. The connection between the climatological intertropical convergence zone position and ENSO we find in the past is expected to operate in our warming world with important ramifications for ENSO variability.

Atmospheric dynamics↗

Imprecise global sensitivity analysis using bayesian multimodel inference and importance sampling

Global Sensitivity Analysis (GSA) aims to understand the relative importance of uncertain input variables to model response. Conventional GSA involves calculating sensitivity (Sobol’) indices for a model with known model parameter distributions. However, model parameters are affected by aleatory and epistemic uncertainty, with the latter often caused by lack of data. In this paper, we propose a new framework to quantify uncertainty in probability model-form and model parameters resulting from small datasets and integrate these uncertainties into Sobol’ index estimates. First, the process establishes, through Bayesian multimodel inference, a set of candidate probability models and their associated probabilities. Imprecise Sobol’ indices are calculated from these probability models using an importance sampling reweighting approach. This results in probabilistic Sobol’ indices, whose distribution characterizes uncertainty in the sensitivity resulting from small dataset size. The imprecise Sobol’ indices thus provide a measure of confidence in the sensitivity estimate and, moreover, can be used to inform data collection efforts targeted to minimize the impact of uncertainties. Through an example studying the parameters of a Timoshenko beam, we show that these probabilistic Sobol’ indices converge to the true/deterministic Sobol’ indices as the dataset size increases and hence, distribution-form uncertainty reduces. The approach is then applied to assess the sensitivity of the out-of-plane properties of an E-glass fiber composite material to its constituent properties. This second example illustrates the approach for an important class of materials with wide-ranging applications when data may be lacking for some input parameters.

42 ENGINEERING↗

Atmospheric effects on earth rotation and polar motion

The variability in the earth's rotation rate not due to known solid body tides is dominated on time scales of about four years and less by variations in global atmospheric angular momentum (M) as derived from the zonal wind distribution. Among features seen in the length of day record produced by atmospheric forcing are the strong seasonal cycle, quasi-periodic fluctuations around 40-50 days, and an interannual signal forced by a strong Pacific warming event known as the El Nino. Momentum variations associated with these time scales arise in different latitudinal regions. Furthermore, winds in the stratosphere make a particularly important contribution to seasonal variability. Other related topics discussed here are: (1) comparisons of the M series from wind fields produced at different weather centers; (2) the torques that dynamically link the atmosphere and earth; and (3) longer-term nonatmospheric effects that can be seen upon removal of the atmospheric signal.an interestigapplication for climatological purposes is the use of the historical earth rotation series as a proxy for atmospheric wind variability prior to the era of upper-air data. Lastly, results pertaining to the role of atmospheric pressure systems in exciting rapid polar motion are presented.

Salstein, David A.↗

Quantifying spatiotemporal variability in occupant exposure to an indoor airborne contaminant with an uncertain source location

Well-mixed zone models are often employed to compute indoor air quality and occupant exposures. While effective, a potential downside to assuming instantaneous, perfect mixing is underpredicting exposures to high intermittent concentrations within a room. When such cases are of concern, more spatially resolved models, like computational-fluid dynamics methods, are used for some or all of the zones. But, these models have higher computational costs and require more input information. A preferred compromise would be to continue with a multi-zone modeling approach for all rooms, but with a better assessment of the spatial variability within a room. Here to do so, we present a quantitative method for estimating a room's spatiotemporal variability, based on influential room parameters. Our proposed method disaggregates variability into the variability in a room's average concentration, and the spatial variability within the room relative to that average. This enables a detailed assessment of how variability in particular room parameters impacts the uncertain occupant exposures. To demonstrate the utility of this method, we simulate contaminant dispersion for a variety of possible source locations. We compute breathing-zone exposure during the releasing (source is active) and decaying (source is removed) periods. Using CFD methods, we found after a 30 minutes release the average standard deviation in the spatial distribution of exposure was approximately 28% of the source average exposure, whereas variability in the different average exposures was lower, only 10% of the total average. We also find that although uncertainty in the source location leads to variability in the average magnitude of transient exposure, it does not have a particularly large influence on the spatial distribution during the decaying period, or on the average contaminant removal rate. By systematically characterizing a room's average concentration, its variability, and the spatial variability within the room important insights can be gained as to how much uncertainty is introduced into occupant exposure predictions by assuming a uniform in-room contaminant concentration. We discuss how the results of these characterizations can improve our understanding of the uncertainty in occupant exposures relative to well-mixed models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The search for, and study of black holes in the galaxy

The study of the short-term (less than 10 s) time variability of Cygnus X-1 is reviewed, with attention given to measurements of the autocorrelation function and of the cross-correlation between different energy bandwidths, as well as to a simple shot noise model for describing all such short-term variations. Observational questions with regard to Cygnus X-1 are discussed, relating to the average shape of the randomly occurring flares, energy dependence in the cross-correlation and the dependence on energy of the X-ray polarization. Concerning the study of the Cyg X-1-like aperiodic X-ray sources, an estimate is made of what can be accomplished at the limits of sensitivity where Poisson fluctuations are at least as important as the intrinsic variability of the source, considering the third moment due to the source fluctuations.

Weisskopf, M. C.↗

Spatial and temporal variability of the extent of chemically processed stratospheric air

Simulations of the spatial and temporal variability of the extent of chemically processed air in the Arctic stratosphere have been carried out using a three-dimensional chemistry-transport model for the winters of 1979 and 1989. Chemically processed air is identified in the model as that in which the amounts of hydrogen chloride (HCl) calculated with parameterized loss for conditions appropriate to polar stratospheric cloud (PSC) formation are substantially smaller than those calculated in a model with gas phase chemistry only. It is seen that chemically processed air may be identified over much of the Arctic lower stratosphere from early January to late February, with HCl depletions being larger in 1989 than in 1979. Near the latitude of the Arctic circle, there is important spatial and temporal variability in the extent of chemically processed air. There is some evidence for transport to midlatitudes of processed air during these winters, but the HCl reductions are much smaller and more sporadic than those near the pole. At 62 and 42N, processed air is calculated to occur preferentially over the longitude regions from 60-120E and 270-330E.

Kaye, Jack A.↗

Impact of a dynamic grid mix and climate on operational carbon emissions modeling for different building typologies and climate zones

Calculating operational carbon emissions through a building’s lifecycle is complex due to the dynamic nature of influencing factors such as climate and energy grid mix. This paper introduces a novel methodology for modeling 30-year operational carbon impacts of buildings and applies this method to mid-rise office and residential typologies across various US climate zones. The method accounts for these temporal variabilities using new and scarcely cited data sources. Key findings indicate that future changes in the climate, while impactful, play a relatively modest role in operational carbon emissions compared to significant reductions with modeling scenarios using the projected decarbonization of the electricity grid. Here, the study also finds that using annual, month-hourly, or hourly grid emission factors have a minimal impact on carbon accounting, except in certain climates and program types where emission patterns do not align with a building’s energy consumption. Warmer climates like Miami, Florida and Tucson, Arizona, which rely heavily on cooling, demonstrate larger variations in carbon emissions when using higher temporal resolution emission factors. Ultimately, this study underscores the critical role of grid decarbonization in reducing long-term emissions and the importance of incorporating this variable in life cycle assessment (LCA) modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reconstruction of Snow on Arctic Sea Ice

Snow on Arctic sea ice is a poorly observed variable that plays an important role in the Arctic climate system and impacts the remote sensing systems that monitor Arctic sea ice. We present and validate a reconstruction of Arctic snow depth based on observed sea ice motion and snowfall derived from reanalysis data. Overall, the reconstruction is in good agreement with direct measurements of snow depth from Operation IceBridge, slightly underestimating mean IceBridge snow depth. At the local scale (10 km), the reconstruction is more skilled than a weighted climatology over first year ice, but underestimates deeper snow over multiyear ice. Reconstructions of single buoy snow depths are unskilled, but the reconstruction shows skill in simulating the mean snow depth across all buoys. Spring snow depths show a lowtohigh crossArctic gradient and tend to be greatest in the Atlantic sector of the eastern Arctic. The relationship between ice type (multiyear or firstyear ice) and snow depth previously documented in the western Arctic is not evident in the eastern Arctic. Using ice type to weight snow depths for satellite ice thickness retrievals may not be justifiable in the eastern Arctic. Reconstructed snow depth across the Arctic shows significant interannual variability, suggesting that use of a fixed snow depth climatology may lead to biases in retrieved ice thickness and its variability. However, interannual variability in panArctic mean snow depth is comparable or smaller than the uncertainty in both the reconstruction and IceBridge, highlighting the need for high accuracy snow depth products and reconstructions.

sea ice↗

Subdecadal Holocene Warm-Season Temperature Variability in Central Europe Recorded by Biochemical Varves

Paleoclimate data provide important information about the character of natural climate variability. However, records with sufficient length and resolution to resolve high-frequency (decadal-scale) variability across the Holocene are scarce. We present a 10,800-year reconstruction of spring and summer temperature at three-year resolution based on biochemical varves from Lake $\dot{Z}abi$$\acute{n}skie$, Poland. The reconstruction is based on Ca/Ti ratio, which are significantly correlated with instrumental spring and summer temperature spanning 240 years. Major climate events of the Holocene period are represented in the reconstruction, including the Holocene Thermal Maximum, 8.2 ka Event, Medieval Climate Anomaly, and Little Ice Age. A low-frequency 8,000-year decreasing trend in warm-season temperatures is driven by declining summer insolation. Temperature variability is highest during the early Holocene, likely related to warmer and drier conditions. The rate of warming during the past 90 years is extremely unusual, if not unprecedented for the Holocene, based on our reconstruction.

58 GEOSCIENCES↗

The Freundlich isotherm equation best represents phosphate sorption across soil orders and land use types in tropical soils of Puerto Rico

Biomass production in the lowland wet tropical forest is greater than in any other biome, and it is typically limited by soil phosphorus (P) availability. However, the mechanisms involved in the P cycle remain poorly represented in Earth System Models (ESMs). Soil P sorption processes are key in the P cycle and for understanding the extent of P limitation for plant productivity. Currently, a few ESMs include isotherm equations to model these processes. Although the Langmuir equation is widely cited, other isotherm equations may better describe sorption in tropical soils. Here, we use a diverse range of soil samples from Puerto Rico to test the validity of the Langmuir, Freundlich, and Temkin equation. We found that across four soil orders (Inceptisols, Mollisols, Oxisols, Ultisols), and forested and cultivated land use types, the Freundlich equation best represented soil P sorption. Furthermore, the Langmuir and the Temkin equations poorly represent soil P adsorption, especially at low P concentrations. Specifically, the Langmuir equation underestimated soil P adsorption by 40% and the Temkin equation overestimated adsorption by 76%. We also found, as expected, that soil clay content and pH were the most important parameters explaining the variability of the Freundlich (K f ) constant. Greater clay content and lower pH, common in highly weathered Ultisols and Oxisols which are abundant in the tropics, led to greater K f values. Overall, our results suggest that a diversity of soils can prompt underestimation of P sorption when using the Langmuir isotherm, which leads to an overestimation of available P that can have repercussions on ESM predictions of the P cycle and tropical forest productivity.

Earth system models↗

Concrete shrinkage and creep under drying/wetting cycles

Concrete shrinkage and creep under variable hydric conditions are important factors for the safety and durability of concrete especially in nuclear reactor or nuclear waste storage background. A large (and long – more than 900 days) experimental campaign, conducted on two different concretes, has therefore been designed to study the strains of non loaded and loaded concretes submitted to drying and liquid water imbibition cycles. For the purposes of comparison, concrete strains and mass variations during drying only (50% RH) and/or following cycles of drying/rewetting were also recorded. This allowed the identification of desiccation shrinkage, basic creep and drying creep at 10 MPa of axial stress. Important results were found and have shown that the final mass and strain are not deeply modified by the introduction of a rewetting phase either for a loaded or a non loaded material.

36 MATERIALS SCIENCE↗

Dynamic energy system modeling using hybrid physics-based and machine learning encoder–decoder models

Three model configurations are presented for multi-step time series predictions of the heat absorbed by the water and steam in a thermal power plant. The models predict over horizons of 2, 4, and 6 steps into the future, where each step is a 5-minute increment. The evaluated models are a pure machine learning model, a novel hybrid machine learning and physics-based model, and the hybrid model with an incomplete dataset. The hybrid model deconstructs the machine learning into individual boiler heat absorption units: economizer, water wall, superheater, and reheater. Each configuration uses a gated recurrent unit (GRU) or a GRU-based encoder–decoder as the deep learning architecture. Mean squared error is used to evaluate the models compared to target values. The encoder–decoder architecture is over 11% more accurate than the GRU only models. The hybrid model with the incomplete dataset highlights the importance of the manipulated variables to the system. The hybrid model, compared to the pure machine learning model, is over 10% more accurate on average over 20 iterations of each model. Automatic differentiation is applied to the hybrid model to perform a local sensitivity analysis to identify the most impactful of the 72 manipulated variables on the heat absorbed in the boiler. The models and sensitivity analyses are used in a discussion about optimizing the thermal power plant.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗

Identifying schools at high-risk for elevated lead in drinking water using only publicly available data

Estimating the risk of lead contamination of schools' drinking water at the State level is a complex, important, and unexplored challenge. Variable water quality among water systems and changes in water chemistry during distribution affect lead dissolution rates from pipes and fittings. In addition, the locations of lead-bearing plumbing materials are uncertain. We tested the capability of six machine learning models to predict the likelihood of lead contamination of drinking water at the schools' taps using only publicly available datasets. The predictive features used in the models correspond to those with a proven correlation to the dominant, but commonly unavailable, factors that govern lead leaching: the presence of lead-bearing plumbing materials and water quality conducive to lead corrosion. By combining water chemistry data from public reports, socioeconomic information from the US census, and spatial features using Geographic Information Systems, we trained and tested models to estimate the likelihood of lead contaminated tap water in over 8,000 schools across California and Massachusetts. Our best-performing model was a Random Forest, with a 10-fold cross validation score of 0.88 for Massachusetts and 0.78 for California using the average Area Under the Receiver Operating Characteristic Curve (ROC AUC) metric. The model was then used to assign a lead leaching risk category to half of the schools across California (the other half was used for training). There was good agreement between the modeled risk categories and the actual lead leaching outcomes for every school; however, the model overestimated the lead leaching risk in up to 17% of the schools. This model is the first of its kind to offer a tool to predict the risk of lead leaching in schools at the State level. Further use of this model can help deploy limited resources more effectively to prevent childhood lead exposure from school drinking water.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗