1 scalar1 = 10; 400 0 8.55e-05 1 SGDH

1 all[] all[] -0.25 2.25 0 0.05 2.25 0 0.05 2.75 0 0.05 0.5 0 0.05 -0.25 0 0.75 -0.25 0 0.75 0.25 0 0.75 0.75 0 -0.5 2.25 0 1 2

300 0 0.05 1 34 8 or 9 or 28 or 29 or 38 or 39 24 and 0.1<=x and 0.5>=x 2 or 3 4 or 7 or 21 or 22 or 23 6 and y>1 6 and y<=1 31 or 33 200 if (0<=t<=3.8){ temperature=20 + 100*t; } else { temperature=400; } 0.5 1 0 0 1 0 1 1 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 25 0 0 0 0 320 0 40 0 1 1 0 A=-4.0668*10^(-3); B=-5.0754*10^(-2); C=1000.9; density=temperature*(A*temperature+B)+C; 725.737 8.951e-05 5483 0.02495 101325 20 0 -9.81 0 0 0 0 0.1 25 5,24,32 1 1 all[] mesh_viscosity_1 = 1; 0 1 temperature = 20.; 400 0 1 SGDH 1 0 1 0 1 3 velocity[0] = 0.; velocity[1] = 0.; velocity[2] = 0.; 2 1 1