verify_lid_cavity_sdflow Namespace Reference

Functions
	run (N=128, Re=100.0, U=1.0, nz=4, max_steps=5000)
	main ()

Variables
	GHIA_Y
	GHIA_U
	GHIA_X
	GHIA_V

Detailed Description

Verification (sdflow): the lid-driven cavity -- the canonical NATIVE-domain-BC benchmark. Three
no-slip walls (-x,+x,-y) + a lid (+y) moving in +x; quasi-2D (periodic z). Driven only by the lid, so
this exercises the whole BC framework: Dirichlet/no-slip velocity ghosts (mac_bc.cuh), a non-periodic
halo on x and y, and Neumann pressure on the walls (boundary-face openness zeroed). We compare the
centreline profiles to the tabulated Ghia, Ghia & Shin (1982) data at Re=100.

Uses the canonical `sdflow` module. NO immersed solid -- the cavity is set up with set_domain_bc +
set_pressure_geometry(all-fluid). Physical units: set_rho/set_mu; Re = U_lid * L / nu, L = N (grid units).

Function Documentation

run()

verify_lid_cavity_sdflow.run	(		N = 128,
			Re = 100.0,
			U = 1.0,
			nz = 4,
			max_steps = 5000
	)

Run the lid-driven cavity to steady state and return the comparison to Ghia et al.

Sets three no-slip walls and a moving lid (Dirichlet velocity U on +y), marches to steady state, and
returns (N, Re, u-rms, v-rms, min centreline u, max divergence, steps) where the rms values are the
centreline-profile errors against the tabulated Ghia, Ghia & Shin (1982) data at the given Reynolds
number. Quasi-2D (periodic z). Returns None on non-root MPI ranks.

Definition at line 31 of file verify_lid_cavity_sdflow.py.

Referenced by main().

main()

verify_lid_cavity_sdflow.main

(

)

Run the default cavity case, print the Ghia comparison, and exit non-zero if it fails the tolerance.

Definition at line 72 of file verify_lid_cavity_sdflow.py.

References main(), and run().

Referenced by main().

Variable Documentation

GHIA_Y

verify_lid_cavity_sdflow.GHIA_Y

Initial value:

=  np.array([0, .0547, .0625, .0703, .1016, .1719, .2813, .4531, .5, .6172, .7344, .8516, .9531,
                   .9609, .9688, .9766, 1])

Definition at line 21 of file verify_lid_cavity_sdflow.py.

GHIA_U

verify_lid_cavity_sdflow.GHIA_U

Initial value:

=  np.array([0, -.03717, -.04192, -.04775, -.06434, -.10150, -.15662, -.21090, -.20581, -.13641,
                   .00332, .23151, .68717, .73722, .78871, .84123, 1])

Definition at line 23 of file verify_lid_cavity_sdflow.py.

GHIA_X

verify_lid_cavity_sdflow.GHIA_X

Initial value:

=  np.array([0, .0625, .0703, .0781, .0938, .1563, .2266, .2344, .5, .8047, .8594, .9063, .9453,
                   .9531, .9609, .9688, 1])

Definition at line 25 of file verify_lid_cavity_sdflow.py.

GHIA_V

verify_lid_cavity_sdflow.GHIA_V

Initial value:

=  np.array([0, .09233, .10091, .10890, .12317, .16077, .17507, .17527, .05454, -.24533, -.22445,
                   -.16914, -.10313, -.08864, -.07391, -.05906, 0])

Definition at line 27 of file verify_lid_cavity_sdflow.py.