Problem with boundary conditions for PDEPE
7 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Shubhanshu Agarwal
am 19 Apr. 2019
Beantwortet: Josh Meyer
am 1 Mai 2019
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?
0 Kommentare
Akzeptierte Antwort
Josh Meyer
am 1 Mai 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.
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu PDE Solvers finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!