EpsilonBoundaryCondition
Posted: Mon Mar 23, 2020 11:35 am
Hi,
I am trying to impose an epsilon boundary condition to the near wall node I that is equal to the one proposed by Chien and Launder (1980) paper. The approach is to tweak the discretised equation of the near wall node by using a big source term and an imposed value equivalent to the equilibrium lenght scale if y+>5.1. The implementation seems to be ok but I am testing the approach on a rectangular cavity and it seems that something goes wrong ( code attached). A few questions:
1) Does the imposed value need to be assigned to the pointer cvar_ep(iel) or to a new dynamic array ?
2) Is there any examples on how to do a weighting when there are 2 boundary faces for a single cell volume?
3) Is epsilon equation solved before k when ikecou=0?
Best regards,
Costas
I am trying to impose an epsilon boundary condition to the near wall node I that is equal to the one proposed by Chien and Launder (1980) paper. The approach is to tweak the discretised equation of the near wall node by using a big source term and an imposed value equivalent to the equilibrium lenght scale if y+>5.1. The implementation seems to be ok but I am testing the approach on a rectangular cavity and it seems that something goes wrong ( code attached). A few questions:
1) Does the imposed value need to be assigned to the pointer cvar_ep(iel) or to a new dynamic array ?
2) Is there any examples on how to do a weighting when there are 2 boundary faces for a single cell volume?
3) Is epsilon equation solved before k when ikecou=0?
Best regards,
Costas