From da942a863ed9eca6572a1f8fddf51de607ed43c3 Mon Sep 17 00:00:00 2001
From: Bruno Seny <bruno.seny@student.uclouvain.be>
Date: Fri, 29 Jan 2010 16:16:56 +0000
Subject: [PATCH] Point(1) = {-0.5, -0.5, 0, .03}; Point(2) = {0, -0.5, 0,
 .03}; Point(3) = {0, 0.5, 0, .03}; Point(4) = {-0.5, 0.5, 0, .03};

Point(5) = {0.5, -0.5, 0, .03};
Point(6) = {0.5, 0.5, 0, .03};

Line(1) = {4, 3};
Line(2) = {3, 2};
Line(3) = {2, 1};
Line(4) = {1, 4};

Line(7) = {3, 6};
Line(8) = {6, 5};
Line(9) = {5, 2};

Line Loop(5) = {2, 3, 4, 1};
Line Loop(6) = {8, 9, -2, 7};
Plane Surface(7) = {5};
Plane Surface(8) = {6};
Transfinite Line {4,8} = 20 ;
Transfinite Line{2}=20;
Transfinite Surface {7};
Recombine Surface {7};
Transfinite Line {3,1,7,9} = 10;
Transfinite Surface{8};

Physical Line("Border") = {1, 3, 4, 7, 8, 9};
Physical Surface("Inside1") = {7};
Physical Surface("Inside2") = {8};
---
 Solver/TESTCASES/CylinderEddies.lua | 6 +++---
 Solver/TESTCASES/cyl.geo            | 1 +
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/Solver/TESTCASES/CylinderEddies.lua b/Solver/TESTCASES/CylinderEddies.lua
index 0a7d9cf3c5..92987c61d0 100644
--- a/Solver/TESTCASES/CylinderEddies.lua
+++ b/Solver/TESTCASES/CylinderEddies.lua
@@ -1,7 +1,7 @@
 MACH = .05;
 RHO  = 1.0;
 PRES = 1./(MACH*RHO*RHO*1.4*1.4) 
-V = 1.0 
+V = 1 
 SOUND = V/MACH
 
 --[[ 
@@ -26,7 +26,7 @@ FS = functionLua(4, 'free_stream', {'XYZ'}):getName()
 
 -- diffusivity
 mu=fullMatrix(1,1);
-mu:set(0,0,0.01)
+mu:set(0,0,0.025)
 kappa=fullMatrix(1,1);
 kappa:set(0,0,0.01)
 
@@ -62,7 +62,7 @@ CFL = 20.1;
 dt = CFL * DG:computeInvSpectralRadius();
 print('DT = ',dt)
 T = 0;
-for i=1,1000000 do
+for i=1,2 do
     dt = CFL * DG:computeInvSpectralRadius();    
     norm = DG:RK44(dt)
     T = T + dt
diff --git a/Solver/TESTCASES/cyl.geo b/Solver/TESTCASES/cyl.geo
index 5779ec0566..8b8595b31e 100644
--- a/Solver/TESTCASES/cyl.geo
+++ b/Solver/TESTCASES/cyl.geo
@@ -73,3 +73,4 @@ Physical Surface("BoundaryLayer") = {100};
 Physical Line("Cylinder") = {4, 1, 2, 3};
 Physical Line("Box") = {5, 6, 7, 8};
 Mesh.CharacteristicLengthExtendFromBoundary=1;
+Recombine Surface {4, 1, 2, 3};
-- 
GitLab