Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scattering problem”

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

Coarse Mesh Finite Difference Acceleration for Pebble Tracking Transport in Griffin

We implemented a coarse mesh finite difference (CMFD) for accelerating transport calculations with PTT (pebble tracking transport) in the Griffin code. More specifically, extensions for transport update with the consideration of scattering operator and CMFD projection were implemented for PTT. The implementation was verified with a simplified PBR (pebble bed reactor) benchmark problem and significant performance improvements in CPU time was observed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Tomographic Reconstruction of the Neutron Time-Energy Spectrum from a Dense Plasma Focus

A dense plasma focus (DPF) nuclear fusion device is an attractive pulsed neutron source in many applications, due to its relatively large neutron yield, produced in a short time duration. In order to design a DPF that generates neutrons within a specified time profile, or generates a neutron energy spectrum with specific properties, it is necessary to be able to characterize and model the results of the fusion process releasing the neutrons. The time-energy spectrum of fusion neutrons is an ideal quantity to use to validate multiphysics codes that simulate the pinch and fusion processes, because it is a quantity that requires the physics of each of the stages leading up to and ending in the fusion reaction to be simulated correctly. In particular, since DPF fusion neutrons are not monoenergetic -- and there can often be several fusion pinches creating neutrons -- a computer simulation matching high quality neutron spectrum measurements provides great confidence in the fidelity of the simulation. In order to make such a comparison, it is first necessary to have quality measurements from which to infer the spectrum. In this work we pose neutron spectroscopy as the classical tomographic inverse problem from neutron time of flight data at multiple distances, but enhanced by using an additional measurement of the time profile of the fusion pinch near its source and by using detector pairs set up in a geometry that allows for scatter background subtraction. The detector pairs enhance the quality of the time of flight measurements, and the additional constraint posed by the measured time profile allows for reconstructions discretized as finely as the time measurements and in energy as finely as 100 keV, without the problem being underdetermined. We present results from a Deuterium-fueled DPF at the U.S. Department of Energy's Nevada National Security Site and show that we can infer the time-energy spectrum from our measurements for both single and multi-pinch fusion reactions with equal fidelity.

Catenacci, Jared↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nonparaxial propagation of an intense relativistic electron beam through dense media

An envelope equation is used to explore the propagation of a focused electron beam through dense conductive media like partially ionized high-pressure gas, metal vapor plumes from beam-target interactions, or low-density solid-phase targets. Envelope equations have been a most useful tool for predicting experimental results of beam propagation in gas, and they are refined for the general problem of propagation through a dense medium. The envelope equation is modified to account for the large beam divergence resulting from multiple scattering. The envelope code results are in qualitative agreement with published experimental data. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Entanglement Witness for Indistinguishable Electrons Using Solid-State Spectroscopy

Characterizing entanglement in quantum materials is crucial for advancing next-generation quantum technologies. Despite recent strides in witnessing entanglement in magnetic materials with distinguishable spin modes, quantifying entanglement in systems formed by indistinguishable electrons remains a formidable challenge. To solve this problem, we introduce a method to extract various four-fermion correlations by analyzing the nonlinearity in resonant inelastic x-ray scattering spectra. These correlations constitute the primary components of the cumulant two-particle reduced density matrix. We further derive bounds for its eigenvalues and demonstrate the linear scaling with fermionic entanglement depth, providing a reliable witness for entanglement. Using the material-relevant strongly correlated models as examples, we show how this entanglement witness can efficiently quantify multipartite entanglement across different phase regions, highlighting its advantage over quantum Fisher information. Published by the American Physical Society 2025

Liu, Tongtong (ORCID:0000000295324061)↗
DOE OSTI2020

Scattering Signatures of Bond-Dependent Magnetic Interactions

