朗肯动力循环

通过朗肯动力循环的实际示例演示Sankey类。

  1. import matplotlib.pyplot as plt
  3. from matplotlib.sankey import Sankey
  5. fig = plt.figure(figsize=(8, 9))
  6. ax = fig.add_subplot(1, 1, 1, xticks=[], yticks=[],
  7.                      title="Rankine Power Cycle: Example 8.6 from Moran and "
  8.                      "Shapiro\n\x22Fundamentals of Engineering Thermodynamics "
  9.                      "\x22, 6th ed., 2008")
  10. Hdot = [260.431, 35.078, 180.794, 221.115, 22.700,
  11.         142.361, 10.193, 10.210, 43.670, 44.312,
  12.         68.631, 10.758, 10.758, 0.017, 0.642,
  13.         232.121, 44.559, 100.613, 132.168]  # MW
  14. sankey = Sankey(ax=ax, format='%.3G', unit=' MW', gap=0.5, scale=1.0/Hdot[0])
  15. sankey.add(patchlabel='\n\nPump 1', rotation=90, facecolor='#37c959',
  16.            flows=[Hdot[13], Hdot[6], -Hdot[7]],
  17.            labels=['Shaft power', '', None],
  18.            pathlengths=[0.4, 0.883, 0.25],
  19.            orientations=[1, -1, 0])
  20. sankey.add(patchlabel='\n\nOpen\nheater', facecolor='#37c959',
  21.            flows=[Hdot[11], Hdot[7], Hdot[4], -Hdot[8]],
  22.            labels=[None, '', None, None],
  23.            pathlengths=[0.25, 0.25, 1.93, 0.25],
  24.            orientations=[1, 0, -1, 0], prior=0, connect=(2, 1))
  25. sankey.add(patchlabel='\n\nPump 2', facecolor='#37c959',
  26.            flows=[Hdot[14], Hdot[8], -Hdot[9]],
  27.            labels=['Shaft power', '', None],
  28.            pathlengths=[0.4, 0.25, 0.25],
  29.            orientations=[1, 0, 0], prior=1, connect=(3, 1))
  30. sankey.add(patchlabel='Closed\nheater', trunklength=2.914, fc='#37c959',
  31.            flows=[Hdot[9], Hdot[1], -Hdot[11], -Hdot[10]],
  32.            pathlengths=[0.25, 1.543, 0.25, 0.25],
  33.            labels=['', '', None, None],
  34.            orientations=[0, -1, 1, -1], prior=2, connect=(2, 0))
  35. sankey.add(patchlabel='Trap', facecolor='#37c959', trunklength=5.102,
  36.            flows=[Hdot[11], -Hdot[12]],
  37.            labels=['\n', None],
  38.            pathlengths=[1.0, 1.01],
  39.            orientations=[1, 1], prior=3, connect=(2, 0))
  40. sankey.add(patchlabel='Steam\ngenerator', facecolor='#ff5555',
  41.            flows=[Hdot[15], Hdot[10], Hdot[2], -Hdot[3], -Hdot[0]],
  42.            labels=['Heat rate', '', '', None, None],
  43.            pathlengths=0.25,
  44.            orientations=[1, 0, -1, -1, -1], prior=3, connect=(3, 1))
  45. sankey.add(patchlabel='\n\n\nTurbine 1', facecolor='#37c959',
  46.            flows=[Hdot[0], -Hdot[16], -Hdot[1], -Hdot[2]],
  47.            labels=['', None, None, None],
  48.            pathlengths=[0.25, 0.153, 1.543, 0.25],
  49.            orientations=[0, 1, -1, -1], prior=5, connect=(4, 0))
  50. sankey.add(patchlabel='\n\n\nReheat', facecolor='#37c959',
  51.            flows=[Hdot[2], -Hdot[2]],
  52.            labels=[None, None],
  53.            pathlengths=[0.725, 0.25],
  54.            orientations=[-1, 0], prior=6, connect=(3, 0))
  55. sankey.add(patchlabel='Turbine 2', trunklength=3.212, facecolor='#37c959',
  56.            flows=[Hdot[3], Hdot[16], -Hdot[5], -Hdot[4], -Hdot[17]],
  57.            labels=[None, 'Shaft power', None, '', 'Shaft power'],
  58.            pathlengths=[0.751, 0.15, 0.25, 1.93, 0.25],
  59.            orientations=[0, -1, 0, -1, 1], prior=6, connect=(1, 1))
  60. sankey.add(patchlabel='Condenser', facecolor='#58b1fa', trunklength=1.764,
  61.            flows=[Hdot[5], -Hdot[18], -Hdot[6]],
  62.            labels=['', 'Heat rate', None],
  63.            pathlengths=[0.45, 0.25, 0.883],
  64.            orientations=[-1, 1, 0], prior=8, connect=(2, 0))
  65. diagrams = sankey.finish()
  66. for diagram in diagrams:
  67.     diagram.text.set_fontweight('bold')
  68.     diagram.text.set_fontsize('10')
  69.     for text in diagram.texts:
  70.         text.set_fontsize('10')
  71. # Notice that the explicit connections are handled automatically, but the
  72. # implicit ones currently are not.  The lengths of the paths and the trunks
  73. # must be adjusted manually, and that is a bit tricky.
  75. plt.show()

朗肯动力循环示例

参考

此示例中显示了以下函数，方法，类和模块的使用：

  1. import matplotlib
  2. matplotlib.sankey
  3. matplotlib.sankey.Sankey
  4. matplotlib.sankey.Sankey.add
  5. matplotlib.sankey.Sankey.finish

下载这个示例