From 3d4480564e11c093deac217041eddc7da3fccc9c Mon Sep 17 00:00:00 2001
From: adfboyd <adfboyd@UbuntuVM.myguest.virtualbox.org>
Date: Sun, 13 Oct 2024 22:11:36 +0100
Subject: [PATCH] added inviscid option if Re<0.0001
---
src/parameters.f90 | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/parameters.f90 b/src/parameters.f90
index 6dab4241..8c9bcc49 100644
--- a/src/parameters.f90
+++ b/src/parameters.f90
@@ -255,7 +255,11 @@ subroutine parameter(input_i3d)
dy2 = dy * dy
dz2 = dz * dz
- xnu=one/re
+ if (re.gt.0.001) then
+ xnu=one/re
+ else
+ xnu=zero
+ endif
!! Constant pressure gradient, re = Re_tau -> use to compute Re_centerline
if (cpg) then
re_cent = (re/0.116_mytype)**(1.0_mytype/0.88_mytype)
--
GitLab