Search NASASearch

Engineering topics

Michael I Mishchenko

Publications and source records attributed to Michael I Mishchenko.

At least 19 records

Scattering By Relatively Small Oblate Spheroidal Drops of Water in the Rainbow Region: T-Matrix Results and Geometric Interpretation

Late in 1999 Michael Mishchenko carried out T-matrix calculations of light scattered by oblate spheroidal drops of water at Marston's request. The purpose of the investigation was to determine if the caustic features, known by that time to be visible in the far-field scattering near the primary rainbow region for horizontally illuminated relatively large drops, would also be present in T-matrix calculations of the scattering by small drops. The axis of rotational symmetry of the drop was taken to be vertical, and the illumination was taken to be horizontal. Scattering into the relevant angular region was evaluated for drop aspect ratios D/H = 1, 1.23, 1.311, and 1.37 where D denotes the equatorial diameter and H the vertical diameter. The results are limited to polarized monochromatic illumination with ka = 140 where a = D/2 and k = 2π/λ with λ being the wavelength of the illumination. As expected from geometrical analysis, a horizontal “V” shaped caustic is visible for D/H = 1.311. The T-matrix results, gray-tone intensity displays provided by Mishchenko, remain instructive.

Scattering by spheroids

Applying Orbital Multi-Angle Photopolarimetric Observations to Study Properties of Aerosols in the Earth's Atmosphere: Implications of Measurements in the 1.378 µm Spectral Channel to Retrieve Microphysical Characteristics and Composition of Stratospheric Aerosols

We analyze the possibilities of orbital photopolarimetric measurements to study properties of aerosols in the Earth's atmosphere. As an example, we consider the case when such measurements are performed within a narrow spectral channel centered at 1.378 µm that allows to retrieve microphysical characteristics of stratospheric aerosols separately from those of tropospheric aerosols. We consider the case of stratospheric aerosols caused by volcanic eruption, and adopt the model of the stratosphere in the form of a homogeneous plane-parallel layer composed of polydisperse spherical particles. We use numerically exact solutions of the vector radiative transfer equation to theoretically simulate measurements carried out at various numbers of scattering angles, including: (i) radiance measurements alone; (ii) polarization measurements alone; and (iii) radiance and polarization measurements together. The results of computations show that the simultaneous use of radiance and polarization measurements at a sufficiently large number of scattering angles enables one to retrieve the optical thickness, effective radius, and refractive index of aerosols with adequate accuracy. We demonstrate how the accuracy of the derived values of the optical parameters of aerosols depends on the accuracy of measurements of the intensity and polarization of the reflected light, optical thickness of aerosol layer itself, effective radius of aerosols, width of the particle size distribution, and number of viewing angles.

Aerosols

Spectral Spherical Harmonics Discrete Ordinate Method

A new method for modeling the radiative transfer in inhomogeneous three-dimensional media illuminated by a Gaussian beam is described. This approach, called the Spectral Spherical Harmonics Discrete Ordinate Method (SSHDOM), uses the Fourier expansion method to transform the three-dimensional radiative transfer into an one-dimensional equation in the spectral domain, and the Spherical Harmonics Discrete Ordinate Method (SHDOM) for its solution. Specifically, (i) the source function is represented in the spectral domain through a spherical harmonic expansion, (ii) the spectral one-dimensional radiative transfer equation is integrated along discrete ordinates through a spatial grid, and (iii) the solution method is based on the Picard iteration. Both SSHDOM and SHDOM algorithms are implemented in a common computer code.

Gaussian beam

Electromagnetic scattering by discrete random media illuminated by a Gaussian beam I: Derivation of the radiative transfer equation

In this paper we present the vector radiative transfer theory for a discrete random medium illuminated by a Gaussian beam. The analysis is based on a plane wave spectrum representation for a Gaussian beam and uses an approach developed previously for a discrete random medium illuminated by a plane electromagnetic wave. Specifically, we establish an integral representation for the coherent field, define an approximate coherent field that satisfies the differential equation fulfilled by the coherent field corresponding to a plane electromagnetic wave and matches the Gaussian beam at the interface of the particulate medium, and finally, derive the vector radiative transfer equation. For weakly focused Gaussian beams, the resulting equation is the traditional radiative transfer equation

