Dimensional analysis of two-phase flow including a rate-dependent capillary pressure-saturation relationship

被引:44
作者
Manthey, Sabine [1 ]
Hassanizadeh, S. Majid [2 ]
Helmig, Rainer [1 ]
Hilfer, Rudolf [3 ]
机构
[1] Univ Stuttgart, Inst Hydraul Engn, Res Grp Hydromech & Modelling Hydrosyst, D-70569 Stuttgart, Germany
[2] Univ Utrecht, Fac Geosci, Environm Hydrogeol Grp, NL-3508 TA Utrecht, Netherlands
[3] Univ Stuttgart, Inst Computat Phys, D-70569 Stuttgart, Germany
关键词
two-phase flow; capillary pressure-saturation relationship; rate dependence; dimensional analysis; dimensionless numbers;
D O I
10.1016/j.advwatres.2008.01.021
中图分类号
TV21 [水资源调查与水利规划];
学科分类号
081501 ;
摘要
The macroscopic modelling of two-phase flow processes in subsurface hydrosystems or industrial applications on the Darcy scale usually requires a constitutive relationship between capillary pressure and saturation, the P-c(S-w) relationship. Traditionally, it is assumed that a unique relation between P-c and S-w exists independently of the flow conditions as long as hysteretic effects can be neglected. Recently, this assumption has been questioned and alternative formulations have been suggested. For example, the extended P-c(S-w) relationship by Hassanizadeh and Gray [Hassanizadeh SM, Gray WG. Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries. Adv Water Resources 1990;13(4):169-86] proposes that the difference between the phase pressures to the equilibrium capillary pressure is a linear function of the rate of change of saturation, thereby introducing a constant of proportionality, the coefficient tau. It is desirable to identify cases where the extended relationship needs to be considered. Consequently, a dimensional analysis is performed on the basis of the two-phase balance equations. In addition to the well-known capillary and gravitational number, the dimensional analysis yields a new dimensionless number. The dynamic number Dy quantifies the ratio of dynamic capillary to viscous forces. Relating the dynamic to the capillary as well as the gravitational number gives the new numbers DyC and DyG, respectively. For given sets of fluid and porous medium parameters, the dimensionless numbers Dy and DyC are interpreted as functions of the characteristic length and flow velocity. The simulation of an imbibition process provides insight into the interpretation of the characteristic length scale. The most promising choice for this length scale seems to be the front width. We conclude that consideration of the extended P-c(S-w) relationship may be important for porous media with high permeability, small entry pressure and high coefficient tau when systems with a small characteristic length (e.g. steep front) and small characteristic time scale are under investigation. (C) 2008 Elsevier Ltd. All rights reserved.
引用
收藏
页码:1137 / 1150
页数:14
相关论文
共 50 条
[1]  
[Anonymous], 1999, NATURAL ATTENUATION
[2]  
[Anonymous], ADV CHEM PHYS
[3]  
[Anonymous], 1978, IAHR S SCAL EFF POR
[4]  
[Anonymous], 1964, HYDRAULIC PROPERTIES
[5]   Trapping and mobilization of residual fluid during capillary desaturation in porous media [J].
Anton, L ;
Hilfer, R .
PHYSICAL REVIEW E, 1999, 59 (06) :6819-6823
[6]  
Barenblatt G. I., 1987, Journal of Engineering Physics, V52, P335, DOI 10.1007/BF00872519
[7]  
Barenblatt G. I., 1997, APPL ANAL, V65, P19, DOI DOI 10.1080/00036819708840547
[8]  
BARENBLATT GI, 2002, SPE DOE 13 S IMPR OI
[9]  
Bialas Z., 1970, ARCH HYDROTECH, P407
[10]   Experimental evidence of nonlocal hydrodynamic influence on the dynamic contact angle [J].
Blake, TD ;
Bracke, M ;
Shikhmurzaev, YD .
PHYSICS OF FLUIDS, 1999, 11 (08) :1995-2007