Spatial discretization in pdepe compared to others
2 次查看(过去 30 天)
显示 更早的评论
The spatial discretization in pdepe seems to be a unique one. Could someone please explain it briefly so that it is understandable in terms of how it differs from other finite differences? Also on the stability and accuracy.
11 个评论
Torsten
2024-1-30
编辑:Torsten
2024-1-30
Did you read the reference for "pdepe" ?
[1] Skeel, R. D. and M. Berzins, "A Method for the Spatial Discretization of Parabolic Equations in One Space Variable," SIAM Journal on Scientific and Statistical Computing, Vol. 11, 1990, pp.1–32.
?
It should be a usual finite-element discretization in 1d, shouldn't it ?
Torsten
2024-1-30
编辑:Torsten
2024-1-30
The authors claim in the introduction that they "propose a simple piecewise nonlinear Galerkin / Petrov-Galerkin method that is second-order accurate in space". These methods are summarized under Finite Element Methods as far as I know. But if you have doubts, you can either read the relevant chapters in a book on the numerical solution of partial differential equations or contact Technical Support to get further details.
feynman feynman
2024-2-1
Thanks. It will be good to get answers from this community because I don't really know how to contact the technical support. If it's a Galerkin / Petrov-Galerkin method, I wonder about the stability and what sorts of hyperbolic equations it can't solve.
Torsten
2024-2-1
编辑:Torsten
2024-2-1
Simply speaking, hyperbolic equations need specific discretization schemes for the flux terms so that the solution remains stable. The schemes used for parabolic-elliptic differential equations don't account for this.
Since there is an interaction of spatial discretization and integration in time with ode45, I think not much about the stability of the integration with "pdepe" can be said.
But why do you need these specific information ? Do you want to compare different pde solvers with respect to their capability in solving pdes of different type ?
feynman feynman
2024-2-1
thanks a lot and this helps. I want to summarize what matlab can do in solving PDEs. The pdepe and solvepde seem to be the two go to methods. Can we conclude that the latter is always better?
Torsten
2024-2-1
编辑:Torsten
2024-2-1
"pdepe" is for one-dimensional, "solvepde" for two- and three-dimensional problems (as far as I can see).
Define "better" !
Broader range of equations that can be solved, more user-friendly, faster, less memory consumption ? I cannot answer this, but "better" is always related to what is required. And since both codes solve problems of different dimension, I don't know exactly how to compare them.
feynman feynman
2024-2-2
I guess solvepde can also solve 1d problems, then in terms of dimensions of solvepde is superior to pdepe.
By better, I mean if solvepde can solve all hyperbolic equations that pdepe can't?
Torsten
2024-2-2
I guess solvepde can also solve 1d problems, then in terms of dimensions of solvepde is superior to pdepe.
If this is the case, my guess is that internally, "pdepe" is called.
By better, I mean if solvepde can solve all hyperbolic equations that pdepe can't?
There is no official MATLAB solver that can solve hyperbolic PDEs.
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Eigenvalue Problems 的更多信息
标签
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!发生错误
由于页面发生更改,无法完成操作。请重新加载页面以查看其更新后的状态。
您也可以从以下列表中选择网站:
如何获得最佳网站性能
选择中国网站(中文或英文)以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。
美洲
- América Latina (Español)
- Canada (English)
- United States (English)
欧洲
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom(English)
亚太
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)