A naive question regarding MATLABs definition of parabolic/elliptical PDEs (pdepe)
2 次查看(过去 30 天)
显示 更早的评论
This is a shortly put question, and possibly very obvious, but MATLAB says in the documentation of pdepe (here) that is solves 1D parabolic and elliptical PDEs. However, any definition I've come across for PDEs' conic classification holds that the coefficients in the PDE expansion should be constant, i.e.
where
For pdepe these coefficients are functions however. How does MATLAB classify them as parabolic or elliptical if the coefficients are variable? Is it that it must follow the inequality below?
0 个评论
采纳的回答
Torsten
2024-4-28
编辑:Torsten
2024-4-28
Second-order linear partial differential equations (PDEs) are classified as either elliptic, hyperbolic, or parabolic. Any second-order linear PDE in two variables can be written in the form ... where A, B, C, D, E, F, and G are functions of x and y.
Loosely speaking, elliptic for "pdepe" means the equation has the form D*d^2u/dx^2 - v*du/dx + s = 0 and parabolic means it has the form du/dt = D*d^2u/dx^2 + v*du/dx + s.
You shouldn't care so much about dependency of the coefficients. The important thing is that the equation doesn't have the form du/dt =- v*du/dx + s or - v*du/dx + s = 0 (hyperbolic type).
2 个评论
Torsten
2024-4-28
编辑:Torsten
2024-4-28
The presence of the second-derivative term is crucial because "pdepe" expects 2 boundary conditions for each PDE that you implement. If you have no second-order term, no (if no spatial derivatives are present) or only one (if only a first-order spatial derivative is present) boundary condition is mathematically necessary ( and sufficient ) to define a well-posed problem. The remaining boundary conditions that you would have to impose to make "pdepe" work would be artificial ones.
Another aspect why a second-order term is necessary is the discretization scheme used for the first-order spatial derivatives. It usually is not stable if no stabilizing second-order term is present.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Boundary Conditions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!