Skip to content
Snippets Groups Projects
Commit f1847cee authored by Bruno Seny's avatar Bruno Seny
Browse files

BackwardFacingStep working in 2d for Euler equations thanks to limiter

parent bff6847f
No related branches found
No related tags found
No related merge requests found
...@@ -41,10 +41,11 @@ DG:setup() ...@@ -41,10 +41,11 @@ DG:setup()
print'*** setting the initial solution ***' print'*** setting the initial solution ***'
DG:L2Projection(FS) DG:L2Projection(FS)
DG:limitSolution()
print'*** export ***' print'*** export ***'
DG:exportSolution('solution_0') DG:exportSolution('output/solution_0')
print'*** solve ***' print'*** solve ***'
...@@ -52,11 +53,11 @@ LC = 0.1 ...@@ -52,11 +53,11 @@ LC = 0.1
dt = .3*LC/(SOUND+V); dt = .3*LC/(SOUND+V);
print('DT=',dt) print('DT=',dt)
for i=1,1000 do for i=1,10000 do
norm = DG:RK44(dt) norm = DG:RK44_limiter(dt)
print('*** ITER ***',i,norm) print('*** ITER ***',i,norm)
if (i % 10 == 0) then if (i % 100 == 0) then
DG:exportSolution(string.format("solution-%03d", i)) DG:exportSolution(string.format("output/solution-%06d", i))
end end
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment