INSFVTKEDSourceSink

The object computes the turbulent source and sink term for the turbulent kinetic energy dissipation rate equation.

Two terms are computed: destruction and production and the term destruction - production is passed to the residual. A different treatment is used for the bulk and the near wall regions.

Bulk formulation:

The production of turbulent kinetic energy dissipation GϵG_\epsilon is modeled as follows:

Gϵ=C1,ϵϵkGk,G_{\epsilon} = C_{1,\epsilon} \frac{\epsilon}{k} G_k \,,

where:

  • C1,ϵ=1.44C_{1,\epsilon} = 1.44 is a closure parameter,

  • GkG_k is the limited turbulent kinetic energy production. For more details please refer to INSFVTKESourceSink.

The destruction of the turbulent kinetic energy dissipation rate is modeled as follows:

ϵϵ=C2,ϵρϵtk,\epsilon_{\epsilon} = \frac{C_{2,\epsilon} \rho \epsilon}{t_k} \,,

where:

  • C2,ϵ=1.92C_{2,\epsilon} = 1.92 is a closure parameter,

  • ϵ\epsilon is the solution variable, i.e., the dissipation rate of the turbulent kinetic energy,

  • kk is the turbulent kinetic energy,

  • tk=kϵt_k = \frac{k}{\epsilon} is the turbulent time scale; if the "linearized_model" is true, this timescale is computed from the previous iteration; if "linearized_model" is false, in a nonlinear solve, this timescale is aded to the Jacobian.

Wall formulation:

All cells in contact with a boundary identified in the "walls" list are applied a different treatment in which the equilibrium value for the ϵ=ϵeq\epsilon = \epsilon_{eq} is set. A separate formulation is used for the sub-laminar and logarithmic boundary layers. The determination of whether the near-wall cell lies in the laminar or logarithmic region is performed via the non-dimensional wall distance y+y^+. The non-dimensional wall distance can be defined differently according to the "wall_treatment" parameter.

The four formulations are described in more detail in INSFVTurbulentViscosityWallFunction.

If an equilibrium "wall_treatment" is defined, i.e. eq_newton,eq_incremental or eq_linearized, the standard wall function formulations are used in which y+y^+ is found:

y+=ρypuτμ,y^+ = \frac{\rho y_p u_{\tau}}{\mu} \,,

where:

  • ρ\rho is the density,

  • ypy_p is the distance from the wall to the centroid of the next-to-wall cell,

  • uτu_{\tau} is the friction velocity, defined as uτ=τwρu_{\tau} = \sqrt{\frac{\tau_w}{\rho}} with τw\tau_w the shear stress at the wall for which the condition is applied,

  • μ\mu is the dynamic molecular viscosity.

If a non-equilibrium "wall_treatment" is defined, i.e. neq, the y+y^+ is defined non-iteratively as follows:

y+=ypCμkμ,y^+ = \frac{y_p \sqrt{\sqrt{C_{\mu}}k}}{\mu} \,,

commentnote

Using non-equilibrium wall functions is recommended for problems with recirculations and boundary layer detachment. However, using non-equilibrium wall functions will deteriorate results for standard problems such as flow developing over walls.

The cells with y+<11.25y^+ < 11.25 belong to sub-laminar boundary layer. The ones belonging to the logarithmic boundary layer are those for which y+11.25y^+ \ge 11.25.

A different value is used for ϵeq\epsilon_{eq} in each of the two regions. For the sub-laminar boundary layer, the equilibrium value is determined as follows:

ϵeq=2kμyp2,\epsilon_{eq} = 2 \frac{k \mu}{y_p^2}\,,

where:

  • μt\mu_t is the turbulent dynamic viscosity.

For the logarithmic boundary layer, the value is determined as follows:

ϵeq=Cμ0.75ρk1.5κyp,\epsilon_{eq} = \frac{C_{\mu}^{0.75} \rho k^{1.5}}{\kappa y_p}\,,

where:

  • κ=0.4187\kappa = 0.4187 is the von Kármán constant.

commentnote

When using wall functions, since the equilibrium value for ϵ\epsilon is set in the cells near the wall, the user is recommended to deactivate advection and diffusion for those near wall cells.

commentnote

When the wall treatment is specified in this kernel, any boundary condition for ϵ\epsilon will be ignored. In other words, there is no need to impose boundary conditions for ϵ\epsilon when the wall treatment is specified in his kernel.

