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 145 records · Page 8

The Effects of an Improved Dynamic Vegetation Phenology Representation in a Global Land Surface Model

Evapotranspiration (ET) is a major driver of the interaction between the land surface and the atmosphere through its component mechanisms, including plant transpiration (T) and soil evaporation. To accurately capture land-atmosphere interactions in global Earth System Models, it is thus critical that the underlying land surface models accurately model both the land hydrology as well as the dynamic response of vegetation to environmental drivers. In an effort to introduce a more realistic vegetation representation, the NASA Catchment land surface model, which is part of the Goddard Earth Observing System (GEOS), has previously been merged with the carbon and nitrogen physics modules of the Community Land Model version 4, resulting in the new Catchment-CN model. Catchment-CN has inherited the advanced treatment of land surface hydrology of Catchment, but is additionally able to dynamically model the response of vegetation to environmental drivers, in contrast to the fixed vegetation climatology that was prescribed in Catchment. Recently, the parameterization of Catchment-CN vegetation has been augmented to better account for variability of vegetation responses to external forcings within existing plant functional types, and vegetation parameters have been calibrated against Moderate Resolution Imaging Spectrometer observations of the fraction of absorbed photosynthetically radiation. These efforts have led to a significant reduction in the RMSE of modeled photosynthetic activity with respect to observations.This presentation investigates the effect of the improved vegetation representation on the partitioning of ET within Catchment-CN. Specifically, we compare global maps of the T:ET ratio across different temporal scales in (1) the original Catchment model, (2) the original Catchment-CN model, and (3) the augmented and calibrated Catchment-CN model. The modeled T and ET estimates are compared against a comprehensive set of ground observations from various field studies, as well as independent global T:ET estimates from different ET algorithms provided in the context of the Water Cycle Observation Multi-mission Strategy ? Evapotranspiration (WACMOS-ET) initiative.

Kolassa, J.↗

Local and Remote Mean and Extreme Temperature Response to Regional Aerosol Emissions Reductions

The climatic implications of regional aerosol and precursor emissions reductions implemented to protect human health are poorly understood. We investigate the mean and extreme temperature response to regional changes in aerosol emissions using three coupled chemistry–climate models: NOAA GFDL CM3, NCAR CESM1, and NASA GISS-E2. Our approach contrasts a long present-day control simulation from each model (up to 400 years with perpetual year 2000 or 2005 emissions) with 14 individual aerosol emissions perturbation simulations (160–240 years each). We perturb emissions of sulfur dioxide (SO2) and/or carbonaceous aerosol within six world regions and assess the statistical significance of mean and extreme temperature responses relative to internal variability determined by the control simulation and across the models. In all models, the global mean surface temperature response (perturbation minus control) to SO2 and/or carbonaceous aerosol is mostly positive (warming) and statistically significant and ranges from +0.17 K (Europe SO2) to -0.06 K (US BC). The warming response to SO2 reductions is strongest in the US and Europe perturbation simulations, both globally and regionally, with Arctic warming up to 1 K due to a removal of European anthropogenic SO2 emissions alone; however, even emissions from regions remote to the Arctic, such as SO2 from India, significantly warm the Arctic by up to 0.5 K. Arctic warming is the most robust response across each model and several aerosol emissions perturbations. The temperature response in the Northern Hemisphere midlatitudes is most sensitive to emissions perturbations within that region. In the tropics, however, the temperature response to emissions perturbations is roughly the same in magnitude as emissions perturbations either within or outside of the tropics. We find that climate sensitivity to regional aerosol perturbations ranges from 0.5 to 1.0 K (W m(exp -2))(exp -1) depending on the region and aerosol composition and is larger than the climate sensitivity to a doubling of CO2 in two of three models. We update previous estimates of regional temperature potential (RTP), a metric for estimating the regional temperature responses to a regional emissions perturbation that can facilitate assessment of climate impacts with integrated assessment models without requiring computationally demanding coupled climate model simulations. These calculations indicate a robust regional response to aerosol forcing within the Northern Hemisphere midlatitudes, regardless of where the aerosol forcing is located longitudinally. We show that regional aerosol perturbations can significantly increase extreme temperatures on the regional scale. Except in the Arctic in the summer, extreme temperature responses largely mirror mean temperature responses to regional aerosol perturbations through a shift of the temperature distributions and are mostly dominated by local rather than remote aerosol forcing.

Westervelt, Daniel M.↗

Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry

To provide theoretical basis for the connection between observed radar scattering and wind-generated waves, a model for the response of surface waves in the gravity-capillary equilibrium region of the spectrum is proposed on the basis of a local (in wavenumber) balance between wind input and dissipation. The wind input function was constructed on the basis of laboratory observations of short-wave growth, while the dissipation function was developed from ideas of viscous dissipation and wave breaking in response to local accelerations and modified by kinematic effects of phase and group velocity differences. The model was exercised at L, C, X, and Ka bands to demonstrate the differences in wind speed and water temperature sensitivity.

Donelan, Mark A.↗

Implementing Geometric Surface Imperfections into Sandwich Composite Cylinder Finite Element Method Models

The buckling responses of certain cylindrical shell structures are extremely sensitive to geometric surface imperfections. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF) is conducting research to develop analysis-based buckling design recommendations. Experiments are used to verify the analysis-based factors, but the sensitivity of the test articles to geometric imperfections requires implementing as-manufactured imperfections into high-fidelity finite element method models. Data collection methods such as structured light scanning are used for all geometric surface data used in this work. Common preprocessing and visualization steps used in SBKF are discussed, and steps on how surface scans are prepared for implementation into a finite element model is described. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS), written specifically for the use with SBKF, is briefly described and uses eight functions to extract, modify, and write geometric imperfections into Abaqus input files. Results of the pre-processing methods and results from Py_TIGIRS are provided and compared for Composite Test Article (CTA) 8.2B. Excellent agreement between the visualized scan data and the FEM-extracted geometry is demonstrated. A brief example of why geometric surface imperfections are significant in nonlinear numerical analyses for thin cylinders in axial compression is provided as motivation to use tools such as Py_TIGIRS. Future development of Py_TIGIRS including expansion to structures of arbitrary geometry is planned.

Sandwich structures↗

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗

Topology optimization of an airfoil fin microchannel heat exchanger using artificial intelligence

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, is analyzed and optimized using a design of experiments with artificial intelligence techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin printed circuit heat exchangers (PCHEs) has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure is developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. Computational fluid dynamics is used to analyze a simplified two-channel configuration where five design parameters are varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods (a 5D polynomial and a regression neural network) are compared for generating surrogate models and the resulting response surface approximation is input to a genetic algorithm that is used to identify a set of optimal parameters. The optimal geometries are found across six channel Reynolds numbers ranging from 1000 to 5000, since inlet conditions affect flow through the heat exchanger. Additionally, a set of optimal designs that maximizes heat transfer and minimizes pressure drop is identified, and a thermal stress analysis is performed on the optimal design. Correlations for the Nusselt number and Darcy friction factor are developed that can be useful for thermal hydraulic analyses using system codes. Thermal stresses are analyzed and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

Correlations Between Sea-Surface Salinity Tendencies and Freshwater Fluxes in the Pacific Ocean

Temporal changes in sea-surface salinity (SSS) from 21 years of a high resolution model integration of the Pacific Ocean are correlated with the freshwater flux that was used to force the integration. The correlations are calculated on a 1 x10 grid, and on a monthly scale to assess the possibility of deducing evaporation minus precipitation (E-P) fields from the salinity measurements to be taken by the upcoming Aquarius/SAC-D mission. Correlations between the monthly mean E-P fields and monthly mean SSS temporal tendencies are mainly zonally-oriented, and are highest where the local precipitation is relatively high. Nonseasonal (deviations from the monthly mean) correlations are highest along mid-latitude storm tracks and are relatively small in the tropics. The response of the model's surface salinity to surface forcing is very complex, and retrievals of freshwater fluxes from SSS measurements alone will require consideration of other processes, including horizontal advection and vertical mixing, rather than a simple balance between the two.

Li, Zhen↗

Transequatorial effects of sea surface temperature anomalies in a global general circulation model

The global response of the atmosphere to a persistent anomalous pool of warm sea surface temperatures in the extratropical Pacific Ocean is examined in terms of the meridional pole-to-pole profile of the zonally-averaged 600 mb surface for periods up to 90 days. Following an initial hydrostatic inflation of the isobaric surface in the latitude of the warm pool, effects spread poleward within the hemisphere, then begin to appear after about two to three weeks in high latitudes of the opposite hemisphere, but with little or no response in the tropics. The same sea temperature anomaly field generates a stronger response in winter than in summer, and a very different reaction when located in the Southern Hemisphere than when in the Northern Hemisphere.

