valvatne_blunt#

valvatne_blunt(phase, pore_viscosity='pore.viscosity', throat_viscosity='throat.viscosity', pore_shape_factor='pore.shape_factor', throat_shape_factor='throat.shape_factor', pore_area='pore.area', throat_area='throat.cross_sectional_area', conduit_lengths='throat.conduit_lengths')[source]#

Calculates the single phase hydraulic conductance of conduits.

Function has been adapted for use with the Statoil imported networks and makes use of the shape factor in these networks to apply Hagen-Poiseuille flow for conduits of different shape classes: triangular, square and circular [2].

Parameters:

  • phase (OpenPNM Phase object) – The phase object to which this model is associated (i.e. attached). This controls the length of the calculated array(s), and also provides access to other necessary properties.

  • pore_viscosity (str) – Name of the dictionary key on the target object pointing to the array containing values of pore viscosity

  • throat_viscosity (str) – Name of the dictionary key on the target object pointing to the array containing values of throat viscosity

  • pore_shape_factor (str) – Name of the dictionary key on the target object pointing to the array containing values of pore geometric shape factor

  • throat_shape_factor (str) – Name of the dictionary key on the target object pointing to the array containing values of throat geometric shape factor

  • pore_area (str) –

    Name of the dictionary key on the target object pointing to the array containing values of pore area The pore area is calculated using following formula:

    \[A_P = \frac{R_P^2}{(4 \cdot SF_P)}\]

    where theoratical value of pore_shape_factor in a circular tube is calculated using following formula:

    \[SF_P = \frac{A_P}{P_P^2} = 1/4π\]

  • throat_area (str) –

    Name of the dictionary key on the target object pointing to the array containing values of throat area. The throat area is calculated using following formula:

    \[T_A = \frac{R_T^2}{(4 \cdot SF_T)}\]

    where theoratical value of throat shape factor in circular tube is calculated using :

    \[SF_T = \frac{T_A}{T_P^2} = 1/4π\]

  • conduit_lengths (str) – Name of the dictionary key on the target object pointing to the array containing values of throat conduit lengths

Returns:

values – A numpy ndarray containing the computed values of

Return type:

ndarray

References

[1] Valvatne, Per H., and Martin J. Blunt. “Predictive pore‐scale modeling of two‐phase flow in mixed wet media.” Water Resources Research 40, no. 7 (2004).

[2] Patzek, T. W., and D. B. Silin (2001), Shape factor and hydraulic conductance in noncircular capillaries I. One-phase creeping flow, J. Colloid Interface Sci., 236, 295–304.

On this page