nan values in throat.hydraulic_conductance

[minzhangut]

I am trying to follow the tutorial here: https://openpnm.org/examples/reference/uncategorized/intro_to_openpnm_advanced.html?highlight=permeability

After I did:
mod = op.models.physics.hydraulic_conductance.hagen_poiseuille phys_water_internal.add_model(propname='throat.hydraulic_conductance', model=mod)

I saw nan values in throat.hydraulic_conductance:

print(phys_water_internal['throat.hydraulic_conductance'])
[nan nan nan ... nan nan nan]

[jgostick]

Can you post more of your script so I can see at which point the process goes astray? Thanks!

[minzhangut]

import warnings
import numpy as np
import scipy as sp
import openpnm as op
%config InlineBackend.figure_formats = ['svg']
import matplotlib.pyplot as plt
%matplotlib inline
np.random.seed(10)
ws = op.Workspace()
ws.settings['loglevel'] = 40
np.set_printoptions(precision=4)
pn = op.network.Cubic(shape=[10, 10, 10], spacing=0.00006, name='net')

pn.add_boundary_pores(labels=['top', 'bottom'])

Ps = pn.pores('*boundary', mode='not')
Ts = pn.throats('*boundary', mode='not')
geom = op.geometry.SpheresAndCylinders(network=pn, pores=Ps, throats=Ts, name='intern')
Ps = pn.pores('*boundary')
Ts = pn.throats('*boundary')
boun = op.geometry.Boundary(network=pn, pores=Ps, throats=Ts, name='boun')

air = op.phases.Air(network=pn)
water = op.phases.Water(network=pn)
water['throat.contact_angle'] = 110
water['throat.surface_tension'] = 0.072

phys_water_internal = op.physics.GenericPhysics(network=pn, phase=water, geometry=geom)
phys_air_internal = op.physics.GenericPhysics(network=pn, phase=air, geometry=geom)
phys_water_boundary = op.physics.GenericPhysics(network=pn, phase=water, geometry=boun)
phys_air_boundary = op.physics.GenericPhysics(network=pn, phase=air, geometry=boun)

def mason_model(target, diameter='throat.diameter', theta='throat.contact_angle',
                sigma='throat.surface_tension', f=0.6667):
    proj = target.project
    network = proj.network
    phase = proj.find_phase(target)
    Dt = network[diameter]
    theta = phase[theta]
    sigma = phase[sigma]
    Pc = 4*sigma*np.cos(f*np.deg2rad(theta))/Dt
    return Pc[phase.throats(target.name)]

mod = op.models.physics.hydraulic_conductance.hagen_poiseuille
phys_water_internal.add_model(propname='throat.hydraulic_conductance',
                              model=mod)
print(phys_water_internal['throat.hydraulic_conductance'])