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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Sarabandi, Kamal↗
NASA NTRS2002

A Two-Dimensional Linear Bicharacteristic FDTD Method

The linear bicharacteristic scheme (LBS) was originally developed to improve unsteady solutions in computational acoustics and aeroacoustics. The LBS has previously been extended to treat lossy materials for one-dimensional problems. It is a classical leapfrog algorithm, but is combined with upwind bias in the spatial derivatives. This approach preserves the time-reversibility of the leapfrog algorithm, which results in no dissipation, and it permits more flexibility by the ability to adopt a characteristic based method. The use of characteristic variables allows the LBS to include the Perfectly Matched Layer boundary condition with no added storage or complexity. The LBS offers a central storage approach with lower dispersion than the Yee algorithm, plus it generalizes much easier to nonuniform grids. It has previously been applied to two and three-dimensional free-space electromagnetic propagation and scattering problems. This paper extends the LBS to the two-dimensional case. Results are presented for point source radiation problems, and the FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗
NASA NTRS1992

A finite element: Boundary integral method for electromagnetic scattering

A method that combines the finite element and boundary integral techniques for the numerical solution of electromagnetic scattering problems is presented. The finite element method is well known for requiring a low order storage and for its capability to model inhomogeneous structures. Of particular emphasis in this work is the reduction of the storage requirement by terminating the finite element mesh on a boundary in a fashion which renders the boundary integrals in convolutional form. The fast Fourier transform is then used to evaluate these integrals in a conjugate gradient solver, without a need to generate the actual matrix. This method has a marked advantage over traditional integral equation approaches with respect to the storage requirement of highly inhomogeneous structures. Rectangular, circular, and ogival mesh termination boundaries are examined for two-dimensional scattering. In the case of axially symmetric structures, the boundary integral matrix storage is reduced by exploiting matrix symmetries and solving the resulting system via the conjugate gradient method. In each case several results are presented for various scatterers aimed at validating the method and providing an assessment of its capabilities. Important in methods incorporating boundary integral equations is the issue of internal resonance. A method is implemented for their removal, and is shown to be effective in the two-dimensional and three-dimensional applications.

Collins, J. D.↗
NASA NTRS1987

New main reflector, subreflector and dual chamber concepts for compact range applications

A compact range is a facility used for the measurement of antenna radiation and target scattering problems. Most presently available parabolic reflectors do not produce ideal uniform plane waves in the target zone. Design improvements are suggested to reduce the amplitude taper, ripple and cross polarization errors. The ripple caused by diffractions from the reflector edges can be reduced by adding blended rolled edges and shaping the edge contour. Since the reflected edge continues smoothly from the parabola onto the rolled surface, rather than being abruptly terminated, the discontinuity in the reflected field is reduced which results in weaker diffracted fields. This is done by blending the rolled edges from the parabola into an ellipse. An algorithm which enables one to design optimum blended rolled edges was developed that is based on an analysis of the continuity of the surface radius of curvature and its derivatives across the junction. Futhermore, a concave edge contour results in a divergent diffracted ray pattern and hence less stray energy in the target zone. Design equations for three-dimensional reflectors are given. Various examples were analyzed using a new physical optics method which eliminates the effects of the false scattering centers on the incident shadow boundaries. A Gregorian subreflector system, in which both the subreflector and feed axes are tilted, results in a substantial reduction in the amplitude taper and cross polarization errors. A dual chamber configuration is proposed to eliminate the effects of diffraction from the subreflector and spillover from the feed. A computationally efficient technique, based on ray tracing and aperture integration, was developed to analyze the scattering from a lossy dielectric slab with a wedge termination.

Pistorius, C. W. I.↗
NASA NTRS1995

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗
DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗
DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Machine-learning-assisted insight into spin ice Dy 2 Ti 2 O 7

Complex behavior poses challenges in extracting models from experiment. An example is spin liquid formation in frustrated magnets like Dy 2 Ti 2 O 7 . Understanding has been hindered by issues including disorder, glass formation, and interpretation of scattering data. Here, we use an automated capability to extract model Hamiltonians from data, and to identify different magnetic regimes. This involves training an autoencoder to learn a compressed representation of three-dimensional diffuse scattering, over a wide range of spin Hamiltonians. The autoencoder finds optimal matches according to scattering and heat capacity data and provides confidence intervals. Validation tests indicate that our optimal Hamiltonian accurately predicts temperature and field dependence of both magnetic structure and magnetization, as well as glass formation and irreversibility in Dy 2 Ti 2 O 7 . The autoencoder can also categorize different magnetic behaviors and eliminate background noise and artifacts in raw data. Our methodology is readily applicable to other materials and types of scattering problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Conservation law for electromagnetic fields in a space-time-varying medium and its implications

