feat : Padulles1 matlab file added #119

Documents/MATLAB/Padulles1.m

@@ -0,0 +1,92 @@
+close all
+
+try
+    python_version = pyversion;
+    fprintf(2,'** Python Version : %s\n',python_version);
+catch e
+    fprintf(2,'** Error : %s\n',e.message);
+end
+
+% Import model
+opem = py.importlib.import_module('opem');
+model = opem.Dynamic.Padulles1;
+test_vector = opem.Params.Padulles_Standard_Vector;
+
+% Model inputs
+test_vector{'T'} = 343; % Fuel cell temperature [K]
+test_vector{'E0'} = 0.6; % No load voltage [V]
+test_vector{'N0'} = 88; % Number of cells
+test_vector{'KO2'} = 0.0000211; % Oxygen valve constant [kmol.s^(-1).atm^(-1)]
+test_vector{'KH2'} = 0.0000422; % Hydrogen valve constant [kmol.s^(-1).atm^(-1)]
+test_vector{'tH2'} = 3.37; % Hydrogen time constant [s]
+test_vector{'tO2'} = 6.74; % Oxygen time constant [s]
+test_vector{'B'} = 0.04777; % Activation voltage constant [V]
+test_vector{'C'} = 0.0136; % Activation constant parameter [A^(-1)]
+test_vector{'Rint'} = 0.00303; % Fuel cell internal resistance [ohm]
+test_vector{'rho'} = 1.168; % Hydrogen-Oxygen flow rate
+test_vector{'qH2'} = 0.0004; % Molar flow of hydrogen [kmol.s^(-1)]
+test_vector{'i-start'} = 0; % Cell operating current start point [A]
+test_vector{'i-stop'} = 100; % Cell operating current end point [A]
+test_vector{'i-step'} = 0.1; % Cell operating current step
+test_vector{'Name'} = 'PadullesI_Test';
+
+% Run simulation
+test_mode = true;
+print_mode = true;
+report_mode = false;
+result = model.Dynamic_Analysis(test_vector,test_mode,print_mode,report_mode);
+
+% Model outputs
+P = cellfun(@vector_filter, cell(result{'P'}));
+I = cellfun(@vector_filter, cell(result{'I'}));
+V = cellfun(@vector_filter, cell(result{'V'}));
+EFF = cellfun(@vector_filter, cell(result{'EFF'}));
+Ph = cellfun(@vector_filter, cell(result{'Ph'}));
+VE = cellfun(@vector_filter, cell(result{'VE'}));
+PO2 = cellfun(@vector_filter, cell(result{'PO2'}));
+PH2 = cellfun(@vector_filter, cell(result{'PH2'}));
+
+% Power-Stack
+figure(1)
+plot(I,P)
+xlabel('I(A)')
+ylabel('P(W)')
+legend('Power-Stack')
+
+% Voltage-Stack & Linear Approximation
+figure(2)
+plot(I,V)
+xlabel('I(A)')
+ylabel('V(V)')
+hold on
+plot(I,VE)
+legend('Voltage-Stack','Linear-Apx')
+
+% Efficiency
+figure(3)
+plot(I,EFF)
+xlabel('I(A)')
+ylabel('EFF')
+legend('Efficiency')
+
+% Power-Thermal
+figure(4)
+plot(I,Ph)
+xlabel('I(A)')
+ylabel('P(W)')
+legend('Power(Thermal)')
+
+% PO2
+figure(5)
+plot(I,PO2)
+xlabel('I(A)')
+ylabel('PO2(atm)')
+legend('PO2')
+
+% PH2
+figure(6)
+plot(I,PH2)
+xlabel('I(A)')
+ylabel('PH2(atm)')
+legend('PH2')
+