Trying to produce non-negative value
조회 수: 2 (최근 30일)
이전 댓글 표시
%Variables
T = 1073; %Temperature
E_thermo = 1; %Thermodynamic voltage
x_H2_Ia = 0.95;%Hydrogyen inlet mole fraction
x_O2_Id = 0.21; %Oxygen inlet mole fraction
p_c = 1;%cathode pressure
p_a = 1;%anode pressure
D_h2 = 1e-4;%Effective hydrogyen diffusivity
D_O2 = 2e-5;%Feective oxygen diffusivity
a = 0.5; %transfer coefficent
j_0 = 0.1; %Exchange current density
A_SOFC = 9e7; %lectrolyte constant
G_act = 100000; %electrolyte activation energy
t_m = 0.00002; %electrolyte thickness
t_a = 50;%anode thickness
t_c = 0.0008;%cathode thickess
R = 8.314;%Gas constant
F = 96485;%Faraday constant
j = 5000; %Current density
%-------------------------------------------------------------------------
%ohmic voltage loss
n_ohmic = j * ((t_m*T)/(A_SOFC*exp(-(G_act/(R*T)))));
%cathode overvoltage
n_cathode = ((R*T)/(4*a*F))*log(j/((j_0*p_c)*((x_O2_Id - t_c)*((j*R*T)/4*F*p_c*D_O2))));
%V = operating voltage of fuel cell
V = E_thermo - n_ohmic - n_cathode ;
The value for n_cathode is meant to be 0.158V
댓글 수: 3
David Hill
2022년 2월 9일
No, but your equation is wrong
n_cathode = R*T/(4*a*F)*log(j/(j_0*p_c*(x_O2_Id - t_c*j*R*T/(4*F*p_c*D_O2))));%division is only for t_c*j*R*t not (x_O2_Id - t_c*j*R*T)
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Thermal Analysis에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!