SEARCH · Search NASA

Results for “Model Coupling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Reduced Complexity Model Intercomparison Project Phase 1: introduction and evaluation of global-mean temperature response

Reduced-complexity climate models (RCMs) are critical in the policy and decision making space, and are directly used within multiple Intergovernmental Panel on Climate Change (IPCC) reports to complement the results of more comprehensive Earth system models. To date, evaluation of RCMs has been limited to a few independent studies. Here we introduce a systematic evaluation of RCMs in the form of the Reduced Complexity Model Intercomparison Project (RCMIP). We expect RCMIP will extend over multiple phases, with Phase 1 being the first. In Phase 1, we focus on the RCMs' global-mean temperature responses, comparing them to observations, exploring the extent to which they emulate more complex models and considering how the relationship between temperature and cumulative emissions of CO2 varies across the RCMs. Our work uses experiments which mirror those found in the Coupled Model Intercomparison Project (CMIP), which focuses on complex Earth system and atmosphere–ocean general circulation models. Using both scenario-based and idealised experiments, we examine RCMs' global-mean temperature response under a range of forcings. We find that the RCMs can all reproduce the approximately 1°C of warming since pre-industrial times, with varying representations of natural variability, volcanic eruptions and aerosols. We also find that RCMs can emulate the global-mean temperature response of CMIP models to within a root-mean-square error of 0.2°C over a range of experiments. Furthermore, we find that, for the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP)-based scenario pairs that share the same IPCC Fifth Assessment Report (AR5)-consistent stratospheric-adjusted radiative forcing, the RCMs indicate higher effective radiative forcings for the SSP-based scenarios and correspondingly higher temperatures when run with the same climate settings. In our idealised setup of RCMs with a climate sensitivity of 3°C, the difference for the ssp585–rcp85 pair by 2100 is around 0.23°C(±0.12°C) due to a difference in effective radiative forcings between the two scenarios. Phase 1 demonstrates the utility of RCMIP's open-source infrastructure, paving the way for further phases of RCMIP to build on the research presented here and deepen our understanding of RCMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Simulations of aerosol pH in China using WRF-Chem (v4.0): sensitivities of aerosol pH and its temporal variations during haze episodes

Aerosol pH is a fundamental property of aerosols in terms of atmospheric chemistry and its impact on air quality, climate, and health. Precise estimation of aerosol pH in chemical transport models (CTMs) is critical for aerosol modeling and thus influences policy development that partially relies on results from model simulations. We report the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem) simulated PM 2.5 pH over China during a period with heavy haze episodes in Beijing, and explore the sensitivity of the modeled aerosol pH to factors including emissions of nonvolatile cations (NVCs) and NH 3 , aerosol phase state assumption, and heterogeneous production of sulfate. We find that default WRF-Chem could predict spatial patterns of PM 2.5 pH over China similar to other CTMs, but with generally lower pH values, largely due to the underestimation of alkaline species (NVCs and NH 3 ) and the difference in thermodynamic treatments between different models. Increasing NH 3 emissions in the model would improve the modeled pH in comparison with offline thermodynamic model calculations of pH constrained by observations. In addition, we find that the aerosol phase state assumption and heterogeneous sulfate production are important in aerosol pH predictions for regions with low relative humidity (RH) and high anthropogenic SO 2 emissions, respectively. These factors should be better constrained in model simulations of aerosol pH in the future. Analysis of the modeled temporal trend of PM 2.5 pH in Beijing over a haze episode reveals a clear decrease in pH from 5.2 ± 0.9 in a clean period to 3.6 ± 0.5 in a heavily polluted period. The increased acidity under more polluted conditions is largely due to the formation and accumulation of secondary species including sulfuric acid and nitric acid, even though being modified by alkaline species (NVCs, NH 3 ). Our result suggests that NO 2 oxidation is unlikely to be important for heterogeneous sulfate production during the Beijing haze as the effective pH for NO 2 oxidation of S(IV) is at a higher pH of ~6.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Xanthos Output Dataset Under ISIMIP3b Selected CMIP6 Scenarios: 1850 - 2100

Overview Xanthos is an open-source hydrologic model, written in Python, designed to quantify and analyze global water availability. Xanthos simulates historical and future global water availability on a monthly time step at a spatial resolution of 0.5 geographic degrees. Xanthos was designed to be extensible and used by scientists that study global water supply and work with the Global Change Analysis Model (GCAM). Details This dataset includes Xanthos outputs forced by the ISIMIP3b bias-adjusted CMIP6 climate forcing (1850 - 2100). There are 30 scenarios consists of 10 Coupled Model Intercomparison Project phase 6 (CMIP6) General Circulation Models (GCM) and 3 Shared Socioeconomic Pathways (SSP). Each folder is named as cmip6_{gcm}_{ensemble-member}_{ssp}. 10 CMIP6 GCMs: GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, UKESM1-0-LL, CanESM5, CNRM-CM6-1, CNRM-ESM2-1, EC-Earth3, and MIROC6 3 SSP-RCP Scenarios: SSP1-2.6, SSP3-7.0, and SSP5-8.5 The table below summarizes the details of output variables within each scenario folder. | Output Variable | File Short Name | Spatial Resolution | Unit | | ----------------------------------------------- | -------------------------------------- | -------------------------------- | --------------------- | | Runoff | q | 0.5x0.5 degree | km3/month | | Runoff | Basin_runoff | GCAM basin | km3/month | | Runoff | Country_runoff | Country | km3/month | | Runoff | GCAMRegion_runoff | GCAM region | km3/month | | Streamflow | avgchflow | 0.5x0.5 degree | m3/second | | Potential Evapotranspiration | pet | 0.5x0.5 degree | km3/month | | Actual Evapotranspiration | aet | 0.5x0.5 degree | km3/month | | Accessible Water | accessible_water | GCAM basin | km3/year | | Hydropower Generation | actual_hydro | GCAM region | EJ/year |

Climate Change↗

DOE OSTI2023

Xanthos Output Dataset Under ISIMIP3b Selected CMIP6 Scenarios: 1850 - 2100

Climate Change↗

DOE OSTI2025

Quantum simulation of charge and exciton transfer in multi-mode models using engineered reservoirs

Quantum simulation enables studies of open-system dynamics in non-perturbative regimes by programming electronic, vibrational, and environmental interactions on comparable energy scales. Trapped ions offer this capability, combining spins, phonons, and tunable dissipation on one platform. We demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Using tailored spin-phonon interactions with reservoir engineering, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor sites and track its non-equilibrium dynamics. We continuously tune the system from the charge transfer regime to the vibrationally assisted exciton transfer regime and find that degenerate modes enhance transfer rates at large energy gaps, while non-degenerate modes activate pathways that reduce the energy-gap dependence. Thus, the presence of one additional vibration introduces interfering pathways and reshapes non-perturbative excitation transfer. Our results establish a scalable, hardware-efficient route to simulate vibronic processes with engineered environments.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

A continuum model for heat and mass transfer in moving-bed reactors for thermochemical energy storage

In this work, a continuum heat and mass transfer model coupling transport phenomena and high-temperature thermochemical reactions is developed for stationary packed-bed and counter-flow moving-bed reactors. After presenting the general modeling framework, we focus on the 2D axisymmetric version of the model for which validation is conducted with experimental results for a packed-bed reactor in the literature for manganese-iron oxide reduction/oxidation and an in-house counter-flow moving-bed reactor for magnesium-manganese oxide reduction up to 1450 °C. Transient simulation results including the local distributions of gas/solid temperatures, oxygen concentration and the extent of reaction, as well as the various energy flow components and energy conversion efficiencies are reported. The results based on the 2D axisymmetric model are also compared with those obtained from a previous 1D model. The comparison shows that capturing the radial variation is critical in reactor modeling and the 2D results demonstrate improved agreement with experiments. Specifically, large temperature variations along the radial direction are observed especially in the reaction zone; this non-uniform radial temperature distribution has a significant effect on the chemical reaction extent due to its strong dependence on temperature; and the overall oxygen concentration at the reactor exit and the predicted system efficiency are slightly lower in the 2D model compared to the 1D model. Finally, the present heat and mass transfer model can provide valuable insights into reactor design, scale-up, and operating conditions selection to maximize system energy storage efficiency.

25 ENERGY STORAGE↗

DOE OSTI2015

Analysis of the S{sub 2}←S{sub 0} vibronic spectrum of the ortho-cyanophenol dimer using a multimode vibronic coupling approach

