freeflow/rans/twoeq/komega/model.hh File Reference

A single-phase, isothermal k-omega 2-Eq. model. More...

#include <dumux/common/properties.hh>
#include <dumux/freeflow/properties.hh>
#include <dumux/freeflow/rans/model.hh>
#include <dumux/freeflow/rans/twoeq/indices.hh>
#include <dumux/freeflow/turbulencemodel.hh>
#include "problem.hh"
#include "fluxvariables.hh"
#include "localresidual.hh"
#include "volumevariables.hh"
#include "iofields.hh"

Go to the source code of this file.

Description

A single-phase, isothermal k-omega 2-Eq. model.

Single-phase Reynolds-Averaged Navier-Stokes flow For a detailed model decription see freeflow/rans/model.hh.

Two additional PDEs, one for the turbulentKineticEnergy (k) and a second for the dissipation (omega) are used to calculate the eddy viscosity for this model. The model is taken from Wilcox, 2008 [70].

Turbulent Kinetic Energy balance:

$$\frac{\partial \left( k \right)}{\partial t} + \nabla \cdot \left( \mathbf{v} k \right) - \nabla \cdot \left[ \left( \nu + \sigma_\textrm{k} \nu_\textrm{t} \right) \nabla k \right] - P + \beta_k^{*} k \omega = 0$$

with $P = 2 \nu_\textrm{t} \mathbf{S} \cdot \mathbf{S}$ and $S_{ij} = \frac{1}{2} \left[ \frac{\partial}{\partial x_i} v_j + \frac{\partial}{\partial x_j} v_i \right]$ based on $a_{ij} \cdot b_{ij} = \sum_{i,j} a_{ij} b_{ij}$.

Dissipation balance:

$$\frac{\partial \left( \omega \right)}{\partial t} + \nabla \cdot \left( \mathbf{v} \omega \right) - \nabla \cdot \left[ \left( \nu + \sigma_{\omega} \nu_\textrm{t} \right) \nabla \omega \right] - \alpha \frac{\omega}{k} P + \beta_{\omega} \omega^2 - \frac{\sigma_d}{\omega} \nabla k \nabla \omega = 0$$

The kinematic eddy viscosity $\nu_\textrm{t}$ is calculated as follows:

$$\nu_\textrm{t} = \frac{k}{\tilde{\omega}}$$

With a limited dissipation:

$$\tilde{\omega} = \textrm{max} \left\{ \omega, 0.875 \sqrt{\frac{P}{\nu_\textrm{t} \beta_\textrm{k}}} \right\}$$

And a cross-diffusion coefficient $\sigma_\textrm{d}$

$$\sigma_\text{d} = \begin{cases} 0 & \mbox{, if } \; \nabla k \cdot \nabla \omega \le 0 \\ 0.125 & \mbox{, if } \; \nabla k \cdot \nabla \omega > 0 \end{cases}.$$

Classes
struct	Dumux::Properties::KOmegaModelTraits< dimension >
	Traits for the k-omega model. More...
struct	Dumux::Properties::TTag::KOmega
	The type tag for the single-phase, isothermal k-omega model. More...
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::KOmega >
	states some specifics of the isothermal k-omega model More...
struct	Dumux::Properties::FluxVariables< TypeTag, TTag::KOmega >
	The flux variables. More...
struct	Dumux::Properties::LocalResidual< TypeTag, TTag::KOmega >
	The local residual. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::KOmega >
	Set the volume variables property. More...
struct	Dumux::Properties::IOFields< TypeTag, TTag::KOmega >
	The specific I/O fields. More...
struct	Dumux::Properties::TTag::KOmegaNI
	The type tag for the single-phase, non-isothermal k-omega 2-Eq. model. More...
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::KOmegaNI >
	The model traits of the non-isothermal model. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::KOmegaNI >
	Set the volume variables property. More...
struct	Dumux::Properties::IOFields< TypeTag, TTag::KOmegaNI >
	The specific non-isothermal I/O fields. More...

Namespaces
namespace	Dumux
namespace	Dumux::Properties
namespace	Dumux::Properties::TTag
	Type tag for numeric models.

Include dependency graph for freeflow/rans/twoeq/komega/model.hh: