def Static_Analysis(
InputMethod=opem.Functions.Get_Input,
TestMode=False,
PrintMode=True,
ReportMode=True):
"""
Run Larminie-Dicks static analysis.
:param InputMethod : input function or input test vector
:param TestMode : test mode flag
:type InputMethod : dict or Get_Input function object
:type TestMode:bool
:param PrintMode : print mode control flag (True : print outputs)
:type PrintMode:bool
:param ReportMode : report mode control flag (True : generate report)
:type ReportMode: bool
:return: Result as dict
"""
OutputFile = None
CSVFile = None
Warning1 = False
Warning2 = False
I_Warning = 0
Overall_Params_Max = {}
Overall_Params_Linear = {}
Simulation_Title = "Larminie-Dicks"
try:
if PrintMode:
print("###########")
print(Simulation_Title + "-Model Simulation")
print("###########")
OutputParamsKeys = sorted(OutputParams.keys())
Output_Dict = dict(
zip(OutputParamsKeys, [None] * len(OutputParamsKeys)))
if not TestMode:
Input_Dict = InputMethod(InputParams)
else:
Input_Dict = InputMethod
if PrintMode:
print("Analyzing . . .")
Input_Dict = opem.Functions.filter_alpha(Input_Dict)
Name = Input_Dict["Name"]
if ReportMode:
OutputFile = opem.Functions.Output_Init(
Input_Dict, Simulation_Title, Name)
CSVFile = opem.Functions.CSV_Init(
OutputParamsKeys,
OutputParams,
Simulation_Title,
Name)
HTMLFile = opem.Functions.HTML_Init(Simulation_Title, Name)
IEnd = Input_Dict["i-stop"]
IStep = Input_Dict["i-step"]
Precision = opem.Functions.get_precision(IStep)
[i, IEnd, IStep] = opem.Functions.filter_range(
Input_Dict["i-start"], IEnd, IStep)
I_List = []
Efficiency_List = []
Power_List = []
Vstack_List = []
Power_Thermal_List = []
B = B_Calc(Input_Dict["T"])
while i < IEnd:
try:
I_List.append(i)
Output_Dict["Vcell"] = Vcell_Calc(
E0=Input_Dict["E0"],
i=i,
i_0=Input_Dict["i_0"],
i_n=Input_Dict["i_n"],
i_L=Input_Dict["i_L"],
R_M=Input_Dict["RM"],
B=B,
A=Input_Dict["A"])
[Warning1, I_Warning] = opem.Functions.warning_check_1(
Output_Dict["Vcell"], I_Warning, i, Warning1)
Output_Dict["PEM Efficiency"] = Efficiency_Calc(
Output_Dict["Vcell"])
Output_Dict["Power"] = Power_Calc(Output_Dict["Vcell"], i)
Output_Dict["VStack"] = VStack_Calc(
Input_Dict["N"], Output_Dict["Vcell"])
Output_Dict["Power-Stack"] = PowerStack_Calc(
Output_Dict["Power"], Input_Dict["N"])
Output_Dict["Power-Thermal"] = Power_Thermal_Calc(
VStack=Output_Dict["VStack"], N=Input_Dict["N"], i=i)
Vstack_List.append(Output_Dict["VStack"])
Efficiency_List.append(Output_Dict["PEM Efficiency"])
Power_List.append(Output_Dict["Power-Stack"])
Power_Thermal_List.append(Output_Dict["Power-Thermal"])
if ReportMode:
opem.Functions.Output_Save(
OutputParamsKeys,
Output_Dict,
OutputParams,
i,
OutputFile,
PrintMode)
opem.Functions.CSV_Save(
OutputParamsKeys, Output_Dict, i, CSVFile)
i = opem.Functions.rounder(i + IStep, Precision)
except Exception as e:
print(e)
i = opem.Functions.rounder(i + IStep, Precision)
if ReportMode:
opem.Functions.Output_Save(
OutputParamsKeys,
Output_Dict,
OutputParams,
i,
OutputFile,
PrintMode)
opem.Functions.CSV_Save(
OutputParamsKeys, Output_Dict, i, CSVFile)
[Estimated_V, B0, B1] = opem.Functions.linear_plot(
x=I_List, y=Vstack_List)
Linear_Approx_Params = Linear_Aprox_Params_Calc(B0, B1)
Max_Params = Max_Params_Calc(Power_List, Efficiency_List, Vstack_List)
Power_Total = Power_Total_Calc(Vstack_List, IStep, Input_Dict["N"])
Overall_Params_Linear["Pmax(L-Approx)"] = Linear_Approx_Params[0]
Overall_Params_Linear["V0"] = B0
Overall_Params_Linear["K"] = B1
Overall_Params_Linear["VFC|Pmax(L-Approx)"] = Linear_Approx_Params[1]
Overall_Params_Max["Pmax"] = Max_Params["Max_Power"]
Overall_Params_Max["VFC|Pmax"] = Max_Params["Max_VStack"]
Overall_Params_Max["Efficiency|Pmax"] = Max_Params["Max_EFF"]
Overall_Params_Max["Ptotal(Elec)"] = Power_Total[0]
Overall_Params_Max["Ptotal(Thermal)"] = Power_Total[1]
if ReportMode:
OutputFile.close()
CSVFile.close()
if PrintMode:
print(Report_Message)
opem.Functions.HTML_Desc(
Simulation_Title, Larminiee_Description, HTMLFile)
opem.Functions.HTML_Input_Table(
Input_Dict=Input_Dict,
Input_Params=InputParams,
file=HTMLFile)
opem.Functions.HTML_Overall_Params_Table(
Overall_Params_Max,
Overall_Params_Max_Description,
file=HTMLFile,
header=True)
opem.Functions.HTML_Chart(
x=str(I_List),
y=str(Power_List),
color='rgba(255,99,132,1)',
x_label="I(A)",
y_label="P(W)",
chart_name="Power-Stack",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(
x=str(I_List), y=[
str(Vstack_List), str(Estimated_V)], color=[
'rgba(99,100,255,1)', 'rgb(238, 210, 141)'], x_label="I(A)", y_label="V(V)", chart_name=[
"Voltage-Stack", "Linear-Apx"], size="600px", file=HTMLFile)
opem.Functions.HTML_Overall_Params_Table(
Overall_Params_Linear,
Overall_Params_Linear_Description,
file=HTMLFile,
header=False)
opem.Functions.HTML_Chart(
x=str(I_List),
y=str(Efficiency_List),
color='rgb(255, 0, 255)',
x_label="I(A)",
y_label="EFF",
chart_name="Efficiency",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(list(map(opem.Functions.rounder,
Power_List))),
y=str(Efficiency_List),
color='rgb(238, 210, 141)',
x_label="P(W)",
y_label="EFF",
chart_name="Efficiency vs Power",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(
x=str(I_List),
y=str(Power_Thermal_List),
color='rgb(255, 0, 255)',
x_label="I(A)",
y_label="P(W)",
chart_name="Power(Thermal)",
size="600px",
file=HTMLFile)
opem.Functions.warning_print(
warning_flag_1=Warning1,
warning_flag_2=Warning2,
I_Warning=I_Warning,
file=HTMLFile,
PrintMode=PrintMode)
opem.Functions.HTML_End(HTMLFile)
HTMLFile.close()
if PrintMode:
print("Done!")
if not TestMode:
if PrintMode:
print(
"Result In -->" +
os.path.join(
os.getcwd(),
Simulation_Title))
else:
return {
"Status": True,
"P": Power_List,
"I": I_List,
"V": Vstack_List,
"EFF": Efficiency_List,
"Ph": Power_Thermal_List,
"V0": B0,
"K": B1,
"VE": Estimated_V}
except Exception:
if TestMode:
return {
"Status": False,
"Message": "[Error] " +
Simulation_Title +
" Simulation Failed!(Check Your Inputs)"}
print(
"[Error] " +
Simulation_Title +
" Simulation Failed!(Check Your Inputs)")
def Dynamic_Analysis(InputMethod=opem.Functions.Get_Input,
TestMode=False,
PrintMode=True,
ReportMode=True):
"""
Run Padulles-Amphlett analysis.
:param InputMethod : input function or input test vector
:param TestMode : test mode flag
:type InputMethod : dict or Get_Input function object
:type TestMode:bool
:param PrintMode : print mode control flag (True : print outputs)
:type PrintMode:bool
:param ReportMode : report mode control flag (True : generate report)
:type ReportMode: bool
:return: result as dict
"""
OutputFile = None
CSVFile = None
Warning1 = False
Warning2 = False
I_Warning = 0
Overall_Params_Max = {}
Overall_Params_Linear = {}
Simulation_Title = "Padulles-Amphlett"
try:
if PrintMode:
print("###########")
print(Simulation_Title + "-Model Simulation")
print("###########")
OutputParamsKeys = sorted(OutputParams.keys())
Output_Dict = dict(
zip(OutputParamsKeys, [None] * len(OutputParamsKeys)))
if not TestMode:
Input_Dict = InputMethod(InputParams, params_default=Defaults)
else:
Input_Dict = InputMethod
Input_Dict = opem.Functions.filter_default(input_dict=Input_Dict,
params_default=Defaults)
Input_Dict = opem.Functions.filter_lambda(Input_Dict)
if PrintMode:
print("Analyzing . . .")
Name = Input_Dict["Name"]
if ReportMode:
OutputFile = opem.Functions.Output_Init(Input_Dict,
Simulation_Title, Name)
CSVFile = opem.Functions.CSV_Init(OutputParamsKeys, OutputParams,
Simulation_Title, Name)
HTMLFile = opem.Functions.HTML_Init(Simulation_Title, Name)
IEndMax = Input_Dict["JMax"] * Input_Dict["A"]
IEnd = min(IEndMax, Input_Dict["i-stop"])
IEnd = Input_Dict["i-stop"]
IStep = Input_Dict["i-step"]
Precision = opem.Functions.get_precision(IStep)
[i, IEnd, IStep] = opem.Functions.filter_range(Input_Dict["i-start"],
IEnd, IStep)
I_List = []
Power_List = []
Vstack_List = []
Efficiency_List = []
PH2_List = []
PO2_List = []
PH2O_List = []
Eta_Ohmic_List = []
Eta_Conc_List = []
Eta_Active_List = []
Power_Thermal_List = []
Kr = Kr_Calc(Input_Dict["N0"])
qH2 = qH2_Calc(Input_Dict["qMethanol"], Input_Dict["CV"],
Input_Dict["t1"], Input_Dict["t2"])
qO2 = qO2_Calc(qH2, Input_Dict["rho"])
while i < IEnd:
try:
I_List.append(i)
Output_Dict["PO2"] = PO2_Calc(Input_Dict["KO2"],
Input_Dict["tO2"], Kr, i, qO2)
Output_Dict["PH2"] = PH2_Calc(Input_Dict["KH2"],
Input_Dict["tH2"], Kr, i, qH2)
PH2_List.append(Output_Dict["PH2"])
PO2_List.append(Output_Dict["PO2"])
Output_Dict["Eta Activation"] = Eta_Act_Calc(
Input_Dict["T"], Output_Dict["PO2"], Output_Dict["PH2"], i,
Input_Dict["A"])
Eta_Active_List.append(Output_Dict["Eta Activation"])
Output_Dict["Eta Ohmic"] = Eta_Ohmic_Calc(
i,
Input_Dict["l"],
Input_Dict["A"],
Input_Dict["T"],
Input_Dict["lambda"],
R_elec=Input_Dict["R"])
Eta_Ohmic_List.append(Output_Dict["Eta Ohmic"])
Output_Dict["Eta Concentration"] = Eta_Conc_Calc(
i, Input_Dict["A"], Input_Dict["B"], Input_Dict["JMax"])
Eta_Conc_List.append(Output_Dict["Eta Concentration"])
Output_Dict["Loss"] = Loss_Calc(
Output_Dict["Eta Activation"], Output_Dict["Eta Ohmic"],
Output_Dict["Eta Concentration"])
Output_Dict["PH2O"] = PH2O_Calc(Input_Dict["KH2O"],
Input_Dict["tH2O"], Kr, i, qH2)
PH2O_List.append(Output_Dict["PH2O"])
Output_Dict["E"] = Enernst_Calc(Input_Dict["E0"],
Input_Dict["N0"],
Input_Dict["T"],
Output_Dict["PH2"],
Output_Dict["PO2"],
Output_Dict["PH2O"])
Output_Dict["FC Voltage"] = Vcell_Calc(Output_Dict["E"],
Output_Dict["Loss"],
Input_Dict["N0"])
[Warning1, I_Warning
] = opem.Functions.warning_check_1(Output_Dict["FC Voltage"],
I_Warning, i, Warning1)
Warning2 = opem.Functions.warning_check_2(
Vcell=Output_Dict["FC Voltage"], warning_flag=Warning2)
Vstack_List.append(Output_Dict["FC Voltage"])
Output_Dict["FC Efficiency"] = Efficiency_Calc(
Output_Dict["FC Voltage"], Input_Dict["N0"])
Efficiency_List.append(Output_Dict["FC Efficiency"])
Output_Dict["FC Power"] = Power_Calc(Output_Dict["FC Voltage"],
i)
Output_Dict["Power-Thermal"] = Power_Thermal_Calc(
VStack=Output_Dict["FC Voltage"], N=Input_Dict["N0"], i=i)
Power_List.append(Output_Dict["FC Power"])
Power_Thermal_List.append(Output_Dict["Power-Thermal"])
if ReportMode:
opem.Functions.Output_Save(OutputParamsKeys, Output_Dict,
OutputParams, i, OutputFile,
PrintMode)
opem.Functions.CSV_Save(OutputParamsKeys, Output_Dict, i,
CSVFile)
i = opem.Functions.rounder(i + IStep, Precision)
except Exception as e:
print(str(e))
i = opem.Functions.rounder(i + IStep, Precision)
if ReportMode:
opem.Functions.Output_Save(OutputParamsKeys, Output_Dict,
OutputParams, i, OutputFile,
PrintMode)
opem.Functions.CSV_Save(OutputParamsKeys, Output_Dict, i,
CSVFile)
[Estimated_V, B0, B1] = opem.Functions.linear_plot(x=I_List,
y=Vstack_List)
Linear_Approx_Params = Linear_Aprox_Params_Calc(B0, B1)
Max_Params = Max_Params_Calc(Power_List, Efficiency_List, Vstack_List)
Power_Total = Power_Total_Calc(Vstack_List, IStep, Input_Dict["N0"])
Overall_Params_Linear["Pmax(L-Approx)"] = Linear_Approx_Params[0]
Overall_Params_Linear["V0"] = B0
Overall_Params_Linear["K"] = B1
Overall_Params_Linear["VFC|Pmax(L-Approx)"] = Linear_Approx_Params[1]
Overall_Params_Max["Pmax"] = Max_Params["Max_Power"]
Overall_Params_Max["VFC|Pmax"] = Max_Params["Max_VStack"]
Overall_Params_Max["Efficiency|Pmax"] = Max_Params["Max_EFF"]
Overall_Params_Max["Ptotal(Elec)"] = Power_Total[0]
Overall_Params_Max["Ptotal(Thermal)"] = Power_Total[1]
if ReportMode:
OutputFile.close()
CSVFile.close()
if PrintMode:
print(Report_Message)
opem.Functions.HTML_Desc(Simulation_Title,
Padulles_Amphlett_Description, HTMLFile)
opem.Functions.HTML_Input_Table(Input_Dict=Input_Dict,
Input_Params=InputParams,
file=HTMLFile)
opem.Functions.HTML_Overall_Params_Table(
Overall_Params_Max,
Overall_Params_Max_Description,
file=HTMLFile,
header=True)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(Power_List),
color='rgba(255,99,132,1)',
x_label="I(A)",
y_label="P(W)",
chart_name="FC-Power",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(
x=str(I_List),
y=[str(Vstack_List), str(Estimated_V)],
color=['rgba(99,100,255,1)', 'rgb(238, 210, 141)'],
x_label="I(A)",
y_label="V(V)",
chart_name=["FC-Voltage", "Linear-Apx"],
size="600px",
file=HTMLFile)
opem.Functions.HTML_Overall_Params_Table(
Overall_Params_Linear,
Overall_Params_Linear_Description,
file=HTMLFile,
header=False)
opem.Functions.HTML_Chart(
x=str(I_List),
y=[
str(Eta_Active_List),
str(Eta_Conc_List),
str(Eta_Ohmic_List)
],
color=[
'rgba(255,99,132,1)', 'rgba(99,100,255,1)',
'rgb(238, 210, 141)'
],
x_label="I(A)",
y_label="V(V)",
chart_name=["Eta Active", "Eta Conc", "Eta Ohmic"],
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(Efficiency_List),
color='rgb(255, 0, 255)',
x_label="I(A)",
y_label="EFF",
chart_name="Efficiency",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(PO2_List),
color=' rgb(0, 255, 128)',
x_label="I(A)",
y_label="PO2(atm)",
chart_name="PO2",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(PH2_List),
color=' rgb(128, 0, 255)',
x_label="I(A)",
y_label="PH2(atm)",
chart_name="PH2",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(PH2O_List),
color=' rgb(165, 185, 112)',
x_label="I(A)",
y_label="PH2O(atm)",
chart_name="PH2O",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(
list(map(opem.Functions.rounder, Power_List))),
y=str(Efficiency_List),
color='rgb(238, 210, 141)',
x_label="P(W)",
y_label="EFF",
chart_name="Efficiency vs Power",
size="600px",
file=HTMLFile)
opem.Functions.HTML_Chart(x=str(I_List),
y=str(Power_Thermal_List),
color='rgb(255, 0, 255)',
x_label="I(A)",
y_label="P(W)",
chart_name="Power(Thermal)",
size="600px",
file=HTMLFile)
opem.Functions.warning_print(warning_flag_1=Warning1,
warning_flag_2=Warning2,
I_Warning=I_Warning,
file=HTMLFile,
PrintMode=PrintMode)
opem.Functions.HTML_End(HTMLFile)
HTMLFile.close()
if PrintMode:
print("Done!")
if not TestMode:
if PrintMode:
print("Result In -->" +
os.path.join(os.getcwd(), Simulation_Title))
else:
return {
"Status": True,
"P": Power_List,
"I": I_List,
"V": Vstack_List,
"EFF": Efficiency_List,
"PO2": PO2_List,
"PH2": PH2_List,
"PH2O": PH2O_List,
"Ph": Power_Thermal_List,
"V0": B0,
"K": B1,
"Eta_Active": Eta_Active_List,
"Eta_Conc": Eta_Conc_List,
"Eta_Ohmic": Eta_Ohmic_List,
"VE": Estimated_V
}
except Exception:
if TestMode:
return {
"Status":
False,
"Message":
"[Error] " + Simulation_Title +
" Simulation Failed!(Check Your Inputs)"
}
print("[Error] " + Simulation_Title +
" Simulation Failed!(Check Your Inputs)")