Skip to content

Instantly share code, notes, and snippets.

@vashu1

vashu1/flame_licker_efficiency.py

Last active November 30, 2022 09:55
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save vashu1/1bbe573f7e22ddf463b11e1fb4d4d088 to your computer and use it in GitHub Desktop.
Save vashu1/1bbe573f7e22ddf463b11e1fb4d4d088 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
flame_licker_efficiency.py
""" Output:
T Efficiency
100 1.6
200 5.4
300 7.8
400 9.6
500 11.0
600 12.2
700 13.1
"""
def efficiency(extra_t):
T0 = 300
T1 = T0 + extra_t
Tend = 273+100
p_air = 1.2
E = p_air * (T1 - T0) * 1000
V = T1 / T0 # m3
P = 1 - T0 / T1
work = 0
for t in range(T1-1, Tend, -1):
new_V = t / T0
new_P = 1 - T0 / t
work += (V - new_V) * P * 1e5
V = new_V
P = new_P
return round(work / E*100,1)
for i in range(100, 700+1, 100):
print(i, efficiency(i))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment