From f1847ceef53ee827aba985b93a90040273d71c31 Mon Sep 17 00:00:00 2001
From: Bruno Seny <bruno.seny@student.uclouvain.be>
Date: Tue, 22 Dec 2009 13:25:18 +0000
Subject: [PATCH] BackwardFacingStep working in 2d for Euler equations thanks
to limiter
---
Solver/TESTCASES/BackwardFacingStep.lua | 11 ++++++-----
1 file changed, 6 insertions(+), 5 deletions(-)
diff --git a/Solver/TESTCASES/BackwardFacingStep.lua b/Solver/TESTCASES/BackwardFacingStep.lua
index 7002a931fe..a6b545b272 100644
--- a/Solver/TESTCASES/BackwardFacingStep.lua
+++ b/Solver/TESTCASES/BackwardFacingStep.lua
@@ -41,10 +41,11 @@ DG:setup()
print'*** setting the initial solution ***'
DG:L2Projection(FS)
+DG:limitSolution()
print'*** export ***'
-DG:exportSolution('solution_0')
+DG:exportSolution('output/solution_0')
print'*** solve ***'
@@ -52,11 +53,11 @@ LC = 0.1
dt = .3*LC/(SOUND+V);
print('DT=',dt)
-for i=1,1000 do
- norm = DG:RK44(dt)
+for i=1,10000 do
+ norm = DG:RK44_limiter(dt)
print('*** ITER ***',i,norm)
- if (i % 10 == 0) then
- DG:exportSolution(string.format("solution-%03d", i))
+ if (i % 100 == 0) then
+ DG:exportSolution(string.format("output/solution-%06d", i))
end
end
--
GitLab