flow 0.4.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
Loading...
Searching...
No Matches
flatwall_displacement.py File Reference

Go to the source code of this file.

Namespaces

namespace  flatwall_displacement
 

Functions

 flatwall_displacement.channel (N, s)
 
 flatwall_displacement.exact_profile (exp, y, w_lo, w_hi)
 
 flatwall_displacement.exact_discharge (exp, w_lo, w_hi)
 
 flatwall_displacement.make_solver (kind, N, s)
 
 flatwall_displacement.set_forces (sol, exp, ny, w_lo, w_hi)
 
 flatwall_displacement.march (sol)
 
 flatwall_displacement.run (kind, exp, N, s)
 

Variables

 flatwall_displacement.MU
 
 flatwall_displacement.F0
 
int flatwall_displacement.NX = 8
 
int flatwall_displacement.SLAB = 8
 
 flatwall_displacement.TOL = float(os.environ.get("MARCH_TOL", "1e-9"))
 
 flatwall_displacement.CAP = int(os.environ.get("MARCH_MAX", "8000"))
 
 flatwall_displacement.DT = float(os.environ.get("DT", "60.0"))
 
 flatwall_displacement.ap = argparse.ArgumentParser()
 
 flatwall_displacement.default
 
 flatwall_displacement.a = ap.parse_args()
 
list flatwall_displacement.Ns = [int(x) for x in a.N.split(",")]
 
list flatwall_displacement.Ss = [float(x) for x in a.s.split(",")]
 
 flatwall_displacement.exps = a.exp.split(",")
 
 flatwall_displacement.kinds = a.solvers.split(",")
 
str flatwall_displacement.hdr = f"{'N':>4} {'s':>5} {'solver':>12} {'steps':>6} {'e_flux':>11} {'e_cell':>11} {'L2prof':>10}"
 
 flatwall_displacement.flush
 
dict flatwall_displacement.base = {}
 
 flatwall_displacement.r = run(kind, exp, N, s)
 
tuple flatwall_displacement.row
 
dict flatwall_displacement.r1 = base.get(kind)
 
 flatwall_displacement.du = float(np.abs(r['U'] - r1['U']).max()) / r['uscale']
 
 flatwall_displacement.dv = float(np.abs(r['V']).max()) / r['uscale']