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 361 records · Page 20

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↗

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↗

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↗

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↗

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.↗

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↗

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↗

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↗

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↗

CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences

We report the development and deployment of a coarse-graining method that is well suited for computer simulations of aggregation and phase separation of protein sequences with block-copolymeric architectures. Our algorithm, named CAMELOT for Coarse-grained simulations Aided by MachinE Learning Optimization and Training, leverages information from converged all atom simulations that is used to determine a suitable resolution and parameterize the coarse-grained model. To parameterize a system-specific coarse-grained model, we use a combination of Boltzmann inversion, non-linear regression, and a Gaussian process Bayesian optimization approach. The accuracy of the coarse-grained model is demonstrated through direct comparisons to results from all atom simulations. We demonstrate the utility of our coarse-graining approach using the block-copolymeric sequence from the exon 1 encoded sequence of the huntingtin protein. This sequence comprises of 17 residues from the N-terminal end of huntingtin (N17) followed by a polyglutamine (polyQ) tract. Simulations based on the CAMELOT approach are used to show that the adsorption and unfolding of the wild type N17 and its sequence variants on the surface of polyQ tracts engender a patchy colloid like architecture that promotes the formation of linear aggregates. These results provide a plausible explanation for experimental observations, which show that N17 accelerates the formation of linear aggregates in block-copolymeric N17-polyQ sequences. The CAMELOT approach is versatile and is generalizable for simulating the aggregation and phase behavior of a range of block-copolymeric protein sequences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

Stochastic Parameterization of Subgrid-Scale Velocity Enhancement of Sea Surface Fluxes

Abstract Subgrid-scale (SGS) velocity variations result in gridscale sea surface flux enhancements that must be parameterized in weather and climate models. Traditional parameterizations are deterministic in that they assign a unique value of the SGS velocity flux enhancement to any given configuration of the resolved state. In this study, we assess the statistics of SGS velocity flux enhancement over a range of averaging scales (as a proxy for varying model resolution) through systematic coarse-graining of a convection-permitting atmospheric model simulation over the Indian Ocean and west Pacific warm pool. Conditioning the statistics of the SGS velocity flux enhancement on 1) the fluxes associated with the resolved winds and 2) the precipitation rate, we find that the lack of a separation between “resolved” and “unresolved” scales results in a distribution of flux enhancements for each configuration of the resolved state. That is, the SGS velocity flux enhancement should be represented stochastically rather than deterministically. The spatial and temporal statistics of the SGS velocity flux enhancement are investigated by using basic descriptive statistics and through a fit to an anisotropic space–time covariance structure. Potential spatial inhomogeneities of the statistics of the SGS velocity flux enhancement are investigated through regional analysis, although because of the relatively short duration of the simulation (9 days) distinguishing true inhomogeneity from sampling variability is difficult. Perspectives for the implementation of such a stochastic parameterization in weather and climate models are discussed.

54 ENVIRONMENTAL SCIENCES↗

Joint Bayesian Inference for Near-surface Explosion Yield

Geophysical observations of an explosion scale with the yield of the explosion. Guided by the principal of yield-scaling and observed relationships in the data, we defined and parameterized seismoacoustic forward models in Ford et al. [2014]. This work extends the model with additional data features of maximum velocity and positive phase duration for the seismic and acoustic data types, respectively, and a new optic data type. Parameters for the seismoacoustic models are estimated via Bayesian inference. The parameterized models can then be used to infer the yield of a new explosion.

58 GEOSCIENCES↗

Target Soil Impact Verification: Experimental Testing and Kayenta Constitutive Modeling.

This report details experimental testing and constitutive modeling of sandy soil deformation under quasi - static conditions. This is driven by the need to understand constitutive response of soil to target/component behavior upon impact . An experimental and constitutive modeling program was followed to determine elastic - plastic properties and a compressional failure envelope of dry soil . One hydrostatic, one unconfined compressive stress (UCS), nine axisymmetric compression (ACS) , and one uniaxial strain (US) test were conducted at room temperature . Elastic moduli, assuming isotropy, are determined from unload/reload loops and final unloading for all tests pre - failure and increase monotonically with mean stress. Very little modulus degradation was discernable from elastic results even when exposed to mean stresses above 200 MPa . The failure envelope and initial yield surface were determined from peak stresses and observed onset of plastic yielding from all test results. Soil elasto - plastic behavior is described using the Brannon et al. (2009) Kayenta constitutive model. As a validation exercise, the ACS - parameterized Kayenta model is used to predict response of the soil material under uniaxial strain loading. The resulting parameterized and validated Kayenta model is of high quality and suitable for modeling sandy soil deformation under a range of conditions, including that for impact prediction.

