Dense Media Scattering Based on Numerical Simulations
L. Tsang,C. H. Chan,J. Guo,B.-K. Huang,M. Y. Xia,C-T Chen
2002-01-01
Abstract:The frequency and polarization dependence of scattering by the densely packed spherical particles is rigorously studied with the Monte Carlo simulations for solving both 2and 3-dimensional Maxwell's equations at microwave frequencies. Numerical simulation solutions indicate that the frequency dependence of densely packed sticky small particles is much weaker than that of independent scattering while the cross polarization is stronger than that of independent scattering. The typical simulation results are illustrated with the snow parameters. To facilitate the Monte Carlo simulations, the Sparse-Matrix Canonical-Grid (SMCG) method has been applied to the simulation of scattering by 2D small particles. We also study the simulations of effective permittivity properties of foam. INTRODUCTION In the volume scattering of geophysical media, wave scattering from the dense media is an important issue [1][2]. Dense media consists of particles densely packed together. Traditional independent scattering theory states that the frequency dependence of the scattering is proportional to the fourth power in 3-D scattering and frequency squared in 2D scattering for the small particles and the scattering of cross polarization by spherical particles is small. With the increase of particle density, the issue to address is whether this will change the behaviors as predicted by classical scattering theory. In this paper, we study rigorously the frequency and polarization dependence of scattering with the Monte Carlo simulations of the solutions of both 2and 3-dimensional Maxwell's equations. We use the deposition and bonding technique to generate the positions of particles and then solve the Maxwell’s equations that include all multiple scattering. The extinction, scattering and absorption properties of dense media are calculated for dense media of sticky and non-sticky particles. The T-matrix method is used. We use the internal field formulation of Foldy-Lax equations that have been shown with much better condition numbers for the matrix equations [2]. It is particularly advantageous for iterative solution of the Foldy-Lax equation. Numerical simulation solutions of Maxwell's equations indicate that the frequency dependence of densely packed sticky small particles is much weaker than that of independent scattering. The results also show strong cross polarization. The simulations are performed with the snow parameters at microwave frequencies from 18GHz to 90GHz. Comparisons are made with experimental data of dry snow. The simulations are able to account for the observed frequency dependence that is much weaker than frequency to the fourth power. To facilitate the Monte Carlo solution, the fast computational method of the sparse matrix canonical grid method (SMCG) is also developed in the paper for cylinder cases. In the SMCG, the impedance matrix elements corresponding to strong interactions among nearby scatterers are computed and stored as a sparse matrix. The sparse matrix-vector multiplication in the iterative solution is computed directly in the conventional manner. Interaction between two cylinders that were far apart was computed indirectly through the use of the addition theorem and a uniformly spaced canonical grid. The memory requirement of the method is O(N) and the computational complexity is O(N log N) as the fast Fourier transforms were used to compute interactions among all the canonical grid points simultaneously. 3-D FORMULATION Consider a medium consisting of N spheres with radii a and permittivities p ε centered at 1 r , 2 r , ... and N r . An incident wave with wave vector i k impinging onto the medium. The solution of Maxwell's equations can be cast into the Foldy-Lax multiple scattering equations of matrix form: ∑ ≠ = ⋅ + = N q p p inc q i p p p q q a r k i w T r r k w 1 ) ( ) ( ) ( ) exp( ) ( σ (1) with N q ,..., 2 , 1 = . Vector ) (q w is the exciting field coefficients, ) ( p q r r k σ is the vector spherical wave transformation matrix, and ) ( p T is the T-matrix representing scattering from the scatterer p . inc a is the incident field coefficients. The physical interpretation of the above is the exciting field at particle q is equal to the incident field plus the scattered field from all other particles except itself. The above matrix equation can be solved with Gaussian elimination method or iterative techniques for the exciting field coefficients ) (q w and then they can be used to calculate the scattering fields from the scatterers [2]. The computational complexity of iterative approach for matrix equation is O( 2 N ). However with the sparse matrix canonical grid (SMCG) method, the computational complexity can be reduced to O( N N log ). This is done for 2-D Maxwell’s equation described as follows. 2-D SMCG The SMCG [3] method decomposes the all interactions into strong interaction (near field) part and weak interaction (non-near field) part. For strong interaction part, field solution are calculated directly based on the Foldy-Lax equations. For the weak interactions, which account for most computational time, the computation is done indirectly with cylinder translation to the canonical grid points by using the translational addition theorem. The formulation for the non-near field interactions with the SMCG is vi in n k i q N n E k e i e w i q i φ ρ − • − = ) )( ( ) )( ( ' ) ( ' ) ' ( ' ) ( ) 1 ( ' ) ( p N n N n n n i q p n N q p p S p n n w T e k H q p ρ ρ φ ρ ρ − − ∞