dolfinx-tutorial
dolfinx-tutorial copied to clipboard
jorgensd
Example for DG poisson
Hi @jorgensd
Are there plans to add a demo for solving the poisson equation using the DG method?
Thanks! Rem
Hi Rem, I made a notebook on this a long time ago (https://github.com/jorgensd/dolfinx-tutorial/compare/main...dokken/dg_tutorial), and it needs some tender love and care to be updated to the main branch, and maybe a review of the mathematics. Maybe @IgorBaratta or @jpdean is willing to help me with this?:)
Sounds good, happy to help. How about DG convection-diffusion though? The convection term doesn't add much difficulty (it's actually easier to understand than the diffusion term), and is often one of the main reasons to use DG. I could update this repo and use it as a starting point. I've also written a tutorial section on DG for convection-diffusion in my thesis which could be used too.
@jorgensd @jpdean Sure a convection diffusion problem would be even more complete! Let me know if I can give you a hand
My personal agenda is to impose discontinuities at interfaces between materials, maybe one day we'll have a tutorial about this too! 😃
My personal agenda is to impose discontinuities at interfaces between materials, maybe one day we'll have a tutorial about this too! smiley
Hopefully when @jpdean merges his MeshView assembly into DOLFINx.
Hopefully when @jpdean merges his MeshView assembly into DOLFINx.
I thought we could "just" specify the jump over some internal boundaries? What is MeshView?
Hopefully when @jpdean merges his MeshView assembly into DOLFINx.
I thought we could "just" specify the jump over some internal boundaries? What is MeshView?
You could, But you would Need to use a DG space for everything. I would like a part mesh Lagrange (continuous) coupled with a continuous solution on another subdomain (with discontinuity between them).
@jorgensd @jpdean I think I'm almost there with my DG case, just having a hard time finding a tutorial for getting the weak form of the interface condition.
Do you think this is something that could be added to a tutorial too?