Multicomponent Adsorption | Three Components Breakhthrough Curve Matlab (PDE ODE Simultaneous Solved
Вставка
- Опубліковано 5 вер 2024
- #matlab #adsorption #multicomponents
Multicomponent Adsorption | Three Components Breakhthrough Curve Matlab (PDE ODE Simultaneous Solved)
Hi Rahmat please dont stop uploading at any cost, I have been following your vidoes , they are so useful that they cant be explained in words
My pleasure . Thank you for your support
Sir. Good work. It has been impactful to me. What about pepde solver for two or three components systems with and without axial dispersion. Tks
Hello Sir. Plese visit at following link for pdepe solver axial dispersion ua-cam.com/video/ov8A5P4RVjc/v-deo.html
Hi Rahmat, what if I have a variable velocity, how do I discretize the mass balance for total mixture
You can see at another videos for variable velocity, following link will guide you : ua-cam.com/video/gH_nkOxIAsY/v-deo.html
Have you ever tried to solve with respect velocity change via space? (du/dz is not equal to zero)
Hello. Thank you for your comment. I never tried with unconstant BC in adsorption, but for heat transfer I have tried with following link ua-cam.com/video/gH_nkOxIAsY/v-deo.html Thank You
@@rahmatsunarya Thanks for your contribution. I am deal with the air separation where most of nitrogen adsorbed on the surface that stands for velocity of gas wrt bed will change dramatically. This accounts for du/dz nonequality to zero. Since dudz term can not be integrated directly as you did at mentioned video, in my code, I wrote dudz(i)=(-dcTdt(i)-u*dcTdz(i)-dqTdt(i))./cT(i) in the "interior" section of yours. Eventually, however, the velocity sometimes drops below zero. I wonder how can this be interpreted?
@@kerembulbul8511 Hmm, I see. Maybe you can use the Ergun equation to see the pressure drop in a long bed. I can't imagine that the velocity tends to zero, which means no flow in the system.
@@rahmatsunarya Here is another question; we already calculate change of concentration wrt column (z). From species summation (dA/dz+dB/dz+..=cT/dz) We know total concentration change where can be translated into dP/dz. What does the Ergun equation represent in this calculation?
@@kerembulbul8511 In ideal gas, we assume C=P/RT. so the effect pressure drop a long bed may be considered.
But, I still can't imagine the gas velocity was drop to zero.
I, have you tried with heat balance coupled ?
and with b1 b2, b3 depending on T (with the exp)?
Hello Sara, Thank you for your comment. I have not tried with non-isothermal process.
@@rahmatsunarya thank you for your response because I have a problem with my non-isothermal process do you have time to take a quick look on it please?
@@sarab297 please send me at rahmatsunarya@gmail.com
@@rahmatsunarya Perfect I sent it to you with my student adress
Hy.. I tried putting the same code on matlab and it's not working..
I am actually trying to solve a system of 4PDE in matlab. ( Absolutely new to this)
I feel like this video is very similar to the problem..so I was trying to draw parallels.. unfortunately I couldn't run the code in matlab..
Can you help me out too?..
Sure. Please let me know
@@rahmatsunarya Error using vertcat
Dimensions of arrays being concatenated are not consistent.
Error in self3>f (line 124)
dydt=[dc1dt;dc2dt;dc3dt;dq1dt;dq2dt;dq3dt];
Error in odearguments (line 92)
f0 = ode(t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ode15s (line 148)
odearguments(odeIsFuncHandle, odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin);
Error in self3 (line 36)
[t, y] =ode15s(@f,t,IC,[],eps,u,LDF,b1,b2,b3,qm1,qm2,qm3,c10,c20,c30,rhop,Nz,dz);
>> these kinda errors i am getting .can i contact you via gmail?
@@rahmatsunarya i am getting erros like
Error using vertcat
Dimensions of arrays being concatenated are not consistent.
Error in self3>f (line 124)
dydt=[dc1dt;dc2dt;dc3dt;dq1dt;dq2dt;dq3dt];
Error in odearguments (line 92)
f0 = ode(t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ode15s (line 148)
odearguments(odeIsFuncHandle, odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin);
Error in self3 (line 36)
[t, y] =ode15s(@f,t,IC,[],eps,u,LDF,b1,b2,b3,qm1,qm2,qm3,c10,c20,c30,rhop,Nz,dz);
>> for the same code youre using. Can i contact you via gmail ?
Surely, rahmatsunarya@gmail.com
@@rahmatsunarya I have mailed you
Hello Rahmat, I tried to type the whole of your code into matlab but it’s not working for me :\ would you be able to send me the exact code so I can paste it into matlab for it to work please. Would really appreciate it
Please email me at rahmatsunarya@gmail.com i will send it
@@rahmatsunarya of course. Thankyou so much
Ok. I just replied your email. Please enjoy