rcp_permeability_ab Namespace Reference

Functions
	sdf_from_pack (Ng, pos, r, side)
	permeability (Ng, sdf, side, colloc, ghost, mode=0, mu=0.1, F=1e-3, dt=80.0, max_steps=3000, tol=1e-6)

Variables
	HERE = os.path.dirname(os.path.abspath(__file__))
	d = np.load(os.path.join(HERE, "rcp_pack_seed3.npz"))
	pos
	r
	side
	phi
	flush
list	variants
	sdf = sdf_from_pack(Ng, pos, r, side)
str	row = f"{Ng:>4} |"
	kref = None
	t0 = time.time()
	out = permeability(Ng, sdf, side, colloc, ghost, mode)

Detailed Description

RCP permeability A/B: {staggered, collocated} x {cutcell, ghost} on the SAME random close
packing (the peclet-examples random-packed-bed geometry: dem pack_bed, N=180, phi~0.63, seed 3;
cached in rcp_pack_seed3.npz — regenerate with the example's pack_bed if absent).

Context: the staggered examples-check found the ghost projection inflates the RCP permeability
(+21/+15/+12% at Ng=32/44/56 vs the cutcell reference) because point-based closure faces cannot
throttle sub-cell throats the way apertures do, and the binary COUPLED graph fragments (guard
decouples the pockets). This study measures whether the collocated ghost inherits the same
defect, with the collocated cutcell (mode 0) and both staggered variants as anchors. Protocol =
the example's steady loop with a tighter stop (tol 1e-6) so the scheme, not the stop criterion,
is compared.

Function Documentation

sdf_from_pack()

rcp_permeability_ab.sdf_from_pack	(		Ng,
			pos,
			r,
			side
	)

Definition at line 23 of file rcp_permeability_ab.py.

permeability()

rcp_permeability_ab.permeability	(		Ng,
			sdf,
			side,
			colloc,
			ghost,
			mode = 0,
			mu = 0.1,
			F = 1e-3,
			dt = 80.0,
			max_steps = 3000,
			tol = 1e-6
	)

Definition at line 38 of file rcp_permeability_ab.py.

Variable Documentation

HERE

rcp_permeability_ab.HERE = os.path.dirname(os.path.abspath(__file__))

Definition at line 20 of file rcp_permeability_ab.py.

d

rcp_permeability_ab.d = np.load(os.path.join(HERE, "rcp_pack_seed3.npz"))

Definition at line 70 of file rcp_permeability_ab.py.

pos

rcp_permeability_ab.pos

Definition at line 71 of file rcp_permeability_ab.py.

r

rcp_permeability_ab.r

Definition at line 71 of file rcp_permeability_ab.py.

side

rcp_permeability_ab.side

Definition at line 71 of file rcp_permeability_ab.py.

phi

rcp_permeability_ab.phi

Definition at line 71 of file rcp_permeability_ab.py.

flush

rcp_permeability_ab.flush

Definition at line 72 of file rcp_permeability_ab.py.

variants

list rcp_permeability_ab.variants

Initial value:

=  [("stag cutcell", False, False, 0), ("stag ghost", False, True, 0),
                ("col cutcell", True, False, 0), ("col ghost", True, True, 0),
                ("col hyb9", True, False, 9)]

Definition at line 75 of file rcp_permeability_ab.py.

sdf

rcp_permeability_ab.sdf = sdf_from_pack(Ng, pos, r, side)

Definition at line 81 of file rcp_permeability_ab.py.

row

rcp_permeability_ab.row = f"{Ng:>4} |"

Definition at line 82 of file rcp_permeability_ab.py.

kref

rcp_permeability_ab.kref = None

Definition at line 83 of file rcp_permeability_ab.py.

t0

rcp_permeability_ab.t0 = time.time()

Definition at line 85 of file rcp_permeability_ab.py.

out

rcp_permeability_ab.out = permeability(Ng, sdf, side, colloc, ghost, mode)

Definition at line 86 of file rcp_permeability_ab.py.