Axial Dispersion Model, How to solve the following set of equations
7 ビュー (過去 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!