WVVerticalDamping

Vertical viscosity and diffusivity


Declaration

WVVerticalDamping < WVForcing

Overview

The damping is designed to mimic the VerticalScalarDiffusivity in Oceananigans to allow for direct comparison between the models. This should probably be used in combination with WVHorizontalScalarDiffusivity. In general, you should be using the WVAdaptiveDamping.

The specific form of the forcing is given by

\[\begin{align} \mathcal{S}_u &= \nu \frac{\partial^2 u}{\partial z^2} \\ \mathcal{S}_v &= \nu \frac{\partial^2 v }{\partial z^2} \\ \mathcal{S}_w &= \nu \frac{\partial^2 w}{\partial z^2} \\ \mathcal{S}_\eta &= \kappa \frac{\partial^2 \eta}{\partial z^2} - \kappa \frac{\partial}{\partial z} \ln N^2 \end{align}\]

with viscosity, \(\nu\), and diffusivity, \(\kappa\). This should be combined with WVHorizontalScalarDiffusivity for a complete closure. For help choosing appropriate values, see the notes in WVAdaptiveDamping.

Usage

Assuming there is a WVTransform instance wvt, to add this forcing,

wvt.addForcing(WVVerticalDamping(wvt,nu=5e-4, kappa=1e-6));

Notes

This is currently implemented in the spatial domain, an is thus highly un-optimized.

Topics



This site uses Just the Docs, a documentation theme for Jekyll.