예제 #1
0
    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)
예제 #3
0
    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)
예제 #4
0
    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)