Eulerian simulation of drift trajectories
Convection consists of advection and diffusion.
Diffusion is given by:
The convection equation is (wiki):
with the assumptions that:
the diffusion constant D is constant for the field,
and that the flow u is incompressible (i.e. has no divergence).
the equation simplifies to:
where \(\nabla^2 = \triangle\) is the Laplacian.
Diffusivity (\(m^2/s\)). E.g. between 0.01 and 0.1 for oil on the surface of the ocean (Matsuzakia et. al., 2017).
Decreasing diffusivity places stricter stability criteria on time step.
Porosity, rate of liquid volume to total volume (fraction of flux)
A simple explicit scheme for integrating the convection-equation.
Forward difference in time
ndimage.laplace and np.gradient for spatial differences.