We introduce a description of passive scalar transport based on a (deterministic and hyperbolic) Liouville master equation. Defining a noise term based on time-independent random coefficients, instead of time-dependent stochastic processes, we circumvent the use of stochastic calculus to capture the one-point space–time statistics of solute particles in Lagrangian form deterministically. To find the proper noise term, we solve a closure problem for the first two moments locally in a streamline coordinate system, such that averaging the Liouville equation over the coefficients leads to the Fokker–Planck equation of solute particle locations. This description can be used to trace solute plumes of arbitrary shape, for any Péclet number, and in arbitrarily defined grids, thanks to the time reversibility of hyperbolic systems. In addition to grid flexibility, this approach offers some computational advantages as compared with particle tracking algorithms and grid-based partial differential equation solvers, including reduced computational cost, no Monte-Carlo-type sampling and unconditional stability. We reproduce known analytical results for the case of simple shear flow and extend the description of mixing in a vortex model to consider diffusion radially and nonlinearities in the flow, which govern the long time decay of the maximum concentration. Finally, we validate our formulation by comparing it with Monte Carlo particle tracking simulations in a heterogeneous flow field at the Darcy (continuum) scale.
Domínguez-Vázquez et al. (Thu,) studied this question.