Personal tools

Difference between revisions of "Biharmonic"

From MohidWiki

Jump to: navigation, search
Line 3: Line 3:
 
<math>\mathrm{J} = -\mathrm{K}\cdot \nabla C + \nabla_{h}\left( \kappa_4 \nabla_{h}^2 C \right) </math>
 
<math>\mathrm{J} = -\mathrm{K}\cdot \nabla C + \nabla_{h}\left( \kappa_4 \nabla_{h}^2 C \right) </math>
  
A ready made formula that would estimate the order of magnitude of the k_4 coefficient would be
+
A ready made formula that would estimate the order of magnitude of the <math>\kappa_4</math> coefficient would be
  
 
<math>\kappa_4 \sim \frac{U}{16} \Delta x^3</math>
 
<math>\kappa_4 \sim \frac{U}{16} \Delta x^3</math>

Revision as of 14:21, 28 June 2011

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.