Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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 199 records · Page 11

Regional and Global Climate Response to Anthropogenic SO2 Emissions from China in Three Climate Models

We use the HadGEM3-GA4, CESM1, and GISS ModelE2 climate models to investigate the global and regional aerosol burden, radiative flux, and surface temperature responses to removing anthropogenic sulfur dioxide (SO2) emissions from China. We find that the models differ by up to a factor of 6 in the simulated change in aerosol optical depth (AOD) and shortwave radiative flux over China that results from reduced sulfate aerosol, leading to a large range of magnitudes in the regional and global temperature responses. Two of the three models simulate a near-ubiquitous hemispheric warming due to the regional SO2 removal, with similarities in the local and remote pattern of response, but overall with a substantially different magnitude. The third model simulates almost no significant temperature response. We attribute the discrepancies in the response to a combination of substantial differences in the chemical conversion of SO2 to sulfate, translation of sulfate mass into AOD, cloud radiative interactions, and differences in the radiative forcing efficiency of sulfate aerosol in the models. The model with the strongest response (HadGEM3-GA4) compares best with observations of AOD regionally, however the other two models compare similarly (albeit poorly) and still disagree substantially in their simulated climate response, indicating that total AOD observations are far from sufficient to determine which model response is more plausible. Our results highlight that there remains a large uncertainty in the representation of both aerosol chemistry as well as direct and indirect aerosol radiative effects in current climate models, and reinforces that caution must be applied when interpreting the results of modelling studies of aerosol influences on climate. Model studies that implicate aerosols in climate responses should ideally explore a range of radiative forcing strengths representative of this uncertainty, in addition to thoroughly evaluating the models used against observations.

sulfates↗

Climatic consequences of observed ozone loss in the 1980s: Relevance to the greenhouse problem

Recently published findings using satellite and ground-based observations indicate a large winter and summertime decrease in the column abundance of ozone at high and middle latitudes during the last decade. Using a simple ozone depletion profile reflecting the observed decrease in ozone column abundance, Ramaswamy et al. (1992) showed that the negative radiative forcing that results from the ozone decrease between 1979 and 1990 approximately balanced the greenhouse climate forcing due to the chlorofluorocarbons emitted during the same period. Here, we extend the forcing analyses by calculating the equilibrium surface temperature response explicitly, using an updated version of the Atmospheric and Environmental Research two-dimensional radiative-dynamical seasonal model. The calculated steady state responses suggest that the surface cooling due to the ozone depletion in the lower stratosphere offsets about 30% of the surface warming due to greenhouse gases emitted during the same decade. The temperature offset is roughly a factor of 2 larger than the corresponding offset obtained from forcing intercomparisons. This result appears to be related to the climate feedback mechanisms operating in the model troposphere, most notably that associated with atmospheric meridional heat transport. Thus a comprehensive assessment of ozone change effects on the predicted greenhouse warming cannot be accomplished based on forcing evaluations alone. Our results also show that calculations adopting a seasonally and latitudinally dependent ozone depletion profile produce a negative forcing about 50% smaller than that calculated for the depletion profile used by Ramaswamy et al. (1992).

Molnar, G. I.↗

Analysis of Atmospheric Conditions Responsible for an Ozone Exceedance Event in Southeast Virginia on June 15, 2022

On June 15, 2022, the Virginia Department of Environmental Quality monitoring site at Suffolk/Holland, a rural site in southeast Virginia, recorded its first ozone exceedance since 2016. An ozone exceedance day occurs when the daily maximum 8-hour average surface ozone concentration is greater than 70 ppb. On this day, it was 75 ppb. This event is also noteworthy due to the rapid increase in surface ozone between 6 am and 9 am EDT as well as the hourly maximum ozone concentration of 82 ppb that was measured at 1 pm EDT. In this analysis, we utilize various observational and model data to determine the atmospheric conditions responsible for this ozone exceedance event. The analysis is conducted in two parts: (1) an evaluation of the accuracy of the GEOS-CF and WRF-Chem model forecasts and (2) an investigation of the mechanisms responsible for the high surface ozone at the Virginia DEQ Suffolk/Holland monitoring site. Comparisons of model forecasted ozone time-height cross sections with measured ozone by lidars from the Tropospheric Ozone Lidar Network (TOLNet) at NASA Langley and NASA Goddard reveal that both models forecasted ozone reasonably well near the surface. Model and observational data reveal that an upper-level ridge was present over much of the eastern United States, with a broad anticyclonic circulation near the surface, resulting in north-northeasterly flow over southeast Virginia. Back trajectory calculations using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) and model data suggest that a low-level plume rich in ozone and precursor species originating from the New York City and New Jersey region was transported along the east coast between June 14 and June 15, arriving in southeast Virginia around 6 am EDT on June 15. In the early morning hours, this plume mixed down to the surface, elevating the background ozone mixing ratio as well as the mixing ratios of several precursor species. Other potential contributions are also explored and discussed. Lastly, the authors note that this work represents the analysis of the authors and not the Virginia Department of Environmental Quality.

Daniel B Phoenix↗

Effects of Forcing Differences and Initial Conditions on Inter-Model Agreement in the VolMIP Volc-Pinatubo-Full Experiment

This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.

metadata↗

Cloud Feedback Uncertainty in the Equatorial Pacific Across CMIP6 Models

Cloud feedback is the largest uncertainty in estimating Equilibrium Climate Sensitivity. In this study we focus on the equatorial Pacific, where CMIP6 model cloud feedback spread is notably large. Cloud radiative effects in this region are relevant for the global climate. Our findings show that models predict a consistent shift towards the ascent regime in response to El Nino-like sea surface warming. Models diverge in terms of the radiative impact due to differences in cloud characteristics in ascent and subsidence regimes. Using the observed relationship between circulation regime and cloud radiative effect, we find a reduction in the regional mean cloud feedback estimate from 0.77 to 0.22 W m -2 K -1 , though this does not substantially lessen the model spread in total feedback. Pathways to reduce this spread include: improving confidence in the regional ocean warming pattern, and using observations and models to understand cloud type and circulation interactions.

CMIP6↗

The HHV-6B U20 glycoprotein binds ULBP1, masking it from recognition by NKG2D and interfering with natural killer cell activation

Human Herpesvirus 6B (HHV-6B) impedes host immune responses by downregulating class I MHC molecules (MHC-I), hindering antigen presentation to CD8+ T cells. Downregulation of MHC-I disengages inhibitory receptors on natural killer (NK) cells, resulting in activation and killing of the target cell if NK cell activating receptors such as NKG2D have engaged stress ligands upregulated on the target cells. Previous work has shown that HHV-6B downregulates three MHC-like stress ligands MICB, ULBP1, and ULBP3, which are recognized by NKG2D. The U20 glycoprotein of the related virus HHV-6A has been implicated in the downregulation of ULBP1, but the precise mechanism remains undetermined. We set out to investigate the role of HHV-6B U20 in modulating NK cell activity. We used HHV-6B U20 expressed as a recombinant protein or transduced into target cells, as well as HHV-6B infection, to investigate binding interactions with NK cell ligands and receptors and to assess effects on NK cell activation. Small-angle X-ray scattering was used to align molecular models derived from machine-learning approaches. We demonstrate that U20 binds directly to ULBP1 with sub-micromolar affinity. Transduction of U20 decreases NKG2D binding to ULBP1 at the cell surface but does not decrease ULBP1 protein levels, either at the cell surface or in toto. HHV-6B infection and soluble U20 have the same effect. Transduction of U20 blocks NK cell activation in response to cell-surface ULBP1. Structural modeling of the U20 – ULBP1 complex indicates some similarities to the m152-RAE1γ complex.

60 APPLIED LIFE SCIENCES↗

A Statistics-Based Material Property Analysis to Support TPS Characterization

Accurate characterization of entry capsule heat shield material properties is a critical component in modeling and simulating Thermal Protection System (TPS) response in a prescribed aerothermal environment. The thermal decomposition of the TPS material during the pyrolysis and charring processes is poorly characterized and typically results in large uncertainties in material properties as inputs for ablation models. These material property uncertainties contribute to large design margins on flight systems and cloud re- construction efforts for data collected during flight and ground testing, making revision to existing models for entry systems more challenging. The analysis presented in this work quantifies how material property uncertainties propagate through an ablation model and guides an experimental test regimen aimed at reducing these uncertainties and characterizing the dependencies between properties in the virgin and charred states for a Phenolic Impregnated Carbon Ablator (PICA) based TPS. A sensitivity analysis identifies how the high-fidelity model behaves in the expected flight environment, while a Monte Carlo based uncertainty propagation strategy is used to quantify the expected spread in the in-depth temperature response of the TPS. An examination of how perturbations to the input probability density functions affect output temperature statistics is accomplished using a Kriging response surface of the high-fidelity model. Simulations are based on capsule configuration and aerothermal environments expected during the Mars Science Laboratory (MSL) entry sequence. We identify and rank primary sources of uncertainty from material properties in a flight-relevant environment, show the dependence on spatial orientation and in-depth location on those uncertainty contributors, and quantify how sensitive the expected results are.

Copeland, Sean R.↗

Application of Neural Networks to Wind tunnel Data Response Surface Methods

The integration of nonlinear neural network methods with conventional linear regression techniques is demonstrated for representative wind tunnel force balance data modeling. This work was motivated by a desire to formulate precision intervals for response surfaces produced by neural networks. Applications are demonstrated for representative wind tunnel data acquired at NASA Langley Research Center and the Arnold Engineering Development Center in Tullahoma, TN.

Lo, Ching F.↗

Sustainable bioleaching of lithium-ion batteries for critical metal recovery: Process optimization through design of experiments and thermodynamic modeling

Recycling spent lithium-ion batteries (LIBs) could alleviate supply risks for critical metals and be less harmful to the environment compared to new production of metals from mining. Developing a cost-effective LIB bioleaching process could be a promising alternative to traditional energy-intensive recycling technologies. Here, this study aimed to optimize bioleaching conditions for maximum economic competitiveness through design of experiments using iterative response surface methodology (RSM), assisted by thermodynamic modeling. The optimal condition was identified as 2.5% pulp density in 75 mM gluconic acid biolixiviant at 55°C for 30 h which could recover 57%–84% of nickel, 71%–86% of cobalt, and 100% of lithium and manganese, yielding a 17%–26% net profit margin. The recommended pulp density and acid concentrations, together with the observed metal solubilization, were supported by thermodynamic modeling predictions. Our study demonstrated that combining RSM with thermodynamic simulations could be a powerful tool for optimizing bioleaching conditions.

60 APPLIED LIFE SCIENCES↗

Constraining the Models' Response of Tropical Clouds to SST Forcings Using CALIPSO Observations

Here we present preliminary results from the analysis of the low cloud cover (LCC) and cloud radiative effect (CRE) interannual changes in response to sea surface temperature (SST) forcings in two GISS climate models, and 12 other climate models. We further classify them as a function of their ability to reproduce the vertical structure of the cloud response to SST change against 10 years of CALIPSO observations: “the constrained models, which match the observation constraint, and the unconstrained models”. The constrained models replicate the observed interannual LCC change particularly well (ΔLCC(sub con)=-3.49 ±1.01 %/K vs. ΔLCC(sub obs)=-3.59 ±0.28 %/K) as opposed to the unconstrained models, which largely underestimate it (ΔLCC(sub unc) = -1.32 ± 1.28 %/K). As a result, the amount of short-wave warming simulated by the constrained models (ΔCRE(sub con)=2.60 ±1.13 W/m2/K) is in better agreement with the observations (ΔCRE(sub obs)=3.05 ±0.28 W/m2/K) than the unconstrained models (ΔCRE(sub con)=0.87 ±2.63 W/m2/K). Depending on the type of low cloud, the observed relationship between cloud/radiation and surface temperature varies. Over the stratocumulus regions, increasing SSTs generate higher cloud top height along with a large decrease of the cloud fraction below as opposed to a slight decrease of the cloud fraction at each level over the trade cumulus regions. Our results suggest that the models must generate sustainable stratocumulus decks and moist processes in the planetary boundary layer to reproduce these observed features. Future work will focus on defining a method to objectively discriminate these cloud types that can be applied consistently in both the observations and the models.

Cesana, Gregory↗

Retrieval of Soil Moisture and Roughness from the Polarimetric Radar Response

The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.

Sarabandi, Kamal↗

Physics-based Entry, Descent and Landing Risk Model

A physics-based risk model was developed to assess the risk associated with thermal protection system failures during the entry, descent and landing phase of a manned spacecraft mission. In the model, entry trajectories were computed using a three-degree-of-freedom trajectory tool, the aerothermodynamic heating environment was computed using an engineering-level computational tool and the thermal response of the TPS material was modeled using a one-dimensional thermal response tool. The model was capable of modeling the effect of micrometeoroid and orbital debris impact damage on the TPS thermal response. A Monte Carlo analysis was used to determine the effects of uncertainties in the vehicle state at Entry Interface, aerothermodynamic heating and material properties on the performance of the TPS design. The failure criterion was set as a temperature limit at the bondline between the TPS and the underlying structure. Both direct computation and response surface approaches were used to compute the risk. The model was applied to a generic manned space capsule design. The effect of material property uncertainty and MMOD damage on risk of failure were analyzed. A comparison of the direct computation and response surface approach was undertaken.

descent and landing↗

Direct Geologic Constraints on the Timing of Late Holocene Ice Thickening in the Amundsen Sea Embayment, Antarctica

Abstract Constraining past West Antarctic Ice Sheet (WAIS) change helps validate numerical models simulating future ice sheet dynamics. Following rapid deglaciation during the mid‐Holocene, ice near Thwaites Glacier was ∼35 m thinner than present; however, the timing of ice regrowth to its present configuration remains unknown. To fill this knowledge gap, we present cosmogenic nuclide exposure ages of cobbles from the surface of a moraine situated between Thwaites and Pope glaciers. We infer that the moraine formed and stabilized in the Late Holocene (∼1.4 ka) when a small glacier thickened. We also present a novel reconstruction of WAIS volume constrained by sea‐level data, which demonstrates that moraine formation coincided with a large‐scale WAIS readvance. Our new geologic constraints will help inform models of the solid Earth response to surface mass loading, improving our understanding of ice sheet dynamics in a vulnerable part of WAIS.

