A Fast Explicit Operator Splitting Method for Passive Scalar Advection

The dispersal and mixing of scalar quantities such as concentrations or thermal energy are often modeled by advection-diffusion equations. Such problems arise in a wide variety of engineering, ecological and geophysical applications. In these situations a quantity such as chemical or pollutant concentration or temperature variation diffuses while being transported by the governing flow. In the passive scalar case, this flow prescribed and unaffected by the scalar. Both steady laminar and complex (chaotic, turbulent or random) time-dependent flows are of interest and such systems naturally lead to questions about the effectiveness of the stirring to disperse and mix the scalar. The development of reliable numerical methods for advection-diffusion equations is crucial for understanding their properties, both physical and mathematical. In this paper, we extend a fast explicit operator splitting method, recently proposed in (A. Chertock, A. Kurganov, G. Petrova, Int. J. Numer. Methods Fluids 59:309–332, 2009), for solving deterministic convection-diffusion equations, to the problems with random velocity fields and singular source terms. A superb performance of the method is demonstrated on several two-dimensional examples.