verify_periodic_spheres_sdflow Namespace Reference

Functions
	packing_sdf (N, radius_frac=0.18)
	run (N, rho=1.0, mu=0.1, dt=60.0, F=1e-3, max_steps=200)
	main ()

Detailed Description

Verification (sdflow): creeping (Stokes) flow through a periodic sphere packing -- the porous-media
target case. A body force drives flow through the pore space; the Robust-Scaled cut-cell IBM enforces
no-slip on the spheres and the cut-cell pressure operator (Galerkin multigrid + CG) enforces
incompressibility. We report the Darcy permeability k = mu*<u>/F and check:
  * the flow is incompressible (small cut-cell flux divergence),
  * no-slip holds (velocity ~ 0 inside the solid),
  * the permeability is finite/positive and increases as the spheres resolve.

Uses the canonical `sdflow` module (one GPU as plain `python`, or multi-rank under `mpirun -np N python`).
Physical units: set_rho/set_mu, body force F is a force per unit volume (= -dp/dx). Grid units (dx = 1).

Function Documentation

packing_sdf()

verify_periodic_spheres_sdflow.packing_sdf	(		N,
			radius_frac = 0.18
	)

2x2x2 sphere packing as a 3-D array sdf[x,y,z]; negative inside the spheres (min-image).

Definition at line 22 of file verify_periodic_spheres_sdflow.py.

Referenced by run().

run()

verify_periodic_spheres_sdflow.run	(		N,
			rho = 1.0,
			mu = 0.1,
			dt = 60.0,
			F = 1e-3,
			max_steps = 200
	)

Definition at line 39 of file verify_periodic_spheres_sdflow.py.

References packing_sdf().

Referenced by main().

main()

verify_periodic_spheres_sdflow.main

(

)

Definition at line 78 of file verify_periodic_spheres_sdflow.py.

References main(), and run().

Referenced by main().