config

[backend]
precision = single
rank-allocator = linear

[backend-cuda]
device-id = local-rank

[constants]
nu = 0.005
Uin = 1.0
Vin = 0.0
Win = 0.0
Pc = 1.0
ac-zeta = 2.5

[solver]
system = ac-navier-stokes
order = 3

[solver-time-integrator]
formulation = dual
scheme = sdirk33
pseudo-scheme = rk45
controller = none
pseudo-controller = local-pi
tstart = 0.0
tend = 375.0
dt = 1.25e-2
pseudo-dt = 1.0e-3
pseudo-niters-min = 5
pseudo-niters-max = 5
pseudo-resid-norm = l2
pseudo-resid-tol = 1e-3
atol = 1e-3
pseudo-dt-max-mult = 2.5

[solver-dual-time-integrator-multip]
pseudo-dt-fact = 1.75
cycle = [(3, 1), (2, 1), (1, 1), (0, 2), (1, 1), (2, 1), (3, 1)]

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

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

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

[solver-interfaces-tri]
flux-pts = williams-shunn

[solver-elements-tri]
soln-pts = williams-shunn

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

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

[solver-elements-pri]
soln-pts = williams-shunn~gauss-legendre

[soln-plugin-nancheck]
nsteps = 50

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

[soln-bcs-in]
type = ac-in-fv
u = Uin
v = Vin
w = Win

[soln-bcs-out]
type = ac-out-fp
p = Pc

[soln-bcs-sym1]
type = slp-wall

[soln-bcs-sym2]
type = slp-wall

[soln-bcs-front]
type = slp-wall

[soln-ics]
u = 1.0 + 0.00675*cos(2*3.1415*x/20.0)*sin(2*3.1415*y/20.0)
v = 0.00675*sin(2*3.1415*x/20.0)*cos(2*3.1415*y/20.0)
w = 0.0
p = Pc

#[solver-source-terms]
#u = 0.015*exp(-2*t)*sin(x)*sin(y)
#v = 0.015*exp(-2*t)*sin(x)*sin(y)

[soln-plugin-writer]
dt-out = 0.075
basedir = results/
basename = soln-{t:.2f}
post-action = echo "Wrote file {soln} at time {t} for mesh {mesh}."
post-action-mode = blocking

[soln-plugin-pseudostats]
flushsteps = 500
file = results/residual.csv
header = true

[soln-plugin-tavg]
nsteps = 10
dt-out = 7.5
mode = continuous
basedir = results/
tstart = 22.5
basename = averages-{t:06.2f}

avg-u = u
avg-v = v
avg-uu = u*u
avg-vv = v*v
avg-uv = u*v

fun-avg-upup = uu - u*u
fun-avg-vpvp = vv - v*v
fun-avg-upvp = uv - u*v
fun-avg-urms = sqrt(uu - u*u + vv - v*v)