{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 241 暖房向け自然冷媒ヒートポンプサイクル\n", "\n", "## 読み込み" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "\n", "import CoolProp.CoolProp as CP\n", "\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "import warnings\n", "warnings.filterwarnings(\"ignore\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## パラメータ定義\n", "### 汎用" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "K = 273.16 # converter from degC to K" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 外生条件" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "params = {\n", "'WF':'CarbonDioxide', # working fluid\n", "'T_bldg': 20 +K, #freezer air temperature degC\n", "'T_amb': -5 +K, # Outdoor air temperature degC\n", "\n", "'DT_cond':5, # condenser approach temperature difference K\n", "'DT_sc':3, # degree of subcooling K\n", "'DT_evap':4, # evaporator approach temprature difference K\n", "'DT_sh':2, # degree of superheat K\n", "\n", "'V_dot_cond_f':2.0, # condenser fan volumetric flow rate m3/s\n", "'V_dot_evap_f':2.5, # evaporator fan volumetric flow rate m3/s\n", "\n", "'W_dot_cond_f':200, # condenser fan power W\n", "'W_dot_evap_f':300, # evaporator fan power W\n", "\n", "'ETA_cond':0.72, # condenser effectiveness\n", "'ETA_evap':0.81, # evaporator effectiveness\n", "\n", "'c':0.025, # clearance volume rate\n", "'V_dot_disp':2*1e-3, # compressor displacement m^3/sec\n", "}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## データフレーム用意\n", "\n", "### Nomenclature\n", "\n", "| Symbol | Description | Unit |\n", "| :--- | :---: | ---: |\n", "| $T$ | Temperature | $K$ |\n", "| $P$ | Pressure | $Pa$ |\n", "| $s$ | Specific entropy | $J/kg-K$ |\n", "| $s_s$ | Specific entropy leaving reversible process | $J/kg-K$ |\n", "| $v$ | Specific volume | $m^3/kg$ |\n", "| $h$ | Specific enthalpy | $J/kg$ |\n", "| $h_s$ | Specific enthalpy leaving reversible process | $J/kg$ |" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "states = pd.DataFrame(columns=['T', 'P', 's', 's_s', 'v', 'h', 'h_s'], \n", " index=list(range(1,4+1)))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## State 4\n", "\n", "\\begin{align}\n", "State4 =& f(T_{cond}, x_4=0) \\\\\n", "\\rm{where} & \\\\\n", "T_4 =& T_{bldg} + \\Delta T_{cond}, \\\\\n", "T_{cond} =& T_4 + \\Delta T_{sc}.\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIkAAAArCAYAAABM8KssAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAABNhJREFUeJztmVuIVlUUx3//GcmokQzTUkhIssZSc8ZeFNQBtR70oTHKCAQfEmXGC9JDqET3FKKLGWQFMTJPRkHlJZLRHmrGFG28jKYVVOYg3iYLKSes1cPeA4fTzOzvm/o+Q9YPNpyz91p77b2//9lrn/PJzHCc/qi40gNw/v+4SJwkLhIniYvESeIicZK4SJwkLhIniYvESVIykUi6T5IlyoJSxR8IkholHZL0ayy7Jc3pxa5B0veSLknaL2larn26pI8ldcZ5LiwgdkE+kkZK2iTpbIx/VNKMfvptkrQ1VzdX0m+SXkiNC0q7k3wOjMyUn4CXc3WbSxh/IJwEngBqgXuBXcCHkib2GEiaD6wHXgRqgDbgE0mjM/1UAR3ACuD3AmMnfSQNBVoBAXOAccAy4EyBMYgP5gfAKjNbU5CTmZW8ADcAfwH15Yj3H4+9C1icud8DvJOz+RZY24f/RWBhkTF79SEIs7XIvpqArfF6BdANLCimj3KdSWoJ6t9fpngASFot6WKiTOvDt1LSI4QnvC3WXQNMBnbkzHcAU0s5l8gDwB5JmyWdkXRA0lJJSjlKeg5YC8wzs+Zigg4a4GCLZTJw3sxOlCOYpE+Bx4GNwHsJ886c7wRgN3At4YmuN7PDsfkmoBI4nevjNDDrXw67EMYADcCrwDpgErAhtr3Rj99sQnqaa2bbig1aLpHUAl+VKRZANXDMzC4T0kUxHCcs/lDgQWCTpDoz68jY5P86Vy91paAC2Gdmq+J9u6SxQCP9i6SDkPKfktRqZheKDVoOaukl1Ui6VdJHktolHZF0m6S7JbVJ6pC0RdJ10XaLpGfjG8dJSbWZfu6S9IWkg5JWAl1mdnkg6cbM/jCz78ys58c4AKyMzeeAP4FbclMZwT93l1JwCjiaq/saGN2Lbd5vBkEoLZJuLCZoyUUiqQoYS24nifl9O/CamdUQcnoXIT00mtl44BtgYXQZD3Sa2RTCAa4+9jMIeBdYYmb3AFOAg9FnI2FX6K/sS0yhAhgMQUAEsc/O2cwmnltKTCtwZ67uDuDHlKOZdQJ1wPXATknDCg1ajnRTQ1jofLqpB/aa2WcAZvZLPCi2mFl7tDkGjJI0BKgws7difSVwPl7PI2zBPengOPBz7LOLItKNpHXANsLr+hDgUcLCZr+VvAI0S9pL+NGWAKMIguzppwq4Pd5WAKMlTSLscCeizVJgqZlVF+pDOIu0SVpD+HxQAywHVhcyPzM7JakO2AnskjTTzM4V4ljqV8gVwAVAufrnybxaxrpngMcy968DDxF2h/cz9W8DszL9LMq0bQdmDnCsTYSnspvw7aEFuL8Xuwbgh2i3H5iea68jnFHypSlj83RY/sJ9ot0cwk55ibDTLs+vbV+vwJm64bGPw8CI5LqUWiT9DH4ZsCFeVwDDgMXAS7FuIiH/DgYWAU9mfL8Ebo7XK4H18Xom4XvM8Cs1r6uxXMn/bpqAMZKOEFLROKAZqJZ0GHgTeNjMuoEJwCGA+E1ghJn1HBSbgamS2gnnl5NmdracE7naUXwCHadP/F9gJ4mLxEniInGSuEicJC4SJ4mLxEniInGSuEicJC4SJ4mLxEniInGSuEicJC4SJ4mLxEniInGSuEicJC4SJ4mLxEniInGSuEicJC4SJ8nfNrUc2onj5bsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "states.loc[4, 'T'] = params['T_bldg'] + params['DT_cond']\n", "T_cond = states.loc[4, 'T'] + params['DT_sc']\n", "states.loc[4, 'P'] = CP.PropsSI('P', 'T', T_cond, 'Q', 0, params['WF']) # Satuation pressure\n", "states.loc[4, 's'] = CP.PropsSI('S', 'T', states.loc[4, 'T'],\n", " 'P', states.loc[4, 'P'], params['WF'])\n", "states.loc[4, 'h'] = CP.PropsSI('H', 'T', states.loc[4, 'T'], \n", " 'P', states.loc[4, 'P'], params['WF'])\n", "ax = plt.axes([0,0,0,0])\n", "ax.axis('off')\n", "plt.text(0, 0, \"$T_{cond}$=\"+'{:.02f} K'.format(T_cond), size=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## State 2\n", "\n", "\\begin{align}\n", "State2 =& f(T_{evap}, x_2=1) \\\\\n", "\\rm{where} & \\\\\n", "T_2 =& T_{amb} - \\Delta T_{evap}, \\\\\n", "T_{evap} =& T_2 - \\Delta T_{sh}.\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIkAAAArCAYAAABM8KssAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAABDRJREFUeJztmV2IVVUUx3//MsqaSu3LgZqHPogJlQyKCEwxsocp+gAJE9FAMGhqICJRQYQpfYisHvIpaEoY0Aj6GIiGqYdGA50sZXQmpCajLDLRErGxr93D2SOb4z1n33Ps3qu39YMN+6yz1l7r7PnvffY9I+cchpHHeY0uwDj7MZEYUUwkRhQTiRHFRGJEMZEYUUwkRpSaiUTSAkku0pbUKn8ZJK2SNCTpmKRfJH0gaUYFv1ZJb3qfcUkjkuYWHSc15t2S3pd00M/Nsgy/3NwZMT2S+lK2+yWdkPRC7qRQ251kEGgN2vfASynblhrmL8M8YBNwFzAf+AsYkDRtwkHSFGA7IKADaAeeAg4VGacCLcBeoAv4vZJDlbmj+MX5DrDKObcmGuCcq3kDLgf+AR6uR77/sO4W4G/ggcC2Hth+puNE/I8DyyrYC+f2cT1An+93ASeBJdXG1+tMchuJ+nfVKd8pJK2WdDzS5mSEX0qy2x4NbA8BOyRtkXRI0m5JnZKUU0alccpQJvcpJHUDG4BHnHObq85apxX5LHC4QbvBNODGSJucEbsV+BI4P7CN+7YBmA08TrLyO3NqOG2ckjtJ4dzBTnIScEBH4Tms0x+qF+hvhEjOoOaNwI/A9Sn7H8BnKdt6YLTIOCVFUih3SiS7gK+BncCUInNRz9fNaa8aSTdJ6pP0uaRBSdMl3Snpw8DnXklbff8xSTskDUv6SNLF3v6u34KHJI1JuiOIL/y6kfQysAiY75wbS5X9EzCSso0CbRWeL2+cMlSdOyN2Lsn5cEDS1Kqz1mFFThzaFqbsFwIfA23+ejnQDVwGfBP4fQrc4vtXBPZNwIO+Pwas9v1FwNuBX6HXDfAq8DPQnvE8vcBgytYNjKRsueOU3Emqyp2xk0wcXFtJhPVFOJ+58XUQyRySd+ENKfujJNvwbt9GgLX+3nfARcACoDeIeQ4YAvYAR4D7vAgPApO8zyxgoGStrwHHSH62Tg9aS+BzO/AnsMYLbCHwG/BkwXE6ga9Si+lW304Aa32/rUjumEj89TUkP7f3AFeeDSLpAn4FlLI/DyzOiOkHZpJ8a7nZ25b6h53sr/cD15F8i+gPYpcCG0vW6jLaupRfh5/gcV/H0+HzVTMOsA5wwfW8jJieIrmrEYm3XeXHGQaubqhIcgrv9MXLX88M7r0CvA68FdheBFb4/grgqO8/4SfrAmAqycHs2kY9VzO2STSON4B7gFFJ48AnwDP+3j6SM0d74L8ZeE/SYu877O2zSL4ebgMuAVY6536offn/HyZW8TmLpG0kXw+/bXQtzUoz/Be4DTjQ6CKamXN+JzFqTzPsJEaNMZEYUUwkRhQTiRHFRGJEMZEYUUwkRhQTiRHFRGJEMZEYUUwkRhQTiRHFRGJEMZEYUUwkRhQTiRHFRGJEMZEYUUwkRhQTiRHlX/z+IBfwuRb1AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "states.loc[2, 'T'] = params['T_amb'] - params['DT_evap']\n", "T_evap = states.loc[2, 'T'] - params['DT_sh']\n", "states.loc[2, 'P'] = CP.PropsSI('P', 'T', T_evap, 'Q', 1, params['WF']) # Satuation pressure\n", "states.loc[2, 's'] = CP.PropsSI('S', 'T', states.loc[2, 'T'], \n", " 'P', states.loc[2, 'P'], params['WF'])\n", "states.loc[2, 'h'] = CP.PropsSI('H', 'T', states.loc[2, 'T'], \n", " 'P', states.loc[2, 'P'], params['WF'])\n", "states.loc[2, 'v'] = 1/CP.PropsSI('D', 'T', states.loc[2, 'T'], \n", " 'P', states.loc[2, 'P'], params['WF'])\n", "ax = plt.axes([0,0,0,0])\n", "ax.axis('off')\n", "plt.text(0, 0, \"$T_{evap}$=\"+'{:.02f} K'.format(T_evap), size=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## State 1\n", "\n", "\\begin{align}\n", "State1 =& f(P_1, h_1) \\\\\n", "\\rm{where} \\\\\n", "P_1 =& P_2 \\text{ which means that pressure loss in the evaporator is neglected,} \\\\\n", "h_1 =& h_4 \\text{ which means that the valve is isenthalpic.}\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "states.loc[1, 'P'] = states.loc[2, 'P']\n", "states.loc[1, 'h'] = states.loc[4, 'h']\n", "states.loc[1, 's'] = CP.PropsSI('S', 'P', states.loc[1, 'P'], \n", " 'H', states.loc[1, 'h'], params['WF'])\n", "states.loc[1, 'T'] = CP.PropsSI('T', 'P', states.loc[1, 'P'], \n", " 'H', states.loc[1, 'h'], params['WF'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## State 3\n", "\n", "\\begin{align}\n", "State3 =& f(P_3, h_3, ) \\\\\n", "\\rm{where} \\\\\n", "P_3 =& P_4 \\text{ which means that pressure loss in the condenser is neglected,} \\\\\n", "\\frac{\\dot{W}_{s,c}}{\\dot{m}} =& h_{s,3} -h_2 \\\\\n", "\\frac{\\dot{W}_{c}}{\\dot{m}} =& \\frac{1}{\\eta_c}\\frac{\\dot{W}_{s,c}}{\\dot{m}}\\\\\n", "h_3 =& h_2+\\frac{\\dot{W}_{c}}{\\dot{m}}\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "states.loc[3, 'P'] = states.loc[4, 'P']\n", "eta_c = 0.79 - 0.0074*(T_cond - T_evap)\n", "states.loc[3, 's_s'] = states.loc[2, 's'] # specific entropy leaving isentropic compressor\n", "states.loc[3, 'h_s'] = CP.PropsSI('H', 'P', states.loc[3, 'P'], \n", " 'S', states.loc[3, 's_s'], params['WF']) # specific enthalpy leaving isentropic compressor\n", "W_dot_s_c_m_dot = states.loc[3, 'h_s'] - states.loc[2, 'h']\n", "W_dot_c_m_dot = W_dot_s_c_m_dot/eta_c\n", "states.loc[3, 'h'] = states.loc[2, 'h'] + W_dot_c_m_dot\n", "states.loc[3, 's'] = CP.PropsSI('S', 'P', states.loc[3, 'P'], \n", " 'H', states.loc[3, 'h'], params['WF'])\n", "states.loc[3, 'T'] = CP.PropsSI('T', 'P', states.loc[3, 'P'], \n", " 'H', states.loc[3, 'h'], params['WF'])\n", "states.loc[3, 'v'] = 1/CP.PropsSI('D', 'P', states.loc[3, 'P'], \n", " 'H', states.loc[3, 'h'], params['WF'])" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAALkAAABICAYAAACqY1TgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAACpxJREFUeJztnX+wVVUVxz9fQLFEESTMSV5CIGoy8CxNK3+WY6WOOlZok2YZmIiWP8bSzHBynEYlc2QGREugmZRKEzWdMJFKUFSEFPyRGDJqiuZvfknQ6o+1r5x33v19z73vvjP7M3PmzVl3n73OOXudfdbe+3vvk5kRieSZPj19ApFIs4lBHsk9McgjuScGeST3xCCP5J4Y5JHcE4M8kntikEdyT26CXNIkSaslbZK0VNIhFcpPkWSp7dVG6420H7kIcknjgeuAK4FOYDFwr6SOCoc+C+ye2MZkVG+kjVCzl/UlLQfmAYOBU4D/Adeb2c8y9LEEeMLMJiRszwF/MLOLSxwzBfiqme2XZb2R9qOpPbmk7YB9gG8CjwAHAdOByyV9tMxxl0haV2E7JJTdHvgUMD9VzXzgsxVOcYSkl0M6cqukEYlzaKTeSBvRr8n17wtsD1xgZncASLoJuAzYBeiWAwdmAL+rUPfL4e8QoC+wNvX5WuCLZY5fApwOPAMMBS4FFkv6pJm90UC9kTaj2UE+DngLuCth+0j4WyrAMbM3gTdr9JXOu1TElvRxb5fC0sPAv4BvAb+ot95I+9HsgedY4DEz25qwjQPWmNnbkoZJmidpmaSVkoZDbekK8B9gK5BOf4bSvRcuiZmtA1YCo7KsN9LztKInfzRl6wSWh5z3HuBcM3tA0kBgUyhTdbpiZpslLQWOAn6f+Pwo4LZqT1TSDsDewANZ1htpA8ysaRveG56csj0I/BQYD/wqIz/jgc3Ad/GB7nXAOuDjiTKTgWcS+9cAhwHDgc8AdwPvpo6pWG/c2n9rWk8uaQ9gV2B5wiY8hfk18Bvgqix8mdlcSbvig8fdgRXAV8xsTaLYEGB0Yn8P4JZgfx14GDgoeUyV9UbanKbPk3dxJu2FL8AMB74NDDGzsyX1AQaZz2pEIpnS6hXPTuANM3sBTxf2lLQSeBxPB1qKpI4wiO1ste9I62j2wDNNJyF9MbP3gGNa7D/Nv/HB8Ys9fB6RJtLSdCUSKYekhcAKM5ucZb3NXtafJSm96LJ/UPxNT9m/IWmjpKEZ+Z4h6dos6qrD9yXhGqel7IdKujNICUzS6SWO313SbEmvB/XjU5IOS3y+k6RfSloT7tliSQfU4yt1TEVlZj31huNmSbq7iL3p7dTsnPxtYGDKdiHwPr6sn+Q8YI6Zvdao0zCLcxwuDGspkg4CJgBPFPl4AD5D831gY4njdwEW4Surx+BjlXOA5H25CTgaX50dg+tp/iLpY7X4KkFZZWYD9XajZe3UzPlJfD58ZWK/Aw/w64F7EvZDcXXi6Iz8Hgi8AfQL+3vgS/HjgQXABuAf+OLPp4G/BdsjQEcDfgcCzwNHAguBaWXKrgNOL2K/ElhU5rgPAVuA41P2pcAVtfgqUm4Kni5Ue71V1RvKzgLurtBOXe4Z8AW8ozwz7O8IzAl+1wIX4+sbs8r5bnVP/gN8lfNRuvbkFwB3mdmzGfk9AfiTmW0J++PC37OAy/HFn/74jb8KuAQ4GJ8zP79QSY3yAoCZuAx3QYPnvkTSXEmvSVouaXLo9cAnC/qybXW4wEbg8w34LVBSmdkE0u30AZJOAv4ITDSzG4J5Kr6AdyLekYwFKn6JpdmzKx8EuaSd8ZXDL+EirYJ9JHAscHiGfo8HfpLYHwu8A4w3s7XB733AycDeFubnJf0Vf0UXqFpeIGkCMBI4tcFzHwFMAq4Ffo4/oNeHz6aZ2XuSHgIulbQCF7qdgj+kqxr0XUmZmTXpdgJA0kTgalzvPz/YBgDfAU4zs/uC7QzgpUpOmh3kbwEDJPUFzgSeNLPFko5gW09+Hi7i+nsWDsNDMwL4c8I8Dn9VJoVVHcDtqcbrIKG1sSrVkJJG42nGIWa2uYHTBx8nPWbbvpSxTNIo4GygMJA9FV81fgkXkT2Or97u34hjq16Z2TAl2gk88M8EDjWzhxL2TwDb4Sll4XzXhwe9LK1IV8CX98/Fn05wjchASYPxnuOaDH2eANxvZusTtrH4sn2STuChlG0ssKywU0O6Ukh1VkjaImkL/lqdFPb713D+rwBPpWxP4w8gAGb2vJkdhg8Ch5nZgXgArK7BT0WsuzIzS4q1E/iA/RXgjESKBj4Qhzpkzq3oycGfzA3AnWH/HXwQMRkfQNxeOCC89ieFc1tlZieWqlzSMLx368C/nHEs3hPMTpTZEe8FksE7GBiWsg3DH8YPbFSfrjwJPJay3ww8h/fwtfTui+iqsQHYC+imlwkBsl7SIHy25aIa/FREKWVmxnRppwSr8dmkhcBMSRPNR52rgP/ig9XV4fw+DOyHD/ZLk8VsRpkR9TD8yXsLH0AU7EMT9nMS9kH4rEdhtL1Lmbq3x4PrCNs2s7FbuBG7JcodjL/SB6RG7ZuA7RK244D3gD4ZXftCUrMreM87Lmwb8G9IjSMxowMcEK7hx3iO/zW8Uzg7UeZo4Mu4BugofBV5Sep6qvFVjzKzYr2VZlfwMVmXdkrfM7xjehG4kW2LljOAF0L77Qv8Ntybm8v6bnKQDwjBvBbYIWEfGezvADsm7DvhPdZUYEyFurtJdfGByeKU7XvJhgy2C/G8N2m7FHgww2svFuSHh+tOb7NS5Y7BH/ZNwD/xVE+Jz7+O917v46/2acDAWn3hU4aW2L8Vlzpsxt9QtwH71nMNRe7HHOCOUu1U7J4lAn0mnq4MwNWr60NM/Qi4H5jeY0Fe5oL7hUDvX+SznYHT8FzwpDJ1XEGYP03Y5gEX9cQ1xa1im88vBGNW7YRPA7+Kf4e4ZLke+d0VM9tiZqvM7P2kXdIoM3vXzObgT+gOic/uT63orcXzMST1keu+FwG3KKoL2wZJQyQdj6dB9wXzInw2qNa6OuXyj5GhbWfjb/+5ZY8LT0RbIGk2/rMVG/DZkMlmtjWMsl/A57Q3hrI74a/XPfH8brKZPRg+6xfsL6YfpEhrkbQAn52ZBVxmDQRcCOwb8YH5FnwscqGZLS17XDsFeSkk7QNMMLPzKxaOtBVqkrKwFnqFCtHMni4EeCtUa6VQVBcWq7vH1IXV0qtUiCFtierCqC6sjSaPqDNVIRLVhVX5KlJuCj2kLkzfL1LKwmCrS11Y7dbbVIhRXVg/Pa4uVHFlIdSpLqyW3qZCjOrC+uhxdaGKKAuDvW51YbX0GhWiorqwbqzn1YWllIXQgLqwWnqTCjGqCzPCWq8uLKUshAbUhdXSjirEbspCM1tNVBdmhlqvLiylLIRG1IXVksXotczIu1YVYjFlYX+KqNaI6sJeoS6kgrIw2OtSF1bdZk0O8lIqxP7B/iZdVYhFfwSUqC6syRdtpC5M3y9SysJEnNSsLqx2a6tlfUlX4HqTG1L2efgccyY/EBrJDknzgefN7Kys2imMddYAV5vZ1EbPsdU/E1eJLspCtv0IaF2qtUjzkDQE+ByeBs0I5rrVhfgK8CO4qvCHVKEurLr+NuvJSyoLI+2F2kBdWHX97RTkkUgzyMU/q41EyhGDPJJ7YpBHck8M8kjuiUEeyT0xyCO5JwZ5JPfEII/knhjkkdwTgzySe2KQR3JPDPJI7olBHsk9McgjuScGeST3xCCP5J4Y5JHcE4M8kntikEdyTwzySO6JQR7JPf8HnsVar1AJE4MAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
TPss_svhh_s
1262.162.57478e+061271.96NaNNaN269834NaN
2264.162.57478e+061914.6NaN0.0147807438320NaN
3367.8456.89339e+062032.341914.60.00827627521470480011
4298.166.89339e+061229.75NaNNaN269834NaN
\n", "
" ], "text/plain": [ " T P s s_s v h h_s\n", "1 262.16 2.57478e+06 1271.96 NaN NaN 269834 NaN\n", "2 264.16 2.57478e+06 1914.6 NaN 0.0147807 438320 NaN\n", "3 367.845 6.89339e+06 2032.34 1914.6 0.00827627 521470 480011\n", "4 298.16 6.89339e+06 1229.75 NaN NaN 269834 NaN" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = plt.axes([0,0,0,0.1])\n", "ax.axis('off')\n", "plt.text(0, 1.0, \"$\\eta_{c}$=\"+'{:.02f}'.format(eta_c), size=14)\n", "plt.text(0, 0.5, \"$\\dot{W}_{s,c}/\\dot{m}$=\"+'{:.02f} J/kg'.format(W_dot_s_c_m_dot), size=14)\n", "plt.text(0, 0, \"$\\dot{W}_{c}/\\dot{m}$=\"+'{:.02f} J/kg'.format(W_dot_s_c_m_dot), size=14)\n", "plt.show()\n", "display(states)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 飽和蒸気曲線の作成" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "P_list = list(range(int(5.2*1e+5), int(7.377*1e+6), int(1e+3)))\n", "saturation_curve = {}\n", "saturation_curve['Ts'] = {'s':[], 'T':[]}\n", "for P in P_list:\n", " saturation_curve['Ts']['T'].append(CP.PropsSI('T', 'P', P, 'Q', 0, params['WF']))\n", " saturation_curve['Ts']['s'].append(CP.PropsSI('S', 'P', P, 'Q', 0, params['WF']))\n", "for P in P_list[::-1]:\n", " saturation_curve['Ts']['T'].append(CP.PropsSI('T', 'P', P, 'Q', 1, params['WF']))\n", " saturation_curve['Ts']['s'].append(CP.PropsSI('S', 'P', P, 'Q', 1, params['WF']))\n", " \n", "saturation_curve['ph'] = {'p':[], 'h':[]}\n", "T_list = list(range(217, 304))\n", "for T in T_list:\n", " saturation_curve['ph']['p'].append(CP.PropsSI('P', 'T', T, 'Q', 0, params['WF'])/1e6)\n", " saturation_curve['ph']['h'].append(CP.PropsSI('H', 'T', T, 'Q', 0, params['WF'])/1e6)\n", "for T in T_list[::-1]:\n", " saturation_curve['ph']['p'].append(CP.PropsSI('P', 'T', T, 'Q', 1, params['WF'])/1e6)\n", " saturation_curve['ph']['h'].append(CP.PropsSI('H', 'T', T, 'Q', 1, params['WF'])/1e6)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Chart\n", "### T-s" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdUAAAEOCAYAAADfQI0bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xl8VOW9+PHPd5KQfWPJIlsEgmHfgmHfd1EUioJ1qUrx/rStVv313l5/va213ttbr9rbqm1ptVZrRVBZVFCpEhAFZIewhyUQCIQlkISQkGSe3x9zggESMpNMZv2+ec1rZs6c55zvkzPkm+ec5zyPGGNQSimlVNPZvB2AUkopFSg0qSqllFJuoklVKaWUchNNqkoppZSbaFJVSiml3ESTqlJKKeUmmlSVUkopN9GkqpRSSrmJJlWllFLKTUK9HYA7tW7d2qSlpTWq7IULF4iOjnZvQD5O6xwctM7BoSl13rRp02ljTJum7H/Tpk1JoaGhfwF6ErgNNjuQU1VVNWfAgAGFda0QUEk1LS2NjRs3NqpsdnY2o0aNcm9APk7rHBy0zsGhKXUWkbym7j80NPQvKSkp3dq0aVNks9kCcvxbu90up06d6n7ixIm/ALfVtU6g/jWhlFLKs3q2adOmOFATKoDNZjNt2rQ5j6M1Xvc6HoxHKaVU4LIFckKtYdWx3typSVUppZRyE02qSimllJtoUlVKKRUQXnjhhdYZGRndMzIyuttstgE1r+fMmdPOUzEEVO9fpZRSwevJJ588/eSTT54+dOhQ2NChQzP27Nmzy9MxaEtVKaVUQNm8eXPkTTfddNEb+9aWqlIqoL376Q9JiGjNxJE/93YowePBB9uTkxPl1m327FnG668fdWbVbdu2RXbr1s0rSVVbqkqpgFVlr+KV49l8fuSf3g5FedDOnTsje/XqdUVSveuuuzrWfj9jxoy02u/ff//9uFdeeaVlU/etLVWlVMDacGIDRTaYENPJ26EEFydblM1lz549kU899dTJmvelpaWSm5sb8cQTT9ywb9++iH/84x+Ho6Ki7OXl5TJ37tz2cXFx1Zs2bYp+8cUXj373u9/tkJiYWL1+/fqYZcuW5bZp06balX1rS1UpFbA+PfARkXY7w1r18nYoykOqq6vJy8uL6Nu3b3nNsrVr10ZNmzat6MUXXzweFxdXvXbt2qj+/fuXPf/8820eeuih0y+//PKxS5cuyWeffRY7Z86c0y+//PKxhISEKlcTKmhSVUoFqCp7FV/kZzOy7CIRCR28HY7ykJ07d4YnJydfioyMvDy609q1a6P79u17EaCsrMy2bt26qCFDhlzYunVr1MCBAy8WFRXZEhMTq7Zv3x6VlZV18fz587akpKTKxuxfT/8qpQLShhMbKLpUzMQLZRCb6u1wlIf07t274sCBAztrL9u1a1fkiRMnwubPn584Z86c02+99VbLp59+unDChAnF9957b8fIyEh7ly5dKvr161c2a9astNjY2Op+/fqVNWb/mlSVUgHps7zPiLSFMexiOcTd4O1wlBfNnz//ill4pk6dWgLw8MMPn3344YfP1ix/+eWXW6WmplYaY3jkkUfONGZfmlSVUgGn2lTzed7njIxsS4Q5oC1V5ZQf/OAHjUqktXnsmqqIRIjINyKyTUR2isgz1nIRkedEZJ+I7BaRH9Va/jsRyRWR7SLS31OxKqX8W255LkUVRUywxUNEPLRw7y2TStXHky3VCmCMMaZURMKANSKyHOgGtAcyjDF2EUmy1p8MpFuPLOAP1rNSSl3XlrItRIZGMqy8CuLaejscFUQ81lI1DqXW2zDrYYD/A/zSGGO31iu01pkGvGmVWwckiIiew1FKXVeVvYptZdsY0W4EkaUn9NSv8iiP3lIjIiEishUoBFYYY9YDnYG7RGSjiCwXkXRr9bZA7RuI861lSilVr40nN1JqL2VCxwlQXABxmlSV53i0o5IxphroKyIJwCIR6QmEA+XGmEwRmQ68DgwHpK5NXL1AROYCcwGSk5PJzs5uVGylpaWNLuuvtM7BIdjqPP/MfMIIQ3KrMaWF5J29xOEgqH+wHWdf5ZXev8aYcyKSDUzC0QJ93/poEfBX63U+jmutNdoBx+vY1jxgHkBmZqYZNWpUo2LKzs6msWX9ldY5OARTnSvtlfxswc/oFdWL8QN7wxo7ab0Gk5Y5ytuhNbtgOs6+zJO9f9tYLVREJBIYB+wBFgNjrNVGAvus10uB+6xewIOA88aYAk/Fq5TyP2uPr+VcxTkyozOhxPp1Eav3qCrP8WRLNRX4m4iE4EjmC4wxH4nIGuBtEfkxUArMsdZfBkwBcoEy4AEPxqqU8kMfHfyI+PB4ukV2g2LrxJZeU1Ue5LGkaozZDvSrY/k54JY6lhvgUQ+EppQKAGWVZWQfzWZqp6mEVoRqS1V5hQ6or5QKCF8c/YKLVRe5pZP1N3rxcbCFQVQr7wamPOaFF15onZGR0T0jI6O7zWYbUPN6zpw57TwVgw5TqJQKCB8f/JjU6FT6JfVj9e7VjpZqbCrYtO0QLJ588snTTz755OlDhw6FDR06NGPPnj27PB2DftuUUn7vzMUzrD2+lsk3TsYm1q+14uN6PTVIbd68OfKmm2666I19a0tVKeX3Psv7jGpT/e2pX3Ak1RSdnNwbHnyQ9jk5uHXA5Z49KXv99SsGBKrXtm3bIrt16+aVpKotVaWU3/v44MekJ6bTNbGrY4ExjtO/OuVbUNq5c2dkr169Gp1U77rrro6NLastVaWUXztacpRtp7bxeP/HLy8LrboAlWWaVL3E2RZlc9mzZ0/kU089dbL2sieffDL17NmzoQkJCdU//vGPC2fMmNF54sSJ53bt2hU5ePDg0s8//zzu2WefPd6tW7eK3NzciCeeeOKGffv2RSxduvSgzYXr8tpSVUr5tY8PfgzAlBunXF7W4pI177QOph90qqurycvLi+jbt295zbJDhw6FVVZWSkJCQvWGDRui161bF3XbbbcV/epXvzpZXFwc8sQTT5yeMWNG0YEDB8LXrl0bNW3atKIXX3zxeFxcXPXZs2dDXNm/JlWllN+yGztLcpeQlZJFasy3CTS8wpprWluqQWfnzp3hycnJlyIjIy+PFf+Tn/yk7TPPPHPi+9///unU1NTKLVu2RE2dOrW4oqJCWrZsWRUSEkJOTk7kgAEDLq5duza6b9++FwHKyspsrVu3rnZl/3r6Vynltzad3ER+aT6P9H3kiuWXk6q2VINO7969Kw4cOLCz9rLu3btf/OUvf5ly5syZ0L59+5Zt27Ytsnfv3uXr16+PzMjIKAfIy8sLT09Pv7Rr167IEydOhM2fPz9xzpw5p13dvyZVpZTfWpy7mJiwGMZ1HHfF8vAKPf2rvvXMM8+crGv5kCFDLg4ZMuQiwOLFiw8BzJ8/P68p+9LTv0opv3Sh8gIr8lYwMW0ikaGRV3zW4tIZiGwJYRFeik4FK02qSim/9Nnhz7hYdZHbu9x+zWfhFWf0eqryCk2qSim/tDh3MWlxafRp0+eaz8IrzuqpX+UVmlSVUn4nrziPzYWbub3L7YjINZ87WqqaVJXnaVJVSvmdJblLsImNWzvfeu2H1ZWEVZ7XKd+UV2hSVUr5lSp7FUsOLGHoDUNJikq6doWSEwhGW6rKKzSpKqX8ypf5X1JYVsiM9Bl1r1B83PGsLVXlBR5LqiISISLfiMg2EdkpIs9c9fnvRaS01vtwEXlXRHJFZL2IpHkqVqWU71qwbwFJkUmMbD+y7hVKrKSqLVXlBZ5sqVYAY4wxfYC+wCQRGQQgIplAwlXrPwQUGWO6AC8B/+3BWJVSPii/JJ+vjn3FjK4zCLXVM3ZNcYHjWVuqygs8llSNQ01LNMx6GBEJAZ4HfnJVkWnA36zX7wFjpa5ufkqpoPHevvcQEaanT69/pZLj2CUMolp6LjDlM2bOnJnWsmXLPunp6T2cLZObmxuWlZXVtVOnTj26dOnS49lnn63jYr1zPHpNVURCRGQrUAisMMasB34ALDXGFFy1eltwTB9kjKkCzgOtPBmvUsp3VFZXsih3ESPbjSQlOqX+FYsLqAhvBfo3eFB68MEHTy9dunS/K2XCwsJ44YUX8g8ePLhzw4YNu1977bWkTZs2NWo4Lo+O/WuMqQb6ikgCsEhERgAzgVF1rF7X/whzzUoic4G5AMnJyWRnZzcqttLS0kaX9Vda5+AQKHXedGETZ8vPklGRcd369D26m+rQeNYHQJ1dESjHuakmT55cunfv3hYNrbd+/frI+Pj46oyMjEsdO3as7NixYyVAYmKivXPnzhePHDnSYsCAAeUNbedqXhlQ3xhzTkSygdFAFyDXOrMbJSK51nXUfKA9kC8ioUA8cLaObc0D5gFkZmaaUaNGNSqm7OxsGlvWX2mdg0Og1PnNT9+kbUxb/mXSv2CT65xk21bGyfB2AVFnV/jScf6/721rv+9ESZQ7t9k1Jbbs+e/0cdvk5xcuXLDde++9Ny5dujQ3IyPjUs3yvXv3tti1a1fUyJEjS69Xvj4eS6oi0gaotBJqJDAO+G9jTEqtdUqthAqwFLgfWAt8B/jCGHNNS1UpFfgOnjvIhhMbeLz/49dPqMZASQGXUq4dulAFr1dffbXlSy+9dE138FOnToXNnDmz844dO3YDnD9/3jZ9+vTOv/71r4+2bNnS3ph9ebKlmgr8zeqYZAMWGGM+us76rwFviUgujhbqLA/EqJTyQW/vfpsWthZ1Dp5/hYtFUFVORbh2UvImd7Yo3eGRRx45+8gjj1xxpnP//v0tbr311i4vvPDCEYCKigq55ZZbOs+cOfPs/ffff66x+/JYUjXGbAf6NbBOTK3X5Tiutyqlgtj5ivMsPbCUWzrdQqvIBvoqljj6O1aEa59GdX05OTkRv//97/PGjx9/wW63M2vWrI5du3Yt/8UvflHn3KvO0hGVlFI+beG+hZRXl3NP93saXrlYk2qwu/XWW28cNmxYxqFDh8KTk5N7v/TSS63rWu+OO+4oHj9+/AWAFStWxCxevLjVmjVrYjMyMrpnZGR0f/fdd+Mbs3+vdFRSSilnVNoreWfPO2SlZtE1sWvDBazRlC610NO/werDDz885GqZiRMnlhpjNrlj/9pSVUr5rBWHV1BYVsi93e51rsDllqomVeUdmlSVUj7JGMNbu96iY1xHhrcb7lyh4mMQ1RpjC2ve4JSqhyZVpZRP2nZqGzlncvhut+9e/zaa2koKdCB95VWaVJVSPumvOX8ltkUs0zpPc75QcYEOpK+8SpOqUsrnHDh3gC+OfsHdGXcTFebCwDwlx7WlqrxKk6pSyue8tuM1IkMj+W637zpfqKoCys5AXNvmC0ypBmhSVUr5lPySfJYdWsZ3un6HxIhE5wuW1Myjqi1V5T2aVJVSPuWNnW9gExv3d7/ftYI1k5Pr6V/lRZpUlVI+4/TF0yzav4hpXaaRHJ3sWmFr4AftqBS8GjvZeFlZmfTq1avbTTfd1L1Lly49fvzjHzf6S6QjKimlfMYbOW9QZap4sMeDrhe+oqVa6Na4lH+omWx82LBhZUVFRbZ+/fp1nzJlSnFD86JGRESYNWvW7I2Pj7dXVFTIwIEDb/r888/Pjx079oKrMWhLVSnlEwrLCpm/dz5TO02lfVx71zdQUgChkRCR4P7glF/o2LFj5bBhw8rgysnG61p3/fr1kXv27GkBYLPZiI+PtwNcunRJqqqqxJrj22XaUlVK+YR52+dRba/mX/r8S+M2UGzdTtPIX4bKjRY/2p7CXW6dpJyk7mXc/orTU8o1NNn41ZOUV1VV0bNnz+5HjhwJv//++wvHjBnjcisVNKkqpXzAsdJjvL//faanT6d9bCNaqeBoqer1VMW1k407M0l5aGgoe/bs2XX69OmQW265pfOGDRsiBg4ceN3TxnXRpKqU8ro/bP0DIRLC3N5zG7+R4uPQ/mb3BaUaz4UWpbvVNdm4M5OU12jdunX1sGHDSj788MP4xiRVvaaqlPKqg+cP8uHBD7nrprtc7/Fbwxirpaq30wQzVyYbr5mkfMKECReOHz8eevr06RCA0tJSyc7OjuvWrZvLCRWcaKmKSDtjTH49n91qjPnQmR2JSASwGgi39vueMebnIvI2kAlUAt8ADxtjKq2rxP8LTAHKgO8ZYzY7sy+llP/47abfEhESwUO9Hmr8RsrOQPUliNPTv8GsZrLx9PT0ixkZGd0BnnnmmWN33XXX+avXveOOO4prXh89ejTse9/73o3V1dUYY2TatGlnZ8+efU0ZZzhz+vdzEZlojDlce6GIPAg8DTiVVIEKYIwxplREwoA1IrIceBu4x1rnH8Ac4A/AZCDdemRZy7Kc3JdSyg+sL1jPyqMreaz/Y7SMaMIcqMU196hqSzWYNXay8aysrIu7d+/e5Y4YnDn9+2NghYik1ywQkZ9ay0c6uyPjUNMLK8x6GGPMMuszg6Ol2s5aZxrwpvXROiBBRPR/jFIBotpezW82/IYbom/g3u5OTkJen5ohCrWlqryswaRqjFkG/AuwXER6ishvganAiPpOC9dHREJEZCuOO7NXGGPW1/osDLgX+MRa1BaofbE731qmlAoAi3IXsa9oH09kPkF4SHjTNqYtVeUjnOr9a4z5XES+B2QDXwNjjTEuX8Q1xlQDfUUkAVgkIj2NMTnWx68Cq40xX1rv67rZzFy9QETmAnMBkpOTyc7OdjUsAEpLSxtd1l9pnYODL9b5ov0iLx57kU7hnWhxqAXZh7ObtL20Q2vpiLB6816M7YBP1rm5BWOdfZEzHZVKcCQzwdHJaCxQaHUkMsaYOFd3aow5JyLZwCQgR0R+DrQBHq61Wj5Q+4a1dsDxOrY1D5gHkJmZaUaNGuVqOABkZ2fT2LL+SuscHHyxzs+te45Seyl/Hv9nerTq0fQNnn8PziQxcsw4wDfr3NyCsc6+yJnTv7HGmDjruYUxJrrWe6cTqoi0sVqoiEgkMA7YIyJzgInAbGOMvVaRpcB94jAIOG+MKXCpdkopn7Pt1Dbe3fsuszNmuyehgt5Oo3yGJwd/SAX+JiIhOJL5AmPMRyJSBeQBa62xFj8wxvwSWIbjdppcHLfUPODBWJVSzaDSXskza5+hTVQbftjvh+7bcHEBJHZ03/aUaiSPJVVjzHagXx3L64zB6g38aHPHpZTynDd3vsn+ov38dvRviWkR474NlxyHDoPctz2lGqnB078i8j+eCEQpFdgOnDvAq1tfZUz7MYztMNZ9G668CBeLdHJy5ROcuU91TLNHoZQKaJXVlfzbl/9GdFg0Pxv8M/duvOYeVR1MP+g1dbLxqqoqunXr1n306NFdGhuDDqivlGp2r257lT1n9/C/o/+X1pGt3bvxKyYnV8GsqZON/+pXv0ru0qXLxdLS0pDGxuBMS7WPiBwSkaUi8p8iMltEelmDNSil1HVtPLGR13NeZ3r6dMZ0aIYTX9pSVRZXJhuvPUk5wIEDB8I+/fTT+O9///unmxKDMy3V7cAtQE+gFzABeBJIF5GjxpieTQlAKRW4TpWd4v+u/r+0j23PTwb+pHl2UjOakrZUfcbPvvpZ+9yiXLdOUt4lsUvZs0OfbXBKOWcnG796kvJHH320/W9+85v88+fPN7qVCs6PqHQcx8ALn9Uss9J/o887K6UCW2V1JU+uepILlReYN34e0WHRzbOj4uMQFg3hLo9DowJQXZONb9iwIep6k5T/+7//+/HWrVtXDR8+vOyjjz6KbdL+nVjnlboWWre87G/KzpVSgev5jc+zpXALvxnxG9IT0xsu0Fglxx2t1HpO8ynPc6ZF2dxqTzb+y1/+8uT1JilfsmRJwooVKxLatm0bX1FRYbtw4YJt2rRpNy5ZsuSQq/t1ZkSlv7i6UaVUcHtr11u8s+cd7ut+H5NvnNy8OyvW0ZSUgyuTjdeepPyVV145dvLkye3Hjh3b8cYbbxwcNGhQSWMSKmjvX6WUm316+FOe3/A84zqM44kBTzT/DksKoOOQ5t+P8nmuTDZee5Jyd9KkqpRym6+OfcVPv/wpfZP68l/D/4sQW5P6fDTMbnckVZ1HVeGeycanTp1aMnXq1JLGlnfmlhrA0TFJRO4Rkf+w3ncQkZsbu2OlVGBZnb+aH33xIzondOb3Y35PRGhE8++07DTYq/R2GuUznE6qOOY7HQzMtt6XUE8nJqVUcPln3j95fOXjdE7ozF8m/IX48HjP7Fhvp1E+xpXTv1nGmP4isgXAGFMkIi0aKqSUClzGGP6686+8tOklerfpzatjX/VcQgUd+EH5HFeSaqU1bZsBx/yogP36RZRSgepC5QWeW/ccHx78kElpk3h26LOeOeVbm7ZUlY9xJan+DlgEJInIc8B3gP/XLFEppa5RXV1NeXk55eXlVFdXA45h2WoeInLF+4iICOoboq2pck7n8K+r/5X80nwe6fMID/d5GJu4cjXJTUoKQGwQneT5fStVB6eSqjV60mpgEzAWEOB2Y8zuZoxNqaB14cIF8vLyyMvLo7CwkDNnzlBS4lqHRJvNRnR0NMYYCgsLadOmDUlJSbRt25aEhIRGxVV8qZiXt7zMu3vfJSkqidcmvEZmSmajtuUWxQUQkwwheiOD8g3ODlNoRGSxMWYAsKeZY1IqKFVUVJCTk0NOTg6HDx8GHEOupaSk0KlTJxISEoiMjCQiIoLQUMd/XbvdfvlhjLn8urq6mrKyMkpLSzly5AgnT55k9+5v/wZOSEggLS2NjIwMOnfufHl79Sm5VMK7e9/lzZ1vcv7See7seic/6PcDz14/rTOw4zrwg/Iprvx5t05EBhpjNjRmRyISgaO1G27t9z1jzM9F5EZgPtAS2Azca4y5JCLhwJvAAOAMcJcx5nBj9q2ULysuLmb9+vVs2rSJiooKWrVqxYgRI0hPTyc1NZWQkKbd65mdnc2oUaOorKzk9OnTHDlyhLy8PPbs2cPWrVsJDw+ne/fuZGVlkZycfLmcMYbtp7fz4YEP+fjgx5RWljKs7TAe6/8YGS0zmlpt9ygugFadvR2FUpe5klRHAw+LSB5wAccpYGOM6e1k+QpgjDGm1Jo2bo2ILAeeAF4yxswXkT8CDwF/sJ6LjDFdRGQW8N/AXS7Eq5RPyz+Xz2drPyMnJ4fq6mpuvPFGBvQcQOvWrRERTnGKwtOFABgMjuG2Ha9rGGMuv79mHeN43n1xN6HHQr/9rK0hpW0KbbLacOLkCfLy8vhk/yd8vOdjkpKTaJvelqP2o6w7vo7Ci4VEhEQwpsMYvtfje3Rr1c2TP6KGlRyHtGHejkL5mKqqKnr16tU9JSXl0sqVK3OdKdO2bdte0dHR1TabjdDQUJOTk9Ooy5uuJNUmDeBpDcBfar0Nsx4GGAPcbS3/G/ALHEl1mvUa4D3gZRERU/ObQSk/dKn6Eh8f/Jh3tr/D7pLdjj9Na/rYlAMbm2nH/2zg89rzhu+HSCK5OeVmftj/h4zrMI6YFjHNFFgTXLoA5ee156+6RmMnG1+1atW+1NTUqqbs2+mkaozJa8qOAKxbcjbhmDLuFeAAcM4YU1OJfKCt9botcNTad5WInAdaAU2aQFYpbyivKuedPe/wt5y/cabiDLGVsQy2DWbygMm0atUKQS731BXrn/Xm8muRb5fXXr/m/TXrWK83b95M//796y0vjp0AUF1VzaHdh9i9fjeVRyqJDY0lvGN48/5wGqtY71FV16qZbPynP/1pwUsvvZRc33rr16+PjI+Pr87IyLjkzv07nVRrhie8mjHml85uwxhTDfQVkQQct+fUdS6ppiVa170A17RSRWQuMBcgOTmZ7OxsZ8O5QmlpaaPL+iutc/MzxrC1bCuLixZztvosyeXJjDw/kqzULNq1a4cUCfai5r3dO6kqiXO7zjm9fjzxDOg/gAMHDrBmzRo2bdpEt27diI1t0jSTbpdQtJ2+wNZDhZw7l33FZ/rd9q7j//50+4r9+906SXl4enrZDf/5XINTyjk72fjVk5QDjB07Nl1EeOCBB0499dRTjWrAuXL6t/bs6RHAVKBR55yNMedEJBsYBCSISKjVWm2HYzJ0cLRa2wP5IhIKxANn69jWPGAeQGZmphk1alRjQrrcmSOYaJ2b15FvvuD9Fb/j9XYHSA1NZcSJEfRt2ZcZc2bQsmVLj8QAja/z+PHjOXjwIEuWLGHbtm1MmDCBgQMHNtu9ry7bdhK2Qd9hk6FN1ys+0u92cHrnnXfi65ps/NVXX215vUnKd+zYsfurr77ak5aWVnns2LHQMWPGdO3Ro0f55MmTS68u0xBXTv++UPu9iPwPsNTZ8tYITJVWQo0ExuHofLQSx0AS84H7gSVWkaXW+7XW51/o9VTlTza//XtGf76XQ0+Mo+XhVmTdnMX48eMbvH3Fl3Tq1ImHH36YxYsXs3z5cvLz87ntttt8ow4lOpqSr3KmRdkc1qxZE1PfZOPXm6QcIC0trRKgbdu2Vbfccsu5tWvXRjcmqTZlCJQooJML66cCK0VkO7ABWGGM+Qj4V+AJEcnFcc30NWv914BW1vIngH9rQqxKeVz3zsMJr4L2eS2Z+Z2ZTJ482TeSkYuioqKYPXs2Y8aMYceOHbz11luUlZV5OyzHNdXwOAj3rdPSyntcmWy89iTlxcXFtqKiIhtAcXGxbeXKlXG9e/e+2JgYXLmmuoNvr2mGAG2AZ50tb4zZDvSrY/lB4Jop5Iwx5cBMZ7evlC/Jzc1l7c593AzMHj+edj16eDukJhERhg8fTmJiIosXL+b111/nvvvuIy4uzntB6cAPqglqT1Ken58fescdd3QBqK6ulhkzZpz5zne+06hJzF35s3lqrddVwMlavXaVUpYtW7bw4Ycf0jnF0fEwwRqnNxD07NmTuLg4/vGPf/DGG29w//33Ex/vpVGVigv01K+qlyuTjXfv3v3S3r17mzS5eQ1XTv8+YozJsx7HrNtc/tsdQSgVKNatW8fSpUvp1KkTdzw0B4DKEye9HJV7dejQgXvuuYeysjLeeOMNzp1zvmexW5UU6O00yue4klTH17GsSQNCKBUojDFkZ2fz6aef0q1bN2bNmkXkDalgs1F18oS3w3O7du3acd9991FeXs4DrmqvAAAgAElEQVTf//53z19jtVdDyQltqSqfIw11qBWR/wM8gqNT0oFaH8UCXxtjvtt84bkmMzPTbNzYiCFpPvmEsrlziYpy621VPq+srEzr7AYG+KxfP9ZlZND34EFu/eYbbNb/q/1ANODN9lRzHucjrVvz1ujRJJ87x31ffEELT53qjqiG2wthYxzkRl/zcTB+t3PHjKHLq682qqyIbDLGNGm6oW3bth3u06dPUAzOs23bttZ9+vRJq+szZ66p/gNYDvwXV/bALTHGXHPfqF9KSKC0SxeikoJrTsbSwkKtcxMZ4NOkJNa3asXNZ88yqaIC6dPn8ueh+/ZTFWKDzt4b9L05j3MHYEZBAQvateO9yZOZlZ/fpFsKnBZxHiiElC6Oqd+uEozf7UsevPe5Hna73S42my2gb3202+0C1DtiS4NJ1RhzHjgPzBaRRCAdx+APiAjGmNVuitV7Bg1i13/8B0lBduP0ruxsrXMTGGP45z//yfqvvyYrK4uJEydeMzBC2A9/RMXBgzB/vlv22RjNfZwzgMkbNrBs2TI+f+ghxo+v60qRm+35GOZ/Dc+8AG37X/NxMH63C7Oz6e7dEHJOnTrVvU2bNucDNbHa7XY5depUPJBT3zqu3FIzB3gMx6hHW3GMhrQWx4D4SgWdlStX8vXXX5OZmVlnQgUITUnhwtdfeyE6zxo4cCCFhYV8/fXXpKam0rNnz+bdYXHNwA/aUclXVFVVzTlx4sRfTpw40ZOmjYHgy+xATlVV1Zz6VnDllprHgIHAOmPMaBHJAJ5pYoBK+aVVq1bx5Zdf0q9fP6ZMmVLv0H1hKcnYL1ygurSUkBgfnOnFjSZNmkRhYSFLliyhdevWpKSkNN/Oio+DhEB0m+bbh3LJgAEDCoHbvB2Ht7ny10S5NSADIhJujNkD3NQ8YSnlu9atW0d2djZ9+vTh1ltvve5YuKHJjsRSdSLwegBfLSQkhJkzZxIZGcnChQu5dMmtk39cqaQAYlPA1rQJ3JVyN1eSar41u8xiYIWILOHbwe+VCgrbt2+/fNvMbbfd1uDg8mHWABCBdq9qfWJiYpg+fTpFRUUsX768+XZUrKMpKd/kVFIVx2+OHxljzhljfgH8DMfYvLc3Y2xK+ZTc3FyWLFlCWloa06dPx2Zr+L9PqHUKNBDvVa1PWloaw4cPZ+vWrezYsaN5dlKioykp3+RUUrVmh1lc6/0qY8xSY0wznt9RynccO3aMBQsW0KZNG+666y6nB8YPtW7rqDwZHC3VGiNHjqR9+/Z89NFHzTPiUnEBxLV1/3aVaiJXTv+uE5GBzRaJUj7q9OnTvP3220RHR3PPPfcQERHhdFlbixaEtGpFVZCc/q1hs9mYPn06AEuXLsWtszZWlMClEj39q3ySK0l1NI7EekBEtovIDmsaN6UCVklJCX//+9+x2Wzce++9xDSiB29YcjKVQXT6t0ZCQgITJkzg0KFDbN682X0bLi5wPOvtNMoHuXJLjY7zq4LKpUuXmD9/PmVlZXzve9+jZSNHrAlNSaHy2DE3R+cf+vfvz86dO/nss8/o0qWLe2a0qZmcXFuqyge50lI9AgwH7jfG5OEYoe3a8cGUCgB2u51FixZRUFDAjBkzuOGGxreKwlKSg+KWmrqICLfeeivGGD766CP3nAbWlqryYa4k1VeBwcBs630J8IrbI1LKB/zzn/9kz549TJgwgZtuatrt2KHJKVSfP4/94kU3RedfEhMTGTt2LLm5uezcubPpG9SWqvJhriTVLGPMo0A5gDGmCGjhbGERaS8iK0Vkt4jsFJHHrOV9RWSdiGwVkY0icrO1XETkdyKSa13DvXaAT6WawcaNG1m7di0DBw4kKyuryduruVe1Ksh6ANc2cOBAUlNT+eyzz6ioqGjaxooLICIeWgTXLDTKP7iSVCtFJATHaV9EpA3XGam/DlXAk8aYbjjGDX5URLoDvwGeMcb0Bf7Deg+Oa7jp1mMu8AcX9qVUo+Tm5rJs2TLS09OZNGlSg4M7OKNmVKVgGQCiLjabjSlTplBSUsKqVauatjGdnFz5MFeS6u+ARUCyiDwHrAH+09nCxpgCY8xm63UJsBtoiyNJx1mrxfPtKE3TgDeNwzogQUT0fI9qNoWFhSxcuJCkpCRmzJjh1OAOzvi2pRqc11VrtGvXjn79+rFu3ToKCwsbv6Hi4zrwg/JZTv/WMMa8DfwERyI9DtxujFnYmJ2KSBrQD1gPPA48LyJHgf8Bfmqt1hY4WqtYvrVMKbe7cOEC//jHP2jRogV333034eHhbtt2aHJwDVV4PePGjSMiIoJly5Y1vtOStlSVD3Nl6rcIYAqOHsB2oIWIHKoZZN+F7cQA7wOPG2OKReRXwI+NMe+LyJ04hj8cB9R13u2a/4UiMhfH6WGSk5PJzs52JZzLSktLG13WX2mdHex2O9u3b6ekpIS+ffu6955KS5voaPI2byLHCz9vXzvO7dq1Y//+/bz33nu0aePaLDNir2ZEyUnyii5x+Dp18rU6e0Iw1tkXuXKf6ps4evz+zno/G3gLmOnsBkQkDEdCfdsY84G1+H4c08oBLAT+Yr3OB9rXKt6OOgbwN8bMA+YBZGZmmlGNnJg4Ozubxpb1V1pnh48++ojz588zffp0evXq1Sz7PdiuHXEhoQzwws/b146z3W7nT3/6EwUFBdxxxx1OD/kIwPljsNpOWs9BpA0cVe9qvlZnTwjGOvsiVy4a3WSMecgYs9J6zAW6OlvYGpT/NWC3MebFWh8dB0Zar8cA+63XS4H7rF7Ag4DzxpgCF+JVqkEbN25k06ZNDBkypNkSKkBoEN+rejWbzcb48eMpKipiw4YNrhUu0XtUlW9zpaW6RUQGWZ2GEJEs4CsXyg8F7gV2iMhWa9m/A98H/ldEQnHcrjPX+mwZjtPNuUAZ8IAL+1KqQXl5eSxfvpz09HTGjh3brPsKS06hPMcN92gGiC5dutC5c2dWr15N3759iYyMdK5gsXWySpOq8lGuJNUsHC3HI9b7DsBuEdmBYyKb3tcrbIxZQ93XSQEG1LG+AR51IT6lnHbu3DkWLFhAYmKi09O4NUVoSjLVZ85gv3QJWwunb+8OaBMmTOCPf/wjq1atYtKkSc4Vqmmpakcl5aNcSapOfuuV8m01Y/pWV1cza9Ysl2adaaww617VqsJCWrRr1+z78wdJSUn069ePDRs2cPPNNzs3tnLxcbCFQVSr5g9QqUZw5ZaavOs9mjNIpdzFGMOSJUs4efIkM2bMoHXr1h7Zb2jNvap6XfUKo0aNwmazOd9rtaTAMTxhM59ZUKqxnP5mikimiCwSkc069ZvyV0eOHGHXrl2MGzeO9PR0j+03rOZe1SAeqrAusbGxZGVlsWPHDucGhNCBH5SPc+XPvbeBvwIzgFuBqdazUn5h3759HD58mF69ejFkyBCP7js0xTr9qwNAXGPo0KGEh4ezcuXKhleuaakq5aNcSaqnjDFLjTGH9LSv8jdFRUUsWrSImJgYbr31VreM6euKkJgYbNHRQTlZeUMiIyMZPHgwe/bs4dj15p01xjGYvvb8VT7MlaT6cxH5i4jMFpHpNY9mi0wpN6msrGTBggUAdO/enbCwMK/EEZqSoi3VegwaNIjIyMjrt1YriqHygrZUlU9zJak+APTF0Qv4Vr49BayUzzLGsGzZMk6cOMEdd9zh/P2QzSAsOVlbqvUIDw9n2LBhHDhwgLy8ek6A6eTkyg+4klT7GGMyjTH3G2MesB4PNltkSrnB5s2b2bp1KyNGjKBrV6cHAGsW2lK9voEDBxITE8MXX3xR92D7Ojm58gOuJNV11vynSvmFY8eOsXz5cjp37szIkSMbLtDMwlKSqTp1ClNV5e1QfFJYWBjDhw/nyJEjdbdWL7dUNakq3+XK4A/DgPtF5BBQgWN0pAZHUvIHJSVw9Ggk+/Z5OxLPCuQ6V1SU8cknCwkPj6FPn+nk5jr+fvRmne2kgN3Ovm9OI61TPLZffzrOMTH9iIj4kuXLVzF2bNoVn7U8dJzWwP4TqZjT19+OP9XZXc6d805fAXUlHVEJ+OQTuO++LG+H4QWBWWcRO9/97gekpZXy+usP8JOfRNX61Ht1HhGdzB/bwZ1jT7C93HNJ1b+OcxiDBw9h4sTPGDfuCEePdrj8yau3HOfO7ol07eHMdXF/qrN7zJ7dnttv93YUypWkegT4LtDJGPNLEekApAB+f1tNVhY8/fQuuncPrrPbu3YFZp3PnFnF2bMHSEqayvPPXzmvvTfrHHEyBf4E//XkSc55MAR/O852+wAOH17DD3+4mrZt77m8fMSxAmyVN/D22w1vw9/q7A5lZYU4hmRX3uRKUn0Vx+TkY4Bf4phb9X1gYDPE5VEdOsC4cYWMGhVc/wmzswOvzvv27eOddxwzn9x2W3+uvh3Vm3WuPpfMvj/BoPQTtLzbc/v1v+PcgjVrBvP5558zcuQx2ra1/jD603GIvoG7nfjZ+V+dmy47u9TbIShc66iUZYx5FMf0bBhjigCdbkP5jJoBHlJSUpgyZYrHB3hoiC0+HomIoFJ7ADdo4MCBREZGsnr16m8XFhdoJyXl81xJqpUiEgIYABFpg6PlqpTX1R7g4c477/TaAA/XIyKEJSdTpfeqNig8PJysrCz27dtHQUEBVFfChVM65Zvyea6c/v0dsAhIEpHngO8A/69ZolLKRTUDPMyePZvExERvh1Ov0JSUgGqpVlbbqaiyU15Z/e1zpZ3yqtrP1fWuU15pp+Ly85XrlF+ykV/Rg6V/3Ej78AssxrCjJIpe3q60UtfRYFIVkVBjTJUx5m0R2QSMxXE7ze3GmN3NHqFSDagZ4GH48OFeH+ChIWEpyZRt2Oj27VbVldys591nqmFv4eUEVlErkdW1/jXJrcp+OTFW1HpfXmWn2l7HIA1OCrEJEaE2wsNCLj+H13qfEB1OaOsYzp0ppEdaa9ZdmEHL9v3d+FNTyv2caal+A/QHMMbsAfY0a0RKueD48eMsW7aMzp07M2rUKG+H06Cy+FZcOnmSpVvyKa82dbfirk5qV7XmKq74zPG6qqHktmFDvR/ZBCKshHb1c3hYCPGRYYTHhtf6zEZ4aAgRYTYiQkMID6u7XM37+tYJC2n46lNZWRm//e1v6Z4Yx6AHXnf1x62UxzmTVN3S20NE2gNv4rgNxw7MM8b8r/XZD4EfAFXAx8aYn1jLfwo8BFQDPzLGfOqOWFRgKCsrY8GCBcTExDB9+nRsfjBx9X4TRVp1NT//25cURcRd8ZkI3yagWs81SSw2IpTWtd5H1JGo6nres3MHN2f2v6Jc7fKhNvG5Tl01oqKi6NevHxs3bmT06NHEx8d7OySlrsuZpNpGRJ6o70NjzItO7qsKeNIYs1lEYoFNIrICSAamAb2NMRUikgRgDYk4C+gB3AD8U0S6GmOqndyfCmB2u51FixZRWlrKAw88QFRUVMOFfMCAe6ZTOHwQC7p0JiLCav1ZyTMspHmSm5wIYUBH373O3JDBgwezceNG1q5dy6RJATsGjQoQziTVECCGJrZYjTEFQIH1ukREdgNtge8DvzbGVFifFVpFpgHzreWHRCQXuBlY25Q4VGBYvXo1ubm53HLLLd/ex+gHkjreQFJH7cHqioSEBHr16sXmzZsZMWKE3/wBpYKT1DkbRO0VRDYbY9zaO0BE0oDVQE/reQmOYRDLgaeMMRtE5GVgnTHm71aZ14Dlxpj3rtrWXGAuQHJy8oD58+c3KqbS0lJiYmIaVdZf+Wudz5w5Q05ODsnJydx0000ute78tc5NEQh1vnDhAhs3bqRjx46kpaU1uH4g1NlVTanz6NGjNxljMt0cUlDy2DXVyxsTicExEtPjxphiEQkFEoFBOEZnWiAinerZ7zV/ARhj5gHzADIzM01jO6tkZ2f7RUcXd/LHOhcVFTFv3jySk5N56KGHXL4f1R/r3FSBUufi4mKOHj3K3XffTYsW1x93JlDq7IpgrLMvcqZnx1h37UxEwnAk1LeNMR9Yi/OBD4zDNzg6MbW2lrevVbwdcNxdsSj/U1VVxcKFCzHG+OwAD6r5DBs2jIsXL7JlyxZvh6JUvRpMqsaYs+7YkTjO0b0G7L6qc9NiHOMJIyJdcQx9eBpYCswSkXARuRFIx3F7jwpSy5Yto6CggOnTp9OyZUtvh6M8rH379nTo0IGvv/6a6mrtr6h8kyfvQRgK3AuMEZGt1mMK8DrQSURygPnA/VardSewANgFfAI8qj1/g9fmzZvZsmWLXwzwoJrPsGHDKC4uJicnx9uhKFUnV4YpbBJjzBrqvz57T10LjTHPAc81W1DKL9QM8NCpUye9ZhTkunTpQlJSEl999RW9e/f22ftrVfDy/bvlVVC7ePEiCxcuJDo6mhkzZvjFAA+q+YgIQ4cO5dSpU+zbt8/b4Sh1Df0NpXyWMYYPPviAkpIS7rzzTr0/UQHQs2dPEhIS+Oqrr7wdilLX0KSqfNaqVavIzc1l0qRJfjXAg2peNpuNwYMHc/ToUY4cOeLtcJS6giZV5ZNyc3NZtWoVffr0YcCAAd4OR/mYfv36ERUVpa1V5XM0qSqfc+7cOd5//32Sk5O55ZZbtDOKukZYWBg333wz+/bto7CwsOECSnmIJlXlU6qqqliwYIEO8KAaNHDgQMLCwvj666+9HYpSl2lSVT6lZoCHO+64Qwd4UNcVFRVF//792bFjB+fPn/d2OEoBmlSVD9myZQtbtmxh2LBh3HTTTd4OR/mBwYMHA7B2rU5epXyDJlXlEwoKCvj444/p1KkTo0eP9nY4yk/Ex8dfnhaurKzM2+EopUlVed/FixdZsGAB0dHRTJ8+XQd4UC4ZMmQIlZWVbNiwwduhKKVJVXmXMYZFixZRXFzMzJkziY6O9nZIys8kJSXRtWtXvvnmGyorK70djgpymlSVV61evZr9+/czadIk2rVr5+1wlJ8aOnQoZWVlOi2c8jpNqspr9u/fT3Z2Nn369CEzM9Pb4Sg/1qFDB9q3b8/atWux2+3eDkcFMU2qyiuKior44IMPdIAH5TZDhw7l3Llz7Ny509uhqCCmSVV5XGVlJQsXLtQBHpRbde3alTZt2vDVV19hjPF2OCpIaVJVHlczwMP06dN1gAflNiLCkCFDOHnyJEVFRd4ORwUpTarKozZv3szWrVsZPnw4Xbt29XY4KsD06tWLuLg4nb1GeY3HkqqItBeRlSKyW0R2ishjV33+lIgYEWltvRcR+Z2I5IrIdhHp76lYVfM4duwYy5Yto3PnzowaNcrb4agAFBISwqBBgzh//jz5+fneDkcFIU+2VKuAJ40x3YBBwKMi0h0cCRcYD9T+83IykG495gJ/8GCsys3KyspYuHAhMTExOsCDalb9+/cnNDRUp4VTXuGx32zGmAJjzGbrdQmwG6iZefol4CdA7d4F04A3jcM6IEFEUj0Vr3Ifu93OBx98QGlpKXfeeSdRUVHeDkkFsPDwcG644Qb27NnD6dOnvR2OCjLijV5yIpIGrAZ6AqOAscaYx0TkMJBpjDktIh8BvzbGrLHKfA78qzFm41XbmoujJUtycvKA+fPnNyqm0tJSYmJiGlXWX3mqzocOHeLIkSN07dqV1FTv/l2kxzk4FBUVkZOTQ1JSUtBMztCU4zx69OhNxhi9WdwNQj29QxGJAd4HHsdxSvhpYEJdq9ax7Jq/AIwx84B5AJmZmaax1+qys7OD7jqfJ+q8b98+Vq1aRd++fZk2bVqz7ssZepyDQ3Z2Nv3792fLli3cfffdxMbGejukZheMx9kXefTCloiE4UiobxtjPgA6AzcC26xWajtgs4ikAPlA+1rF2wHHPRmvapqzZ8/ywQcfkJqaypQpU7wdjgoyQ4YMwW63s27dOm+HooKIJ3v/CvAasNsY8yKAMWaHMSbJGJNmjEnDkUj7G2NOAEuB+6xewIOA88aYAk/Fq5qmsrKSBQsWICLMnDlTB3hQHpeYmEiPHj3YuHEj5eXl3g5HBQlPtlSHAvcCY0Rkq/W4XvNlGXAQyAX+DDzigRiVGxhj+Oijjzh58iTTp08nMTHR2yGpIDVkyBAuXbrExo0bG15ZKTfw2DVVq8PRdQd4tVqrNa8N8Ggzh6WawcaNG9m+fTujRo0iPT3d2+GoIJaamkrnzp1Zt24dgwYNIjTU491IVJDRmwWVWx09epRPPvmE9PR0RowY4e1wlGLo0KFcuHCBrVu3ejsUFQQ0qSq3KSkpYcGCBcTHx3PHHXfozDPKJ6SlpXHDDTfotHDKIzSpKreorq5m4cKFVFRUcNdddxEZGentkJQCHAPtDx06lLNnz7J7925vh6MCnCZV5RbLly/n6NGjTJs2jeTkZG+Ho9QVMjIyaNWqFWvWrNFp4VSz0qSqmmzz5s1s2rSJoUOH0qNHD2+Ho9Q1bDYbw4YN48SJE+zbt8/b4agApklVNUl+fv7lmWfGjBnj7XCUqlevXr1ITExk9erV2lpVzUaTqmq00tJSFixYQGxsLDNmzNCZZ5RPCwkJYdiwYRw/fpzc3Fxvh6MClP4WVI1SXV3NggULKC8vZ9asWdoxSfmFPn36EB8fz6pVq7S1qpqFJlXVKJ988glHjx7ltttu045Jym+EhIQwfPhwjh07xsGDB70djgpAmlSVy7Zs2cLGjRsZMmQIPXv29HY4SrmkT58+xMXFaWtVNQtNqsolx44d4+OPP6ZTp06MHTvW2+Eo5bLQ0FCGDRvG0aNHOXz4sLfDUQFGk6pyWmlpKe+++652TFJ+r1+/fsTGxrJq1Spvh6ICjP5WVE6pGTHp4sWL3HXXXURFRXk7JKUaLTQ0lKFDh5KXl6etVeVWmlSVUz755BOOHDnCbbfdRkpKirfDUarJ+vfvT0xMjLZWlVtpUlUN2rBhw+WOSb169fJ2OEq5RVhYGEOHDuXw4cPaE1i5jSZVdV2HDx++PJWbdkxSgSYzM5O4uDi++OIL7Qms3EKTqqpXUVERCxYsoGXLlkyfPl07JqmAExoaysiRIzl27JiOCazcwmO/JUWkvYisFJHdIrJTRB6zlj8vIntEZLuILBKRhFplfioiuSKyV0QmeipWBRUVFcyfPx9jDLNnzyYiIsLbISnVLPr06UPLli1ZuXKltlZVk3my6VEFPGmM6QYMAh4Vke7ACqCnMaY3sA/4KYD12SygBzAJeFVEQjwYb9AyxrBo0SJOnTrFzJkzadmypbdDUqrZhISEMGrUKE6ePMnOnTu9HY7ycx5LqsaYAmPMZut1CbAbaGuM+cwYU2Wttg5oZ72eBsw3xlQYYw4BucDNnoo3mK1cuZK9e/cyceJEOnXq5O1wlGp2PXv2JCkpiezsbOx2u7fDUX7MKxfJRCQN6Aesv+qjB4Hl1uu2wNFan+Vby1QzysnJ4csvv6Rfv37cfLP+DaOCg4gwevRozpw5w7Zt27wdjvJjoZ7eoYjEAO8DjxtjimstfxrHKeK3axbVUfyaCx4iMheYC5CcnEx2dnaj4iotLW10WX91dZ1LSkrYunUrcXFxAXv/nh7n4NCYOhtjiI2N5dNPP6WoqMjvOuYF43H2RR5NqiIShiOhvm2M+aDW8vuBqcBY821PgXygfa3i7YDjV2/TGDMPmAeQmZlpRo0a1ajYsrOzaWxZf1W7zqWlpfz5z38mJiaGuXPnEh0d7d3gmkmwH+dg0dg6d+jQgb///e9ER0eTlZXl/sCaUTAeZ1/kyd6/ArwG7DbGvFhr+STgX4HbjDFltYosBWaJSLiI3AikA994Kt5gUlVVxbvvvsvFixeZPXt2wCZUpRrSqVMn0tLSWL16NRUVFd4OR/khT57fGArcC4wRka3WYwrwMhALrLCW/RHAGLMTWADsAj4BHjXGVHsw3qBgjOHjjz8mPz+f22+/XYcgVEFNRBg3bhxlZWWsXbvW2+EoP+Sx07/GmDXUfZ102XXKPAc812xBKQBatWrFyJEj6d69u7dDUcrr2rZty4ABA4iNjfV2KMoPebyjkvItIsKwYcO8HYZSPmXq1KneDkH5Kf/q3qaUUkr5ME2qSimllJtoUlVKKaXcRJOqUkop5SaaVJVSSik30aSqlFJKuYkmVaWUUspNNKkqpZRSbiKBNNO9iJwC8hpZvDVw2o3h+AOtc3DQOgeHptS5ozGmjTuDCVYBlVSbQkQ2GmMyvR2HJ2mdg4PWOTgEY519kZ7+VUoppdxEk6pSSinlJppUvzXP2wF4gdY5OGidg0Mw1tnn6DVVpZRSyk20paqUUkq5iSZVQEQmicheEckVkX/zdjzuJCKHRWSHiGwVkY3WspYiskJE9lvPidZyEZHfWT+H7SLS37vRO0dEXheRQhHJqbXM5TqKyP3W+vtF5H5v1MVZ9dT5FyJyzDrWW0VkSq3PfmrVea+ITKy13C+++yLSXkRWishuEdkpIo9ZywP2OF+nzgF7nAOCMSaoH0AIcADoBLQAtgHdvR2XG+t3GGh91bLfAP9mvf434L+t11OA5YAAg4D13o7fyTqOAPoDOY2tI9ASOGg9J1qvE71dNxfr/AvgqTrW7W59r8OBG63ve4g/ffeBVKC/9ToW2GfVK2CP83XqHLDHORAe2lKFm4FcY8xBY8wlYD4wzcsxNbdpwN+s138Dbq+1/E3jsA5IEJFUbwToCmPMauDsVYtdreNEYIUx5qwxpghYAUxq/ugbp54612caMN8YU2GMOQTk4vje+8133xhTYIzZbL0uAXYDbQng43ydOtfH749zINCk6viSHq31Pp/rf3H9jQE+E5FNIjLXWpZsjCkAx39cIMlaHkg/C1frGCh1/4F1uvP1mlOhBFidRSQN6AesJ0iO81V1hiA4zv5Kk6rj9NDVAqlL9FBjTH9gMvCoiIy4zrqB/rOA+usYCHX/A9AZ6AsUAC9YywOmziISA7wPPG6MKb7eqnUsC5Q6B/xx9meaVB1/tbWv9b4dcNxLsWbWyokAAAavSURBVLidMea49VwILMJxKuhkzWld67nQWj2Qfhau1tHv626MOWmMqTbG2IE/4zjWECB1FpEwHMnlbWPMB9bigD7OddU50I+zv9OkChuAdBG5UURaALOApV6OyS1EJFpEYmteAxOAHBz1q+n1eD+wxHq9FLjP6jk5CDhfc2rND7lax0+BCSKSaJ1Om2At8xtXXf++A8exBkedZ4lIuIjcCKQD3+BH330REeA1YLcx5sVaHwXsca6vzoF8nAOCt3tK+cIDR0/BfTh6yD3t7XjcWK9OOHr6bQN21tQNaAV8Duy3nltaywV4xfo57AAyvV0HJ+v5Do7TYJU4/ip/qDF1BB7E0bkjF3jA2/VqRJ3fsuq0HccvzdRa6z9t1XkvMLnWcr/47gPDcJyy3A5stR5TAvk4X6fOAXucA+GhIyoppZRSbqKnf5VSSik30aSqlFJKuYkmVaWUUspNNKkqpZRSbqJJVSmllHITTapKKaWUm2hSVUoppdxEk6ryOBF52pofcrs1H2SWm7f/da3XP7Lmo3y79nJ3E5EEEXmkubbvLiJSetX7P4nIUBFJk1pzszZluyIyxZqrtENTYlXKH4V6OwAVXERkMDAVxzyRFSLSGsccj25jjBlS6+0jOEaWOeTOfdQhwdrXq1d/YA03J8YxVquvycIRd/uGVnSGiIwFfg9MMMYcccc2lfIn2lJVnpYKnDbGVAAYY04bY45bLaU9IvI3qwX7nohEAYjIPSLyjdWq/ZOIhFjL77PW3SYib9XsoKbVJCJ/xDFU41IR+fFVrak6y9ZW136tOHeLyJ+t1vZnIhIJ/BrobK37fK31XgU2A+1F5AkRybEej1v7uF69nxWRx2rF85yI/KiOOO8Xx9R+20XkS2cPhIh0A/YZY6qvWt5JRLaIyEAR+ZkV3woReUdEnrrO9objGOD9FmPMAWfjUCqgeHucRH0E1wOIwTGG6T4crbqR1vI0HOOcDrXevw48BXQDPgTCrOWvAvcBPXCMb9raWt6y1j5Ka70+XGudUuu53rK1ytW33zSgCuhrLV8A3GMtz6lVPg2wA4Os9wNwjNcabf0MduKYH7POetfaxmbrtQ3HuK2troozFtgFtLDeJzTw86/9s3kCeLDWvnKAm4AtOKYVy7SOVaS1n/01sdWx3Uock6b39vZ3TB/68OZDW6rKo4wxpTgSzP9v7w5erKziMI5/n2ghMgspatGiJBAhhCzSFunCWmRUYJRCFjQtksD+ghJ0YctWIpG5KAgqsEQQiiLMoRGhuk6LalFgIoigI4bORMrwtDjn4u1y37nXfJlkfD4w8N7zvuc958zm957znnt/24BzwKeSxuvp07Yn6/FHlB8Uf6Je/72kqfr5fuBx4IDt8/W+F66jG6PUbWoX4KTtqXr8IyUgDXLK9vF6vA44aHum/g8+B9bPM25s/wFMS3qIkk3lhO3pvjbmKEHvHUmP2L44dPTXPAl82fP5LkqWl5fr+NYBh2z/ZfsS5SGjyVXgGOWH/SNuWQmqseBcckF+a3sn8AbwfPdU/6WUbCMf2l5d/1ba3lXL/2s2iFHqNrUL8HfPdXM0702Y6btfk0Hj7toPjAOvUmax/77QngVWAZPAvlE3S9Ul5mWu+XarP4HTwGPD+ixpe13qnpJ0D2VWvgVYI+nNUfoQsRglqMaCkrRS0oqeotXAqXp8b93IBPAi8B0lndcLku6u9e+QdF8t3yLpzm75dXRjlLpN7Ta5RFkibTIBbJK0VCW37XNA9/3noHF3HQQ2AmsYkPdT0oo6+/0EOAwsmacPvTYAR/rKrgCbKHlIt9Z+PCtpiaQx4Onuhbb39jxwnKlls5RNaC9Jyow1bknZ/RsLbQzYI2kZ5d3k75Sl4DHgV+AVSe9R3t+9a3tW0g7gK0m3UZYZt9s+Lult4KikOcp7wPFROmD752F1bf8yqF3gbMM9pyVN1q+lfEHJ5dl7viPpA0rSaID9tk9IWj5o3D31rkg6Alx034ai6q0akGco72lfaxq3pNu5Nst+CjgwYBwzkp4BvgZ2U/J1/kR58PmBMpttZPuCpI3AhKTztg/Nd33EYpN8qnFTqMHlsO1V/3NXFtSwcdeA3gE22/7tBtt6EHjf9lpJHeBR21eH1BmzfbkuF08A22x3bqQfEYtZln8jblKSHqDM5L9pIaC+DnwM7ACw/fCwgFrtqxu1OsBnCagR88tMNSIioiWZqUZERLQkQTUiIqIlCaoREREtSVCNiIhoSYJqRERESxJUIyIiWpKgGhER0ZIE1YiIiJb8A6uwMkTGEhFBAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "s_max = 2500\n", "s_list = range(0, s_max)\n", "plt.plot(s_list, [params['T_bldg']]*s_max, label='$T_{bldg}$', color='red')\n", "plt.plot(s_list, [params['T_amb']]*s_max, label='$T_{amb}$', color='blue')\n", "plt.plot(saturation_curve['Ts']['s'], saturation_curve['Ts']['T'], color='gray')\n", "\n", "nb_points = 1000\n", "s12 = np.linspace(states.loc[1, 's'], states.loc[2, 's'], nb_points)\n", "s23 = np.linspace(states.loc[2, 's'], states.loc[3, 's'], nb_points)\n", "s34 = np.linspace(states.loc[3, 's'], states.loc[4, 's'], nb_points)\n", "s41 = np.linspace(states.loc[4, 's'], states.loc[1, 's'], nb_points)\n", "\n", "T12 = np.linspace(states.loc[1, 'T'], states.loc[2, 'T'], nb_points)\n", "T23 = np.linspace(states.loc[2, 'T'], states.loc[3, 'T'], nb_points)\n", "T34 = []\n", "for s in s34:\n", " T34.append(CP.PropsSI('T', 'S', s, 'P', states.loc[3, 'P'], params['WF']))\n", "T41 = np.linspace(states.loc[4, 'T'], states.loc[1, 'T'], nb_points)\n", "\n", "L_T = [T12, T23, T34, T41]\n", "L_s = [s12, s23, s34, s41]\n", "\n", "for i in range(len(L_s)):\n", " plt.plot(L_s[i],L_T[i],label='{}$\\\\to${}'.format(i+1,(i+1)%5+1))\n", "\n", "plt.legend(bbox_to_anchor=(1.05, 1), loc='upper left', borderaxespad=0)\n", "plt.xlabel('Specific entropy $s$ J/kg-K')\n", "plt.ylabel('Temperature $T$ K')\n", "plt.grid()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### P-h" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcUAAAEQCAYAAADBMm8YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8lOd97/3PT/suIRACxCIQqxDIGMmsdhCrMdgGAw5OvGXDrp02T5O053TJ6ZKkT5+euj5NDG1wcOs8jo0JTVgcGwwYsUlmEbvYMSCQQAgQIAmkQdJ1/piRq4BAM9LMfc/ye79e80Izuu+5vroY+OlerusSYwxKKaWUgjC7AyillFL+QouiUkop5aJFUSmllHLRoqiUUkq5aFFUSimlXLQoKqWUUi5aFJVSSikXLYpKKaWUixZFpZRSykWLolJKKeUSYXeAtnTr1s1kZmZa1l5dXR3x8fGWtRcItE/upX1yL+2Te9nVJyUlJVeMMWmdfI/uERERvwRyCN6DpmbgcGNj47dHjx59+e5v+mVRzMzMZM+ePZa1V1hYyKRJkyxrLxBon9xL++Re2if3sqtPRORcZ98jIiLilz169BiWlpZWHRYWFpQTYzc3N0tVVVX2pUuXfgk8dff3g/U3AaWUUp7LSUtLuxmsBREgLCzMpKWl3cB5NHzv9y3Oo5RSyn+FBXNBbOH6Gdusf1oUlVJKKRctikoppZSLFkWllFJ+ZcGCBZmpqam5gwYNGu7uPqdOnYocM2bM4AEDBgwfOHDg8B//+MfdO9K2FkWllFJ+5Zvf/OaVNWvWnPRkn8jISN54440LX3zxRenu3buPLlu2rHtJSUmMp2375ZAMpfyVw+Hg7NmzGBM49yLExsbSp08fRMTuKEq5ZebMmbXHjx+Pam+7nTt3xiYnJzcNHTrU0a9fvzv9+vW7A9ClS5fmrKys22VlZVGjR4+u96RtLYpKuamhoYF33nmHyspKu6N4LDc3l1mzZhEZGWl3FBUg/mzlgT4nLtXEefM9B/dIvPW/5+ee99b71dXVhb3wwgv916xZc2ro0KGOltePHz8edeTIkbivfOUrtZ6+pxZFpdxw6dIl9u7dC8D8+fNJTU21OZH7jh07xtatW6msrOSrX/0qKSkpdkdSymNLlixJffPNN3ve/XpVVVXkggULsg4dOnQU4MaNG2HPPPNM1j/+4z+eT01Nbfa0HZ8XRREZAnzY6qUBwP8yxvwfX7etlDecPHmSlStXIiJ84xvfID093e5IHunZsycZGRn89re/ZenSpcybN4+srCy7Yyk/580jOm947bXXrr322mvXWr928uTJqCeffHLgG2+8UQbQ0NAgs2bNylqwYMG1l1566XpH2vH5jTbGmOPGmIeMMQ8Bo4FbwO983a5S3lBSUsIHH3xAamoqo0aNCriC2GLw4MEsWrSIxMRE3nvvPbZt2xZQ10WVasvhw4djfv7zn5+bPn16XXNzMwsXLuw3ePDg+r/927/t8DUOq+8+nQKcNsZ0eo4+pXzJGMPGjRv56KOPGDhwIC+//DLR0dF2x+qU1NRUvvWtb5GTk8Nnn33GihUrcDgc7e+olMWefPLJ/hMnThx65syZ6PT09JFvvvlmt7a2mzt37s1p06bVAWzYsCFh1apVXbdv3544dOjQ7KFDh2Z/+OGHyZ62bfU1xYXABxa3qZTHCgsL2bFjB6NHj+aJJ54gLCw4Ri9FRUXxzDPPkJGRwaeffsqKFSt47rnnCA8PtzuaUl9au3btGU/3mTFjRq0xpqSzbYtVp1BEJAqoAIYbY+45tBWRRcAigPT09NHLly+3JBdAbW0tCQkJlrUXCEK5Ty5cuMDp06fp0aMHgwcP/nIoQ7D1ycWLFzlx4gRpaWkMGzasQ0M2gq1PvMGuPikoKCgxxuR15j0OHDhwNjc394q3MvmzAwcOdMvNzc28+3UrjxRnAnvbKogAxpilwFKAvLw8Y+XSK7r8zb1CtU8OHjzIli1bGDp0KAsWLPiDI8Rg7JMdO3awceNG+vbty6xZszwujMHYJ52lfRLYrCyKz6GnTpUFPvqTZ6iur+azZzI92u/27dtcvnyZ6MxoziSdYd2n6/7g+9evX+fdde96Mal/qB5cTWF5Ib9c/ku6dOni0b7B2ied0eE+uXWNoWGx/I/Zv4KIdsetKx+xpCiKSBwwDXjFivZUaEsuu0bknVse7dPQ0MDlqstERUXRPa17SM3+0qVLF5qbm7lx8wbh4eEkJSXZHSk03SiDOwbCdPi4nSzpfWPMLaCrFW0p1TexLwD/8fh/uLV9TU0Nb7/9NuHh4Xz7698mPj6+ze2C+bRYc3MzK1as4OTxk7z88sv06dPHrf2CuU86qkN9cnY7fD4LZr8JQXJTV6DS3lchrbGxkRUrVlBfX8/ChQvvWxCDXVhYGHPmzCEpKYmVK1dy65ZnR9qqk4regriukPuc3UlCnhZFFbKMMXz88cdcuHCBOXPmBOzAfG+JiYlhwYIF1NXVsWrVKh3cb5UrJ+HEJ5D/bYiMtTtNyNOiqELWnj172LdvHxMnTiQ7O9vuOH6hV69eTJ8+nZMnT7Jjxw6744SGz5dAeLSzKCrbaVFUIam8vJx169YxaNAgCgoK7I7jV/Lz8xk+fDifffYZFy9etDtOcKu7Cvvfh5HPQkKH1sQNOh1dLPjWrVsyYsSIYUOGDMkeOHDg8D/90z/t1ZH2tSiqkFNfX8/KlStJTExk7ty5QTNbjbeICLNnzyY+Pp41a9bQ1NRkd6TgtWcZNNbDuO/ancRvdHSx4JiYGLN9+/bjx48fP1JaWnpk06ZNSZs2bfL4JgH930CFFGMMa9eu5caNG8ybN4/YWL2G05aYmBieeOIJLl26RHFxsd1xgtOdeti1FAZOg+5D7U7jN/r163dn4sSJt+APFwtua9udO3fGHjt2LAqcN4slJyc3AzgcDmlsbJSODK3SATEqpJSUlHDkyBGmTJni9rCDUDVs2DCGDRtGYWEhQ4cOpVu3NudkVh116DdQVwXjXrc7SdtWvd6Hy0e8usgw3bNvMWex20tStbdY8N2LDDc2NpKTk5NdVlYW/dJLL12ePHlynacR9UhRhYzLly+zfv16srKymDBhgt1xAsITTzxBZGQka9eu1btRvckYKF4M6TkwYJLdafzS3YsFL1myJHXQoEHDWz9eeeWVzMuXL0ctWLAgCyAiIoJjx44dKSsrO7h379743bt3t3va9W56pKhCQmNjI7/97W+Jjo5mzpw5ITVjTWckJCQwbdo01q5dS2lpKTk5OXZHCg6nNkHVUZjz7+Cvn0UPjui8ra3Fgt1ZZLhFt27dmiZOnFizdu3a5Pz8/HpP2tYjRRUSCgsLqays5Mknn9RVHTw0atQoevTowYYNG7hz547dcYJD8c8hoQfkzLM7id/xZLHg1osMV1RURFy5ciUcoLa2VgoLC5OGDRvmUUEELYoqBJSVlVFUVMSoUaMYMmSI3XECjogwY8YMbt68qTfdeMOlw/BFIYxZpBN/t8GTxYJbLzJ8/vz5yEcffXTI4MGDs0eNGpVdUFBw87nnnrvhaft6+lQFtYaGBlatWkVycjIzZsywO07AyszMZNiwYWzfvp1Ro0aRmJhod6TAVbwYIuNg9DfsTuKXOrpY8JgxY24fPXr0SGfb1yNFFdQ2btxIdXU1c+bMITo62u44AW3atGk0NzezefNmu6MErpsXnXedjnoe4lLtTqPaoEVRBa0zZ86wZ88exo4dS79+/eyOE/C6dOnC6NGj2b9/P9euXWt/B3Wv3W9DcyOMedXuJOo+tCiqoNTsGqTfpUsXJk+ebHecoDFx4kTCw8PZunWr3VECj6MOdi+DobOga5bdadR9aFFUQel6dTXV1dU89dRTREZG2h0naCQmJpKfn8/Bgwe5cuWK3XECy/73of46jP9ju5OoB9CiqIJOfUMDN2tukp+fT2Zmpt1xgs6ECROIiIhgy5YtdkcJHM1NzhtsMvKgzxi706gH0KKogkpTUxPXrl4lIjyCKVOm2B0nKMXHx5Ofn09paSm3b9+2O05gOP4JVJ9xTunmr4P1FaBFUQWZ4uJiHHccpHZN1btNfWjs2LGEhYVx/rxtk54EluK3ILkvDHvK7iSqHVoUVdCorq5my5YtxMXFExfr3XmM1R9KTEwkNzeXS5cuUVvb5lzNqsWFEigrhrF/BOE6NNzfaVFUQcEYw+9//3vCwsJITdXxX1YYP348xhg+//xzu6P4t+KfQ3QSPPyC3UkCQmcXC25sbGTYsGHZBQUFAzvSviVFUURSRGSliBwTkaMiMs6KdlXoOHz4MKdPn2by5MlEhIfbHSckdO3alW7dulFSUqJzot5P9Tk4shpGvwTROguQOzq7WPBPfvKT9IEDB3b4YrdVR4r/CqwzxgwFcoGjFrWrQsDt27dZv349vXr1Ij8/3+44ISUjI4P6+noOHjxodxT/tPMXIGE6WN8DniwW3HqRYYDTp09Hrl+/Pvk73/lOh8cL+fwEt4gkAY8BLwMYYxyAw9ftqtCxceNGbt26xfPPP09YmF4RsFJycjI9evRg586dPPzww7okV2v1N2Dvr2D4XEjubXcaj/1ox4/6nKo+5dWL8wO7DLz14wk/bvfuLHcXC757keHXX3+9zz/90z9duHHjRodPF1lx1XcAUAX8h4jkAiXA94wxHq+IrNTdysrK2Lt3L+PGjaNHjx52xwk5IsIjjzzCmjVrOHPmDAMGDLA7kv8oeRccNc5hGMojLYsFX7lyJXzWrFlZu3fvjtm9e3fcm2++2fPubauqqiIXLFiQ9Zd/+ZcV3bp1a3z00UdvffTRRx0+Vy2+Xk1bRPKAz4EJxpidIvKvwE1jzI/u2m4RsAggPT199PLly32aq7Xa2lpdY+8ugdAnxhhKSkpobGwkPz+fcNe1xC5v/AsA1T/4vlfbC4Q+sVptbS1xcXEUFxeTkpLC8OHD7Y5ku9raWhLjYhiz8xXqY3qwf9RPLWm3oKCgxBiT15n3OHDgwNnc3Fy/mqroBz/4Qc/4+Pjmv//7v79nbcWWRYZ/9rOfnVu9enXKypUru0ZERJiGhoawurq6sOnTp19fvXr1mbbe98CBA91yc3Mz737diiPFC8AFY8xO1/OVwP+8eyNjzFJgKUBeXp6ZNGmSBdGcCgsLsbK9QBAIfbJr1y7q6upYsGAB2dnZX75+btk7AOR6OX8g9InVWvrE4XCwe/du8vPziY93+56IoFRYWMhXul6BhivEPPMWk4ZMsjtSQKmoqIiIiooy3bp1a2pZLPiHP/zhpba2bVlkeNq0aXXTp0+vW7x4cTnARx99lPjGG2+k368gPojPi6Ix5pKInBeRIcaY48AUoNNrXqnQVldXx+bNmxkwYADDhg2zO07IGzVqFDt37uTgwYOMGxfiN5cbA0U/h66DYJCu4emp8+fPR7788sv9m5qaMMbI008/fe1+iwXPnTv3prfbt2ok6R8DvxaRKOALQFfXVJ2yadMmHA4Hjz/+uN7c4QfS09PJyMhg3759jB07NqT/TpJvHIGL+2H2m6A3fnnMG4sFz549u2b27Nk1HdnXkr8xY8x+Y0yeMWakMWaOMabainZVcCovL2ffvn2MGTOGtLQ0u+Mol1GjRlFVVUV5ebndUWzV5/wqiE2FkQvtjqI6QH+NUQHFGMPHH39MQkICX/nKV+yOo1rJyckhMjKSffv22R3FPldO0fXqbsj/NkTpVIOBSIuiCij79u2joqKCadOm6YTffiY6Oprs7GwOHz6MwxGiQ5E/X4yRCHjkO3YnUR2kRVEFjNu3b7Nx40b69u3LiBEj7I6j2jBq1CgcDgdHjoTgvXR1V2H/B1SmfwUSutudRnWQFkUVMDZv3kx9fT0zZ84M6Rs5/Fnfvn1JTU3lwIEDdkex3p53oPE25/s8bXcS1QlaFFVAuHTpEnv27CEvL09nrvFjIkJOTg5nz56lpqZDN/8Fpjv1sGspDJzKrfi+dqdRnaBFUfm9lptrYmNjKSgosDuOakdOTg5AaJ1CPfQbqLsM475rdxLVSVoUld87dOgQ58+fZ8qUKcTGxtodR7UjLS2N7t27U1paancUaxgDxYshPQcGTLI7jeokLYrKrzkcDjZs2ECvXr0YNWqU3XGUm4YPH8758+e5caPNiUiCy+lNUHXUOfG3Xuv2mo4sFpyRkTFi8ODB2UOHDs3Oycnp0FRXWhSVX9u+fTu1tbU6c02ACalTqEVvQUIPyJlvd5Kg0tHFgrds2XLi2LFjRw4fPtyhdXu1KCq/df36dYqLixkxYgR9+vSxO47yQGpqKj179gz+U6iXDsMXm2HMIoiIan975RZ3Fwu+e5Fhb7Bq7lOlPLZx40YApkyZYnMS1RHDhw9n48aNVFdX06VLF7vj+MbnSyAyDkYH33TOFX/5V30aTp706rQ80YMG3er1Dz9td5FhdxcLvnuRYYApU6YMEhG+8Y1vVP3whz/0eBksPVJUfqmsrIzS0lLGjx9PcnKy3XFUB7SsrXjs2DGbk/hIzSU4uAIe+jrEpdqdJmh88MEHyS2LBbd+fcmSJamDBg0a3vrxyiuvZF6+fDlqwYIFWQA7duw4duTIkaOffvrpybfffrv7J5984vECqHqkqPyOMYb169eTmJjIhAkT7I6jOiglJYX09HSOHz8enMtJ7VoKzY0w9o/sTuIT7hzR+cL27dsTNmzYkJKRkZHcsljw008/3X/16tVnXnvttWutt21ZZPiNN94oA8jMzLwDkJGR0Thr1qzrxcXF8TNnzqz1pH09UlR+58CBA1RUVDB16lSiovQ6TSAbPHgwZWVl3Lp1q/2NA4mjDnYvg6GzoGuW3WmCyuLFi8srKysPlpeXH/rP//zPL8aOHVtzv8WCWxYZnj59et3NmzfDqqurwwBu3rwZtnnz5qSRI0d6fKOOHikqv+JwONi0aRMZGRk6v2kQGDJkCNu2bePkyZPk5ubaHcd79r8P9dd1sL7NWi8yfOHChYi5c+cOBGhqapJ58+ZdnT9/vseLEGtRVH6lZQjGs88+q0MwgkCvXr1ISEjgxIkTwVMUm5ucN9hkjIa+Y+1OE9Q8WSw4Ozvbcfz48U6PAdLTp8pv6BCM4CMiDBkyhFOnTtHY2Gh3HO84/glc+8J5lKi/uAUdLYrKb+gQjOA0ZMgQHA4HZ8+etTuKdxQvhuS+MOwpu5MoH9CiqPyCDsEIXv379ycyMpLjx4/bHaXzykugrAjGvgrhevUpGGlRVLYzxrBu3TodghGkIiIi6N+/P6dPn7Y7SucVvQXRSTDqBbuT+Epzc3Nz0J8Tdv2MzW19T4uist2BAwe4ePGiDsEIYllZWVRXV3Pt2rX2N/ZX18vgyGoY/RLEJNmdxlcOV1VVJQdzYWxubpaqqqpk4HBb37fk+F9EzgI1QBPQaIzJs6Jd5f90CEZoyMpyjuU7ffo0qakBOvvLzl84/xzzqr05fKixsfHbly5d+uWlS5dyCN6DpmbgcGNj47fb+qaVJ8ULjDEez0OngpsOwQgNqamppKSkcPr0afLz8+2O47n6G1DyLgyfC8m97U7jM6NHj74MhPQdRB4VRRHpAgwCYlpeM8Zs9XYoFRquX79OUVGRDsEIASJCVlYWhw4doqmpifDwB87z7H/2/gocNTBeB+sHOzHGuLehyLeB7wG9gf3AWKDYGDPZjX3PANWAAX5hjFnaxjaLgEUA6enpo5cvX+7uz9BptbW1JCR4PG9sULOiT44cOcLVq1fJz88nJiam/R3c1OWNfwGg+gff99p7gn5O2uJJn1RVVXHkyBFyc3NJSUnxcTLvkeZGxux8hfqYHuwf9dN2t7frc1JQUFCil6Y6z5Mjxe8B+cDnxpgCERkK/J2b+04wxlSISHdgg4gcu/sI01UolwLk5eWZSZMmeRCtcwoLC7GyvUDg6z4pKytjy5YtPPbYYxQUFHj1vc8teweAXC/n18/JvTzpk/r6eo4ePUpSUlJg9eOhldBwhZhn3mLSkEntbq6fk8DmyYXUemNMPYCIRBtjjgFD3NnRGFPh+vMy8DvgEU+DquChQzBCU0xMDBkZGYE1iN8YKH4Lug6EQTPsTqMs4ElRvCAiKcAqnEd7q4GK9nYSkXgRSWz5GpjOfW6FVaHh4MGDXLx4kSlTpugQjBDTr18/ysvLuXPnjt1R3HOuCCr2wdjXICxYb8ZUrbn9t2yMmWuMuW6M+VvgR8AyYI4bu6YD20XkALAL+L0xZl1HwqrA1zIEo1evXowcOdLuOMpi/fr1o7m5mfPnbVmqz3PFb0FsKuQ+Z3cSZZF2rymKSAzwKjAQOAQsM8ZscbcBY8wXQJBMj686a8eOHdTU1LBgwQIdghGC+vbti4hw7tw5BgwYYHecB7tyyjn592N/BlFxdqdRFnHnSPFdIA9nQZwJvOHTRCpo3bhxg6KiInJycnQIRoiKjo6mZ8+enDt3zu4o7ft8CYRHQn6bY7xVkHLn7tNsY8wIABFZhvMUqFIe01UwFDhPoe7atYvGxkYiIvx0Uu1b15wLCY98FhLT7U6jLOTOkeKXV8SNMUGyIJqy2vnz5zl8+DDjxo0LqDFqyvv69etHU1MTFy5csDvK/e1ZBo23nWsmqpDiTlHMFZGbIlIjIjXAyJavReSmrwOqwGeMYf369SQkJDBx4kS74yib9evXD8B/T6E2NsCut2HgVOg+zO40ymLtnrswxgTYfEzK3xw6dIjy8nKefvppHYKhiImJoUePHv5bFA/9BmorYdy/251E2cCdu0/XPOj7xpiQnjxWPZjD4WDjxo307NmT3Fy9CVk59evXj5KSEv+bB9UYKF4M6TkwwLszLanA4M5V7nHABeB9YCeg99ErtxUVFVFTU8O8efN0CIb6Ur9+/di5cycVFRX+dSfy6c/g8hGY82+gn9eQ5M41xR7AXwA5wL8C04ArxpgtnoxXVKHn5s2bFBUVkZ2d/eV1JKUAevd2Lr9UXl5uc5K7FL8FCT0gZ77dSZRN2i2KxpgmY8w6Y8xLOFfGOAUUisgf+zydCmibNm2iubmZqVOn2h1F+ZnExESSkpL8qyhWljqPFB/5DkTote9Q5dYgIRGJBmYBzwGZwM+A3/oulgp05eXlHDx4kAkTJtClSxe74yg/1Lt3b/8allG8GCLjIO+bdidRNnLnRpt3cZ46/QT4O2OMTuatHqhlFYz4+HgeffRRu+MoP5WRkcGRI0f8Y53KmktwcAWMfhniUu3NomzlzjXFF4DBONdTLHKNWbyp4xTV/ZSWlnLhwgUmT55MdHS03XGUn/Kr64q73obmRhj7R3YnUTZz55pimDEm0fVIavVINMYkWRFSBY47d+6wceNGevTowUMPPWR3HOXHevbsSVhYmP2nUB11zhlshs6Crln2ZlG20wXClFcVFxdz48YNZsyYQZiuP6ceIDIykvT0dPuPFPe/D7erdUo3BWhRVF5UU1PD9u3bGTp0KJmZmXbHUQEgIyOD8vJympub7QnQ3OxcDSNjNPQda08G5Ve0KCqv+eyzz2hqamLatGl2R1EBonfv3jgcDq5cuWJPgBOfwLUvnEeJOlhfoUVReUl5eTn79+9n7NixpKbq3XvKPRkZGQBUVFTYE6DoLUjuC8N0tkrl1G5RFJGXROSKiFwTkV+JSKIVwVTgMMbwySefEB8fz2OPPWZ3HBVAUlNTiYyM5OLFi9Y3Xl4CZUUw9lUI99N1HZXl3DlS/BHOqd2GAueAf/BpIhVwDhw4QHl5OVOnTtUhGMojYWFhpKenc+nSJesbL14M0Ukw6gXr21Z+y52ieNMYs88Yc9kY8yPgEV+HUoGjoaGBTZs2kZGRoatgqA7p0aMHly5dwhhjXaPXz0PpKhj9EsToyDL139wpij1FZJGIPCoiaUBkRxoSkXAR2SciH3Vkf+Wftm7dSm1tLTNnztRVMFSH9OzZE4fDQXV1tXWN7nStlfjIK9a1qQKCOyfS/wYYCXwdGAEkiMjHwAHgoDHmAzfb+h5wFNBfy4LE1atX+fzzz3nooYe+vGFCKU/16NEDgIsXL1pzk1b9TSh5F4bPhRQ/WrZK+QV3ZrRZaoz5rjHmK8aYVGAA8BZwHXjCnUZEpDfOCcV/2Zmwyr+sX7+eiIgIpkyZYncUFcC6d+9OWFiYddcV9/4KHDUw7nVr2lMBxeNbrowxF3AuOvyxB7v9H+DPAb1zNUicOHGCkydPMm3aNPsnc1YBLSIigrS0NGuKYlOj89RpvwmQ8bDv21MBx+f3IYvIbOCyMaZERCY9YLtFwCKA9PR0CgsLfR3tS7W1tZa2Fwge1CfNzc3s2bOH2NhY6uvr/a7vuly/DsAZL+fSz8m9vNkn586d83n/dq/cSvaN8xzq8wJXfdSWfk4CmxWDcyYAT4nIE0AMkCQi7xljnm+9kTFmKbAUIC8vz0yaNMmCaE6FhYVY2V4geFCf7Nixg9u3b/P1r3+dgQMHWhvMDeeWvQNArpf/TvVzci9v9UlMTAzr168nLy/Pd2cejIG3/xa6DmTEM38GPpqbVz8ngc3tT4U4PS8i/8v1vK+ItDs8wxjzF8aY3saYTGAh8NndBVEFjpqaGrZu3crgwYP9siCqwNS9e3cALl++7LtGyoqhYh+Mfc1nBVEFPk8+GUuAccBzruc1wGKvJ1J+bdOmTTQ1NTFjxgy7o6gg0lIUq6qqfNdI0VsQmwq5z7W/rQpZnhTFMcaY14F6AGNMNRDlSWPGmEJjzGxP9lH+o6ysjAMHDuj8psrr4uPjiY2N9d2R4tXTcPxjyP8WRMX5pg0VFDwpindEJBwwAK6B/Dat96Ks1tzczMcff0xSUpLOb6q8TkTo3r27744UP18C4ZGQ/x3fvL8KGp4UxZ8BvwO6i8hPge3oPKghY9euXVRWVjJjxgyiojw6QaCUW9LS0rh8+bL3p3u7dQ32/RpGPguJ6d59bxV03Lr7VJzzd20FSoApgABzjDFHfZhN+YmamhoKCwvJyspi2LBhdsdRQap79+40NDRQU1NDUpIXJ77aswwab8NYHayv2udWUTTGGBFZZYwZDRzzcSblZzZs2EBjY6POb6p8Ki0tDXDegeq1otjYALvehqwpkJ4cJgUZAAAbO0lEQVTtnfdUQc2T06efi0i+z5Iov3T27FkOHTrE+PHj6dq1q91xVBDzyR2oh1ZCbSWM/6733lMFNU8G7xcAr4jIOaAO5ylUY4wZ6ZNkynZNTU18/PHHpKSk8Oijj9odRwW5uLg44uLivFcUjXGumdh9OAwo8M57qqDnSVGc6bMUyi/t3LmTqqoqFi5cSGRkh1YMU8ojXbt25dq1a955s9OfweVSeHoJ6Gl/5Sa3i6Ix5pwvgyj/0tDQQFFREYMHD2bIkCF2x1EhomvXrpw6dco7b1a8GBLSYcR877yfCgntXlMUke2uP2tE5KbrUdPy3PcRlR1Onz6NMYbHH3/c7igqhKSmplJbW4vD4ejcG1UegdOb4JFFEBHtnXAqJLiznuJE15+Jxpgk1yOx5bnvIyqrnT59mqqqKiZOnEiXLl3sjqNCSMvNXFevXu3cGxUvhsg4yPumF1KpUOLJhOALRCTR9fVfi8hvRWSU76IpOzQ2NvLJJ58QExPDhAkT7I6jQkzL9IGduq5YUwmHVsBDX4M4nY5QecaTIRk/MsbUiMhEYAbwLvDvvoml7LJt2zauXr3KoEGDiIiwYmUxpf5bS1Hs1JHi7reh6Y5zNQylPORJUWxy/TkL+DdjzGo8nBBc+beqqiq2b9/OiBEjdMJvZYuoqCgSExM7fqTouAW7l8HQWdA1y7vhVEjwpCiWi8gvgGeBj0Uk2sP9lR8zxrB27Vqio6N1WShlq65du3b8SPHA+3D7GozTwfqqYzwpas8C64HHjTHXgS7An/kklbLc3r17OX/+PNOmTSM+Pt7uOCqEpaamdqwoNjdD8RLo9TD0Hev9YCokeFIUZwEbjDEnReSvcS46fMU3sZSVampq2LBhA5mZmTz00EN2x1EhrmvXrty+fZvbt297tuOJdXDttHNKNx2srzqoszfa/JtvYikrrV+/nsbGRmbPnq0TfivbtVzPrq6u9mzH4rcguS8Me9oHqVSo0BttQtyJEycoLS3lscce0wm/lV9ITk4G4Pr16+7vVL4Xzu2Asa9CuN41rTquIzfafBW90SYoOBwOPv74Y9LS0nRMovIbKSkpgIdFsfgtiE6CUS/4KJUKFR250WaG60abVPRGm4C2efNmbty4wezZswkPD7c7jlIAxMTEEBUVxY0bN9zb4fp5KF0FD78IMTrJluocT4ribSAeeM71PBLw4Fc55U8qKirYuXMno0ePpm/fvnbHUepLIkJKSor7RXGnaw6RMa/6LpQKGZ4UxSXAWP67KNYAi9vbSURiRGSXiBwQkVIR+bsO5FRe1NTUxJo1a4iPj2fq1Kl2x1HqHikpKe6dPq2/CXt/BcPnQkof3wdTQc+TojjGGPM6UA9gjKnGvRttGoDJxphc4CHgcRHRQUQ22rZtG5WVlcyePZuYmBi74yh1j+TkZPeK4t5fQcNNGPe670OpkOBJUbwjIuGAARCRNKC5vZ2MU63raaTrYTwNqrzj0qVLbNu2jREjRug6icpvJScn09DQQH19/f03amp0njrtNwEyHrYunApqnhTFnwG/A7qLyE+B7cA/uLOjiISLyH7gMs4JAHZ6nFR1WlNTE6tWrSI2NlbXSVR+reUO1AdeVzy6Gm6c1yndlFeJMe0ftIlzRHdvnDfaTAEE2GSMOepRYyIpOAvrHxtjDt/1vUXAIoD09PTRy5cv9+StO6W2tpaEhATL2rPL2bNnOXfuHMOHD6dbt24P3DaQ+6TLG/8CQPUPvu/V9w3kPvEVX/XJzZs32bdv3/0/q8bw8N4/I6Kxjl2PLAbxn9Fhdn1OCgoKSowxeZY3HGTcGuVqjDEissoYMxo41tHGjDHXRaQQeBw4fNf3lgJLAfLy8sykSZM62ozHCgsLsbI9O1RWVrJt2zZycnKYN29eu9sHcp+cW/YOALlezh/IfeIrvuqT2tpa9u3bR9++fXnkkUfu3eBcMWw5CbPeYFL+ZK+33xn6OQlsnvx69bmI5HvagIikuY4QEZFYYCqdKKzKc01NTaxevZrY2Fhmzpxpdxyl2hUXF4eIUFtb2/YGxW9BbCrkfs3aYCroeTIfUgHwqoicBepwnkI1xpiR7ezXE3jXdZNOGLDCGPNRR8KqjtmxYwcXL15kwYIFxMXF2R1HqXaFhYURHx9PTU3Nvd+8ehqO/R4e+yFE6edZeZcnRbFDhxjGmIPAqI7sqzrv8uXLbNmyheHDh5OdnW13HKXclpiY2PaR4udLIDwS8r9jfSgV9NotiiISA7wKDAQOAcuMMY2+DqY6r+W0aUxMjJ42VQEnMTHx3rtPb12Dfb+GEc9CYro9wVRQc+ea4rtAHs6COBN4w6eJlNds3bqViooKZs2apQsHq4CTkJBw75Hinneg8bYO1lc+487p02xjzAgAEVkG7PJtJOUNFy5cYNu2bYwcOVJPm6qAlJiYSF1dHU1NTc4J6xsbYNdSyJoC6fqZVr7hzpHinZYv9LRpYHA4HPzud78jKSlJT5uqgNUy1q+urs75wqGVUFupR4nKp9wpirkictP1qAFGtnwtIjd9HVB57tNPP+XatWvMmTNH5zZVASsxMRHAeQeqMVC8GLoPhyz/Gpeogku7p0+NMbrQXgA5ceIEJSUljBs3jszMTLvjKNVhLUeKtbW18MVmuFwKTy8BEZuTqWDmyZAM5efq6upYs2YN3bt3Z/Lk0P1t+uzVOm45mvjzXxR79X2vX7/Nvx337nsGOl/2SWNjI+UNQ9j1URk/dfy/ZEV0JXHEfJ+0pVQL/5kwUHWKMYa1a9dSX1/PM888Q0RE6P6+c71Xfy6l6dp6gS483PnfU0pjFTkNB9jdfR5ERNucSgW70P2fM8js2bOH48ePM336dNLTQ3v81tO//GcAvuHl93XOaTnOy+8a2HzZJ8YYfvKTnzBu9DgiHjnE5EidvUb5nhbFIFBZWcn69esZOHAgY8fq+s0qOIgIcXFx3Lp1C5J62R1HhQg9fRrg7ty5w8qVK4mNjWXOnDmI3oSggkhcXBy3b9+2O4YKIVoUA9y6deu4cuUKc+fO1VlrVND58khRKYtoUQxgpaWl7N27lwkTJjBgwAC74yjldbGxsVoUlaW0KAao69evs3btWnr37k1BQYHdcZTyidjYWD19qiylRTEANTY28pvf/AaAZ555xjkvpFJBSIuispoWxQC0fv16KioqmDt3Ll26dLE7jlI+Ex0dTXNzM42NOu2ysoYWxQBz6NAh9uzZw/jx4xkyZIjdcZTyqeho52D9hoYGm5OoUKFFMYBUVVWxdu1a+vbty5QpU+yOo5TPRUVFAc6VX5SyghbFAOFwOFixYgVRUVHMnz+fsDD9q1PBT48UldX0f9YA0DKv6dWrV5k3b96XS+ooFey0KCqraVEMAMXFxRw+fJiCggL69+9vdxylLKNFUVnN50VRRPqIyGYROSoipSLyPV+3GUxOnz7Nxo0byc7OZuLEiXbHUcpSLdcUtSgqq1gxIXgj8ANjzF4RSQRKRGSDMeaIBW0HtGvXrrFy5UrS0tJ4+umndV5TFXIiIyMBdEiGsozPjxSNMReNMXtdX9cAR4EMX7cb6BwOBx9++CEACxcu/PI3ZqVCScvEFE1NTTYnUaFCjDHWNSaSCWwFcowxN+/63iJgEUB6evro5cuXW5artraWhIQEy9prjzGGI0eOcOXKFUaMGEFqaqrlGfytT/yB9sm9fN0nd+7coaioiKysLHr37u2zdrzJrs9JQUFBiTEmz/KGg4xl6ymKSALwX8D/c3dBBDDGLAWWAuTl5ZlJkyZZFc21UKp17bWnsLCQK1euMG3aNMaPH29bBn/qE3+gfXIvX/eJw+GgqKiI/v37M2HCBJ+14036OQlsltx9KiKROAvir40xv7WizUB16NAhtmzZQm5uLuPG6SrvKrTp6VNlNSvuPhVgGXDUGPMvvm4vkJWVlbF69Wr69evHk08+qTfWqJDXMkmFFkVlFSuOFCcALwCTRWS/6/GEBe0GlOrqaj788EOSk5N59tlndeULpQARITw8XIuisozPrykaY7YDesjzAPX19bz//vsYY/ja175GXFyc3ZGU8htaFJWVdEYbmzU2NrJixQquXbvGs88+S9euXe2OpJRfCQ8P13GKyjJaFG1kjGHVqlWcOXOGp556iszMTLsjKeV39Nq6spIWRZsYY1i3bh2lpaVMnTqV3NxcuyMp5ZeMMVoYlWW0KNpkx44d7Nq1i7Fjx9o2FlGpQKBFUVlJi6IN9u/fz6ZNmxgxYgTTp0/Xf/BKPYAWRWUlLYoWO3bsGGvWrGHAgAE6ybdSbtCiqKykRdFCp06dYuXKlfTq1UvHIirlJi2KykpaFC1y9uxZPvzwQ9LS0vj617/+5eKpSqkH06KorKRF0QIXLlzggw8+ICUlheeff57Y2Fi7IykVMJqbm7UoKstoUfSxixcv8utf/5r4+HhefPFF4uPj7Y6kVMBobm6mubn5y8WGlfI1y5aOCkXl5eW89957REdH8+KLL5KYmGh3JKUCisPhANBFtpVltCj6yIULF3jvvfeIjY3lpZdeIiUlxe5ISgUcLYrKaloUfaCsrOzLU6YvvfQSycnJdkdSKiBpUVRW06LoZWfPnuX9998nKSmJF198kaSkJLsjKRWw7ty5A2hRVNbRouhFx48fZ+XKlaSkpOg1RKW8QI8UldW0KHrJ/v37WbNmDT179uRrX/ua3mWqlBc0NDQA6N2nyjJaFL2gqKiIDRs20L9/f7761a/qwHylvOTWrVsAuvC2sowWxU4wxrBx40aKiorIzs5m7ty5RERolyrlLXV1dQAkJCTYnESFCv0fvIPu3LnD6tWrKS0tZfTo0TzxxBOEhelcCEp5U21tLREREXr6VFlGi2IH1NXVsXz5ci5cuMDUqVMZP368TkOllA/U1dWRkJCg/76UZXxeFEXkHWA2cNkYk+Pr9nytqqqK999/n9raWhYsWEB2drbdkZQKWnV1dXrTmrKUFef7/hN43IJ2fO7UqVMsW7aMO3fu8PLLL2tBVMrHtCgqq/n8SNEYs1VEMn3dji8ZY9i+fTufffYZ6enpLFy4UKdtU8rHjDFcv36dvn372h1FhRC9ptiOhoYGVq1axbFjxxgxYgRPPvmkXvRXygL19fU0NDTQpUsXu6OoECLGGN834jxS/OhB1xRFZBGwCCA9PX308uXLfZ6rRW1tbZu3fNfV1XHkyBFu3bpFVlYWGRkZIXPB/359Esq0T+7lyz65efMm+/btY/jw4XTr1s0nbfiCXZ+TgoKCEmNMnuUNBxm/OVI0xiwFlgLk5eWZSZMmWdZ2YWEhrdszxrB//36KioqIjIzkxRdfpH///pbl8Qd394nSPmmLL/vk8OHD7Nu3j0cffZT09HSftOEL+jkJbH5TFP1FQ0MDv//97zl06BD9+/dn7ty5OoepUja4evUqgJ4+VZayYkjGB8AkoJuIXAD+xhizzNftdkRFRQX/9V//RXV1NQUFBUycOFEH5Ctlk8rKSlJTU3UycGUpK+4+fc7XbXRWc3MzmzdvZtu2bSQmJvLyyy/rHW9K2ayysjKgTpuq4BDyp08rKyvZt28ftbW15Obm8vjjjxMTE2N3LKVCmsPh4Nq1a4wcOdLuKCrEhGxRbGpqoqioiC1bthAWFsZXv/pVhg4dancspRRw8eJFAHr27GlzEhVqQrIolpWV8dFHH1FVVUV2djbJyclaEJXyI+fOnQOgT58+NidRoSakiuLt27fZuHEje/fuJSkpiYULFzJkyBAKCwvtjqaUauX8+fOkpaURGxtrdxQVYkKiKDY3N7N37142b97M7du3GTt2LAUFBXpXm1J+qKmpifPnz5OTE/DrB6gAFPRF8cyZM6xfv57Kykr69u3L448/rtcplPJjZWVlNDQ0kJWVZXcUFYKCtihWVVXx2WefcezYMZKTk5k/fz7Z2dkhM02bUoHqxIkThIeHa1FUtgi6olhdXc2WLVs4ePAgkZGRFBQUMG7cOJ3EW6kAYIzh+PHj9O/fXy9vKFsETVF0OBxs3LiRkpISwsLCGDt2LBMnTiQuLs7uaEopN1VUVFBdXc2jjz5qdxQVooKmKEZERHDu3DlGjRrFY489RlJSkt2RlFIeqq6uJiEhgWHDhtkdRYWooCmKYWFhLFq0iPDwcLujKKU6KCcnh+zsbJ1zWNkmqD55WhCVCnxaEJWd9NOnlFJKuWhRVEoppVy0KCqllFIuWhSVUkopFy2KSimllIsWRaWUUspFi6JSSinlIsYYuzPcQ0SqgHMWNtkNuGJhe4FA++Re2if30j65l1190s8Yk2ZDu0HFL4ui1URkjzEmz+4c/kT75F7aJ/fSPrmX9klg09OnSimllIsWRaWUUspFi6LTUrsD+CHtk3tpn9xL++Re2icBTK8pKqWUUi56pKiUUkq5hFRRFJHHReS4iJwSkf/ZxvcfE5G9ItIoIvPtyGg1N/rk+yJyREQOisgmEelnR04rudEnr4rIIRHZLyLbRSTbjpxWaq9PWm03X0SMiAT93ZdufE5eFpEq1+dkv4h8246cykPGmJB4AOHAaWAAEAUcALLv2iYTGAn8Cphvd2Y/6ZMCIM719R8BH9qd2w/6JKnV108B6+zObXefuLZLBLYCnwN5due2u0+Al4G37M6qD88eoXSk+AhwyhjzhTHGASwHnm69gTHmrDHmINBsR0AbuNMnm40xt1xPPwd6W5zRau70yc1WT+OBYL8w326fuPwY+Ceg3spwNnG3T1SACaWimAGcb/X8guu1UOZpn3wL+MSnieznVp+IyOsichpnEfgTi7LZpd0+EZFRQB9jzEdWBrORu/925rkuPawUkT7WRFOdEUpFUdp4Ldh/w2+P230iIs8DecD/9mki+7nVJ8aYxcaYLOB/AH/t81T2emCfiEgY8CbwA8sS2c+dz8laINMYMxLYCLzr81Sq00KpKF4AWv+m1huosCmLv3CrT0RkKvBXwFPGmAaLstnF08/JcmCOTxPZr70+SQRygEIROQuMBdYE+c027X5OjDFXW/17eRsYbVE21QmhVBR3A4NEpL+IRAELgTU2Z7Jbu33iOi32C5wF8bINGa3mTp8MavV0FnDSwnx2eGCfGGNuGGO6GWMyjTGZOK89P2WM2WNPXEu48znp2erpU8BRC/OpDoqwO4BVjDGNIvJdYD3OO8feMcaUisjfA3uMMWtEJB/4HdAFeFJE/s4YM9zG2D7lTp/gPF2aAPxGRADKjDFP2Rbax9zsk++6jp7vANXAS/Yl9j03+ySkuNknfyIiTwGNwDWcd6MqP6cz2iillFIuoXT6VCmllHogLYpKKaWUixZFpZRSykWLolJKKeWiRVEppZRy0aKolFJKuWhRVEoppVxCZvC+Uso3XHOf/hhIwjlwXef4VAFLjxTVHxCRvxKRUtfM/vtFZIyX37+o1dd/IiJHReTXrV/3NhFJEZHXWj3PFJHDHr5HrRfzTBWR/7+jbYvIL0RkQkd+Djfbe0dELt/93g9YVPdpnCtE3ME5J6hSAUuLovqSiIwDZgMPu2b2n8ofLo/TacaY8a2evgY8YYz5+l2ve1uKqy1/kQvs68T+Y3DOL+or/wk83voFEQkHFgMzgWzgORHJdn17CFBsjPk+zoWolQpYWhRVaz2BKy0z+xtjrhhjKlxHJMdE5N1Wa8PFgXNJKRHZ5Tqq/IXrP09E5EXXtgdaHxW1HPWIyL/jXLV8jYj8aeujofvt21pb7bpyHhWRt11Hu5+KSCzwj0CWa9uWpa/C29gOEVklIiWu1xe10e6D+uLHIvK9Vtv+VETaWmsxF+ghIttE5JJrHlW3iMgw4IQxpumu1weIyD7X/L2IyI9cOTeIyAci8sO7ts935Y8RkXjXz5sDYIzZinOuztYetKjuBZxzwAI0oVQgM8boQx8YY8A58fd+4ASwBPiK6/VMnGvFTXA9fwf4ITAM55pxka7XlwAvAsOB40A31+uprdqobfX12Vbb1Lr+vO++rfa7X7uZOCdffsj1+grgedfrh1vt3+Z2rdsDYoHDQNe78rXZF62+t9f1dRhwumX/u/IfAP7c9fUzwH+08/fSus++D3yzVXuHcR6p7Wv18+S5/h5jcS7rdLIl413v+xPgn3EeAf7FXd+7u8/mA79s9fwF4C3X13HAMuDnwOt2f471oY/OPPRGG/UlY0ytiIwGHgUKgA9d144KgfPGmB2uTd/Dudp8Pc414na7VtCIBS4DycBKY8wV1/vefdTxIJPd2HfKfdrdCpwxxux3bVeC8z/37W28R1vbgXNlg7mur/sAg4Crd+3bVl/8szHmrIhcFedyW+nAPmPMH+wrIpFAKs5iBM6b3a63ke9+ZgDfaPU8DVgNzDPGlLpemwisNsbcdrW59j7v9fc4l0Cqd/0MD3LfRXWNMbeAb7mVXik/p0VR/QHjPC1XiHPB2EM4l0Uq5N5VxQ3O/yjfNcb8RetvuE4ZdnT5FXFj3/u1mwm0XgS5CWfBbMs924nIJJzXUccZY26JSCEQ08a+bfVFi1/iXCKoB86jyLtlAweMMc2u5yNxHu21y3WaNsUY03ox2xs4r/tOAFqKYlsFrC2pOM8OROL8OesesK0u0q1Cgl5TVF8SkSHyhwvoPgScc33d13UjDsBzOI++NgHzRaS7a/9UEennev1ZEena8roHMdzZ937t3k8NztOI7UkGql0FcSjOFeTb0lZftPgdzptU8nGutXe3XJynT1uMBA66kQ2cR++b73rNAcwBXhSRr7le245zPdAYEUnAuRByW5YCPwJ+Dfx/7bSti3SrkKBHiqq1BODnIpKC85rbKWCR6/WjwEsi8guc16j+zVU8/hr4VJxj1e7gvKb0uYj8FNgiIk04r3e97E4A41yo9YH7GmOOtNUucOk+73lVRHaIc4jBJzivobVlHfCqiBzEeV3zfnd43tMXrdpyiMhm4Lq562YYl1xgV6vnOTzgSFFEIvjvo9qZwMo2fr46EZkNbBCROmPMahFZg7P4ngP24DyibP2+LwKNxpj3XTdHFYnIZGPMZyLyATAJ6CYiF4C/McYskzYW1b1fbqUClS4yrNrlOi35kTEmx+YotmuvL1xFei+wwBhz0gvt5QJvG2MeEZG9wBhjzB039ktwXSOOw3mtdZExZm9n8ygV7PT0qVJeIs5xe6eATV4qiK8CHwB/DWCMedidguiyVET24yzQ/6UFUSn36JGiUkop5aJHikoppZSLFkWllFLKRYuiUkop5aJFUSmllHLRoqiUUkq5aFFUSimlXLQoKqWUUi5aFJVSSikXLYpKKaWUy/8FpuI9c1ZLnCIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(saturation_curve['ph']['h'], saturation_curve['ph']['p'], color='gray')\n", "\n", "nb_points = 1000\n", "h12 = np.linspace(states.loc[1, 'h'], states.loc[2, 'h'], nb_points)\n", "h23 = np.linspace(states.loc[2, 'h'], states.loc[3, 'h'], nb_points)\n", "h34 = np.linspace(states.loc[3, 'h'], states.loc[4, 'h'], nb_points)\n", "h41 = np.linspace(states.loc[4, 'h'], states.loc[1, 'h'], nb_points)\n", "\n", "P12 = np.linspace(states.loc[1, 'P'], states.loc[2, 'P'], nb_points)\n", "P23 = np.linspace(states.loc[2, 'P'], states.loc[3, 'P'], nb_points)\n", "P34 = np.linspace(states.loc[3, 'P'], states.loc[4, 'P'], nb_points)\n", "P41 = np.linspace(states.loc[4, 'P'], states.loc[1, 'P'], nb_points)\n", "\n", "L_h = [h12, h23, h34, h41]\n", "L_P = [P12, P23, P34, P41]\n", "\n", "for i in range(len(L_h)):\n", " plt.plot(L_h[i]/1e6, L_P[i]/1e6,label='{}$\\\\to${}'.format(i+1,(i+1)%5+1))\n", "\n", "plt.legend(bbox_to_anchor=(1.05, 1), loc='upper left', borderaxespad=0)\n", "plt.xlabel(\"Specific enthalpy $h$ J/kg x$10^6$\")\n", "plt.ylabel('Pressire $P$ MPa')\n", "#plt.yscale(\"log\")\n", "plt.grid()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Perfoemance of heat pump cycle\n", "### Power required by the compressor\n", "\\begin{align}\n", "\\eta_{vol} = 1-c(\\frac{v_2}{v_3}-1) \\\\\n", "\\dot{m} = \\frac{\\dot{V}_{disp}\\eta_{vol}}{v_2} \\\\\n", "\\dot{W}_c = \\frac{\\dot{W}_c}{\\dot{m}}\\dot{m} \\\\\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIgAAABlCAYAAACWYKsOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAACxhJREFUeJztnXuwVVUdxz9fXr4uYZpaGmhKiqYCZphO4kUjbVLRMcVxiChNybDSrNRR09HSfCXJJOKYF50xH5OpmI74Ih+Y+EAFDEokAydxFNMARaRff/zWgcXm3nXOuZd7zwbXZ2bPvft31mvv/d17rX3O+u4tMyOTaYtujW5AptxkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIEUkHSqpAWSPpD0nKQDq6TvLelqSa9Jel/SdElfKqTpLumiqNwFki6W1KNzt2Y9YGZ5CQswElgJfA/YHbgGWAr0S+S5Dfgb0Az0By4A3gV2iNKcAywBjgB2Ao4E3gHOa/Q2V90njW5A2IEvABeGA/IW8GYjdh7wNHB9IfYP4JI20m8GfASMKMSfAy6O1u8FJhfSTAbubfS+r7Y0vIuR1BM/W0cBM4AvA9cCF0r6dAfLPkfS0irLgSFtL+CLwNRCMVOBA9qoogfQHfigEH8f+Eq0/gQwTNKAUNcewMHAfR3Zvi6h0QoFBgIGHBXF+obYgHaUdzgwIfy/FX7ZTy2bhbTbhzqHFso7H5iXqG868DiwAy6WUcCqOA8g4JfA//AuzIiuMGVeyjBIGoT3x1Oi2Dbh7xvtKG9v4EUAM1uC9/31UPx5W63EYr4F/B5YhAvjeeAPwD5RmpHAaOAEYA6+zeMlLTCzG+psX5fS8C4Gv4I8a2arotgg4DVggKT7K0FJwyXdHv7/QrhjmC1piqTNo/JeCmlq7mLwsc8qoNitbQssbqvxZjbfzA4CmoC+ZjYE6AksiJJdDlxhZrea2Swzuxm4Cji7jv3UEMpyBXmmEBuMD1xfBnaN4ucBYyVtCtwOjDKzmZKuBMYAvwP2BGaH9BNDuhSvA5jZh5KeA4YDd0SfDwf+WG0jzGwZsEzSJ4FDgZ9FH2+Oiy9mFeU4QZOUQSB7A5MKscHAg2b2nqQeQRBDgUVm9rKk44GHzGxmSD8X2D6kUzhY7elirgJuljQDeBIYi49NJlYSSBoHjDOzyoDzUPxAz8XHNJcD84Abo3KnAGdJWoB3MYOBM4Cb6mhbY2jwAPWzFAajeJ//X+DIsD4V2AsfCO4WYhcCJ0V5fgscC+wL3NHBNp0K/BNYgd+uFgetF/huW71+HDA/pP83MAHoU8jTG7ga7zbfB14FfgVs2uhBaLVFYQNKg6Rd8TOwn5ktlHQ13r/3MrPRIc0pQH8z+6mkvYFb8bNyFLC9mV3UoOZvdJShiykyGHjbzBaG9Tn42GL3KM3NwG2SZgHvAceZ2Yoglke6tLUbOWW8glwK7GtmX210WzIlFEimXJT+NivTWLJAMkmyQDJJskAySbJAMkmyQDJJskAySbJAMkmyQDJJskAySbJAMklKKRBJ/cJ0wMGNbsvHnVL+WBccZzsBC81sRYOb87GmlALZmJE0DZhtZuMa3ZZaKGUX0xm0w3M7VNI9kl6XZJLGtJLmB5JekvReWJ6S9I1O24h1658o6TedWcfHQiCSRgLj8Xmgg3Gz0/2S+iWyNeGz43+EzyNtjUXAz3EPzL74bLa7wsy2TkWScK/v3Z1aUaMnxbJm4vJIfAcvx41PA/Cd/liIzSBhoq5SR12e21byLwXG1Jh2CXBK4vNpBOdfWD8E+E+cB9gCn/G+FPfknI37e1uiNEOAt4EeYX0o8NeQ592wzXt29PiU4QoyKPz9Pj5bfT9gE6AFuAx3xu8PfAq3CnSF57ZuwiMejsevPNNrzHMM8CfgZDO7LvroSuAg4GjcwzsQKHaJRwF/NrOPwqD+btwDPBDfh+NZ14tTN2WYtDwQV/xIM1sMIOlB4HjcDvF2iP0F+EzIU7MhChdWd9Z1xy0GOjzvVdJewFPApvjZe7SZzaoh38m4h+abZjY1ijcB3wVGm9mDIXYi3p3FjMCNZACfALYEppjZ/BCb2+6NiiiDQAbhj0GID2A/4M6KOKLYM9BlnttamYdvw5bAMcBkSc1mNjuRZwRwCu65earw2S64dXNGJWBmyyStLk9Sf2Bn4IHw+RJJLcADkh4GHsb9QQvpIGXoYgbifWfMYPysLKabCV3jua0VM/vQzF4xs2fN7GzcMnp6lWwv4SarE8NgM6aynhLvUcDDFhyEoR3fwbuWx/AH1Pw9uP46REMFImkL/IyZGcW2wh//EMf6AltHsYn4WZtangU/gLhDbnih+uHUOFaok274GCrFAvyJRF8DJhVE8gr+iIghlUAwpu8ZpRkB3FUs1MxeNLNfm1kzPhj+dv3NX7fQRt7B7I+f3U2FUf0HQM8odgRux+zWznpGAh8CJ+EGrPH4eGHHKM04YG603sQasS3HnxMyiOhOCrgUHzzuhNtDL8GfAfL1Wu5i8JNjIXA94UvLEJ+I2z8PAfYAbsHHaTfij8ZYCWwXpf9caMsBwI7AMHwMdm6Hj1GDBTI2Pighdib+OIg4di7wRAfrqtdz24xf5otLS5SmBffbrsAfm/UQcGiVdqwWSEEkk1jzzXYT7h5chneDZ+HjimvxAez0QpnbAXcGUawA/oXfAfasZx+1tuSv2jcAJG2CC/Fy/PuOJ83ssq6ouwx3MZkC4Vfs3fE7md74t7W98ScqdsefYNQlZIGUlzOA3fCnKL6Ad4mL8K6jy8hdTCZJGb4HyZSYLJBMkiyQTJIskEySLJBMkiyQTJIskEySLJBMkiyQTJLSCiS768pBab9qz+66clBagWTKQWm7mEw5yALZwJA0TdKErqqv4QKR1BK/VSrE9gl+2GsL8RPk76bddj3VXYv/tmqaNspOeoE7w9craaykZcEsVon1krQ8vPggTvv5sD0Hp8psuEBw22GfQuxMfG7lloX46cBNZvbmeqq7Fv9tLWnWokYvcGf4eh/B3241JIrth0943lXSNlG8Gd/H6Zn9jZy0HAbIvwDmROv9QsOvAe6L4kPxGeO7dVI7qvpva0kT0rXLC0wVX28bk57X8vbiE5fPiz4/H58A/SRwbBS/BXi02raU8QryY/x9ss+w9hXkJ7i1cF5bBdVpqOoU2uMFbo+vN+Rrzdv7KG57qDAMF9W0Qrw5pE1ShjmpqwUi6RO4d+Uw3P9RiffH34fbXKWsejy7nUXNXuD2+npD3la9vfhBnxBmwgt/UfVJuLVifMg7APc5V335UhkE8g7QJKk77ledZWbTJQ1jzRXkdNwr83iqIGufZ7ezqMUL3B5fL6S9vY/igts/1PmWmc2X9Aawi/xt5sNwM9jT1TaiLF0MuLXyh/hZAf6qsT7BijkGuKJaQWXoYqjDC2zt8/VCwttrZq/iHprmsEwL8WW4YawSf8LMVlarqCxXEPAzYjlwT1h/F3+Qyjh8x95ZyRC8uhPwAW0v4HAzW0AJuhjr2Pt3a/H1gnt7T8MP/iRJJ9vaX4lXxiECJkfxafjzRprxV8BWpwR3MX3xS+87+GCrEt82ip8WxXsBs4BhYb0PsEk7667Ff1tLmqKvtxYvcN2+3uJdDG17e0fjd4IrgF2i+GH4ldmAITXtoxIIpCk0eDHRe2TxM8nwMcUWhZ1/w3qqu5nq/tta0lxA5OsNsWpe4Bbq9PUWBVIQSeztrZx0C1vZ1yvxq3P3WvbRBvdjnaSL8Q2/rmriTIcpwyC1XhYTnpUhqZukrRvcno2aDVEgLcDOkuYAz7P2C5cz65kNrovJdC0b4hUk04VkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJJkgWSSZIFkkmSBZJL8H7MP452NxHk9AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "eta_vol = 1- params['c']*(states.loc[2, 'v']/states.loc[3, 'v'] - 1)\n", "m_dot = params['V_dot_disp']*eta_vol/states.loc[2, 'v']\n", "W_dot_c = W_dot_c_m_dot*m_dot\n", "\n", "ax = plt.axes([0,0,0,0.2])\n", "ax.axis('off')\n", "plt.text(0, 1, \"$\\eta_{vol}$=\"+'{:.02f}'.format(eta_vol), size=14)\n", "plt.text(0, 0.5, \"$\\dot{m}$=\"+'{:.02f} kg/s'.format(m_dot), size=14)\n", "plt.text(0, 0, \"$\\dot{W}_{c}$=\"+'{:.02f} kW'.format(W_dot_c/1e3), size=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Condenser heat transfer rate\n", "\\begin{align}\n", "\\dot{Q}_{cond} = \\dot{m}(h_3-h_4) \\text{ or } \\eta_{cond}\\dot{m}_{a,cond}c_{p,a,cond}(T_{cond}-T_{bldg}) \\\\\n", "\\dot{m}_{a,cond} = \\frac{\\dot{V}_{cond,f}}{v_{a,cond}} \n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMMAAACGCAYAAABg6j51AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFlZJREFUeJztnXm0XFWVh79fwphBFIMjCQFpCUgmZBIZAmlARIU4RemF0oahmYwDLStp1BZQRAUNKyIEXSZq0yB0N4EAHTTJE6NhNCZ5QVRIgDAYaSMgSQgJ7v5jn8q776aGW1Wv6t56nG+tu17Vufued+6wz1Rn/67MjEgkAgPyLkAkUhSiM0QigegMkUggOkMkEojOEIkEojNEIoGOcAZJIyS9KGl83mWJ9F/UCb8zSNoOGAmsMbNNORcn0k/pCGeIRNpBR3STIpF2EJ0hgqQuSTPzLkfeFN4ZJI2W9BNJT0t6WdJjki6XtHPeZauEpHMlLZf0QtiWSDqxXpsG8x0o6RJJqyW9FP5eGsZdjZzLv0haL2mHRNoOkjZIWpGy/QdJJumYRv5X3hTaGSSdCjwIvABMAkYB04DTgFvyK1lNngQuBA4ADgQWArdIGlOnTSP5XgicC3wav15Tw/dpDZ7LQmAQcHAi7RDgeeDtknZLpE8ANgG/bvB/5YuZFXIDDge2AGeW2fchwIDD8y5nHeezDjirWZtaxwDzgDkpmznAvCp5dAEzE98nAs+V8gWeAr6Y2P8l4MfAr4CPJNKvBxblfa0b3YrcMswAusxsVpl9i8Lfse0oiKTp4XeOatsRFY4dKOljwBAq1JhZbOo4ZjFwtKRRwW4/4Bjgjoz5fgj4H7wSujYkLwKOTpgdjTtQVyp9Aj33pvPI2xsr1FSj8Zp/UoX9bw3766pFmyjPrsDeNbady5zDi3jr9hxwYoXzrGpT7zGAgK8Cfwc2h+t0aY08u4CZwJl49+e41P4pwEZgR2Cn8PltwHHA74LNKDqstd7mOuRdgAo359RwYUdU2P/+sP/QFpZhPrB/E8fvEJzkQOAy4P/S+WWxqTdf4GPAmvB3dLiW64ApVfLsCsdsBt5VZv9e4XpPwFuCNSF9MPAy8CbgbGA9sH3ez0/D9yzvAlS4OVPCxd+twv7/BR4m/GjYojI8DmwXPk8PtXG17Yga+f0c+EGzNrWOCQ/11JTNRcAjVfLoAm4HngC+X+66Ao8B/w58BfhxIn1JcLwbgfl5PzvNbA1Nt7WBZeHvUcDNyR2SpuB94OMt3A1Jw/FmfgRec74PnwG5DngNsBqYbGYbJN0GLAWOBYYDHzCz34S+9SxgKDAbWGdmW8K/vQb4aY0yP1Vj/wC8m9GsTa1jBgGvpGxeofbM4WrgfNwxZkk6s3R9A6Vxg/ABeYku/H5MAK6sr+gFI29vrFJb3Q48g9c6e+FTe98FXgI+mrDbAVgBHB2+7xK2lcD4kHYFcE74vJqeWZJzgEuA7YB7CN0N/MGf3UTZvw4cga+nGo13Z/4OnJDVBjgPeLiBfGfjU7AnBrtJwLPAFVXK20WYTcLHAmvwikQJm0/g06abgLcl0t+DT30bcHDez01Tz1zeBahyg3YON3s1Plg04BfAnim7yaS6FsGBZiS+nwF8Ga/1H0+knw98BvgovacWLwE+10TZZ+PdrE3An/GuzPH12OBdEmsg36HAd4LdRmAV8DVgpyrl3eoM4XvJIWbRs35teLgHa1LHDsHHGs8DA/N+bprZOmahnqRv4w/5AWb2TCL9UvwGXZtI+0pI+374fhXwS7zG/LyZfTikz8JbgQm4k1wX0u/Aa9IF7Ti3SDEo8u8Maf4Vr/EOSKWvBfYHkDRA0uuBp4F9QtoY4B+BW4PdssSxY/Au1l8SeUzEm/7lrTqRSDHpmJahEpKGAjfg/ePNeF/7N/jsxki8P3uWmXWHFmKBmc2VJOBRM9tL0jDgTnzs0A0cZWYj2n4ykVzpeGeIRPqKTuomRSItJTpDJBKIzhCJBKIzRCKBjnCGKBUTaQcdMZsUpWIi7aAjWgYz22Jmj9TrCJLeJ2mmpJGSTmpV+SL9g45whiYYg//ifDzhF+ZIpBL9whkknSLpXkkrJM2XNCjsGgv8AV/wd5qk30raJb+SFpMoFeMU3hkySsXMN7NDzGw08CgeqwDeGjyArzM6zszGmdnzbSjzkZJulfRUkE45rYxNI1IxQyV9R9LjkjZK+rWkg6rYTw//v+EHPUrFFIQ6pGKmSLpf0jJ8ZetLknbCJwjW44Pvx9pXcobga5ym4suoy9GIVMz38S7fJ/F4hruAn0t6a9pQ0qH40vVmFxxGqZi8NzJKxeAPxmxCQD7eLRqOP2A3AbsD9+R4Hi8Cp2W0rSgVg8d3bAFOSqU/SCrgHw9uehSPQOsiEatQIe9eNkSpmMKRVSpmf2CJmW2UdBYeN70m7FuOO8bTzRSkGamYjPlnkYrZDhiIR/ol2YhXCklmATeb2cIGyvKqlYopZAy0pNF49+GDFUxK44UteA01V9I/4U16qR87Jnx/CNgj9G9PMbMV6cwy0Bcx0NsQznMJLr/yIi6NU7Z8ZvY3SUuAiyR1A38CPg68C3gkkecZuHrGqQ2U50zgm8CHzeyuxK5FwExJO+Ix0IcCp+PRcDPCsaOAN+PXvDPJu2mq0GznLhXTh+dSsZtEnVIxeDjmL8K5bwHuA34CPBT274PHO49KHNNFtm7Sq14qpqjdpNLMRaXB57nA74F7W1WAMEVbin5rSTfJzF42/zHxATObBvwW+GwV+0fN7Ci8OzXczA4GtsfjxMFbiWFAt6QtkrbgCiPnhO/VlDeW4wIMU0LgU/L/rsJjqieErSukr8fHLKX0xWa2OfMFKBiF7CZRp1RMixiFazNBi7pJZcgkFRMewvWSXofPLn0h7LoFn0pO8kPgj7gowMtVso1SMXk3TVWa7kxSMcF2ODAX10NaCewJvAMfjHYDtwGDgu1twMV4X/1JXGAAYD9cp3QZXjsvbaLsQ4BxYduAz76MI9HtozGpmOOBE8L5HYu3JPdSpWtCnbNJRKmY4m1kl4opom7ShFDe9DY7YTOb+qViPopPmW7CK4qZwC5ZH/SsNkSpmGJTRSpmMv7r8pRE2sfwgeDU8P0M4C14M95tZnuE9PPx6cqngSPN7LyQfgnwVzPr7GY/UhdFHUCXo5JUzGh8ViXJvvRMsZZsHsJ/k7g/ld5Nz4K+Eu9MfY+8CugYZzBfxn25md2e2hV1kyJ9Qsc4QxVmA3tJWonrJe2L/xA3KvzQ9j18wL0JbwmWA4TpwzeY2dpgf5ikpfi6pyfN7Nl2n0gkXzpmzBCJtJr+0DJEIn1CdIZIJBCdIRIJFN4ZokxMpF0UfgAdZWIi7aLwLYM1KBMT6Xv6QjigyOIDhXeGVztZxAUqHHeOpNWSXpL0YLkl5lls6izrbEnzyqRfE5bT9CnJfMv9b7lu1gZJX82SX3SG4pNFXKAXYb3WDHzZ9nh89e6dkkbUY9MXhB8334+vKm5bvkFM4r+AaWb2b1nyLLQzSHq9pO9J+lOovbolHZd3udqJmd1hZtPN7GZ8iXcWPoevkL3OzH5nZufjq1zPrtOmKpImSnouxJ5X4iA8rHVxljwkDZb0ozBpslbSNEnzJM3Omq+kqbiSyOlmNiPr+RTWGSTtjq8jejMeC70/8C187XzH0gZxgR3whYZ3pXbdBRyW1SbD/yknHFCOk4Hbreed2rXyuAIP6pqEBw2NxeM+MuUbVhxfBnzQzH6c5VxKFDXSDeBafB3RJOuZ8nqkin2n0OqouWH4svS1qfS1+ILFrDYVqSIcUI6TgC9myUPSEOBTwCfM7GchbQoehJUl32Px91+/r8yCzpoU0hlCv/W9wEFW9LnfDEiaj79yt9vM1uH6SK0mfd1UJi2LTZqTgLPw+I8l1Qwl7Y1HKc7PmMfb8JjurUvyzWx9UAPJkm83Htj1ZUm/MrPnapxLL4raTRqPR7c9mHdB+oit8dSt7ibhChuv4IoVSd5AT0uQxaYSFYUDynAy/nbV9RnzKH2u5ZCV8n0G72LtgisNvq5GPr0oqjNsxlutoeV2SrpF0o1ySclVkg4uZ5ewryRMXM52uKS5kpZKWilpT0nvkOuadku6rXR8+HyxXCf1SUkHhPT9JC2WtEzSZ4F1ib7tNfTER1fa0kH9mTGzl/FK5NjUrmMJAmVZbKqwGg9rPQ4XDqjmECfRWwa0Vh6P4Pd+6/0M1zqtoF4pX8zsqZD3YGCBPL4lG3nHnZbbgNfjXYnr8cD+Ubho1diwfxUwPXz+OHBTrfwSn68mJdGY2FfEeOqq4gKUFw6YjCthnI7Hd8zA9Zv2qMemTFm6qCIcgMeWzAufd8Mf7DfWmcc1uC7uRFyk4Xo8vvqHNfLd+r/D9zfi3aZlwLAs17qQLYOZ/QWfQ94Df7DuwW/e2jDI2hH4RjBfCdRqDrcRJq5gNwm4z8wWhXI8j6tR/NzMlgabh4Hd5C9jH2A9syAD8Yi5DwIPmFmpn/t7mouaOxBX/ViKiyR8JXy+OOwfRojqK2FmNwKfAS7CFTQOB95rZo/XY1MNM3sUr4HfA1wbavcBePcW/P7dbx48VU8eFwC/xCMTF+HX7gF67lnNfEPea+mRvlwk6Q1ZTqqjNnzq767E908CV1axLytMXMH2UlLCv/jDd3ri+1XAR3DBrpsT6bPwmZhLgTMS6XcAE/O+bm26N3cB3wuf5wJf6IM8d8SlND/fl/mW2wrZMtRgDDBS0vZhgHQuCfEqSQvUW6K9kjBxOdsYT90AkobJXxN2FPCzkPwr4D8byGt8GOPtLV+pPAcfO97YTL6ZyLs2aaCmuBr/UeVevE94YmKfcC2inRNpY/C+/d24FtHdVWyH4uJlK+npOgzChcdWhBtRGgtcRRh7hLxWhc/DcAWOpXhs9RN5X7M23JOFeN//EhLCYw3mNR7vFv0N+CveVXpnO86j8Eu400haDJxqZqvL7NsX76J8LkM+mW0jrw460RmewGc8OqvgkcLTcc4QibSKThxARyItITpDJBKIzhCJBKIzRBomrNdaFtZsfSnv8jRLHEBHGkbSa8zsBUkD8Yizs83st3mXq1FiyxBpGDMrRR3uQM97+MqiAqtilIjOkBOqU5lC/q7oSxLHrJZ0qVxXqmRzrqTlkl4I2xJJJ7b4PO4lvHko3SqoRaoYrSI6Qw6oMWWKC/F1WJ/Gl7RPDd+nJWyeDHYH4KtdFwK3hDVVLcHMDgHeCoxTeDsqtE4Vo6Xkva6lyhqVPwPnpdJG4Wvwx+VdvibP7V7gulTaH4HLqhwzD5iTSptDYg1/hePWkVqJm9q/Ox5ZNhl3ng34AsRRuEPdHdLuo8J7uUM+FwIXJL4fjC9a3C5876L3e+MmAs/REw8yGPgRHlOxFnfyeTQRC1LvVuSWYQUe3JHk68ANVrBBWj2hnGpcmWIxcLSkUSGf/XD1iDsqlGmg/N12Q6gevTYu/D0bX65+CL5sejYeMzKdnvdLb13HJem1koaFzzvhUWsP00OrVDFaRiEFAQIr6LlRSDoMXw69T8Uj8qMexYtGlSkux1fVPiTpFfzefdXMrk4aSRqNv9Z3J7yWnWRmK9KZJRiLR5JNthAwI+lneBDUKPNAKyT9ApftKbErcJOk7fHu9k/NLKlo1ypVjJZRdGc4JfH9cryZzRSJ1WrUvOJFvcoUk/F3MZ+CLzEfB8yQtNrMfpCw+33Y91rgQ8AcSROsJ/IuzTi8q5V0zhHAf5ccIZG29eWQZrYKb+G2QS1SxWg1Re8m7SZX1TsJD5j5Ws5lStKo4kWjyhTfBL5lZjeY2Qpzgawr6T2AxsxeNhdqfsDMpuFxGZ+tku9YPKw2yXi8dUnbLSUbrVLFaClFdoaVuJziaNwJLgs1MNC5ihfWuDLFINyJkrxC7Xs4AB8DlLsug/FaeWkibVf8BejJtOG4SENWZ2iVKkZryXtmpcZMyCP4A7IG2Cm1r5MVL7KoV/RSvcAHtE/iinEj8YHms8AVCZuv44POkXglchleoZxQoRzvwh1qSCJtIh58v30i7f145NmADOfWElWMdmxFHjOAP8QnA58ys7SixRRJH8Ef9uFUns/eRvEizLKkFS/eovoUL/7a6EmZ2Y3y+OqL8EFpN9sqU6RVL87HHfNqvEv1DP4wXZyweRPwk/D3ebxbcoKZpfvuJcYCfzSzFxNp44FuM9ucsltmZlmEjzOpYkiagDvItfLY9Avw6dVb8Yrh27jcSyUlk74njxq/2Y2oeFHYjRapYrRjK/KYoRpR8aK4tEoVo+V05KrV8ADPxfudC4FjzOzIMCB7DJ8f3xhshwI34P3ozXhf/Df4RR6JS9yfZWbdkq7CZ0HmhrweNbO9wo9Ld+Jjh27gKDPr05d6vNoJDnAdXkltwWfBLjCztuntdqQzVEJR8SLSBP3KGSKRZujUMUMk0udEZ4hEAtEZIpFAdIZIJNBRziBpRFj0Nj7vskT6Hx01mxTifUcCa8xsU87FifQzOsoZIu1DUhe+Rum8vMvSLjqqm9SfkTRN/qqtFyQ9G5aS17WEOcRVmFKSLJKOlHSrpKfC/tP6tPC1y9URKhnRGYrDBHxF6mF4DPAW/PWtu2Y5WNKhwBmUXzM1BF9GMhXY2BeFzUpY1tIRKhmFdQZJu4dabLKkhZI2hMCaUZIOlHR3SLtP1SVWOgIzO97MfmgeRroCOBWPDXh3rWMl7QL8BzCFMkvLzewOM5tuZjfj8Q11I2mipOfCwkgkDZb0ozChsTa0bPMkzU4dehAej704HHekpHvCcc+HAK32BvFUoLDOQIOqDUWiznDQNEPx+5MlbmIWvvR8YZ8VPoGaU7PYqpIRJkDm4o4xFr+nM9g2gi8Xihzc06hqQ5GoRzUjzQx85WY6FrkXks4A9sZbkj5HzatZJFUyXoMLFdxm/tpb6C0vkytFdoaGVBvyRs2rZiDpSvzlioebWcVaU9I+eHz4Eeax1X1NU2oWSqlkmNm60I2aL2kBsAB/of2aFpS9borcTWqFakM7aFQ1g3DMt4GP4zEaq2r8r1I3sVvSFklb8K7LOeF7WSGAOmhWzWIblQwz+2e8e3Q38AHgD5KOb7KcfUPeYYIVQv4G4/3IdyfSdsUv/thE2vCQ9vbw/RRcunEFXhsNqvI/SnHTS3FxgD1D+jtwEYJu/JW3g0L6bXi88RK8O3BASN8P7wMvwyVZlqbKvHeNLfnq3Rl4ZN6+Ga/Ta/EIvOR2Px5Mvz9lXkOLxxefliHvLmAmXrM/Qe/g/SG4oMHkhP0gfHwzO5H2SxLhtRX+z53A9Xk/c2ZWWGdoSLWB5hQzdsRnPXJRzQC+i0fdHYMH9Je25DXopZhR7SFOpQ2hR7JmA/Cl8LmadurWfGhAzYIyKhnAnriCx2HAHsDR+JjporyfObPixkA3qtowJfxwtQwfaFdSVthGMcN8ecfJbKuasVudqhmNxkafg88gLcC7JqXtgoRNWjEjKwfiLeBSYGd8dm4pvZU1KmI+2J2Ax35fG7pMF+A1/634i8uX49pQpWteTiVjA/B24CZcxGEOPiV8eQPn1Pfk7Y19tdGkYkZIj6oZjV//XmoW9JFKRju3orYMjdCsYgZE1YzMZFCzaEglI0/6zUK9ZhUzzGyxXNIwqmZkoAhqFn1Nv3GGSkTFjEhW+r0zRCJZ6U9jhkikKaIzRCKB6AyRSCA6QyQSiM4QiQSiM0QigegMkUggOkMkEojOEIkEojNEIoHoDJFIIDpDJBKIzhCJBKIzRCKB6AyRSCA6QyQSiM4QiQSiM0QigegMkUjg/wFaTflPbVzmlAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Q_dot_cond = m_dot*(states.loc[3, 'h'] -states.loc[4, 'h'])\n", "cP_a_cond = CP.PropsSI('Cp0mass', 'T', params['T_bldg'], 'P', 101325, 'Air')\n", "v_a_cond = 1/CP.PropsSI('D', 'T', params['T_bldg'], 'P', 101325, 'Air')\n", "m_dot_a_cond = params['V_dot_cond_f'] / v_a_cond\n", "\n", "Q_dot_cond_alt = params['ETA_cond']*m_dot_a_cond*cP_a_cond*(T_cond - params['T_bldg'])\n", "\n", "ax = plt.axes([0,0,0,0.3])\n", "ax.axis('off')\n", "plt.text(0, 1.0, \"$\\dot{Q}_{cond}$=\"+'{:.02f} kW'.format(Q_dot_cond/1e3), size=14)\n", "plt.text(0, 0.75, \"$\\dot{Q}_{cond}^{alt}$=\"+'{:.02f} kW'.format(Q_dot_cond_alt/1e3), size=14)\n", "plt.text(0, 0.5, \"$\\dot{c}_{p,a,cond}$=\"+'{:.02f} kJ/kg/K'.format(cP_a_cond/1e3), size=14)\n", "plt.text(0, 0.25, \"$v_{a,cond}$=\"+'{:.02f} $m^3$/kg'.format(v_a_cond), size=14)\n", "plt.text(0, 0.0, \"$\\dot{m}_{a,cond}$=\"+'{:.02f} kg/s'.format(m_dot_a_cond), size=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Evaporator heat transfer rate\n", "\\begin{align}\n", "\\dot{Q}_{evap} = \\dot{m}(h_2-h_1) \\text{ or } \\eta_{evap}\\dot{m}_{a,evap}c_{p,a,evap}(T_{amb}-T_{evap}) \\\\\n", "\\dot{m}_{a,evap} = \\frac{\\dot{V}_{evap,f}}{v_{a,evap}}\n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMMAAACGCAYAAABg6j51AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFXRJREFUeJztnXm0HHWVxz9fViFBkc0w7IEwoGRDCYsQIkzIIDIJMAwIIoxhGRbHY0bFsCgekUUdMEdAEI8EQYZFZgLCcIgQnqxmCMMSdiFBAyKoCDGEEHDu/HF/nVTqVXVXd7/uqn78PufUea9/9evbv6quW7+l7/2WzIxIJAKrld2ASKQqRGeIRALRGSKRQHSGSCQQnSESCURniEQCPeEMkraUtETS2LLbEhm8qBd+Z5C0BrA1sMjM3i65OZFBSk84QyTSDXpimBSJdIPoDBEk9Um6qOx2lE3lnUHSSElXS/qdpOWSXpB0vqR1ym5bHpKmS3pQ0mJJf5D0c0k7NVsnw+7Jkh4L71ks6QFJB6TqnCXJUtvv2ziWf5H0pqS1EmVrSVoqaX6q7ojwefu0+nllUmlnkHQU8BCwGDgI2AGYDhwDzCqvZQ2ZAFwC7AHsA7wL3CFpgybrpHkROBXYGfgYMAeYJWlUqt4zwKaJbWQbxzIHWBcYlyjbFXgD2F7SxonyCcDbwP1tfF55mFklN2BP/AI5PmPfIYABe5bdzoLHMhT4K3BgO3Vy3vcacELi9VnA403a6AMuSrzeF3i9Zhd4CTgzsf9rwFXAfcChifJrgLvKPt+tblXuGWYAfWb2w4x9d4W/o7vREEmnhd856m171TGxHt4L/7nNOsk2rS7pcNyJ0nfi4ZJekrRQ0rWShhexGeweAvwXfhO6LBTfBXwiUe0TuAP1pconsPK76TnWKLsBWUgaiQ8FDs6pUpsvvNudFnEpcH2DOi/V2TcDeAR4oM06tXPzAPA+YAlwkJklx+5z8WHk08AmwBnA/ZI+YmZ/amD7eOA7wD+a2ezErruAiyStDQjYDTgWWBTajaQd8CHZnHqfUWnK7ppyuu2j8GHQljn7Dwz7dyu7rQWO5QLgd8Dwduok6q4FbIfPGc4F/gjsVKf+UOBVYFqdOn34hf0OsHvG/uHhfE/Ae4JFoXwIsBwYBpwIvAmsWfY5b3Wr6jCptnLxVs7+k/FJ4txuNKbVYZKkC4FPA/uY2YIc2w3rJDGz5Wb2nJnNM7PpeG/yxTr1lwBPACMamH4MeBmYKkkpGwuA3+DOMAF3HszsTXyBo1Z+r5m90+gYqkpVneHR8Hfv9A5JU/HVlxMt3J7Ckt4tkuZJukfSMEm7Sbot8b6Jkq6XdISkuZLmS7pd0rph/yxJ14XlzgWSkqsnlwJjGmzzUu2cARyBX+RPZx1kkToFWA1YO2+npPfhq3AvN7CzEL+g9wN+mHYIVs4bavOFGn349zGBXh4iQTWHSeEavxX/Ag/Hu+lxwMXAMuCfEvXWBu4kDKnwsew3gfcDzyfq3Q18GNgwUXYJMDn8vwA4Lfz/aeCGNtp+Mb4cvA8+hKhtQ5uscwrwdOL1ecBeeJzWSHyY9H/A/ok638VvItvgS6C3hM/Zqk57+wirScC2+JDpckK4Tij/LL5s+jawbaL874N9A8aVfd20dc2V3YA6X9A64cteiE+UDfglsE2q3mH4ePuRsD0JfC3s+w0+0dwPuCaUfQV4EO99XgMm4ePql4A1Qp1RwB1ttN1ytrOarHMWYInXM8MxvY3PA+4AJqU++9pwPpaHY7oR+HCD9q5whpRD/JCV8WtbhPYtSr13KD7XeANYvezrZlA6Q8YXdiHeU2yaKj8bODLnPbPDHfQe4G+Bo8MFtU7Y/2z4kvcAZifedzRwQdnHHLfublWdM2TxZeB7+JJrkt8DE2tj3LD0WONJ4AvAQjN7BtgJeMDM3pJ0ArCxmS3Ce4KtJa0p6YP4BP2Czh5OpGr0fAi3pCHA1cCO+HxijplNC/uOw+cFO5rZcyFs4SZ8CDAHn7iOl3QJ3s3vgy8Xnmpmt3b/aCJl0vPOMBBIuhc4yswWlt2WSHn00jCpk2wJvFB2IyLlEnuGSCQQe4ZIJNATzhDVMSLdoCeGSVEdI9INeqJnMLN3zYPTmnIESZ+SdJGkrSVN7lT7IoODnnCGNhiFh11Mwn9wi0RyGRTOkBeJimfCPYvHOB0j6RFJHyivpdVEUR0D6AFnUDF1jNvNbFczGwk8D0wM5TvhodWPAfuZ2Rgze6MLbR4v6eaQemmSjsmpt6mkK+XqGMskPSmpX9h66j0nhXTOZZIeysmjaNpunc+L6hhVQMXVMaaGPIRH8ZDvZSGOX+YJKFvT3R/VhgKP43FRmQlKktbHE+oFHICHk3wej0bNRNJheJrlOcBYPPf5NklbtmO3AVEdo+yNguoY5Eeifgy4Adgc+FWJx7EEOCaj/BzgviZtzQUuT5X9Gji3Tbt9RHWMSvcMRdUx8iJRR+PDoy3w+P6WaTXtswFTgLkhu+7VMJ85JSPDrNaGtYCP4mHpSWbjIegt2c34nPesOkYlnUEr1TEuzqmSVMe4CviqpLtxdYbaOHYU7gxPAluFyXWrYlpNp30WYDhwEp5hNwl3/vPw8PEsNgJWB15Jlb+CZ8i1ancFcnWMH+PqGEk1kLuA3SWtHYafu+GO8EuCMyiqY3RsaDGY1DHyhknLgftTZecAT+XY+ZtwzHulyr/OqqmhTdkN+/uI6hjV7BkYJOoYDXgZ77WSPIVH0GbxR1xxb1iqfBNW7S2atVvjPa+OUUkRMVZVx/hZcodWqmNMsnB76gLtiohlcR+eippke/yi64eZLZf0EL5sfENi10Q8z7kluwkW4qtOfbg6xvGp81ubNwi4MlHex0p1jN7ODiy7a6rTdRdSxwh1R+AqEPPwfOdh+Lj2tkSdifgFfQTeo8wHbgfWDftnAdfhYgELaEPpAV9arc0lluKrL2NIDPuAXfBhyem4KNih+HLlyYk6aXWMw/BhybH4kukMfBi2VTN2M9rbR1THqLQzFFXHqKJUzASylS9mpuodgPeCy/Al4X9NXYBnkVDHCGUn4b+ZvI0PUcZnfH5duxn1VzhDeB3VMaq8ka+OUTmpmLj15lbVOUMWX8aVMHbGh1A1RgJfNrOfZrznGXwIdSZwrKSj8d5hvPnvEs/izjMKeMLMakLGY/EJZeQ9RFVXk/phHsZ9vvVXrYhSMZEBoZd6hjyuwMMHnpK0DP/RZ1rY9wRBKia8vgq4SdKRoV7yB7obgXtZKRXzYneaH6kKPZHp1mmiVEwEemiY1GGiVEwk9gyRSI3YM0Qigco7Q5SJiXSLyg+TokxMpFtUvmewFmViIgPPQAgHVFl8oPLO8F6nqLhAxvuKCAc0rNNCe2dKuiWj/FL5wxwHjKTNrM+V62YtlfStIvaiM1SfhuICaQoKBzSsM1CE6IAD8WdjdMVmEJO4EZhuZqcXsVlpZ5C0oaQfSPp9uHs9Lmm/stvVTczsv83sNDP7Gf4wwyJMwyNkLzezp8zs83iQ44lN1qmLpH0lvR5CW+qxCx4weW8jG5KGSPpJWDR5RdJ0+ZNcZzZh8wvAj4BjzWxG0eOprDNI2hyPLN0UOBiPK/ouHjvfs3Qoay5pv6FwQBPiAvU+J0s4II8pwK2JQMh6Nv4dT+o6CE8aGo0/4bSozW/iof8Hm9lVRY6lRpVjky7DI0cPspVLXs+V2J6BohNZc0nqCQf8XRN1cgnCAd/BhQPSDpXFZDxyuK4NSUOBzwGfNbNfhLKpQFacWD+beALXAcCnMgI6G1JJZwjj1k8Cu1jV136bxMxew/MoOv5RqdfKKCtSJ81k4AQ8DP6BRo2QtB2eqXh7ARvbAmsC/7OigWZvSnq8gE3wudUHgK9Lus/MXm/UviRVHSaNxbPbHiq7IQNNp4dJFBMOKCoukEWucEAOU4A7zcUDGtmo/d/IIbNsEmzujTvEHSEcvzBVdYZ38F5rvaydkmbJRbIelLRA0riseon6ecLEWXVHhAnbPEn3SBomaTdJtyXqTJR0fZ7dBu3rhAbTCsxsOX4TmZjaNZEg+1ikTh0W4mmt++HCAY0cYjKrSoHWs/Ec/t2vOF/hnKYV1LNsAmBmLwXbQ4A7JW3YoH2rvLlyG7AhPpS4BvgIrrF6LDA67G8qX5mcvOeMelXMpy4iLtCKcEDDOhlt6aOxcMBM4Jbw/8b4xf2hojbwm8ULeI7Kh8M18AZwRZ7N9OeG1x/Ch02PAhsVOtdlX/h1TvzHcdmTv+C6n7/Au/Wm85XJyHvOqVe5fGoKiAvQunBAwzp5zpC6mJPCAT8BZoX/P0d/QbO6NsL5uwoXJHsF+Cp+g/pBns0sZ0g4zqN4EtcmPesMdb6QPYDZiddHAxfUqX80GcLEOXXPBo7M2Tcbz7e+B9clyrTbbPsG2xbOU+3CvQn4Spv21sZTe/9toGzmbVWdM9Sjbr6ypDslbZaon5f3nFU35lO3iKSN5I8K2xvvxcF79v9o0s7YMBfbTh6pfCU+d7yuVZuFKftO0sKd4hL8R5W5+JjwgMQ+4cOZdRJlo/AJ2934cOLuOnWH4D8CPY0Pky5I7DsOH6tu18BubvsG84bnlC/C51i5Gk0F7IzFFxD+AvwZV/L7aDeOofIh3GlUJ19Z0o7AcWY2rf87W687UO2LVJtedIbf4iselWx41dsXyafnnCES6RS9OIGORDpCdIZIJBCdIRIJRGeIRALRGSItI+l+SY+GDMSvld2edomrSZGWkfR+M1ssaXU8/fJEM3uk7Ha1SuwZIi1jZrUU3LVY+VDKTFRhiZga0RlKQk3KtEh6QS4Vk95uTdXbVNKVkv4QbD8pae8OHsdc4FU8Mrdfr6AOSMR0iugMJaDWZFp2wcURatvOeCj3inxqSevjgWzCc4F3xJ/g+erAH4VjZrsCmwFjJK2ShBMCHgdUIqajlB3gVSdg61XglFTZDnhCypiy29fmsc0FLk+V/Ro4twkbp+N5Husmys4B7muhPZvjjnUYHnC3FM8D2AH4GB6MuBTPTc57UP2pwJdSZeOAP+FZi32smsewb2j/CeH1EDwXYgmexzAdf4LrzG59L1XuGebjmU5JzgOutYpN0prJa9bAyLQImApcbWZLE7umAHNDyumrkh6RdEqB1Mwx4e+JwDeAXfE8gpnAt4HTgN1xVY1poQ3rS9oo/F9Leno6ZbcTEjEdo5LqGIH5rPySkLQH/szh9AO/q0Az8i9tybQEJgLb4EJZSYbj2WsX4jeOMcD3w756k9fReGrlYWb2CoCkX+DP4N7BzP4Uyn6JD9EANgBukLQmPty+3szSspKdkIjpGFV3hiMSr8/Hu9lGT7rvOtaa/EsrMi01jgMezOghVwPmmdn08PphSSPwBKN6zjAGT5lMOuiWwH/WHCFR9iCAmS3Ae7hM1AGJmE5T9WHSxnKJycl4Ztk5Jbcpk2aGSbQn04KkTfCL6vKM3S/jGXlJnsIv4nqMBn6VKhsLpHWRRgMPN2pjoBMSMR2lys7wBK4tOhJ3gnPDHRjonPxLPdvKl4ApLP9i7cm0AByDJ/Bfm7HvPvoPI7fHM/oykTQEvzM/nCjbAM/nTpZtgauWFHWGTkjEdJayV1YarHI8h18gi4D3pfZ1RP6lnm0GSAKGYlIuq8i/hDLhwgOX59jdBb+oTge2Aw7F5wIn12nL7nhPNTRRti+wDFgzUXYgnoq5WoHjG3CJmG5sVZ4zgA+VpgCfM7NlqX1TJR2K//K5Bflr2VPwC+7mcCNaC1+RWixpjbASMh540cxqQ4x+tsMkb218dQW855raykGZ2XVycasz8Anp48AnbdX50Eb0v8tPAEYAn8mx+6CkKXhPeibw2/D3kjrNGQ382syWJMrGAo+b2Tupeo+aWREl8APxOU3msM/Mnpc0AXeQy+SCCl/Cl1dvxm8MF+LaR+nvvXN0y+sGcqND8i/1bPMel4Bp8vsZcImYbmxVnjPUo1PyL/VsvyclYFqkExIxHafqw6Q8rsKHLkfiv5jOhxU/Rm3HqsucV+Dj0KckLQv1a4oYT+BDiB0b2cad4UY8OnMIcKqZdXUdvFcws283rpXJNHxo+C4u1TO+m+d4UIVwq0PyL8F2lIAZ5AwqZ+gkihIwg57oDJFIoFcn0JHIgBOdIRIJ9JQzSNoyxPmMLbstkcFHT80ZJK0BbA0sMrO3S25OZJDRU84Q6R6S+vCQjFPKbku36Klh0mBG0smSHpO0OGwPSDqgwPumhyjaxXIRgJ9n5CKvJ+l7kn4j6S253tEunTuafm3sCVGA6AzV4UU8j3hnPO94DjBL0qgG75uA/4q+B54u+S7+2NcNEnV+hD9v7mg8Fmt2qLMZHSZEBfSGKEDZQV11ArXaTlLv9Q0PKzmhyfcMxUOyDwyv18EdZHKq3kPA2XXs9DEACfwkRAHC6/F4ItESPER7LrBT2efarNqBek0nqVeNJjPgku9bXdLh+IVdJOEnyXp4j//n8HoNPOc6HQr9FrBnweNoJ4F/hShAWAC5CY/vGo1/pzNw5y2fsr2xzp2pJoWSTBD5PvAHVk2+uQK4ruz25hzDBnjgYL0t+Uy5kfgd891w7E0/Dw4XJngYWD1Rdj8epr4Z7hifwS/AZ+rY6cPzpo/H7+D7JfYNxZOTDk+UDcEdcGbKzhPAwYnzYcDeZX83WVuVo1abTlKvGta8UMAz+HGvDxwCXClpgpkVSoyXdAF+t9/TzJJ326OAH+Pzkr8C/4uHWO/cwGRbCfxpUQAze03STOB2SXfiGYg3WAi/L5sqD5M6kaTeVZodJpnZcjN7zsxqChePAF8s+FkX4qmo+5grV6zAzJ43s73xO/oWZjYOv5gbReC2m8DfTxTAzP4ZHx7dDfwD8KykSQ3sdIeyu6acLnoIfgf7eGrIYcDoRNkWoWz78PoIfEI2H78brVvnM0bgE755+BBiGLAbcFuizkRcDyjXNp70fh3eOy0AxqXaXHiYlNHGObhQWKPzNQOfxO5Y8Px+EB+GHV+nTh8+TBqOp48m85Vrw6TDEvXXJTVMCuf12AZtuQ24puxrzswq6wwtJanT2yIB5+ET0K3xucO5uDrI/ok6WSIBFwOL8UnssMSWPHeTgP1x4bGJeI8zN3ku85wh/N90Aj/ZogDbhOPcA9gK+AQurnZG2decWXXnDK0mqfesSAB+AV8d/r6BD1H2N7OkCFeWSMBJ4e+dqfJv4A9rB/gA7lyb43OYG4HTU+cyF2stgT9LFGApLl1zQziWV4Cf4gJx5VO2Nw7URhQJKPv8r5LAzwCIAnR7q/IEulmiSEAXKZDA37QoQNlUdZjUClEkoPvkJvBb66IApTHoo1ajSECkKIPeGTpJFAkYXERniEQCg2kCHYm0RXSGSCQQnSESCURniEQC0RkikUB0hkgkEJ0hEglEZ4hEAtEZIpFAdIZIJBCdIRIJRGeIRALRGSKRQHSGSCQQnSESCURniEQC0RkikUB0hkgkEJ0hEgn8P1PHQbzRGsorAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Q_dot_evap = m_dot*(states.loc[2, 'h']-states.loc[1, 'h'])\n", "cP_a_evap = CP.PropsSI('Cp0mass', 'T', params['T_amb'], 'P', 101325, 'Air')\n", "v_a_evap = 1/CP.PropsSI('D', 'T', params['T_amb'], 'P', 101325, 'Air')\n", "m_dot_a_evap = params['V_dot_evap_f'] / v_a_evap\n", "\n", "Q_dot_evap_alt = params['ETA_evap']*m_dot_a_evap*cP_a_evap*(params['T_amb'] - T_evap)\n", "\n", "ax = plt.axes([0,0,0,0.3])\n", "ax.axis('off')\n", "plt.text(0, 1.0, \"$\\dot{Q}_{evap}$=\"+'{:.02f} kW'.format(Q_dot_evap/1e3), size=14)\n", "plt.text(0, 0.75, \"$\\dot{Q}_{evap}^{alt}$=\"+'{:.02f} kW'.format(Q_dot_evap_alt/1e3), size=14)\n", "plt.text(0, 0.5, \"$\\dot{c}_{p,a,evap}$=\"+'{:.02f} kJ/kg/K'.format(cP_a_evap/1e3), size=14)\n", "plt.text(0, 0.25, \"$v_{a,evap}$=\"+'{:.02f} $m^3$/kg'.format(v_a_evap), size=14)\n", "plt.text(0, 0.0, \"$\\dot{m}_{a,evap}$=\"+'{:.02f} kg/s'.format(m_dot_a_evap), size=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### COP of the heat pump system\n", "\\begin{align}\n", "COP = \\frac{\\dot{Q}_{cond}}{\\dot{W}_{c} + \\dot{W}_{cond,f} + \\dot{W}_{evap,f}} \n", "\\end{align}" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "2.895050928533937" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "COP = Q_dot_cond / (W_dot_c + params['W_dot_cond_f'] + params['W_dot_evap_f'])\n", "display(COP)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [default]", "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.6" } }, "nbformat": 4, "nbformat_minor": 2 }