Problem with boundary conditions for PDEPE

2 views (last 30 days)
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?

Accepted Answer

Josh Meyer
Josh Meyer on 1 May 2019
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.

More Answers (0)




Community Treasure Hunt

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

Start Hunting!