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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Electromagnetic Scattering by Spheroidal Volumes of Discrete Random Medium

We use the superposition T-matrix method to compare the far-field scattering matrices generated by spheroidal and spherical volumes of discrete random medium having the same volume and populated by identical spherical particles. Our results fully confirm the robustness of the previously identified coherent and diffuse scattering regimes and associated optical phenomena exhibited by spherical particulate volumes and support their explanation in terms of the interference phenomenon coupled with the order-of-scattering expansion of the far-field Foldy equations. We also show that increasing non-sphericity of particulate volumes causes discernible (albeit less pronounced) optical effects in forward and backscattering directions and explain them in terms of the same interference/multiple-scattering phenomenon.

Electromagnetic scattering; Multi-particle groups;↗
NASA NTRS2017

Retrieval of Microphysical Characteristics of Particles in Atmospheres of Distant Comets from Ground-Based Polarimetry

We summarize unique aperture data on the degree of linear polarization observed for distant comets C/2010 S1, C/2010 R1, C/2011 KP36, C/2012 J1, C/2013 V4, and C/2014 A4 with heliocentric distances exceeding 3 AU. Observations have been carried out at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences (Nizhnij Arkhyz, Russia) during the period from 2011 to 2016. The measured negative polarization proves to be significantly larger in absolute value than what is typically observed for comets close to the Sun. We compare the new observational data with the results of numerical modeling performed with the T-matrix and superposition T-matrix methods. In our computer simulations, we assume the cometary coma to be an optically thin cloud containing particles in the form of spheroids, fractal aggregates composed of spherical monomers, and mixtures of spheroids and aggregate particles. We obtain a good semi-quantitative agreement between all polarimetric data for the observed distant comets and the results of numerical modeling for the following models of the cometary dust: (i) a mixture of submicrometer water-ice oblate spheroids with aggregates composed of submicrometer silicate monomers; and (ii) a mixture of submicrometer water-ice oblate spheroids and aggregates consisting of both silicate and organic monomers. The microphysical parameters of these models are presented and discussed.

Numerical modeling↗
NASA NTRS2018

Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study

This paper provides a thorough modeling-based overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols. Using the numerically-exact superposition T-matrix method, we examine the absorption enhancement, absorption Angstroem exponent (AAE), backscattering linear depolarization ratio (LDR), and scattering matrix elements of black-carbon aerosols with 11 different model morphologies ranging from bare soot to completely embedded soot-sulfate and soot-brown carbon mixtures. Our size-averaged results show that fluffy soot particles absorb more light than compact bare-soot clusters. For the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Absorption increases as soot accumulates more coating material and can become saturated. The absorption enhancement is affected by particle size, morphology, wavelength, and the amount of coating. We refute the conventional belief that all carbonaceous aerosols have AAEs close to 1.0. Although LDRs caused by bare soot and certain carbonaceous particles are rather weak, LDRs generated by other soot-containing aerosols can reproduce strong depolarization measured by Burton et al. for aged smoke. We demonstrate that multi-wavelength LDR measurements can be used to identify the presence of morphologically complex carbonaceous particles, although additional observations can be needed for full characterization. Our results show that optical constants of the host/coating material can significantly influence the scattering and absorption properties of soot-containing aerosols to the extent of changing the sign of linear polarization. We conclude that for an accurate estimate of black-carbon radiative forcing, one must take into account the complex morphologies of carbonaceous aerosols in remote sensing studies as well as in atmospheric radiation computations.

carbonaceous aerosols; scattering matrix; polariza↗
NASA NTRS2001

Complex Correlation Kohn-T Method of Calculating Total and Elastic Cross Sections: Electron-Hydrogen Elastic Scattering - Part 1

