Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moderator”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Influence of support on Rh-Co bimetallic catalysts for ethylene hydroformylation

Bimetallic Rh- and Co-based catalysts are promising materials for the heterogenous ethylene hydroformylation reaction. Here, in this study, the influence of a mesoporous silica (SBA-15) support on Rh and Co bimetallic interactions was investigated through a comparison with silica gel and alumina supports. The bimetallic catalyst supported on mesoporous silica (RhCo 3 /SBA-15) showed the best C 3 oxygenate yield among the three bimetallic catalysts. In-situ vibrational studies suggested moderate binding of the gem-dicarbonyl and Rh(CO)(C 2 H 4 ) intermediates due to this bimetallic interaction that led to improved hydroformylation performance. Kinetic studies revealed a lower hydroformylation barrier for the bimetallic compared to a Rh monometallic catalyst, and in-situ X-ray absorption spectroscopy investigations showed clear Rh-Co alloy formation on RhCo 3 /SBA-15. The bimetallic enhancement effect from the interaction with the mesoporous silica support shown here can be further optimized to design olefin hydroformylation catalysts.

58 GEOSCIENCES↗
DOE OSTI2025

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Probing the shear-induced microstructure of a smectite clay aqueous suspensions using rheo-USANS and rheo-SIPLI measurements

Hypothesis: The static microstructure of aqueous sodium-montmorillonite (Na-Mt) suspensions at low ionic strengths (where Particle Size/Debye Length ≈1) exhibits both the particle–particle ordering as well as aggregation with repulsive ordered domains having characteristic optical birefringence and attractive aggregated entities larger than 20 µm resulting in ever-increasing yield stresses also known as physical aging-rejuvenation behavior. We hypothesize that the attractive particle–particle aggregation is the underlying cause behind the physical aging-rejuvenation behavior observed in Na-Mt suspensions with no contribution from structural dynamics driven by repulsive particle–particle ordering or jamming. Experiments: We investigate the shear-induced microstructure of aqueous Na-Mt suspensions in the sol and gel state using rheo-ultra-small angle neutron scattering (rheo-USANS) experiments at shear rates of 1, 50, 500, and 2000 s −1 . We also perform rheo-shear-induced polarization light imaging (rheo-SIPLI) experiments to relate ordering with shearing and aging. Findings: Shearing the suspensions at low to moderate shear rates induces particle–particle aggregation and shearing at high shear rates induces the breakage of particle–particle aggregation in the sol and gel states, suggesting the microstructural aggregation in the sol and gel state is shear sensitive and a full rejuvenation or breakage of particle–particle aggregation is only achieved at a minimum critical shear rate. The rheo-SIPLI experiments reveal that the sol and gel state exhibited strong Maltese cross patterns at a shear rate of 1000 s −1 , indicating particle–particle ordering. Post shearing, the gel exhibited temporal evolution of storage modulus without any noticeable influence on the appearance of the Maltese cross pattern indicating physical aging and particle ordering are distinct length scale phenomena in Na-Mt suspensions and the physical aging-rejuvenation behaviour is a feature of particle–particle aggregation as opposed to ordering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Home charging for all: Techno-economic and life cycle assessment of multi-unit dwelling electric vehicle charging hubs

Ubiquitous electric vehicle adoption can drastically reduce greenhouse gas emissions (GHG) but will require equitable home charging infrastructure for all residences. Unlike single-family homes, Multi-Unit Dwellings (MUD) currently lack market share and access to nearby charging infrastructure partly due to expensive capital costs for reluctant residential property-owners. Further, this work evaluates the levelized cost of charging (LCOC) for Battery Electric Vehicles (BEVs) at MUD community charging hubs through a techno-economic analysis (TEA) that leverages real-world charging data and costs. Three different MUD charging types are investigated for a baseline and optimistic case: 1.9-kW Level 1 (L1), 6.6-kW Level 2 (L2), and 50-kW Direct Current Fast Charging (DCFC) stations under three different ownership models: resident, utility, and private company. Results demonstrate L1 and L2 chargers to be less expensive than the gasoline equivalent across the United States. Further, utility and private company ownership models, which avoid initial costs for residential property owners, result in a large LCOC premium for L1 and moderate LCOC premium for L2 relative to the resident ownership model. In contrast, DCFC is shown to be expensive for baseline scenarios but economical for optimistic scenarios especially under private company ownership. This work also performs a cradle to grave (C2G) life cycle assessment (LCA) of an average passenger BEV and gasoline conventional vehicle (CV) using yearly (grid mix & vehicle parameters), hourly (grid mix), and state-level (grid mix) resolution. Results show BEVs to have lower GHG emissions (-86% to-10%) than gasoline CVs in the contiguous U.S. Next, the system boundary of the TEA is extended to the total cost of ownership (TCO) of BEVs. The TCO is then coupled with the C2G GHG emissions to calculate the cost of GHG emissions reduction. Ultimately, the cost of GHG emissions reduction from MUD BEVs relative to gasoline CVs is shown to be negative for every scenario except baseline DCFC, meaning MUD BEV charging infrastructure can be a cost-effective endeavor to reduce GHG emissions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Neural-network based collision operators for the Boltzmann equation

