summaryrefslogtreecommitdiff
path: root/Fluidization_Engineering_by_K_Daizo_And_O_Levenspiel/ch5.ipynb
blob: 791c6af29d83649166d5c683541db1c51e708521 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
{
 "metadata": {
  "name": "",
  "signature": "sha256:502044f71f0f4f73ada9ff3f4a68fa509018810bcbdaf16eaa17d381758e9755"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 5 : Bubbles in Dense Beds"
     ]
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 1, Page 126\n"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "dt=60.0;        #ID of tube in cm \n",
      "dp=300;         #Size of particles of bed in micrometers\n",
      "umf=3;          #Velocity at minimum fluidization condition in cm/s\n",
      "ephsilonmf=0.5; #Void fraction of bed at minimum fluidization condition\n",
      "db=5;           #Diameter of bubble in cm\n",
      "g=980;          #Acceleration due to gravity in cm/s**2\n",
      "\n",
      "#CALCULATION\n",
      "#Computation of rise velocity of bubble\n",
      "if (db/dt)<0.125:\n",
      "    ubr=(0.711*((g*db)**0.5));#Rise velocity by Eqn.(3)\n",
      "elif (db/dt)<0.6:\n",
      "    ubr=(0.711*((g*db)**0.5))*1.2*exp(-1.49*(db/dt));#Rise velocity by Eqn.(4)   \n",
      "\n",
      "#Computation of cloud thickness\n",
      "Rb=db/2.0;          #Radius of bubble\n",
      "uf=umf/ephsilonmf;  #Velocity of emulsion gas\n",
      "Rc=Rb*((ubr+(2*uf))/(ubr-uf))**(1/3.0);#Radius of cloud by Eqn.(6)\n",
      "\n",
      "#OUTPUT\n",
      "print 'The rise velocity of the bubble=%.1f cm/s'%ubr\n",
      "print 'The cloud thickness=%.2f cm'%(Rc-Rb)\n",
      "print 'From Fig.8(page 124)comparing fw vs dp, for dp = %.0f micrometer, wake fraction = 0.24'%dp\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The rise velocity of the bubble=49.8 cm/s\n",
        "The cloud thickness=0.30 cm\n",
        "From Fig.8(page 124)comparing fw vs dp, for dp = 300 micrometer, wake fraction = 0.24\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 2, Page 132\n"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "uo=15;     #Superificial gas velocity in cm/s\n",
      "umf=1;     #Velocity at minimum fluidization condition in cm/s\n",
      "lor=2.0;   #Pitch of perforated plate in cm\n",
      "g=980;     #Acceleration due to gravity in cm/s**2\n",
      "#CALCULATION\n",
      "#Case(a) For porous plate\n",
      "dbo1=(2.78/g)*(uo-umf)**2;#Initial bubble size using Eqn.(19)\n",
      "\n",
      "#Case(b) For Perforated plate\n",
      "Nor=(2/math.sqrt(3))*(1/lor)**2;#Number of orifices in cm**-2\n",
      "dbo2=(1.30/(g**0.2))*((uo-umf)/Nor)**0.4;#Initial bubble size using Eqn.(15) assuming inital bubble size is smaller than hole spacing\n",
      "\n",
      "#OUTPUT\n",
      "print 'Case(a) For porous plate'\n",
      "print '\\tInitial bubble size=%.2fcm'%dbo1\n",
      "print 'Case(b) For Perforated plate'\n",
      "print '\\tInitial bubble size=%.2fcm'%dbo2\n",
      "print '\\tSince %f<%f, the equation used is correct.'%(dbo2,lor)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Case(a) For porous plate\n",
        "\tInitial bubble size=0.56cm\n",
        "Case(b) For Perforated plate\n",
        "\tInitial bubble size=1.55cm\n",
        "\tSince 1.548765<2.000000, the equation used is correct.\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}