Bond-dependent magnetic interactions can generate exotic phases such as Kitaev spin-liquid states. Experimentally determining the values of bond-dependent interactions is a challenging but crucial problem. In this paper, I show that each symmetry-allowed nearest-neighbor interaction on triangular and honeycomb lattices has a distinct signature in paramagnetic neutron-diffraction data, and that such data contain sufficient information to determine the spin Hamiltonian unambiguously via unconstrained fits. Moreover, I show that bond-dependent interactions can often be extracted from powder-averaged data. These results facilitate experimental determination of spin Hamiltonians for materials that do not show conventional magnetic ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fourier analysis of continuous fractional diffusion synthetic acceleration schemes in slab geometry

We propose two fractional extensions of continuous diffusion synthetic acceleration (DSA) with fractional derivative order α varying over the interval 2 ≥ α ≥ 1 . We investigate the spectral properties of the corresponding continuous families of fractional preconditioners by performing Fourier analysis for a model infinite homogeneous medium problem in slab geometry. The first family results in a fractional acceleration scheme, FrDSAo, that reduces to traditional DSA for .α = 2 and scattering ratio c limiting to a unit value (c → 1) but is otherwise optimized via the Fourier analysis, to obtain the smallest possible spectral radius, for c < 1 and 2 ≥ α ≥ 1. The second family corresponds to a fractional acceleration scheme, FrDSAs, that reduces to traditional DSA for α = 2 for all values of c. The latter scheme is not optimized but has the advantage of lending itself to a more straightforward implementation. For high values of c, the results of the Fourier analysis point to the existence of an interval 2 > α > ∼1.8 where both FrDSAo and FrDSAs can achieve a lower spectral radius than DSA. For example, DSA has a spectral radius of ∼0.2246 for c = 0.9999 while FrDSAo produces a value of ∼0.1616 at α = 1.92 and FrDSAs results in ∼0.2116 at α =1.93. (author)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Qualitative Conditions for Reasonable Estimates of Kerma Using Total and Kerma Cross Section Attenuation of Photons

