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

Multiple and dependent scattering by densely packed discrete spheres: Comparison of radiative transfer and Maxwell theory

Dependent Scattering and Absorption by Densely Packed Discrete Spherical Particles: Effects of Complex Refractive Index

Effect of Dependent Scattering on Light Absorption in Highly Scattering Random Media

Comparison Between the Pseudo-Spectral Time Domain Method and the Discrete Dipole Approximation for Light Scattering Simulations

Dependent scattering criterion modification of disordered dispersed particulate medium with the consideration of particle random distribution and high complex refractive index effects

Radiative Property Investigation of Disperse D Particulate Me Dium with the Consideration of Non-Uniform Particle Size Distribution and Dependent Scattering Effects

Simultaneously Simulating the Scattering Properties of Nonspherical Aerosol Particles with Different Sizes by the Mrtd Scattering Model

Dependent scattering and plasmon coupling in concentrated suspensions of optically hard nanoparticles

A modified recursive transfer matrix algorithm for radiation and scattering computation of multilayer spheres

Dense Media Scattering Based on Numerical Simulations

Light Scattering by Densely Packed Optically Soft Particle Systems, with Consideration of the Particle Agglomeration and Dependent Scattering.

Study of Radiative Transfer in 1D Densely Packed Bed Layer Containing Absorbing–scattering Spherical Particles

Beyond mie theory

Poynting-Stokes tensor and radiative transfer in discrete random media: the microphysical paradigm

ESTIMATION OF RAY EFFECT AND FALSE SCATTERING IN APPROXIMATE SOLUTION METHOD FOR THERMAL RADIATIVE TRANSFER EQUATION

Investigation of Dependent Scattering Criterion for Dispersed Particulate Medium: Effects of Metal Particle Plasmon Resonances and Host Medium Absorption

Microwave emission from snowpacks: modeling the effects of volume scattering, surface scattering and layering

The dependent scattering effect on radiative properties of micro/nanoscale discrete disordered media

Numerical Simulation on Radiative Heat Transfer Through A Two-Dimensional Rectangular Domain with Inhomogeneous, Absorbing and Isotropic/Anisotropic Scattering Media Exposed to Collimated Irradiation

Evaluation of Numerical Scattering in Finite Volume Method for Solving Radiative Transfer Equation by a Central Laser Incidence Model

Electromagnetic Scattering from Multiple Cylinders