Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeled”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nuclear Energy in Long-Term System Models: A Multi-Model Perspective

Long-term energy system models–including electric sector capacity expansion models–are widely used tools for informing planning, technology assessment, and policy analysis. Recent decarbonization goals and rapid technological change have increased the need to appropriately represent economic characteristics and technical details of energy system resources, including variable renewable energy, energy storage technologies, carbon-capture-equipped capacity, and nuclear energy. Nuclear power represents about 20% of electricity generation and 50% of carbon-free electricity in the United States as of 2021. However, there are many perspectives on the role of existing and new nuclear in the future U.S. energy system, which is reflected in the broad range of potential contributions reported in the literature. This project aims to understand how issues central to nuclear energy are represented in long-term energy models. Building on earlier collaborations that focused on variable renewable energy and energy storage, this project convenes four modeling teams that use national-scale long-term energy system models from the Electric Power Research Institute, the National Renewable Energy Laboratory, the U.S. Energy Information Administration, and the U.S. Environmental Protection Agency to share methods and data, update models, run coordinated scenarios, and identify research needs. Improving tools can provide more insightful analyses and ensure that methods are more transparent.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2020

GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere. Part 2: Validation of Large-Scale Transport and Evaluation of Climate Response

Here we examine the large‐scale transport characteristics of the new ``Middle Atmosphere" NASA Goddard Institute for Space Studies climate model (E2.2). First we evaluate the stratospheric transport circulation in historical atmosphere‐only simulations integrated with interactive trace gas and aerosol chemistry. Compared to lower vertical resolution model versions, E2.2 exhibits improved tropical ascent and older stratospheric mean ages that are more consistent with observed values. In the troposphere, poleward transport to the Arctic and interhemispheric mean ages in E2.2 are comparable to models participating in the Chemistry Climate Modeling Initiative. In addition to validating E2.2 we also assess its ``transport sensitivity" using the coupled atmosphere‐ocean abrupt 4xCO2 and transient 1%CO2 simulations submitted to the Coupled Model Intercomparison Project, Phase 6, along with a 2xCO2 simulation used to evaluate the linearity of the transport circulation's response to increased CO2. We show that decreases (increases) in a stratospheric mean age (idealized surface loss) tracer scale linearly with increased lower stratospheric upwelling, which also increases linearly with warming tropical sea surface temperatures (SSTs). Abrupt 2xCO2 and 4xCO2 experiments constrained with (fixed) pre‐industrial SSTs are also used to quantify the relative importance of rapid adjustments versus SST feedbacks to the transport circulation responses in the model. Finally, sensitivity experiments are presented to illustrate the impact of changes in the convective parameterization on stratospheric transport.

GISS Model E2.2↗
DOE OSTI2018

Coupling methods between finite element-based Boussinesq-type wave and particle-based free-surface flow models [Coupling methods between finite element-based Boussinesq-typewave and particle-based free-surface flow models]

We develop one-way coupling methods between a Boussinesq-type wave model based on the discontinuous Galerkin finite element method and a free-surface flow model based on a mesh-free particle method to strike a balance between accuracy and computational cost. In our proposed model, computation of the wave model in the global domain is conducted first, and the nonconstant velocity profiles in the vertical direction are reproduced by using its results. Computation of the free-surface flow is performed in a local domain included within the global domain with interface boundaries that move along the reproduced velocity field in a Lagrangian fashion. To represent the moving interfaces, we used a polygon wall boundary model for mesh-free particle methods. Verification and validation tests of our proposed model are performed, and results obtained by the model are compared with theoretical values and experimental results to show its accuracy and applicability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

North American extreme precipitation events and related large-scale meteorological patterns: a review of statistical methods, dynamics, modeling, and trends

This paper surveys the current state of knowledge regarding large-scale meteorological patterns (LSMPs) associated with short-duration (less than 1 week) extreme precipitation events over North America. In contrast to teleconnections, which are typically defned based on the characteristic spatial variations of a meteorological feld or on the remote circulation response to a known forcing, LSMPs are defned relative to the occurrence of a specifc phenomenon—here, extreme precipitation—and with an emphasis on the synoptic scales that have a primary infuence in individual events, have medium-range weather predictability, and are well-resolved in both weather and climate models. For the LSMP relationship with extreme precipitation, we consider the previous literature with respect to defnitions and data, dynamical mechanisms, model representation, and climate change trends. There is considerable uncertainty in identifying extremes based on existing observational precipitation data and some limitations in analyzing the associated LSMPs in reanalysis data. Many diferent defnitions of “extreme” are in use, making it difcult to directly compare diferent studies. Dynamically, several types of meteorological systems—extratropical cyclones, tropical cyclones, mesoscale convective systems, and mesohighs—and several mechanisms—fronts, atmospheric rivers, and orographic ascent—have been shown to be important aspects of extreme precipitation LSMPs. The extreme precipitation is often realized through mesoscale processes organized, enhanced, or triggered by the LSMP. Understanding of model representation, trends, and projections for LSMPs is at an early stage, although some promising analysis techniques have been identifed and the LSMP perspective is useful for evaluating the model dynamics associated with extremes.

Barlow, Mathew↗
DOE OSTI2019

Investigation of goodness of model data fit using PLSR and PCR regression models to determine informative wavelength band in NIR region for non-invasive blood glucose prediction

Diabetes mellitus is a metabolic disorder characterized by the deficiency in insulin secretion or insulin action on target tissues. Diabetes being diagnosed, requires frequent monitoring. Current measurement techniques are invasive in nature, which is painful. Non-invasive measurement methods for the prediction of blood glucose is a big challenge. Hence this paper focuses on a comparative study to analyze the performance of the two regression models: partial least square regression (PLSR) and principal component regression (PCR) models to arrive at the best regression model for the prediction of blood glucose non-invasively based on NIR diffuse reflectance spectroscopy. The regression model developed involves spectral data collected from 32 subjects from a diabetic center using diffuse reflectance spectrometer (DRS). The spectroscopic data were acquired from the subjects both invasively using clinical method and non-invasively using DRS setup in the NIR spectral range whose wavelength ranges from 750 to 1040 nm. Multivariate analysis using PLSR is found to be a better regression model in terms of the fitted response and estimated mean square prediction error. It is observed that the mean square error value is 0.04 mg/dl for the tenth component in the PLSR model. The developed regression model also proves to be a useful tool in identifying the informative wavelength bands for blood glucose measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Improving the Modeling of Slag and Steel Bath Chemistry in an Electric Arc Furnace Process Model

An improvement to the thermochemical module of the electric arc furnace (EAF) process model developed by Meier based on the work of Logar, Dovžan, and Škrjanc is presented. Different models for the calculation of activities in melt and slag are implemented, and separate reaction zones are defined for the interaction of slag and melt, the injection of oxygen, and the injection of carbon. For each zone, equilibrium compositions and reaction rates are calculated. Furthermore, diffusion of species is considered as a rate-limiting factor in reactions between slag and melt, and diffusion rates are calculated from the bulk melt to the reaction zone where reactions with the slag take place. Oxygen and sulfur dissolved in the melt and CaS in the slag are added as new species not previously considered in the EAF model. The treatment of carbon is revised to reduce model complexity and improve accuracy. The improved model is validated using extensive data from an industrial EAF, and results are compared to measured data as well as results obtained with the unmodified model. The different models for the determination of thermodynamic activities and their impacts on the duration of the simulation as well as its results are evaluated.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A heuristic model inversion for coupled thermo-hydro-mechanical modelling of triaxial experiments

