From b73b4a4812ae4cac64fe9a5b77265a42078de4dd Mon Sep 17 00:00:00 2001 From: lumlauf Date: Tue, 28 May 2024 13:57:54 +0200 Subject: [PATCH] bug fix in k-omega model --- src/turbulence/turbulence.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/turbulence/turbulence.F90 b/src/turbulence/turbulence.F90 index a4bff3fc..41346c10 100644 --- a/src/turbulence/turbulence.F90 +++ b/src/turbulence/turbulence.F90 @@ -2058,7 +2058,7 @@ subroutine analyse_model ! compute model propeties sig_k = sig_kw gen_d = - 1./cw2 - gen_alpha = -4.*sqrt(sig_k)/(3.*sqrt(sig_k)-sqrt(sig_k+24.*sig_w*cw2)) + gen_alpha = 4.*sqrt(sig_k)/(3.*sqrt(sig_k)-sqrt(sig_k+24.*sig_w*cw2)) gen_l = cm0*sqrt(rcm)*sqrt( (5.*sig_k+12.*sig_w*cw2 - & 3.*sqrt(sig_k*(sig_k+24.*sig_w*cw2 ) ) ) / 12. ) @@ -3558,7 +3558,7 @@ REALTYPE function omega_bc(bc,type,zi,ki,z0,u_tau) if (bc.eq.Dirichlet) then omega_bc = K**0.5/(cm0*gen_l)*(zi+z0)**(0.5*gen_alpha-1.) else - omega_bc = -cmsf*K*(0.5*alpha-1.)/(sig_w*cm0) & + omega_bc = -cmsf*K*(0.5*gen_alpha-1.)/(sig_w*cm0) & *(zi+z0)**(gen_alpha-1.) endif case default @@ -3698,7 +3698,7 @@ REALTYPE function psi_bc(bc,type,zi,ki,z0,u_tau) psi_bc = - (gen_m*gen_alpha+gen_n)*cmsf*cm0**gen_p/sig_psi & *K**(gen_m+0.5)*gen_l**(gen_n+1.) & *(zi+z0)**((gen_m+0.5)*gen_alpha+gen_n) - endif + endif case default end select end function psi_bc