Search NASA⌕ Search

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.

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

At least 1,099 records · Page 61

Tropical Pacific response to continental ice sheet topography

The last glacial maximum was marked by maximum land ice extent and lowest greenhouse gases concentration during the last ice age. We investigate the effect of glacial continental ice sheet topography on the large-scale tropical ocean–atmosphere climate, in particular the tropical Pacific, in an intermediate complexity coupled model. Increasing the thickness of continental ice sheets causes a southward displaced Pacific Intertropical Convergence Zone (ITCZ) and a strengthening (weakening) of northern (southern) hemisphere winter Hadley cell. The equatorial zonal sea surface temperature gradient weakened with an increased continental ice sheets thickness, the reduction being caused by cooling in the western equatorial Pacific and warming in the eastern equatorial Pacific. The evolution of the tropical climate with changing ice thickness has distinct quasi-linear and nonlinear parts. While the linear part is a direct response to the ice topographic changes, the nonlinear part was a finding of the tropical thermocline adjustment. Our analysis of a fully-coupled transient deglacial simulation strongly indicates the dominant role of ice sheet topography in determining the deglacial evolution of the simulated Pacific climate. The thickness of continental ice sheet, separate from ice albedo effect, has significant impact on the tropical ocean–atmosphere climate in particular with the meridional displacement in the Pacific ITCZ. The altered circulation states seen in the model may aid understanding of the relationship between tropical and high-latitude climate records in glacial-interglacial cycles.

54 ENVIRONMENTAL SCIENCES↗

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Steady-state radiochemical transport model of the molten salt reactor experiment

A radiochemical transport model of the Molten Salt Reactor Experiment (MSRE) was formulated based on the material derivative featuring system level 1-D flow providing spatial dependence accounting for multiple phenomena including radioactive decay, chemical removal, and fuel depletion across custom volume regions. Further, the steady-state version, which is based on a system of ODE's, can be solved efficiently and is useful for sensitivity studies utilizing Monte Carlo sampling. Accounting for multiple phenomena in a fully coupled model is necessary to understand the complexities of transport and distribution of radionuclides throughout the primary loop of a flowing fuel reactor.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Impacts of Suppressing Excessive Light Rain on Aerosol Radiative Effects and Health Risks

Global climate models (GCMs) have been used widely to study radiative forcing and health risks of aerosols. A recent study using two GCMs found that light rain plays a dominant role in controlling aerosol loading. However, “too much light rain and too little heavy rain” is a longstanding bias in GCMs. It is unclear how much light rain affects aerosol-cloud-radiation interactions and health risks from air pollution. Here, in this work, we show that, with the correction of the rainfall intensity spectrum in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 by introducing a stochastic deep convection scheme, the reduced frequency of light rain (1–20 mm d -1 ) results in changes of aerosol direct radiative effects (DRE) of up to -0.5 ± 0.03 W/m 2 and aerosol cloud radiative effects (CRE) of up to -0.9 ± 0.03 W/m 2 . The total (CRE + DRE) radiative effects of light rain-mediated aerosol changes exceed the present-day anthropogenic forcing of aerosols relative to preindustrial levels from the Coupled Model Intercomparison Project (CMIP5&6) models. However, the correction of the rainfall intensity spectrum has little effect on anthropogenic aerosol forcing (defined as the radiative perturbation due to changes in aerosol concentrations between the industrial era and preindustrial levels). Due to increased exposure to fine particulates (PM 2.5 ), the estimated global total premature mortality is much higher than previously estimated, by 300,000 ± 60,000 deaths per year, and is more severe in populous regions such as India and China. The findings in this study highlight the need to understand uncertainties in radiative effects and health risks of aerosols due to simulation biases of precipitation in GCMs.

54 ENVIRONMENTAL SCIENCES↗

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗

Impacts of air pollutants from fire and non-fire emissions on the regional air quality in Southeast Asia

