From 41a197f17f2ad6dbb5144d86d54edd2399338475 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Gizem=20=C3=96zler?= <gizemozler95@gmail.com>
Date: Mon, 28 Aug 2023 09:24:16 +0000
Subject: [PATCH] Update particlesTransport.f90
---
src/particlesTransport.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/particlesTransport.f90 b/src/particlesTransport.f90
index b8428d4..6cdcfe6 100755
--- a/src/particlesTransport.f90
+++ b/src/particlesTransport.f90
@@ -81,7 +81,7 @@
do 1 n=1,np
urel= sqrt((uf(n)-up(n))**2+(vf(n)-vp(n))**2+(wf(n)-wp(n))**2)
Reyp= max(1.e-10_pr,2._pr*radp(n)*urel/nuf)
- Reyf= radp(n)**2/nuf*sqrt(dufdy(n)**2+dufdz(n)**2)
+ Reyf= 4._pr*radp(n)**2/nuf*sqrt(dufdy(n)**2+dufdz(n)**2)
beta= Reyf/Reyp/2._pr
if (Reyp.le.40._pr) then
--
GitLab