File:Energy density.svg

#!/usr/bin/env python
import matplotlib
matplotlib.use('agg')
 
from matplotlib.font_manager import FontProperties
import pylab
import sys
 
# (name, MJ/Kg, MJ/L, [label-adjustments])
data = (
    ('Liquid Hydrogen', 143, 10.1),
    ('Hydrogen Gas (700 bar)', 143, 5.6),
    ('Hydrogen Gas', 143, 0.01079),
#    ('Beryllium', 67.6, 125.1),
    ('Lithium Borohydride', 65.2, 43.4),
#    ('Boron', 58.9, 137.8, {'valign': 'top'}),
    ('Methane', 55.6, 0.0378),
    ('Natural Gas', 53.6, 10),
    ('LPG Propane', 49.6, 25.3, {'valign': 'center', 'x': 8}),
    ('LPG Butane', 49.1, 27.7, {'valign': 'center', 'x': 7}),
    ('Gasoline', 46.4, 34.2, {'valign': 'center', 'x': 5}),
    ('Diesel', 46.2, 37.3, {'valign': 'center', 'x': 4}),
    ('Polyethylene', 46.3, 42.6, {'valign': 'center', 'x': 7.5}),
#    ('Polypropylene', 46.4, 41.7),
#    ('Gasohol E10', 43.54, 33.72, {'valign': 'center', 'x': -8}),
#    ('Gasohol E85', 33.1, 25.65, {'valign': 'center', 'x': -8}),
    ('Lithium', 43.1, 23.0, {'valign': 'top'}),
    ('Kerosene', 42.8, 33, {'valign': 'top'}),
#    ('Biodiesel', 42.20, 33),
#    ('DMF (2,5-Dimethylfuran)', 42, 37.8),
#    ('Crude oil', 46.3, 37),
    ('Polystyrene', 41.4, 43.5),
    ('Fat Metabolism', 38, 35, {'x': -0.5}),
    ('Butanol', 36.6, 29.2, {'valign': 'top'}),
#    ('Graphite', 32.7, 72.9),
    ('Anthracite', 32.5, 72.4, {'valign': 'top'}),
    ('Silicon', 32.2, 75.1),
    ('Aluminum', 31.0, 83.8),
    ('Ethanol', 30, 24, {'valign': 'center', 'x': 5}),
    ('Polyester', 26.0, 35.6, {'valign': 'top'}),
    ('Magnesium', 24.7, 43.0),
    ('Bituminous', 24, 20, {'valign': 'center', 'x': 6.5}),
    ('Methanol', 19.7, 15.6, {'valign': 'center', 'x': 6}),
    ('Hydrazine', 19.5, 19.3, {'valign': 'center', 'x': -6.5}),
    ('Liquid Ammonia', 18.6, 11.5, {'valign': 'top'}),
#    ('PVC', 18.0, 25.2),
    ('Sugar Metabolism', 17, 26.2),
#    ('Calcium', 15.9, 24.6),
    ('Glucose', 15.55, 23.9, {'valign': 'center', 'x': 5}),
    ('Sodium', 13.3, 12.8), #  (burned to wet NaOH)
#    ('Sodium (burned to dry Na2O)', 9.1, 8.8),
    ('Zinc', 5.3, 38.0, {'valign': 'top'}),
#    ('Teflon', 5.1, 11.2),
    ('Iron', 5.2, 40.68), # (burned to Fe2O3)
#    ('Iron (burned to FeO)', 4.9, 38.2),
 
    # 0.54-0.72, 0.9-1.9 (values used are from the Li-Ion battery article)
    ('Lithium Ion Battery', 0.576, 0.972, {'valign': 'center', 'x': 11}),
)
 
font12 = FontProperties(family=['sans-serif'], size=12)
font10 = FontProperties(family=['sans-serif'], size=10)
font8 = FontProperties(family=['sans-serif'], size=8)
 
pylab.rcParams.update({'figure.figsize': [12, 7.5]})
 
fig=pylab.figure(1)
fig.figurePatch.set_alpha(0.0)
ax=pylab.axes([0.05, 0.05, 0.95 - 0.05, 0.95 - 0.05])
ax.axesPatch.set_alpha(0.0)
 
for v in data:
    # (name, MJ/Kg, MJ/L, label-adjustments)
    label = v[0]
    x = v[1]
    y = v[2]
 
    if len(v) < 4:
        halign = 'center'
        valign = 'bottom'
        x_text = x
        y_text = y + 0.9
    else:
        halign = v[3].get('halign', 'center')
        valign = v[3].get('valign', 'bottom')
        x_text = x + v[3].get('x', 0)
        y_text = y + v[3].get('y', 0)
        if valign == 'bottom':
            y_text += 0.9
        elif valign == 'top':
            y_text -= 0.9
 
    pylab.scatter(x=[x], y=[y])
    pylab.text(x=x_text, y=y_text, s=label,
               horizontalalignment=halign,
               verticalalignment=valign,
               fontproperties=font10)
 
pylab.title('Selected Energy Densities',
            fontproperties=font12)
 
pylab.ylabel('MJ/L', fontproperties=font10)
pylab.xlabel('MJ/Kg', fontproperties=font10)
 
pylab.xticks(fontproperties=font8)
pylab.yticks(fontproperties=font8)
 
ax.xaxis.grid(True)
ax.yaxis.grid(True)
 
pylab.xlim(0, 160)
pylab.ylim(0, 90)
 
pylab.savefig('energy_density.svg')