54 ENVIRONMENTAL SCIENCES↗

Time-series machine-learning error models for approximate solutions to parameterized dynamical systems

This work proposes a machine-learning framework for modeling the error incurred by approximate solutions to parameterized dynamical systems. In particular, we extend the machine-learning error models (MLEM) framework proposed in Ref. Freno and Carlberg (2019) to dynamical systems. The proposed Time-Series Machine-Learning Error Modeling (T-MLEM) method constructs a regression model that maps features – which comprise error indicators that are derived from standard a posteriori error-quantification techniques – to a random variable for the approximate-solution error at each time instance. The proposed framework considers a wide range of candidate features, regression methods, and additive noise models. We consider primarily recursive regression techniques developed for time-series modeling, including both classical time-series models (e.g., autoregressive models) and recurrent neural networks (RNNs), but also analyze standard non-recursive regression techniques (e.g., feed-forward neural networks) for comparative purposes. Finally, numerical experiments conducted on multiple benchmark problems illustrate that the long short-term memory (LSTM) neural network, which is a type of RNN, outperforms other methods and yields substantial improvements in error predictions over traditional approaches.

42 ENGINEERING↗

Atomistically resolved hot exciton relaxation dynamics in CdSe quantum dots: Experiment and theory

Semiconductor quantum dots (QDs) are well known to give rise to a quantum confined structure of excitons. Because of this quantum confinement, new physics of hot exciton relaxation dynamics arises. Decades of work using transient absorption (TA) spectroscopy have yielded initial simple observations, such as estimates of the cooling rate from single pump photon energy experiments. More detailed TA experiments employed variable pump photon energies to measure excitonic state-resolved transition rates. These TA measurements, usually the simplest form, have been employed to characterize QDs and their relaxation dynamics to this day. Yet, these TA measurements are fundamentally lacking in their ability to measure energy-resolved hot exciton cooling, which requires observation of the full cooling history through the real excitonic manifold. Here, we employ coherent multi-dimensional spectroscopy (CMDS) to perform an atomistically directed study of hot exciton cooling in CdSe QDs, revealing energy resolved relaxation dynamics. CMDS experiments are compared with simulations and prior TA measurements and simpler theories. Our findings reveal a hot exciton relaxation dynamics landscape. This relaxation dynamics landscape is a linear or sub-linear function of excess energy for different structures of QDs, with a strong size dependence. Our model simulations parameterized by the empirical pseudopotential model reproduces the experimental functional form and the dependence upon QD diameter and shell.

Atomistic simulations↗

Testing to Evaluate Processes Expected to Occur during MSR Salt Spill Accidents

Obtaining a license for a new nuclear reactor requires the identification and assessment of the potential consequences of specified accident scenarios, which are achieved using accident progression models. Those models need to be parameterized and validated using experimental data, but existing experimental data addressing processes relevant to molten salt reactor accidents are sparse. Specifically, experimental data that quantify the sensitives of important processes to the initial conditions of the spill, the ambient environment, and the containment features are needed to parameterize individual process models. Integrated experiments that simulate accident scenarios are also needed to provide data for model validation, but these experiments will require the use of proven methods to quantify the processes under evaluation. The overarching objectives of this work are to develop the methods for simulating the targeted processes, to determine the effectiveness of the methods in producing the data required for model development, to generate data that can be used to parameterize individual process models, and to provide key insights into the behavior of spilled molten salt that should be considered in models. Experimental methods were designed to quantify aspects of individual processes expected to occur during or after a molten salt spill accident that will affect the fate of spilled molten salt and the radionuclides within. These processes include 1) molten salt spreading and heat transfer, 2) molten salt flowing and freezing in tubing, 3) stainless steel corrosion kinetics in molten salt, and 4) molten salt splashing and aerosol generation. The initial tests described in this report were conducted using eutectic FLiNaK to demonstrate the test methods, the data that are generated, and the analyses of the data to derive values needed for modeling. The primary variables that were tested include initial salt temperature and the presence of volatile surrogate fission products (e.g., cesium and iodine). The developed methods are shown to be effective in quantifying the desired processes and can be applied to study more complex salt compositions of interest to molten salt reactor developers, a wide range of environmental conditions of interest to modelers, and additional variables relevant to salt spill accidents. The developed methods and insights gained from laboratory tests can also be incorporated in future large-scale integrated tests used to simulate molten salt spill accidents.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