Skip to content
Snippets Groups Projects
Commit a1deba0e authored by Anna Wellmann's avatar Anna Wellmann
Browse files

Add omega modification to sponge layer kernel

parent 63f02aca
No related branches found
No related tags found
1 merge request!186Fix boundary condition edge diagonals
......@@ -342,6 +342,24 @@ __global__ void LB_Kernel_CumulantK17Sponge(
//!
real omega = omega_in;
if(turbulenceModel != TurbulenceModel::None){ omega /= (c1o1 + c3o1*omega_in*turbulentViscosity[k_000]); }
////////////////////////////////////////////////////////////////////////////////////
// Calculate modified omega for sponge layer
// sponge layer inflow
real startXsponge = 1507.0f;
real endXsponge = 1537.0f;
real sizeSponge = endXsponge - startXsponge;
if (coordX[k_000] > startXsponge) {
real spongeFactor = (((endXsponge - coordX[k_000]) / sizeSponge) * c1o2) + c1o2;
omega = spongeFactor * omega;
}
//sponge layer outflow
endXsponge = 30.0f;
if (coordX[k_000] < endXsponge) {
real spongeFactor = (((coordX[k_000]) / endXsponge) * c1o2) + c1o2;
omega = spongeFactor * omega;
}
////////////////////////////////////////////////////////////
// 2.
real OxxPyyPzz = c1o1;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment