Trying to produce non-negative value

4 次查看（过去 30 天）
Nathaniel Porter
Nathaniel Porter 2022-2-9
评论： Nathaniel Porter 2022-2-9
%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 个评论
显示 1更早的评论
David Hill
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 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by