Search NASA⌕ Search

SEARCH · Search NASA

Results for “T-matrix method”

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.

63 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Ab initio method for calculating total cross sections

A method for calculating total cross sections without formally including nonelastic channels is presented. The idea is to use a one channel T-matrix variational principle with a complex correlation function. The derived T matrix is therefore not unitary. Elastic scattering is calculated from T-parallel-squared, but total scattering is derived from the imaginary part of T using the optical theorem. The method is applied to the spherically symmetric model of electron-hydrogen scattering. No spurious structure arises; results for sigma(el) and sigma(total) are in excellent agreement with calculations of Callaway and Oza (1984). The method has wide potential applicability.

Bhatia, A. K.↗
NASA NTRS2010

Extinction-to-Backscatter Ratios of Saharan Dust Layers Derived from In-Situ Measurements and CALIPSO Overflights During NAMMA

We determine the extinction-to-backscatter (Sa) ratios of dust using (1) airborne in-situ measurements of microphysical properties, (2) modeling studies, and (3) the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) observations recorded during the NASA African Monsoon Multidisciplinary Analyses (NAMMA) field experiment conducted from Sal, Cape Verde during Aug-Sept 2006. Using CALIPSO measurements of the attenuated backscatter of lofted Saharan dust layers, we apply the transmittance technique to estimate dust Sa ratios at 532 nm and a 2-color method to determine the corresponding 1064 nm Sa. This method yielded dust Sa ratios of 39.8 plus or minus 1.4 sr and 51.8 plus or minus 3.6 sr at 532 nm and 1064 nm, respectively. Secondly, Sa at both wavelengths is independently calculated using size distributions measured aboard the NASA DC-8 and estimates of Saharan dust complex refractive indices applied in a T-Matrix scheme. We found Sa ratios of 39.1 plus or minus 3.5 sr and 50.0 plus or minus 4 sr at 532 nm and 1064 nm, respectively, using the T-Matrix calculations applied to measured size spectra. Finally, in situ measurements of the total scattering (550 nm) and absorption coefficients (532 nm) are used to generate an extinction profile that is used to constrain the CALIPSO 532 nm extinction profile and thus generate a stratified 532 nm Sa. This method yielded an Sa ratio at 532 nm of 35.7 sr in the dust layer and 25 sr in the marine boundary layer consistent with a predominantly seasalt aerosol near the ocean surface. Combinatorial simulations using noisy size spectra and refractive indices were used to estimate the mean and uncertainty (one standard deviation) of these Sa ratios. These simulations produced a mean (plus or minus uncertainty) of 39.4 (plus or minus 5.9) sr and 56.5 (plus or minus 16.5) sr at 532 nm and 1064 nm, respectively, corresponding to percent uncertainties of 15% and 29%. These results will provide a measurements-based estimate of the dust Sa for use in backscatter lidar inversion algorithms such as CALIOP.

Omar, Ali H.↗
NASA NTRS1997

Particle Scattering in the Resonance Regime: Full-Wave Solution for Axisymmetric Particles with Large Aspect Ratios

Reliable descriptions of the optical properties of clouds and aerosols are essential for studies of radiative transfer in planetary atmospheres. The scattering algorithms provide accurate estimates of these properties for spherical particles with a wide range of sizes and refractive indices, but these methods are not valid for non-spherical particles (e.g., ice crystals, mineral dust, and smoke). Even though a host of methods exist for deriving the optical properties of nonspherical particles that are very small or very large compared with the wavelength, only a few methods are valid in the resonance regime, where the particle dimensions are comparable with the wavelength. Most such methods are not ideal for particles with sharp edges or large axial ratios. We explore the utility of an integral equation approach for deriving the single-scattering optical properties of axisymmetric particles with large axial ratios. The accuracy of this technique is shown for spheres of increasing size parameters and an ensemble of randomly oriented prolate spheroids of size parameter equal to 10.079368. In this last case our results are compared with published results obtained with the T-matrix approach. Next we derive cross sections, single-scattering albedos, and phase functions for cylinders, disks, and spheroids of ice with dimensions extending from the Rayleigh to the geometric optics regime. Compared with those for a standard surface integral equation method, the storage requirement and the computer time needed by this method are reduced, thus making it attractive for generating databases to be used in multiple-scattering calculations. Our results show that water ice disks and cylinders are more strongly absorbing than equivalent volume spheres at most infrared wavelengths. The geometry of these particles also affects the angular dependence of the scattering. Disks and columns with maximum linear dimensions larger than the wavelength scatter much more radiation in the forward and backward directions and much less radiation at intermediate phase angles than equivalent volume spheres.

Zuffada, Cinzia↗
NASA NTRS1996

T-Matrix Computations of Light Scattering by Nonspherical Particles: A Review

We review the current status of Waterman's T-matrix approach which is one of the most powerful and widely used tools for accurately computing light scattering by nonspherical particles, both single and composite, based on directly solving Maxwell's equations. Specifically, we discuss the analytical method for computing orientationally-averaged light-scattering characteristics for ensembles of nonspherical particles, the methods for overcoming the numerical instability in calculating the T matrix for single nonspherical particles with large size parameters and/or extreme geometries, and the superposition approach for computing light scattering by composite/aggregated particles. Our discussion is accompanies by multiple numerical examples demonstrating the capabilities of the T-matrix approach and showing effects of nonsphericity of simple convex particles (spheroids) on light scattering.

Mischenko, Michael I.↗
NASA NTRS1989

Kaon-nucleus scattering

