{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Chapter 13 : Non newtonian fluid flow in circular pipes" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Example 13.1 page no : 432" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The pressure gradient is 61.241859 Pa/m\n" ] } ], "source": [ "#Calculate the pressure gradient\n", "import math\n", "\n", "# variables\n", "v=1. #ft/s\n", "d=0.5 #ft\n", "\n", "# calculation\n", "A=(math.pi)/4*d**2 #ft**2\n", "Q=v*A #ft**3/s\n", "#Let DP denote the pressure gradient\n", "n=0.41 #dimentionless\n", "K=0.66 #kg/m/s\n", "#1 m = 3.281 ft \n", "Q1=Q/3.281**3 #m**3/s\n", "d1=d/3.281 #m\n", "DP=(Q1*8*(3*n+1)/(n*(math.pi)*d1**3))**n*(4*K/d1) #Pa/m\n", "\n", "# result\n", "print \"The pressure gradient is %f Pa/m\"%DP" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Example 13.3 page no : 437" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The fanning friction factor is 0.050081\n", "The reynolds number is 318.531771\n", "The flow is Laminar\n" ] } ], "source": [ "#Calculate the fanning friction factor and reynolds number by power law\n", "\n", "# variables\n", "DP=61.3 #Pa/m (pressure gradient)\n", "D=0.152 #m\n", "V_avg=0.305 #m/s\n", "rho=1000. #kg/m**3\n", "\n", "# calculation\n", "f=DP*D/(4*rho*V_avg**2/2.) #dimentionless\n", "print \"The fanning friction factor is %f\"%f\n", "n=0.41 #dimentionless\n", "K=0.66 #dimentionless\n", "R_pl=8*rho*V_avg**(2-n)*D**n/(K*(2*(3*n+1)/n)**n) #dimentionless\n", "\n", "# result\n", "print \"The reynolds number is %f\"%R_pl\n", "if (R_pl<2000):\n", " print \"The flow is Laminar\"\n", "else:\n", " print \"The flow is turbulent\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Example 13.4 page no : 438" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The pressure gradient is 0.486400 KPa/m\n" ] } ], "source": [ "#Calculate the pressure gradient\n", "\n", "# variables\n", "D=0.152 #m\n", "V_avg=3.04 #m/s\n", "rho=1000. #kg/m**3\n", "n=0.41 #dimentionless\n", "K=0.66 #dimentionless\n", "\n", "# calculation\n", "R_pl=8*rho*V_avg**(2-n)*D**n/(K*(2*(3*n+1)/n)**n) #dimentionless\n", "#print \"The reynolds number is %f\"%R_pl\n", "f=0.004 #dimentionless (fanning friction factor)\n", "#Let DP denote the pressure gradient\n", "DP=4*f*(rho/D)*(V_avg**2/2)/1000 #KPa/m\n", "\n", "# result\n", "print \"The pressure gradient is %f KPa/m\"%DP" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Example 13.5 page no : 440" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The headstrom number is 52427.752042\n", "The reynolds number is 15942.778426\n", "The fanning friction factor is 0.006189\n" ] } ], "source": [ "#Calculate the headstrom ,reynold numbers and the fanning friction factor\n", "\n", "# variables\n", "tow_yield=3.8 #Pa\n", "mew=0.00686 #Pa.s\n", "D=0.0206 #m\n", "rho=1530.0 #kg/m**3\n", "V=3.47 #m/s\n", "\n", "# calculation and Result\n", "He=tow_yield*D**2*rho/mew**2 #dimentionless (headstrom number)\n", "print \"The headstrom number is %f\"%He\n", "R=D*V*rho/mew #dimentionless (reynolds number)\n", "print \"The reynolds number is %f\"%R\n", "dP=11069. #Pa/m\n", "f=dP*D/(4*rho*V**2/2) #dimentionless (fanning friction factor)\n", "print \"The fanning friction factor is %f\"%f" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.6" } }, "nbformat": 4, "nbformat_minor": 0 }