Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiation absorption”

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

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

H α emission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Direct Measurement of the Selective Microwave-Induced Heating of Agglomerates of Dipolar Molecules: The Origin of and Parameters Controlling a Microwave Specific Superheating Effect

Agglomerates of polar molecules in nonpolar solvents are selectively heated by microwave radiation. The magnitude of the selective heating was directly measured by using the temperature dependence of the intensities of the Stokes and anti-Stokes bands in the Raman spectra of p-nitroanisole (pNA) and mesitylene. Under dynamic heating conditions, a large apparent temperature difference (ΔT) of over 100 °C was observed between the polar pNA solute and the nonpolar mesitylene solvent. This represents the first direct measurement of the selective microwave heating process. The magnitude of the selective microwave heating was affected by the properties of the agglomerated pNA. As the concentration of the pNA increases, the magnitude of the selective heating of the pNA was observed to decrease. This is explained by the tendency of the pNA dipoles to orient in an antiparallel fashion in the aggregates as measured by the Kirkwood g value, which decreased with increasing concentration. This effect reduces the net dipole moment of the agglomerates, which decreases the microwave absorption. After the radiation was terminated, the effective temperature of the dipolar molecules returned slowly to that of the medium. The slow heat transfer was modeled successfully by treating the solutions as a biphasic solvent/solute system. Based on modeling and the fact that the agglomerate can be heated above the boiling temperature of the solvent, an insulating layer of solvent vapor is suggested to form around the heated agglomerate, slowing convective heat transfer out of the agglomerate. Furthermore, this is an effect unique to microwave heating.

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

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Orbital Hybridization and Magnetic Coupling at Cuprate–Manganite Interfaces Driven by Manganite Doping

Fundamental understanding of interface mechanisms that generate unexpected physical properties in the cuprate/manganite heterostructures is essential for possible applications as spintronic devices. In this study, CaCuO 2 /La 1–x Sr x MnO 3 (CCO/LSMO) superlattices are investigated, where the infinite-layer cuprate CCO does not have apical oxygen for Cu in the Ca-plane so that the Mn–O–Cu superexchange coupling can be only present at the MnO 2 –(La,Sr)O–CuO 2 interface. Two different doping states for the manganites are studied, namely x = 0.1 and x = 0.3, corresponding to the ferromagnetic insulating and ferromagnetic metallic states, respectively, in the manganite phase diagram. Linear and circular dichroism in X-ray absorption by synchrotron radiation at Cu and Mn L-edges clearly demonstrate that in the absence of apical oxygen in the cuprate block, the magnetic coupling between the LSMO and CCO is weakly ferromagnetic when LSMO is metallic, while it is antiferromagnetic when LSMO is insulating, along with an increased Cu 3d(3z 2 –r 2 ) and Mn 3d(3z 2 –r 2 ) orbital occupation. It is proposed that the Mn 3d valence band upward shift driven by the enhancement of Mn 3+ content in the underdoped sample enhances the orbital hybridization. The stronger hybridization at the interface MnO 2 –(La,Sr)O–CuO 2 gives rise to the antiferromagnetic coupling between cuprate and underdoped manganite.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strong Local Evaporative Cooling Over Land Due to Atmospheric Aerosols

Abstract Aerosols can enhance terrestrial productivity through increased absorption of solar radiation by the shaded portion of the plant canopy—the diffuse radiation fertilization effect. Although this process can, in principle, alter surface evaporation due to the coupling between plant water loss and carbon uptake, with the potential to change the surface temperature, aerosol‐climate interactions have been traditionally viewed in light of the radiative effects within the atmosphere. Here, we develop a modeling framework that combines global atmosphere and land model simulations with a conceptual diagnostic tool to investigate these interactions from a surface energy budget perspective. Aerosols increase the terrestrial evaporative fraction, or the portion of net incoming energy consumed by evaporation, by over 4% globally and as much as ∼40% regionally. The main mechanism for this is the increase in energy allocation from sensible to latent heat due to global dimming (reduction in global shortwave radiation) and slightly augmented by diffuse radiation fertilization. In regions with moderately dense vegetation (leaf area index >2), the local surface cooling response to aerosols is dominated by this evaporative pathway, not the reduction in incident radiation. Diffuse radiation fertilization alone has a stronger impact on gross primary productivity (+2.18 Pg C y −1 or +1.8%) than on land evaporation (+0.18 W m −2 or +0.48%) and surface temperature (−0.01 K). Our results suggest that it is important for land surface models to distinguish between quantity (change in total magnitude) and quality (change in diffuse fraction) of radiative forcing for properly simulating surface climate.

