Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A Hierarchical Evaluation of Regional Climate Simulations

Global climate models (GCMs) are the primary tools for predicting the evolution of the climate system. Through decades of development, GCMs have demonstrated useful skill in simulating climate at continental to global scales. However, large uncertainties remain in projecting climate change at regional scales, which limit our ability to inform decisions on climate change adaptation and mitigation. To bridge this gap, different modeling approaches including nested regional climate models (RCMs), global stretch-grid models, and global high-resolution atmospheric models have been used to provide regional climate simulations (Leung et al. 2003). In previous efforts to evaluate these approaches, isolating their relative merits was not possible because factors such as dynamical frameworks, physics parameterizations, and model resolutions were not systematically constrained. With advances in high performance computing, it is now feasible to run coupled atmosphere-ocean GCMs at horizontal resolution comparable to what RCMs use today. Global models with local refinement using unstructured grids have become available for modeling regional climate (e.g., Rauscher et al. 2012; Ringler et al. 2013). While they offer opportunities to improve climate simulations, significant efforts are needed to test their veracity for regional-scale climate simulations.

climate↗
DOE OSTI2014

Determination of the Fermi level position in dilute magnetic Ga{sub 1-x}Mn{sub x}N films

We report on a combined theoretical and experimental determination of the Fermi level position in wurtzite Ga{sub 1−x}Mn{sub x}N films with x=4% and x=10% as grown by molecular beam epitaxy. By means of ellipsometric measurements, the real part of the frequency-dependent conductivity is determined. An electronic model in the framework of the effective bond-orbital model is parameterized in order to theoretically reproduce the measured transport properties. Predictions for the long-wavelength behaviour as a function of the Fermi level are made. The corresponding density of states obtained in this model is in qualitative agreement with first-principle calculations. The absence of a significant experimental peak in the AC conductivity for small frequencies indicates that the Fermi level lies in a gap between two Mn-related impurity bands in the host band gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Stoichiometrically coupled carbon and nitrogen cycling in the MIcrobial-MIneral Carbon Stabilization model version 1.0 (MIMICS-CN v1.0)

Explicit consideration of microbial physiology in soil biogeochemical models that represent coupled carbon–nitrogen dynamics presents opportunities to deepen understanding of ecosystem responses to environmental change. The MIcrobial-MIneral Carbon Stabilization (MIMICS) model explicitly represents microbial physiology and physicochemical stabilization of soil carbon (C) on regional and global scales. Here we present a new version of MIMICS with coupled C and nitrogen (N) cycling through litter, microbial, and soil organic matter (SOM) pools. The model was parameterized and validated against C and N data from the Long-Term Inter-site Decomposition Experiment Team (LIDET; six litter types, 10 years of observations, and 13 sites across North America). The model simulates C and N losses from litterbags in the LIDET study with reasonable accuracy (C: R 2 =0.63; N: R 2 =0.29), which is comparable with simulations from the DAYCENT model that implicitly represents microbial activity (C: R 2 =0.67; N: R 2 =0.30). Subsequently, we evaluated equilibrium values of stocks (total soil C and N, microbial biomass C and N, inorganic N) and microbial process rates (soil heterotrophic respiration, N mineralization) simulated by MIMICS-CN across the 13 simulated LIDET sites against published observations from other continent-wide datasets. We found that MIMICS-CN produces equilibrium values in line with measured values, showing that the model generates plausible estimates of ecosystem soil biogeochemical dynamics across continental-scale gradients. MIMICS-CN provides a platform for coupling C and N projections in a microbially explicit model, but experiments still need to identify the physiological and stoichiometric characteristics of soil microbes, especially under environmental change scenarios.

58 GEOSCIENCES↗
DOE OSTI2026

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Centralized Activities in Support of the Radiative Forcing Model Intercomparison Project

A new paradigm in benchmark absorption-scattering radiative transfer is presented that enables both the globally-averaged and spatially-resolved testing of climate model radiation parameterizations using a novel benchmark computational strategy to uncover persistent sources of biases in the aerosol Instantaneous Radiative Effect (IRE). A proof-of-concept is demonstrated with the GFDL AM4 and CESM 1.2.2 climate models. Instead of prescribing atmospheric conditions and aerosols, as in prior intercomparisons, snapshots of the atmospheric state and aerosol optical properties from the participating models are used as inputs to a highly accurate radiation solver to uncover relevant biases unique to each climate model. These diagnostic results show that the models’ aerosol IRE bias is of the same magnitude as the persistent intermodel spread (~1 W/m2), and varies spatially and with intrinsic aerosol optical properties. The findings underscore the significance of native model error analysis and its ability to diagnose global biases, confirming its fundamental value for the Radiative Forcing Model Intercomparison Project (RFMIP).

58 GEOSCIENCES↗
DOE OSTI2017

