def __init__(self,
shape,
spacing=1,
label_1='primary',
label_2='secondary',
**kwargs):
super().__init__(**kwargs)
spacing = sp.array(spacing)
shape = sp.array(shape)
# Deal with non-3D shape arguments
shape = sp.pad(shape, [0, 3 - shape.size],
mode='constant',
constant_values=1)
net = Cubic(shape=shape, spacing=1)
net['throat.' + label_1] = True
net['pore.' + label_1] = True
single_dim = shape == 1
shape[single_dim] = 2
dual = Cubic(shape=shape - 1, spacing=1)
faces = [['front', 'back'], ['left', 'right'], ['top', 'bottom']]
faces = [faces[i] for i in sp.where(~single_dim)[0]]
faces = sp.array(faces).flatten().tolist()
dual.add_boundary_pores(faces)
# Add secondary network name as a label
dual['pore.' + label_2] = True
dual['throat.' + label_2] = True
# Shift coordinates prior to stitching
dual['pore.coords'] += 0.5 * (~single_dim)
topotools.stitch(net,
dual,
P_network=net.Ps,
P_donor=dual.Ps,
len_max=1)
net['throat.interconnect'] = net['throat.stitched']
del net['throat.stitched']
# Clean-up labels
net['pore.surface'] = False
net['throat.surface'] = False
for face in faces:
# Remove face label from secondary network since it's internal now
Ps = net.pores(labels=[face, label_2], mode='xnor')
net['pore.' + face][Ps] = False
Ps = net.pores(labels=[face + '_boundary'])
net['pore.' + face][Ps] = True
Ps = net.pores(face)
net['pore.surface'][Ps] = True
Ts = net.find_neighbor_throats(pores=Ps, mode='xnor')
net['throat.surface'][Ts] = True
net['throat.' + face] = net.tomask(throats=Ts)
[net.pop(item) for item in net.labels() if 'boundary' in item]
# Label non-surface pores and throats as internal
net['pore.internal'] = True
net['throat.internal'] = True
# Transfer all dictionary items from 'net' to 'self'
[self.update({item: net[item]}) for item in net]
ws.close_project(net.project)
# Finally, scale network to requested spacing
net['pore.coords'] *= spacing
np.random.seed(10)
# Define network, geometry and constituent phases
shape = np.array([100, 100, 1])
net = Cubic(shape=shape)
geom = StickAndBall(network=net, pores=net.Ps, throats=net.Ts)
air = Air(network=net, name="air")
water = Water(network=net, name="water")
water["pore.diffusivity"] = air["pore.diffusivity"] * 0.05
# Define the regions to be occupied by the two phases (air and water)
x, y, z = net["pore.coords"].T
ps_water = net.Ps[(y >= 75) + (y <= 25)]
ps_air = np.setdiff1d(net.Ps, ps_water)
ts_water = net.find_neighbor_throats(pores=ps_water, mode="xnor")
ts_air = net.find_neighbor_throats(pores=ps_air, mode="xnor")
ts_interface = net.find_neighbor_throats(pores=ps_water, mode="xor")
# Define multiphase and set phase occupancy
mphase = MultiPhase(network=net, phases=[air, water], name="mphase")
mphase._set_automatic_throat_occupancy()
mphase.set_occupancy(air, pores=ps_air, throats=ts_air)
mphase.set_occupancy(water, pores=ps_water, throats=ts_water)
# Define physics
phys = Standard(network=net, phase=mphase, geometry=geom)
# Assign a partition coefficient (concentration ratio)
K_water_air = 0.5 # c @ water / c @ air
mphase.set_binary_partition_coef(phases=[water, air],
model=constant,