Abstract. Severe haze events in Southeast Asia caused by particulate pollution have become more intense and frequent in recent years. Widespread biomass burning occurrences and particulate pollutants from human activities other than biomass burning play important roles in degrading air quality in Southeast Asia. In this study, numerical simulations have been conducted using the Weather Research and Forecasting (WRF) model coupled with a chemistry component (WRF-Chem) to quantitatively examine the contributions of aerosols emitted from fire (i.e., biomass burning) versus non-fire (including fossil fuel combustion, and road dust, etc.) sources to the degradation of air quality and visibility over Southeast Asia. These simulations cover a time period from 2002 to 2008 and are driven by emissions from (a) fossil fuel burning only, (b) biomass burning only, and (c) both fossil fuel and biomass burning. The model results reveal that 39 % of observed low-visibility days (LVDs) can be explained by either fossil fuel burning or biomass burning emissions alone, a further 20 % by fossil fuel burning alone, a further 8 % by biomass burning alone, and a further 5 % by a combination of fossil fuel burning and biomass burning. Analysis of an 24 h PM 2.5 air quality index (AQI) indicates that the case with coexisting fire and non-fire PM 2.5 can substantially increase the chance of AQI being in the moderate or unhealthy pollution level from 23 to 34 %. The premature mortality in major Southeast Asian cities due to degradation of air quality by particulate pollutants is estimated to increase from ∼ 4110 per year in 2002 to ∼ 6540 per year in 2008. In addition, we demonstrate the importance of certain missing non-fire anthropogenic aerosol sources including anthropogenic fugitive and industrial dusts in causing urban air quality degradation. An experiment of using machine learning algorithms to forecast the occurrence of haze events in Singapore is also explored in this study. All of these results suggest that besides minimizing biomass burning activities, an effective air pollution mitigation policy for Southeast Asia needs to consider controlling emissions from non-fire anthropogenic sources.

54 ENVIRONMENTAL SCIENCES↗

Changes in the dependence between global precipitation and temperature from observations and model simulations

Precipitation and temperature are physically related to each other with significant dependences that vary with seasons at regional and global scales. The dependence between precipitation and temperature (P–T dependence) plays a central role in characterizing the joint behaviour of the two variables. Along with extensive studies on the variation in the mean, variance, or extremes of precipitation and temperature under global warming, it is of particular interest to understand the climate–change impacts on the covariability of precipitation and temperature. The aim of this study therefore is to assess climate change impacts on the P–T dependence over global land areas from observations, and the performance of the current–generation Coupled Model Intercomparison Project (CMIP5) models in simulating these patterns. In general, the CMIP5 models correctly simulate the overall pattern of the P–T dependence (around 87–94% of global land areas). Furthermore, they capture both the sign and change of P–T dependence for relatively small regions (around 41–47% of global land areas). Results from this study may provide useful insights for the future development of climate models.

54 ENVIRONMENTAL SCIENCES↗

Impacts of anthropogenic aerosols on fog in North China Plain

Fog poses a sever environmental problem in the North China Plain (NCP), China, which has been witnessing increases anthropogenic emission since the early 1980s. Here, this work first uses the WRF/Chem model coupled with the local anthropogenic emissions to simulate and evaluate a severe fog event occurring in NCP. Comparison of the simulations against observations shows that WRF/Chem well reproduces the general features of temporal evolution of PM 2.5 mass concentration, fog spatial distribution, visibility, and vertical profiles of temperature, water vapor content and relative humidity in the planetary boundary layer throughout the whole period of the fog event. Sensitivity studies are then performed with five different levels of anthropogenic emission as model inputs to systematically examine the comprehensive impacts of aerosols on fog microphysical, macrophysical, radiative, and dynamical properties. The results show that as aerosol concentration increases, fog droplet number concentration and liquid water content all increase nonlinearly; but effective radius decreases. Macrophysical properties (fog fraction, fog duration, fog height, and liquid water path) also increase nonlinearly with increasing aerosol concentration, with rates of changes smaller than microphysical properties. Further analysis reveals distinct aerosol effects on thermodynamic and dynamical conditions during different stages of fog evolution: increasing aerosols invigorate fog formation and development by enhancing longwave-induced instability, fog droplet condensation accompanying latent heat release, and thus turbulence, but delay fog dissipation by reducing surface solar radiation, surface sensible and latent heat fluxes, and thus suppressing turbulence during the dissipation stage.

54 ENVIRONMENTAL SCIENCES↗

Decoherence and energy relaxation in the quantum-classical dynamics for charge transport in organic semiconducting crystals: An instantaneous decoherence correction approach

We explore an instantaneous decoherence correction (IDC) approach for the decoherence and energy relaxation in the quantum-classical dynamics of charge transport in organic semiconducting crystals. These effects, originating from environmental fluctuations, are essential ingredients of the carrier dynamics. The IDC is carried out by measurement-like operations in the adiabatic representation. While decoherence is inherent in the IDC, energy relaxation is taken into account by considering the detailed balance through the introduction of energy-dependent reweighing factors, which could be either Boltzmann (IDC-BM) or Miller-Abrahams (IDC-MA) type. For a non-diagonal electron-phonon coupling model, it is shown that IDC tends to enhance diffusion while energy relaxation weakens this enhancement. As expected, both the IDC-BM and IDC-MA achieve a near-equilibrium distribution at finite temperatures in the diffusion process, while in the Ehrenfest dynamics the electronic system tends to infinite temperature limit. The resulting energy relaxation times with the two kinds of factors lie in different regimes and exhibit different dependences on temperature, decoherence time, and electron-phonon coupling strength, due to different dominant relaxation processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasi-stationary Planetary Scale Waves in Modern Climate Models

We evaluate climate model simulations of geopotential height at the largest space-time scales, considering the first Fourier components (in longitude) of long-term climatological means. Quasi-stationary planetary scale waves emerge from this analysis. Variations of these waves are important components of climate change, for example associated with regional drought. Systematic examination of their long-term climatology provides context for understanding their evolution in time. We compare four reanalysis datasets with “historical runs” from the latest version of the Coupled Model Intercomparison Project (CMIP) and find general agreement, but we also identify model errors outside the range of both observational uncertainty and the uncertainty implied by interannual variability. Errors in wave phase have possible implications for climate model projections of regional climate and of shifting weather patterns. One model in the set we examine, originally created nearly 20 years before the others, exhibits noticeably greater errors.

54 ENVIRONMENTAL SCIENCES↗

The Southern Ocean as a constraint to reduce uncertainty in future ocean carbon sinks

Earth system model (ESM) simulations exhibit large biases compares to observation-based estimates of the present ocean CO 2 sink. The inter-model spread in projections increases nearly 2-fold by the end of the 21st century and therefore contributes significantly to the uncertainty of future climate projections. In this study, the Southern Ocean (SO) is shown to be one of the hot-spot regions for future uptake of anthropogenic CO 2 , characterized by both the solubility pump and biologically mediated carbon drawdown in the spring and summer. Here, we show, by analyzing a suite of fully interactive ESMs simulations from the Coupled Model Intercomparison Project phase 5 (CMIP5) over the 21st century under the high-CO 2 Representative Concentration Pathway (RCP) 8.5 scenario, that the SO is the only region where the atmospheric CO 2 uptake rate continues to increase toward the end of the 21st century. Furthermore, our study discovers a strong inter-model link between the contemporary CO 2 uptake in the Southern Ocean and the projected global cumulated uptake over the 21st century. This strong correlation suggests that models with low (high) carbon uptake rate in the contemporary SO tend to simulate low (high) uptake rate in the future. None the less, our analysis also shows that none of the models fully capture the observed biophysical mechanisms governing the CO 2 fluxes in the SO. The inter-model spread for the contemporary CO 2 uptake in the Southern Ocean is attributed to the variations in the simulated seasonal cycle of surface p CO 2 . Two groups of model behavior have been identified. The first one simulates anomalously strong SO carbon uptake, generally due to both too strong a net primary production and too low a surface p CO 2 in December–January. The second group simulates an opposite CO 2 flux seasonal phase, which is driven mainly by the bias in the sea surface temperature variability. Furthermore, we show that these biases are persistent throughout the 21st century, which highlights the urgent need for a sustained and comprehensive biogeochemical monitoring system in the Southern Ocean to better constrain key processes represented in current model systems.

58 GEOSCIENCES↗

Isolating the effect of biomass burning aerosol emissions on 20th century hydroclimate in South America and Southeast Asia

Abstract Biomass burning is a significant source of aerosol emissions in some regions and has a considerable impact on regional climate. Earth system model simulations indicate that increased biomass burning aerosol emissions contributed to statistically significant decreases in tropical precipitation over the 20th century. In this study, we use the Community Earth System Model version 1 Large Ensemble (CESM1-LENS) experiment to evaluate the mechanisms by which biomass burning aerosol contributed to decreased tropical precipitation, with a focus on South America and Southeast Asia. We analyze the all-but-one forcing simulations in which biomass burning aerosol emissions are held constant while other forcings (e.g., greenhouse gas concentrations) vary throughout the 20th century. This allows us to isolate the influence of biomass burning aerosol on processes that contribute to decreasing precipitation, including cloud microphysics, the radiative effects of absorbing aerosol particles, and alterations in regional circulation. We also show that the 20th century reduction in precipitation identified in the CESM1-LENS historical and biomass burning experiments is consistent across Coupled Model Intercomparison Project Phase 5 models with interactive aerosol schemes and the CESM2 single-forcing experiment. Our results demonstrate that higher concentrations of biomass burning aerosol increases the quantity of cloud condensation nuclei and cloud droplets, limiting cloud droplet size and precipitation formation. Additionally, absorbing aerosols (e.g., black carbon) contribute to a warmer cloud layer, which promotes cloud evaporation, increases atmospheric stability, and alters regional circulation patterns. Corresponding convectively coupled circulation responses, particularly over the tropical Andes, contribute to further reducing the flow of moisture and moisture convergence over tropical land. These results elucidate the processes that affect the water cycle in regions prone to biomass burning and inform our understanding of how future changes in aerosol emissions may impact tropical precipitation over the 21st century.

Environmental Sciences & Ecology↗

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗

Temporal and Spatial Variation in Peatland Carbon Cycling and Implications for Interpreting Responses of an Ecosystem-Scale Warming Experiment

Here, we are conducting a large-scale, long-term climate change response experiment in an ombrotrophic peat bog in Minnesota to evaluate the effects of warming and elevated CO 2 on ecosystem processes using empirical and modeling approaches. To better frame future assessments of peatland responses to climate change, we characterized and compared spatial vs. temporal variation in measured C cycle processes and their environmental drivers. We also conducted a sensitivity analysis of a peatland C model to identify how variation in ecosystem parameters contributes to model prediction uncertainty. High spatial variability in C cycle processes resulted in the inability to determine if the bog was a C source or sink, as the 95% confidence interval ranged from a source of 50 g C m –2 yr –1 to a sink of 67 g C m –2 yr –1 . Model sensitivity analysis also identified that spatial variation in tree and shrub photosynthesis, allocation characteristics, and maintenance respiration all contributed to large variations in the pretreatment estimates of net C balance. Variation in ecosystem processes can be more thoroughly characterized if more measurements are collected for parameters that are highly variable over space and time, and especially if those measurements encompass environmental gradients that may be driving the spatial and temporal variation (e.g., hummock vs. hollow microtopographies, and wet vs. dry years). Together, the coupled modeling and empirical approaches indicate that variability in C cycle processes and their drivers must be taken into account when interpreting the significance of experimental warming and elevated CO 2 treatments.

54 ENVIRONMENTAL SCIENCES↗

Impact of I/O and Data Management in Ensemble Large Scale Climate Forecasting Using EC-Earth3

The EC-Earth climate model is a seamless Earth System Model (ESM) used to carry out climate research in 24 academic institutions and meteorological services from 11 countries in Europe. This model couples several components and it is continuously under development. In this work we present a study regarding the impact of the I/O and data management when using EC-Earth in well-known supercomputing environments. Most large-scale and long-term climate simulators have been developed bearing in mind the paramount importance of its scalability. However, the computational capabilities of the High Performance Computing (HPC) environments increase at so great speed that it is almost impossible to re-implement the whole models so that they are able to exploit efficiently the new features. Therefore, it is necessary to design different strategies to take advantage of them. In this work we present an operational framework to run ensemble simulations in HPC platforms. A set of experiments are presented in order to validate the suitability of this technique. Moreover, the derived impact regarding the I/O and data management aspects is analyzed.

97 MATHEMATICS AND COMPUTING↗

Predicting Near-Term Changes in the Earth System: A Large Ensemble of Initialized Decadal Prediction Simulations Using the Community Earth System Model

The objective of near-term climate prediction is to improve our fore-knowledge, from years to a decade or more in advance, of impactful climate changes that can in general be attributed to a combination of internal and externally forced variability. Predictions initialized using observations of past climate states are tested by comparing their ability to reproduce past climate evolution with that of uninitialized simulations in which the same radiative forcings are applied. A new set of decadal prediction (DP) simulations has recently been completed using the Community Earth System Model (CESM) and is now available to the community. This new large-ensemble (LE) set (CESM-DPLE) is composed of historical simulations that are integrated forward for 10 years following initialization on 1 November of each year between 1954 and 2015. CESM-DPLE represents the “initialized” counterpart to the widely studied CESM Large Ensemble (CESM-LE); both simulation sets have 40-member ensembles, and they use identical model code and radiative forcings. Comparing CESM-DPLE to CESM-LE highlights the impacts of initialization on prediction skill and indicates that robust assessment and interpretation of DP skill may require much larger ensembles than current protocols recommend. CESM-DPLE exhibits significant and potentially useful prediction skill for a wide range of fields, regions, and time scales, and it shows widespread improvement over simpler benchmark forecasts as well as over a previous initialized system that was submitted to phase 5 of the Coupled Model Intercomparison Project (CMIP5). Finally, the new DP system offers new capabilities that will be of interest to a broad community pursuing Earth system prediction research.

54 ENVIRONMENTAL SCIENCES↗

Vibrationally resolved high-resolution NEXAFS and XPS spectra of phenanthrene and coronene

We performed a combined experimental and theoretical study of the C1s Near-Edge X-ray Absorption Fine-Structure (NEXAFS) spectroscopy and X-ray Photoelectron Spectroscopy in the gas phase of two polycyclic aromatic hydrocarbons (phenanthrene and coronene), typically formed in combustion reactions. In the NEXAFS of both molecules, a double-peak structure appears in the C1s → LUMO region, which differ by less than 1 eV in transition energies. The vibronic coupling is found to play an important role in such systems. It leads to weakening of the lower-energy peak and strengthening of the higher-energy one because the 0 − n (n > 0) vibrational progressions of the lower-energy peak appear in nearly the same region of the higher-energy peak. Vibrationally resolved theoretical spectra computed within the Frank-Condon (FC) approximation and linear coupling model agree well with the high-resolution experimental results. We find that FC-active normal modes all correspond to in-plane vibrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