|
flow 0.4.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
|
Functions | |
| lattice_sdf (N, phi=0.125) | |
| drag (N, mode, mu=0.1, F=1e-3, dt=80.0, warm_tol=1e-7, tail=40, max_steps=4000) | |
Variables | |
| float | kref = 4.2920 |
| dict | M0 = {32: +1.004, 48: +0.685, 64: +0.598, 96: +0.397, 128: +0.299} |
| dict | MG = {32: -0.175, 48: -0.084, 64: -0.056, 96: -0.029, 128: -0.018} |
| flush | |
| dict | prev = {} |
| t0 = time.time() | |
| K9 | |
| i9 | |
| K10 | |
| i10 | |
| int | e9 = 100 * (K9 - kref) / kref |
| int | e10 = 100 * (K10 - kref) / kref |
| dict | o9 = np.log(abs(prev["9"]) / abs(e9)) / np.log(N / prev["N"]) if prev else float("nan") |
| dict | o10 = np.log(abs(prev["10"]) / abs(e10)) / np.log(N / prev["N"]) if prev else float("nan") |
Z&H drag for the collocated cutcell-ghost HYBRID modes (set_face_interp 9/10): mode-0's aperture projection with the directional gpCenterGrad -grad(P)/cell correction (9), plus the open-centroid wall-aware constraint quadrature (10). Question: how much of the mode-0 first-order drag was the O(1/h) gradient defect alone (9), and does the a-priori-O(h^2) open-centroid flux quadrature finally pay off once paired with a telescoping 2nd-order force (10)? Baselines: mode-0 +1.00/+0.69/+0.60/+0.40/+0.30 %, ghost (1,2) -0.175/-0.084/-0.056/-0.029/-0.018 % at N=32..128.
| collocated_zh_hybrid.lattice_sdf | ( | N, | |
phi = 0.125 |
|||
| ) |
Definition at line 17 of file collocated_zh_hybrid.py.
Referenced by drag().
| collocated_zh_hybrid.drag | ( | N, | |
| mode, | |||
mu = 0.1, |
|||
F = 1e-3, |
|||
dt = 80.0, |
|||
warm_tol = 1e-7, |
|||
tail = 40, |
|||
max_steps = 4000 |
|||
| ) |
Definition at line 30 of file collocated_zh_hybrid.py.
References lattice_sdf().
| float collocated_zh_hybrid.kref = 4.2920 |
Definition at line 61 of file collocated_zh_hybrid.py.
| dict collocated_zh_hybrid.M0 = {32: +1.004, 48: +0.685, 64: +0.598, 96: +0.397, 128: +0.299} |
Definition at line 62 of file collocated_zh_hybrid.py.
| dict collocated_zh_hybrid.MG = {32: -0.175, 48: -0.084, 64: -0.056, 96: -0.029, 128: -0.018} |
Definition at line 63 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.flush |
Definition at line 65 of file collocated_zh_hybrid.py.
| dict collocated_zh_hybrid.prev = {} |
Definition at line 68 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.t0 = time.time() |
Definition at line 70 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.K9 |
Definition at line 71 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.i9 |
Definition at line 71 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.K10 |
Definition at line 72 of file collocated_zh_hybrid.py.
| collocated_zh_hybrid.i10 |
Definition at line 72 of file collocated_zh_hybrid.py.
Definition at line 73 of file collocated_zh_hybrid.py.
Definition at line 74 of file collocated_zh_hybrid.py.
| dict collocated_zh_hybrid.o9 = np.log(abs(prev["9"]) / abs(e9)) / np.log(N / prev["N"]) if prev else float("nan") |
Definition at line 75 of file collocated_zh_hybrid.py.