Neko
0.8.0rc2
A portable framework for highorder spectral element flow simulations

Neko's solves the NavierStokes equations 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.$$
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 treated as constants.
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 s}{\partial t} + u_j \frac{\partial s}{\partial x_j} \right) = \lambda \frac{\partial^2 s }{\partial x_j \partial x_j} + \rho c_p \sum_j f_j^s.$$
Here, \( s \) is the scalar field, \( c_p \) is the specific heat capacity, \( \lambda \) is the thermal conductivity, and \( f_j^s \) is the active source terms. Both \( c_p \) and \( \lambda \) are treated as constants.