{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"Solar Energy Storage.ipynb","version":"0.3.2","provenance":[]},"kernelspec":{"display_name":"Python [default]","language":"python","name":"python2"}},"cells":[{"metadata":{"id":"TG6JHTaUoLa3","colab_type":"text"},"cell_type":"markdown","source":["# Analysis of Thermochemical System for Solar Energy Storage"]},{"metadata":{"id":"dxHfRagnoLa5","colab_type":"text"},"cell_type":"markdown","source":["## Antoine Equation and Ammonia"]},{"metadata":{"id":"hi1jHIqzoLa6","colab_type":"code","colab":{}},"cell_type":"code","source":["import numpy as np\n","import matplotlib.pylab as plt\n","%matplotlib inline"],"execution_count":0,"outputs":[]},{"metadata":{"id":"oZD_qG-yoLa9","colab_type":"text"},"cell_type":"markdown","source":["#### Part a. \n","The ammonia will be distributed as liquid at ambient temperatures which may range up to 50 $^\\circ$C. What would be the minimum operating pressure of the piping network in psig?\n","\n","#### Solution\n","The solution to this problem is found using Antoine's equation to predict ammonia vapor pressure as a function of temperature."]},{"metadata":{"id":"4YBCh44BoLa-","colab_type":"code","colab":{},"outputId":"544a5b2f-524b-4bcd-8f91-cb1074c91fe8"},"cell_type":"code","source":["A = 7.36050\n","B = 926.132\n","C = 240.17\n","\n","def Psat(T):\n","    return 10.0**(A - B/(T + C))\n","\n","T = np.linspace(-83,60)\n","plt.plot(T,Psat(T))\n","plt.xlabel('Temperature [deg C]')\n","plt.ylabel('Pressure [mmHg]')\n","plt.grid()"],"execution_count":0,"outputs":[{"output_type":"display_data","data":{"image/png":"iVBORw0KGgoAAAANSUhEUgAAAZYAAAEPCAYAAABhkeIdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcFNW5//HPI4iCgqgJ6gUNLhBxRQyo0ShGxSWiGBPE\nuIGaeF2i0RijUS/GmIsa43JjMPGnESQaolFZIoLruEQEFMcNBYyiYAQTlDUIDDy/P06N05BZepqq\nrqqe7/v16td0n67u/k7NMA91TtU55u6IiIjEZaO0A4iISGVRYRERkVipsIiISKxUWEREJFYqLCIi\nEisVFhERiVWihcXMupjZ02b2lpm9YWYXRu1bmtnjZjbTzCaZ2RYFr7nCzGab2dtm1q+gvZeZvW5m\ns8zs1oL2NmY2OnrNZDPbIcnvSUREGpf0EUsNcIm77w4cAJxvZrsClwNPuvtXgaeBKwDMbDdgINAD\nOBoYbmYWvdcdwFnu3h3obmZHRu1nAZ+6ezfgVuDGhL8nERFpRKKFxd3nu3t1dH8Z8DbQBTgeGBlt\nNhIYEN0/Dhjt7jXuPgeYDfQxs22B9u4+Ldru3oLXFL7XX4DDkvuORESkKWUbYzGzrkBP4CVgG3df\nAKH4AJ2izToDcwte9lHU1hmYV9A+L2pb5zXuvgZYZGZbJfJNiIhIk8pSWMxsc8LRxEXRkcv688jE\nOa+MNb2JiIgkpXXSH2BmrQlFZZS7j42aF5jZNu6+IOrm+iRq/wjYvuDlXaK2htoLX/MPM2sFdHD3\nT+vJoUnRRERK4O7N+g97OY5Y/gDMcPfbCtrGAYOj+2cAYwvaB0Vneu0I7AJMjbrLFptZn2gw//T1\nXnNGdP+7hJMB6uXumb8NHTo09QzKqYzKqZy1t1IkesRiZgcCpwBvmNmrhC6vnwE3AA+Y2ZnAB4Qz\nwXD3GWb2ADADWA2c53Xf2fnACGBTYIK7T4za7wZGmdlsYCEwKMnvKWlz5sxJO0JRlDM+ecgIyhm3\nvOQsRaKFxd3/BrRq4OnDG3jNMGBYPe2vAHvW076SqDCJiEj6dOV9xgwePDjtCEVRzvjkISMoZ9zy\nkrMUVmofWt6YmbeU71VEJC5mhmdw8F6aoaqqKu0IRVHO+OQhIyhn3PKSsxQqLCIiEit1hYmISIPU\nFSYiIqlTYcmYvPS7Kmd88pARlDNueclZChUWERGp18qVpb1OYywiIlKvgQPhwQc1xiIiIjGYOxee\nfLK016qwZExe+l2VMz55yAjKGbes5xw+HE4/vbTXqrCIiMg6VqyAu+6CCy4o7fUaYxERkXXcdReM\nHQvjx+s6FhER2UDu8H//BxdeWPp7qLBkTNb7XWspZ3zykBGUM25Zzfnss1BTA4fXu7BJcVRYRETk\nC7fdFo5WrFmdX+vSGIuIiADw/vvwta/Bhx/CZpuFNo2xiIhIyYYPhyFD6opKqVRYMiar/a7rU874\n5CEjKGfcspZz+XK45x44//wNfy8VFhERYdQo+MY3YMcdN/y9NMYiItLCucPuu8NvfwuHHrrucxpj\nERGRZnvqKWjVCvr2jef9VFgyJmv9rg1RzvjkISMoZ9yylPPWWzf8FONCreN5GxERyaO334Zp0+DB\nB+N7T42xiIi0YN//Pmy/PfzP/9T/fCljLCosIiIt1Pz50KMHzJ4NX/pS/dto8L4CZKnftTHKGZ88\nZATljFsWcv72t3DyyQ0XlVJpjEVEpAVavhx+9zt48cX431tdYSIiLdDtt8PTT8PDDze+ncZYGqHC\nIiISrFkD3brBfffBAQc0vq3GWCpAFvpdi6Gc8clDRlDOuKWZ85FHYLvtmi4qpVJhERFpQdzhV7+C\nSy9N7jPUFSYi0oK88EKYGv+dd8I0Lk1RV5iIiDTqppvgkkuKKyqlUmHJGPUPxysPOfOQEZQzbmnk\nnDkTJk+GM85I9nNUWEREWoibb4Zzz4V27ZL9HI2xiIi0AB9/HNZceecd6NSp+NdpjEVEROp1881w\n2mnNKyqlUmHJGPUPxysPOfOQEZQzbuXM+emncPfdyZ5iXEiFRUSkwv3mN3DCCWF6/HLQGIuISAVb\nuhR22ilMNtmtW/NfrzEWERFZx+9/D4cdVlpRKZUKS8aofzheeciZh4ygnHErR87PPw+D9ldckfhH\nrUOFRUSkQo0YAb16wd57l/dzNcYiIlKBVq+G7t3h/vs3bBZjjbGIiAgAo0dD167JTY3fGBWWjFH/\ncLzykDMPGUE545ZkzrVrYdgw+NnPEvuIRiVaWMzsbjNbYGavF7QNNbN5ZjY9uh1V8NwVZjbbzN42\ns34F7b3M7HUzm2Vmtxa0tzGz0dFrJpvZDkl+PyIieTBmDGy2GRx+eDqfn+gYi5kdBCwD7nX3vaK2\nocBSd795vW17APcDvYEuwJNAN3d3M5sCXODu08xsAnCbu08ys3OBPd39PDM7CTjB3Qc1kEVjLCJS\n8dyhd2+46ioYMGDD3y9zYyzu/gLwWT1P1RfyeGC0u9e4+xxgNtDHzLYF2rv7tGi7e4EBBa8ZGd3/\nC3BYXNlFRPLo0Udh5Uo47rj0MqQ1xnKBmVWb2V1mtkXU1hmYW7DNR1FbZ2BeQfu8qG2d17j7GmCR\nmW2VaPKEqX84XnnImYeMoJxxSyKnO1xzTbhtlOIIeusUPnM4cG3UxXUd8Gvg7Jjeu9HDtcGDB9O1\na1cAOnbsSM+ePenbty9Q90NO+3GtrORp6HF1dXWm8uR9f+bhcXV1daby5P1xEvtzyZK+rF4NW25Z\nRVVVae9XVVXFiBEjAL74e9lciV/HYmZfAcbXjrE09JyZXQ64u98QPTcRGAp8ADzj7j2i9kHAIe5+\nbu027j7FzFoBH7t7vZNCa4xFRCqZO+y7L1x9dZhwMi6ZG2OJGAVHEtGYSa1vA29G98cBg6IzvXYE\ndgGmuvt8YLGZ9TEzA04Hxha8pnaRze8CTyf3bYiIZNfYsaG4xDFgv6GSPt34fuBFoLuZfWhmQ4Ab\no1OHq4FDgIsB3H0G8AAwA5gAnFdwiHE+cDcwC5jt7hOj9ruBL5nZbOBHwOVJfj/lsH4XTlYpZ3zy\nkBGUM25x5nSHn/88jK1Ys44tkpHoGIu7f6+e5nsa2X4YMKye9leAPetpXwkM3JCMIiJ5N2ZMGKxP\n80ywQporTEQkx9auhX32geuug/7943//rI6xiIhIQh55BNq0gWOPTTtJHRWWjGmJ/cNJykPOPGQE\n5YxbHDnXrq27biULYyu1VFhERHLqoYegXTs45pi0k6xLYywiIjm0di3stRf86ldw9NHJfY7GWERE\nWoj774cOHeCoo5rettxUWDKmJfUPl0MecuYhIyhn3DYk56pV4Qr766/P1thKLRUWEZGcufNO6NED\nDj447ST10xiLiEiOLFsG3brBxImw997Jf57GWEREKtwtt8A3v1meolKqBguLmY0r4jaijFlbhJbQ\nP1xOeciZh4ygnHErJee//gW33QbXXht/njg1NldYDxpfJ8WA38YbR0REGjJsGJx0Euy8c9pJGtfg\nGIuZDXT3Bxp9cRHbZIXGWEQkzz78MMwJ9uabsN125fvcUsZYNHgvIpIDZ54ZCsovf1nez01k8N7M\nxtcztjLKzC4ys01Ljyv1qeT+4TTkIWceMoJyxq05OWfMgL/+FX7yk+TyxKmYs8LeA5YB/y+6LQGW\nAt2jxyIikqCrroLLLoOOHdNOUpwmu8LMbJq7966vzczecvfdE00YE3WFiUgeTZ4MAwfCrFnQtm35\nPz+p61g2N7MdCj5kB2Dz6OGq5nyYiIgUzx0uvjiMq6RRVEpVTGH5MfCCmT1jZlXA88ClZrYZMDLJ\ncC1RJfYPpykPOfOQEZQzbsXk/POfYfVqOPXU5PPEqck17919gpl1A3aNmma6++fR/VsTSyYi0oKt\nWAGXXw4jR4b17POksetYvt3YC9394UQSJURjLCKSJ9dfD1OnwsMp/6WN9ToWM7un4GF/YHzBY3f3\nM5sfMT0qLCKSFwsWwO67h4H7bt3SzRLr4L27D6m9AXMLH+etqORJJfUPZ0EecuYhIyhn3BrLOXQo\nnH56+kWlVE2OsUT0X30RkTJ4883Q/TVzZtpJSlfUlC5mNt3de5UhT2LUFSYieXDUUXDMMXDhhWkn\nCUrpCmvwiMXMxlN3pLKTmY0rfN7dj2t+RBERacjEifD++3DuuWkn2TCNncR2E/Dr6HZ8wf3amySg\nEvqHsyQPOfOQEZQzbuvnrKmBH/8YfvUr2HjjdDLFpcEjFnd/tpxBRERast/9DrbZBvr3TzvJhitm\nrrBjgV8AXyEUIiOcbtwh+Xjx0RiLiGTVJ5+E04ufeQb22CPtNOtKZD0WM3sX+DbwRp7/MquwiEhW\nnXUWbLEF3Hxz2kn+U1KTUM4F3tRf5fLIa/9wVuUhZx4ygnLGrTbnlCnw2GNwzTWpxolVMdexXAZM\nMLNngZW1je6ewdoqIpIfa9bA+efDDTdAh1wNLjSumK6wxwkLfb0BrK1td/efJxstXuoKE5GsufNO\nGDUKnnsOrFmdTeWT1BjLm+6eseGk5lNhEZEsWbgQdtsNHn8c9t477TQNS2qMZYKZ9SsxkzRT3vqH\nsy4POfOQEZQzbkOGVDFwYLaLSqmKGWM5l7Cw10pgNTk93VhEJCumT4cXXghrrVSiouYKqwTqChOR\nLFi7Fg48EL7/fTgzB/PExzpX2HpvvBfQtXD7vC30JSKSBSNGhOIyeHDaSZLT5BiLmf0B+ANwImHB\nr/7AsQnnarHy0j+snPHJQ0ZQzjh88glccQX8/vfw3HNVacdJTDFHLPu7+26JJxERqXAXXwxnnAE9\ne0KG698GK+Z047uBX7v7jPJESobGWEQkTRMnwnnnwRtvwGabpZ2meEmNsdwLTDaz+YQr72vPCtur\nhIwiIi3O8uVhjZXf/S5fRaVUxVzHcjdwGnAUdeMrFTCxczZluX+4kHLGJw8ZQTk3xDXXwNe/Dkce\nWdeWxZxxKeaI5Z/uPq7pzUREZH2vvgr33hu6wFqKYsZYhgMdgfGsOwllrk431hiLiJTbmjWw335h\noskhQ9JOU5qkxljaEgpK4bQuDuSqsIiIlNtvfgPt21f2NSv1aXKMxd2H1HPLwfWi+ZSXflfljE8e\nMoJyNtcHH8B114VrVuqbuTgrOZPQYGExsx809eKmtjGzu81sgZm9XtC2pZk9bmYzzWySmW1R8NwV\nZjbbzN4unPjSzHqZ2etmNsvMbi1ob2Nmo6PXTDazHZrKLCKSNHc45xz40Y+ge/e005Rfg2MsZvYe\ncGljrwWudffdG9zA7CDCWi731p6ebGY3AAvd/UYz+ymwpbtfbma7AfcBvYEuwJNAN3d3M5sCXODu\n08xsAnCbu08ys3OBPd39PDM7CTjB3Qc1kEVjLCJSFnfdBXfcAS+9BBtvnHaaDRP3GMuzNH1a8RON\nPenuL5jZV9ZrPh44JLo/EqgCLgeOA0a7ew0wx8xmA33M7AOgvbtPi15zLzAAmBS919Co/S/A7U3k\nFRFJ1AcfhGlbnnkm/0WlVA12hTUwtrL+7UclfGYnd18QfcZ8oFPU3hmYW7DdR1FbZ2BeQfu8qG2d\n17j7GmCRmW1VQqbMyEu/q3LGJw8ZQTmL4Q5nnw2XXAJ7NLE8Yl72ZymKmt04YXH2TzV6uDZ48GC6\ndu0KQMeOHenZsyd9+/YF6n7IaT+ulZU8DT2urq7OVJ687888PK6urs5Uniw+njmzL4sXQ58+VVRV\n5XN/VlVVMWLECIAv/l42V+LrsURdYeMLxljeBvq6+wIz2xZ4xt17mNnlhKliboi2m0jo5vqgdpuo\nfRBwiLufW7uNu08xs1bAx+7e6T9TaIxFRJL1/vvQpw88+2xYcrhSJLU08YYy1j2SGAcMju6fAYwt\naB8Unem1I7ALMDXqLltsZn3MzIDT13vNGdH97wJPJ/ZdiIg0YO1aOOss+MlPKquolKqY9Vi2iU4b\nfix6vJuZnVXMm5vZ/cCLQHcz+9DMhgDXA0eY2UzgsOgx0ezJDwAzgAnAeQWHGOcT5iybBcx294lR\n+93Al6KB/h8RTgLItfW7cLJKOeOTh4ygnI254w5YsQJ+/OPiX5OX/VmKYsZYRgD3AFdGj2cBfyb8\nUW+Uu3+vgacOb2D7YcCwetpfAfasp30lMLCpHCIiSXnvvTDJ5AsvQKtWaafJhmLmCpvm7r3N7FV3\n3ydqq3b3nmVJGBONsYhI3Gpq4NBD4YQTwplglSipMZblZrY10dlbZrY/sLiEfCIiFWXYMNhkk3CF\nvdQpprBcQhgk39nM/ka4QPGHiaZqwfLS76qc8clDRlDO9U2eDLffDiNHwkYlnAaVl/1ZikbHWMxs\nI2BTwpXyXyWc3TXT3VeXIZuISCYtWQKnnhpWhOzcuentW5pixli+GFvJM42xiEhczjgjdIHdeWfa\nSZKX1HosT5nZicDD+sssIi3d6NFhcsnp09NOkl3F9AyeAzwIrDSzJWa21MyWJJyrxcpLv6tyxicP\nGUE5IUwweeGFcP/9sNlmG/ZeedmfpWjyiMXd25cjiIhIlq1ZA6edBpdeCvvum3aabCtmjOXg+trd\n/blEEiVEYywisiGuuy5Mhf/EE6WdBZZXpYyxFFNYxhc83BToA7zi7t9sfsT0qLCISKmqquDkk+Hl\nl1veWWCJXCDp7v0LbkcAewCflRpSGpeXflfljE8eMkLLzTl/PpxySrheJc6ikpf9WYpSDujmAT3i\nDiIikjU1NeFI5eyzoV+/tNPkRzFdYb+hbjGujYCewBx3PzXhbLFSV5iINNdVV8GUKTBxYsudYDKp\n61heLrhfA/zJ3f/WrGQiIjnz2GMwYkS4XqWlFpVSFTPGMrL2RlgnZWnysVquvPS7Kmd88pARWlbO\nDz+EIUPgT3+CTvWuSbvh8rI/S1HMQl9VZtbBzLYCpgP/z8xuST6aiEj5rVoFJ50UpsH/xjfSTpNP\nRc8VZmZnA9u7+1Aze712Dfu80BiLiBTj4ovh3Xdh7NiWdb1KQ5IaY2ltZtsRVmq8sqmNRUTyatQo\nGD8epk5VUdkQxey6a4FJwLvuPs3MdgJmJxur5cpLv6tyxicPGaHyc06bFrq/xoyBrbaKN1N98rI/\nS1HMXGEPEiahrH38HnBikqFERMpp/nw48cQwDf4ee6SdJv+KGWO5EbgOWAFMBPYCLnb3PyYfLz4a\nYxGR+qxcCd/8ZrgAcujQtNNkT1Jr3vdz9yXAscAcYBfgJ82PJyKSLe5wwQWwzTZw9dVpp6kcxRSW\n2u6ybwEPuvviBPO0eHnpd1XO+OQhI1RmzjvuCGvXl7pu/YbIy/4sRTFnhf3VzN4hdIWda2ZfBj5P\nNpaISLKefRZ+/nN48UVor1WnYtXkGAtAdHHkYndfY2abAe3dfX7i6WKkMRYRqfX3v8NBB4XTiw8/\nPO002ZbIGIuZtQPOA+6Imv4L+Frz44mIpG/hQjjmGLjmGhWVpBTTq3gPsAr4evT4I8JZYpKAvPS7\nKmd88pARKiPn55/DgAHhds455ctUn7zsz1IUU1h2dvcbgdUA7v5voFmHRSIiaVu7Nkwsud12MGxY\n2mkqWzHXsbwIHAb8zd17mdnOhKnz+5QjYFw0xiLSsl15ZViz/qmnoG3btNPkR1JzhQ0lXBi5vZnd\nBxwIDG5+PBGRdNx1F/z5z+HUYhWV5DXaFWZmBrwDfJtQTP4EfM3dqxJP1kLlpd9VOeOTh4yQ35yP\nPx5WgpwwAb785XQy1Scv+7MUjR6xuLub2QR33xN4tEyZRERiMX06nHoqPPwwdO+edpqWo5gxlpHA\n7e4+rTyRkqExFpGWZeZM6Ns3XF0/YEDaafKrlDGWYgrLO0A3wjxhywlnhLkW+hKRrJo7N1wAec01\n4UwwKV1Sk1AeCewEfBPoT5iMsn/z40kx8tLvqpzxyUNGyE/OMWOqOOIIuOiibBeVvOzPUjQ4xmJm\nmwL/TZjN+A3gbnevKVcwEZHmWrIELrsMBg4Mi3ZJOhrsCjOzPxMuinweOBr4wN0vKmO2WKkrTKSy\nrVgBRx8NPXrA8OFguow7FrGOsZjZG9HZYJhZa2Cqu/fa8JjpUGERqVyrV8N3vgPt2sEf/witWqWd\nqHLEPcayuvaOusDKJy/9rsoZnzxkhOzmXL0avvc9qKkJ66o8/3xV2pGKktX9GYfGrmPZ28yWRPcN\naBs9rj0rrEPi6UREGrF6NZxyCixfHq5VadMm7UQCRa7HUgnUFSZSWWpqwpHK0qXwyCOw6aZpJ6pM\nSZ1uLCKSKTU14UhFRSWbVFgyJi/9rsoZnzxkhOzkrKkJ07QsXlx/UclKzqbkJWcpipndWEQkE2pq\n4LTT4LPPYOxYHalklcZYRCQXVq4MYyrLlsGYMZr+vlw0xiIiFWn5cujfP1z0OG6cikrWqbBkTF76\nXZUzPnnICOnl/OwzOOII6NIFRo+GTTZpfHvtz/SlVljMbI6ZvWZmr5rZ1KhtSzN73MxmmtkkM9ui\nYPsrzGy2mb1tZv0K2nuZ2etmNsvMbk3jexGRZMyfH6a+33//sApka40K50JqYyxm9h6wr7t/VtB2\nA7DQ3W80s58CW7r75Wa2G3Af0BvoAjwJdIsWIpsCXODu08xsAnCbu0+q5/M0xiKSI3PmhCOV008P\nK0Bq7q905G2Mxer5/OOBkdH9kUDt8jzHAaPdvcbd5wCzgT5mti3QvmARsnsLXiMiOTVjBhx8MFx4\nIVx9tYpK3qRZWBx4wsymmdnZUds27r4AwN3nA52i9s7A3ILXfhS1dQbmFbTPi9pyKy/9rsoZnzxk\nhPLlfOYZOPRQ+N//hR/+sPmv1/5MX5o9lge6+8dm9mXgcTObSSg2hWLtuxo8eDBdu3YFoGPHjvTs\n2ZO+ffsCdT/ktB/Xykqehh5XV1dnKk/e92ceHldXVyf+eR9+2JdLL4XLL6+iSxeA7Hz/cT8ux/4s\n5XFVVRUjRowA+OLvZXNl4joWMxsKLAPOBvq6+4Kom+sZd+9hZpcTJr68Idp+IjAU+KB2m6h9EHCI\nu59bz2dojEUko9zhF7+Ae+6BRx+F3XZLO5HUys0Yi5m1M7PNo/ubAf0Iq1SOAwZHm50BjI3ujwMG\nmVkbM9uRsKrl1Ki7bLGZ9TEzA04veI2I5MCqVXDmmTB+PEyerKJSCdIaY9kGeMHMXgVeAsa7++PA\nDcARUbfYYcD1AO4+A3gAmAFMAM4rOPw4H7gbmAXMdveJZf1OYrZ+F05WKWd88pARksm5aBEccwws\nXAhVVbDtthv+ni15f2ZFKmMs7v4+0LOe9k+Bwxt4zTBgWD3trwB7xp1RRJL1zjswYAD06we33KJV\nHytJJsZYykFjLCLZMX48nHUWDBsWvkp2lTLGoutYRaRs1q6FX/4Sfv/7MOfX/vunnUiSoLnCMiYv\n/a7KGZ88ZIQNz7l0KXz3u/DYYzBtWnJFpaXszyxTYRGRxL37LhxwAGy1VbgAcrvt0k4kSdIYi4gk\n6oEH4IIL4Npr4ZxzND1L3miMRUQyY8UKuPhiePLJ0P21775pJ5JyUVdYxuSl31U545OHjNC8nDNn\nhjGURYtg+vTyFpVK3J95o8IiIrH64x/hoIPg/PPhT3+CDh3STiTlpjEWEYnFkiVw0UVhWpYHHoC9\n9ko7kcQhN3OFiUhlqaoKhaRNG3j5ZRWVlk6FJWPy0u+qnPHJQ0aoP+eKFXDJJXDKKTB8eLjwcfPN\ny5+tUJ73Z6VQYRGRkrzyShiU/+gjeP31MJmkCGiMRUSaadWqMMfX8OFw220waFDaiSRJuo5FRBL1\nt7/BD34AO+8cTiPunOuFwCUp6grLmLz0uypnfPKQcdEiOO64KgYODFfQjx2b3aKSh/0J+clZChUW\nEWmQOzz0EOy+e3j81ltw4omalkUapzEWEanX3/8epmR59124885w0aO0PLqORUQ22LJl8LOfwX77\nhWlZXn1VRUWaR4UlY/LS76qc8clKRvcwHcuuu8LcufDaa6HAbLJJeD4rOZuinOnTWWEiwssvw4UX\nwurV8OCDYe0UkVJpjEWkBXv3Xbj6anj22bBk8BlnwEbqx5ACGmMRkaJ8/DGcd14YQ9ljD5g9G4YM\nUVGReOjXKGPy0u+qnPEpZ8bFi+HKK0MxadcurJty5ZWw2WZNvzYP+xKUMwtUWERagEWL4LrroFs3\nmD8/nOl1002w9dZpJ5NKpDEWkQq2cCHceivccQd861vhLK+vfjXtVJInGmMREQAWLICf/hS6dw/3\np06FkSNVVKQ8VFgyJi/9rsoZnzgzzpoVBuV79IDly0OX1513wk47bfh752FfgnJmgQqLSM65hxUc\njzsuXCH/pS/BjBlw++2www5pp5OWSGMsIjm1alVYW/7mm8NKjhdfDKedBm3bpp1MKkkpYywqLCI5\n8/e/h+6tESPC2vKXXAJHHqlrUCQZGryvAHnpd1XO+BSTcfVqePhh6NcvXNS4Zg08/zw88QQcfXR5\nikoe9iUoZxZorjCRDJsxA0aNCmd07bILnHMOjBsHm26adjKRhqkrTCRj/vlPGD0a7r0XPvoITj01\nzOFVu9iWSDlpjKURKiySZcuWwaOPwn33wXPPQf/+YSD+sMOgVau000lLpjGWCpCXflfl3HBLl4Yj\nk4MPrqJzZ7jnHvj2t8NaKKNGhfGULBWVLO/LQsqZPo2xiJTRwoXw2GNhHfmnnoIDDwyD8WPGwFZb\npZ1OJB7qChNJkDu88w6MHx9ur70Ghx4KAwbA8cermEj2aYylESosUi5Ll4aFs554Av7613AhY//+\ncOyxoajoAkbJE42xVIC89LsqZ52aGpg8Ga69Fr7xDfiv/4JbboFttw3Xnnz4IQwfDsccU39R0b6M\nl3KmT2MsIs20ciVMmxbO3nr+eXjxRdhxRzj8cLjqqlBc2rVLO6VIetQVJtKEf/4TpkyBl14KheSV\nV2DXXeHgg0MROegg+PKX004pkgyNsTRChUWK8e9/hwH2KVPqbp9+Cr17w377hUJywAHQoUPaSUXK\nQ2MsFSAv/a6VkHPhwnDK7003wSmnwG67hSnnL7ggrGty1FEwYUIoLE88EZb2PfLI+ItKJezLLFHO\n9GmMRSo4gOXZAAALxklEQVTekiXw1lt1tzffDF+XLYOePWGffeCII+Cyy8ICWW3apJ1YJN/UFSYV\nYdUqmDMnHGnMmgUzZ9bdX7QoFIzdd4c99qj7uv32YM06wBdpeTTG0ggVlnyrqYF//CNMd/LBB/De\ne+H2/vvh6/z50KVLWNO9e/d1b126aK0SkVK12MJiZkcBtxLGjO529xvq2SYXhaWqqoq+ffumHaNJ\nceZctiwUjY8/Dl9r78+bF64BmTsXFiwIZ17tsEO47bxzOMV3p53CbfvtoXU9Hbt52J95yAjKGbe8\n5CylsOR+jMXMNgJuBw4D/gFMM7Ox7v5OuslKU11dnYtftoZyuocxjYUL173961/htN1PPglF4pNP\n6u6vXRsuKtxuu7qv220HvXqFgrH99tC5M2y8cXw5syQPGUE545aXnKXIfWEB+gCz3f0DADMbDRwP\n5LKwLFq0KJXPdYfPPw9HD4W3pUvDbcmSdW9PP72IF18M4xeLFsFnn9Xdb9sWtt76P2+dOsG++4av\n22wTvnbqBO3bJzfWkdb+bI48ZATljFtecpaiEgpLZ2BuweN5hGKTWStXhjGD+m6ffhomLVy9Otxq\naurur1pVd6t9vHJlw7cVK0Kx+PzzuvsrVoRrNeq7tWkDm2++7m2zzWCLLcIptrW3Tp3CEcWJJ8KW\nW0LHjuG25ZZhW51VJdKyVUJhyZ1evcIZTK1br3tr1Qo++2wOkyaFbp/aW+vW4Y91mzbhceH9TTb5\nz1v79uF6jLZtwxK2tV9r77drV/+tvjGKhrz11hxOOimxXRSbOXPmpB2hSXnICMoZt7zkLEXuB+/N\nbH/gGnc/Knp8OeDrD+CbWb6/URGRlLS4s8LMrBUwkzB4/zEwFTjZ3d9ONZiISAuV+64wd19jZhcA\nj1N3urGKiohISnJ/xCIiItlScdcjm9l3zOxNM1tjZr3We+4KM5ttZm+bWb+C9l5m9rqZzTKzW1PI\nvLeZTTazV81sqpl9ranMaTGzH0ZZ3jCz6wvaM5UTwMx+bGZrzWyrgrbM5DSzG6Mc1Wb2kJl1KHgu\nMzmjPEeZ2TvRv5Gfpp2nlpl1MbOnzeyt6Hfywqh9SzN73MxmmtkkM9siA1k3MrPpZjYuwxm3MLMH\no9+7t8xsv5JyuntF3YCvAt2Ap4FeBe09gFcJ3X9dgXepO2KbAvSO7k8Ajixz5klAv+j+0cAz0f3d\nGsqc0r7tS+hybB09/lJT+zbFrF2AicD7wFZZzAkcDmwU3b8eGJbRn/tGUYavABsD1cCuaf58C7Jt\nC/SM7m9OGG/dFbgBuCxq/ylwfQayXgz8ERgXPc5ixhHAkOh+a2CLUnJW3BGLu89099nA+mcxHA+M\ndvcad58DzAb6mNm2QHt3nxZtdy8woGyBg7WEHyBAR+Cj6P5x1JO5zNkKnUv4paoBcPd/Re317tt0\nIn7hFuAn67VlKqe7P+nua6OHLxGKIWTv5/7FRcjuvhqovQg5de4+392ro/vLgLcJ+/F4YGS02UjK\n/296HWbWBTgGuKugOWsZOwDfcPd7AKLfv8WUkLPiCksj1r+Q8qOorTPhospa86K2croYuMnMPgRu\nBK6I2hvKnJbuwMFm9pKZPWNm+0btmcppZscBc939jfWeylTO9ZxJOFqG7OWs7yLkrOy3L5hZV6An\noUhv4+4LIBQfoFN6yYC6/+gUDmpnLeOOwL/M7J6oy+5OM2tHCTlzeVaYmT0BbFPYRPiBXenu49NJ\n1bjGMhO6RC5y9zFm9h3gD8AR5U/ZaM6rCL8vW7r7/mbWG3gQ2Kn8KZvM+TNS2n/rK+Z31cyuBFa7\n+59SiFgRzGxz4C+Ef0fL6rluLbWzlMzsW8ACd682s76NbJr2mVStgV7A+e7+spndAlzOf+ZqMmcu\nC4u7l/JH4yNg+4LHXaK2htpj1VhmMxvl7hdF2/3FzGoPl8uSrVATOf8beDjablp0gsTWUaYdspDT\nzPYgjEu8ZmYWZZluZn2ylLOWmQ0mdJF8s6C57D/3JpR9vzWHmbUmFJVR7j42al5gZtu4+4Kou/uT\n9BJyIHCcmR0DtAXam9koYH6GMkI4Ep3r7i9Hjx8iFJZm78tK7worHGcZBwwyszZmtiOwCzA1OrRb\nbGZ9oj9EpwNj63mvJH1kZocAmNlhhD71BjOXOVuhMUR/AM2sO9DG3RcScp6UhZzu/qa7b+vuO7n7\njoR/LPu4+ydZyglfLPfwE+A4d19Z8FTWfu7TgF3M7Ctm1gYYFGXMij8AM9z9toK2ccDg6P4ZlP/f\n9Bfc/WfuvoO770TYd0+7+2nAeDKSESDq7pob/duGcNH5W5SyL9M+CyHuG2FgaS6wgnAl/mMFz11B\nOLvlbaKzsKL2fYE3CH/Qb0sh89eBlwlnAk0m/CFsNHNK+3ZjYFS0r14GDslizvUyv0d0VljWcka/\nbx8A06Pb8CzmjPIcRTjjajZwedp5CnIdCKwhnKn2arQfjwK2Ap6MMj8OdEw7a5T3EOrOCstcRmBv\nwn8kqgm9E1uUklMXSIqISKwqvStMRETKTIVFRERipcIiIiKxUmEREZFYqbCIiEisVFhERCRWKiyS\nO2a2lYUlBqab2cdmNq/gcSZnkzCzIWaWyFxQZrazma0ws3ovojSzX9ROJx/DZ21nZn+Ops+fZmbj\nzWwnM+se/Qw+jeNzJN8y+Y9QpDHu/imwD4CZ/Q+wzN1vTjdVWG/D62YrXt+ZhIv3ip62w8xaufua\nIjd/293LMQPyGOB37n4ShLWECJMUTrawjtA/y5BBMk5HLJJ36yyPYGanm9mU6Ojl9qitlZl9Zma/\ntrAI3GPRFD5VZvZuNLUKZnaWmT0ctc+MJocs5n1vMbNqoLeZXWNhsbbXzWx4tN1Awqy7o6PXb2xm\nc6NpyrGwmNIT0f1fmNlIM3sBuCf6jF9bmFG62szOLGqnmP1P9D08R1ifqLZ9FzObGB1tVJnZLgXt\nL5nZa2Z2nZl9Vs97HgEs9WhadQB3f83dJxeTSVoOFRapGGa2O3ACcIC79wI2NrNB0dNbAI+6+x7A\namAoYd6zgcAvCt6mN2E9lH2A75nZXkW8b5W793T3KcCt7t7H3fcCOprZke7+AGGKjIHu3svDmiaN\nzRj7VeBQdz8d+AFhZtz9CeuiXGBhbY/G9kPvKO+ewLGsu5bLncC57t6bMAv0b6P23wA3uvvehKmQ\n6puSYw/glcY+WwTUFSaV5XDga8DL0YSimxLm4gL4t7s/Hd1/A1jk7mvN7A3Cyoi1Jrn7EgAzewQ4\niDBHWkPvu9LrZtQFOMLMLo222Zowp9qk6LnCo6v1F6IrNDYqPgD9gF3N7OTocQfCEci8el8ZHAw8\n5O6rgFVmVjs9/xbA/sBD0fcBdf+53M/dj47u38+6xVakWVRYpJIY8Ad3H7pOo1krYFVB01pgZcH9\nwn8Hhf9Tt4LHDb3vioLHbQn/8+/p7vPN7BeEAlOfGur+qK+/zfL1Mpzn7s808D7NYcA/o6Ou9RUz\naeBbhCMgkUapK0wqyZPAQAtrxNSePVbbbdTYEULhc/3MrIOFlfOOB/4GPFXk+7YlzLS70MzaAycW\nPLeUcLRR633CrNqst936JgHnR0WM6OyrTRrZHuA54AQz2yQaxzkWwN0XAR+b2YDovczM9opeM9XM\nvh3dP/k/3jG8/nHCWiKDa9uirsIDmsgjLYwKi1QMd38T+DnwpJm9RvijXLt6Y2P/Iy98bhph/YlX\ngfvc/fVi3zc6W20kYar7RwlL5Na6B7jL6k6J/jlwh5lNoe7oqT6/J0xVX21mrwPDaaKnwd2nAY8A\nrxPW/JhS8PTJwH9HJxu8CXwrar8I+GnU3hVY3MDbHw98Kzrp4Q1Cl9n8xvJIy6Np80UiZnYWsLu7\nX5J2luYws52Bv7j7PhvwHu3c/d/R/VOAAe7+3Wa+R2tCV9uWpeaQyqAjFpH8qwG2tgYukCxSbwsX\nOL4GnE1Y2bJoFlYdnEY4o0xaOB2xiIhIrHTEIiIisVJhERGRWKmwiIhIrFRYREQkViosIiISKxUW\nERGJ1f8HJ4xn2SNWmXgAAAAASUVORK5CYII=\n","text/plain":["<matplotlib.figure.Figure at 0x10b78b550>"]},"metadata":{"tags":[]}}]},{"metadata":{"id":"J8HM3uQ2oLbE","colab_type":"text"},"cell_type":"markdown","source":["At the operating conditions"]},{"metadata":{"id":"dnZS_eQboLbE","colab_type":"code","colab":{},"outputId":"56fdb21f-147a-4b70-cde2-f720a8597048"},"cell_type":"code","source":["print \"Absolute Pressure at 50 deg C = \", Psat(50), \"mmHg\"\n","print \"   Gauge Pressure at 50 deg C = \", (Psat(50)-760)*14.696/760.0, \"psig\"\n"],"execution_count":0,"outputs":[{"output_type":"stream","text":["Absolute Pressure at 50 deg C =  14750.6911081 mmHg\n","   Gauge Pressure at 50 deg C =  270.5357849 psig\n"],"name":"stdout"}]},{"metadata":{"id":"Q-RRmWH9oLbI","colab_type":"code","colab":{},"outputId":"c4c02efc-4438-4ea7-b320-e41d51575d5f"},"cell_type":"code","source":["print 14.696*34.47/1.01325"],"execution_count":0,"outputs":[{"output_type":"stream","text":["499.946824574\n"],"name":"stdout"}]},{"metadata":{"id":"zCrpW7lUoLbL","colab_type":"text"},"cell_type":"markdown","source":["So the minimum operating pressure would be 271 psig in order to keep ammonia in the liquid phase."]},{"metadata":{"id":"vXJlmLj2oLbM","colab_type":"text"},"cell_type":"markdown","source":["#### Part b.\n","The liquid ammonia is delivered to the solar collector at \\unit[50]{$^\\circ$C} and \\unit[500]{psig}. The solar collector contains a pumping unit to bring the pressure  to \\unit[15]{MPa}. What is the pump requirement in watts for a flowrate of 1 kilogram/min of ammonia? \n"]},{"metadata":{"id":"VZmp9kvdoLbO","colab_type":"code","colab":{},"outputId":"0198471f-358c-4300-e99d-6b648321861f"},"cell_type":"code","source":["rho = 0.60  #kg/liter\n","rho = 600.0   # kg/m3\n","\n","Vspecific = 1.0/rho\n","\n","Pin = (500+14.696)*101325/14.696\n","Pout = 15.0*1.0e6\n","\n","W = Vspecific*(Pout-Pin)/60.0\n","print W"],"execution_count":0,"outputs":[{"output_type":"stream","text":["318.091899156\n"],"name":"stdout"}]},{"metadata":{"id":"VIneM3muoLbR","colab_type":"code","colab":{}},"cell_type":"code","source":["dGf = -16600.0\n","dHf = -46150.0\n","\n","Grxn = -2.0*dGf\n","Hrxn = -2.0*dHf\n","\n","T = 600.0 + 273.15\n","lnK = -((Grxn-Hrxn)/298.15 + Hrxn/T)/8.314\n","K = np.exp(lnK)\n","\n"],"execution_count":0,"outputs":[]},{"metadata":{"id":"FL-mSXGqoLbV","colab_type":"text"},"cell_type":"markdown","source":["#### Part c.\n","\n","\n","#### Solution\n","We need to compute the enthalpy change in going from liquid ammonia at 500 psia and 50 deg C to nitrogen and oxygen at 15 MPa absolute and 600 deg C, assuming ideal gas behavior."]},{"metadata":{"id":"oT_wtF5coLbV","colab_type":"code","colab":{},"outputId":"6e17af22-6770-4352-9d29-21cef0438877"},"cell_type":"code","source":["rho = 600.0\n","MW = \n","\n","\n","Hf_NH3 = -46150.0\n","Hf_N2 = 0.0\n","Hf_H2 = 0.\n","\n","Cp_NH3 = 35.6\n","Cp_N2 = 29.1\n","Cp_H2 = 29.1\n","\n","dH_0 = \n","dH_1 = 2.0*Cp_NH3*(25-50)\n","dH_2 = 2.0*Hf_NH3 - Hf_N2 - 3.0*Hf_H2\n","dH_3 = 1.0*Cp_N2*(600-25)\n","dH_4 = 3.0*Cp_H2*(600-25)\n","\n","Hrxn = dH_1 + dH_2 + dH_3 + dH_4\n","print dH_1, dH_2, dH_3, dH_4\n","\n","print \"Net heat of reaction at operating conditions = \", Hrxn, \"Joules/gmole\""],"execution_count":0,"outputs":[{"output_type":"stream","text":["-1780.0 -92300.0 16732.5 50197.5\n","-27150.0\n"],"name":"stdout"}]},{"metadata":{"id":"h3ypgnkfoLbZ","colab_type":"code","colab":{},"outputId":"d5287af0-9c51-4c3e-8a9a-a9d0e41241d4"},"cell_type":"code","source":["15.0*14.696"],"execution_count":0,"outputs":[{"output_type":"execute_result","data":{"text/plain":["2175.5736491487787"]},"metadata":{"tags":[]},"execution_count":15}]},{"metadata":{"id":"asVWdKJaoLbc","colab_type":"code","colab":{}},"cell_type":"code","source":[""],"execution_count":0,"outputs":[]}]}