Gaussian beam

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

An Overview of the Null-Field Method. I. Formulation and Basic Results

In this paper we revisit the fundamentals of the null-field method with discrete sources. We prove the unique solvability of the null-field equations for the total field inside the particle and the internal field outside the particle, at all wavenumbers. For this purpose, we use the equivalence between the null-field and surface integral equations methods. Furthermore, we discuss the completeness property of different systems of discrete sources for a surface field approximation, and derive an infinite set of integral equations for the surface fields in a variety of discrete sources. Finally, we formulate the null-field scheme as an approach aiming to construct an approximate solution to the scattering problem. The way in which we introduce the matrix of a particle is different from the standard approach relying on the assumption that the system of regular vector spherical wave functions for the interior problem is a basis.

Null-field method

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method

Comprehensive Thematic T-Matrix Reference Database: a 2017-2019 Update

Since its inception in the mid-1960s, the T-matrix method has been one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It has been widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This compilation is the ninth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2017. It also includes a few earlier publications that have so far been overlooked.

Michael I Mishchenko

Co- and Counter-Propagating Wave Effects in an Absorbing Medium

In this semi-tutorial paper, we revisit the interference phenomena caused by pairs of copropagating or counter-propagating transverse electromagnetic waves by letting the host medium be absorbing. We first consider plane waves in an unbounded medium, summarize the standingwave solution of the Maxwell equations, and discuss specific effects caused by nonvanishing absorption. We then consider the superposition of plane and spherical waves in the context of fartfield electromagnetic scattering by a particle. To this end we modify the classical Jones lemma by allowing nonzero absorption in the host medium and consider its most obvious consequences such as forward- and backscattering interference. The physical similarity of the two scenarios (superpositions of plane waves and superpositions of plane and spherical waves) is discussed.

Interference

Electromagnetic Scattering by Discrete Random Media. IV: Coherent Backscattering

The problem of backscattering of light by a discrete random medium illuminated by an obliquely incident plane electromagnetic wave is considered.The analysis is performed in a linear-polarization basis and includes a complete derivation of the cross reflection matrix for a layer with densely and sparsely distributed particles, the design of an approximate method for computing the ladder and cross reflection matrices in the case of a semi-infinite medium with a sparse distribution of particles, the derivation of the relations between the elements of the ladder and cross reflection matrices in the exact backscattering direction for dense and sparse media, and the development of practical algorithms for solving the underlying integral equations by the method of Picard iterations and the discrete ordinate method. Simulation results for particles with large size parameters are also presented.

Adrian Doicu

Electromagnetic Scattering by Discrete Random Media. III: The Vector Radiative Transfer Equation

A vector radiative transfer equation with an additional source term typical of dense media is obtained. The analysis includes (i) the derivation of an integral equation for the correlation matrix of the exciting field coefficients accounting for the correlation between the particles, (ii) the derivation of an integral representation for the specific coherency dyadicin terms of this matrix, and (iii) the simplification of the integral equation for the correlation matrix and of the integral representation for the specific coherency dyadic by employing a series of approximations which are characteristic of sparse media.

Adrian Doicu

Electromagnetic Scattering by Discrete Random Media. II: The Coherent Field

The computation of the coherent field in the case of a plane electromagnetic wave obliquely incident on a discrete random layer with non-scattering boundaries is addressed. For dense media, the analysis is based on a special-form solution for the conditional configuration-averaged exciting field coefficients, and is restricted to the computation of the so-called zeroth-order fields without a special treatment of the boundary regions. In this setting, we calculate the coherent fields reflected and transmitted by the layer, and the coherent field inside the layer. We found that these fields are analytically equivalent to plane electromagnetic waves, and investigated the fulfillment of the boundary conditions for the electric fields at the layer interfaces. The results are then particularized to the cases of normal incidence and a semi-infinite discrete random medium. For sparsely distributed particles, we present a self-consistent derivation of the coherent field and discuss the Twersky and Foldy approximations.

Adrian Doicu

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu

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

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;

Electromagnetic Scattering and Emission by a Fixed Multi-Particle Object in Local Thermal Equilibrium: General Formalism

The majority of previous studies of the interaction of individual particles and multi-particle groups with electromagnetic field have focused on either elastic scattering in the presence of an external field or self-emission of electromagnetic radiation. In this paper we apply semi-classical fluctuational electrodynamics to address the ubiquitous scenario wherein a fixed particle or a fixed multi-particle group is exposed to an external quasi-polychromatic electromagnetic field as well as thermally emits its own electromagnetic radiation. We summarize the main relevant axioms of fluctuational electrodynamics, formulate in maximally rigorous mathematical terms the general scattering emission problem for a fixed object, and derive such fundamental corollaries as the scattering emission volume integral equation, the Lippmann Schwinger equation for the dyadic transition operator, the multi-particle scattering emission equations, and the far-field limit. We show that in the framework of fluctuational electrodynamics, the computation of the self-emitted component of the total field is completely separated from that of the elastically scattered field. The same is true of the computation of the emitted and elastically scattered components of quadratic linear forms in the total electromagnetic field. These results pave the way to the practical computation of relevant optical observables.

Electromagnetic scattering; Thermal emission; Fluc

First-Principles Modeling Of Electromagnetic Scattering By Discrete and Discretely Heterogeneous Random Media

A discrete random medium is an object in the form of a finite volume of a vacuum or a homogeneous material medium filled with quasi-randomly and quasi-uniformly distributed discrete macroscopic impurities called small particles. Such objects are ubiquitous in natural and artificial environments. They are often characterized by analyzing theoretically the results of laboratory, in situ, or remote-sensing measurements of the scattering of light and other electromagnetic radiation. Electromagnetic scattering and absorption by particles can also affect the energy budget of a discrete random medium and hence various ambient physical and chemical processes. In either case electromagnetic scattering must be modeled in terms of appropriate optical observables, i.e., quadratic or bilinear forms in the field that quantify the reading of a relevant optical instrument or the electromagnetic energy budget. It is generally believed that time-harmonic Maxwell's equations can accurately describe elastic electromagnetic scattering by macroscopic particulate media that change in time much more slowly than the incident electromagnetic field. However, direct solutions of these equations for discrete random media had been impracticable until quite recently. This has led to a widespread use of various phenomenological approaches in situations when their very applicability can be questioned. Recently, however, a new branch of physical optics has emerged wherein electromagnetic scattering by discrete and discretely heterogeneous random media is modeled directly by using analytical or numerically exact computer solutions of the Maxwell equations. Therefore, the main objective of this Report is to formulate the general theoretical framework of electromagnetic scattering by discrete random media rooted in the Maxwell- Lorentz electromagnetics and discuss its immediate analytical and numerical consequences. Starting from the microscopic Maxwell-Lorentz equations, we trace the development of the first principles formalism enabling accurate calculations of monochromatic and quasi-monochromatic scattering by static and randomly varying multiparticle groups. We illustrate how this general framework can be coupled with state-of-the-art computer solvers of the Maxwell equations and applied to direct modeling of electromagnetic scattering by representative random multi-particle groups with arbitrary packing densities. This first-principles modeling yields general physical insights unavailable with phenomenological approaches. We discuss how the first-order-scattering approximation, the radiative transfer theory, and the theory of weak localization of electromagnetic waves can be derived as immediate corollaries of the Maxwell equations for very specific and well-defined kinds of particulate medium. These recent developments confirm the mesoscopic origin of the radiative transfer, weak localization, and effective-medium regimes and help evaluate the numerical accuracy of widely used approximate modeling methodologies.

Michael I Mishchenko

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