blast_3d_amr.in

# AthenaPK - a performance portable block structured AMR MHD code
# Copyright (c) 2020, Athena Parthenon Collaboration. All rights reserved.
# Licensed under the BSD 3-Clause License (the "LICENSE");

<comment>
problem = spherical blast wave

<job>
problem_id = blast

<parthenon/mesh>
refinement = adaptive
numlevel = 3
nghost = 2

nx1 = 32
x1min = -0.5
x1max = 0.5
ix1_bc = periodic
ox1_bc = periodic

nx2 = 32
x2min = -0.5
x2max = 0.5
ix2_bc = periodic
ox2_bc = periodic

nx3 = 32
x3min = -0.5
x3max = 0.5
ix3_bc = periodic
ox3_bc = periodic

<parthenon/meshblock>
nx1=8 
nx2=8 
nx3=8 

<parthenon/time>
integrator = vl2
cfl = 0.3
tlim = 0.10
nlim = 100000
perf_cycle_offset = 2 # number of inital cycles not to be included in perf calc
ncycle_out_mesh = -100

<hydro>
eos = adiabatic
riemann = hlle
reconstruction = plm
gamma = 1.666666666666667 # gamma = C_p/C_v
scratch_level = 0 # 0 is actual scratch (tiny); 1 is HBM

<refinement>
type                         = pressure_gradient
threshold_pressure_gradient  = 0.1

<problem/blast>
pressure_ambient  = 0.001   # ambient pressure
pressure_ratio    = 1.6e8   # Pressure ratio initially
radius_outer      = 0.03125 # Radius of the outer sphere
radius_inner      = 0.0     # Radius of the inter sphere (with ramp between inner and outer sphere)
density_ambient   = 1.0     # ambient density
density_ratio     = 1.0     # density ratio in sphere

<parthenon/output0>
file_type = hdf5
dt = 0.01
variables = cons