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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Machine learning on neutron and x-ray scattering and spectroscopies

Neutron and x-ray scattering represent two classes of state-of-the-art materials characterization techniques that measure materials structural and dynamical properties with high precision. These techniques play critical roles in understanding a wide variety of materials systems from catalysts to polymers, nanomaterials to macromolecules, and energy materials to quantum materials. In recent years, neutron and x-ray scattering have received a significant boost due to the development and increased application of machine learning to materials problems. This article reviews the recent progress in applying machine learning techniques to augment various neutron and x-ray techniques, including neutron scattering, x-ray absorption, x-ray scattering, and photoemission. We highlight the integration of machine learning methods into the typical workflow of scattering experiments, focusing on problems that challenge traditional analysis approaches but are addressable through machine learning, including leveraging the knowledge of simple materials to model more complicated systems, learning with limited data or incomplete labels, identifying meaningful spectra and materials representations, mitigating spectral noise, and others. We present an outlook on a few emerging roles machine learning may play in broad types of scattering and spectroscopic problems in the foreseeable future.

Chen, Zhantao↗

DOE OSTI2024

ROSE: A reduced-order scattering emulator for optical models

A new generation of phenomenological optical potentials requires robust calibration and uncertainty quantification, motivating the use of Bayesian statistical methods. These Bayesian methods usually require calculating observables for thousands or even millions of parameter sets, making fast and accurate emulators highly desirable or even essential. Emulating scattering across different energies or with interactions such as optical potentials is challenging because of the nonaffine parameter dependence, meaning the parameters do not all factorize from individual operators. Here we introduce and demonstrate the reduced-order scattering emulator ( ROSE ) framework, a reduced basis emulator that can handle nonaffine problems. ROSE is fully extensible and works within the publicly available band framework software suite for calibration, model mixing, and experimental design. As a demonstration problem, we use ROSE to calibrate a realistic nucleon-target scattering model through the calculation of elastic cross sections. This problem shows the practical value of the ROSE framework for Bayesian uncertainty quantification with controlled trade-offs between emulator speed and accuracy as compared to high-fidelity solvers. Finally, planned extensions of ROSE are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Fabrication of a Monolithic 5-Meter Aluminum Reflector for Millimeter-Wavelength Observations of the Cosmic Microwave Background

We have demonstrated the fabrication of a monolithic, 5-meter diameter, aluminum reflector with 17.4 $\mu$m RMS surface error. The reflector was designed to avoid the problem of pickup due to scattering from panel gaps in a large, millimeter-wavelength telescope that will be used for measurements on the cosmic microwave background.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Determining population densities in bimodal micellar solutions using contrast-variation small angle neutron scattering

Self-assembly of amphiphilic polymers in water is of fundamental and practical importance. Significant amounts of free unimers and associated micellar aggregates often coexist over a wide range of phase regions. The thermodynamic and kinetic properties of the microphase separation are closely related to the relative population density of unimers and micelles. Although the scattering technique has been employed to identify the structure of micellar aggregates as well as their time-evolution, the determination of the population ratio of micelles to unimers remains a challenging problem due to their difference in scattering power. In this work, using small-angle neutron scattering (SANS), we present a comprehensive structural study of amphiphilic n-dodecyl-PNIPAm polymers, which shows a bimodal size distribution in water. By adjusting the deuterium/hydrogen ratio of water, the intra-micellar polymer and water distributions are obtained from the SANS spectra. The micellar size and number density are further determined, and the population densities of micelles and unimers are calculated to quantitatively address the degree of micellization at different temperatures. Our method can be used to provide an in-depth insight into the solution properties of microphase separation, which are present in many amphiphilic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A hybrid Monte Carlo, discontinuous Galerkin method for linear kinetic transport equations

Here we present a hybrid method for time-dependent particle transport problems that combines Monte Carlo (MC) estimation with deterministic solutions based on discrete ordinates. For spatial discretizations, the MC algorithm computes a piecewise constant solution and the discrete ordinates use bilinear discontinuous finite elements. From the hybridization of the problem, the resulting problem solved by Monte Carlo is scattering free, resulting in a simple, efficient solution procedure. Between time steps, we use a projection approach to “relabel” collided particles as uncollided particles. In conclusion, from a series of standard 2-D Cartesian test problems we observe that our hybrid method has improved accuracy and reduction in computational complexity of approximately an order of magnitude relative to standard discrete ordinates solutions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

