Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative diffusion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Consistent solutions of the radiation diffusion equation in spherical and cylindrical geometries

We have extended the radiation diffusion model of Hammer and Rosen [Phys. Plasmas 10, 1829 (2003)] to diverging spherical and cylindrical geometries. The effect of curvilinear geometry on the supersonic, expanding wavefront increases as the internal radius of a spherical or cylindrical shell approaches zero. Small spherical geometries are important for modeling systems at the size scale of inertial confinement fusion capsules, at these scales existing quasi-analytic models for planar geometry significantly disagree with the results of simulation. With this method, the benefits of rapid iteration can be applied to common spherical systems at much smaller length scales. We present comparisons between numerical diffusion solutions and the analytic model to give ranges of applicability for the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

A Case Study of View-Factor Rectification Procedures for Diffuse-Gray Radiation Enclosure Computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Robert P Taylor↗
DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Rethinking agrivoltaic incentive programs: A science-based approach to encourage practical design solutions

Agrivoltaic systems are promising solutions to address global food and energy challenges by combining agriculture and solar photovoltaics. However, the lack of appropriate regulations to define and guide their implementation constrains the growth of agrivoltaic systems in the U.S. This study uses a shading and radiation tool to evaluate an existing agrivoltaic incentive program that defines agrivoltaic designs based on shading reduction limits and panel height requirements. Our analysis indicates that structuring policy requirements around shading, and not light availability, may lead to an underestimation of crop suitability by neglecting diffuse radiation. Furthermore, we show that agrivoltaic systems can avoid increasing panel height if policy acknowledges use-case scenarios where farming only occurs between rows. In light of these insights, this study proposes two key policy recommendations: (1) benchmark crop suitability based on daily light integral (DLI) requirements for a shade-intolerant crop selected to represent a prevalent crop in the region, and (2) include an incentive scenario where agriculture is only required between rows. Furthermore, these two recommendations can potentially incentivize designs that are practical and closer in cost to conventional solar farms, thereby accelerating the adoption of cost-effective agrivoltaic systems.

14 SOLAR ENERGY↗
DOE OSTI2025

Flux effects in precipitation under irradiation simulation of Fe-Cr alloys

Radiation-enhanced precipitation of Cr-rich α’ in irradiated Fe-Cr alloys, which results in hardening and embrittlement, depends on the irradiating particle and the displacement per atom (dpa) rate. Here, we utilize a Cahn-Hilliard phase-field based approach, that includes simple models for nucleation, irradiating particle and rate dependent radiation-enhanced diffusion and cascade mixing to simulate α’ evolution under neutrons, heavy ions, and electron irradiations. Different irradiating particles manifest very different cascade mixing efficiencies. The model was calibrated using neutron data. For cascade inducing neutron/heavy-ion dpa rates at 300 °C between 10-8 and 10-6 dpa/s the model predicts approximately constant number density, decreasing radius, decreasing α’ Cr composition, and lower α’ volume fraction. The model then predicts a dramatic transition to no α’ formation above approximately 10-5 dpa/s, while electron irradiation, with weak mixing, had little effect at dpa rates up to 10-3 dpa/s. These model predictions are consistent with experiments. We explain the results in terms of the flux dependence of the radiation enhanced diffusion, cascade mixing, and their ratio, which all vary significantly in relevant flux ranges for neutron and cascade inducing ion irradiations. These results show that both cascade mixing and radiation enhanced diffusion must be accounted for when attempting to emulate neutron-irradiation effects using accelerated ion irradiations.

Ke, Jia-Hong↗
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2025

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1997

Radiation Temperature and Extinction of Transient Gaseous Diffusion Flames in Microgravity

The absence of buoyancy-induced flows in micro-g and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between 1-g and micro-g flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment which changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation (and its interaction with flame chemistry), that are often ignored under normal gravity, become very important and sometimes even controlling. This is particularly true for conditions at extinction of a micro-g diffusion flame.

Arvind Atreya↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1995

Numerical Solution of Fluid Flow and Heat Tranfer Problems with Surface Radiation

This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

S Ahuja↗
DOE OSTI2025

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