Transport in Random Media With Inhomogeneous Mixing Statistics
The Levermore-Pomraning (LP) model has been investigated for use in implicit Monte Carlo (IMC) algorithms to describe thermal photon transport computations in stochastically mixed cells with material energy coupling. For purely absorbing non-emitting materials, the LP-model is known to predict the exact material averaged fluxes when the mixing statistics are Markovian. In the presence of scattering or reemmission, or if the mixing is non-Markovian, the LP-model is approximate and has greatly diminished accuracy for optically thick, highly scattering materials. Nevertheless, the model enables a robust extension of the simple atomic mix approximation to accommodate spatial correlation effects through well characterized parameters such as mean chord lengths. In this work, we construct exact analytical solutions of the LP-equations for piecewise constant spatial variation of the material chord lengths. It is not possible in general to obtain closed form solutions when the chord lengths are continuous functions of space, although we briefly examine special conditions under which this is feasible. However, by allowing an arbitrarily fine spatial segmentation of the medium it is possible to construct highly refined solutions to fit desired spatial chord length profiles. Material volume fractions that are consistent with the chord lengths are first constructed by solving the associated Chapmann-Kolmogorov equations and applied to obtain the material averaged fluxes in two settings: when the materials are purely absorbing, and when scattering is present in rod model geometry. Illustrative numerical results for the volume fractions and material averaged fluxes are obtained by implementing the algorithms in MATLAB for spatially linear and quadratic variations of the chord lengths.