Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiation Scattering Equations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Radiative classical gravitational observables at $ \mathcal{O} $(G 3 ) from scattering amplitudes

We compute classical gravitational observables for the scattering of two spinless black holes in general relativity and N =8 supergravity in the formalism of Kosower, Maybee, and O’Connell (KMOC). We focus on the gravitational impulse with radiation reaction and the radiated momentum in black hole scattering at $ \mathcal{O} $(G 3 ) to all orders in the velocity. These classical observables require the construction and evaluation of certain loop-level quantities which are greatly simplified by harnessing recent advances from scattering amplitudes and collider physics. In particular, we make use of generalized unitarity to construct the relevant loop integrands, employ reverse unitarity, the method of regions, integration-by-parts (IBP), and (canonical) differential equations to simplify and evaluate all loop and phase-space integrals to obtain the classical gravitational observables of interest to two-loop order. The KMOC formalism naturally incorporates radiation effects which enables us to explore these classical quantities beyond the conservative two-body dynamics. From the impulse and the radiated momentum, we extract the scattering angle and the radiated energy. Finally, we discuss universality of the impulse in the high-energy limit and the relation to the eikonal phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect-induced phonon-resonant scattering and its influence on thermal transport of irradiated thorium-dioxide

Thermal transport in proton irradiated thorium-dioxide (ThO 2 ) is investigated. Using a combination of experiments and first-principles computational framework, the role of lattice defects on thermal conductivity is analyzed. A resonant-phonon scattering mechanism beyond the traditionally considered Rayleigh scattering is found to significantly influence low-temperature thermal transport in the presence of irradiation-induced point defects. The existence of localized phonon modes associated with irradiation-induced defects is suggested by the inability of the first-principles based thermal conductivity model—which considers only three-phonon interactions and phonon-defects scattering using the Tamura formalism—to predict the experimental results, unless a resonant scattering mechanism is included. The emergence of additional peaks in the Raman spectra in the proximity of phonon-resonant frequency provides further evidence for the existence of localized modes. Coupled with a microstructure evolution model, this analysis enables more accurate analysis for contrasting the contributions of different phonon scattering mechanisms across all irradiation doses and temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Stretched and Filtered Multigroup Pn Transport for Improved Positivity and Accuracy

One of the common methods for solving the radiation transport equation is to use a polynomial expansion for the angle variable(s). Recent research that reduces the oscillations and improves the positivity of the gray transport equation solutions is here applied to the multigroup transport equations. Constant scale factors that stretch the time axis and constant scattering opacities that filter the solution greatly increase the accuracy of the solution with no added nonlinearities. No new solution techniques are required. Test problems are presented in one and two dimensions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Machine learning moment closure models for the radiative transfer equation I: Directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for the radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. Finally, a variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

Computer science↗
DOE OSTI2020

Analytic Insights into the Neutronic Characteristics of Neutron Moderators from MCNP Calculations

In the course of re-processing the ENDF/B VIII.0 thermal scattering law data, formulas and data checking insights were developed from MCNP calculations into the neutronic characteristics of moderators. Infinite medium test calculations were carried out in MCNP as part of the verification effort for the ACE formatted thermal scattering files. Edit results from MCNP were combined with derived results to formulate several equations which give analytic insights into the neutronic behavior of moderators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The Numerical Computation of Clear-Sky Transmittance and Reflectance of Solar Radiation for the Modeling of PV System Performance

Modeling of photovoltaics (PV) system performance requires numerical computation of plane-of-array (POA) irradiances that are influenced by atmospheric absorption and scattering. Clear-sky transmittance and reflectance of diffuse irradiance are essential for the understanding of solar radiation related to multiple reflection between land surface and clouds. This study analytically solves the radiative transfer equation for clear-sky conditions and numerically computes the transmittance and reflectance of diffuse irradiance as functions of the optical thickness of the atmosphere. The computation will be implemented in all-sky radiative transfer models to improve their computational efficiency.

14 SOLAR ENERGY↗
DOE OSTI2025

QCD–Gravity Double Copy in Regge Asymptotics: From \(2\rightarrow n\) Amplitudes to Radiation in Shockwave Collisions

This paper discusses multi-particle production in QCD and in gravity at ultrarelativistic energies, their double-copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double-copy relations; in gravity, Weinberg’s soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. Renormalization group equations in the Color Glass Condensate (CGC) EFT describe wee parton correlations and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron–hadron collisions follows, with multi-particle production described by Cutkosky’s rules in strong time-dependent fields. Gluon radiation in the CGC EFT has a double copy in gravitational shockwave collisions, with a similar correspondence applicable between gluon and graviton shockwave propagators. Possible extensions of this semi-classical double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
DOE OSTI2023

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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)↗
DOE OSTI2022

General Relativistic Implicit Monte Carlo Radiation-hydrodynamics

Abstract We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ : an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation–matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Understanding parton evolution in matter from renormalization group analysis

We perform a renormalization group (RG) analysis of collinear hadron production in deep inelastic scattering on nuclei. We consider the limit where the parent parton energy E is large, while the medium opacity remains small. We identify the fixed order and leading enhanced medium contributions to the semi-inclusive cross sections and derive RG equations that resum multiple emissions near the endpoints of the splitting functions at first order in opacity. These evolution equations treat the same type of radiation enhancement in matter as the modified Dokshitzer-Gribov-Lipatov-Altarelli-Parisi approach, but differ in the way one regulates the collinear divergences. They provide a unique analytic insight into the problem of resummation and a faster and more efficient path to phenomenology. The new RG evolution framework is applied to study fragmentation in eA reactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Neutron Next-Event Estimators Kinematics (Rev.2)

This paper reviews the kinematics of neutron elastic and inelastic scattering with moving and stationary targets for contributions to Neutron Next-Event Estimators (NEEs). NEEs are often used for simulating the response of detectors in locations that have few random walk particles. Contributions to NEEs from collisions differ from sampling the outgoing particle state in standard random walk collisions as the location of the estimator is fixed, thus fixing the scattering angle between the incoming and outgoing directions. To calculate contributions to NEEs, the outgoing energies and probabilities of scatter toward the estimator must be calculated. This paper presents the most general case that encompasses elastic and inelastic scattering with both moving and stationary targets. The moving target case is required for elastic scattering with thermal motion due to the free gas thermal treatment approximation. The stationary target equations are presented for comparison with equations presented in canonical Monte Carlo texts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Implementation of TSL Evaluations Beyond the Incoherent Approximation [Slides]

This presentation discusses the motivation behind the project, which is to facilitate the production of TSL libraries that contain high fidelity physics including coherent inelastic effects, such as graphite and beryllium. High fidelity TSL libraries are of interest in various applications, for example in relation to neutron scattering instrument response analysis and directional biases (reflection) in neutronic analysis. In summation, the removal of typical approximations, in terms of incoherent, cubic, and atom site are examined. The application of generalized equations to crystalline graphite and beryllium TSL evaluations and the validation of data and data formats is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