We report on the first part of a study of electron-hydrogen scattering, using a method which allows for the ab initio calculation of total and elastic cross sections at higher energies. In its general form the method uses complex 'radial' correlation functions, in a (Kohn) T-matrix formalism. The titled method, abbreviated Complex Correlation Kohn T (CCKT) method, is reviewed, in the context of electron-hydrogen scattering, including the derivation of the equation for the (complex) scattering function, and the extraction of the scattering information from the latter. The calculation reported here is restricted to S-waves in the elastic region, where the correlation functions can be taken, without loss of generality, to be real. Phase shifts are calculated using Hylleraas-type correlation functions with up to 95 terms. Results are rigorous lower bounds; they are in general agreement with those of Schwartz, but they are more accurate and outside his error bounds at a couple of energies,

Bhatia, A. K.↗
NASA NTRS2014

Comprehensive Thematic T-matrix Reference Database: a 2013-2014 Update

This paper is the sixth update to the comprehensive thematic database of peer-reviewedT-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2013. It also lists several earlier publications not incorporated in the original database and previous updates.

T-matrix method↗
NASA NTRS2016

Applicability of the Effective-Medium Approximation to Heterogeneous Aerosol Particles.

The effective-medium approximation (EMA) is based on the assumption that a heterogeneous particle can have a homogeneous counterpart possessing similar scattering and absorption properties. We analyze the numerical accuracy of the EMA by comparing superposition T-matrix computations for spherical aerosol particles filled with numerous randomly distributed small inclusions and Lorenz-Mie computations based on the Maxwell-Garnett mixing rule. We verify numerically that the EMA can indeed be realized for inclusion size parameters smaller than a threshold value. The threshold size parameter depends on the refractive-index contrast between the host and inclusion materials and quite often does not exceed several tenths, especially in calculations of the scattering matrix and the absorption cross section. As the inclusion size parameter approaches the threshold value, the scattering-matrix errors of the EMA start to grow with increasing the host size parameter and or the number of inclusions. We confirm, in particular, the existence of the effective-medium regime in the important case of dust aerosols with hematite or air-bubble inclusions, but then the large refractive-index contrast necessitates inclusion size parameters of the order of a few tenths. Irrespective of the highly restricted conditions of applicability of the EMA, our results provide further evidence that the effective-medium regime must be a direct corollary of the macroscopic Maxwell equations under specific assumptions.

aerosols↗
NASA NTRS2016

Ice Cloud Backscatter Study and Comparison with CALIPSO and MODIS Satellite Data

An invariant imbedding T-matrix (II-TM) method is used to calculate the single-scattering properties of 8-column aggregate ice crystals. The II-TM based backscatter values are compared with those calculated by the improved geometric-optics method (IGOM) to refine the backscattering properties of the ice cloud radiative model used in the MODIS Collection 6 cloud optical property product. The integrated attenuated backscatter-to-cloud optical depth (IAB-ICOD) relation is derived from simulations using a CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite) lidar simulator based on a Monte Carlo radiative transfer model. By comparing the simulation results and co-located CALIPSO and MODIS (Moderate Resolution Imaging Spectroradiometer) observations, the non-uniform zonal distribution of ice clouds over ocean is characterized in terms of a mixture of smooth and rough ice particles. The percentage of the smooth particles is approximately 6 percent and 9 percent for tropical and mid-latitude ice clouds, respectively.

MODIS↗
NASA NTRS2001

Microwave Polarized Signatures Generated within Cloud Systems: SSM/I Observations Interpreted with Radiative Transfer Simulations