Chakraborty, TC↗
DOE OSTI2020

Non-invasive measurement of leaf water content and pressure–volume curves using terahertz radiation

In this paper we describe a non-invasive method of measuring leaf water content using THz radiation and combine this with psychrometry for determination of leaf pressure–volume relationships. In contrast to prior investigations using THz radiation to measure plant water status, the reported method exploits the differential absorption characteristic of THz radiation at multiple frequencies within plant leaves to determine absolute water content in real-time. By combining the THz system with a psychrometer, pressure–volume curves were generated in a completely automated fashion for the determination of leaf tissue water relations parameters including water potential at turgor loss, osmotic potential at full turgor and the relative water content at the turgor loss point. This novel methodology provides for repeated, non-destructive measurement of leaf water content and greatly increased efficiency in generation of leaf PV curves by reducing user handling time.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Valence shell electronically excited states of norbornadiene and quadricyclane

The absolute photoabsorption cross sections of norbornadiene (NBD) and quadricyclane (QC), two isomers with chemical formula C 7 H 8 that are attracting much interest for solar energy storage applications, have been measured from threshold up to 10.8 eV using the Fourier transform spectrometer at the SOLEIL synchrotron radiation facility. The absorption spectrum of NBD exhibits some sharp structure associated with transitions into Rydberg states, superimposed on several broad bands attributable to valence excitations. Sharp structure, although less pronounced, also appears in the absorption spectrum of QC. Assignments have been proposed for some of the absorption bands using calculated vertical transition energies and oscillator strengths for the electronically excited states of NBD and QC. Natural transition orbitals indicate that some of the electronically excited states in NBD have a mixed Rydberg/valence character, whereas the first ten excited singlet states in QC are all predominantly Rydberg in the vertical region. In NBD, a comparison between the vibrational structure observed in the experimental 1 1 B 1 –1 1 A 1 (3sa 1 ← 5b 1 ) band and that predicted by Franck–Condon and Herzberg–Teller modeling has necessitated a revision of the band origin and of the vibrational assignments proposed previously. Similar comparisons have encouraged a revision of the adiabatic first ionization energy of NBD. Simulations of the vibrational structure due to excitation from the 5b 2 orbital in QC into 3p and 3d Rydberg states have allowed tentative assignments to be proposed for the complex structure observed in the absorption bands between ~5.4 and 7.0 eV.

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

Probing the optical properties and toxicological profile of zinc tungstate nanorods

Zinc tungstate is a semiconductor known for its favorable photocatalytic, photoluminescence, and scintillation properties, coupled with its relatively low cost, reduced toxicity, and high stability in biological and catalytic environments. In particular, zinc tungstate evinces scintillation properties, namely the ability to emit visible light upon absorption of energetic radiation such as x rays, which has led to applications not only as radiation detectors but also for biomedical applications involving the delivery of optical light to deep tissue, such as photodynamic therapy and optogenetics. Here, we report on the synthesis of zinc tungstate nanorods generated via an optimized but facile method, which allows for synthetic control over the aspect ratio of the as-synthesized anisotropic motifs via rational variation of the solution pH. Additionally, we investigate the effect of aspect ratio on their resulting photoluminescent and radioluminescent properties. We further demonstrate the potential of these zinc tungstate nanorods for biomedical applications, such as photodynamic therapy for cancer treatment, by analyzing their toxicological profile within cell lines and neurons.

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

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
DOE OSTI2022

Use of diamond sensors for a high-flux, high-rate X-ray pass-through diagnostic

X-ray free-electron lasers (XFELs) deliver pulses of coherent X-rays on the femtosecond time scale, with potentially high repetition rates. While XFELs provide high peak intensities, both the intensity and the centroid of the beam fluctuate strongly on a pulse-to-pulse basis, motivating high-rate beam diagnostics that operate over a large dynamic range. The fast drift velocity, low X-ray absorption and high radiation tolerance properties of chemical vapour deposition diamonds make these crystals a promising candidate material for developing a fast (multi-GHz) pass-through diagnostic for the next generation of XFELs. A new approach to the design of a diamond sensor signal path is presented, along with associated characterization studies performed in the XPP endstation of the LINAC Coherent Light Source (LCLS) at SLAC. Qualitative charge collection profiles (collected charge versus time) are presented and compared with those from a commercially available detector. Quantitative results on the charge collection efficiency and signal collection times are presented over a range of approximately four orders of magnitude in the generated electron–hole plasma density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Emerging low-cloud feedback and adjustment in global satellite observations

From mid-2003 to mid-2024, a global decrease in low-cloud amount enhanced the absorption of solar radiation by 0.22±0.07 W m −2 per decade (±1σ range), accelerating the energy imbalance trend during that period (0.44 W m −2 per decade). Through controlling factor analysis, here we show that the low-cloud trend is due to a combination of cloud feedback and adjustments to greenhouse gases and aerosols (respectively 0.09±0.02, 0.05±0.03, and 0.03±0.03 W m −2 per decade), which jointly account for 74 % of the trend. The contribution of natural climate variability is weak but uncertain (0.01±0.08 W m −2 per decade), owing to a poorly constrained trend in boundary-layer inversion strength. Importantly, the observed low-cloud radiative trend lies well within the range of values simulated by contemporary global climate models under conditions close to present day. Any systematic model error in the representation of present-day global energy imbalance trends is thus likely to originate in processes unrelated to low clouds.

Geosciences↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2023

Impact of Wildfires on Solar Resource Availability in California in a Changing Climate

Wildfires can emit large amounts of atmospheric particulate matters and influence not only air quality but also availability of photovoltaic (PV) generation due to scattering and absorption of solar radiation. Under anthropogenic changing climate, wildfire activity is projected to increase over western North America due to drier and warmer climate, implying increasing impact on solar resource and larger uncertainty in solar generation especially in regions with faster PV penetration. This study focuses on quantifying the impact of wildfires on aerosol optical depth (AOD) and thus solar resource over California using National Solar Radiation Database (NSRDB), developed by the National Renewable Energy Laboratory (NREL). This assessment includes historical analysis and estimation of solar resource under wildfire scenarios (2020 wildfire and an enhanced wildfire scenario based on 2020 wildfires). Historical analysis for the period of July-October 2019/2020 (low/high fire activity period) shows that the averaged global horizontal irradiance (GHI) and direct normal irradiance (DNI) are reduced by about 30 and 90 W m-2 (5.6 and 15.3%), respectively, during high fire activity period. To create AOD dataset for enhanced fire scenario, 165% increase in burn area (in around 2050) is selected based on comprehensive literature review, which is further applied to the Fire INventory from NCAR (FINN) fire emission. WRF-Chem model with the enhanced FINN emissions is used to simulate and represent a preserving spatial distribution of burned areas and wildfire-emitted aerosols. Our initial analysis suggests that the"enhanced 2020 wildfire" AOD can increase by a factor of 1.3 - 2, which can significantly reduce solar irradiance and increase uncertainty in generation and reliability of power system in extreme wildfire events under a high solar penetration scenario. Overall, this study provides an estimate of the impacts of wildfires on solar resource to make informed decisions on reserve planning, generation scheduling, and reliability investments.

aerosol optical depth↗
DOE OSTI2023

Impact of Actinide Complexation on Ligand Radiolysis

Although the actinides boast many unique physical and chemical properties, their inherent susceptibility to radioactive decay are what make them truly interesting elements to study. The absorption of ionizing radiation from actinide decay leads to the formation of a variety of transient and steady-state radicals, ions, and molecular radiolysis products that can lead to significant changes in the surrounding environment, and ultimately dictate steady-state actinide redox distributions and the longevity of molecules designed for actinide complexation. Radiolysis of the latter leads to complexant destruction and the concomitant formation of degradation products that can complicate actinide studies and processes. However, the radiation chemistry of most actinide complexants have been studied in the absence of the actinides they were designed to complex, which can lead to inaccurate conclusions on longevity and degradation product distributions, as metal ion complexation has been historically shown to influence a ligand’s radiolytic behavior. Consequently, bridging this knowledge gap is important for actinide science. Here, I will discuss the impact of actinide complexation on the steady-state and time-resolved radiation-induced reactivity of a variety of complexants, including, tributyl phosphate (TBP) and N,N,N',N'-tetraoctyl diglycolamide (TODGA) .

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

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