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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

(U) Average Scattering Angle Cosine Sensitivities

The Evaluated Nuclear Data File (ENDF) provides covariances of the average scattering angle cosine $\overline{μ}$ for 117 nuclides, and the NJOY nuclear data processing code can process these covariances. To propagate them to the uncertainty of a response due to uncertainties in $\overline{μ}$ requires sensitivities of the response with respect to $\overline{μ}$. This report uses NJOY’s definition of $\overline{μ}$ and derives a relationship between the sensitivity of a response to $\overline{μ}$ and the sensitivity of the response to the anisotropic P 1 elastic scattering cross section. This relationship has been implemented in the multigroup neutron sensitivity code SENSMG, which relies on the discrete ordinates code PARTISN for the neutron transport calculations. Kodeli has also calculated sensitivities with respect to $\overline{μ}$ and associated uncertainties. However, the formal derivation of this paper shows that the definition of the sensitivity with respect to $\overline{μ}$ is ambiguous because there are multiple ways to change $\overline{μ}$ by changing scattering cross sections. Here we derive two such ways. Both are implemented in SENSMG. This report is organized as follows. Section II discusses $\overline{μ}$ and its significance. Section III presents the equations used for the sensitivity of a response with respect to $\overline{μ}$, and Section IV discusses how to compute these sensitivities with SENSMG and MCNP6.3. Section V is a fast-spectrum test problem, and Sec. VI is a thermal-spectrum problem. Section VII is a summary. The example problem inputs are listed in the appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Analytic Sensitivity Coefficients for General Multigroup Infinite Medium k-Eigenvalue Problems

The general multigroup infinite medium k-eigenvalue neutron transport equation is used to derive analytic expressions for the infinite medium k-eigenvalue, the scalar neutron flux and adjoint, and the sensitivity of $k$ ∞ to perturbations in the multigroup nuclear data of a single species isotropic elastic scattering material. In the appendix, we present the multigroup nuclear data for U-235 and U-238 along with the corresponding k-eigenvalue, flux, adjoint, and sensitivity profiles, which include the sensitivity of $k$ ∞ to the total, fission, capture, and scattering macroscopic cross sections as well as to the group-to-group scattering cross section matrix, group neutron production, and the unconstrained and constrained fission neutron energy distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Analytic Sensitivity Coefficients for General Multigroup Infinite Medium k-Eigenvalue Problems

This work presents a general set of equations that can be used to rapidly generate new benchmarks to verify nuclear data sensitivity calculations. The general multigroup infinite medium k-eigenvalue neutron transport equation is used to derive analytic expressions for the infinite medium k-eigenvalue, the scalar neutron flux and adjoint flux, and the sensitivity of k∞ to perturbations in the multigroup nuclear data of a single species, isotropic and elastic scattering, material. The multigroup nuclear data for U-235 and U-238 is presented along with their corresponding k-eigenvalues, forward flux, adjoint flux, and sensitivity profiles, which include the sensitivity of k∞ to the total, fission, capture, and scattering macroscopic cross sections as well as to the group-to-group scattering cross section matrix, group-wise fission neutron production, and the unconstrained and constrained fission neutron energy distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quark and gluon helicity evolution at small x: revised and updated

We revisit the problem of small Bjorken-x evolution of the gluon and flavor-singlet quark helicity distributions in the shock wave (s-channel) formalism. Earlier works on the subject in the same framework resulted in an evolution equation for the gluon field-strength F 12 and quark “axial current” $\overline{\psi}\gamma$ + γ 5 ψ operators (sandwiched between the appropriate light-cone Wilson lines) in the double-logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant). In this work, we observe that an important mixing of the above operators with another gluon operator, ${}_D{}^{\leftarrow i} {}_D{}^{i}$, also sandwiched between the light-cone Wilson lines (with the repeated transverse index i = 1, 2 summed over), was missing in the previous works. This operator has the physical meaning of the sub-eikonal (covariant) phase: its contribution to helicity evolution is shown to be proportional to another sub-eikonal operator, ${}_D{}^{i} {}_D{}^{\leftarrow i}$, which is related to the Jaffe-Manohar polarized gluon distribution. In this work we include this new operator into small-x helicity evolution, and construct novel evolution equations mixing all three operators (${}_D{}^{i} {}_D{}^{\leftarrow i}$, F 12 , and $\overline{\psi}\gamma$ + γ 5 ψ), generalizing the results of . We also construct closed double-logarithmic evolution equations in the large-N c and large-N c &N f limits, with N c and N f the numbers of quark colors and flavors, respectively. Solving the large-N c equations numerically we obtain the following small-x asymptotics of the quark and gluon helicity distributions ΔΣ and ΔG, along with the g 1 structure function, $$ \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {g}_1\left(x,{Q}^2\right)\sim {\left(\frac{1}{x}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2\pi }}} $$ in complete agreement with the earlier work by Bartels, Ermolaev and Ryskin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Cross section of neutrons from the H 2 ( n , 2 n ) reaction at E n = 15 MeV

