Skip to content
Snippets Groups Projects
Commit c5ebad2b authored by Jean-François Remacle's avatar Jean-François Remacle
Browse files

No commit message

No commit message
parent f1847cee
No related branches found
No related tags found
No related merge requests found
MACH = 0.1; MACH = 3;
RHO = 1.0; RHO = 1.0;
PRES = 1./(MACH*RHO*RHO*1.4*1.4) PRES = 1./(MACH*RHO*RHO*1.4*1.4)
V = 1.0 V = 1.0
...@@ -49,7 +49,7 @@ DG:exportSolution('output/solution_0') ...@@ -49,7 +49,7 @@ DG:exportSolution('output/solution_0')
print'*** solve ***' print'*** solve ***'
LC = 0.1 LC = 0.1*.1
dt = .3*LC/(SOUND+V); dt = .3*LC/(SOUND+V);
print('DT=',dt) print('DT=',dt)
......
Point(1) = {0, 0, 0, .1}; Point(1) = {0, 0, 0, .1};
Point(2) = {3, 0, 0, .1}; Point(2) = {.5, 0, 0, .1};
Point(3) = {3, 1, 0, .1}; Point(3) = {.5, .25, 0, .1};
Point(4) = {5, 1, 0, .1}; Point(4) = {3, .25, 0, .1};
Point(5) = {5, 2, 0, .1}; Point(5) = {3, 1, 0, .1};
Point(6) = {0, 2, 0, .1}; Point(6) = {0, 1, 0, .1};
Line(1) = {6, 5}; Line(1) = {6, 5};
Line(2) = {5, 4}; Line(2) = {5, 4};
Line(3) = {4, 3}; Line(3) = {4, 3};
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment