skip to main content


Title: Vector spherical wave function truncation in the invariant imbedding T-matrix method
Both the computational costs and the accuracy of the invariant-imbedding T-matrix method escalate with increasing the truncation number N at which the expansions of the electromagnetic fields in terms of vector spherical harmonics are truncated. Thus, it becomes important in calculation of the single-scattering optical properties to choose N just large enough to satisfy an appropriate convergence criterion; this N we call the optimal truncation number. We present a new convergence criterion that is based on the scattering phase function rather than on the scattering cross section. For a selection of homogeneous particles that have been used in previous single-scattering studies, we consider how the optimal N may be related to the size parameter, the index of refraction, and particle shape. We investigate a functional form for this relation that generalizes previous formulae involving only size parameter, a form that shows some success in summarizing our computational results. Our results indicate clearly the sensitivity of optimal truncation number to the index of refraction, as well as the difficulty of cleanly separating this dependence from the dependence on particle shape.  more » « less
Award ID(s):
1826936
NSF-PAR ID:
10380070
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Optics Express
Volume:
30
Issue:
17
ISSN:
1094-4087
Page Range / eLocation ID:
30020
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Scattering codes are used to study the optical properties of polar stratospheric clouds (PSCs). Particle backscattering and depolarization coefficients can be computed with available scattering codes once the particle size distribution (PSD) is known and a suitable refractive index is assumed. However, PSCs often appear as external mixtures of supercooled ternary solution (STS) droplets, solid nitric acid trihydrate (NAT) and possibly ice particles, making the assumption of a single refractive index and a single morphology to model the scatterers questionable.Here we consider a set of 15 coincident measurements of PSCs above McMurdo Station, Antarctica, using ground-based lidar, a balloon-borne optical particle counter (OPC) and in situ observations taken by a laser backscattersonde and OPC during four balloon stratospheric flights from Kiruna, Sweden. This unique dataset of microphysical and optical observations allows us to test the performances of optical scattering models when both spherical and aspherical scatterers of different composition and, possibly, shapes are present. We consider particles as STS if their radius is below a certain threshold value Rth and NAT or possibly ice if it is above it. The refractive indices are assumed known from the literature. Mie scattering is used for the STS, assumed spherical. Scattering from NAT particles, considered spheroids of different aspect ratio (AR), is treated with T-matrix results where applicable. The geometric-optics–integral-equation approach is used whenever the particle size parameter is too large to allow for a convergence of the T-matrix method.The parameters Rth and AR of our model have been varied between 0.1 and 2 µm and between 0.3 and 3, respectively, and the calculated backscattering coefficient and depolarization were compared with the observed ones. The best agreement was found for Rth between 0.5 and 0.8 µm and for AR less than 0.55 and greater than 1.5.To further constrain the variability of AR within the identified intervals, we have sought an agreement with the experimental data by varying AR on a case-by-case basis and further optimizing the agreement by a proper choice of AR smaller than 0.55 and greater than 1.5 and Rth within the interval 0.5 and 0.8 µm. The ARs identified in this way cluster around the values 0.5 and 2.5.The comparison of the calculations with the measurements is presented and discussed. The results of this work help to set limits to the variability of the dimensions and asphericity of PSC solid particles, within the limits of applicability of our model based on the T-matrix theory of scattering and on assumptions on a common particle shape in a PSD and a common threshold radius for all the PSDs. 
    more » « less
  2. The index of refraction (n) of particles is an important parameter in optical models that aims to extract particle size and carbon concentrations from light scattering measurements. An inadequate choice ofncan critically affect the characterization and interpretation of optically-derived parameters, including those from satellite-based models which provide the current view of how biogeochemical processes vary over the global ocean. Yet, little is known about hownvaries over time and space to inform such models. Particularly, in situ estimates ofnfor bulk water samples and at diel-resolving time scales are rare. Here, we demonstrate a method to estimatenusing simultaneously and independently collected particulate beam attenuation coefficients, particle size distribution data, and a Mie theory model. We apply this method to surface waters of the North Pacific Subtropical Gyre (NPSG) at hourly resolution. Clear diel cycles innwere observed, marked by minima around local sunrise and maxima around sunset, qualitatively consistent with several laboratory-based estimates ofnfor specific phytoplankton species. A sensitivity analysis showed that the daily oscillation innamplitude was somewhat insensitive to broad variations in method assumptions, ranging from 11.3 ± 4.3% to 16.9 ± 2.9%. Such estimates are crucial for improvement of algorithms that extract the particle size and production from bulk optical measurements, and could potentially help establish a link betweennvariations and changes in cellular composition of in situ particles.

     
    more » « less
  3. Summary

    We explore a truncation error criterion to steer adaptive step length refinement and coarsening in incremental‐iterative path following procedures, applied to problems in large‐deformation structural mechanics. Elaborating on ideas proposed by Bergan and collaborators in the 1970s, we first describe an easily computable scalar stiffness parameter whose sign and rate of change provide reliable information on the local behavior and complexity of the equilibrium path. We then derive a simple scaling law that adaptively adjusts the length of the next step based on the rate of change of the stiffness parameter at previous points on the path. We show that this scaling is equivalent to keeping a local truncation error constant in each step. We demonstrate with numerical examples that our adaptive method follows a path with a significantly reduced number of points compared to an analysis with uniform step length of the same fidelity level. A comparison with Abaqus illustrates that the truncation error criterion effectively concentrates points around the smallest‐scale features of the path, which is generally not possible with automatic incrementation solely based on local convergence properties.

     
    more » « less
  4. The aerosol single scattering albedo (SSA) is the dominant intensive particle parameter determining aerosols direct radiative forcing. For homogeneous spherical particles and a complex refractive index in- dependent of wavelength, the SSA is solely dependent on size parameter (ratio of particle circumference and wavelength) and complex refractive index of the particle. Here, we explore this dependency for the small and large particle limits with size parameters much smaller and much larger than one. We show that in the small particle limit of Rayleigh scattering, a novel, generalized size parameter can be introduced that unifies the SSA dependence on particle size parameter independent of complex refractive index. In the large particle limit, SSA decreases with increasing product of imaginary part of the refractive index and size parameter, another generalized parameter, until this product becomes about one, then stays fairly constant until the imaginary part of the refractive index becomes comparable with the real part minus one. Beyond this point, particles start to acquire metallic character and SSA quickly increases with the imaginary part of the refractive index and approaches one. 
    more » « less
  5. Bleed air is brought into aircraft cabins in order to maintain the quality of the air for passenger and crew health and comfort. The bleed air can be contaminated by oil due to oil seal leaks in the compressor which have been reported randomly and generated significant public interest. Previous studies have measured the particulate size distribution in the bleed air entering the cabin, but never distinguished the type and material of the particulate matter (PM). The particulates could be potentially hazardous oil droplets from the oil seal leaks, water droplets due the presence of fog generated by the cooling system, and so on. In this study we propose a novel technique using light scattering technology to distinguish between contaminant types. This technique uses size and complex index of refraction as the measure. Since each material has a distinct index of refraction, by determining the index of refraction, our proposed low-cost detector could distinguish the compound in the aerosol as well as determine the particle size simultaneously. 
    more » « less