diff --git a/Solver/TESTCASES/ForwardFacingStep.lua b/Solver/TESTCASES/ForwardFacingStep.lua index d1bbe4aeab5fee353e59554e65febac1b6a93c44..c1bb142ad0f143e6d36fb0d91e044175f6054a99 100644 --- a/Solver/TESTCASES/ForwardFacingStep.lua +++ b/Solver/TESTCASES/ForwardFacingStep.lua @@ -25,7 +25,7 @@ end --[[ Example of a lua program driving the DG code --]] - +sv order = 1 print'*** Loading the mesh and the model ***' myModel = GModel () @@ -55,14 +55,12 @@ print'*** export ***' DG:exportSolution('output/solution_0') print'*** solve ***' +CFL = 2.0; -LC = 0.1*.1 -dt = .3*LC/(SOUND+V); -print('DT=',dt) - -for i=1,10 do +for i=1,1000 do + dt = CFL * DG:computeInvSpectralRadius(); norm = DG:RK44_limiter(0.1*dt) - print('*** ITER ***',i,norm) + print('*** ITER ***',i,dt, norm) if (i % 1 == 0) then DG:exportSolution(string.format("output/solution-%06d", i)) end