The S{sub 2}←S{sub 0} vibronic spectrum of the ortho-cyanophenol dimer (oCP){sub 2} is analyzed in a joint experimental and theoretical investigation. Vibronic excitation energies up to 750 cm{sup −1} are covered, which extends our previous analysis of the quenching of the excitonic splitting in this and related species [Kopec et al., J. Chem. Phys. 137, 184312 (2012)]. As we demonstrate, this necessitates an extension of the coupling model. Accordingly, we compute the potential energy surfaces of the ortho-cyanophenol dimer (oCP){sub 2} along all relevant normal modes using the approximate second-order coupled cluster method RI-CC2 and extract the corresponding coupling constants using the linear and quadratic vibronic coupling scheme. These serve as the basis to calculate the vibronic spectrum. The theoretical results are found to be in good agreement with the experimental highly resolved resonant two-photon ionization spectrum. This allows to interpret key features of the excitonic and vibronic interactions in terms of nodal patterns of the underlying vibronic wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Linearized Boltzmann transport model for jet propagation in the quark-gluon plasma: Heavy quark evolution

A linearized Boltzmann transport (LBT) model coupled with hydrodynamical background is established to describe the evolution of jet shower partons and medium excitations in high energy heavy-ion collisions. We extend the LBT model to include both elastic and inelastic processes for light and heavy partons in the quark-gluon plasma. A hybrid model of fragmentation and coalescence is developed for the hadronization of heavy quarks. Within this framework, we investigate how heavy flavor observables depend on various ingredients, such as different energy loss and hadronization mechanisms, the momentum and temperature dependences of the transport coefficients, and the radial flow of the expanding fireball. Our model calculations show good descriptions of the D meson suppression and elliptic flow observed at the Larege Hadron Collider and the Relativistic Heavy-Ion Collider. In addition, the prediction for the Pb-Pb collisions at s NN =5.02 TeV is provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Retrospective Analysis and Bayesian Model Averaging of CMIP6 Precipitation in the Nile River Basin

The Nile River basin is one of the global hotspots vulnerable to climate change impacts because of a fast-growing population and geopolitical tensions. Previous studies demonstrated that general circulation models (GCMs) frequently show disagreement in the sign of change in annual precipitation projections. Here, we first evaluate the performance of 20 GCMs from phase six of the Coupled Model Intercomparison Project (CMIP6) benchmarked against a high-spatial-resolution precipitation dataset dating back to 1983 from Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks–Climate Data Record (PERSIANN-CDR). Next, a Bayesian model averaging (BMA) approach is adopted to derive probability distributions of precipitation projections in the Nile basin. Retrospective analysis reveals that most GCMs exhibit considerable (up to 64% of mean annual precipitation) and spatially heterogenous bias in simulating annual precipitation. Moreover, it is shown that all GCMs underestimate interannual variability; thus, the ensemble range is underdispersive and is a poor indicator of uncertainty. The projected changes from the BMA model show that the value and sign of change vary considerably across the Nile basin. Specifically, it is found that projected changes in the two headwaters basins, namely, the Blue Nile and Upper White Nile, are 0.03% and -1.65%, respectively; both are statistically insignificant at α = 0.05. The uncertainty range estimated from the BMA model shows that the probability of a precipitation decrease is much higher in the Upper White Nile basin whereas projected change in the Blue Nile is highly uncertain both in magnitude and sign of change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

The E3SM Diagnostics Package (E3SM Diags v2.7): a Python-based diagnostics package for Earth system model evaluation

Abstract. The E3SM Diagnostics Package (E3SM Diags) is a modern, Python-based Earth system model (ESM) evaluation tool (with Python module name e3sm_diags), developed to support the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM). E3SM Diags provides a wide suite of tools for evaluating native E3SM output, as well as ESM data on regular latitude–longitude grids, including output from Coupled Model Intercomparison Project (CMIP) class models. E3SM Diags is modeled after the National Center for Atmospheric Research (NCAR) Atmosphere Model Working Group (AMWG, 2022) diagnostics package. In its version 1 release, E3SM Diags included a set of core essential diagnostics to evaluate the mean physical climate from model simulations. As of version 2.7, more process-oriented and phenomenon-based evaluation diagnostics have been implemented, such as analysis of the quasi-biennial oscillation (QBO), the El Niño–Southern Oscillation (ENSO), streamflow, the diurnal cycle of precipitation, tropical cyclones, ozone and aerosol properties. An in situ dataset from DOE's Atmospheric Radiation Measurement (ARM) program has been integrated into the package for evaluating the representation of simulated cloud and precipitation processes. This tool is designed with enough flexibility to allow for the addition of new observational datasets and new diagnostic algorithms. Additional features include customizable figures; streamlined installation, configuration and execution; and multiprocessing for fast computation. The package uses an up-to-date observational data repository maintained by its developers, where recent datasets are added to the repository as they become available. Finally, several applications for the E3SM Diags module were introduced to fit a diverse set of use cases from the scientific community.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Evaluating precipitation distributions at regional scales: a benchmarking framework and application to CMIP5 and 6 models

