Processing math: 100%
Neko 0.9.99
A portable framework for high-order spectral element flow simulations
All Classes Namespaces Files Functions Variables Typedefs Enumerator Macros Pages
Governing equations

Fluid

When possible, Neko integrates the Navier-Stokes equations in time formulated as follows:

\rho \left( \frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j} \right) = -\frac{\partial p}{\partial x_i} + \mu \frac{\partial u_i }{\partial x_j \partial x_j} + \rho \sum_j f^u_{i,j} , \quad i=1,2,3,

together with the continuity equation

\frac{\partial u_i}{\partial x_i} = 0 \ .

Here, u_i are the components of the velocity vector, p is the pressure, \rho is the density, \mu is the dynamic viscosity, and f^u_{i,j} are the components of the source terms. Both \rho and \mu are usually treated as constants.

When the viscosity is varying in space (e.g. for LES modelling), the viscous stress tensor cannot be simplified and the following equations are solved in a coupled manner.

\rho \left( \frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j} \right) = -\frac{\partial p}{\partial x_i} + \frac{\partial}{\partial x_j} \left(\mu_{tot} \left(\frac{\partial u_i }{\partial x_j} + \frac{\partial u_j }{\partial x_i} \right) \right) + \rho \sum_j f^u_{i,j} , \quad i=1,2,3.

Here, \mu_{tot} is the total viscosity field, potentially including the contribution of turbulence modelling.

Scalar

Optionally, and additional equation for scalar transport can be solved. Here we define it as an equation for temperature, but the physical meaning can of course differ from case to case.

\rho c_p \left( \frac{\partial T}{\partial t} + u_j \frac{\partial T}{\partial x_j} \right) = \frac{\partial}{ \partial x_j} \left(\lambda_{tot} \frac{\partial T }{ \partial x_j} \right)+ \rho c_p \sum_j f_j^s.

Here, T is the scalar temperature field, c_p is the specific heat capacity, \lambda is the total thermal conductivity, and f_j^s is the active source terms.

Non-dimensionalisation

A non-dimensional for of the Navier-Stokes equations may be found by defining the Reynolds number Re = \frac{\rho U L}{\mu} \ , with reference velocity U, reference length L, density \rho and dynamic viscosity \mu; the latter two may be combined into the kinematic viscosity \nu=\mu/\rho. In order to get a properly scaled problem, one can define the mesh such that the reference length is unity, and set the velocity boundary conditions (or forcing) such that the reference velocity is unity. Then specifying the Reynolds number (e.g. via the case file) will set the (now non-dimensional) density \rho=1 and the non-dimensional viscosity to \nu = 1/Re. Similarly for a scalar, a Prandtl number Pr may be specified.