Discontinuous Galerkin Method 2D

Hello everyone,
I try to find a discontinuous galerkin method solver of the simple equation :
- div(p(nabla(u))= f on omega u=g on the boundary
Where omega is a square [-1 1]*[-1 1] here with triangular meshes!
Thank you for your help! :)

2 个评论

Hello, can you share why you need the discontinuous Galerkin method for this problem? What kind of discontinuity are you planning to model?
Hello, thank you for your answer,
I forgot to give some extra indications : p=1 if y>=0 ans p=100 if y<0. Also, f= exp(x)-5*exp(x+2y)
I'm planning to model the IP model
Thank you for your help! :)

请先登录,再进行评论。

回答(1 个)

Mukul Rao
Mukul Rao 2017-12-5

0 个投票

Hi Bilal, you could model your domain as two rectangles sharing a common face at y = 0. You can then use the specifyCoefficients function to set different values of parameter "c" ( which corresponds to "p" in your equation), by face value.
See the specify coefficients by sub-domain documentation page for examples.

类别

帮助中心File Exchange 中查找有关 Mathematics 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by