Kinetic gas dynamics in rarefied and moderate-density regimes have complex behavior associated with collisional processes. These processes are generally defined by convolution integrals over a high-dimensional space (as in the Boltzmann operator), or require evaluating complex auxiliary variables (as in Rosenbluth potentials in Fokker-Planck operators) that are challenging to implement and computationally expensive to evaluate. In this work, we develop a data-driven neural network model that augments a simple and inexpensive BGK collision operator with a machine-learned correction term, which improves the fidelity of the simple operator with a small overhead to overall runtime. The composite collision operator has a tunable fidelity and, in this work, is trained using and tested against a direct-simulation Monte-Carlo (DSMC) collision operator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Efficient particle control in systems with large density gradients

Simulations of large density gradients present a number of challenges for direct Monte Carlo methods, since they lead to too few particles in dilute regions and prohibitively many in the dense regions. Here, we propose a particle control methodology that gives the user more control of the number of particles per cell by introducing a variable weight for each particle. The proposed scheme is based on the stochastic weighted particle method, requires only small modifications to DSMC, and exactly conserves mass, momentum, and energy. In validation tests of systems with density ratios of order 102-104, we observe 101-102 times less variance in the dilute region compared to a DSMC solution with the same number of system particles, while introducing a moderate additional computational cost.

DSMC↗
DOE OSTI2022

Morphological and Phenotypic Characteristics of the Bovine Nasopharyngeal Mucosa and Associated Lymphoid Tissue

The mammalian nasopharynx is an anatomically complex region of the upper respiratory tract that directly communicates with the nasal cavity, laryngopharynx, oesophagus and trachea. The nasopharyngeal mucosa contains moderate quantities of mucosa-associated lymphoid tissue (MALT) that is appropriately located for immunological sampling but also creates vulnerability to pathogens. In recent years, the nasopharynx has been inculpated in the pathogenesis of important diseases of cattle (foot-and-mouth disease) and humans (COVID-19), yet the tissue has never been described in detail in any species. In order to characterize the morphology and cellular composition of the bovine nasopharynx, samples of mucosa were collected from the nasopharynx of five 8–13-month-old steers and examined using light microscopy, immunohistochemistry and multichannel immunofluorescence. Morphologically, the nasopharyngeal epithelium was highly heterogeneous, with a continuum ranging from stratified squamous epithelium to highly attenuated, follicle-associated epithelium (FAE). Distribution of MALT was similarly regionally variable ranging from absent to clusters of multiple lymphoid follicles. Phenotypic characterization demonstrated dense distributions of dendritic cells and T lymphocytes surrounding lymphoid follicles, which comprised mostly B lymphocytes. The FAE overlaying the lymphoid follicles also contained higher numbers of dendritic cells and lymphocytes compared with the adjacent non-lymphoid epithelium, although cytotoxic T cells were notably scarce in the FAE. The bovine nasopharyngeal lymphoid tissue had comparable elements to other MALTs with specific differences that may help to elucidate the pathogenesis of infectious agents that have specific tropism for this tissue.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Advanced high-performance large diameter Cs 2 HfCl 6 (CHC) and mixed halides scintillator

This paper reports on successful growth and performance evaluation of two large diameter Cs 2 HfCl 6 (CHC) and Cs 2 HfCl 4 Br 2 (CHCB), both recently developed scintillator crystals. The discovery of Cs 2 HfCl 6 (CHC) as a scintillator has lately generated much interest in this material and its family, which belongs to the K 2 PtCl 6 cubic crystal structure. CHC is an intrinsic scintillator that is non-hygroscopic (slightly deliquescent, i.e., slightly moisture sensitive but does not deteriorate much when left in air), has no self-radioactivity, provides excellent energy resolution, and has excellent non-proportionality. CHC has a moderate density of 3.9 g/cm 3 and an effective atomic number of 58. Reported in this paper are transparent crack-free single crystal CHC and CHCB boules of one inch in diameter, both grown using the vertical Bridgman method. Samples retrieved from the boules, sized 23 mm × 30 mm and 23 mm × 26 mm, respectively, are characterized for their optical and scintillation properties. Energy resolutions of 3.5% and 3.7% (FWHM) at 662 keV, respectively, are reported. Light yields of 23,000 ph/MeV and 20,000 ph/MeV have been calculated for the large diameter CHC and CHCB crystals, respectively. Results comparable to previously reported results for smaller crystals have been obtained. Studies on decay time, non-proportionality, as well as detector characterization are also reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Characterization of Spent Purolite A530E Resin with Implications for Long-Term Radioactive Contaminant Removal

