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 829 records · Page 46

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↗

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↗

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗

Computing Resonant Inelastic X-Ray Scattering Spectra Using The Density Matrix Renormalization Group Method

We present a method for computing the resonant inelastic x-ray scattering (RIXS) spectra in one-dimensional systems using the density matrix renormalization group (DMRG) method. By using DMRG to address this problem, we shift the computational bottleneck from the memory requirements associated with exact diagonalization (ED) calculations to the computational time associated with the DMRG algorithm. This approach is then used to obtain RIXS spectra on cluster sizes well beyond state-of-the-art ED techniques. Using this new procedure, we compute the low-energy magnetic excitations observed in Cu L-edge RIXS for the challenging corner shared CuO4 chains, both for large multi-orbital clusters and downfolded t-J chains. We are able to directly compare results obtained from both models defined in clusters with identical momentum resolution. In the strong coupling limit, we find that the downfolded t-J model captures the main features of the magnetic excitations probed by RIXS only after a uniform scaling of the spectra is made.

97 MATHEMATICS AND COMPUTING↗

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↗

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)↗

Unsteady temperature distribution in volume reflectors

The development of unsteady temperature distributions in semitransparent scattering materials under the influence of arbitrary incident radiative flux with other than specified temperature boundary conditions is considered. Specifically, such a problem is considered for an approximate radiative transfer model which allows analytic solutions. The model represents the unsteady, one-dimensional, radiative heating of a weakly absorbing, semitransparent medium with specified conductive heat flux at one boundary. The radiation field is modeled using an approximation to the radiative heat flux obtained by solution of the Kubelka-Munk differential equations.

Weston, K. C.↗

Resonance-line transfer with partial redistribution. VII Angle-dependent redistribution

A method is presented for treating radiative transfer in resonance lines, allowing for the full angle and frequency dependence of redistribution in the scattering process, as seen in the laboratory frame. The case of an equivalent-two-level-atom source function is considered; the problem to be treated is then linear in the radiation field. We apply this method to the Ca II lines in the solar atmosphere, using a redistribution function which takes into account a mixture of coherence in the atom's frame, with Doppler redistribution in the laboratory frame (for atoms which have not suffered an elastic collision), and of complete redistribution in the laboratory frame (for atoms that are collisionally perturbed during the emission process). Both the angle-averaged approximation and the full angle-dependent solution were obtained, and were compared to assess, differentially, the effects of angular redistribution upon the computed line profile and its center-to-limb behavior. For the Ca II line in a homogeneous solar chromosphere the angle-dependent effects are found to be negligible, indicating that one may use angle-averaged redistribution functions when studying partial redistribution effects in line profiles.

Milkey, R. W.↗

Thermal microwave emission from a random inhomogeneous layer over a homogeneous medium using the method of invariant imbedding

The paper studies thermal microwave emission from an inhomogeneous slab of a random medium, with possible nonuniform absorption, scattering, and temperature profiles, bounded by different dielectrics on both sides. The invariant imbedding method is used to cast the boundary value problem of the radiative transfer equations into an initial value problem at zero slab thickness. As a numerical example, the angular and polarization variations of brightness temperatures for ice over water are considered.

Tsang, L.↗

Super-alfvenic propagation of cosmic rays: The role of streaming modes

Numerous cosmic ray propagation and acceleration problems require knowledge of the propagation speed of relativistic particles through an ambient plasma. Previous calculations indicated that self-generated turbulence scatters relativistic particles and reduces their bulk streaming velocity to the Alfven speed. This result was incorporated into all currently prominent theories of cosmic ray acceleration and propagation. It is demonstrated that super-Alfvenic propagation is indeed possible for a wide range of physical parameters. This fact dramatically affects the predictions of these models.

Morrison, P. J.↗

Speed-dependent collision effects on radar back-scattering from the ionosphere

A computer code to accurately compute the fluctuation spectrum for linearly speed dependent collision frequencies was developed. The effect of ignoring the speed dependence on the estimates of ionospheric parameters was determined. It is shown that disagreements between the rocket and the incoherent scatter estimates could be partially resolved if the correct speed dependence of the i-n collision frequency is not ignored. This problem is also relevant to the study of ionospheric irregularities in the auroral E-region and their effects on the radio communication with satellites.

Theimer, O.↗

Modeling of vegetation canopy reflectance: Status, issues and recommended future strategy

Various technical issues related to mapping of vegetative type, condition and stage of maturity, utilizing remotely sensed spectral data are reviewed. The existing knowledge base of models, especially of radiative properties of the vegetation canopy and atmosphere, is reviewed to establish the state of the art for addressing the problem of vegetation mapping. Activities to advance the state of the art are recommended. They include working on canopy reflectance and atmospheric scattering models, and field measurements of canopy reflectance as well as of canopy components. Leaf area index (LAI) and solar radiation interception (SRI) are identified as the two most important vegetation variables requiring further investigation. It is recommended that activities related to sensing them or understanding their relationships with measurable variables, should be encouraged and supported.

Goel, N. S.↗

RIMSTAR: Rapid image scanning technology and archiving

Conventional schemes for digitizing large volumes of photographic data are far too costly and time consuming to encourage the undertaking of the ambitious projects using traditional technology. There is a distinct need for fast digitizing systems such as could be derived from development of large format, optically multiplexed CCD systems designed to address these problems. The use of CCD chips for data handling and the capability of using polaroids in the optical path for reduced light scattering are discussed.

Craine, E. R.↗

On the electromagnetic scattering from infinite rectangular conducting grids

The study and development of two numerical techniques for the analysis of electromagnetic scattering from a rectangular wire mesh are described. Both techniques follow from one basic formulation and they are both solved in the spectral domain. These techniques were developed as a result of an investigation towards more efficient numerical computation for mesh scattering. These techniques are efficient for the following reasons: (a1) make use of the Fast Fourier Transform; (b2) they avoid any convolution problems by converting integrodifferential equations into algebraic equations; and (c3) they do not require inversions of any matrices. The first method, the SIT or Spectral Iteration Technique, is applied for regions where the spacing between wires is not less than two wavelengths. The second method, the SDCG or Spectral Domain Conjugate Gradient approach, can be used for any spacing between adjacent wires. A study of electromagnetic wave properties, such as reflection coefficient, induced currents and aperture fields, as functions of frequency, angle of incidence, polarization and thickness of wires is presented. Examples and comparisons or results with other methods are also included to support the validity of the new algorithms.

Christodoulou, C.↗

The Tukey algorithm for enhancing MST radar data

One of the most troublesome features in MSR (mesosphere stratosphere troposphere) velocity measurements is the determination of unwanted scatterer whose velocity is different from that of the surrounding atmosphere. Aircraft seen in the sidelobes of the antenna are the principal problem. Because coherent integration essentially eliminates echoes with line of sight velocities greater than 10 or 20 m/s, aircraft are seen only when their flight path is almost perpendicular to the line of sight. Then, they give large returns whose velocities may be positive or negative, and certainly different from that of the surrounding air. The glitches in the minute by minute velocity records are quite troublesome in that they may distort the statistics of the velocity. An objective way is therefore needed to remove sporadic points of this kind. For this purpose, the Tukey algorithm is appropriate and has some advantages over averaging. The Tukey algorithm, applied to a data array, uses for each data point the median of it and the two points surrounding it. If the three points form a monotonically increasing or decreasing sequence, the original point is copied with change. However, if the central point is remote from the other two, it is replaced by whichever of the two surrounding points is closest in value. The greatest effect of the Tukey algorithm is on data where the successive points are uncorrelated. Examples are given.

Bowhill, S. A.↗

Particle acceleration at shocks - A Monte Carlo method

A Monte Carlo method is presented for the problem of acceleration of test particles at relativistic shocks. The particles are assumed to diffuse in pitch angle as a result of scattering off magnetic irregularities frozen into the fluid. Several tests are performed using the analytic results available for both relativistic and nonrelativistic shock speeds. The acceleration at relativistic shocks under the influence of radiation losses is investigated, including the effects of a momentum dependence in the diffusion coefficient. The results demonstrate the usefulness of the technique in those situations in which the diffusion approximation cannot be employed, such as when relativistic bulk motion is considered, when particles are permitted to escape at the boundaries, and when the effects of the finite length of the particle mean free path are important.

Kirk, J. G.↗

D-wave Resonances in Positronium Hydride

We have recently reexamined the problem of calculating the positions and widths of the lowest-lying S- and P-state resonances in the Ps+H scattering system. The resonances are generated by bound Rydberg states in the closed rearranged channel (e(+)+H(-)), shifted somewhat by coupling with the open scattering channel. This picture had been successfully applied for S-waves, but more recent results showed unexpectedly large shifts in both P- and D-wave resonances. Our recalculation for the S- and P-waves showed that the series of resonances begins with the second Rydberg state, giving rise to an apparent large shift. Now we have extended the calculations to D states. Here we find the expected results: The 3D Rydberg state produces a resonance very close to the result of Yan and Ho, and inclusion of the 4D state gives a second resonance only slightly shifted.

Drachman, Richard J.↗

Analysis of blended rolled edge reflectors using numerical UTD

The uniform geometrical theory of diffraction (UTD) concept is used to predict the scattered fields in the target zones of compact-range blended rolled-edge reflectors. Since the necessary diffraction coefficients are not known in a closed form, a numerical method to calculate the diffraction coefficients is described. In the numerical method, the problem is reduced to two dimensions, and physical optics line integration is used to compute the diffraction coefficients. Thus, the method is computationally efficient. The method is used to analyze two compact-range reflectors. The results obtained using the numerical UTD show good agreement with the scattered fields obtained using a corrected physical optics surface integration.

Ellingson, Steven W.↗