Square Cylinder M1.5 blows up

General
1

Dear All-
I was trying to simulate the flow past over a square cylinder for Mach 1.5. It starts but afte a few time steps it just blows up. I have attached my .ini file. I will appreciate if you can provide any helpful pointer. sq_geo.ini

Thanks!

[constants]
gamma = 1.4
mu = 0.0001
Pr = 0.71
M = 1.5


[solver]
system = navier-stokes
order = 4


[solver-time-integrator]
scheme = rk4
controller = none
tstart = 0.0
tend = 200
dt = 1e-5


[solver-interfaces]
riemann-solver = hllc
ldg-beta = 0.5
ldg-tau = 0.1

[solver-interfaces-line]
flux-pts = gauss-legendre

[solver-elements-quad]
soln-pts = gauss-legendre

[soln-plugin-sampler]
nsteps = 10
samp-pts = [(1.5, 0.0), (7.0, 0.0), (12.0, 0.0)]
format = primitive
file = point-data.csv
header = true

[soln-plugin-nancheck]
nsteps = 50

[soln-plugin-fluidforce-wall]
nsteps = 10
quad-deg = 4
file = square-forces.csv
header = true

[soln-plugin-writer]
dt-out = 0.5
basedir = .
basename = sq-geo-{t:.2f}

[soln-bcs-inflow]
type = char-riem-inv
rho = 1.0
u = 1.0
v = 0
w = 0
p = 0.31746031746031744

[soln-bcs-outflow]
type = char-riem-inv
rho = 1.0
u = 1.0
v = 0
w = 0
p = 0.31746031746031744

[soln-bcs-slip]
type = char-riem-inv
rho = 1.0
u = 1.0
v = 0
w = 0
p = 0.3174603174603175

[soln-bcs-wall]
type = no-slp-adia-wall

[soln-ics]
rho = 1
u = 1
v = 0
w = 0
p = 0.31746031746031744
2

Try starting with a lower order e.g. order = 1, and try reducing the dt. Also, given the flow is supersonic you will need to turn on one of the shock capturing approaches (e.g. entropy filtering [solver-entropy-filter] — Documentation).

3

Thanks @p.vincent I did as you suggested but it still blows up.Thanks!