Special Sensor Microwave /Imager (SSM/I) observations in cloud systems are studied over the tropics. Over optically thick cloud systems, presence of polarized signatures at 37 and 85 GHz is evidenced and analyzed with the help of cloud top temperature and optical thickness extracted from visible and IR satellite observations. Scattering signatures at 85 GHz (TbV(85) less than or = 250 K) are associated with polarization differences greater than or = 6 K, approx. 50%, of the time over ocean and approx. 40% over land. In addition. over thick clouds the polarization difference at 37 GHz is rarely negligible. The polarization differences at 37 and 85 GHz do not stem from the surface but are generated in regions of relatively homogeneous clouds having high liquid water content. To interpret the observations, a radiative transfer model that includes the scattering by non-spherical particles is developed. based on the T-matrix approach and using the doubling and adding method. In addition to handling randomly and perfectly oriented particles, this model can also simulate the effect of partial orientation of the hydrometeors. Microwave brightness temperatures are simulated at SSM/I frequencies and are compared with the observations. Polarization differences of approx. 2 K can be simulated at 37 GHz over a rain layer, even using spherical drops. The polarization difference is larger for oriented non-spherical particles. The 85 GHz simulations are very sensitive to the ice phase of the cloud. Simulations with spherical particles or with randomly oriented non-spherical ice particles cannot replicate the observed polarization differences. However, with partially oriented non-spherical particles, the observed polarized signatures at 85 GHz are explained, and the sensitivity of the scattering characteristics to the particle size, asphericity, and orientation is analyzed. Implications on rain and ice retrievals are discussed.

Prigent, Catherine↗
NASA NTRS2011

Drop Size Distribution Measurements Supporting the NASA Global Precipitation Measurement Mission: Infrastructure and Preliminary Results

Global Precipitation Measurement Mission (GPM) retrieval algorithm validation requires datasets that characterize the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (5 -- 50 km). Key to this process is the combined use of disdrometer and polarimetric radar platforms. Here the disdrometer measurements serve as a reference for up-scaling dual-polarimetric radar observations of the PSD to the much larger volumetric sampling domain of the radar. The PSD observations thus derived provide a much larger data set for assessing DSD variability, and satellite-based precipitation retrieval algorithm assumptions, in all three spatial dimensions for a range of storm types and seasons. As one component of this effort, the GPM Ground Validation program recently acquired five 3rd generation 2D Video disdrometers as part of its Disdrometer and Radar Observations of Precipitation Facility (DROP), currently hosted in northern Alabama by the NASA Marshall Space Flight Center and the University of Alabama in Huntsville. These next-generation 2DVDs were operated and evaluated in different phases of data collection under the scanning domain of the UAH ARMOR C-band dual-polarimetric radar. During this period approximately 7500 minutes of PSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSDs were then used as input to T-matrix code for computation of polarimetric radar moments at C-band. The combined dataset was then analyzed with two basic objectives in mind: 1) the investigation of seasonal variability in the rain PSD parameters as observed by the 2DVDs; 2) the use of combined polarimetric moments and observed gamma distribution parameters in a functional form to retrieve PSD parameters in 4-D using the ARMOR radar for precipitation occurring in different seasons and for different rain system types. Preliminary results suggest that seasonal variations in the DSD parameters do occur, but are most pronounced when comparing tropical PSDs to either winter or summer convective precipitation. For example the previously documented shift to relatively smaller drop diameters in higher number concentrations for equivalent rain rate bins was observed in tropical storm rainbands occurring over Huntsville. On a more inter seasonal basis empirical fits between parameters such as D0 and ZDR do not appear to exhibit robust seasonal biases- i.e., one fit seems to work for all seasons within acceptable standard error (O[10%]) for estimates of D0. In polarimetric retrievals of the vertical variability in PSD (rain layer) for a tropical rainband we find that the Do varies with height when partitioned by specified precipitation categories (e.g., convective or stratiform, heavy and light stratiform etc.) but this variation is of order 10-20% and is smaller than the difference in D0 observed between the basic delineation of convective and stratiform precipitation types. Currently we are expanding our analysis of the vertical structure of the PSD to include several seasonally and/or dynamically-different storm system types (e.g., winter convection and stratiform events; summer mid-latitude convective etc.) sampled by ARMOR. The study will present the results of our combined analyses.

Petersen, Walter A.↗
NASA NTRS2011

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗
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.↗