demo_network_flow_packing.py File Reference

Go to the source code of this file.

Namespaces
namespace	demo_network_flow_packing

Variables
int	demo_network_flow_packing.vti = sys.argv[1] if len(sys.argv) > 1 else "../flow/data/packing_128.vti"
	demo_network_flow_packing.sdf_zyx
	demo_network_flow_packing.origin_zyx
	demo_network_flow_packing.spacing_zyx
	demo_network_flow_packing.nz
	demo_network_flow_packing.ny
	demo_network_flow_packing.nx
	demo_network_flow_packing.MU
	demo_network_flow_packing.FX
	demo_network_flow_packing.sdf_xyz = np.asfortranarray(sdf_zyx.T.astype(np.float64))
	demo_network_flow_packing.s = peclet.flow.Solver(nx, ny, nz)
	demo_network_flow_packing.cutcell_pressure
	demo_network_flow_packing.t0 = time.time()
	demo_network_flow_packing.u
	demo_network_flow_packing.v
	demo_network_flow_packing.w
	demo_network_flow_packing.p
	demo_network_flow_packing.ox
	demo_network_flow_packing.oy
	demo_network_flow_packing.oz
	demo_network_flow_packing.F = np.array([(ox[i] * u[i]).sum() for i in range(nx)]) * spacing_zyx[1] * spacing_zyx[0]
	demo_network_flow_packing.net
	demo_network_flow_packing.Q = np.array(net["throat_flow"])
	demo_network_flow_packing.dp = np.array(net["throat_dp"])
	demo_network_flow_packing.A = np.array(net["throat_area"])
	demo_network_flow_packing.res = np.array(net["pore_residual"])
	demo_network_flow_packing.g = Q / dp
int	demo_network_flow_packing.ok = np.abs(res).max() < 1e-6 * np.abs(Q).max() and np.isfinite(g).all()