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

Functions

 inlet_profile (H, S, nz, U_in)
 
 reattachment (u_bottom)
 
 run (Re, S=16, Lr=12, U_in=1.0, nz=4, dt=0.4, max_steps=8000)
 
 main ()
 

Detailed Description

Phase-5b verification (collocated grid): the backward-facing step (BFS) -- the separated-flow benchmark
and the showcase for the collocated per-face INLET PROFILE. Gartling (1990) expansion-ratio-2 step: a
channel of full height H=2S with no-slip walls top/bottom; flow enters ONLY the upper half (y in [S,2S])
with the developed parabola and is zero over the lower half (the step face); it leaves through the outflow.
The step is realized purely as the inlet condition (set_domain_bc_profile -> the collocated profile inlet
ghost, bcVelocityColocated with a per-position value), no immersed solid.

Laminar de-risk (Re_S = U_in*S/nu = 100, 200): a recirculation bubble must form behind the step, the
reattachment length x_r/S must grow with Re, mass is conserved, and the flow stays incompressible (the open
outflow leaves the approximate projection's O(h^2) residual, looser than the staggered 1e-6). The staggered
solver reaches x_r/S 5.3 (Re=100) -> 8.3 (Re=200) on the Armaly/Biswas curve.

Function Documentation

◆ inlet_profile()

verify_colocated_bfs.inlet_profile (   H,
  S,
  nz,
  U_in 
)

Definition at line 24 of file verify_colocated_bfs.py.

Referenced by run().

◆ reattachment()

verify_colocated_bfs.reattachment (   u_bottom)

Definition at line 33 of file verify_colocated_bfs.py.

Referenced by run().

◆ run()

verify_colocated_bfs.run (   Re,
  S = 16,
  Lr = 12,
  U_in = 1.0,
  nz = 4,
  dt = 0.4,
  max_steps = 8000 
)

Definition at line 43 of file verify_colocated_bfs.py.

References inlet_profile(), and reattachment().

Referenced by main().

◆ main()

verify_colocated_bfs.main ( )

Definition at line 76 of file verify_colocated_bfs.py.

References main(), and run().

Referenced by main().