We discuss the propagation of electromagnetic waves in a space-time modulated medium, for which neither energy nor momentum remains conserved. We have found a new conservation law when modulation moves at a constant speed v in a traveling-wave fashion. We use this law to study how light is scattered inside such a medium. We consider both the subluminal (v < c) andsuperluminal (v > c) cases and find that they differ considerably in terms of what remains conserved. We show that the total number of photons is conserved in the subluminal case. In contrast, it is the difference between the forward and backward moving photons that remains conserved in the superluminal case. We also study the reciprocity issue for a space-time modulated medium and found that reciprocity holds in all situations. Here, we develop a quantum formulation of our scattering problem and show that the subluminal case is similar to the action of a beam splitter. In the superluminal case, photon pairs can be generated from vacuum in a fashion analogous to the two-mode frequency down-conversion process.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hamiltonian parameter inference from resonant inelastic x-ray scattering with active learning

Identifying model Hamiltonians is a vital step toward creating predictive models of materials. Here, in this study, we combine Bayesian optimization with the EDRIXS numerical package to infer Hamiltonian parameters from resonant inelastic x-ray scattering (RIXS) spectra within the single atom approximation. To evaluate the efficacy of our method, we test it on experimental RIXS spectra of NiPS 3 , NiCl 2 , Ca 3 ⁢LiOsO 6 , and Fe 2⁢ O 3 , and demonstrate that it can reproduce results obtained from hand-fitted parameters to a precision similar to expert human analysis while providing a more systematic mapping of parameter space. Our work provides a key first step toward solving the inverse scattering problem to extract effective multi-orbital models from information-dense RIXS measurements, which can be applied to a host of quantum materials. We also propose atomic model parameter sets for two materials, Ca 3⁢ LiOsO 6 and Fe 2⁢ O 3 , that were previously missing from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis

In this work, a butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N log 1.5 N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

42 ENGINEERING↗
DOE OSTI2021

Electron ptychography achieves atomic-resolution limits set by lattice vibrations

Transmission electron microscopes use electrons with wavelengths of a few picometers, potentially capable of imaging individual atoms in solids at a resolution ultimately set by the intrinsic size of an atom. However, owing to lens aberrations and multiple scattering of electrons in the sample, the image resolution is reduced by a factor of 3 to 10. By inversely solving the multiple scattering problem and overcoming the electron-probe aberrations using electron ptychography, we demonstrate an instrumental blurring of less than 20 picometers and a linear phase response in thick samples. The measured widths of atomic columns are limited by thermal fluctuations of the atoms. As a result, our method is also capable of locating embedded atomic dopant atoms in all three dimensions with subnanometer precision from only a single projection measurement.

42 ENGINEERING↗
DOE OSTI2020

Lattice phase shifts and mixing angles for an arbitrary number of coupled channels

We present a lattice method for determining scattering phase shifts and mixing angles for the case of an arbitrary number of coupled channels. The proposed method combines a spherical wall boundary condition and a channel-mixing auxiliary potential to extract the full-rank S-matrix from the radial wave functions. We consider the scattering problem of two spin-1 bosons interacting with a test potential involving up to four coupled channels. For this benchmark system, the phase shifts and mixing angles are shown to agree on the lattice and in the continuum. Our method should allow to extend previous two-channel nuclear lattice EFT simulations to mixing of more than two partial waves.

Bovermann, Lukas↗
DOE OSTI2021

Toward nonthermal control of excited quantum materials: framework and investigations by ultrafast electron scattering and imaging

Quantum material systems upon applying ultrashort laser pulses provide a rich platform to access excited material phases and their transformations that are not entirely like their equilibrium counterparts. The addressability and potential controls of metastable or long-trapped out-of-equilibrium phases have motivated interests both for the purposes of understanding the nonequilibrium physics and advancing the quantum technologies. Thus far, the dynamical spectroscopic probes eminently focus on microscopic electronic and phonon responses. For characterizing the long-range dynamics, such as order parameter fields and fluctuation effects, the ultrafast scattering probes offer direct sensitivity. Bridging the connections between the microscopic dynamics and macroscopic responses is central toward establishing the nonequilibrium physics behind the light-induced phases. Here, we present a path toward such understanding by cross-examining the structure factors associated with different dynamical states obtained from ultrafast electron scattering, imaging, and modeling. We give the basic theoretical framework on describing the non-equilibrium scattering problems and briefly describe how such framework relates to the out-of-equilibrium phenomena. We give effective models outlining the emergences of nonthermal critical points, hidden phases, and non-equilibrium relaxational responses from vacuum-suspended rare-earth tritellurides, tantalum disulfides thin films, and vanadium dioxide nanocrystalline materials upon light excitations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Generalized pseudopotential approach for electron-atom scattering.

A generalized many-electron pseudopotential approach is presented for electron-neutral-atom scattering problems. A calculation based on this formulation is carried out for the singlet s-wave and p-wave electron-hydrogen phase shifts with excellent results. We compare the method with other approaches as well as discuss its applications for inelastic and rearrangement collision problems.

Zarlingo, D. G.↗