{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NOTEBOOK_HEADER-->\n",
    "*This notebook contains course material from [CBE30338](https://jckantor.github.io/CBE30338)\n",
    "by Jeffrey Kantor (jeff at nd.edu); the content is available [on Github](https://github.com/jckantor/CBE30338.git).\n",
    "The text is released under the [CC-BY-NC-ND-4.0 license](https://creativecommons.org/licenses/by-nc-nd/4.0/legalcode),\n",
    "and code is released under the [MIT license](https://opensource.org/licenses/MIT).*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NAVIGATION-->\n",
    "< [Process Variables](http://nbviewer.jupyter.org/github/jckantor/CBE30338/blob/master/notebooks/02.01-Process-Variables.ipynb) | [Contents](toc.ipynb) | [Blending Tank Simulation](http://nbviewer.jupyter.org/github/jckantor/CBE30338/blob/master/notebooks/02.03-Blending-Tank-Simulation.ipynb) ><p><a href=\"https://colab.research.google.com/github/jckantor/CBE30338/blob/master/notebooks/02.02-Gravity-Drained-Tank.ipynb\"><img align=\"left\" src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open in Colab\" title=\"Open in Google Colaboratory\"></a><p><a href=\"https://raw.githubusercontent.com/jckantor/CBE30338/master/notebooks/02.02-Gravity-Drained-Tank.ipynb\"><img align=\"left\" src=\"https://img.shields.io/badge/Github-Download-blue.svg\" alt=\"Download\" title=\"Download Notebook\"></a>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "VzFcLHpl3OPe"
   },
   "source": [
    "# Gravity Drained Tank"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "VzFcLHpl3OPe"
   },
   "source": [
    "## Summary\n",
    "\n",
    "In the example we study the problem of constructing nonlinear process model for the liquid level in a gravity drained tank."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "xSIqWFm93OPg"
   },
   "source": [
    "## Torricelli's law\n",
    "\n",
    "Torricelli's law states the velocity of an incompressible liquid stream exiting a liquid tank at level $h$ below the surface is \n",
    "\n",
    "$$v = \\sqrt{2gh}$$ \n",
    "\n",
    "This is the same velocity as an object dropped from a height $h$. The derivation is straightforward. From Bernoulli's principle,\n",
    "\n",
    "$$\\frac{v^2}{2} + gh + \\frac{P}{\\rho} = \\text{constant}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "zGDdqtSj3OPj"
   },
   "source": [
    "Applying this principle, compare a drop of water just below the surface of the water at distance $h$ above the exit, to a drop of water exiting the tank\n",
    "\n",
    "$$gh + \\frac{P_{atm}}{\\rho} = \\frac{v^2}{2} + \\frac{P_{atm}}{\\rho}$$\n",
    "\n",
    "$$\\implies v^2 = 2gh$$\n",
    "$$\\implies v = \\sqrt{2gh}$$\n",
    "\n",
    "Torricelli's law is an approximation that doesn't account for the effects of fluid viscosity, the specific flow geometry near the exit, or other flow non-idealities. Nevertheless it is a useful first approximation for flow from a tank."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "H8dmoxAJ3OPk"
   },
   "source": [
    "## Mass Balance for Tank with Constant Cross-Sectional Area\n",
    "\n",
    "For a tank with constant cross-sectional area, such as a cylindrical or rectangular tank, the liquid height is described by a differential equation\n",
    "\n",
    "$$A\\frac{dh}{dt} = q_{in}(t) - q_{out}(t)$$\n",
    "\n",
    "where $q_{out}$ is a function of liquid height. Torricelli's law tells the outlet flow from the tank is proportional to square root of the liquid height\n",
    "\n",
    "$$ q_{out}(h) = C_v\\sqrt{h} $$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "rcxaOAy73OPl"
   },
   "source": [
    "Dividing by area we obtain a nonlinear ordinary differential equation \n",
    "\n",
    "$$ \\frac{dh}{dt} = - \\frac{C_V}{A}\\sqrt{h} + \\frac{1}{A}q_{in}(t) $$\n",
    "\n",
    "in our standard form where the LHS derivative appears with a constant coefficient of 1."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "56j29DiG3OPm"
   },
   "source": [
    "## Step-by-Step Approach to Nonlinear Simulation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "C-BiES7y3OPn"
   },
   "source": [
    "### Step 1. Initialize Jupyter and Python\n",
    "\n",
    "The first step in any Python application is to initialize graphics display (if you will be creating graphic output), and import any needed libraries.  \n",
    "\n",
    "The `%matplotlib inline` is a command telling the Jupyter notebook to display any graphics inside this notebook.\n",
    "\n",
    "We then import a standard graphics library (`matplotlib.pyplot`) that will be referred to with the prefix `plt`, a standard numerical library `numpy` that will be referred to with the prefix `np`, and the function `odeint` from the `scipy.integrate` library.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "lAtqLB373OPp"
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from scipy.integrate import odeint"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "Rz2LNFOP3OPu"
   },
   "source": [
    "### Step 2. Define parameters\n",
    "\n",
    "Provide values for all relevant parameters.  Use comments to describe the parameters and units."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "qYEqVEsU3OPv"
   },
   "outputs": [],
   "source": [
    "Cv  = 0.1     # Outlet valve constant [cubic meters/min/meter^1/2]\n",
    "A   = 1.0     # Tank area [meter^2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "QqTpi1bf3OPy"
   },
   "source": [
    "### Step 3. Write Functions for the RHS of the Differential Equations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "nuIhX1VO3OPy"
   },
   "outputs": [],
   "source": [
    "# inlet flow rate in cubic meters/min\n",
    "def qin(t):\n",
    "    return 0.15\n",
    "\n",
    "def deriv(h,t):\n",
    "    return qin(t)/A - Cv*np.sqrt(h)/A"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "n1EJHyQd3OP1"
   },
   "source": [
    "### Step 4. Choose an Initial Condition, a Time Grid, and Integrate the Differential Equation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "HiA2HZcb3OP1"
   },
   "outputs": [],
   "source": [
    "IC = [0.0]\n",
    "t = np.linspace(0,200,101)\n",
    "h = odeint(deriv,IC,t)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "UJ5008ZA3OP6"
   },
   "source": [
    "### Step 5. Visualize and Analyze Results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 294
    },
    "colab_type": "code",
    "id": "DyQyXfpp3OP8",
    "outputId": "be349d5d-4806-4249-e3c5-85f14d31ff28"
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEVCAYAAADpbDJPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8XXWd//HXzd4sbdI23TdK208L\nZS0ge8suiIqy+Rt0VHAYFxz8Oe6i44yoI44wCPMb1FFER1FkEZCOKKAwZa1QKIX2U0rXNG2TtNn3\n3Ht/f5yTNg1Jm4Te3Juc9/Px6CP3nuXed05vzuec7zn3+40lk0lERCR6stIdQERE0kMFQEQkolQA\nREQiSgVARCSiVABERCJKBUBEJKJy0h1ARj8ziwGfBq4B8oBcYD3wNXd/8RC+zzpgqbvvMrO/c/cf\nD3L9zUAMaAWKw4z/5u7Lh5DlO8AWd79jsOv283r/BVS4+zd6TCsFngufFgJTgI3h88fc/bohvtcK\n4HZ3//XQE8tIoAIgw+FbwFnAO919h5llAx8DHjOzBe5efSjexN0XApjZFOALwKAKQOgqd18Rvs45\nwJ1m9kV3v3uQWb48hPceFHevA7p/52XAf3VvA5GBUAGQlDKz8cBngGPcfQeAu8eBH5rZr9y9MVzu\nL8DTwPsJzhTeBO4C5gD5wG3ufrOZ3QSMcfdPh+tNBLYA04A6YCbwFDAjPCP4HyC3+2jYzMqA7cAs\nd685UHZ3f9zMrgZ+Zma/Bj4MvAcYB7zo7l8ws68BHyT4W1oLfNDd68zsZ8AGd78xPLP4Tvh7zQR+\n5e7/GOZ5L3AjUARsAP7G3WvMbAJwNzAfeB1oASoGs+3D1/974P+G+SqAD7n7NjP7GHAu0AacCnQC\nl7r7ul7rvwu4GTjV3XcP9v0ls+kagKTaycBWd3+j94zunX8PS4Aj3f0Z4AZgU3hEew7wHTObCdwL\nvLvHOu8GHnf3+h7Trg7fcyHwC+ByM+s+2LkYeOpgO/8engBKgQXh8/OBj4c7/yXAdcCJBDvq/PB5\nX84ETgl/x0+b2Qwzmxvm+z/uPhf4M9DdZPRFoNrdDwM+BVwwwLx7mdk04N+Bs919HrAV+GqPRS4G\n/t3dFwArgOt7rb8IuB14t3b+o5MKgKRaGbC3icfMSs1sXfivwsy+0GPZ5e6eCB//A8F1A9x9I7AT\nOMzdXwBiZnZMuNz7gHv6e3N3f4ngzOCcHsv/ZqDhwzxNBEf9AOu7i1l4/WKmuzeEyz0DzO3npX7l\n7nF3rwR2EZwJvBP4i7uvCZe5A3hP2ER2Zvfv5e6bgScHmrlH9kpgXPgT4H975XvV3V8OH78EzOox\nr4yg2F7t7usH+94yMqgJSFKtmqB5BnhLu/V/EVy87Lanx+MTCY76ZwFxYCr7DljuI9hRbgBOB646\nSIa7gb8xs6eAZcDVZnYS8PNw/gP9tdmb2RhgElAFHNEzo5kVAreE7e8A44FH+snQ8wwlDmQTnFmc\nGTZV9VxuQvhaPdepPfCv2Gf2HOBGM3s3wbYbC7x2kEzdvkVwQXzHYN9XRg4VAEm1Z4FJZnacu68a\nxHr/DdwC3OHuSTPb3mPevcCtBDuzJ/toSurtbuB5gusBT4dF6AXCQnQQlxK05W82s97zPkPQ9LPE\n3ZvM7FvA9AG8ZrdKgrt1Lus9w8xq2XfWAVDOvjt8BupvgAuB0919t5l9guD3GYibCYrOXWZ2anjd\nRkYZNQFJSoU7528CvzCzeQBmlmVmHwCuILjw2ZdJBBdak2b2YYKLpMXhvGeBycBH6Lv5pxMo7m73\nd3cnuKj8rwyi+Sc8sr8J+NwBMq4Ld/6zgYt6ZByIR4EzwmsBmNlJZnZrOO9ZguYqzOxwgjOdwZpE\ncB1ld3ix/PJB5NsA/D+C5q8vDeG9ZQRQAZCUc/ebCC5G3hs2d2wEPgpc5u7/3c9qXwMeMLPVBDut\nHwI/NrPD3T0J/I7gLpaH+1h3NUFTzc6wCQmCs4DJwIMHifvL8PrEdoKCcY27P9TPsncAS83Mge8D\nnwXOMbPPHOQ9AAjvivq78PdcS3DBtbtAfQeYbWabgNuA+wfymr1/F2Bq2FT2S+ArwNzwTqqB5EsS\nXFD/vz2uucgoEtN4ABIFZnYFQcG5It1ZRDKFzgBk1Asv1n4R+EG6s4hkEhUAGdXM7GJgHfBw9zd8\nRSSgJiARkYjSGYCISESNmO8BVFc3DvlUpayskNralkMZ55DI1FyQudmUa3CUa3BGY67y8pJYf/Mi\ncQaQk5N98IXSIFNzQeZmU67BUa7BiVquSBQAERF5KxUAEZGIUgEQEYkoFQARkYhSARARiSgVABGR\niFIBEBGJqBHzRTARiY5EMkk8niSeSJBIJIknkvv/TO7/PJkkeBxOT3Q/TiZJJILXS4bTkkn2Lrf3\ncTi9qCifhsa2cF4wLZlMktj7c9+0ZBKS9HzcPe+t0xLhwnuXh+B5+DiZ7DWdfcvHgPedPZ/xhbmH\nfDurALwNO3ZUcsMNX+QnP/nFftNvvfX7XH75B5g2beCDQy1f/jBFRcUsXXrWftPf9a5zeOSRx/eb\ndt111/LZz36BuXPnDT18L/39LjL6JZNJOrsSdHQl6OiM7/ezs/tnV4LOePizx/OuHj+74sHjeDy5\n92dXPJjelUgSjyfoCqcRi9HREacrkQh39D129vFwRyh7zZg6lguWzDjkr6sCkALXX/+Pg17noove\nnYIkMholk0naOuK0tnfR0t5Fa3sXre1x2jr2f9zeGYesLOob2mjriNPeGd/7s70jTkdX8LizM5Hy\nHW4sBjnZWeRkx8jOyiIvN4sYMfJzs8nOD6ZlZ8XIzo4FP7v/ZWeRFds3PSsW/uz5Lxb+y4LsrCyy\nstg7LZYVIytGr+VixGLsewzEwvXHjSuksbF13/rhcrHux+G6McKfse712bvM3p/sW5+eywD0tW73\n9HDhWPCUrKwYC+ZOpKam6ZD/v6gApED3EXpxcQk33PBF8vLyWLjwCNate53bb//Rfkf1N9zwBd7/\n/itYtepFSktLee97L+Wf//kGqqp2sWjREQN+z3g8zk03fYvKyu10dXXxsY99nLFjx3HbbTfzgx/c\nAcBPf/ojSkrGcsIJJ3HLLTcRi8UoLCzkK1/5Rio2gxxE9468saWDxtZOGls6aW7tpCn819zWRXNr\nJ81tnTS3dtHS3klLW7DTH2onvjnZWeTnZpGfl03xmFzGjy0gPyeLvNzs8F8WeTlZ5OZ0P84mLyeL\nnJwscnP2zcvJjpGbk0VudhY52cG87h18TnawfE5WFrk5sXCnvH93NOXlJVRXH2wo5+GXqblisX67\n83lbRk0BuOeJDaxcV9XnvOzsGPH44P9iTlw4iSvOHnozy29/+2suuOBCLr30Sn7xi58NaJ2VK5+j\nq6uLH/7wTl57bQ333juwIWz/9Kc/MGHCRL785a9TV1fH9dd/nLvu+jU1NdU0NjZSUlLCihVP8d3v\n3syNN/4Tn//8V5g5cxb33/9b7r//Hs4//8Ih/56yv86uOHVNHdQ2tlPf3EFdUzv1TR3UN7fT0NxJ\nQ3MHDS0dNLZ00DXAz2VebhZFBbmUFuczdWIRhfk5FBbkMCY/h8L8HAryshmTn8OYvOBxQV42Bfk5\nTJ08lpamNgrycsjPyyI7S/d9yD6jpgBkoi1bNnHOOecBcNxxS3j++WcOus6mTZs46qijATjyyMXk\n5+cP6L3WrFnNK6+sYvXqlwFob2+ns7OT0047k+eff4bFi48hPz+P8vJJvP76a3z3uzcC0NnZOagz\njahLJJPUN3Wwu76NmoZW9jS0s7uhjea2ODtrmtjT2E5Ta+cBXyMvJ4uxRXnMnFRCSWEuYwvzKC7M\npWRMLsXhv6LwX3FBDoUFueTmDG3HXV5eQrX2+dKPUVMArjh7Xr9H6+k6rQsG2wlO3frrza+rq6v3\nWsRi+/5iBzpgT05OLn/7t1dz3nnv3G/60qVncd9991BfX8fSpWcDUFBQwG23/XC/08odOyoH9D5R\n0BVPUFPfxq49LeyqbaW6tpWqulaq61qpqW8LLmL2IS83i/ElBcycVExpcT6lJXnBz+J8xhXlMa44\nj7GFeRTkZafslF5kMEZNAchEs2bN4fXX17Bw4SJWrnxh7/RYLEZbWxutrTmsX++91pnNn/70KACv\nvvoKHR0dA3qvI45YzIoVT3Leee+ktnYP99xzN3//95/iyCOP4vvf/1caGhr4/Oe/AsC8efN57rln\nOOWU03jssUcpLS1j+vRDf4dBpmvviFO5u5nKmmYqdzezc3cLlbtbqK5tDW7b66WoIIcZ5UVMLB3D\nxHEFTBxXwPixBUwYW8CCwybQ0tSmHbuMKCoAb9PWrVu47rpr9z7/5Cf/Ye/jyy67kq9//Uv8+c+P\nMW/egr3TL7nkMq699sOYLcBs0X6vd/LJp/HIIw9x3XXXMm/efMrLJ/X5vt/+9r9QUFAAwJIlJ/Kh\nD32Ul15aycc/fjXxeJyrrw4yxWIxFi8+hjfecKZMmQLA9dd/jptu+ha//OVd5OXl841v3Ehzc/Oh\n2SAZKJlMsru+jS27mthW1UhFdTMVVU1U1bW+ZdmighwOm1bClLJCJo8vZFLZGCaXFVJeOobCgv7/\nXIoL82htbk/lryFyyI2YMYHfzohgmXBlf+PGDdx8803cfvuP9k7LhFz9ydRsA8lV29jOxsp6Nu1o\nZPPOBrbsbKS5bf+mtuIxucwoL2LaxCKmTwx+Tp1YRMmY3CEdxY/k7ZUOyjU4byfXgUYE0xmAjGjx\nRIKtu5p4o6KeDRV1vFnZQG3j/kfik0rHsGjOeGZPLmbW5BJmlBdTWpyn5hqJPBWAYTJ37rz9jv5l\naLriCTZsr8e31rJuSy0btjcEX3gKjS3M5bj5E5k7bSyHTR3L7CklFBUc+q/Qi4wGKgCS8XbuaWHN\nxt28vrkW31ZHa/u+5pxpE4tYMGMc82eUMm/GOCaOK9CRvcgAqQBIxumKJ1i3tZbVG3azeuNuqmr3\nXaydNrGIk4+YzMLZZdjMUsYW5aUxqcjIpgIgGaGto4vVb+7mpfXVvLpxN63tQbNOQV42SxaUs3ju\neI48bDyL5k3KyIt0IiORCoCkTXtnnFc21LBybRWrN+6msyv4gtXEcQWcdtRUjps3kfkzS8nJ1ldZ\nRVJBBUCGVSKRZO2WWp59bScvrq+mvSM40p86oZATbBJLrJyZk4rVji8yDFQAZFhU1bXy9OodPL1m\nB3sagts0J44r4NwlM3jHEZOZUV6c5oQi0aMCICkTTyRYvWE3f161nTWb9gBBm/6Zx0zjtKOmMG/6\nOB3pi6SRCoAccs1tnTz1ciWPv1Sx92h/3vRxLD12GifYJPLz+u4YT0SGV0oLgJndBJwRvs933P3+\nHvPOBb4NxIHl7v7NVGaR1Kuqa+XRF7by9Ks76OhMkJ+bzbLjpnPWcdOZOUlNPCKZJmUFwMzOAha7\n+ylmNgFYBdzfY5EfABcA24Enzew+d389VXkkdSqqmlj+3BZeWFtFIplkwth8zjl9JmccM1XfwhXJ\nYKk8A3gK6O4DuQ4oMrNsd4+b2Vxgj7tvAzCz5cA5gArACFJR1cSDT2/iRa8GYEZ5ERedPJsTF03S\nyFMiI0DKCoC7x4HuPoavIWjm6e60ZQpQ3WPxKuDwA71eWVlhv4OqDER5ecmQ102lTM0F/WfbXt3E\nf//PWla8EgwiY7PKuPK8BZywaPKwXNTN1G2mXIOjXIOTilwpvwhsZu8lKADnH2Cxg+41amtbhpxh\nNHbxmmp9Zatv7uChFZt48uVKEskkc6aUcMkZczlq7nhisRg1NU1pyZUJlGtwlGtw3mZ30P3OS/VF\n4AuArwLvdPf6HrMqCc4Cuk0Pp0kG6uyK88eV2/j9M1to74wzeXwhly2dy/ELynUbp8gIlsqLwOOA\n7wHnuvuenvPcfbOZjTWzOUAFcDFwVaqyyNAkk0lefqOGXz/xBtV1bZQU5nLF2fM44+ip6p5BZBRI\n5RnAlcBE4B4z6572BPCquz8AfAK4O5z+G3dfn8IsMkhVe1q49d7VrH5zN9lZMc4/cSbvOW0Ohbqr\nR2TUSOVF4B8B/Y6A4u5PAaek6v1laOKJBH9aWcGDT2+ivSPOotllXHXeAqZNLEp3NBE5xPRNYNmr\nsqaZnzyylk07GigpzOOD5y3g1MVT1M4vMkqpAAiJRJI/rtzG/U9tpCue4OQjJ3PdFcfR0dqR7mgi\nkkIqABFX29jOjx9+jXVb6xhbmMuHLjiSJVbOuOJ8qlUAREY1FYAIW7W+mp8uX0tzWxfHzZ/IRy5c\nSEmhhlgUiQoVgAjqiie4588beOyvFeTmZPGhC4xlx05TW79IxKgARExtYzv/+eAaNlTUM3VCIZ+4\nZLEGYxGJKBWACHmjoo7/eGANDc0dnLhwEh+9aCEFefoIiESV/voj4ulXd3DXH9aRSMAHzpnPeSfM\nUJOPSMSpAIxyiWSS+5/cyPLntlCYn8Mn37eYI+aMT3csEckAKgCjWGdXnB//fi1/XVfF5LIxXH/5\nMUwZX5juWCKSIVQARqmWti5uv38167bWsWBmKde9/yiKx6gfHxHZRwVgFKpraueWe15hW1UTSxaU\nc+17jiD3bQymIyKjkwrAKLOnoY2b7l5FVW0ry46bzgfPW0BWli72ishbqQCMItV1rXzv7lXU1Lfx\nrlNm8/4z5+pOHxHplwrAKLGrtoXv3b2KPQ3tXHLGYbzntMPSHUlEMpwKwChQU9fKTb9aRW1jO5cv\nO5wLT56d7kgiMgJoXL8RrraxnX/79cva+YvIoKkAjGANLR18/zcvU1XXyrtPnaOdv4gMigrACNXW\n0cW/3/MKlTXNnH/iTC45Q23+IjI4KgAjUDyR4I4HX2PzzkZOO2oKV549T3f7iMigqQCMMMlkkl88\nup7Vb+5m8WHj+fA7F2rnLyJDogIwwvz+2S089UolsyYX84lLFpOTrf9CERka7T1GkJfWV/PAUxuZ\nMDafz1x+DGPydReviAydCsAIUVHdxI9//zp5uVl8+tKjKS3OT3ckERnhVABGgKbWTn5w72raO+Jc\n864jmDW5JN2RRGQUUAHIcIlEkjseXENNfRvvPnUOJy6clO5IIjJKqABkuIee3sTrm2s5dt5E3qt7\n/UXkEFIByGCvbdrDw09vZsLYAq65eBFZut1TRA4hFYAMVdvYzo8efo2srBifuGQxRQUazUtEDi0V\ngAyUSCT50UOv0djSyZVnz2PutLHpjiQio5AKQAZ69IWt+LY6jps/kXOWzEh3HBEZpVQAMszWXY3c\n/9RGxhXl8ZEL1c2DiKSOCkAG6eyK8+OHXyeeSPLRixZRUpiX7kgiMoqpAGSQ+57cyPaaZs46fjpH\nHz4h3XFEZJRTAcgQG7bX86eV25g8vpArzpqX7jgiEgEp7U3MzBYDDwK3uPvtveZtBrYB8XDSVe6+\nPZV5MlVnV4I7l68lCXz0woXk52anO5KIREDKCoCZFQG3AY8fYLEL3b0pVRlGikee3cyO3S2cdfx0\nFswsTXccEYmIVDYBtQMXAZUpfI8Rr6K6iUee3UJZST6XLT083XFEJEJSdgbg7l1Al5kdaLE7zGwO\nsAL4srsn+1uwrKyQnJyhN42Ul2deD5qJRJJfPvYG8USST19xLLNmlKU70n4ycZuBcg2Wcg1OlHKl\nc0SRrwN/APYAvwMuBe7tb+Ha2pYhv1F5eQnV1Y1DXj9VXt64B99Sy4kLJzGnvCijMmbqNlOuwVGu\nwRmNuQ5UONJWANz9592PzWw5cBQHKACjTXNbJ3ctf5383GyuPFt3/YjI8EvLbaBmNs7MHjWz7m86\nLQXWpCNLuvzufzdR39TBu0+bw/ixBemOIyIRlMq7gJYA3wfmAJ1mdhnwELDJ3R8Ij/qfM7NWYBUR\nOvrfVtXEEy9VML28iPNOmJnuOCISUam8CPwisOwA828Fbk3V+2eqZDLJL/+0nmQSrr3kaHJz9F08\nEUkP7X2G2ao3ali/rY5j503keA3vKCJppAIwjOKJBPf+5U2yYjEuP0v3/ItIeqkADKP/fWUHO/e0\ncOYxU5k6oSjdcUQk4lQAhklbRxe/W7GJ/Nxs3nu6BncXkfRTARgmj76wjYbmDi44aSbjivPTHUdE\nRAVgODS0dPCH57cytiiPC06ale44IiKACsCwePSFrbR3xnnXKbMZk5/O3jdERPZRAUixhpYOnnhx\nO+OK81h27LR0xxER2UsFIMUefT48+j95NrlvozdTEZFDTQUghRpaOnj8pQpKi/NYqqN/EckwKgAp\n9OjzW+noTPCuU+bo6F9EMo4KQIo0tXbuPfo/85ip6Y4jIvIWKgAp8sRLFXR0JnjnSbN09C8iGWlA\n9ySaWS4w2d0rzOxo4BjgPncf+jBdo1hHZ5zH/lpBYX4OZ6rtX0Qy1EDPAO4CTjaz6cD9BKN3/SxV\noUa6p1/dQVNrJ2cdP52CPN33LyKZaaAFYLq73wtcCfw/d/8CMD51sUauRCLJoy9sIyc7i3OXzEh3\nHBGRfg20AOSbWQx4H/D7cFpxaiKNbC+tr6aqrpVTF09Rnz8iktEGWgD+AtQDO9x9vZl9BvCUpRqh\nkskk//P8FmLABSdpqEcRyWwDLQC3ArPc/Yrw+YPAx1ITaeR6c3sDm3Y0cuz8iervX0Qy3kCvUP7S\n3c/ufuLum1KUZ0R74qUKALX9i8iIMNACsN7Mfg48A3R0T3T3n6Yk1QhU39zBynVVTJ1QyMLZZemO\nIyJyUAMtAPlAHHhHj2lJQAUg9NTL24knkpx9/AxisVi644iIHNSACoC7f9TMsoBJ7r4zxZlGnHgi\nwV9erqQgL5tTF09JdxwRkQEZ0EVgMzsbeJPgbiDM7BYze1cKc40oq9bXUNvYzqmLp2jAFxEZMQZ6\nF9C3gZOBHeHzbwFfS0miEaj74u/Zx+vir4iMHAMtAE3uvqv7ibvX0ONicJTt2N3Muq11LJpdxrSJ\nuvVTREaOgbZXtJrZUiBmZmXAB4C21MUaOVasDk6KNOCLiIw0Ay0AnwT+EzgR2ACsAP4uVaFGiq54\ngqfX7KSoIIfj5k9MdxwRkUEZaAE43N0v7jnBzC4Bthz6SCPHmo17aGju4JzjZ6jPfxEZcQ5YAMxs\nDnA48G9m9lmg+wb3XODfgd+lNF2GW/Fq0Pxz+tEa8UtERp6DnQFMJegCeg7w9R7TE8AdKco0IjQ0\nd/DKhhpmTipm1mR1jCoiI88BC4C7Pws8a2bL3T3SR/u9PfvaTuKJJKcfPVXf/BWREWmgt4G+bGb3\nmtmfAczsY2Y2P4W5MloymWTF6h3kZMc45Uh981dERqaBFoAfAT/vsfz6cFokbdnVyPaaZo6dN5Hi\nMbnpjiMiMiQDLQC57v4QQds/7v5U6iJlvudfD74Tp6N/ERnJBtxxjZmVEvQAipkdCYwZwDqLCQaP\nucXdb+8171yCLibiwHJ3/+YgcqdNIpnkhbVVjMnPYfHcCemOIyIyZAM9A/gX4DlgiZmtBv4EfOVA\nK5hZEXAb8Hg/i/wAuBQ4DTjfzI4YYJa0emNbHbWN7SxZUE5uzkA3n4hI5hnoHsyBu4DvE3wT+OfA\n6QdZpx24CKjsPcPM5gJ73H2buyeA5cA5Aw2dTi+srQLgHUdMTnMSEZG3Z6BNQP8DvAhsB14Lpx3w\n6qe7dwFdZtbX7ClAdY/nVQRfOOtXWVkhOW/j27bl5SVDXrdbVzzBi+urKS3O54wlM8nOfvtnAIci\nV6pkajblGhzlGpwo5RpoAdjt7lcf8nff56A30tfWtgz5xcvLS6iubhzy+t1e3bh7b9cPe/Y0v+3X\nO1S5UiFTsynX4CjX4IzGXAcqHAMtAA+Y2VXAs0BX90R33zqkREGzUM9baKbTR1NRpum+++ekIyal\nOYmIyNs30AJwNHAVsLvHtCQwayhv6u6bzWxs2NdQBXBx+PoZq7Mrzkvrq5kwNp/Dp49LdxwRkbdt\noAXgZKDM3dsH+sJmtoTgovEcoNPMLgMeAja5+wPAJ4C7w8V/4+7rB5w6DdZs2kNbR5xlx00nS10/\niMgoMNACsBIoILizZ0Dc/UVg2QHmPwWcMtDXS7dV62sAWLKgPM1JREQOjYEWgBnAZjNby/7XAM5M\nSaoME08keHlDDeOK8jhs2th0xxEROSQGWgC+ldIUGW5DRT1NrZ0sO3aamn9EZNQYUAFw9ydTHSST\nrXojaP45Ts0/IjKKqC+Dg0gmk7y0vpqCvGwWzipLdxwRkUNGBeAgKqqbqalv4+jDJ6jvHxEZVbRH\nO4hV64MeK46br+YfERldVAAO4qU3qsnOinGUun4WkVFGBeAAaupb2bqriUWzyygsGPDQCSIiI4IK\nwAG8unEPAMfMm5jmJCIih54KwAGs2Rh0fXTU4Wr+EZHRRwWgH13xBK9vqWVy2RgmlR509EsRkRFH\nBaAfGyrqae+Ia9xfERm1VAD68eqmsPln7vg0JxERSQ0VgH68tnEPOdlZ2Ex9+1dERicVgD7UNbWz\ntaoJmzmO/Lyhj0MsIpLJVAD68Nqm4PbPIw9T+7+IjF4qAH14daPa/0Vk9FMB6CWRSPL65lrKSvKZ\nNrEo3XFERFJGBaCXzTsbaWrt5Ki544lp8BcRGcVUAHpZt7UWgCPmqPlHREY3FYBeuguAzSxNcxIR\nkdRSAeghnkjwRkU9UycUMq44P91xRERSSgWghy07m2jviGMa+lFEIkAFoAdX84+IRIgKQA/rttYB\nYLNUAERk9FMBCMUTCdZX1DFlfCGlav8XkQhQAQjta//X0b+IRIMKQMi3he3/KgAiEhEqACHvbv9X\n988iEhEqAITt/9vqmFw2hrIStf+LSDSoAABbdzXRpvv/RSRiVACANyrqAVgwc1yak4iIDB8VAGBj\nZVAADp+uAiAi0aECAGysbKB4TC6TSsekO4qIyLCJfAGob2qnpr6NudPGqv9/EYmUnFS+uJndApwM\nJIHr3X1lj3mbgW1APJx0lbtvT2WevmysbADg8Gljh/utRUTSKmUFwMyWAvPd/RQzWwT8FDil12IX\nuntTqjIMxJthAZir9n8RiZhUNgGdA/wOwN3XAmVmlnGH2Rsr64kBh03JuGgiIimVyiagKcCLPZ5X\nh9Maeky7w8zmACuAL7t7sr8cui6rAAALiklEQVQXKysrJCcne8hhystL3jItnkiyeWcjMyaXMDtN\n3wDuK1emyNRsyjU4yjU4UcqV0msAvfS+wvp14A/AHoIzhUuBe/tbuba2ZchvXF5eQnV141umb6sK\nvgA2e3Jxn/NTrb9cmSBTsynX4CjX4IzGXAcqHKksAJUER/zdpgE7up+4+8+7H5vZcuAoDlAAUuHN\n7vv/dQFYRCIoldcA/ghcBmBmxwOV7t4YPh9nZo+aWV647FJgTQqz9Gnj9u47gHQBWESiJ2VnAO7+\njJm9aGbPAAngU2b2EaDe3R8Ij/qfM7NWYBXDfPQPsHFHA/l52UybWDTcby0iknYpvQbg7l/qNemV\nHvNuBW5N5fsfSEtbJ5U1zSyaXUZWlr4AJiLRE9lvAm/aEVxQmav2fxGJqMgWgM07g/b/Obr/X0Qi\nKrIFYFtV8AXk2ZOL05xERCQ9Il0AxuTnMGFcQbqjiIikRSQLQHtHnJ27W5g5qVg9gIpIZEWyAFTU\nNJEEZk1S84+IRFckC8C2XUH7/0y1/4tIhEWzAIQXgGdNysxOn0REhkMkC8DWqkays2L6BrCIRFrk\nCkAikaSiqpmpEwrJzYncry8islfk9oBVda20d8aZqeYfEYm4yBWAve3/ugAsIhEXuQKwdVfQB9BM\n3QIqIhEXuQLQfQagAiAiURe5ArB1VyNlJfmUFOYdfGERkVEsUgWgoaWDuqYOHf2LiBCxAqALwCIi\n+0SqAGwPC8CMchUAEZFIFYCdta0ATJ2gbwCLiESqAOza0wLApLIxaU4iIpJ+kSoAO/e0MH5sPvm5\n2emOIiKSdpEpAO0dcWob25lcVpjuKCIiGSEyBWBXbdD8M2W8CoCICESqAAQXgCerAIiIABEqADv3\ndJ8B6AKwiAhEqAB03wGkMwARkUCkCkB2VoyJ4wrSHUVEJCNEpgDs3NNCeekYsrMi8yuLiBxQJPaG\nDc0dNLd16Q4gEZEeIlEAKquDPoD0DWARkX0iUQC2hwVAZwAiIvtEqgDoDiARkX0iUQAqq5sBnQGI\niPQUiQKwvbqJvNwsSos1DKSISLdRXwASySSVNc1MKSskFoulO46ISMbISeWLm9ktwMlAErje3Vf2\nmHcu8G0gDix392+mIkNdYzsdnXG1/4uI9JKyMwAzWwrMd/dTgGuAH/Ra5AfApcBpwPlmdkQqcuxU\nFxAiIn1KZRPQOcDvANx9LVBmZmMBzGwusMfdt7l7AlgeLn/IVYW9gKoTOBGR/aWyCWgK8GKP59Xh\ntIbwZ3WPeVXA4Qd6sbKyQnJyBj+S14lHTcMr6ll24mzGFecPev1UKy8vSXeEfmVqNuUaHOUanCjl\nSuk1gF4OdAX2oFdna8MBXQarKCfGDVe/g+rqRqpbO4b0GqlSXl5CdXVjumP0KVOzKdfgKNfgjMZc\nByocqWwCqiQ40u82DdjRz7zp4TQRERkmqSwAfwQuAzCz44FKd28EcPfNwFgzm2NmOcDF4fIiIjJM\nUtYE5O7PmNmLZvYMkAA+ZWYfAerd/QHgE8Dd4eK/cff1qcoiIiJvldJrAO7+pV6TXukx7ynglFS+\nv4iI9G/UfxNYRET6pgIgIhJRKgAiIhGlAiAiElGxZDKZ7gwiIpIGOgMQEYkoFQARkYhSARARiSgV\nABGRiFIBEBGJKBUAEZGIUgEQEYmo4RwQJi0ONDB9GrLcBJxBsN2/A7wHWALsDhf5nrs/MsyZlgG/\nBV4LJ70K3AT8AsgmGMPhQ+7ePsy5rgE+1GPSCcBfgSKgOZz2j+7+Yu91U5hpMfAgcIu7325mM+lj\nO5nZVcBnCHrB/ZG7/yQNue4EcoFO4IPuvtPMOoGne6x6jrvHhzHXz+jj854B2+u3QHk4ezzwHPBt\ngr+F7s9XtbtfnuJcvfcPK0nx52tUF4CeA9Ob2SLgp6SpB1IzOwtYHGaZAKwCngC+7O6/T0emHp50\n98u6n5jZncB/uPtvzezbwNXAfw5noPBD/ZMwz1LgCuBI4KPuvmY4s4QZioDbgMd7TP4Xem0nM/s5\n8HXgJKADWGlmD7j7nmHMdSPBjuEeM/sU8FngCwRdsS9LRY4B5oJen/dwubRur547djP7KfBf+2YN\n2/bqa//wOCn+fI32JqB+B6ZPg6eA7g9aHcGR7OAHOR4ey4CHwscPA+emLwoQfOC/meYM7cBF7D9y\n3TLeup3eAax093p3byU44j5tmHN9ErgvfFwNTEjh+/enr1x9yYTtBYCZGVDq7i+k8P3709f+YRkp\n/nyN6jMADjww/bAKT7W7my6uAZYDceA6M/ssUAVc5+41w50NOMLMHiI4/f1noKhHk08VMDUNmQAw\nsxOBbWETBsC/mNlEYC3wmfCPIOXcvQvoCjN062s7TSH4nNFr+rDlcvdmADPLBj5FcKYCUGBmvwJm\nA/e5+83DmSu03+edDNhePVxPcHbQbYqZ3UswnO1/uPsvU5irr/3DBan+fI32M4DeDjr4fKqZ2XsJ\n/oOvI2jf+5K7nw28DHwjDZHeINjpvxf4MEGzS88Dg3Rvs48BPwsf3wp83t3PJBxlLl2h+tDfdkrL\n9gt3/r8AnnD37uaOzwHXAucDV5nZCcMcayCf93RtrzzgdHf/czhpN/A14P8QXKv7ppml/ECo1/6h\np5R8vkb7GcCBBqYfdmZ2AfBV4J3uXs/+7aMPMczt7ADuvh34Tfj0TTPbCZxoZmPCo+vpHPw0PpWW\nAZ8GCIcS7fYwcGU6AvXQ1Md26v2Zm05wUXG43Qm84e7/3D3B3e/ofmxmjwNHEVxYHxY9ChHs+7zf\nS2Zsr6XA3qafcPzyO8OnNWb2V2AhKdx/9N4/mFnKP1+j/Qyg34Hph5uZjQO+B1zcfcHGzO4zs7nh\nIsuAdFzcvMrMPhc+ngJMJvjgXxoucinwh+HOFeaZBjS5e4eZxczsMTMrDWcvIw3bq5fHeOt2ep6g\ngJaaWTFB++z/Dmeo8C6RDnf/px7TzMx+FW7HnDDXa/2+SGpy9fV5T/v2Cp1IjyFrzewsM7s5fFwE\nHAukbNzyvvYPDMPna9R3B21m/wrsbTJw91cOskqqclxLcMrb80N0J8GpXgvQRHCHS9Uw5yoBfgWU\nAnkEzUGrgJ8DBcCWMFfncOYKsy0BbnT3C8PnVwBfJGgr3Q5c4+4tw5jl+8AcglsrtwNXETRP7bed\nzOwy4PMEtx7flsq2435yTQLa2Het63V3/6SZfRc4m+Bv4SF3/9Yw57oN+BK9Pu8ZsL3eT/C5X+Hu\nvwmXyyG4G8gIbtb4T3e/s6/XPES5+to/fDjMkLLP16gvACIi0rfR3gQkIiL9UAEQEYkoFQARkYhS\nARARiSgVABGRiFIBEDkAM/tL+K3aoaz7wUOdR+RQ0m2gIikQFo217r4g3VlE+qMCIJEVjoXwVaCC\n4JugzwGrgfcBE4ELgW0E/erfQNCr5gxgPvBnd/+0mX0EONfdPxi+5l8IumP+EPABgq62zw+/xPZp\ngr5bqgn6OKpn35eNksAqd8+k/o1klFMTkETdScA/Egw4cxVQ5+5nEfQie1mvZY8Lp50IfNTMyg7w\nuv9EMIjI+eEALV8lKBSnA38BvkLQF8873P0Udz8VeDnsEkBkWIz2zuBEDmZtj76ZdgPPhNMrgN47\n4xVht72tZlZD0H32QJxC0GXvo2E3xPnAJoIurWvMbDlB53b3hJ0EigwLFQCJuq4DPO/d1W7vZWME\nTTc95fXxHu3AC+5+cR/zzgg7KryYYHSn09w9bT3WSrSoAIi8PQ3ATAAzm0QwbCUEHa7lho9XAj82\nsynhwDaXEwzntx040t3vAl4ys6OABaSxy3KJFhUAkbfnj8DnzOw5giad7iakSmCnmb1I0Bvt9cDv\nzayFoDfMDxMUgX8ys78n6L3zTfYftF0kpXQXkIhIROkuIBGRiFIBEBGJKBUAEZGIUgEQEYkoFQAR\nkYhSARARiSgVABGRiPr/YZwA5Y19DxcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "tags": []
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(t,h)\n",
    "\n",
    "plt.xlabel('minutes')\n",
    "plt.ylabel('meters')\n",
    "plt.title('Gravity-Drained Tank')\n",
    "plt.legend(['Liquid Level']);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "taamgWsu3OQB"
   },
   "source": [
    "### Step 6.  Adjust and Repeat as Needed\n",
    "\n",
    "Rarely does a simulation provide you with the information you seek on the first try.  You will almost certainly need to revise parameter values, choose different time grids, and add visualization features to obtain the most informative results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "fuPuM5fV3OQC"
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NAVIGATION-->\n",
    "< [Process Variables](http://nbviewer.jupyter.org/github/jckantor/CBE30338/blob/master/notebooks/02.01-Process-Variables.ipynb) | [Contents](toc.ipynb) | [Blending Tank Simulation](http://nbviewer.jupyter.org/github/jckantor/CBE30338/blob/master/notebooks/02.03-Blending-Tank-Simulation.ipynb) ><p><a href=\"https://colab.research.google.com/github/jckantor/CBE30338/blob/master/notebooks/02.02-Gravity-Drained-Tank.ipynb\"><img align=\"left\" src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open in Colab\" title=\"Open in Google Colaboratory\"></a><p><a href=\"https://raw.githubusercontent.com/jckantor/CBE30338/master/notebooks/02.02-Gravity-Drained-Tank.ipynb\"><img align=\"left\" src=\"https://img.shields.io/badge/Github-Download-blue.svg\" alt=\"Download\" title=\"Download Notebook\"></a>"
   ]
  }
 ],
 "metadata": {
  "colab": {
   "name": "Gravity Drained Tank.ipynb",
   "provenance": [],
   "version": "0.3.2"
  },
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}