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 OSTI2016

EVOLUTIONARY TRACKS OF THE CLIMATE OF EARTH-LIKE PLANETS AROUND DIFFERENT MASS STARS

The climatic evolution of the Earth depends strongly on the evolution of the insolation from the Sun and the amount of the greenhouse gasses, especially CO{sub 2} in the atmosphere. Here, we investigate the evolution of the climate of hypothetical Earths around stars whose masses are different from the solar mass with a luminosity evolution model of the stars, a mantle degassing model coupled with a parameterized convection model of the planetary interiors, and an energy balance climate model of the planetary surface. In the habitable zone (HZ), the climate of the planets is initially warm or hot, depending on the orbital semimajor axes. We found that, in the inner HZ, the climate of the planets becomes hotter with time owing to the increase in the luminosity of the central stars, while, in the outer HZ, it becomes colder and eventually globally ice-covered owing to the decrease in the CO{sub 2} degassing rate of the planets. The orbital condition for maintaining the warm climate similar to the present Earth becomes very limited, and more interestingly, the planet orbiting in the outer HZ becomes globally ice-covered after a certain critical age (∼3 Gyr for the hypothetical Earth with standard parameters), irrespective of the mass of the central star. This is because the critical age depends on the evolution of the planets and planetary factors, rather than on the stellar mass. The habitability of the Earth-like planet is shown to be limited with age even though it is orbiting within the HZ.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗
NASA NTRS2003

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗
DOE OSTI2026

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
NASA NTRS2020

Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective56equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climatecsensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique amongst intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than un-aggregated simulations. Plain Language Summary This study investigates tropical clouds and climate using results from more than 30 different numerical models set up in a simplified framework. The dataset of model simulations is unique in that it includes a wide range of model types configured in a consistent manner. We address some of the biggest open questions in climate science, including how cloud properties change with warming and the role that the tendency of clouds to form clusters plays in determining the average climate and how climate changes. While there are large differences in how the different models simulate average temperature, humidity, and cloudiness, in a majority of models, the amount of high clouds decreases as climate warms. Nearly all models simulate a tendency for clouds to cluster together. There is agreement that when the clouds are clustered, the atmosphere is drier with fewer clouds overall. We don’t find a conclusive result for how cloud clustering changes as the climate warms.

Allison A. Wing↗
DOE OSTI2014

Simulating weather regimes: impact of model resolution and stochastic parameterization

The simulation of quasi-persistent regime structures in an atmospheric model with horizontal resolution typical of the Intergovernmental Panel on Climate Change fifth assessment report simulations, is shown to be unrealistic. A higher resolution configuration of the same model, with horizontal resolution typical of that used in operational numerical weather prediction, is able to simulate these regime structures realistically. The spatial patterns of the simulated regimes are remarkably accurate at high resolution. A model configuration at intermediate resolution shows a marked improvement over the low-resolution configuration, particularly in terms of the temporal characteristics of the regimes, but does not produce a simulation as accurate as the very-high-resolution configuration. It is demonstrated that the simulation of regimes can be significantly improved, even at low resolution, by the introduction of a stochastic physics scheme. At low resolution the stochastic physics scheme drastically improves both the spatial and temporal aspects of the regimes simulation. These results highlight the importance of small-scale processes on large-scale climate variability, and indicate that although simulating variability at small scales is a necessity, it may not be necessary to represent the small-scales accurately, or even explicitly, in order to improve the simulation of large-scale climate. It is argued that these results could have important implications for improving both global climate simulations, and the ability of high-resolution limited-area models, forced by low-resolution global models, to reliably simulate regional climate change signals.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evaluation of PBL Parameterization Schemes in WRF Model Predictions during the Dry Season of the Central Amazon Basin

Planetary Boundary Layer (PBL) parameterization schemes are employed to handle subgrid-scale processes on atmospheric models, playing a key role in accurately representing the atmosphere. Recent studies have shown that PBL schemes are particularly fundamental to the depiction of PBL height (PBLH), especially over the Amazon. In the present study, we investigated the performance of PBL schemes on the representation of meteorological variables, turbulent fluxes, PBL vertical structures, and PBLH over the central Amazon basin under dry conditions, taking advantage of observations from the Observations and Modeling of the Green Ocean Amazon campaign (GoAmazon2014/5) for validation and evaluation. Numerical experiments were carried out within the WRF model using eight PBL schemes for two dry periods from 2014 (typical year) and 2015 (El-Niño year), and results from the 1-km resolution domain were directly compared to hourly in situ observations. In general, all PBL schemes present good performance to reproduce meteorological variables, with nonlocal (local) PBL schemes producing better performance in the 2014 (2015) study period. All PBL schemes in general overestimate (>100%) daytime turbulent fluxes. Thermodynamic (daytime) vertical structures are better predicted than mechanical (nocturnal) ones. The local MYNN2.5 scheme showed the overall best performance for PBLH prediction, mainly at night.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Implementation and Evaluation of Physics-Driven Dynamic Entrainment-Mixing Parameterization in a Climate Model and Its Impact on Low-Cloud Simulation