The PRad windowless gas flow target

In this work, we report on a windowless, high-density, gas flow target at Jefferson Lab that was used to measure $r_p$, the root-mean-square charge radius of the proton. To our knowledge, this is the first such system used in a fixed-target experiment at a (non-storage ring) electron accelerator. The target achieved its design goal of an areal density of 2$\times$10$^{18}$ atoms/cm$^2$, with the gas uniformly distributed over the 4 cm length of the cell and less than 1% residual gas outside the cell. This design eliminated scattering from the end caps of the target cell, a problem endemic to previous measurements of the proton charge radius in electron scattering experiments, and permitted a precise, model-independent extraction of $r_p$ by reaching unprecedentedly low values of $Q^2$, the square of the electron's transfer of four-momentum to the proton.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Localization and coherent imaging of hidden moving objects using laser speckle

Imaging and sensing of moving objects through opaque scattering media is a challenging but important problem in a variety of applications, including environmental sensing, biomedical imaging, and material inspection. We have previously demonstrated a technique to coherently image a moving object through thick, heavily scattering random media using correlations of speckle images as a function of the object’s spatial translation. Here, we demonstrate that this technique can be combined with localization to achieve imaging without prior knowledge of the object’s motion, greatly extending the application domain. This method is effective beyond the thin or weakly scattering regime and, rather than motion being deleterious, exploits the information available when the hidden object is moving, as could be the case in a cluttered terrestrial environment or through substantial levels of biological tissue scatter.

Hastings, Ryan L. (ORCID:0009000095977807)↗

DOE OSTI2021

Simulating low-energy neutrino interactions with MARLEY

Monte Carlo event generators are a critical tool for the interpretation of data obtained by neutrino experiments. Several modern event generators are available which are well-suited to the GeV energy scale used in studies of accelerator neutrinos. However, theoretical modeling differences make their immediate application to lower energies difficult. In this paper, I present a new event generator, MARLEY, which is designed to better address the simulation needs of the low-energy (tens of MeV and below) neutrino community. The code is written in C++14 with an optional interface to the popular ROOT data analysis framework. The current release of MARLEY (version 1.2.0) emphasizes simulations of the reaction 40 Ar ( ν e , e − ) 40 K ⁎ but is extensible to other channels with suitable user input. This paper provides detailed documentation of MARLEY's implementation and usage, including guidance on how generated events may be analyzed and how MARLEY may be interfaced with external codes such as Geant4. Further information about MARLEY is available on the official website at http://www.marleygen.org. Program title:MARLEY 1.2.0 CPC Library link to program files:https://doi.org/10.17632/4v7zxnc8j3.1 Developer's respository link:http://github.com/MARLEY-MC/marley Code Ocean capsule:https://codeocean.com/capsule/9868179 Licensing provisions: GNU General Public License 3.0 Programming language: C++14 External routines/libraries used: GNU Scientific Library [1,2] (required), ROOT [3,4] (optional) Nature of problem: Simulation of neutrino-nucleus scattering events at energies of tens-of-MeV and below Solution method: Initial two-to-two scattering kinematics are sampled using the allowed approximation differential cross section and tables of precomputed nuclear matrix elements. Subsequent de-excitations of the remnant nucleus are simulated using a Monte Carlo implementation of the Hauser-Feshbach statistical model and tabulated γ-ray decay schemes for discrete nuclear levels. Additional comments including restrictions and unusual features: Input data are provided with the code that are suitable for producing simulations of the charged-current reaction 40 Ar ( ν e , e − ) 40 K ⁎ , coherent elastic neutrino-nucleus scattering on spin-zero target nuclei, and neutrino-electron elastic scattering on any atomic target. Preparation of new reaction input files (whose format is documented in Appendix B) would enable other reaction channels and nuclear targets to be handled by the existing code framework. Although there is no maximum neutrino energy enforced by the code itself, realistic neutrino-nucleus scattering events may be generated up to roughly 50 MeV. Above this energy, the effects of forbidden nuclear transitions, which are neglected in the current treatment of the cross sections (see section 2.1), become increasingly important.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

X-ray multiphoton-induced Coulomb explosion images complex single molecules

Following structural dynamics in real time is a fundamental goal towards a better understanding of chemical reactions. Recording snapshots of individual molecules with ultrashort exposure times is a key ingredient towards this goal, as atoms move on femtosecond (10 -15 s) timescales. For condensed-phase samples, ultrafast, atomically resolved structure determination has been demonstrated using X-ray and electron diffraction. Pioneering experiments have also started addressing gaseous samples. However, they face the problem of low target densities, low scattering cross sections and random spatial orientation of the molecules. Therefore, obtaining images of entire, isolated molecules capturing all constituents, including hydrogen atoms, remains challenging. Here we demonstrate that intense femtosecond pulses from an X-ray free-electron laser trigger rapid and complete Coulomb explosions of 2-iodopyridine and 2-iodopyrazine molecules. We obtain intriguingly clear momentum images depicting ten or eleven atoms, including all the hydrogens, and thus overcome a so-far impregnable barrier for complete Coulomb explosion imaging—its limitation on molecules consisting of three to five atoms. In combination with state-of-the-art multi-coincidence techniques and elaborate theoretical modelling, this allows tracing ultrafast hydrogen emission and obtaining information on the result of intramolecular electron rearrangement. Our work represents an important step towards imaging femtosecond chemistry via Coulomb explosion.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Piezoelectric Electron-Phonon Interaction from Ab Initio Dynamical Quadrupoles: Impact on Charge Transport in Wurtzite GaN

First-principles calculations of $\textit{e}$-ph interactions are becoming a pillar of electronic structure theory. However, the current approach is incomplete. The piezoelectric (PE) $\textit{e}$-ph interaction, a long-range scattering mechanism due to acoustic phonons in noncentrosymmetric polar materials, is not accurately described at present. Current calculations include short-range $\textit{e}$-ph interactions (obtained by interpolation) and the dipolelike Frölich long-range coupling in polar materials, but lack important quadrupole effects for acoustic modes and PE materials. Here we derive and compute the long-range $\textit{e}$-ph interaction due to dynamical quadrupoles, and apply this framework to investigate $\textit{e}$-ph interactions and the carrier mobility in the PE material wurtzite GaN. In this work, we show that the quadrupole contribution is essential to obtain accurate $\textit{e}$-ph matrix elements for acoustic modes and to compute PE scattering. Our work resolves the outstanding problem of correctly computing $\textit{e}$-ph interactions for acoustic modes from first principles, and enables studies of $\textit{e}$-ph coupling and charge transport in PE materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

A global analysis strategy to resolve neutrino NSI degeneracies with scattering and oscillation data

Neutrino non-standard interactions (NSI) with the first generation of standard model fermions can span a parameter space of large dimension and exhibit degeneracies that cannot be broken by a single class of experiment. Oscillation experiments, together with neutrino scattering experiments, can merge their observations into a highly informational dataset to combat this problem. We consider combining neutrino-electron and neutrino-nucleus scattering data from the Borexino and COHERENT experiments, including a projection for the upcoming coherent neutrino scattering measurement at the CENNS-10 liquid argon detector. We extend the reach of these data sets over the NSI parameter space with projections for neutrino scattering at a future multi-ton scale dark matter detector and future oscillation measurements from atmospheric neutrinos at the Deep Underground Neutrino Experiment (DUNE). In order to perform this global anal- ysis, we adopt a novel approach using the copula method, utilized to combine posterior information from different experiments with a large, generalized set of NSI parameters. We find that the contributions from DUNE and a dark matter detector to the Borexino and COHERENT fits can improve constraints on the electron and quark NSI parameters by up to a factor of 2 to 3, even when relatively many NSI parameters are left free to vary in the analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Desmearing Bonse–Hart USANS data using Bayesian Gaussian process regression

Ultra-small-angle neutron scattering (USANS) enables access to micrometer-scale structures but is intrinsically affected by strong, anisotropic resolution smearing arising from slit-geometry optics. As a result, recovery of the intrinsic scattering intensity constitutes an ill-posed inverse problem, and commonly used iterative desmearing methods lack rigorous uncertainty quantification. We present a Bayesian desmearing framework for slit-geometry USANS based on Gaussian process regression. In this approach, the scattering intensity is modeled as a smooth random function, and the instrumental point spread function is incorporated explicitly as a forward operator. The resulting formulation yields a closed-form maximum a posteriori solution with well-defined credibility intervals. Computational benchmarks and experimental validation using combined USANS and small-angle neutron scattering (SANS) measurements demonstrate that the framework enables stable desmearing, suppresses experimental noise, and preserves physically meaningful structural features under realistic conditions.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

