Axial Dispersion Model, How to solve the following set of equations
2 次查看(过去 30 天)
显示 更早的评论
here dq/dt = k*(q*-q)
q* = HC
H= henry constant
Initial conditions
at t=0
z=0 to l
q=0
c=0
boundary condition
at t>0
z=0
c=2
How can i solve the following equation, i am stuck with the above problem
0 个评论
回答(1 个)
Torsten
2019-6-26
Use MATLAB's "pdepe".
As boundary condition for q, set
pl = 0, ql = 1, pr = 0, qr = 1.
4 个评论
Torsten
2019-6-26
m = 2;
Why ? Do you work in a sphere ?
function [u,f,s,p] = pdex4pde(x,t,c,DcDx)
Why 4 return parameters ? pdepe needs 3.
u = 1;
f = 0.006.* DcDx;
s = -1.061;
u, f and s must be vectors of length 2 since you have equations for C and q.
function [c0,q] = pdex4ic(x)
Why 2 return parameters ? pdepe needs one array of length 2.
pl = cl;
ql = 2;
pr = 0;
qr = 0;
These settings don't make sense at all.
My advice:
Read the documentation of "pdepe".
Hasti
2020-9-30
Hi Torsten,
In another question here
I see that you commented ""pdepe" is not designed to solve mixtures of partial differential equations and ordinary differential equations (the equation for u is a partial differential equation, the equation for v is an ordinary differential equation)."
In my idea these two cases are similar, so I am wondering if such a system can be solved with pdepe?
Appreciate your reply.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!