problem = Brio & Wu shock tube author = M. Brio & C. C. Wu journal = J. Comp. Phys. 75, 400-422 (1988) config = --with-problem=shkset1d problem_id = Brio-Wu # problem ID: basename of output filenames maxout = 3 # Output blocks number from 1 -> maxout out_fmt = tab # Tabular data dump dt = 0.0025 # time increment between outputs out_fmt = hst # History data dump dt = 0.0025 # time increment between outputs out_fmt = bin # Binary data dump dt = 0.0025 # time increment between outputs cour_no = 0.8 # The Courant, Friedrichs, & Lewy (CFL) Number nlim = 10000 # cycle limit tlim = 1. # time limit Nx1 = 800 # Number of zones in X1-direction x1min = 0.0 # minimum value of X1 x1max = 1.0 # maximum value of X1 ibc_x1 = 2 # inner (X1) boundary flag obc_x1 = 2 # outer (X1) boundary flag Nx2 = 1 # Number of zones in X2-direction x2min = 0.0 # minimum value of X2 x2max = 1.0 # maximum value of X2 ibc_x2 = 2 # inner (X2) boundary flag obc_x2 = 2 # outer (X2) boundary flag Nx3 = 1 # Number of zones in X3-direction x3min = 0.0 # minimum value of X3 x3max = 1.0 # maximum value of X3 ibc_x3 = 2 # inner (X3) boundary flag obc_x3 = 2 # outer (X3) boundary flag NGrid_x1 = 1 NGrid_x2 = 1 NGrid_x3 = 1 gamma = 2.0 # gamma = C_p/C_v shk_dir = 1 # Shock Direction -- (1,2,3) = (x1,x2,x3) dl = 1.0 # density on left half of grid pl = 1.0 # pressure v1l = 0.0 # X-velocity v2l = 0.0 # Y-velocity v3l = 0.0 # Z-velocity b1l = 0.75 # X-magnetic field b2l = 1.0 # Y-magnetic field b3l = 0.0 # Z-magnetic field s0l = 1.0 # "color" dr = 0.125 # density on right half of grid pr = 0.1 # pressure v1r = 0.0 # X-velocity v2r = 0.0 # Y-velocity v3r = 0.0 # Z-velocity b1r = 0.75 # X-magnetic field b2r = -1.0 # Y-magnetic field b3r = 0.0 # Z-magnetic field s0r = 0.0 # "color"