def test_Erhg(self):
vls = 3.0
Dp = 0.5
d = 0.075 / 1000
epsilon = DHLLDV_constants.steel_roughness
rhol = DHLLDV_constants.water_density[20]
nu = DHLLDV_constants.water_viscosity[20]
rhos = 2.65
Cvs = 0.25
self.assertAlmostEqual(
homogeneous.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
0.01187623153019)
def test_Erhg_med_sand(self):
vls = 3.0
Dp = 0.5
d = 0.4 / 1000
epsilon = DHLLDV_constants.steel_roughness
rhol = DHLLDV_constants.water_density[20]
nu = DHLLDV_constants.water_viscosity[20]
rhos = 2.65
Cvs = 0.25
self.assertAlmostEqual(
homogeneous.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
0.00633356768663)
def test_Erhg_sf(self):
vls = 3.0
Dp = 0.5
d = 8.0 / 1000
epsilon = DHLLDV_constants.steel_roughness
rhol = DHLLDV_constants.water_density[20]
nu = DHLLDV_constants.water_viscosity[20]
rhos = 2.65
Cvs = 0.25
self.assertAlmostEqual(
homogeneous.Erhg(vls,
Dp,
d,
epsilon,
nu,
rhol,
rhos,
Cvs,
use_sf=False), 0.00465260820521)
self.assertAlmostEqual(
homogeneous.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
0.0302993)
def Cvs_Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs, get_dict=False):
"""
Cvs_Erhg - Calculate the Erhg for the given slurry, using the appropriate model
vls = average line speed (velocity, m/sec)
Dp = Pipe diameter (m)
d = Particle diameter (m)
epsilon = absolute pipe roughness (m)
nu = fluid kinematic viscosity in m2/sec
rhol = density of the fluid (ton/m3)
rhos = particle density (ton/m3)
Cvs = insitu volume concentration
get_dict: if true return the dict with all models.
"""
Erhg_obj = {
'il':
homogeneous.fluid_head_loss(vls, Dp, epsilon, nu, rhol),
'FB':
stratified.fb_Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
'SB':
stratified.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
'He':
heterogeneous.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs, use_sf,
use_sqrtcx),
'Ho':
homogeneous.Erhg(vls, Dp, d, epsilon, nu, rhol, rhos, Cvs),
}
if Erhg_obj['FB'] < Erhg_obj['SB']:
regime = 'FB'
else:
regime = 'SB'
if Erhg_obj[regime] > Erhg_obj['He']:
regime = 'He'
if Erhg_obj[regime] < Erhg_obj['Ho']:
regime = 'Ho'
Erhg_obj['regime'] = regime
if get_dict:
return Erhg_obj
else:
return Erhg_obj[Erhg_obj['regime']]
