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