![\"Open](\"https://colab.research.google.com/assets/colab-badge.svg\" "\\"Open"){kind=icon}
![\"Download\"](\"https://img.shields.io/badge/Github-Download-blue.svg\" "\\"Download"){kind=icon}
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Model Library"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 185,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from scipy.integrate import odeint,ode\n",
"\n",
"class GravityTank():\n",
" \n",
" def __init__(self, A=1, Cv=1, h = 0):\n",
" self.A = A\n",
" self.Cv = Cv\n",
" self.Qin = 0\n",
" self._h = h\n",
" self.IC = self._h\n",
" \n",
" def deriv(self,h,t):\n",
" h = max(0,h)\n",
" self._h = h\n",
" return (self.Qin/self.A - self.Cv*np.sqrt(h)/self.A)\n",
" \n",
" @property\n",
" def Qin(self):\n",
" return self._Qin\n",
" \n",
" @Qin.setter\n",
" def Qin(self,Qin):\n",
" self._Qin = Qin\n",
" \n",
" @property\n",
" def Qout(self):\n",
" return self.Cv*np.sqrt(self._h)\n",
" \n",
" "
]
},
{
"cell_type": "code",
"execution_count": 186,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": [
"["
]
}
],
"metadata": {
"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.6.8"
}
},
"nbformat": 4,
"nbformat_minor": 2
}