Commit 0d0354a7 authored by Jan Schäfer's avatar Jan Schäfer Committed by Jakub Klinkovský
Browse files

bracket in Prandtl model for 3D turbulent viscosity

parent e60319fe

src/Examples/turbulent-flows/zero-equation-model/Getter.h

+24 −22
Original line number Diff line number Diff line
@@ -253,6 +253,7 @@ class TurbulentViscosityGetter< Meshes::Grid< 3, MeshReal, Device, MeshIndex >, 
           const RealType& velocity_z_center    = this->velocity.template getData< TNL::Devices::Host >()[ 2 ].template getData< DeviceType >()[ center ]; 



           return rho * this->mixingLength * this->mixingLength *

                  ( 

                     (  ( velocity_y_east - velocity_y_west ) * hxInverse * 0.5

                     -  ( velocity_x_north - velocity_x_south ) * hyInverse * 0.5

                     )
@@ -275,6 +276,7 @@ class TurbulentViscosityGetter< Meshes::Grid< 3, MeshReal, Device, MeshIndex >, 
                     *

                     (  ( velocity_y_up - velocity_y_down ) * hzInverse * 0.5

                     -  ( velocity_z_north - velocity_z_south ) *hyInverse * 0.5

                     )

                  );

           }

           else return 0;