Nichols, Keir A. [Imperial College London London U↗

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗

Modeling Time-Dependent Surrogates of Additive-Manufactured Nuclear Fuels Processes

Additive manufacturing (AM) technology is being increasingly adopted in a wide variety of application areas for its ability to rapidly produce, prototype, and customize designs. Recently, a hybrid AM technique was successfully developed at Idaho National Laboratory (INL) to manufacture nuclear fuels [1]. Despite the advantages, this AM technique needs optimization due to defects from a highly complex melting and sintering process. The complex metallurgical phenomena during AM processes are strongly related to parameters such as applied laser power, traveling speed, and scan style, which could lead to differences in density, residual stress, crystallographic texture, and mechanical properties. In addition, stochastic variations in laser energy interaction and associated multiscale/multiphysics phenomena cause variations in microstructure evolution and mechanical properties. Currently, researchers at INL are focusing on developing a comprehensive modeling framework, leveraging INL’s simulation tools MOOSE/MARMOT/BISON/RAVEN [2-4] to describe all steps of this AM process across multiple length scales. Although this advanced framework plays a critical role in enabling enhancements to traditional trial and error approaches for design and optimization of nuclear fuel materials, it remains computationally intense, limiting its use in sensitivity and optimization analysis. In this case, an accurate and inexpensive surrogate becomes an effective tool for providing a tractable approximation of the underlying underline physics. Surrogate models generally not based on the physics of a system are purely mathematical models used to capture the relationships between specific system inputs and outputs. Popular approaches, including neural networks [5], response surfaces [6], and subspace-based reduced order models [7], have been applied to a wide range of disciplines, such as nuclear reactor design, aerospace design and automotive design. In this summary, we employ advanced time-dependent surrogate models such as high-dimensional model representation (HDMR) [8] and physics-informed deep neural network (PINNs) [9] to accelerate the design and optimization of AM process.

42 ENGINEERING↗

Progressive Response Surfaces

Response surface functions are often used as simple and inexpensive replacements for computationally expensive computer models that simulate the behavior of a complex system over some parameter space. Progressive response surfaces are ones that are built up progressively as global information is added from new sample points in the parameter space. As the response surfaces are globally upgraded based on new information, heuristic indications of the convergence of the response surface approximation to the exact (fitted) function can be inferred. Sampling points can be incrementally added in a structured fashion, or in an unstructured fashion. Whatever the approach, at least in early stages of sampling it is usually desirable to sample the entire parameter space uniformly. At later stages of sampling, depending on the nature of the quantity being resolved, it may be desirable to continue sampling uniformly over the entire parameter space (Progressive response surfaces), or to switch to a focusing/economizing strategy of preferentially sampling certain regions of the parameter space based on information gained in early stages of sampling (Adaptive response surfaces). Here we consider Progressive response surfaces where a balanced indication of global response over the parameter space is desired.We use a variant of Moving Least Squares to fit and interpolate structured and unstructured point sets over the parameter space. On a 2-D test problem we compare response surface accuracy for three incremental sampling methods: Progressive Lattice Sampling; Simple-Random Monte Carlo; and Halton Quasi-Monte-Carlo sequences. We are ultimately after a system for constructing efficiently upgradable response surface approximations with reliable error estimates.

Romero, V. J.↗

Response of U.S. West Coast Mountain Snowpack to Local Sea Surface Temperature Perturbations: Insights from Numerical Modeling and Machine Learning

Sea surface temperature (SST) significantly modulates the precipitation and temperature over land, with important consequences on land surface processes such as snowpack. Compared to the impact of remote SST, the effect of nearshore/local SST is less well understood. In this study, the impact of local SST on the mountain snowpack of the U.S. West Coast is investigated using two 6-km regional climate simulations driven by the same lateral boundary conditions but with time-varying versus time-invariant and warmer local SSTs during 2003–15. Results show that local SST warming leads to warmer winter with more precipitation over the mountains. Meanwhile, the removal of SST temporal variability results in reduced temperature variability but increased precipitation variability. As a result, winter snow accumulation decreases by ~200 mm per season in the Cascade Mountains in the north but increases by ~100 mm per season in the Sierra Nevada in the south. Such a dipole response results from the competing effects of precipitation and temperature change at different elevations and are amplified by the enhanced atmospheric river moisture transport. To further delineate the relative contributions of different meteorological factors to the snowpack response, two neural network models were developed to predict the snow behaviors at seasonal and monthly scales. These models reveal the dominant influence of the total amount and the average temperature of precipitation on the snowpack response. Furthermore, these findings highlight the sensitivity of mountain snowpack to local SST in the western United States and underscore the importance of local SST and atmospheric rives to accurate snowpack estimations for water management.

54 ENVIRONMENTAL SCIENCES↗

Effects of surface anomalies on certain model generated meteorological histories

The Mintz-Arakawa 2-level general circulation model has been used in a series of experiments to compute the response of the model atmosphere to: (1) a positive sea surface temperature anomaly in the North Pacific Ocean in summer and in winter, (2) an identical anomaly in the South Pacific Ocean in the Southern Hemisphere winter, and (3) anomalous northward and southward displacements of the Northern Hemisphere snow line over the continents. In each case computations were carried out for 90 forecast days. Results are shown in terms of the differences between anomaly and control histories. Time series of certain regional response indices, including area-average sea level pressure and 600 mb circulation indices, as well as 30-day mean sea level pressure maps are used in the analysis. Of particular interest is the evidence of significant interhemisphere influence.

Spar, J.↗