Direct removal of pertechnetate anion (TcO4-) from the subsurface contaminated plumes at legacy nuclear weapons production sites, is effective remediation strategy to prevent their spread to adjacent aquifers but is also challenging because TcO4- is a trace component of the contaminated groundwater typically containing large access of other anionic constituents, such nitrate, sulfate, chloride and others. It can be achieved through ion exchange treatment using Purolite A530E resin which to date however has only been evaluated under laboratory conditions and many questions regarding its long-term operational performance remain. To address this need, here we report comprehensive characterization of the spent Purolite A530E resin which processed over 5.38x109 L of contaminated groundwater and successfully removed about 3.78 Ci of Tc-99 during four years of operation at the 200 West Pump & Treat facility at the US DOE Hanford site. This Tc-99 loading constitutes however only about 0.034% of the theoretical capacity of the resin which retained significant amounts of sulfate. Among other radioactive contaminants, small quantities of U-238, Co-60, and I-129 were retained by the resin. The total loading of iodine (combined 1-127 and I-129) on the spent resin exceeded that of Tc-99. To elucidate the mechanism of iodine retention, ion exchange behavior of iodide and iodate was investigated. Purolite A530E resin exhibited highly efficient uptake of iodide and only moderate affinity for iodate in accord with their Gibbs energy of hydration. Sorption isotherms for both anions obeyed Freundlich model.

Levitskaia, Tatiana G.↗
DOE OSTI2021

Dual-atom active sites embedded in two-dimensional C 2 N for efficient CO 2 electroreduction: A computational study

Double-atom catalysts (DACs) have emerged as an enhanced platform of single-atom catalyst for promoting electrocatalytic CO 2 reduction reaction (CO 2 RR). Herein, we present a density-functional theory study on CO 2 RR performance of seven C 2 N-supported homo- and heteronuclear DACs, denoted as M 2 @C 2 N. Our results demonstrate that there exists substantial synergistic effect of dual-metal-atom N 2 M 2 N 2 active site and C 2 N matrix on O&#x2550 C &#x2550O bond activation. The dual-atom M 2 sites are able to drive CO 2 RR beyond C 1 products with low limiting potential (U L ). Specifically, C 2 H 4 formation is preferred on FeM@C 2 N (M = Fe, Co, Ni, Cu) versus CH 4 formation on CuM@C 2 N (M = Co, Ni, Cu). Furthermore, *CO+*CO co-binding strength can serve as a descriptor for CO 2 RR activity for making C 2 products such that the moderate binding results in the lowest U L . Remarkably, C-affinity matters most to C—C bond coupling and C 2 H 4 formation while both C- and O-affinity control CH 4 formation. Furthermore, our results provide theoretical insight into rational design of DACs for efficient CO 2 RR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Revisiting the capacity-fading mechanism of P2-type sodium layered oxide cathode materials during high-voltage cycling

We report P2-type sodium layered oxide cathode (Na 2/3 Ni 1/3 Mn 2/3 O 2 , P2-NNMO) has attracted great attention as a promising cathode material for sodium ion batteries because of its high specific capacity. However, this material suffers from a rapid capacity fade during high-voltage cycling. Several mechanisms have been proposed to explain the capacity fade, including intragranular fracture caused by the P2-O2 phase transion, surface structural change, and irreversible lattice oxygen release. Here we systematically investigated the morphological, structural, and chemical changes of P2-NNMO during high-voltage cycling using a variety of characterization techniques. It was found that the lattice distortion and crystal-plane buckling induced by the P2-O2 phase transition slowed down the Na-ion transport in the bulk and hindered the extraction of the Na ions. The sluggish kinetics was the main reason in reducing the accessible capacity while other interfacial degradation mechanisms played minor roles. Our results not only enabled a more complete understanding of the capacity-fading mechanism of P2-NNMO but also revealed the underlying correlations between lattice doping and the moderately improved cycle performance.

25 ENERGY STORAGE↗
DOE OSTI2022

Tracking the environmental impacts of ecological engineering on coastal wetlands with numerical modeling and remote sensing

