|
flow 0.4.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
|
Functions | |
| channel_sdf (nx, ny, nz, ylo, yhi) | |
| run (N, ghost, rho=1.0, mu=0.1, dt=50.0, F=0.01, max_steps=400) | |
| main () | |
Verification (flow, ghost projection): plane Poiseuille through an SDF channel with the EXPERIMENTAL directional ghost-cell projection active (set_ghost_projection + cutcell_pressure). The steady parabola u(y) = F/(2 mu)(y-ylo)(yhi-y) is exactly quadratic, so a second-order scheme must reproduce it pointwise to solver tolerance. For this unidirectional flow the projection RHS is identically zero (u* depends on y only, v = w = 0, and every ghost closure sees zero data), so this gates that the ghost-projection machinery — overlay build, binary surrogate, BiCGStab — does NOT corrupt an exact solution (classification/indexing bugs inject spurious divergence and break the parabola). The closure-sign/theta physics is gated by tests/study/staggered_zh_ghostproj.py. Single-rank only (ghost projection v1). Pattern: scripts/verify_poiseuille_flow.py.
| verify_poiseuille_ghostproj.channel_sdf | ( | nx, | |
| ny, | |||
| nz, | |||
| ylo, | |||
| yhi | |||
| ) |
Definition at line 24 of file verify_poiseuille_ghostproj.py.
Referenced by run().
| verify_poiseuille_ghostproj.run | ( | N, | |
| ghost, | |||
rho = 1.0, |
|||
mu = 0.1, |
|||
dt = 50.0, |
|||
F = 0.01, |
|||
max_steps = 400 |
|||
| ) |
Definition at line 31 of file verify_poiseuille_ghostproj.py.
References channel_sdf().
Referenced by main().
| verify_poiseuille_ghostproj.main | ( | ) |
Definition at line 65 of file verify_poiseuille_ghostproj.py.
Referenced by main().