examples/notebooks/chp-fixed-heat-power-ratio.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Multi-Link usage: CHP with fixed heat-power ratio\n",
    "\n",
    "This notebook if a demonstration of Link with multiple outputs: Combined-Heat-and-Power (CHP) with fixed heat-power ratio\n",
    "\n",
    "For a CHP with a more complicated heat-power feasible operational area, see the [power-to-gas example](https://github.com/PyPSA/PyPSA/blob/master/examples/notebooks/power-to-gas-boiler-chp.ipynb)\n",
    "\n",
    "This example demonstrates a Link component with more than one bus output (\"bus2\" in this case). In general links can have many output buses.\n",
    "\n",
    "In this example a CHP must be heat-following because there is no other supply of heat to the bus \"Frankfurt heat\"."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pypsa"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network = pypsa.Network()\n",
    "\n",
    "network.add(\"Bus\", \"Frankfurt\", carrier=\"AC\")\n",
    "network.add(\"Load\", \"Frankfurt\", bus=\"Frankfurt\", p_set=5)\n",
    "\n",
    "\n",
    "network.add(\"Bus\", \"Frankfurt heat\", carrier=\"heat\")\n",
    "network.add(\"Load\", \"Frankfurt heat\", bus=\"Frankfurt heat\", p_set=3)\n",
    "\n",
    "\n",
    "network.add(\"Bus\", \"Frankfurt gas\", carrier=\"gas\")\n",
    "network.add(\"Store\", \"Frankfurt gas\", e_initial=1e6, e_nom=1e6, bus=\"Frankfurt gas\")\n",
    "\n",
    "network.add(\n",
    "    \"Link\",\n",
    "    \"OCGT\",\n",
    "    bus0=\"Frankfurt gas\",\n",
    "    bus1=\"Frankfurt\",\n",
    "    p_nom_extendable=True,\n",
    "    capital_cost=600,\n",
    "    efficiency=0.4,\n",
    ")\n",
    "\n",
    "\n",
    "network.add(\n",
    "    \"Link\",\n",
    "    \"CHP\",\n",
    "    bus0=\"Frankfurt gas\",\n",
    "    bus1=\"Frankfurt\",\n",
    "    bus2=\"Frankfurt heat\",\n",
    "    p_nom_extendable=True,\n",
    "    capital_cost=1400,\n",
    "    efficiency=0.3,\n",
    "    efficiency2=0.3,\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network.optimize();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network.loads_t.p"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network.links_t.p0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network.links_t.p1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "network.links_t.p2"
   ]
  }
 ],
 "metadata": {
  "@webio": {
   "lastCommId": null,
   "lastKernelId": null
  },
  "kernelspec": {
   "display_name": "",
   "language": "python",
   "name": ""
  },
  "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.10.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}