Axial Dispersion Model, How to solve the following set of equations
3 views (last 30 days)
Show older comments
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 Comments
Answers (1)
Torsten
on 26 Jun 2019
Use MATLAB's "pdepe".
As boundary condition for q, set
pl = 0, ql = 1, pr = 0, qr = 1.
4 Comments
Torsten
on 26 Jun 2019
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
on 30 Sep 2020
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.
See Also
Categories
Find more on Mathematics and Optimization in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!