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 ρ ρ φ ρ ρ − − ∞

Uniaxial Effective Permittivity Extracted from Anisotropic Bicontinuous Media Using Numerical Solution of Maxwell Equation in 3d and Strong Permittivity Fluctuations

Modeling Snow Anisotropy and Backscattering Co-Polarization Phase Difference Using Bicontinuous Media and Numerical Solutions of Maxwell Equations

Efficient Method for Scattering from Cylindrical Components of Vegetation and Its Potential Application to the Determination of Effective Permittivity

Effective Permittivity and Scattering of Bicontinuous Random Medium with Strong Permittivity Fluctuation Theory

Uniaxial Effective Permittivity of Anisotropic Bicontinuous Random Media Using NMM3D.

MIXTURE EFFECTIVE PERMITTIVITY SIMULATIONS USING IMLMQRF METHOD ON PRECONDITIONED EFIE

The Simulation Of Ice Anisotropy Feature Based On 3-D Fdtd Method

Numerical Simulation of Maxwell's Equation in 3D (NMM3D) Applied to Active and Passive Remote Sensing of Terrestrial Snow and Snow on Sea Ice

Efficient Implementation of the Invariant Imbedding T-matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Scattering and Absorption of Light by Ice Particles: Solution by a New Physical-Geometric Optics Hybrid Method

Scattering and Emission Models for Microwave Remote Sensing of Snow Using Numerical Solutions of Maxwell Equations

Microwave Signatures of Snow Cover Using Numerical Maxwell Equations Based on Discrete Dipole Approximation in Bicontinuous Media and Half-Space Dyadic Green's Function.

Retrieving Material Permittivity Under Oblique Incidence

Rough Surface and Volume Scattering of Soil Surfaces, Ocean Surfaces, Snow, and Vegetation Based on Numerical Maxwell Model of 3-D Simulations.

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

CHARACTERISTICS OF SNOW SCATTERING WITH BICONTINUOUS MEDIA DISCRETE DIPOLE APPROXIMATION AND MEMLS

Full Wave Simulation of Snowpack Applied to Microwave Remote Sensing of Sea Ice

Stochastic Electromagnetic Scattering From Stratified Media of Inhomogeneous Layers With Random Constitutive Parameters

Anisotropic homogenized composite mediums arising from truncated spheres, spheroids, and ellipsoids

Dense Media Scattering Based on Numerical Simulations

Scattering Characteristics of an Anisotropic Medium Sphere