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

Aberration corrected RF flipper for high resolution neutron spectroscopy

Project Summary Company: Adelphi Technology, Inc. Title: Aberration-corrected High Frequency RF Flipper for High-Resolution Neutron Spectroscopy PI: Dr. Jay Theodore Cremer Topic: C55-11 Enhancement of Scattering Instrumentation Technology Used at Pulsed and Continuous Sources Subtopic: d. Other Statement of the problem or situation that is being addressed. The quest to understand heterogeneous and hierarchical materials is gathering momentum, as described in a 2015 report by the Basic Energy Sciences Advisory Committee on Challenges at the Frontiers of Matter and Energy. For the past 40 years a technique called neutron spin echo (NSE) has been used to probe molecular motions in such non-crystalline materials over time scales from 10’s of picoseconds to 100’s of nanoseconds. The method has provided unique information about the dynamics of soft heterogeneous materials, including confirmation of the de Gennes model of polymer reptation and quantitative measurement of bending constants of biologically relevant lipid membranes. However, scientists continue to clamor for even higher resolution than NSE can provide. Biomaterials, polymers, glasses, and artificially nanostructured materials all manifest slow molecular motions because of weak or competing interactions between subunits and are amenable to study with neutrons, provided sufficiently long dynamical correlation times can be achieved. All these materials have important applications to advanced technologies so understanding them is key to technological progress. General statement of how this problem is being addressed. We will address the need for high-resolution neutron spectroscopy by using a technique called Neutron Resonance Spin Echo (NRSE). While similar to NSE in many respects, this method has the potential to exceed the NSE capabilities, if 2 technical hurdles can be overcome. The major impediments to successful high-resolution NRSE are the availability of two technologies: a very high frequency, efficient, radiofrequency (rf) flipper for neutrons and a method to correct certain magnetic aberrations. Based on previous STTR support and follow-on research we have developed a suitable rf flipper and we have invented a method to correct the magnetic aberrations. Both technologies need refinement to make them suitable for implementation at a neutron source such as the Oak Ridge National Laboratory nuclear reactor. In this proposal we seek to perfect the two technologies and to combine them into a single, operationally convenient device. Commercial Applications and Other Benefits In view of the increasing demand for the unique scientific information that high resolution neutron spectroscopy can provide, we expect several major instrumentation upgrades at both U.S. and foreign neutron centers will require make use of the NRSE method over the coming decade, creating a market for the devices we will design. These components will enhance scientists’ abilities to probe the time dependence of density fluctuations in a wide range of hierarchical and heterogeneous materials many of which are vital to existing and future technologies. Key Words – Polarized Neutrons, Neutron Spin Echo, Neutron Scattering, advanced materials. Summary for Members of Congress Neutron beams are a powerful materials-science probe that provide unique information about the structure of matter. The proposed devices will accelerate scientific discoveries required to achieve national goals for new technological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Revisiting the K-edge X-ray absorption fine structure of Si, Ge–Si alloys, and the isoelectronic series: CuBr, ZnSe, GaAs, and Ge

Extended X-ray absorption fine structure (EXAFS) has evolved into an unprecedented local-structure technique that is routinely used to study materials’ problems in the biological, chemical, and physical sciences. Like many other experimental techniques, EXAFS also requires that several key atomic parameters must be known a priori before structural information can be quantitatively determined. Utilizing current analytical methods, we revisit the isoelectronic series CuBr, ZnSe, GaAs, and Ge originally studied by Stern et al. during the early development of EXAFS. We demonstrate that the ab initio EXAFS code FEFF accurately predicts the atomic phase shifts and backscattering amplitudes that are primarily functions of the sum of atomic numbers Z along an EXAFS scattering path. We also investigate quantitative fitting and first- and second-shell phase transferability together with problems that arise if a backscattering atom is identified incorrectly in an EXAFS fitting model. Features in the near-edge region, on the other hand, are shown to require a comprehensive treatment of the band structure and density-of-states, including effects of the screened Coulomb interaction between the photoelectron and core hole. Here, we demonstrate that the Bethe–Salpeter equation (BSE) accurately captures the NEXAFS (or XANES) portion of the spectrum for the isoelectronic series in addition to Si and Ge–Si alloys, including within a few eV of the absorption edge, where band structure and excitonic effects are most important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Bi-cross validation of spectral clustering hyperparameters