Abstract. As the resolution of global Earth system models increases, regional-scale evaluations are becoming ever more important. This study presents a framework for quantifying precipitation distributions at regional scales and applies it to evaluate Coupled Model Intercomparison Project (CMIP) 5 and 6 models. We employ the Intergovernmental Panel on Climate Change (IPCC) sixth assessment report (AR6) climate reference regions over land and propose refinements to the oceanic regions based on the homogeneity of precipitation distribution characteristics. The homogeneous regions are identified as heavy-, moderate-, and light-precipitating areas by K-means clustering of Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (GPM) version 6 final run product (IMERG) precipitation frequency and amount distributions. With the global domain partitioned into 62 regions, including 46 land and 16 ocean regions, we apply 10 established precipitation distribution metrics. The collection includes metrics focused on the maximum peak, lower 10th percentile, and upper 90th percentile in precipitation amount and frequency distributions; the similarity between observed and modeled frequency distributions; an unevenness measure based on cumulative amount; average total intensity on all days with precipitation; and number of precipitating days each year. We apply our framework to 25 CMIP5 and 41 CMIP6 models, as well as six observation-based products of daily precipitation. Our results indicate that many CMIP5 and 6 models substantially overestimate the observed light-precipitation amount and frequency, as well as the number of precipitating days, especially over midlatitude regions outside of some land regions in the Americas and Eurasia. Improvement from CMIP5 to 6 is shown in some regions, especially in midlatitude regions, but it is not evident globally, and over the tropics most metrics point toward degradation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

A new WRF-Chem treatment for studying regional-scale impacts of cloud processes on aerosol and trace gases in parameterized cumuli

A new treatment of cloud effects on aerosol and trace gases within parameterized shallow and deep convection, and aerosol effects on cloud droplet number, has been implemented in the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) version 3.2.1 that can be used to better understand the aerosol life cycle over regional to synoptic scales. The modifications to the model include treatment of the cloud droplet number mixing ratio; key cloud microphysical and macrophysical parameters (including the updraft fractional area, updraft and downdraft mass fluxes, and entrainment) averaged over the population of shallow clouds, or a single deep convective cloud; and vertical transport, activation/resuspension, aqueous chemistry, and wet removal of aerosol and trace gases in warm clouds. These changes have been implemented in both the WRF-Chem chemistry packages as well as the Kain–Fritsch (KF) cumulus parameterization that has been modified to better represent shallow convective clouds. Testing of the modified WRF-Chem has been completed using observations from the Cumulus Humilis Aerosol Processing Study (CHAPS). The simulation results are used to investigate the impact of cloud–aerosol interactions on regional-scale transport of black carbon (BC), organic aerosol (OA), and sulfate aerosol. Based on the simulations presented here, changes in the column-integrated BC can be as large as –50% when cloud–aerosol interactions are considered (due largely to wet removal), or as large as +40% for sulfate under non-precipitating conditions due to sulfate production in the parameterized clouds. The modifications to WRF-Chem are found to account for changes in the cloud droplet number concentration (CDNC) and changes in the chemical composition of cloud droplet residuals in a way that is consistent with observations collected during CHAPS. Efforts are currently underway to port the changes described here to the latest version of WRF-Chem, and it is anticipated that they will be included in a future public release of WRF-Chem.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Equation of state for neutron stars with hyperons and quarks in the relativistic Hartree–Fock approximation

We construct the equation of state (EoS) for neutron stars explicitly including hyperons and quarks. Using the quark–meson coupling model with the relativistic Hartree–Fock approximation, the EoS for hadronic matter is derived by taking into account the strange (σ* and ϕ) mesons as well as the light non-strange (σ, ω, π, and ρ) mesons. Relevant coupling constants are determined to reproduce the experimental data of nuclear matter and hypernuclei in SU(3) flavor symmetry. For quark matter, we employ the MIT bag model with a one-gluon-exchange interaction, and Gibbs criteria for chemical equilibrium in the phase transition from hadrons to quarks. We find that the strange vector (ϕ) meson and the Fock contribution make the hadronic EoS stiff, and that the maximum mass of a neutron star can be consistent with the observed mass of heavy neutron stars even if the coexistence of hadrons and quarks takes place in the core. However, in the present calculation, the transition to pure quark matter does not occur in stable neutron stars. Furthermore, the lower bound of the critical chemical potential of the quark–hadron transition at zero temperature turns out to be around 1.5 GeV in order to be consistent with the recent observed neutron-star data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A time-dependent chloride diffusion model for predicting initial corrosion time of reinforced concrete with slag addition

