Discontinuous Galerkin Method 2D

10 次查看(过去 30 天)
Bilal ALAUDDIN
Bilal ALAUDDIN 2017-12-1
回答: Mukul Rao 2017-12-5
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 个评论
Mukul Rao
Mukul Rao 2017-12-4
Hello, can you share why you need the discontinuous Galerkin method for this problem? What kind of discontinuity are you planning to model?
Bilal ALAUDDIN
Bilal ALAUDDIN 2017-12-5
编辑:Bilal ALAUDDIN 2017-12-5
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
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.

类别

Help CenterFile Exchange 中查找有关 Ordinary Differential Equations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by