commentnote

When using near-wall treatment, we assume that the μt\mu_t functor is computed by an object that performs near-wall treatment. Otherwise, the results obtained won't be physically correct.

Input Parameters

  • muDynamic viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Dynamic viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • mu_tTurbulent viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Turbulent viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • rhofluid density. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:fluid density. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • tkeCoupled turbulent kinetic energy. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Coupled turbulent kinetic energy. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • uThe velocity in the x direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The velocity in the x direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • variableThe name of the variable that this residual object operates on

    C++ Type:NonlinearVariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:The name of the variable that this residual object operates on

Required Parameters

  • C1_eps1.44First epsilon coefficient

    Default:1.44

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:First epsilon coefficient

  • C2_eps1.92Second epsilon coefficient

    Default:1.92

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Second epsilon coefficient

  • C_mu0.09Coupled turbulent kinetic energy closure.

    Default:0.09

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coupled turbulent kinetic energy closure.

  • C_pl10Production limiter constant multiplier.

    Default:10

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Production limiter constant multiplier.

  • blockThe list of blocks (ids or names) that this object will be applied

    C++ Type:std::vector<SubdomainName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The list of blocks (ids or names) that this object will be applied

  • linearized_modelTrueBoolean to determine if the problem should be used in a linear or nonlinear solve

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Boolean to determine if the problem should be used in a linear or nonlinear solve

  • prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

  • use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

  • vThe velocity in the y direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The velocity in the y direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • wThe velocity in the z direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The velocity in the z direction. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • wall_treatmentneqThe method used for computing the wall functions 'eq_newton', 'eq_incremental', 'eq_linearized', 'neq'

    Default:neq

    C++ Type:MooseEnum

    Unit:(no unit assumed)

    Options:eq_newton, eq_incremental, eq_linearized, neq

    Controllable:No

    Description:The method used for computing the wall functions 'eq_newton', 'eq_incremental', 'eq_linearized', 'neq'

  • wallsBoundaries that correspond to solid walls.

    C++ Type:std::vector<BoundaryName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Boundaries that correspond to solid walls.

Optional Parameters

  • absolute_value_vector_tagsThe tags for the vectors this residual object should fill with the absolute value of the residual contribution

    C++ Type:std::vector<TagName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The tags for the vectors this residual object should fill with the absolute value of the residual contribution

  • extra_matrix_tagsThe extra tags for the matrices this Kernel should fill

    C++ Type:std::vector<TagName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The extra tags for the matrices this Kernel should fill

  • extra_vector_tagsThe extra tags for the vectors this Kernel should fill

    C++ Type:std::vector<TagName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The extra tags for the vectors this Kernel should fill

  • matrix_tagssystemThe tag for the matrices this Kernel should fill

    Default:system

    C++ Type:MultiMooseEnum

    Unit:(no unit assumed)

    Options:nontime, system

    Controllable:No

    Description:The tag for the matrices this Kernel should fill

  • vector_tagsnontimeThe tag for the vectors this Kernel should fill

    Default:nontime

    C++ Type:MultiMooseEnum

    Unit:(no unit assumed)

    Options:nontime, time

    Controllable:No

    Description:The tag for the vectors this Kernel should fill

Tagging Parameters

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:Yes

    Description:Set the enabled status of the MooseObject.

  • implicitTrueDetermines whether this object is calculated using an implicit or explicit form

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Determines whether this object is calculated using an implicit or explicit form

  • newton_solveFalseWhether a Newton nonlinear solve is being used

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Whether a Newton nonlinear solve is being used

  • seed0The seed for the master random number generator

    Default:0

    C++ Type:unsigned int

    Unit:(no unit assumed)

    Controllable:No

    Description:The seed for the master random number generator

  • use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

  • ghost_layers2The number of layers of elements to ghost.

    Default:2

    C++ Type:unsigned short

    Unit:(no unit assumed)

    Controllable:No

    Description:The number of layers of elements to ghost.

  • use_point_neighborsFalseWhether to use point neighbors, which introduces additional ghosting to that used for simple face neighbors.

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Whether to use point neighbors, which introduces additional ghosting to that used for simple face neighbors.

Parallel Ghosting Parameters