One challenge impeding the analysis of terabyte scale X-ray scattering data from the Linac Coherent Light Source (LCLS) is determining the number of clusters required for the execution of traditional clustering algorithms. Here, we demonstrate that the previous work using bi-cross validation to determine the number of singular vectors directly maps to the spectral clustering problem of estimating both the number of clusters and hyperparameter values. Here by applying this method to LCLS X-ray scattering data enables the identification of dropped shots without manually setting boundaries on detector fluence and provides a path toward identifying rare and anomalous events.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A New Galerkin Quadrature Method Not Requiring a Matrix Inverse

We derive a new Galerkin quadrature (GQ) method for S 𝑛 calculations that differs from the two methods preceding it in that a matrix inverse for an 𝑁 𝑑 × 𝑁 𝑑 matrix, where 𝑁𝑑 is the number of directions in the quadrature set, is no longer required. Galerkin quadrature methods are designed for calculations with highly anisotropic scattering. Such methods are not simply special angular quadratures but also are methods for representing the S 𝑛 scattering source that offers several advantages relative to the standard scattering source representation when highly truncated Legendre cross-section expansions must be used. Galerkin quadrature methods are also useful when the scattering is moderately anisotropic, but the quadrature being used is not sufficiently accurate for the order of the scattering source expansion that is required. Furthermore, we derive the new method and present computational results showing that its performance for two challenging problems is comparable to those of the two GQ methods that preceded it.

Galerkin quadrature↗
DOE OSTI2023

Material Identification From Radiographs Without Energy Resolution

We propose a method for performing material identification from radiographs without energy-resolved measurements. Material identification has a wide variety of applications, including in biomedical imaging, nondestructive testing, and security. While existing techniques for radiographic material identification make use of dual energy sources, energy-resolving detectors, or additional (e.g., neutron) measurements, such setups are not always practical— requiring additional hardware and complicating imaging. We tackle material identification without energy resolution, allowing standard X-ray systems to provide material identification information without requiring additional hardware. Assuming a setting where the geometry of each object in the scene is known and the materials come from a known set of possible materials, we pose the problem as a combinatorial optimization with a loss function that accounts for the presence of scatter and an unknown gain and propose a branch and bound algorithm to efficiently solve it. We present experiments on both synthetic data and real, experimental data with relevance to security applications— thick, dense objects imaged with MeV X-rays. We show that material identification can be efficient and accurate, for example, in a scene with three shells (two copper, one aluminum), our algorithm ran in six minutes on a consumer-level laptop and identified the correct materials as being among the top 10 best matches out of 8,000 possibilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three-Dimensional Imaging Lidar for Characterizing Particle Fields and Organisms in the Mesopelagic Zone

The ocean’s mesopelagic zone is largely uncharacterized despite its vital role in sustaining ocean ecosystems. The composition, cycling, and fate of particle fields in the mesopelagic lacks an integrative multi-scale understanding of organism migration patterns, distribution, and diversity. This problem is addressed by combining complementary technologies with overlapping size spectra, including profiler mounted optical scattering sensors, profiler, and ship mounted acoustic devices, and a custom Unobtrusive Multi-Static Lidar Imager (UMSLI). This unique sensor suite can observe distributions of particles including organisms over a six order of magnitude dynamic size range, from microns to meters. Overlapping size ranges between different methods allows for cross-validation. This work focuses on the lidar imaging measurements and optical backscattering and attenuation, covering a combined particle size range of 0.1 mm to several cm. Particles at the small end of this range are sized using an existing backscattering time series inversion method after Briggs et al. (2013). Larger particles are resolved with UMSLI over an expanding volume using three-dimensional photo-realistic laser serial imaging. UMSLI’s image rectifying ability over time allows for derivation of particle concentration, size, and spatial distribution. Technical details on the development and post-processing methods for the novel UMSLI system are provided. Image resolved particle size distributions (PSDs) revealed a size shift from smaller to larger particles (>0.5 mm) as indicated by flatter slopes from dawn (slope = 2.6) to dusk (slope = 3.0). PSD trends are supported by an optical backscatter and transmissometer time series inversion analysis. Size shifts in the particle field are largely attributed to aggregation effects. Images support evidence of temporal variation between dusk and dawn stations through statistical analysis of particle concentrations for particle sizes 0.50–5.41 mm. Spatial analysis of the particle field revealed a dominantly uniform distributed marine snow background. The importance and potential of integrated approaches to studying particle and organism dynamics in ocean environments are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evolution of Bose–Einstein condensate systems beyond the Gross–Pitaevskii equation

While many phenomena in cold atoms and other Bose–Einstein condensate (BEC) systems are often described using the mean-field approaches, understanding the kinetics of BECs requires the inclusion of particle scattering via the collision integral of the quantum Boltzmann equation. A rigorous approach for many problems in the dynamics of the BEC, such as the nucleation of the condensate or the decay of the persistent current, requires, in the presence of factors making a symmetry breaking possible, considering collisions with thermal atoms via the collision integral. These collisions permit the emergence of vorticity or other signatures of long-range order in the nucleation of the BEC or the transfer of angular momentum to thermal atoms in the decay of persistent current, due to corresponding terms in system Hamiltonians. Here, we also discuss the kinetics of spin–orbit-coupled BEC. The kinetic equation for the particle spin density matrix is derived. Numerical simulations demonstrate significant effects of the collision integral on the dynamics of the spin–orbit-coupled BEC upon quenching of the Raman coupling that generates synthetic electric and magnetic fields.

Bose–Einstein condensation↗
DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

The ABCs of phase retrieval: Connecting the acronyms of scanning transmission electron microscopy

High-resolution scanning transmission electron microscopy (S/TEM) is an indispensable tool for characterizing the structure and properties of materials down to the atomic scale. Conventional S/TEM imaging, however, is limited by the phase problem, whereby the phase of the electron exit wave is lost upon detection. Recent advances in diffractive imaging and 4D-STEM have enabled a range of phase-retrieval techniques that computationally reconstruct the missing information encoded in the phase of the transmission function. These approaches offer improved dose efficiency and enhanced sensitivity to weakly scattering signals, extending quantitative imaging to beam-sensitive materials composed of light elements. In this work, we introduce the phase problem in electron microscopy and survey the diverse landscape of phase-retrieval techniques used in the field. Despite their many acronyms and algorithmic variations, these techniques share a common physical and mathematical foundation. We present a unified framework that connects these seemingly distinct methods, from parallax imaging and tilt-corrected bright-field (tcBF-STEM), to aberration-corrected bright-field (acBF-STEM), optimum bright-field (OBF-STEM) and single-sideband (SSB) ptychography, as well as first-moment integrated center of mass techniques (iCOM) and iterative ptychographic algorithms. Based on these insights, we discuss the opportunities and practical limitations of applying these methods across different materials systems, detector designs, and microscope configurations.Graphical abstractRepresentative electron microscopy configurations used for phase retrieval and diffractive imaging in S/TEM: (a) Zernike phase-contrast transmission electron microscopy (TEM), (b) small-convergence-angle four-dimensional scanning transmission electron microscopy (4D-STEM) for nanobeam-based phase reconstruction methods, and (c) large-convergence-angle 4D-STEM for ptychographic and related diffractive imaging techniques reviewed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Laser scattering by submicron droplets formed during the electrical explosion of thin metal wires

