staggered_zh_ghostproj_modes Namespace Reference

Functions
	lattice_sdf (N, phi=0.125)
	drag (N, morder, rorder, mu=0.1, F=1e-3, dt=80.0, warm_tol=1e-7, tail=40, max_steps=4000)

Variables
float	kref = 4.2920
dict	CUT = {32: -0.314, 48: -0.082, 64: -0.018, 96: +0.009, 128: +0.013}
dict	G22 = {32: (-0.609, 11), 48: (-0.256, 13), 64: (-0.144, 13), 96: (-0.066, 16), 128: (-0.039, 15)}
	flush
dict	prev = {}
	K11
	s11
	i11
	d11
	t11
	K12
	s12
	i12
	d12
	t12
int	e11 = 100 * (K11 - kref) / kref
int	e12 = 100 * (K12 - kref) / kref
dict	o11 = np.log(abs(prev["e11"]) / abs(e11)) / np.log(N / prev["N"]) if prev else float("nan")
dict	o12 = np.log(abs(prev["e12"]) / abs(e12)) / np.log(N / prev["N"]) if prev else float("nan")
	e22
	i22

Detailed Description

Ghost-projection closure-order comparison on staggered Zick & Homsy drag:
  (2,2) full quadratic  — 13-point nonsymmetric matrix (baseline ghost mode, cached below)
  (1,1) linear closure  — 7-point matrix, 1st-order closure everywhere
  (1,2) MIXED/deferred  — quadratic RHS (2nd-order steady constraint) on the 7-point linear
        matrix; the operator mismatch converges through the time stepping (rate ~0.4 measured
        a-priori). Expectation: (1,2) matches (2,2)'s accuracy with fewer BiCGStab iterations;
        (1,1) shows whether the linear closure's larger O(h^2) constant (or 1st-order term)
        degrades the drag. Cached columns from tests/study/staggered_zh_ghostproj.py runs.

Function Documentation

lattice_sdf()

staggered_zh_ghostproj_modes.lattice_sdf	(		N,
			phi = 0.125
	)

Definition at line 18 of file staggered_zh_ghostproj_modes.py.

Referenced by drag().

drag()

staggered_zh_ghostproj_modes.drag	(		N,
			morder,
			rorder,
			mu = 0.1,
			F = 1e-3,
			dt = 80.0,
			warm_tol = 1e-7,
			tail = 40,
			max_steps = 4000
	)

Definition at line 31 of file staggered_zh_ghostproj_modes.py.

References lattice_sdf().

Variable Documentation

kref

float staggered_zh_ghostproj_modes.kref = 4.2920

Definition at line 61 of file staggered_zh_ghostproj_modes.py.

CUT

dict staggered_zh_ghostproj_modes.CUT = {32: -0.314, 48: -0.082, 64: -0.018, 96: +0.009, 128: +0.013}

Definition at line 63 of file staggered_zh_ghostproj_modes.py.

G22

dict staggered_zh_ghostproj_modes.G22 = {32: (-0.609, 11), 48: (-0.256, 13), 64: (-0.144, 13), 96: (-0.066, 16), 128: (-0.039, 15)}

Definition at line 64 of file staggered_zh_ghostproj_modes.py.

flush

staggered_zh_ghostproj_modes.flush

Definition at line 67 of file staggered_zh_ghostproj_modes.py.

prev

dict staggered_zh_ghostproj_modes.prev = {}

Definition at line 70 of file staggered_zh_ghostproj_modes.py.

K11

staggered_zh_ghostproj_modes.K11

Definition at line 72 of file staggered_zh_ghostproj_modes.py.

s11

staggered_zh_ghostproj_modes.s11

Definition at line 72 of file staggered_zh_ghostproj_modes.py.

i11

staggered_zh_ghostproj_modes.i11

Definition at line 72 of file staggered_zh_ghostproj_modes.py.

d11

staggered_zh_ghostproj_modes.d11

Definition at line 72 of file staggered_zh_ghostproj_modes.py.

t11

staggered_zh_ghostproj_modes.t11

Definition at line 72 of file staggered_zh_ghostproj_modes.py.

K12

staggered_zh_ghostproj_modes.K12

Definition at line 73 of file staggered_zh_ghostproj_modes.py.

s12

staggered_zh_ghostproj_modes.s12

Definition at line 73 of file staggered_zh_ghostproj_modes.py.

i12

staggered_zh_ghostproj_modes.i12

Definition at line 73 of file staggered_zh_ghostproj_modes.py.

d12

staggered_zh_ghostproj_modes.d12

Definition at line 73 of file staggered_zh_ghostproj_modes.py.

t12

staggered_zh_ghostproj_modes.t12

Definition at line 73 of file staggered_zh_ghostproj_modes.py.

e11

int staggered_zh_ghostproj_modes.e11 = 100 * (K11 - kref) / kref

Definition at line 74 of file staggered_zh_ghostproj_modes.py.

e12

int staggered_zh_ghostproj_modes.e12 = 100 * (K12 - kref) / kref

Definition at line 75 of file staggered_zh_ghostproj_modes.py.

o11

dict staggered_zh_ghostproj_modes.o11 = np.log(abs(prev["e11"]) / abs(e11)) / np.log(N / prev["N"]) if prev else float("nan")

Definition at line 76 of file staggered_zh_ghostproj_modes.py.

o12

dict staggered_zh_ghostproj_modes.o12 = np.log(abs(prev["e12"]) / abs(e12)) / np.log(N / prev["N"]) if prev else float("nan")

Definition at line 77 of file staggered_zh_ghostproj_modes.py.

e22

staggered_zh_ghostproj_modes.e22

Definition at line 78 of file staggered_zh_ghostproj_modes.py.

i22

staggered_zh_ghostproj_modes.i22

Definition at line 78 of file staggered_zh_ghostproj_modes.py.