Homec4science
Diffusion Convection Heat Transfer Notebook (master) Browse

Convection/Python_Code/ExHydrogen.py
2f22c81e3f9amaster

Last Change
Raw File

ExHydrogen.py
View Options

import numpy as np
import ipywidgets as widgets
import matplotlib.pyplot as plt
from IPython.display import display, clear_output
def EvaluateNumericalSolution(inputstring, solution, tolerance):
notaNumber = "Result is not a number"
correct = "Your answer is correct!"
wrong = "Wrong. Please try again."
clear_output();
try:
result = float(inputstring)
except ValueError:
print(notaNumber)
return
if np.abs(result-solution)<=tolerance:
print(correct)
else:
print(wrong)
"""Function that clears and updates the displayed value"""
def echo(string, padding=40):
padding = " " * (padding - len(string)) if padding else ""
print(string + padding, end='\r')
"""Nusselt correlation passed to students"""
def Nusselt(Re,Pr):
result = 0.3+(0.62*Re**(0.5)*Pr**(1/3))/ ((1+(0.4/Pr)**(2/3))**(0.25))*(1+(Re/282000)**(5/8))**(4/5)
resultstr = "The Nusselt Number is:%.2f \n" % float(result)
return resultstr
"""PLOT FUNCTIONS"""
def plotTemperature():
def Twall(r,r1=0.05,r2=0.055,T1=275.2,Tinterface=275.7):
"""calculates the temperature inside the wall for r1<r2 without heat sources"""
if r<=r1:
T = T1
else:
T = (T1-Tinterface)/np.log(r1/r2)*np.log(r/r2)+Tinterface
return T
rwall = np.linspace(0.045,0.055,100)
Tin =[Twall(x) for x in rwall]
rout = np.linspace(0.055,0.060,50)
Tout = 296*np.ones((len(rout),1))
fig = plt.figure(figsize=(15,7))
plt.rcParams.update({'font.size': 20})
plt.plot(rwall,Tin,color='r')
plt.plot(rout,Tout,color='r')
plt.plot([0.05,0.05],[273,296],'k')
plt.plot([0.055,0.055],[273,296],'k')
plt.ylabel("Temperature")
plt.xlabel("Radius [m]");
plt.grid(True)
plt.annotate("Metal hydride \n- Canister interface",xy=(0.05,285),xytext=(0.045,290),arrowprops=dict(facecolor='black', shrink=0.05))
plt.annotate("Canister - \n Air interface",xy=(0.055,285),xytext=(0.051,290),arrowprops=dict(facecolor='black', shrink=0.05))
return fig
"""Question 4"""
def EvalEnergyGen():
text = widgets.Text(placeholder='Type your answer')
display(text)
checkbutton = widgets.Button(description="Check answer")
cheatbutton = widgets.Button(description="Show solution")
buttons = HBox([checkbutton, heatbutton])
output = widgets.Output()
# output = widgets.Output()
# display(button, output)
solution = -3982
tolerance = 20
solstr = str(solution)
def evaluateEnGen(b):
with output:
if text.value!="":
resultstring = text.value
# solution = -3982
# tolerance = 20
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
def display_solution_EnGen(s):
with output:
clear_output()
display(Latex(r"""The result is $E_{gen}=""" + solstr + """ m^2 K / W$."""))
# button.on_click(evaluateEnGen)
display(buttons, output)
checkbutton.on_click(evaluateEnGen)
cheatbutton.on_click(display_solution_EnGen)
"""Question 5.2"""
def EvalFlowConditions():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateFlowCondition(b):
with output:
if text.value!="":
resultstring = text.value
solution = 173760
tolerance = 300
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateFlowCondition)
"""Question 6.1"""
def NusseltCorr():
def eval_nuCorr(Nusselt):
if Nusselt =="":
print("")
return
elif Nusselt=="A":
answer = True
else:
answer = False
return widgets.Valid(value=answer)
widgets.interact(eval_nuCorr, Nusselt=widgets.Dropdown(
options = ["","A", "B","C","D"]
))
"""Question 6.2"""
def EvalNusselt():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateNusselt(b):
with output:
if text.value!="":
solution = 315.8
resultstring = text.value
tolerance = 1
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateNusselt)
"""Question 6.3"""
def EvalConvectionCoeff():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateConvectionCoeff(b):
with output:
if text.value!="":
solution = 72.6
resultstring = text.value
tolerance = 3
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateConvectionCoeff)
"""Question 7.1"""
def EvalRwall():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateRwall(b):
with output:
if text.value!="":
solution = 0.001415
resultstring = text.value
tolerance = 0.0001
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateRwall)
"""Question 7.2"""
def EvalRconv():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateRconv(b):
with output:
if text.value!="":
solution = 0.053
resultstring = text.value
tolerance = 0.01
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateRconv)
"""Question 7.3"""
def EvalTout():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateRconv(b):
with output:
if text.value!="":
solution = 275.73
resultstring = text.value
tolerance = 0.2
EvaluateNumericalSolution(resultstring, solution, tolerance)
try:
result = float(resultstring)
except ValueError:
return
if abs(solution-result)<=tolerance:
display(plotTemperature())
else:
clear_output();
button.on_click(evaluateRconv)
"""Quesiton 8.1"""
def EvalHeatFlux():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateHeatFlux(b):
with output:
if text.value!="":
solution = 406
resultstring = text.value
tolerance = 3
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateHeatFlux)
"""Question 8.2"""
def EvalAddHeatFlux():
text = widgets.Text(placeholder='Type your answer')
display(text)
button = widgets.Button(description="Check answer")
output = widgets.Output()
display(button, output)
def evaluateAddHeatFlux(b):
with output:
if text.value!="":
solution = 3576
resultstring = text.value
tolerance = 10
EvaluateNumericalSolution(resultstring, solution, tolerance)
else:
clear_output();
button.on_click(evaluateAddHeatFlux)
c4science · Help