Coastal wetlands are the most valuable ecosystems on the earth but facing severe degradation and losses owing to climate change and anthropogenic activities. Many ecological engineering projects (EEP) have been conducted to mitigate the degradation of coastal wetlands. However, the geomorphological impacts of EEP on coastal wetlands have not been well documented. In this study, a method employed a process-based hydrodynamic model and remote sensing (RS) was developed to evaluate the impacts of EEP on the geomorphological change of a prototype Ramsar site. Results demonstrated that RS has great potential in improving the quality of bathymetry data for the numerical model with a decrease of RMSE from 0.52 m to 0.3 m. It also showed good capacity in trend detection of geomorphological change spatially. Results showed the Chongming Dongtan wetland experienced erosion with an annual rate of -0.035 m/yr from 2013 to 2016 after the implementation of EEP. The deposition rate changed significantly in the area within 200 m of the EEP. It is found that the EEP modified the composition of vegetation, sediment transportation, as well as substrate stability, affecting the geomorphological change of coastal wetlands. It is suggested that the EEP with moderate anthropogenic disturbance is a direct and effective way to recover the coastal habitats for waterbirds. However, the modification of the coastal wetland ecosystem by EEP will lead to the potential vulnerability to global climate change. Therefore, how to mitigate the advantages and disadvantages of the EEP is needed to be further studied to find a more sustainable way for coastal management.

geomorphological change, coastal wetlands, Ecologi↗
DOE OSTI2022

Mitigating life cycle GHG emissions of roads to be built through 2030: Case study of a Chinese province

We report the expansion of road networks in emerging economies such as China causes significant greenhouse gas (GHG) emissions. This development is conflicting with China's commitment to achieve carbon neutrality. Thus, there is a need to better understand life cycle emissions of road infrastructure and opportunities to mitigate these emissions. Existing impact studies of roads in developing countries do not address recycled materials, improved pavement maintenance, or pavement-vehicle interaction and electric vehicle (EV) adoption. Combining firsthand information from Chinese road construction engineers with publicly available data, this paper estimates a comprehensive account of GHG emissions of the road pavement network to be constructed in the next ten years in the Shandong province in Northern China. Further, we estimate the potential of GHG emission reductions achievable under three scenario sets: maintenance optimization, alternative pavement material replacement, and EV adoption. Results show that the life cycle GHG emissions of highways and Class 1–4 roads to be constructed in the next 10 years amount to 147 Mt CO2-eq. Considering the use phase in our model reveals that it is the dominant stage in terms of emissions, largely due to pavement-vehicle interaction. Vehicle electrification can only moderately mitigate these emissions. Other stages, such as materials production and road maintenance and rehabilitation, contribute substantially to GHG emissions as well, highlighting the importance of optimizing the management of these stages. Surprisingly, longer, not shorter maintenance intervals, yield significant emission reductions. Another counter-intuitive finding is that thicker and more material-intensive pavement surfaces cause lower emissions overall. Taken together, optimal maintenance and rehabilitation schedules, alternative material use, and vehicle electrification provide GHG reduction potentials of 11%, 4%–16% and 2%–6%, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Third international challenge to model the medium- to long-range transport of radioxenon to four Comprehensive Nuclear-Test-Ban Treaty monitoring stations

In 2015 and 2016, atmospheric transport modeling challenges were conducted in the context of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) verification, however, with a more limited scope with respect to emission inventories, simulation period and number of relevant samples (i.e., those above the Minimum Detectable Concentration (MDC)) involved. Therefore, a more comprehensive atmospheric transport modeling challenge was organized in 2019. Stack release data of Xe-133 were provided by the Institut National des Radioéléments/IRE (Belgium) and the Canadian Nuclear Laboratories/CNL (Canada) and accounted for in the simulations over a three (mandatory) or six (optional) months period. Best estimate emissions of additional facilities (radiopharmaceutical production and nuclear research facilities, commercial reactors or relevant research reactors) of the Northern Hemisphere were included as well. Model results were compared with observed atmospheric activity concentrations at four International Monitoring System (IMS) stations located in Europe and North America with overall considerable influence of IRE and/or CNL emissions for evaluation of the participants’ runs. Participants were prompted to work with controlled and harmonized model set-ups to make runs more comparable, but also to increase diversity. It was found that using the stack emissions of IRE and CNL with daily resolution does not lead to better results than disaggregating annual emissions of these two facilities taken from the literature if an overall score for all stations covering all valid observed samples is considered. A moderate benefit of roughly 10% is visible in statistical scores for samples influenced by IRE and/or CNL to at least 50% and there can be considerable benefit for individual samples. Effects of transport errors, not properly characterized remaining emitters and long IMS sampling times (12–24 h) undoubtedly are in contrast to and reduce the benefit of high-quality IRE and CNL stack data. Complementary best estimates for remaining emitters push the scores up by 18% compared to just considering IRE and CNL emissions alone. Despite the efforts undertaken the full multi-model ensemble built is highly redundant. An ensemble based on a few arbitrary runs is sufficient to model the Xe-133 background at the stations investigated. The effective ensemble size is below five. An optimized ensemble at each station has on average slightly higher skill compared to the full ensemble. However, the improvement (maximum of 20% and minimum of 3% in RMSE) in skill is likely being too small for being exploited for an independent period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)↗