Does it make sense my results of a pde equation from a metalurgy problem?

2 visualizaciones (últimos 30 días)
Raúl Rivera
Raúl Rivera el 18 de Oct. de 2023
Editada: Torsten el 19 de Oct. de 2023
Hello , I have a metalurgy problem in which I want to model a CSTR reactor , but whose equation is a diffusion-convection equation that is being addresed next, by the way I tried to solve it with the pdepe solver. I adjointed the pdf file so can be better understandable.
As assumption I am working only with the gold in ore equation because cyanide and oxygen equation are supposedly constant in time (not real but ideally) ,so the following equation is left :
kinetic rate is:
which "t" means residence time.
Initial conditions are : t=0 ; Cau(z,t)=0
Boundary conditions are:
z=0 ; Cau(0,t)=Cau_0
z=L ; dC/dt =0
The code is the following.
x = 0:1:100;
t = 0:1:100;
m = 0;
sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
u1 = sol(:,:,1);
surf(x,t,u1)
zlabel('Mineral conc (x,t)')
xlabel('Reactor length')
ylabel('Residence Time')
function [c,f,s] = pdefun(x,t,u,dudx) % Equation to solve
U=10;
E=0.1e-8;
c = 1;
k=0.0017;
a=1.3703;
b=1.071;
cc=0.1195;
Ox=8;
cn=159;
f = E.* dudx;% - U.*u;
y =k*((u-0.17)^a)*((cn)^b)*((Ox)^cc) ;
s = -y -U.*dudx ;
end
% ---------------------------------------------
function u0 = pdeic(x) % Initial Conditions
u0 = 1;
end
% ---------------------------------------------
function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,t) % Boundary Conditions
Cai=1;
pl = ul-Cai;
ql = 0;
pr = 0;
qr = 1;
end
% ---------------------------------------------
My question arise from the fact that when using a low mineral concentration feed, my reactor length spans a large length , on the contrary when feeding with a high mineral concentration the residence time and reactor length reduces. Should´t be the other way around? Don´t know if the code is broke or its not properly written or there is a reason thereby. Thanks in advance for the help.
  5 comentarios
Raúl Rivera
Raúl Rivera el 19 de Oct. de 2023
Yes , C_Au^inf does not depend on the concentration C_Au, because that is the concentration after a infinite leaching time , it is took from lab experiments , as much the other exponential values, just the gold concentration is a variable.
Torsten
Torsten el 19 de Oct. de 2023
Editada: Torsten el 19 de Oct. de 2023
Ok, then you must accept from the form of your equation for r_1 that the reaction rate in case of high concentrations of C_Au is much higher than for low concentrations. Thus there is no reason not to believe that the time to reach equilibrium of 0.17 is almost independent of the concentration of C_Au at the inlet.

Iniciar sesión para comentar.

Respuestas (0)

Productos


Versión

R2023a

Community Treasure Hunt

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

Start Hunting!

Translated by