0

我正在寻找以速率 r 模拟矿物随时间的沉淀。当这种矿物沉淀时，它会降低溶液中 Fe(II) 的浓度，这会降低溶液相对于矿物的饱和状态，从而减慢矿物沉淀的速度。

cla;
clf;

%%
%variables and rate calculation
Fe(1)=0.01*0.26/1000; %mM, * activity coefficient, convert to M
Si=1.8/1000; %mM, convert to M
pH=8;
H=10^(-pH);
k=1.11*10^(-9);
b=1.81;
Ksp=4.68*10^(22);
IAP=(Fe^(3)*Si^(2))/H^(6);
logIAP=log10(IAP);
logKsp=log10(Ksp); 

rate=k*exp(b*(logIAP-logKsp)); %mm Fe Kg-1 hour-1

%% 
%Fe vs time 
dt=0.01;
t=0:dt:10;

for ind=2:length(t)
Fe=Fe(1)-(k*exp(b*(logIAP-logKsp))/1000*t);
end

我最终想要一个 Fe 与时间的关系图。上面的代码不起作用，因为 Fe 依赖于 logIAP，它依赖于上一个时间步的 Fe...我需要 Matlab 迭代计算每个时间步的两者。

任何帮助表示赞赏，谢谢！

4

2 回答 2

0

正如其他人指出的那样，您需要制作Fe一个数组。此外，您需要t在方程式中正确索引，因为这是一个向量。

Fe=zeros(size(t)); % now Fe is an array
Fe(1)=0.01*0.26/1000; % first value assigned
...
Ksp=4.68*10^(22);
IAP=(Fe(1)^(3)*Si^(2))/H^(6);
logIAP=log10(IAP);
logKsp=log10(Ksp); 
....
for ind=2:length(t)

    logIAP=log10((Fe(ind-1)^(3)*Si^(2))/H^(6));
    Fe(ind)=Fe(ind-1)-(k*exp(b*(logIAP-logKsp))/1000*t(ind));
end
于 2017-10-12T14:32:47.430 回答
-1

使用 Fe 作为数组：

Fe=[0.01*0.26/1000]; % Fe is now an array

....

for ind=2:length(t)
    Fe(ind)=Fe(ind-1)-(k*exp(b*(logIAP-logKsp))/1000*t); % We fill Fe
end

然后在 t 的函数中绘制 Fe

于 2017-10-12T14:11:13.150 回答