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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , and α Using PARTISN

SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff , and α using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every isotope and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations. The tool can be used for one-dimensional spherical (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used for fixed-source and eigenvalue problems. The tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists all of the input files whose results are presented in Sec. VI.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

S -matrix positivity without Lorentz invariance: a case study

We investigate the analytic structure of scattering amplitudes in theories in which Lorentz invariance is spontaneously broken. We do so by computing and studying the S-matrix for a simple example: a superfluid described by a complex scalar with quartic interactions. The computation is confined to tree-level, for there are no absolutely stable single-particle states, though the lifetime can be made long by lowering the chemical potential. For the 2 → 2 amplitude in center-of-mass configurations, not only is crossing symmetry violated, there appears a tree level branch cut for unphysical kinematics. Its appearance is a consequence of non-analyticity in the dispersion relation. The branch point defines a new scale in the problem, which scales inversely with the chemical potential. In this example, even derivatives of the forward amplitude are positive while odd derivatives are negative. This pattern can be understood in a general way in the limit of a small chemical potential, or weak Lorentz breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Minimum possible laser pulse duration for SRS

The minimum laser pulse duration for transient stimulated Raman scattering is shown to be limited as a consequence of nonlinear phase modulation of the laser and Stokes waves. It is determined by the ratio of the nonlinear refractive index to the spontaneous Raman scattering cross section. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Theoretical studies of Resonance Enhanced Stimulated Raman Scattering (RESRS) of frequency doubled Alexandrite laser wavelengths in cesium vapor

This work focused on understanding the effects of arbitrary transverse and longitudinal relaxation rates on the susceptibilities of coherently driven three-level systems. The approximation of a single relaxation rate often made in previous work is strongly invalidated by the variation in the spontaneous emission lifetime between various atomic level pairs in systems such as cesium. It is of great importance to the problem of nonlinear infrared generation to determine the dependence of both real and imaginary susceptibility on relaxation rates. The imaginary susceptibility on the pump transition determines the absorption of pump photons and the imaginary susceptibility on the laser transition determines the spectral dependence of the gain. This is of particular importance for pure Raman emission (i.e., absorption at linecenter of the gain transition) as it determines the tunability characteristics we are aiming to predict. The real susceptibility is important when cavities are used at the signal field as this will determine the loaded resonance of the Raman oscillator. Researchers show that in some cases which result from having different relaxation rates mode splitting may result, allowing more than one frequency to have the same Raman wavelength, possibly resulting in a temporal instability.

Lawandy, N. M.↗
DOE OSTI2020

3-3-1 Self interacting dark matter and the galaxy core-cusp problem

The core–cusp problem remains a challenging discrepancy between observations and simulations in the standard Λ CDM model for the formation of galaxies. The problem is that Λ CDM simulations predict a steep power-law mass density profile at the center of galactic dark matter halos. However, observations of dwarf galaxies in the Local Group reveal a density profile consistent with a nearly flat distribution of dark matter near the center. A number of solutions to this dilemma have been proposed. In this paper we summarize investigations into the possibility that the dark matter particles themselves self interact and scatter. Such self-interacting dark matter (SIDM) particles can smooth out the dark-matter profile in high-density regions. We also review the theoretical proposal that self-interacting dark matter may arise as an additional Higgs scalar in the 3–3–1 extension of the Standard Model (SM). We present new simulations of galaxy formation and evolution for this formulation of self-interacting dark matter. Current constraints on this self-interacting dark matter are then summarized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Smart Scattering Scanning Near-Field Optical Microscopy

Scattering scanning near-field optical microscopy (s-SNOM) provides spectroscopic imaging from molecular to quantum materials with few nanometer deep subdiffraction limited spatial resolution. However, in its conventional implementation s-SNOM is slow to effectively acquire a series of spatio-spectral images, especially with large fields of view. This problem is further exacerbated for weak resonance contrast or when using light sources with limited spectral irradiance. Indeed, the generally limited signal-to-noise ratio prevents sampling a weak signal at the Nyquist sampling rate. Here, we demonstrate how acquisition time and sampling rate can be significantly reduced by using compressed sampling, matrix completion, and adaptive random sampling, while maintaining or even enhancing the physical or chemical image content. We use fully sampled real data sets of molecular, biological, and quantum materials as ground-truth physical data and show how deep under-sampling with a corresponding reduction of acquisition time by 1 order of magnitude or more retains the core s-SNOM image information. We demonstrate that a sampling rate of up to 6× smaller than the Nyquist criterion can be applied, which would provide a 30-fold reduction in the data required under typical experimental conditions. Furthermore, our smart s-SNOM approach is generally applicable and provides systematic full spatio-spectral s-SNOM imaging with a large field of view at high spectral resolution and reduced acquisition time.

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

Spin texture induced by non-magnetic doping and spin dynamics in 2D triangular lattice antiferromagnet h -Y(Mn,Al)O 3

Novel effects induced by nonmagnetic impurities in frustrated magnets and quantum spin liquid represent a highly nontrivial and interesting problem. A theoretical proposal of extended modulated spin structures induced by doping of such magnets, distinct from the well-known skyrmions has attracted significant interest. Here, we demonstrate that nonmagnetic impurities can produce such extended spin structures in h-YMnO 3 , a triangular antiferromagnet with noncollinear magnetic order. Using inelastic neutron scattering (INS), we measured the full dynamical structure factor in Al-doped h-YMnO 3 and confirmed the presence of magnon damping with a clear momentum dependence. Our theoretical calculations can reproduce the key features of the INS data, supporting the formation of the proposed spin textures. As such, our study provides the first experimental confirmation of the impurity-induced spin textures. It offers new insights and understanding of the impurity effects in a broad class of noncollinear magnetic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1974

Influence of rotation on the maximum mass of pulsationally stable stars

The critical mass for stability against radial pulsations in rotating, homogeneous main-sequence stars is found to be greater than in the case of no rotation. Analytic and detailed numerical models show that the critical mass rises steeply with increasing concentration of angular momentum to the center of the star. For uniform rotation near breakup velocity at the star's equator the critical mass is about 850 solar masses if an electron-scattering opacity is used, or about 5000 solar masses if the opacities of Cox and Stewart are used. For nonuniform rotation with a constant ratio of centrifugal force to gravity in the star, the critical mass becomes 'infinite' long before breakup velocity is attained. The relevance of the present results to several observational problems is noted.

Stothers, R.↗
NASA NTRS1985

Electron impact excitation of atomic oxygen - Revised cross sections

Revised cross-section values for the excitation of three O I resonance transitions at 1304, 1027, and 989 A, by electron impact on atomic oxygen are presented from threshold to 300 eV. These results are smaller than the excitation cross sections used in some airglow models by a factor of about 2.8. The revised values are in good agreement with recent quantum-scattering calculations. The downward revision is required by new laboratory studies in which the direct and dissociative cross sections for 1304 A excitation were normalized with small probable error to the O and O2 ionization cross sections. The results also reflect new advances in VUV optical calibration techniques. A number of outstanding airglow problems are simplified by these revisions.

Zipf, E. C.↗
NASA NTRS1993

An efficient model for coupling structural vibrations with acoustic radiation

The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.

Frendi, Abdelkader↗
NASA NTRS2004

Overcoming Geometry-Induced Stiffness with IMplicit-Explicit (IMEX) Runge-Kutta Algorithms on Unstructured Grids with Applications to CEM, CFD, and CAA

My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.

Kanevsky, Alex↗
NASA NTRS2016

A Possible Explanation for the Z-R Parameter Inconsistency when Comparing Stratiform and Convective Rainfall

The well-known Z-R power law Z A R(sup b) uses two parameters, A and b, in order to relate rainfall rate R to measured weather radar reflectivity Z. A common method used by researchers is to compute Z and R from disdrometer data and then extract the A-b parameter pair from a log-linear line fit to a scatter plot of Z-R pairs. Even though it may seem far more truthful to extract the parameter pair from a fit of radar Z(sub R) versus gauge rainfall rate R(sub G), the extreme difference in spatial and temporal sampling volumes between radar and rain gauge creates a slew of problems that can generally only be solved by using rain gauge arrays and long sampling averages. Disdrometer derived A-b parameters are easily obtained and can provide information for the study of stratiform versus convective rainfall. However, an inconsistency appears when comparing averaged A-b pairs from various researchers. Values of b range from 1.26 to 1.51 for both stratiform and convective events. Paradoxically the values of A fall into three groups: 150 to 200 for convective; 200 to 400 for stratiform; and 400 to 500 again for convective. This apparent inconsistency can be explained by computing the A-b pair using the gamma DSD coupled with a modified drop terminal velocity model, v(D) alpha D(sup beta) - w, where w is a somewhat artificial constant vertical velocity of the air above the disdrometer. This model predicts three regions of A, corresponding to w less than 0, w = 0, and w greater than 0, which approximately matches observed data.

gamma DSD↗
NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗
NASA NTRS2020

Electromagnetic Scattering by Discrete Random Media Illuminated by a Gaussian Beam II: Solution of the Radiative Transfer Equation

In this paper, we present numerical methods for solving the phenomenological scalar radiative transfer equation for a discrete random medium illuminated by a Gaussian beam. These rely on the Fourier transform method for the horizontal variables and the discrete ordinate method with matrix exponential for solving the underlying one-dimensional radiative transfer equation in the wavenumber domain. The problem of a Gaussian beam at oblique and normal incidence, as well as, the searchlight problem are treated. A complete description of the methods and the numerical algorithms is provided.

Gaussian beam↗
DOE OSTI2019

A composite neural network that learns from multi-fidelity data: Application to function approximation and inverse PDE problems

Presently the training of neural networks relies on data of comparable accuracy but in real applications only a very small set of high-fidelity data is available while inexpensive lower fidelity data may be plentiful. We propose a new composite neural network (NN) that can be trained based on multi-fidelity data. It is comprised of three NNs, with the first NN trained using the low-fidelity data and coupled to two high-fidelity NNs, one with activation functions and another one without, in order to discover and exploit nonlinear and linear correlations, respectively, between the low-fidelity and the high-fidelity data. We first demonstrate the accuracy of the new multi-fidelity NN for approximating some standard benchmark functions but also a 20-dimensional function that is not easy to approximate with other methods, e.g. Gaussian process regression. Subsequently, we extend the recently developed physics-informed neural networks (PINNs) to be trained with multi-fidelity data sets (MPINNs). MPINNs contain four fully-connected neural networks, where the first one approximates the low-fidelity data, while the second and third construct the correlation between the low- and high-fidelity data and produce the multi-fidelity approximation, which is then used in the last NN that encodes the partial differential equations (PDEs). Specifically, by decomposing the correlation into a linear and nonlinear part, the present model is capable of learning both the linear and complex nonlinear correlations between the low- and high-fidelity data adaptively. By training the MPINNs, we can: (1) obtain the correlation between the low- and high-fidelity data, (2) infer the quantities of interest based on a few scattered data, and (3) identify the unknown parameters in the PDEs. In particular, we employ the MPINNs to learn the hydraulic conductivity field for unsaturated flows as well as the reactive models for reactive transport. The results demonstrate that MPINNs can achieve relatively high accuracy based on a very small set of high-fidelity data. Despite the relatively low dimension and limited number of fidelities (two-fidelity levels) for the benchmark problems in the present study, the proposed model can be readily extended to very high-dimensional regression and classification problems involving multi-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Efficient structure-factor modeling for crystals with multiple components

Diffraction intensities from a crystallographic experiment include contributions from the entire unit cell of the crystal: the macromolecule, the solvent around it and eventually other compounds. These contributions cannot typically be well described by an atomic model alone, i.e . using point scatterers. Indeed, entities such as disordered (bulk) solvent, semi-ordered solvent ( e.g. lipid belts in membrane proteins, ligands, ion channels) and disordered polymer loops require other types of modeling than a collection of individual atoms. This results in the model structure factors containing multiple contributions. Most macromolecular applications assume two-component structure factors: one component arising from the atomic model and the second one describing the bulk solvent. A more accurate and detailed modeling of the disordered regions of the crystal will naturally require more than two components in the structure factors, which presents algorithmic and computational challenges. Here an efficient solution of this problem is proposed. All algorithms described in this work have been implemented in the computational crystallography toolbox ( CCTBX ) and are also available within Phenix software. These algorithms are rather general and do not use any assumptions about molecule type or size nor about those of its components.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1975

Calculations of the cosmic ray modulation in interplanetary space taking into account the possible dependence of the transport travel for the scattering of the particles and of the velocity of the solar winds on the angles they make with the helioequator plane: The case of isotropic diffusion

The modulation of galactic cosmic rays is studied by the magnetic heterogeneities stream on the assumption that the diffusion coefficient is reduced whereas the solar wind velocity is increased with the growth of the angle between the sun's rotation axis and the direction of solar plasma motion. The stationary plane problem of isotropic diffusion is solved as it applies to two cases: (1) with due account of particle retardation by the antiphermium mechanism; and (2) without an account of the above mechanism. This problem is solved by the grid method in the polar coordinate system. The results of the calculations are followed by a discussion of the method of solution and of the errors.

Dorman, L. I.↗
DOE OSTI2022

A machine learning photon detection algorithm for coherent x-ray ultrafast fluctuation analysis

X-ray free electron laser experiments have brought unique capabilities and opened new directions in research, such as creating new states of matter or directly measuring atomic motion. One such area is the ability to use finely spaced sets of coherent x-ray pulses to be compared after scattering from a dynamic system at different times. This enables the study of fluctuations in many-body quantum systems at the level of the ultrafast pulse durations, but this method has been limited to a select number of examples and required complex and advanced analytical tools. By applying a new methodology to this problem, we have made qualitative advances in three separate areas that will likely also find application to new fields. As compared to the “droplet-type” models, which typically are used to estimate the photon distributions on pixelated detectors to obtain the coherent x-ray speckle patterns, our algorithm achieves an order of magnitude speedup on CPU hardware and two orders of magnitude improvement on GPU hardware. We also find that it retains accuracy in low-contrast conditions, which is the typical regime for many experiments in structural dynamics. Finally, it can predict photon distributions in high average-intensity applications, a regime which up until now has not been accessible. Our artificial intelligence-assisted algorithm will enable a wider adoption of x-ray coherence spectroscopies, by both automating previously challenging analyses and enabling new experiments that were not otherwise feasible without the developments described in this work.

47 OTHER INSTRUMENTATION↗