現在この質問をフォロー中です
- フォローしているコンテンツ フィードに更新が表示されます。
- コミュニケーション基本設定に応じて電子メールを受け取ることができます。
why does pdepe adopt Petrov-Galerkin?
2 ビュー (過去 30 日間)
古いコメントを表示
feynman feynman
2024 年 2 月 13 日
pdepe is meant to solve parabolic and elliptic PDEs. Petrov-Galerkin seems to be designed to solve convection dominated ones, why would matlab use this algorithm for solving non-hyperbolic ones?
19 件のコメント
Torsten
2024 年 2 月 13 日
Address your question to the creators of "pdepe". Most probably, they found in their tests that it worked well for the problem class. Or they wanted to use a non-standard method to have a basis for a new publication.
Torsten
2024 年 2 月 13 日
編集済み: Torsten
2024 年 2 月 13 日
creators of pdepe are from mathworks, so I asked here :)
No. Responsible for the theoretical approach are the authors of the article:
[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.
TMW won't be able to help in this respect.
feynman feynman
2024 年 2 月 13 日
right, I meant there must be good reason for mathworks to decide to adopt their Petrov-Galerkin?
feynman feynman
2024 年 2 月 13 日
it just makes one wonder why pdepe and solvepde adopt different FEMs (assuming solvepde adopts a more regular FEM)
Torsten
2024 年 2 月 13 日
Different university chairs propagate different methods - just consider all the methods used in codes for ordinary differential equations (BDF, Runge-Kutta, Extrapolation, Multistep,...).
"pdepe" is a stand-alone program to solve parabolic-elliptic PDEs in one spatial dimension. So why should the method used not differ from the one of the PDE Toolbox (or to whatever "solvepde" belongs) ?
feynman feynman
2024 年 2 月 15 日
編集済み: feynman feynman
2024 年 2 月 15 日
makes sense, thanks! I wonder if there's any tests on pdepe regarding if it's dissipative or non-dissipative for conservative PDEs (though it's not designed for hyperbolic ones)?
feynman feynman
2024 年 2 月 15 日
It actually solves some hyperbolic PDEs correctly but am just not sure in which situations it can't solve well
Torsten
2024 年 2 月 15 日
編集済み: Torsten
2024 年 2 月 15 日
It cannot solve hyperbolic pdes because the f-term should not equal 0. Further, you have to specify two boundary conditions for each equation, and for equations of hyperbolic type you need only one. So one boundary condition will be wrong or at most artificial.
Why do you want a code force to solve hyperbolic equations if its name already indicates that it is created for the parabolic-elliptic type ?
If you are in need to solve hyperbolic PDEs, use CLAWPACK:
feynman feynman
2024 年 2 月 16 日
編集済み: feynman feynman
2024 年 2 月 16 日
Many thanks for the suggestion of CLAWPACK, which I'll check out. For pdepe, I don't know if periodic boundary conditions are allowed but if so PDEs having only first order spatial derivatives can also be solved. I think other finite element software packages can't solve first order hyperbolic ones well either, except when periodic boundary conditions are used.
Torsten
2024 年 2 月 16 日
編集済み: Torsten
2024 年 2 月 16 日
For pdepe, I don't know if periodic boundary conditions are allowed but if so PDEs having only first order spatial derivatives can also be solved.
Periodic boundary conditions are not possible with pdepe.
As said, setting up a problem with only first-order derivatives is technically possible. But you have to assume a second boundary condition that is mathematically incorrect. And usually - because the first derivative is in essence approximated by a central difference quotient - the results won't be stable.
Torsten
2024 年 2 月 26 日
Why not ? But it's not commercial - so it works with FORTRAN or python, not MATLAB.
feynman feynman
2024 年 2 月 26 日
but their Installation Prerequisites says it only works in
- Linux
- Mac OS X
回答 (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!エラーが発生しました
ページに変更が加えられたため、アクションを完了できません。ページを再度読み込み、更新された状態を確認してください。
Web サイトの選択
Web サイトを選択すると、翻訳されたコンテンツにアクセスし、地域のイベントやサービスを確認できます。現在の位置情報に基づき、次のサイトの選択を推奨します:
また、以下のリストから Web サイトを選択することもできます。
最適なサイトパフォーマンスの取得方法
中国のサイト (中国語または英語) を選択することで、最適なサイトパフォーマンスが得られます。その他の国の 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 (한국어)