이 질문을 팔로우합니다.
- 팔로우하는 게시물 피드에서 업데이트를 확인할 수 있습니다.
- 정보 수신 기본 설정에 따라 이메일을 받을 수 있습니다.
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!오류 발생
페이지가 변경되었기 때문에 동작을 완료할 수 없습니다. 업데이트된 상태를 보려면 페이지를 다시 불러오십시오.
웹사이트 선택
번역된 콘텐츠를 보고 지역별 이벤트와 혜택을 살펴보려면 웹사이트를 선택하십시오. 현재 계신 지역에 따라 다음 웹사이트를 권장합니다:
또한 다음 목록에서 웹사이트를 선택하실 수도 있습니다.
사이트 성능 최적화 방법
최고의 사이트 성능을 위해 중국 사이트(중국어 또는 영어)를 선택하십시오. 현재 계신 지역에서는 다른 국가의 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 (한국어)