Problem with boundary conditions for PDEPE

5 ビュー (過去 30 日間)
Shubhanshu Agarwal
Shubhanshu Agarwal 2019 年 4 月 19 日
回答済み: Josh Meyer 2019 年 5 月 1 日
I'm trying to apply following boundary equations for 8 PDEs but the solver is assuming zero flux boundary conditions for all the PDEs at the left boundary.
pr= [-rad1;2-ur(2);0-ur(3);0-ur(4);0-ur(5);0-ur(6);0-ur(7);1-ur(8)];
qr = [((Tf2-Ti)/r);0;0;0;0;0;0;0];
ql = [((Tf2-Ti)/r);0;0;0;0;0;0;0];
pl = [rad2;2-ul(2);0-ul(3);0-ul(4);0-ul(5);0-ul(6);0-ul(7);1-ul(8)];
Is there anything wrong with the way I have specified the boundary conditions?


Josh Meyer
Josh Meyer 2019 年 5 月 1 日
From the doc for pdepe,
"When m > 0 and a = 0, boundedness of the solution near x = 0 requires that the flux f vanish at a = 0. pdepe imposes this boundary condition automatically and it ignores values returned in pl and ql."
So if you are using a value of m greater than 0 (i.e. cylindrical or spherical problem symmetry), then the solver always enforces zero flux on the left boundary no matter what you specify.

その他の回答 (0 件)




Community Treasure Hunt

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

Start Hunting!

Translated by