def __init__( self, M=10, Nstar=4.16, # [rps] impeller revolutions phi=0.15, # [1] holdup v0=0.05): self.D = 10 # [cm] impeller diameter self.Nstar = Nstar self.phi = phi # Water self.muc = 0.89e-2 # [P = g * cm^-1 s^-1] self.rhoc = 0.1 # [g/cm3] # Kerosene-dicholorebenzene self.rhod = 0.972 # [g/cm3] self.sigma = 42.82 time = arange(0.0, 3600, 1) mm3_to_cm3 = 0.1**3 vmax = 0.06 * mm3_to_cm3 # Feed distribution self.v0 = v0 * mm3_to_cm3 self.sigma0 = 0.005 * mm3_to_cm3 # Feed theta = 600 self.Vt = 12 * 10**3 self.n0 = self.Vt / theta Ninit = zeros(M) MOCSolution.__init__( self, Ninit, time, vmax / M, beta=beta, gamma=self.g, Q=self.Q, theta=theta, n0=self.n0, A0=self.A0)
def __init__( self, dispProperties, contProperties, domainProperties, modelParameters=None): self.dispProperties = dispProperties self.contProperties = contProperties self.phi = dispProperties['phi'] self.muc = contProperties['mu'] self.rhoc = contProperties['rho'] self.epsilon = contProperties['epsilon'] self.rhod = dispProperties['rho'] self.sigma = dispProperties['sigma'] time = arange(0.0, 3600, 1) vmax = dispProperties['vMax'] # Feed distribution self.v0 = dispProperties['v0'] self.sigma0 = dispProperties['sigma0'] # Feed theta = domainProperties['theta'] M = domainProperties['M'] self.Vt = domainProperties['V'] if theta is None: self.n0 = None else: self.n0 = self.Vt / theta Ninit = zeros(M) if modelParameters is None: self.C1 = 0.4 self.C2 = 0.08 self.C3 = 2.8e-6 self.C4 = 1.83e9 else: self.C1 = modelParameters['C1'] self.C2 = modelParameters['C2'] self.C3 = modelParameters['C3'] self.C4 = modelParameters['C4'] MOCSolution.__init__( self, Ninit, time, vmax / M, beta=beta, gamma=self.g, Q=self.Q, theta=theta, n0=self.n0, A0=self.A0)
def __init__( self, M=10, Nstar=4.16, # [rps] impeller revolutions phi=0.15, # [1] holdup vmax=0.08): self.D = 10 # [cm] impeller diameter self.Nstar = Nstar self.phi = phi # Water self.muc = 0.89e-2 # [P = g * cm^-1 s^-1] self.rhoc = 1.0 # [g/cm3] # Kerosene-dicholorebenzene self.rhod = 0.972 # [g/cm3] self.sigma = 42.82 time = arange(0.0, 3600, 0.5) mm3_to_cm3 = 0.1**3 vmax = vmax * mm3_to_cm3 # vmax = 0.06 * mm3_to_cm3 # Feed distribution self.v0 = vmax / 2 self.sigma0 = (vmax - self.v0) / 3.3 vmin = None # 5.23e-7 # Feed theta = 600 self.Vt = 12 * 10**3 self.n0 = self.Vt / theta MOCSolution.__init__( self, M, time, vmax / M, N0=self.N0, xi0=vmin, beta=beta, gamma=self.g, Q=self.Qf, theta=theta, n0=self.n0, A0=self.A0)