DOE OSTI2020

Comparison of Generation of Higher-Order Neutron Scattering Cross Sections

The generation of high-order neutron scattering cross sections consistent with high-fidelity simulations remains an area of active research. Popular options include generating cross sections from continuous energy Monte Carlo calculations or generating cross sections from a deterministic neutron transport calculation with high-fidelity tabulated cross sections. Both options present challenges. Monte Carlo simulations can naturally process continuous energy cross section data and allow for the general description of anisotropic neutron scattering given a scattering law. However, Monte Carlo simulations are inherently related to particle weighting and it has been suggested that this may be unacceptable for generating high-order neutron scattering cross sections. Deterministic neutron transport calculations can easily calculate high-order moments of the flux to appropriately calculate higher-order neutron scattering cross sections but are generally limited by discretization of space, energy, and angle. In this work, the trade-offs between generation of high-order neutron scattering cross sections via Monte Carlo and deterministic neutron transport methods are investigated. The methods implemented in the Monte Carlo computer program Serpent and the deterministic fast reactor neutron cross section generator MC2-3 are compared. Cross sections resulting from these methods are used in Rattlesnake, a deterministic neutron transport code developed by Idaho National Laboratory (INL), and results are compared to a reference continuous energy Monte Carlo calculation. Whereas previous work investigating the effects of anisotropic neutron scattering has focused on light water reactor (LWR) simulations, this work focuses on high-order neutron scattering cross sections as they relate to fast reactor simulations. To investigate the consequences of the Serpent and MC2-3 methodologies, a test problem is developed. The test problem is a one-dimensional geometry with fast reactor materials designed to demonstrate deep penetration and exacerbate the effects of high-order neutron scattering. Based on the results of the deep penetration test problem, it is concluded that P3 neutron scattering cross sections are sufficient to describe anisotropic scattering in fast reactor materials. Any neutron scattering of order higher than P3 offers negligible change in the eigenvalue of the test problem. Additionally, it is determined that the methodology as implemented in Serpent is applicable for generating high-order neutron scattering cross sections through at least P3 in fast reactor materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Relations between anomalous dimensions in the Regge limit

We extend the recent formalism developed for computing rapidity anomalous dimension of form factors using unitarity to the problem of high-energy near forward scattering. By combining the factorization of 2 → 2 scattering in the effective field theory (EFT) for Glauber operators with definite signature amplitudes, we derive an expression that relates anomalous dimensions (including Regge trajectories) to cut amplitudes, leading to significant computational simplifications. We demonstrate this explicitly by computing the one and two-loop Regge trajectories. Our formalism can also be used to bootstrap anomalous dimensions of operators not related by symmetries. As an example, we show that the full anomalous dimensions (including both the Regge pole and cut pieces) of the two Glauber exchange anti-symmetric octet operator, can be determined from the anomalous dimension of the single Glauber exchange operator. Many other such relations exist between other color channels at each order in α.

Effective Field Theories↗

DOE OSTI2021

Dark matter haloes in the multicomponent model. III. From dwarfs to galaxy clusters

A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a ‘velocity kick’ parameter. The simplest 2cDM model with cross-section 0.01 ≲ σ/m < 1cm 2 g -1 and the kick velocity V k ≃100kms -1 have been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in N-body cosmological simulations tested on Milky Way (MW)-like haloes of a virial mass ~5 × 10 11 M ⊙ (Papers I & II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of ~10 7 -10 8 and ~10 13 -10 14 M ⊙ ⁠, respectively. We find that σ 0 /m ≳ 0.1 cm 2 g -1 is preferentially disfavoured for both dwarfs and galaxy cluster haloes in comparison with observations, while σ 0 /m = 0.001 cm 2 g -1 causes little perceptible difference from that of the CDM counterpart for most of the cross-section’s velocity dependence studied in this work. Our main result is that within the reasonable set of parameters, the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes, and the model leaves us an open window for other possible alternative DM models.

79 ASTRONOMY AND ASTROPHYSICS↗