Spar, J.↗

Comparison of Entry Descent and Landing Aerodynamic Databases with Uncertainty Quantification Developed Using Machine Learning Techniques

When developing the aerodynamic databases for use in trajectory simulations, it is important to develop a system of metrics to qualify which aerodynamic models are best to use. Since aerodynamics are just one input into trajectory simulations, the results of these simulations do not reflect on the quality of the aerodynamic database used. This means that aerodynamic database comparisons must be done offline. While traditional metrics that focus on mean/nominal predictions are a good first step, more robust estimates of the prediction interval become important as more focused uncertainty models are developed. We explore the limitations of evaluating aerodynamic models based purely on nominal-centered response surfaces. Before elaborating and evaluating metrics based on distributed models, the value of evaluating prediction interval and confidence interval are discussed to conclude that prediction intervals are more relevant to the use of trajectory analysis. Several metrics to evaluate the prediction interval are introduced with a focus on the standard calibration metric. Finally, we compare candidate models using both mean and distributed metrics. A finalized candidate model developed using state of the art machine learning methods is compared to a baseline model developed using traditional aerodynamic database modeling techniques.

Aerodynamic Database↗

Epidermal Homeostasis and Radiation Responses in a Multiscale Tissue Modeling Framework

The surface of skin is lined with several thin layers of epithelial cells that are maintained throughout life time by a small population of stem cells. High dose radiation exposures could injure and deplete the underlying proliferative cells and induce cutaneous radiation syndrome. In this work we propose a multiscale computational model for skin epidermal dynamics that links phenomena occurring at the subcellular, cellular, and tissue levels of organization, to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis. Incorporating experimentally measured histological and cell kinetic parameters, we obtain results of population kinetics and proliferation indexes comparable to observations in unirradiated and acutely irradiated swine experiments. At the sub-cellular level, several recently published Wnt signaling controlled cell-cycle models are applied and the roles of key components and parameters are analyzed. Based on our simulation results, we demonstrate that a moderate increase of proliferation rate for the survival proliferative cells is sufficient to fully repopulate the area denuded by high dose radiation, as long as the integrity of underlying basement membrane is maintained. Our work highlights the importance of considering proliferation kinetics as well as the spatial organization of tissues when conducting in vivo investigations of radiation responses. This integrated model allow us to test the validity of several basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and enhance our understanding of the pathophysiological effects of ionizing radiation on skin.

Hu, Shaowen↗

Interactive Biogenic Emissions and Drought Stress Effects on Atmospheric Composition in NASA GISS ModelE

Drought is a hydroclimatic extreme that causes perturbations to the terrestrial biosphere and acts as a stressor on vegetation, affecting emissions patterns. During severe drought, isoprene emissions are reduced. In this paper, we focus on capturing this reduction signal by implementing a new percentile isoprene drought stress (y d ) algorithm in NASA GISS ModelE based on the MEGAN3 (Model of Emissions of Gases and Aerosols from Nature Version 3) approach as a function of a photosynthetic parameter (V c,max ) and water stress (β). Four global transient simulations from 2003–2013 are used to demonstrate the effect without y d (Default_ModelE) and with online y d (DroughtStress_ModelE). DroughtStress_ModelE is evaluated against the observed isoprene measurements at the Missouri Ozarks AmeriFlux (MOFLUX) site during the 2012 severe drought where improvements in the correlation coefficient indicate it is a suitable drought stress parameterization to capture the reduction signal during severe drought. The application of y d globally leads to a decadal average reduction of ∼ 2.7 %, which is equivalent to ∼14.6 Tg yr -1 of isoprene. The changes have larger impacts in regions such as the southeastern US. DroughtStress_ModelE is validated using the satellite ΩHCHO column from the Ozone Monitoring Instrument (OMI) and surface O 3 observations across regions of the US to examine the effect of drought on atmospheric composition. It was found that the inclusion of isoprene drought stress reduced the overestimation of ΩHCHO in Default_ModelE during the 2007 and 2011 southeastern US droughts and led to improvements in simulated O 3 during drought periods. We conclude that isoprene drought stress should be tuned on a model-by-model basis because the variables used in the parameterization responses are relative to the land surface model hydrology scheme (LSM) and the effects of y d application could be larger than seen here due to ModelE not having large biases of isoprene during severe drought.

