def __init__(
self,
M=10,
U=1.1, # [rps] impeller revolutions
phi=0.091, # [1] holdup
v0=5e-10, # [cm^3]
model_parameters=None,
theta=600.):
# pipe diamter and length
self.D = 24.0e-03
self.L = 9.5
contProperties = dict()
# oil
contProperties['mu'] = 1.6e-3 # [P = kg * m^-1 s^-1]
contProperties['rho'] = 801. # [kg/cm3]
# calculate turbulent properties
Re = U * self.D / contProperties['mu'] * contProperties['rho']
I = 0.16 * Re**(-1. / 8.)
u_rms = U * I
k = 3. / 2. * u_rms**2
L_t = 0.038 * self.D
contProperties['epsilon'] = 0.09 * k**(3. / 2.) / L_t
contProperties['Re'] = Re
# water solution
dispersion = DispersionProperties(
phi=phi,
rho=1000., # [kg/m3]
sigma=1.7e-2, # [P = kg * m^-1 s^-1]
v_max=v0 * 3,
v0=v0,
sigma0=v0 / 10)
# Feed distribution
# Feed
domain = DomainProperties(theta=theta,
V=pi * self.L * (self.D / 2)**2,
M=M)
CaseSolution.__init__(self,
dispersion,
contProperties,
domain,
model_parameters=model_parameters)
def __init__(
self,
M=10,
U=2.71, # [rps] impeller revolutions
phi=0.117, # [1] holdup
v0=4e-11, # [cm^3]
model_parameters=None,
theta=600.):
# pipe diamter andlength
self.D = 0.063 # [m] impeller diameter
self.L = 4.5 # [m] impeller diameter
contProperties = dict()
# oil
contProperties['mu'] = 1.8e-3 # [P = kg * m^-1 s^-1]
contProperties['rho'] = 797. # [kg/cm3]
# calculate turbulent properties
Re = U * self.D / contProperties['mu'] * contProperties['rho']
I = 0.16 * Re**(-1. / 8.)
u_rms = U * I
k = 3. / 2. * u_rms**2
L_t = 0.038 * self.D
contProperties['epsilon'] = 0.09 * k**(3. / 2.) / L_t
contProperties['Re'] = Re
# water solution
dispersion = DispersionProperties(
phi=phi,
rho=1166., # [kg/m3]
sigma=1.e-2, # [P = kg * m^-1 s^-1]
v_max=6e-11,
v0=v0,
sigma0=v0 / 10)
domain = DomainProperties(theta=theta,
V=pi * self.L * (self.D / 2)**2,
M=M)
CaseSolution.__init__(self,
dispersion,
contProperties,
domain,
model_parameters=model_parameters)
def __init__(
self,
M=10,
U=1.1, # [rps] impeller revolutions
phi=0.091, # [1] holdup
v0=5e-10, # [cm^3]
model_parameters=None,
theta=600.):
# pipe diamter and length
self.D = 24.0e-03
self.L = 9.5
contProperties = dict()
# oil
contProperties['mu'] = 1.6e-3 # [P = kg * m^-1 s^-1]
contProperties['rho'] = 801. # [kg/cm3]
# calculate turbulent properties
Re = U * self.D / contProperties['mu'] * contProperties['rho']
I = 0.16 * Re ** (-1. / 8.)
u_rms = U * I
k = 3. / 2. * u_rms ** 2
L_t = 0.038 * self.D
contProperties['epsilon'] = 0.09 * k ** (3. / 2.) / L_t
contProperties['Re'] = Re
# water solution
dispersion = DispersionProperties(
phi=phi, rho=1000., # [kg/m3]
sigma=1.7e-2, # [P = kg * m^-1 s^-1]
v_max=v0 * 3,
v0=v0,
sigma0=v0 / 10)
# Feed distribution
# Feed
domain = DomainProperties(
theta=theta, V=pi * self.L * (self.D / 2) ** 2, M=M)
CaseSolution.__init__(
self, dispersion, contProperties, domain,
model_parameters=model_parameters)
def __init__(
self,
M=10,
Nstar=4.16, # [rps] impeller revolutions
phi=0.15, # [1] holdup
v0=0.03):
self.D = 10 # [cm] impeller diameter
contProperties = dict()
dispProperties = dict()
domainProperties = dict()
# Water
contProperties['mu'] = 0.89e-2 # [P = g * cm^-1 s^-1]
contProperties['rho'] = 0.1 # [g/cm3]
contProperties['epsilon'] = 0.407 * Nstar * self.D
# Kerosene-dicholorebenzene
dispProperties['sigma'] = 42.82 # [P = g * cm^-1 s^-1]
dispProperties['rho'] = 0.972 # [g/cm3]
dispProperties['phi'] = phi
mm3_to_cm3 = 0.1**3
dispProperties['vMax'] = 0.08 * mm3_to_cm3
# Feed distribution
dispProperties['v0'] = v0 * mm3_to_cm3
dispProperties['sigma0'] = 0.005 * mm3_to_cm3
# Feed
domainProperties['theta'] = 600.
domainProperties['V'] = 12 * 10**3
domainProperties['M'] = M
Ninit = zeros(M)
CaseSolution.__init__(
self, dispProperties, contProperties, domainProperties)