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! :)

[Mukul Rao]

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]

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! :)

Answer 1

[Mukul Rao]

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.