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

Go to the source code of this file.

Namespaces

namespace  three_solver_study
 

Functions

 three_solver_study._grid (N)
 
 three_solver_study._mm (d, N)
 
 three_solver_study.sdf_zh (N, phi=0.216)
 
 three_solver_study.sdf_random (N, n=8, r_frac=0.18, jit=0.06, seed=12345)
 
 three_solver_study.run (case, N, solver="staggered", re=0.0, mu=0.1, F=1e-3, dt=60.0, max_steps=500, tol=1e-6)
 

Variables

list three_solver_study.ZH_PHI = [0.000125,0.001,0.008,0.027,0.064,0.125,0.216,0.343,0.45,0.5236]
 
list three_solver_study.ZH_K = [1.096,1.212,1.525,2.008,2.810,4.292,7.442,15.4,28.1,42.1]
 
float three_solver_study.zh_ref = lambda phi(np.interp(phi, ZH_PHI, ZH_K))
 
list three_solver_study.cases = [("zh",[16,24,32,48]), ("random",[24,32,48])]
 
 three_solver_study.flush
 
float three_solver_study.ref = "zh" else None
 
 three_solver_study.r = run(case,N,solver,re=0.0); gc.collect()
 
str three_solver_study.err = f"{100*(r['val']-ref)/ref:+.2f}%" if ref else ""