Error modeling for surrogates of dynamical systems using machine learning: Machine-learning-based error model for surrogates of dynamical systems

A machine learning–based framework for modeling the error introduced by surrogate models of parameterized dynamical systems is proposed. The framework entails the use of high-dimensional regression techniques (eg, random forests, and LASSO) to map a large set of inexpensively computed “error indicators” (ie, features) produced by the surrogate model at a given time instance to a prediction of the surrogate-model error in a quantity of interest (QoI). This eliminates the need for the user to hand-select a small number of informative features. The methodology requires a training set of parameter instances at which the time-dependent surrogate-model error is computed by simulating both the high-fidelity and surrogate models. Using these training data, the method first determines regression-model locality (via classification or clustering) and subsequently constructs a “local” regression model to predict the time-instantaneous error within each identified region of feature space. We consider 2 uses for the resulting error model: (1) as a correction to the surrogate-model QoI prediction at each time instance and (2) as a way to statistically model arbitrary functions of the time-dependent surrogate-model error (eg, time-integrated errors). We then apply the proposed framework to model errors in reduced-order models of nonlinear oil-water subsurface flow simulations, with time-varying well-control (bottom-hole pressure) parameters. The reduced-order models used in this work entail application of trajectory piecewise linearization in conjunction with proper orthogonal decomposition. Moreover, when the first use of the method is considered, numerical experiments demonstrate consistent improvement in accuracy in the time-instantaneous QoI prediction relative to the original surrogate model, across a large number of test cases. When the second use is considered, results show that the proposed method provides accurate statistical predictions of the time- and well-averaged errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗
DOE OSTI2024

Non-standard neutrino interactions mediated by a light scalar at DUNE

Abstract We investigate the effect on neutrino oscillations generated by beyond-the-standard-model interactions between neutrinos and matter. Specifically, we focus on scalar-mediated non-standard interactions (NSI) whose impact fundamentally differs from that of vector-mediated NSI. Scalar NSI contribute as corrections to the neutrino mass matrix rather than the matter potential and thereby predict distinct phenomenology from the vector-mediated ones. Similar to vector-type NSI, the presence of scalar-mediated neutrino NSI can influence measurements of oscillation parameters in long-baseline neutrino oscillation experiments, with a notable impact on CP measurement in the case of DUNE. Our study focuses on the effect of scalar NSI on neutrino oscillations, using DUNE as an example. We introduce a model-independent parameterization procedure that enables the examination of the impact of all non-zero scalar NSI parameters simultaneously. Subsequently, we convert DUNE’s sensitivity to the NSI parameters into projected sensitivity concerning the parameters of a light scalar model. We compare these results with existing non-oscillation probes. Our findings reveal that the region of the light scalar parameter space sensitive to DUNE is predominantly excluded by non-oscillation probes, especially when considering all nonzero parameters simultaneously for DUNE.

Physics↗
DOE OSTI2015

A preliminary regional PBPK model of lung metabolism for improving species dependent descriptions of 1,3-butadiene and its metabolites

1,3-Butadiene (BD), a volatile organic chemical (VOC), is used in synthetic rubber production and other industrial processes. It is detectable at low levels in ambient air as well as in tobacco smoke and gasoline vapors. Inhalation exposures to high concentrations of BD have been associated with lung cancer in both humans and experimental animals, although differences in species sensitivity have been observed. Metabolically active lung cells such as Pulmonary Type I and Type II epithelial cells and club cells (Clara cells) 1 are potential targets of BD metabolite-induced toxicity. Metabolic capacities of these cells, their regional densities, and distributions vary throughout the respiratory tract as well as between species and cell types. Here we present a physiologically based pharmacokinetic (PBPK) model for BD that includes a regional model of lung metabolism, based on a previous model for styrene, to provide species-dependent descriptions of BD metabolism in the mouse, rat, and human. Since there are no in vivo data on BD pharmacokinetics in the human, the rat and mouse models were parameterized to the extent possible on the basis of in vitro metabolic data. Where it was necessary to use in vivo data, extrapolation from rat to mouse was performed to evaluate the level of uncertainty in the human model. A kidney compartment and description of downstream metabolism were also included in the model to allow for eventual use of available urinary and blood biomarker data in animals and humans to calibrate the model for estimation of BD exposures and internal metabolite levels. Results from simulated inhalation exposures to BD indicate that incorporation of differential lung region metabolism is important in describing species differences in pulmonary response and that these differences may have implications for risk assessments of human exposures to BD.

1,3-Butadiene↗
DOE OSTI2026

On the Prospect of Chemically Transferable Coarse-Grained Electronic Models for Soft Materials

