.. _doc.basic.parametrization: Parametrization and closures =============================== .. _turbulence: Turbulence models ----------------- Three classes of turbulent closure models have been included in the MARS modelling system. * Numerous models based on **algebraic relations** exist. The following formulations are actually implemented in the MARS model code. * **Constants coefficient** | Regarding nuV and kapaV , in the MARS modelling system we recommend a constant viscosity | between 10e−4 (kinematic viscosity) and 10e−1 m2/s (the largest value measured in situ) * **Prandtl’s law** | Prandtl derived the logarithmic velocity distribution near a solid surface by using a mixing length theory * **Quetin** | based on density stratification * **Pacanowski & Philander** * **Single equation models** | With this approach, the eddy viscosity (resp. the diffusivity) is set as a function of a turbulent velocity V and a mixing length lm * **Turbulent Kinetic Energy** * **Mixing length** * **Two equations models** * **Horizontal diffusion** Air_sea interactions -------------------- The thermal structure of the upper ocean layer is quite complex due to the time and space variability of wind stress and heat fluxes. The heat fluxes at the sea surface, combined with the wind, induced turbulent mixing near surface layer cause temperature variations within the upper layers and consequently, density stratification. .. _vertical_obc: Vertical boundary conditions ---------------------------- The flux at the upper and lower boundaries of the domain together with kinematic conditions leads to the following sets of equations : * **at surface** (:math:`\sigma=0`) +--------------------------------------------------------------------------------------------------------------------+----------------------------------------------------------------------+ | * :math:`\frac{nz}{D}\frac{nz}{D}\frac{\partial u}{\partial \sigma }=\frac{\tau_{sx}}{\rho_0}` | | D= height of water column, D = H + :math:`\xi(x,y)` | | * :math:`\frac{nz}{D}\frac{nz}{D}\frac{\partial v}{\partial \sigma }=\frac{\tau_{sy}}{\rho_0}` | | nz=vertical eddy viscosity ; kz=vertical eddy diffusivity | | * :math:`\frac{kz}{D}\frac{nz}{D}\frac{\partial T}{\partial \sigma }=\frac{Q_T}{\rho_0 C_p}` | | u,v = zonal, meridian velocity components | | * :math:`kz \frac{\partial S}{\partial \sigma }=0` | | :math:`\tau_{sx},\tau_{sy}` = surface stress components | | * :math:`w^{\stackrel{0}{a}}=0` = vertical velocity components | | :math:`\rho_0` = water density ; :math:`\rho_a` = air density | | | | | T,S = Temperature & Salinity | | | | | Q_T = heat flux at the air_sea interface | | * :math:`(\tau_{sx}` , :math:`\tau_{sy})=\rho_a Cd_s \parallel \overrightarrow{W} \parallel (W_x` , :math:`W_y)` | | Cd_s = surface drag coefficient | | | | | :math:`(W_x` , :math:`W_y)` = wind velocity vector | | | | | at 10 m above the sea surface | +--------------------------------------------------------------------------------------------------------------------+----------------------------------------------------------------------+ * **at bottom** (:math:`\sigma=-1`) +--------------------------------------------------------------------------------------------------------------------+----------------------------------------------------------------------+ | * :math:`\frac{nz}{D}\frac{nz}{D}\frac{\partial u}{\partial \sigma }=\frac{\tau_{bx}}{\rho_0}` | | D= height of water column, D = H + :math:`\xi(x,y)` | | * :math:`\frac{nz}{D}\frac{nz}{D}\frac{\partial v}{\partial \sigma }=\frac{\tau_{by}}{\rho_0}` | | nz=vertical eddy viscosity ; kz=vertical eddy diffusivity | | * :math:`\frac{kz}{D}\frac{nz}{D}\frac{\partial T}{\partial \sigma }=0` | | u,v = zonal, meridian velocity components | | * :math:`kz \frac{\partial S}{\partial \sigma }=0` | | :math:`\tau_{bx},\tau_{by}` = bottom stress components | | * :math:`w^{\stackrel{0}{a}}=0` | | :math:`\rho_0` = water density ; | | * :math:`(\tau_{bx}` , :math:`\tau_{by})=\rho_a Cd_B \parallel \overrightarrow{u} \parallel(u` , :math:`v)` | | T,S = Temperature & Salinity | | * :math:`Cd_B=(\frac{\kappa}{ln{\frac{z+h+z_0}{z_0}}})^2` | | :math:`\kappa` =Von Karman constant ; z_0 = bed roughness | +--------------------------------------------------------------------------------------------------------------------+----------------------------------------------------------------------+ Open boundary conditions ------------------------ Numerical models of finite coastal regions inevitably involve the treatment of open boundaries where the numerical grid ends. Ideal open boundaries are ”transparent” to motions which are generated within the computational domain and which are ”felt” at the open boundaries. There is no numerical treatment of open boundaries which achieves this ideal. A large number of open boundary conditions have been proposed in the past (Palma & Matano, 1998 [29]). Due to the essentially hyperbolic nature of the incompressible, hydrostatic primitive equations, external data is required for only inward boundary fluxes. The outward fluxes may be treated with an algorithm for two-dimensional radiation. In the MARS model the prescription of the Open Boundary Conditions (OBCs) is based on the method of characteristics. .. image:: FIG/PHYS/schema_discretisation_obc.jpg Moving boundaries ------------------ In shallow coastal areas, some areas may dry up at low water level, which means that we have moving boundaries. The definition of these boundaries would appear to be rather straightforward : the instantaneous water depth h = 0 (h is defined by eq. (1.31)). Several techniques have been developed to track the moving boundaries. However, there are some difficulties involved in applicable numerical models. In the MARS model, the wetting and drying processes have been cleared up by introducing a water depth H on each side of a grid cell .. _forcing: Forcings ------------------ * Momentum equations * Winds * Tides * River discharges (mass) * Thermodynamics * Heat fluxes (LW, SH, latent, sensible) * River discharges (salinity, temperature)