Monte Carlo propagation in a fluctuating layer

A Monte Carlo approach is developed to study the characteristics of a propagating wave through a randomly inhomogeneous layer. An exact analysis of such problems is not generally available. Previous computer modelings have provided a partial solution, giving only the reflection and transmission coefficients. In contrast, the present simulation approach yields a complete wave solution throughout the nodes, P, of the layer. The technique adopts an efficient numerical procedure for the wave solution across each member of the ensemble, Q. Had a Gauss elimination been used, the increase by a factor, QP, in solution time would have been forbidding. Furthermore, the present implementation extends the simulation to physical circumstances not previously considered. It accounts for energy leakage near a vertexing ray. Also, it yields correct statistical descriptors of the wave field with the development of functional equations in stochastic variables, prior to the application of any averaging operation.