Personal tools

Biharmonic

From MohidWiki

Revision as of 15:21, 28 June 2011 by Guillaume (talk | contribs)
Jump to: navigation, search

The biharmonic coefficient is an extra term that is added to the horizontal turbulent viscous flux, J:

\mathrm{J} = -\mathrm{K}\cdot \nabla C + \nabla_{h}\left( \kappa_4 \nabla_{h}^2 C \right)

A ready made formula that would estimate the order of magnitude of the \kappa_4 coefficient would be

\kappa_4 \sim \frac{U}{16} \Delta x^3

where U is a typical flow velocity and \Delta x is the numerical model horizontal resolution.

MOHID

In mohid, you need to add in the Hydrodynamic_x.dat the following keywords:

BIHARMONIC        : [0/1] <-- Activates de biharmonic viscous flux
BIHARMONIC_COEF   : [Real] <-- Sets the value of \kappa_4

References

Delhez, E., & Deleersnijder, E. (2007). Overshootings and spurious oscillations caused by biharmonic mixing. 
Ocean Modelling, 17(3), 183-198. doi: 10.1016/j.ocemod.2007.01.002.