Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Version 5 (modified by Giersch, 9 years ago) (diff)

--

This site is currently under construction!

Governing equations

By default, PALM has six prognostic quantities: the velocity components u, v, w on a Cartesian grid, the potential temperature θ, specific humidity qv or a passive scalar s, and the SGS turbulent kinetic energy (SGS-TKE) e. The separation of resolved scales and SGS is implicitly achieved by averaging the governing equations over discrete Cartesian grid volumes as proposed by Schumann (1975). Moreover, it is possible to run PALM in a direct numerical simulation mode by switching off the prognostic equation for the SGS-TKE and setting a constant eddy diffusivity. For a list of all symbols and parameters, that we will introduce in this section, see the following Table 1 and 2.

List of genereal model parameters List of genereal symbols

The model is based on the non-hydrostatic, filtered, incompressible Navier-Stokes equations in Boussinesq-approximated form. In the following set of equations, angle brackets denote a horizontal domain average. A subscript 0 indicates a surface value. Note that the variables in the equations are implicitly filtered by the discretization (see above), but that the continuous form of the equations is used here for convenience. A double prime indicates SGS variables. The overbar indicating filtered quantities is omitted for readability, except for the SGS flux terms. The equations for the conservation of mass, energy and moisture, filtered over a grid volume on a Cartesian grid, then read as

\begin{align*} \frac{\partial u_i}{\partial t}&= - \frac{\partial u_i u_j}{\partial x_j} -\varepsilon_{ijk}f_j u_k + \varepsilon_{i3j}f_3 {u_{\mathrm{g},j}} - \frac{1}{\rho_0} \frac{\partial \pi^\ast}{\partial x_i} + g \frac{\theta_\mathrm{v} - \langle\theta_{\mathrm{v}}\rangle}{\langle\theta_{\mathrm{v}}\rangle}\delta_{i3}-\frac{\partial}{\partial x_j} \left(\overline{u_i^{\prime\prime} u_j^{\prime\prime}} - \frac{2}{3}e\delta_{ij}\right), \\ \frac{\partial u_j}{\partial x_j}&=0, \\ \frac{\partial \theta}{\partial t} &= - \frac{\partial u_j \theta}{\partial x_j} -\frac{\partial}{\partial x_j}\left(\overline{u_j^{\prime\prime}\theta^{\prime\prime}}\right) - \frac{L_\mathrm{V}}{c_p \Pi} \Psi_{q_\mathrm{v}}, \\ \frac{\partial q_\mathrm{v}}{\partial t} &= - \frac{\partial u_j q_\mathrm{v}}{\partial x_j} - \frac{\partial}{\partial x_j}\left(\overline{u_j^{\prime\prime}q^{\prime\prime}_\mathrm{v}}\right) + \Psi_{q_\mathrm{v}},\\ \frac{\partial s}{\partial t} &= - \frac{\partial u_j s}{\partial x_j} - \frac{\partial}{\partial x_j}\left(\overline{u_j^{\prime\prime}s^{\prime\prime}}\right) + \Psi_s. \end{align*}

Here, i, j, k ∈ {1, 2, 3}. ui are the velocity components (u1 = u, u2 = v, u3 = w) with location xi (x1 = x, x2 = y, x3 = z), t is time, fi = (0, 2 Ω cos(φ), 2 Ω sin(φ)) is the Coriolis parameter with Ω being the Earth's angular velocity and φ being the geographical latitude. ug,k are the geostrophic wind speed components, ρ0 is the density of dry air, π∗ = p∗ + 2/3 ρ0 e is the modified perturbation pressure with p∗ being the perturbation pressure and the SGS-TKE e = 1/2 u″iu″i

References

  • Schumann U. 1975. Subgrid scale model for finite difference simulations of turbulent flows in plane channels and annuli. J. Comput. Phys. 18: 376–404.

Attachments (2)

Download all attachments as: .zip

                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |