Engineering topics

Adrian Doicu

Publications and source records attributed to Adrian Doicu.

11 records

NASA NTRS2020

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↗

NASA NTRS2020

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↗

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↗

NASA NTRS2020

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↗

NASA NTRS2020

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↗

NASA NTRS2019

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↗

NASA NTRS2019

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↗

NASA NTRS2019

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↗

NASA NTRS2019

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↗

NASA NTRS2018

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗

NASA NTRS2018

Radiative Transfer in a Discrete Random Medium Adjacent to a Half-Space with a Rough Interface

For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.

Adrian Doicu↗