Electronic coarse-graining (ECG) methods predict quantum-mechanical electronic properties directly from coarse-grained (CG) molecular configurations, enabling electronic predictions at mesoscale length scales. Here, we present a diagnostic assessment of the feasibility of chemically transferable ECG models across a broad polymer-relevant chemical space using all-atom, united-atom, and Martini-scale representations. While high-resolution ECG models achieve near-quantitative accuracy, we show that chemically transferable ECG at the Martini resolution fails because the CG force field does not sample the same configurational distribution of local molecular structure as that underlying the DFT-parameterized ECG model. We demonstrate that our proposed Element-Count-Label (ECL) representation, which augments Martini beads with explicit stoichiometric data, significantly improves chemical generalization across diverse polymer chemistries. However, we find that even with improved chemical resolution, the model cannot recover electronic property distributions that are absent from the configurational space sampled by the CG force field. These results demonstrate that chemically transferable ECG requires future Martini-like force fields to explicitly preserve quantum chemistry–compatible local molecular structure in addition to thermodynamic and structural fidelity.

Kidder, Katherine M [Department of Chemistry; Univ↗
DOE OSTI2021

Towards Realistic Models of Heterogeneous Catalysis: Simulations of Oxidation Catalysis from First Principles

Among the most significant developments in heterogeneous catalysis in the last 20 years is the emergence of microkinetic models, often parameterized from density functional theory (DFT) calculations, used to quantify observed catalyst performance and to guide the discovery of new catalytic materials. While DFT directly reports binding energies and elementary step activation energies, the free energies that enter into microkinetic models must be computed from additional approximations. Frequently these free energy approximations assume ideal behavior—that adsorbates do not interact with one another, or that adsorbates are immobile or vibrate harmonically about a binding site. These assumptions can and do have an impact on predicted catalyst performance and potentially even on predicted trends. Our work relates to three categories of non-idealities, explored in the context of nitrogen oxidation and reduction catalysis on metal surfaces. One component relates to adsorbate translational free energy, the contribution least well described by conventional models, in which we develop modeling approaches that improve accuracy with limited increase in computational expense. A second component relates adsorbate-adsorbate interactions, in which we develop benchmark on-lattice interaction models and compare kinetic predictions with conventional mean-field models, in an effort to develop more robust coverage-dependent mean field modeling approaches. The last relates the most fundamental assumption of all—that of energy equipartition—in an effort to rationalize and guide plasma-enhanced catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Microkinetic modeling of aqueous phase biomass conversion: Application to ethylene glycol reforming

In this communication, we develop automatically a reaction network for ethylene glycol reforming on Pt. We employ a recently introduced group additivity scheme for estimation of thermochemistry over metal catalysts in vapor phase and in solution along with Brønsted-Evans-Polanyi (BEP) relations for vapor phase kinetics to parameterize microkinetic models for vapor and aqueous phases. Unlike vapor-phase reforming, we show that solvent occupies a significant fraction of surface sites and by doing this, it shifts the water-gas shift reaction to CO2 and H2 and avoids poisoning of the catalyst by CO. The model is found to be in reasonable agreement with the experimental data.

bio-inspired↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2020

Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models

Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Improved physiologically based pharmacokinetic model for oral exposures to chromium in mice, rats, and humans to address temporal variation and sensitive populations

A physiologically based pharmacokinetic (PBPK) model for hexavalent chromium [Cr(VI)] in mice, rats, and humans developed previously (Kirman et al., 2012, 2013), was updated to reflect an improved understanding of the toxicokinetics of the gastrointestinal tract following oral exposures. Improvements were made to: (1) the reduction model, which describes the pH-dependent reduction of Cr(VI) to Cr(III) in the gastrointestinal tract under both fasted and fed states; (2) drinking water pattern simulations, to better describe dosimetry in rodents under the conditions of the NTP cancer bioassay; and (3) parameterize the model to characterize potentially sensitive human populations. Important species differences, sources of non-linear toxicokinetics, and human variation are identified and discussed within the context of human health risk assessment. - Highlights: • An improved version of the PBPK model for Cr(VI) toxicokinetics was developed. • The model incorporates data collected to fill important data gaps. • Model predictions for specific age groups and sensitive subpopulations are provided. • Implications to human health risk assessment are discussed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Evaluating the Future of Lightning in Cloud‐Resolving Models

Abstract Two proxies for lightning predict very different responses to global warming: the CAPE times precipitation proxy predicts a large increase in lightning over both the continental United States and the tropical oceans, while the ice flux proxy predicts a small increase over the United States and a decrease over the tropical oceans. To date, however, these proxies have been studied only in global climate models with parameterized convection. Here, cloud‐resolving simulations are used to assess their predictions of future lightning rates. Over the United States, all proxies predict a large increase in the lightning rate in the range of 8–16%/K. On the other hand, in the tropics as modeled by radiative convective equilibrium, half of the proxies predict an increase (of 5–12%/K), while the other half predict a decrease (of 1–4%/K). The reasons for the different responses of these proxies is explored, but it remains unclear which proxy is best suited to predicting future lightning rates.

Romps, David M.↗