The turbulent entrainment-mixing process in the Community Earth System Model version 1.2 (CESM1.2) is assumed to follow the extremely inhomogeneous entrainment-mixing. However, different entrainment-mixing scenarios can occur in real clouds. To address this deficiency, a unifying parameterization that represents different entrainment-mixing processes is implemented and evaluated in CESM1.2. The results indicate that the homogeneous mixing degree values simulated by the new parameterization in CESM1.2 are predominantly greater than 50%, suggesting a tendency toward homogeneous mixing. Compared to the extremely inhomogeneous mixing mechanism, the new parameterization increases the cloud droplet number concentration (Nc). More importantly, the new parameterization improves low-cloud fraction (CLDLOW) simulation in Northwest Pacific (NWP) and Southeast Pacific (SEP) regions, with relative improvements of 2.95% and 4.17%, respectively. Furthermore, the improvements reach up to 44.6% and 16.2% in the NWP and SEP regions, respectively, when considering the relationship between N c and CLDLOW. Further analysis reveals that the new parameterization enhances cloud optical depth, longwave radiative cooling effect, net condensation rate, cloud water mixing ratio, lower-troposphere stability, and CLDLOW by increasing N c . Additionally, these results underscore the importance of improving entrainment-mixing parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Parameterizing deep convection using the assumed probability density function method

Due to their coarse horizontal resolution, present-day climate models must parameterize deep convection. This paper presents single-column simulations of deep convection using a probability density function (PDF) parameterization. The PDF parameterization predicts the PDF of subgrid variability of turbulence, clouds, and hydrometeors. That variability is interfaced to a prognostic microphysics scheme using a Monte Carlo sampling method.The PDF parameterization is used to simulate tropical deep convection, the transition from shallow to deep convection over land, and mid-latitude deep convection.These parameterized single-column simulations are compared with 3-D reference simulations. The agreement is satisfactory except when the convective forcing is weak. The same PDF parameterization is also used to simulate shallow cumulus and stratocumulus layers. The PDF method is sufficiently general to adequately simulate these five deep, shallow, and stratiform cloud cases with a single equation set. This raises hopes that it may be possible in the future, with further refinements at coarse time step and grid spacing, to parameterize all cloud types in a large-scale model in a unified way.

Storer, R. L.↗
DOE OSTI2020

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

ExaSat: An Exascale Static Analysis Tool (ExaSat) v1.0

The ExaSAT framework automates the extraction of parameterized performance models directly from source code using compiler analysis. Its parameterized analytic model enables quantitative evaluation of a broad range of hardware design trade-offs and software optimizations on a variety of different performance metrics, with a primary focus on data movement as a metric. ExaSAT analysis provides insights in the hardware and software tradeoffs and lays the groundwork for exploring a more targeted set of design points using cycle-accurate architectural simulators and emulators.

Chan, Cy↗
DOE OSTI2015

ExaSAT: An exascale co-design tool for performance modeling

One of the emerging challenges to designing HPC systems is understanding and projecting the requirements of exascale applications. In order to determine the performance consequences of different hardware designs, analytic models are essential because they can provide fast feedback to the co-design centers and chip designers without costly simulations. However, current attempts to analytically model program performance typically rely on the user manually specifying a performance model. Here we introduce the ExaSAT framework that automates the extraction of parameterized performance models directly from source code using compiler analysis. The parameterized analytic model enables quantitative evaluation of a broad range of hardware design trade-offs and software optimizations on a variety of different performance metrics, with a primary focus on data movement as a metric. Finally, we demonstrate the ExaSAT framework’s ability to perform deep code analysis of a proxy application from the Department of Energy Combustion Co-design Center to illustrate its value to the exascale co-design process. ExaSAT analysis provides insights into the hardware and software trade-offs and lays the groundwork for exploring a more targeted set of design points using cycle-accurate architectural simulators.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Conceptual aspects of a statistical-dynamical approach to represent landscape subgrid-scale heterogeneities in atmospheric models

A parameterization of land surfaces based on a statistical-dynamical approach is presented. With this approach, the most important characteristics of the soil-plant-atmosphere system that affect the partition of energy (e.g., plant stomatal conductance, soil humidity, and surface roughness) are represented by a pdf rather than by a single 'representative' value. A primary simplified version of this parameterization is used to estimate the land-surface energy fluxes that are produced at the grid scale by various distributions of stomatal conductance under a broad range of environmental conditions. To demonstrate the approach's potential, results are compared with the same fluxes calculated with a big leaf model using the mean stomatal conductance that corresponds to the distributions. Large absolute and relative differences are obtained between the two schemes for many combinations of stomatal conductance pdfs and environmental conditions.

Avissar, Roni↗
DOE OSTI2024

Implementation and Exploration of Parameterizations of Large-Scale Dynamics in NCAR's Single Column Atmosphere Model SCAM6

A single column model with parameterized large-scale (LS) dynamics is used to better understand the response of steady-state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large-scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR's Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies -Weak Temperature Gradient (GASS-WTG) Intercomparison project. Using the same parameter settings and simulation configuration as in the GASS-WTG Intercomparison project, SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the large-scale circulation's damping of tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS-WTG Intercomparison project. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid-scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column's reference radiative-convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively.

54 ENVIRONMENTAL SCIENCES↗