Lagrangian Numerical Methods for Ocean Biogeochemical Simulations

Francesco Paparella,Marina Popolizio
DOI: https://doi.org/10.1016/j.jcp.2018.01.031
2017-07-25
Abstract:We propose two closely--related Lagrangian numerical methods for the simulation of physical processes involving advection, reaction and diffusion. The methods are intended to be used in settings where the flow is nearly incompressible and the Péclet numbers are so high that resolving all the scales of motion is unfeasible. This is commonplace in ocean flows. Our methods consist in augmenting the method of characteristics, which is suitable for advection--reaction problems, with couplings among nearby particles, producing fluxes that mimic diffusion, or unresolved small-scale transport. The methods conserve mass, obey the maximum principle, and allow to tune the strength of the diffusive terms down to zero, while avoiding unwanted numerical dissipation effects.
Computational Physics
What problem does this paper attempt to address?