drought stress effects↗

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

Using experimental design modules for process characterization in manufacturing/materials processes laboratories

Modules dealing with statistical experimental design (SED), process modeling and improvement, and response surface methods have been developed and tested in two laboratory courses. One course was a manufacturing processes course in Mechanical Engineering and the other course was a materials processing course in Materials Science and Engineering. Each module is used as an 'experiment' in the course with the intent that subsequent course experiments will use SED methods for analysis and interpretation of data. Evaluation of the modules' effectiveness has been done by both survey questionnaires and inclusion of the module methodology in course examination questions. Results of the evaluation have been very positive. Those evaluation results and details of the modules' content and implementation are presented. The modules represent an important component for updating laboratory instruction and to provide training in quality for improved engineering practice.

Ankenman, Bruce↗

Quantifying Parameter Sensitivity, Interaction and Transferability in Hydrologically Enhanced Versions of Noah-LSM over Transition Zones

We use sensitivity analysis to identify the parameters that are most responsible for shaping land surface model (LSM) simulations and to understand the complex interactions in three versions of the Noah LSM: the standard version (STD), a version enhanced with a simple groundwater module (GW), and version augmented by a dynamic phenology module (DV). We use warm season, high-frequency, near-surface states and turbulent fluxes collected over nine sites in the US Southern Great Plains. We quantify changes in the pattern of sensitive parameters, the amount and nature of the interaction between parameters, and the covariance structure of the distribution of behavioral parameter sets. Using Sobol s total and first-order sensitivity indexes, we show that very few parameters directly control the variance of the model output. Significant parameter interaction occurs so that not only the optimal parameter values differ between models, but the relationships between parameters change. GW decreases parameter interaction and appears to improve model realism, especially at wetter sites. DV increases parameter interaction and decreases identifiability, implying it is overparameterized and/or underconstrained. A case study at a wet site shows GW has two functional modes: one that mimics STD and a second in which GW improves model function by decoupling direct evaporation and baseflow. Unsupervised classification of the posterior distributions of behavioral parameter sets cannot group similar sites based solely on soil or vegetation type, helping to explain why transferability between sites and models is not straightforward. This evidence suggests a priori assignment of parameters should also consider climatic differences.

Rosero, Enrique↗

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗

Many‐body response of benzene at monolayer MoS 2 : Van der Waals interactions and spectral broadening

Abstract Models of surface enhancement of molecular electronic response properties are challenging for two reasons: (a) molecule‐surface interactions require a simultaneous solution of the molecular and the surface dynamic response (a daunting task), and (b) when solving for the electronic structure of the combined molecule + surface system, it is not trivial to single out the particular physical effects responsible for enhancement. To tackle this problem, in this work, we apply a formally exact decomposition of the system's response function into subsystem contributions by using subsystem density functional theory (DFT), which grants access to dynamic polarizabilities and optical spectra. In order to access information about the interactions between the subsystems, we extend a previously developed subsystem‐based adiabatic connection fluctuation‐dissipation theorem of DFT to separate the additive from the nonadditive correlation energy and identify the nonadditive correlation as the van der Waals interactions. As an example, we choose benzene adsorbed on monolayer MoS 2 . We isolate the contributions to benzene's dynamic response arising from the interaction with the surface, and for the first time, we evaluate the enhancements to the effectiveness of C 6 coefficients as a function of benzene‐MoS 2 distance and adsorption site. We also quantify the spectral broadening of the benzene's electronic excited states due to their interaction with the surface. We find that the broadening has a similar decay law with the molecule‐surface distance as the leading van der Waals interactions (ie, R −6 ) and that the surface enhancement of dispersion interactions between benzene molecules is less than 5% but is still large enough (0.5 kcal/mol) to likely play a role in the prediction of interface morphologies.

Umerbekova, Alina↗

Transequatorial effects of sea-surface temperature anomalies in a global general circulation model

An effort is made to describe the development of the model atmosphere's response to the sea-surface temperature (SST) anomalies in terms of their effect on the pole-to-pole meridional profile of the zonal mean 600-mb surface, which represents approximately the middle level of the model. The initial effect of a positive SST anomaly located in extratropical latitudes of the Pacific Ocean is found to be an inflation of the isobaric surfaces, relative to a control run, within the latitude band of the warm pool. In experiments with a North Pacific warm pool, the inflation continued for about 3 weeks in both winter and summer.

Spar, J.↗

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