{ "metadata": { "name": "", "signature": "sha256:89c91579d721ed0f833b399593203d4eb19421ab1695bc830f1be612e46bd826" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Chapter 10 : Gas Dispersion and Gas Interchange in Bubbling Beds" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\n", "Welcome to pylab, a matplotlib-based Python environment [backend: module://IPython.zmq.pylab.backend_inline].\n", "For more information, type 'help(pylab)'.\n" ] } ], "prompt_number": 2 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Example 1, Page 253\n" ] }, { "cell_type": "code", "collapsed": false, "input": [ "\n", "from numpy import *\n", "%matplotlib inline\n", "#Variable declaration\n", "umf=[0.01,0.045]; #Velocity at minimum fluidization condition in m/s\n", "ephsilonmf=[0.5,0.5]; #Void fraction at minimum fluidization condition\n", "D=[2E-5,7E-5]; #Diffusion coefficient of gas in m**2/s\n", "g=9.81; #Acceleration due to gravity in m/s**2\n", "\n", "#CALCULATION\n", "db=[5.,10.,15.,20.];\n", "n=len(umf);\n", "m=len(db)\n", "Kbc = zeros((n,m))\n", "Kce = zeros((n,m))\n", "Kbe = zeros((n,m))\n", " \n", "for i in range(n):\n", " for j in range(m):\n", " Kbc[i][j]=4.5*(umf[i]/db[j])+5.85*((D[i]**0.5*g**0.25)/db[j]**(5.0/4));#Gas interchange coefficient between bubble and cloud from Eqn.(27)\n", " Kce[i][j]=6.77*((D[i]*ephsilonmf[i]*0.711*(g*db[j])**0.5)/db[j]**3)**0.5;#Gas interchange coefficient between emulsion and cloud from Eqn.(34)\n", " Kbe[i][j]=(Kbc[i][j]*Kce[i][j])/(Kbc[i][j]+Kce[i][j]);#Gas interchange coefficient between bubble and emulsion from Eqn.(14)\n", "\n", "#OUTPUT\n", "i=0;\n", "j=0;\n", "k=0;\n", "while k