The effect of granulated blast furnace slag (GBFS) addition on the threshold chloride concentration (TCC) for rebar corrosion was investigated. A modified diffusion model, coupled with a time-dependent effective diffusion coefficient and surface chloride concentration, was proposed to predict the chloride profile. The corrosion initiation time was estimated based on the model predictions and the measured TCC. The results indicate that adding GBFS decreases the TCC by lowering the pH value of the pore solution. The evolution of corrosion potential and current density is found to obey a 3-parameter Weibull distribution. MIP tests show that adding GBFS contributes to refinement of pore structure by decreasing the fraction of large capillary pores. The time-dependent model exhibits good predictive strength and helps understand how GBFS addition delays the corrosion initiation by retarding the chloride diffusion, though a lower TCC is obtained.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Characterizing and quantifying uncertainty in projections of climate change impacts on air quality

Climate change can aggravate air pollution, with important public health and environmental consequences. While major sources of uncertainty in climate change projections - greenhouse gas (GHG) emissions scenario, model response, and internal variability - have been investigated extensively, their propagation to estimates of air quality impacts has not been systematically assessed. Here, we compare these uncertainties using a coupled modeling framework that includes a human activity model, an Earth system model of intermediate complexity, and a global atmospheric chemistry model. Uncertainties in projections of U.S. air quality under 21st century climate change are quantified based on a climate-chemistry ensemble that includes multiple initializations, representations of climate sensitivity, and climate policy scenarios, under constant air pollution emissions. We find that climate-related uncertainties are comparable at mid-century, making it difficult to distinguish the impact of variations in GHG emissions on ozone and particulate matter pollution. While GHG emissions scenario eventually becomes the dominant uncertainty based on the scenarios considered, all sources of uncertainty are significant through the end of the century. The results provide insights into intrinsically different uncertainties in projections of air pollution impacts and the potential for large ensembles to better capture them.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Humidity Bias and Effect on Simulated Aerosol Optical Properties during the Ganges Valley Experiment

The radiosonde humidity profiles available during the Ganges Valley Experiment were compared to those simulated from the regional Weather Research and Forecasting (WRF) model coupled with a chemistry module (WRF -Chern) and the global reanalysis datasets. Large biases were revealed. On a monthly mean basis at Nainital, located in northern India, the WRFChern model simulates a large moist bias in the free troposphere (up to +20%) as well as a large dry bias in the boundary layer (up to -30%). While the overall pattern of the biases is similar, the magnitude of the biases varies from time to time and from one location to another. At Thiruvananthapuram, the magnitude of the dry bias is smaller, and in contrast to Nainital, the higher-resolution regional WRF -Chern model generates larger moist biases in the upper troposphere than the global reanalysis data. Furthermore, the humidity biases in the upper troposphere, while significant, have little impact on the model estimation of column aerosol optical depth (AOD). The frequent occurrences of the dry boundary-layer bias simulated by the large-scale models tend to lead to the underestimation of AOD. It is thus important to quantify the humidity vertical profiles for aerosol simulations over South Asia.

aerosol extinction↗

DOE OSTI2026

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

DOE OSTI2016

Comparison of fluid neutral models for one-dimensional plasma edge modeling with a finite volume solution of the Boltzmann equation

We derive fluid neutral approximations for a simplified 1D edge plasma model, suitable to study the neutral behavior close to the target of a nuclear fusion divertor, and compare its solutions to the solution of the corresponding kinetic Boltzmann equation. The plasma is considered as a fixed background extracted from a detached 2D simulation. We show that the Maxwellian equilibrium distribution is already obtained very close to the target, justifying the use of a fluid approximation. We compare three fluid neutral models: (i) a diffusion model; (ii) a pressure-diffusion model (i.e., a combination of a continuity and momentum equation) assuming equal neutral and ion temperatures; and (iii) the pressure-diffusion model coupled to a neutral energy equation taking into account temperature differences between neutrals and ions. Partial reflection of neutrals reaching the boundaries is included in both the kinetic and fluid models. We propose two methods to obtain an incident neutral flux boundary condition for the fluid models: one based on a diffusion approximation and the other assuming a truncated Chapman-Enskog distribution. The pressure-diffusion model predicts the plasma sources very well. The diffusion boundary condition gives slightly better results overall. Although including an energy equation still improves the results, the assumption of equal ion and neutral temperature already gives a very good approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