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

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")
 

Detailed Description

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.

Function Documentation

◆ lattice_sdf()

collocated_zh_hybrid.lattice_sdf (   N,
  phi = 0.125 
)

Definition at line 17 of file collocated_zh_hybrid.py.

Referenced by drag().

◆ 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().

Variable Documentation

◆ kref

float collocated_zh_hybrid.kref = 4.2920

Definition at line 61 of file collocated_zh_hybrid.py.

◆ M0

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.

◆ MG

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.

◆ flush

collocated_zh_hybrid.flush

Definition at line 65 of file collocated_zh_hybrid.py.

◆ prev

dict collocated_zh_hybrid.prev = {}

Definition at line 68 of file collocated_zh_hybrid.py.

◆ t0

collocated_zh_hybrid.t0 = time.time()

Definition at line 70 of file collocated_zh_hybrid.py.

◆ K9

collocated_zh_hybrid.K9

Definition at line 71 of file collocated_zh_hybrid.py.

◆ i9

collocated_zh_hybrid.i9

Definition at line 71 of file collocated_zh_hybrid.py.

◆ K10

collocated_zh_hybrid.K10

Definition at line 72 of file collocated_zh_hybrid.py.

◆ i10

collocated_zh_hybrid.i10

Definition at line 72 of file collocated_zh_hybrid.py.

◆ e9

int collocated_zh_hybrid.e9 = 100 * (K9 - kref) / kref

Definition at line 73 of file collocated_zh_hybrid.py.

◆ e10

int collocated_zh_hybrid.e10 = 100 * (K10 - kref) / kref

Definition at line 74 of file collocated_zh_hybrid.py.

◆ o9

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.

◆ o10

dict collocated_zh_hybrid.o10 = np.log(abs(prev["10"]) / abs(e10)) / np.log(N / prev["N"]) if prev else float("nan")

Definition at line 76 of file collocated_zh_hybrid.py.