清空筛选条件
清空筛选条件

the degradation of ozone in aqueous solutions (water). When I plot the all three components, the plot stays at the 0 value for all components.

1 次查看（过去 30 天）
MERT
MERT 2024-5-27
回答： Torsten 2024-5-27
function ProjectOzone
clc
ki = 7e-3;
kj = 2e-5;
k1 = 1e3;
k2 = 1e3;
k3 = 5e7;
OH = 1e-6;
%constants
t0=0;
tn=1000000;
% time iteration interval is set to a high value to reach steady state
Ca0 = 0.002292;
Cb0 = 0;
Cc0 = 0;
% initial concentrations
[tsoln, Csoln] = ode45(@Ozone,[t0 tn],[Ca0 Cb0 Cc0]);
plot(tsoln, Csoln)
xlabel('Time')
ylabel('Concentration')
legend('Cb')
function dC = Ozone(~, C)
dC = zeros(3,1);
dC(1) = -C(1)*(3*ki + kj*OH + 2*k2*sqrt((C(1)/(2*k3))*(2*ki+kj*OH)));
dC(2) = (2*ki)*C(1) - k1*C(1)*C(2) + k2*C(1)*C(3);
dC(3) = kj*C(1)*OH + k1*C(1)*C(2) - k2*C(1)*C(3) - (2*k3)*(C(3)^2);
end %dC(1) is evaluated by steady state approximation method
end
%HERE IS THE CODE, I POST IT IF YOU WANT TO CHECK

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2024-5-27
Ran in:
There seems to be a singularity in the solution at t approximately 320.
ProjectOzone()
function ProjectOzone
clc
ki = 7e-3;
kj = 2e-5;
k1 = 1e3;
k2 = 1e3;
k3 = 5e7;
OH = 1e-6;
%constants
t0=0;
tn=321;%1000000;
% time iteration interval is set to a high value to reach steady state
Ca0 = 0.002292;
Cb0 = 0;
Cc0 = 0;
% initial concentrations
[tsoln, Csoln] = ode15s(@Ozone,[t0 tn],[Ca0 Cb0 Cc0]);
plot(tsoln, Csoln)
xlabel('Time')
ylabel('Concentration')
legend('Cb')
function dC = Ozone(~, C)
dC = zeros(3,1);
dC(1) = -C(1)*(3*ki + kj*OH + 2*k2*sqrt((C(1)/(2*k3))*(2*ki+kj*OH)));
dC(2) = (2*ki)*C(1) - k1*C(1)*C(2) + k2*C(1)*C(3);
dC(3) = kj*C(1)*OH + k1*C(1)*C(2) - k2*C(1)*C(3) - (2*k3)*(C(3)^2);
end %dC(1) is evaluated by steady state approximation method
end

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by