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

Functions

 channel_sdf (nx, ny, nz, ylo, yhi)
 
 run (SolverCls, N, rho=1.0, mu=0.1, dt=50.0, F=0.01, max_steps=400)
 
 main ()
 

Detailed Description

Phase-2 verification (collocated grid): plane Poiseuille flow through an SDF channel, driven by a body
force, with Robust-Scaled cut-cell IBM no-slip walls and NO pressure projection (cutcell_pressure=False).

Runs the SAME setup on both grids: sdflow.Solver (staggered) and sdflow.SolverColocated (cell-centered).
The steady centreline velocity must match the analytic parabola U_max = F*H^2/(8*mu) and the error must
shrink with resolution on the collocated grid too. The u-component's wall-normal (y) location is the cell
centre on BOTH grids, so the two solvers should agree closely. This exercises the collocated {0,0,0} IBM
offset + implicit diffusion path (advection is ~0 for unidirectional flow; it is stressed in phase 3).

Function Documentation

◆ channel_sdf()

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

Definition at line 20 of file verify_colocated_poiseuille.py.

Referenced by run().

◆ run()

verify_colocated_poiseuille.run (   SolverCls,
  N,
  rho = 1.0,
  mu = 0.1,
  dt = 50.0,
  F = 0.01,
  max_steps = 400 
)

Definition at line 27 of file verify_colocated_poiseuille.py.

References channel_sdf().

Referenced by main().

◆ main()

verify_colocated_poiseuille.main ( )

Definition at line 56 of file verify_colocated_poiseuille.py.

References main(), and run().

Referenced by main().