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

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

Detailed Description

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.

Function Documentation

◆ channel_sdf()

verify_poiseuille_ghostproj.channel_sdf (   nx,
  ny,
  nz,
  ylo,
  yhi 
)

Definition at line 24 of file verify_poiseuille_ghostproj.py.

Referenced by run().

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

◆ main()

verify_poiseuille_ghostproj.main ( )

Definition at line 65 of file verify_poiseuille_ghostproj.py.

References main(), and run().

Referenced by main().