As multiphysics geomechanical models get developed, their increasing complexity and number of parameters make it particularly difficult to calibrate against experimental data. In this contribution, we present a heuristic workflow to invert for parameters of a coupled Thermo-Hydro-Mechanical (THM) model in a way that helps the theoretical modellers refine their definition of the underlying elasto-visco-plastic model itself. We apply this workflow to the calibration of deviatoric and volumetric data for two sets of triaxial experiments on mudstone and sandstone. Here, we show that beyond the calibration of well-defined parameters of the THM model, one or two scaling factors of the experimentally obtained yield surface and two more factors to capture the confinement and pore pressure dependency of the flow law lead to the satisfactory matching of series of experiments at varying confinements for two different rock types. Using this physical model, we also show that tracking the volumetric component of the mechanical power in the numerical simulations might allow reducing the number of experiments required to calibrate the model.

58 GEOSCIENCES↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An automatic performance model-based scheduling tool for coupled climate system models

The prediction of the climate system is highly depended on the efficient integration of observations and simulations of the Earth. This is regarded as a canonical example of the cyber-physical system. The climate system model, the simulation engine in this cyber-physical system, is one of most challenging applications in scientific computing. It utilizes the multi-physics simulation that couples multiple components, conducts decadal to millennium simulations, and has long been an important application on supercomputers. However, current climate system models suffer from the inefficient task scheduling methods resulting in an intolerable simulation time. Take the Community Earth System Model (CESM), the most widely used climate system model, as an example, one major reason that CESM suffers from bad performances is the huge overhead to rationally distribute processes among the coupled heterogeneous components. According to the report of NCAR, every percent improvement in CESM performance frees up to the equivalent of $250,000 in computing resources in their scientific experiments. To address such challenge, our paper firstly constructs a lightweight and accurate performance model for effectively capturing and predicting the heterogeneous time-to-solution performance of end-to-end CESM components with a given simulation configuration. Then, based on the performance model, we further propose an efficient scheduling strategy based on rectangular packing method to determine the best process layout among different components, and the process numbers assigned to each component. Our evaluations show that we can achieve 58% average run time reductions on CESM comparing to the widely used sequential process layout for a scale of 144 to 480 cores on typical CPU clusters. And we can save 4 million CPU hours when we conduct one standard scientific experiment (a 2870-year simulation), which equals to save $40,089 with a charge of $0.01 per CPU hour. Meanwhile, 26% extra performance improvements could also be gained in our methods comparing to the heuristic branch and bound algorithm with the guidance of the known curve-fitting performance model.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the Forward Modeling of Radar Doppler Spectrum Width From LES: Implications for Model Evaluation

Large-eddy simulations of an observed single-layer Arctic mixed-phase cloud are analyzed to study the value of forward modeling of profiling millimeter wave cloud radar Doppler spectral width for model evaluation. Individual broadening terms and their uncertainties are quantified for the observed spectral width and compared to modeled broadening terms. Modeled turbulent broadening is narrower than the observed values when the turbulent kinetic energy dissipation rate from the subgrid scale model is used in the forward model. The total dissipation rates, estimated with the subgrid scale dissipation rates and the numerical dissipation rates, agree much better with both the retrieved dissipation rates and those inferred from the power spectra of the simulated vertical air velocity. The comparison of the microphysical broadening provides another evaluative measure of the ice properties in the simulation. To accurately retrieve dissipation rates as well as each broadening term from the observations, we suggest a few modifications to previously presented techniques. First, we show that the inertial subrange spectrum filtered with the radar sampling volume is a better underlying model than the unfiltered -5/3 law for the retrieval of the dissipation rate from the power spectra of the mean Doppler velocity. Second, we demonstrate that it is important to filter out turbulence and remove the layer-mean reflectivity-weighted mean fall speed from the observed mean Doppler velocity to avoid overestimation of shear broadening. Finally, we provide a method to quantify the uncertainty in the retrieved dissipation rates, which eventually propagates to the uncertainty in the microphysical broadening.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2023

Towards fast, accurate predictions of RF simulations via data-driven modeling: Forward and lateral models

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