The derivations of the Lippmann-Schwinger equation and Watson multiple scattering are given. A simple optical potential is found to be the first term of that series. The number density distribution models of the nucleus, harmonic well, and Woods-Saxon are used without t-matrix taken from the scattering experiments. The parameterized two-body inputs, which are kaon-nucleon total cross sections, elastic slope parameters, and the ratio of the real to the imaginary part of the forward elastic scattering amplitude, are presented. The eikonal approximation was chosen as our solution method to estimate the total and absorptive cross sections for the kaon-nucleus scattering.

Hong, Byungsik↗
NASA NTRS1989

Kaon-nucleus scattering

Two kinds of number density distributions of the nucleus, harmonic well and Woods-Saxon models, are used with the t-matrix that is taken from the scattering experiments to find a simple optical potential. The parameterized two body inputs, which are kaon-nucleon total cross sections, elastic slope parameters, and the ratio of the real to imaginary part of the forward elastic scattering amplitude, are shown. The eikonal approximation was chosen as the solution method to estimate the total and absorptive cross sections for the kaon-nucleus scattering.

Hong, Byungsik↗
NASA NTRS2010

Extinction-to-Backscatter Ratios of Saharan Dust Layers Derived from In-Situ Measurements and CALIPSO Overflights During NAMMA

We determine the aerosol extinction-to-backscatter (Sa) ratios of dust using airborne in-situ measurements of microphysical properties, and CALIPSO observations during the NASA African Monsoon Multidisciplinary Analyses (NAMMA). The NAMMA field experiment was conducted from Sal, Cape Verde during Aug-Sept 2006. Using CALIPSO measurements of the attenuated backscatter of lofted Saharan dust layers, we apply the transmittance technique to estimate dust Sa ratios at 532 nm and a 2-color method to determine the corresponding 1064 nm Sa. Using this method, we found dust Sa ratios of 39.8 plus or minus 1.4 sr and 51.8 plus or minus 3.6 sr at 532 nm and 1064 nm, respectively. Secondly, Sa ratios at both wavelengths is independently calculated using size distributions measured aboard the NASA DC-8 and estimates of Saharan dust complex refractive indices applied in a T-Matrix scheme. We found Sa ratios of 39.1 plus or minus 3.5 sr and 50.0 plus or minus 4 sr at 532 nm and 1064 nm, respectively, using the T-Matrix calculations applied to measured size spectra. Finally, in situ measurements of the total scattering (550 nm) and absorption coefficients (532 nm) are used to generate an extinction profile that is used to constrain the CALIPSO 532 nm extinction profile.

Omar, Ali H.↗
NASA NTRSDate not supplied

Retrieving DSD Moments from GPM-DPR: A Simulation Study Based on the Full DSD Spectra Characterized by the Generalized Gamma Model

Recently, a method to retrieve rain drop size distribution (DSD) moments from X-band dual polarization radar measurements has been developed using copolar reflectivity, differential reflectivity and specific attenuation. Two reference moments are retrieved first, followed by reconstructing the DSDs and calculating other moments using a generalized gamma model to represent the underlying shape, h(x), corresponding to the pair of chosen reference moments. Here we explore a similar approach but the retrieval method in this study uses dual-frequency radar measurements. In the case of GPM-DPR, the two frequencies are 13.8and 35 GHz, and their products include attenuation corrected reflecivities at both frequencies (Z and Z) and the specific attenuation (k, k, if available). Our approach is to use these products to determine two reference moments, namely M3 andM6 representing the third and the sixth moments respectively. Then, as with the polarimetric radar retrievals, we use the most probable h(x) to reconstruct the full DSD spectra, from which other moments are calculated. The best two DPR products for estimating M3 and M6 appear to be A and Z. Figure 1 shows the retrievals versus the ‘true’ moments. 2930 three-minute DSDs were used for the (T-matrix) scattering calculations at Ku and Ka bands but only cases with A > 0.5 dB/km were chosen. For the retrieved moments, the integration was performed only up to 6 mm drop diameter. The [1:1] line is included in Fig. 1. Even with just two DPR products as input the results seem promising, although the lower order moments show somewhat more scatter, especially the zeroth moment, M0. We will quantify the retrieval errors, and additionally examine the stability of h(x). For the latter, data from (i) Greeley, Colorado, (ii) Huntsville, Alabama, and (iii)Wallops, Virginia, will be used. Finally, a GPM overpass case over Huntsville, during a widespread rain event on 11 April 2016, will be considered as an initial test case.

precipitation↗
DOE OSTI2026

Active learning emulators for nuclear two-body scattering in momentum space

In this work we extend the active learning emulators for two-body scattering in coordinate space with error estimation, recently developed by Maldonado et al. [Phys. Rev. C 112, 024002], to coupled-channel scattering in momentum space. Our full-order model (FOM) solver is based on the Lippmann-Schwinger integral equation for the scattering t-matrix as opposed to the radial Schrödinger equation. We use (Petrov-)Galerkin projections and high-fidelity calculations at a few snapshots across the parameter space of the interaction to construct efficient reduced-order models (ROMs), trained by a greedy algorithm for locally optimal snapshot selection. Both the FOM solver and the corresponding ROMs are implemented efficiently in Python using Google's JAX library. We present results for emulating scattering phase shifts in coupled and uncoupled channels and cross sections, and assess the accuracy of the developed ROMs and their computational speedup factors. We also develop emulator error estimation for both the t-matrix and the total cross section. The software framework for reproducing and extending our results is publicly available. Together with our recent advances in developing active-learning emulators for three-body scattering, these emulator frameworks set the stage for full Bayesian calibrations of chiral nuclear interactions and optical models against scattering data with quantified emulator errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