Abstract This paper presents the results of studying dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators (~1–10 kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and shadow imaging at two wavelengths (1.064 µ m and 0.532 µ m). Our findings indicate the important role in the visualization of the explosion products that belongs to multiple scattering by submicron droplets of dense condensed matter, which are mixed with metal vapor. The hypothesis on the existence of submicron droplets in the products of exploding metal wires correlates with the results obtained by soft x-ray radiography combined with a laser probing technique. Taking into account the multiple scattering by submicron droplets, it is possible to significantly clarify the parameters of the explosion products visualized via laser probing techniques as well as to gain a deeper insight into the physics behind the electrical wire explosion.

Physics↗
DOE OSTI2024

Adaptive Methods for Radial Basis Functions

Radial basis functions (RBFs) are a powerful tool for constructing high-order accurate reduced representations of scattered data in arbitrary dimension and on manifolds. We present a method of constructing data approximations in which we utilize a functional tail to capture a global background profile and a RBF neural network (NN) to capture the smaller-scale features. In the RBF NN the RBF centers, matrix shape parameters were selected adaptively for each RBF. We also utilized a geodesic notion of distance on the manifold on which the data lies, e.g., the spherical geodesic for data on the sphere. Although each of these ideas have been been investigated separately in previous works, their combination into a single algorithm is novel. We defined a machine learning problem in which these properties are learned to minimize the data reduction error. We demonstrate the algorithm for applications of scattered data reduction in the plane and on the sphere.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Data Mining and Visualization of High-Dimensional ICME Data for Additive Manufacturing

Integrated computational materials engineering (ICME) methods combining CALPHAD with process-based simulations can produce rich, high-dimensional data for alloy and process design. In ICME methods for metallurgical applications, the visualization and interpretation of such high-dimensional data has previously been through heat maps represented in 2 or 3 dimensions. While such an approach is ideal when one variable is varied at a time, in the case of high-dimensional data with multiple variables varied simultaneously, as is the case in additive manufacturing, interpreting the trends through two- or three-dimensional heat maps becomes challenging. Here, we propose a strategy of mixed visual data mining and quantitative analysis for high-dimensional metallurgical and process data using high-throughput thermodynamic calculations. Two case studies show the application of the proposed approach. The first case study investigated the effects of feedstock chemistry on the δ ferrite formation in 316L stainless steel powders used for binder jet additive manufacturing. The second case study linked Scheil–Gulliver calculations to a process model for dissimilar joining of aluminum alloys 5356 and 6111 during laser hot-wire additive manufacturing. Both cases contained thousands of calculated data points, showcasing the utility of visual data analysis through parallel coordinate plotting, Pearson correlation coefficient matrices, and scatter matrices compared to traditional process maps. These visualization techniques can be extended to many additive manufacturing problems to capture process–structure–property relationships for additively manufactured components.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiphysics FISPACT-II and TENDL-2019 simulation: neutron-induced damage metrics

Nuclear interactions can be the source of atomic displacement, irradiation-induced defects and transmutation in structural materials. Such quantities are derived from, or can be correlated to, nuclear kinematic simulations of primary atomic energy distributions spectra and the quantification of the numbers of secondary defects produced per primary as a function of the available recoils, residual and emitted, energies. Recoil kinematics of neutral, residual, gas, proton, alpha particle emissions are now more rigorously treated based on recent, complete and enhanced nuclear data parsed in state-of-the-art processing tools. Defect production metrics are the starting point in the complex problem of correlating and simulating the behaviour of materials under irradiation, as experimental information is rare or scattered. Detailed, segregated primary knock-on-atom metrics are now becoming available as the starting point of further simulation processes of isolated and clustered defects in material lattices. This allows more materials, lattices, neutron incident energy ranges, and irradiation conditions to be explored with sufficient data to adequately cover both standard and novel applications and materials: the broader reactor applications landscape. The damage metrics of of materials are systematically explored under typical but different reactor's type environment. The inventory code FISPACT-II combined with the enhanced nuclear data forms of the TENDL-2019 libraries allow one to not only calculate dpa from mostly scattering events but to also properly predict gas production, nuclear heating and transmutation under the same conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Impulse Response Results from FOPEN Simulations

Foliage penetration (FOPEN) radar at lower frequencies (VHF, UHF) is a well-studied area with many contributions. However, there is growing interest in using higher Ku-band frequencies (12-18 GHz) for FOPEN. Specifically, the reduced wavelength sizes provide some key saliencies for developing more optimized detection solutions. The disadvantage is that exploiting Ku-band for FOPEN is complicated because higher frequencies have pronounced scattering effects due to their smaller wavelengths. A methodology has been developed to model and simulate FOPEN problems that characterize the phenomenology of Ku-band EM wave transmissions through moderate foliage. The details of this research are documented in multiple reports. The main focus of this report is to describe the FOPEN model simulation scene setup, validation and results.

42 ENGINEERING↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