In this work, the double-differential cross section of the deuteron breakup reaction 2 H(n, 2n) has been studied experimentally with a neutron beam energy of 15 MeV. Special attention has been devoted to estimation of background condition and multiple scattering effect in the scattering sample. Experimental data have been compared with models based on phase-space approximation used in the ENDF/B-VIII.0 data library and in the MCNP neutron transport code, as well as with rigorous model based on Faddeev equations used for cross section evaluations in JENDL data library. It was found that experimental data are better reproduced by Faddeev model, however, the model overestimates data in the low-energy region of the neutron spectrum (<4 MeV)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A combined theoretical-experimental investigation of thermal transport in low-dose irradiated thorium dioxide

During reactor operation, nuclear fuels are subject to extreme temperature and irradiation conditions which can significantly degrade the fuel's thermal transport properties. The reduction in thermal conductivity of the fuel as a result of irradiation-induced lattice defects is arguably the most important fuel performance metric in regard to reactor efficiency and safety. Because thorium dioxide (ThO2) is suitable as a model system for more complex materials such as UO2 and its mixed oxides, we present a theoretical investigation of thermal conductivity reduction seen in defect-bearing thorium dioxide and compare directly to experimental measurements. Phonon-mediated thermal transport of the fuel is modeled by a solution to the Boltzmann transport equation (BTE) for phonons. A cluster dynamics (CD) model for lattice defect evolution during irradiation predicts defect densities which are used as input to the BTE for modeling phonon-defect scatterings. Phonon scatterings by lattice defects include those from point defects and vacancy clusters and interstitial clusters of various sizes. The CD model is benchmarked against structural defect characterization of irradiated thorium dioxide using electron microscopy. Thermal conductivity predicted by the BTE model is compared to measured values for irradiated thorium dioxide specimens below room temperature to isolate effects of phonon-defect scattering from intrinsic 3-phonon processes, which dominate at higher temperatures. The computed conductivity values are in partial agreement at temperatures close to room temperature while slight deviations are observed at the lowest measured temperatures, suggesting that implemented phonon-defect scattering cross-section expressions may not be adequate for low temperatures. The presented work provides a necessary investigation of the influence of irradiation induced defects on fuel performance and represents a first step toward a full characterization of phonon mediated thermal transport in irradiated materials with complex defect microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nuclear fragmentation at the future electron-ion collider

We explore the potential of conducting low-energy nuclear physics studies, including nuclear structure and decay, at the future Electron-Ion Collider (EIC) at Brookhaven. By comparing the standard theory of electron-nucleus scattering with the equivalent photon method applied to Ultraperipheral Collisions (UPC) at the Large Hadron Collider (LHC) at CERN. In the limit of extremely high beam energies and small energy transfers, very transparent equations emerge. We apply these equations to analyze nuclear fragmentation in UPCs at the LHC and 𝑒𝐴 scattering at the EIC, demonstrating that the EIC could facilitate unique photonuclear physics studies. Furthermore, we have also shown that the fragmentation cross-sections at the EIC are about 1,000 times smaller than those at the LHC. At the LHC, the fragmentation of uranium nuclei displays characteristic double-hump mass distributions from fission events, while at the EIC, fragmentation is dominated by neutron emission and fewer few fission products, about 10,000 smaller number of events.

Bertulani, Carlos A. [East Texas A&M University, C↗
DOE OSTI2020

Spin-dependent Pomeron and Odderon in elastic proton–proton scattering

We introduce a new model of near-forward elastic proton-(anti)proton scattering at high energy based on the modern formulation of Pomeron and Odderon in terms of Wilson lines and generalized TMDs (GTMDs). We compute the helicity-dependent elastic amplitudes $\phi$1,2,3,4,5 in this model and study their energy dependence from the nonlinear small-x evolution equations. While both Pomeron and Odderon contribute to helicity-flip processes in general, in the forward limit t = 0 only the double helicity-flip amplitude $\phi$2, dominated by the spin-dependent Odderon, survives. This may effect the extraction of the ρ parameter as well as the total cross section in the LHC energy domain and beyond

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗
DOE OSTI2021

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗
DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

QCD evolution of entanglement entropy

Entanglement entropy has emerged as a novel tool for probing nonperturbative quantum chromodynamics (QCD) phenomena, such as color confinement in protons. While recent studies have demonstrated its significant capability in describing hadron production in deep inelastic scatterings, the QCD evolution of entanglement entropy remains unexplored. Here, in this work, we investigate the differential rapidity-dependent entanglement entropy within the proton and its connection to final-state hadrons, aiming to elucidate its QCD evolution. Our analysis reveals a strong agreement between the rapidity dependence of von Neumann entropy, obtained from QCD evolution equations, and the corresponding experimental data on hadron entropy. These findings provide compelling evidence for the emergence of a maximally entangled state, offering new insights into the nonperturbative structure of protons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗
DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