Performing a full adjoint simulation in the Integrated Tiger Series (ITS) can be time consuming to obtain statistically significant results. The ray-trace capability in ITS allows for rapid scoping calculations of the uncollided kerma from photon sources without needing to run a full adjoint simulation. However, the uncollided estimate will always underestimate the full-physics kerma since it neglects scattered radiation, and under certain conditions, the result from the capability may not provide a sufficiently accurate estimate of the full-physics kerma. To exemplify the conditions in which the capability provides reasonable estimates, two problem geometries with different materials are simulated using the full adjoint capability as well as the ray-trace capability with a total cross section treatment and a new kerma-attenuation cross-section treatment. The results are then compared to show under which conditions the estimates are accurate. Reasonable estimates are provided from the ray-trace feature when there is minimal scattering into a region of interest occurring, such as with low-Z material and when the photons travel through small amounts of material.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Rapid 3D nanoscale coherent imaging via physics-aware deep learning

Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programing framework trained to predict 3D structure and strain, solely from input 3D X-ray coherent scattering data. Our networks are designed to be “physics-aware” in multiple aspects; in that the physics of the X-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods. Our integrated machine learning and differential programing solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Tycho: Realistically Simulating Exoplanets within Stellar Clusters. I. Improving the Monte Carlo Approach

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses >0.5M {sub ⊙} (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale N-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, Tycho, which employs a variety of different computational techniques, including multiple N-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1% of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which that perturb long-period planets lead to 38.3% of solar-system-like planetary systems becoming long-term unstable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Using ultrasonic attenuation in cortical bone to infer distributions on pore size

Here, in this work we infer the underlying distribution on pore radius in human cortical bone samples using ultrasonic attenuation data. We first discuss how to formulate polydisperse attenuation models using a probabilistic approach and the Waterman Truell model for scattering attenuation. We then compare the Independent Scattering Approximation and the higher-order Waterman Truell models’ forward predictions for total attenuation in polydisperse samples. Following this, we formulate an inverse problem under the Prohorov Metric Framework coupled with variational regularization to stabilize this inverse problem. We then use experimental attenuation data taken from human cadaver samples and solve inverse problems resulting in nonparametric estimates of the probability density function on pore radius. We compare these estimates to the “true” microstructure of the bone samples determined via microCT imaging. We find that our methodology allows us to reliably estimate the underlying microstructure of the bone from attenuation data.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2025

Pulling Order Back from the Brink of Disorder: Observation of a Nodal-Line Spin Liquid and Fluctuation Stabilized Order in K 2 ⁢IrCl 6

Competing interactions in frustrated magnets can give rise to highly degenerate ground states from which correlated liquidlike states of matter often emerge. The scaling of this degeneracy influences the ultimate ground state, with extensive degeneracies potentially yielding quantum spin liquids, while subextensive or smaller degeneracies yield static orders. A long-standing problem is to understand how ordered states precipitate from this degenerate manifold and what echoes of the degeneracy survive ordering. Here, we use neutron scattering to experimentally demonstrate a new “nodal-line” spin liquid, where spins collectively fluctuate within a subextensive manifold spanning one-dimensional lines in reciprocal space. Realized in the spin-orbit-coupled, face-centered-cubic iridate K 2 ⁢IrCl 6 , we show that the subextensive degeneracy is robust, but remains susceptible to fluctuations or longer-range interactions which cooperate to select a magnetic order at low temperatures. Proximity to the nodal-line spin liquid influences the ordered state, enhancing the effects of quantum fluctuations that in turn act to stabilize the sublattice magnetization through the self-consistent opening of a large spin-wave gap. Our results demonstrate how quantum fluctuations can act counterintuitively in frustrated materials: Even in a case where fluctuations are ineffective at selecting an ordered state from a degenerate manifold, at the brink of the nodal spin liquid, they can act to protect the ordered state and dictate its low-energy physics.

36 MATERIALS SCIENCE↗
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 OSTI2022

Updraft dynamics and microphysics: on the added value of the cumulus thermal reference frame in simulations of aerosol–deep convection interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Multimessenger measurements of the static structure of shock-compressed liquid silicon at 100 GPa

The ionic structure of high-pressure, high-temperature fluids is a challenging theoretical problem with applications to planetary interiors and fusion capsules. Here we report a multimessenger platform using velocimetry and angularly and spectrally resolved x-ray scattering to measure the thermodynamic conditions and ion structure factor of materials at extreme pressures. We document the pressure, density, and temperature of shocked silicon near 100 GPa with uncertainties of 6%, 2%, and 20%, respectively. The measurements are sufficient to distinguish between and rule out some ion screening models. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Measuring Neutron Polarisation in Deuteron Photo-disintegration with the CLAS Start Counter [Thesis]

Deuteron photo-disintegration (γd → γp) is a reaction that represents the simplest case in which nuclear and hadron physics models can be tested. Despite this, associated polarization analyses are limited in terms of angular coverage and energy ranges, especially in observables related to the recoil neutron. This is largely due to a lack in dedicated polarimetry equipment, and represents a roadblock in global progress to understand high-energy phenomena such as hexaquarks, and quark-gluon degrees of freedom. To address this problem, this PhD thesis pioneers a new methodology for the parasitic measurement of nucleon polarization using kinematic reconstruction of (spin-dependent) nucleon-nucleus scattering of reaction products, prior to their detection in large acceptance particle detector apparatus. Following this novel approach, which requires no dedicated polarimeter, a determination of the double polarization observable, $C^n_{x'}$, from deuteron photo-disintegration is presented, using Jefferson Lab’s CLAS detector. The analysis utilizes the (n,p) charge exchange reaction in CLAS’s "start counter" (plastic scintillator) to determine the final state neutron polarizations. The results present the first ever data for this observable above 0.7 GeV (photon beam energy) and significantly extend the angular range of the world data set. This new data is largely statistically consistent with the previous measurement of $C^n_{x'}$ by Bashkanov et al . in the overlapping energy range of 0.4-0.7 GeV. It is planned for the statistical accuracy of the presented result to be increased by the inclusion of additional data. The analysis herein serves as a key proof of concept for future applications, including a recommended similar analysis to be implemented with data from the more modern CLAS12 detector. This paves the way for a plethora of additional analyses using existing data sets that would provide crucial new constraints for hadron and nuclear physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