Terminal Settling Velocity of a Glass Cube

 

Illustrative Example 4.4: Terminal Settling Velocity of a Glass Cube

Problem Statement:
Calculate the terminal settling velocity of a glass cube with an edge dimension of 0.1 mm in a fluid with the following properties:

  • Density of fluid (ρf\rho_f) = 982 kg/m³
  • Viscosity of fluid (μ\mu) = 0.0013 kg/m·s
  • Density of glass (ρs\rho_s) = 2820 kg/m³

Additionally, calculate:

  1. The equivalent volume diameter (dvd_v) of the cube.
  2. The sphericity factor (ϕ\phi) of the cube.

Solution

Step 1: Calculate the Equivalent Volume Diameter (dvd_v)

The equivalent volume diameter (dvd_v) is the diameter of a sphere with the same volume as the cube.

  1. Volume of the cube:

    Vcube=a3=(0.0001)3=1×1012m3V_{\text{cube}} = a^3 = (0.0001)^3 = 1 \times 10^{-12} \, \text{m}^3
  2. Volume of a sphere:

    Vsphere=π6dv3V_{\text{sphere}} = \frac{\pi}{6} d_v^3
  3. Equating the two volumes:

    1×1012=π6dv31 \times 10^{-12} = \frac{\pi}{6} d_v^3 dv3=61×1012π=1.91×1012m3d_v^3 = \frac{6 \cdot 1 \times 10^{-12}}{\pi} = 1.91 \times 10^{-12} \, \text{m}^3 dv=1.91×10123=0.000125m(or 0.125 mm)d_v = \sqrt[3]{1.91 \times 10^{-12}} = 0.000125 \, \text{m} \, \text{(or 0.125 mm)}

Step 2: Calculate the Sphericity Factor (ϕ\phi)

Sphericity (ϕ\phi) is the ratio of the surface area of a sphere with the same volume to the surface area of the cube.

  1. Surface area of the cube:

    Acube=6a2=6(0.0001)2=6×108m2A_{\text{cube}} = 6 a^2 = 6 (0.0001)^2 = 6 \times 10^{-8} \, \text{m}^2
  2. Surface area of the sphere:

    Asphere=πdv2=π(0.000125)2=4.91×108m2A_{\text{sphere}} = \pi d_v^2 = \pi (0.000125)^2 = 4.91 \times 10^{-8} \, \text{m}^2
  3. Sphericity:

    ϕ=AsphereAcube=4.91×1086×108=0.818\phi = \frac{A_{\text{sphere}}}{A_{\text{cube}}} = \frac{4.91 \times 10^{-8}}{6 \times 10^{-8}} = 0.818

Step 3: Calculate the Terminal Settling Velocity (vtv_t)

Stokes' law for terminal settling velocity is:

vt=2(ρsρf)grv29μv_t = \frac{2 (\rho_s - \rho_f) g r_v^2}{9 \mu}

Where rv=dv2r_v = \frac{d_v}{2}:

rv=0.0001252=0.0000625mr_v = \frac{0.000125}{2} = 0.0000625 \, \text{m}

Substitute the known values:

vt=2(2820982)(9.81)(0.0000625)29(0.0013)v_t = \frac{2 (2820 - 982) (9.81) (0.0000625)^2}{9 (0.0013)}
  1. Density difference:

    ρsρf=2820982=1838kg/m3\rho_s - \rho_f = 2820 - 982 = 1838 \, \text{kg/m}^3
  2. Square the radius:

    rv2=(0.0000625)2=3.91×109m2r_v^2 = (0.0000625)^2 = 3.91 \times 10^{-9} \, \text{m}^2
  3. Numerator:

    218389.813.91×109=1.41×1042 \cdot 1838 \cdot 9.81 \cdot 3.91 \times 10^{-9} = 1.41 \times 10^{-4}
  4. Denominator:

    90.0013=0.01179 \cdot 0.0013 = 0.0117
  5. Terminal velocity:

    vt=1.41×1040.0117=0.0121m/sv_t = \frac{1.41 \times 10^{-4}}{0.0117} = 0.0121 \, \text{m/s}

Results

  1. Equivalent Volume Diameter (dvd_v): 0.000125m(or 0.125 mm)0.000125 \, \text{m} \, \text{(or 0.125 mm)}
  2. Sphericity Factor (ϕ\phi): 0.8180.818
  3. Terminal Settling Velocity (vtv_t): 0.0121m/s(or 12.1 mm/s)

Comments

Popular posts from this blog

Simulation & Modeling Comprehensive Notes

Question bank: Modeling and Simulation of Mineral Processing Systems

Rosin–Rammler equation and Numerical Example 1