[petsc-users] FVM using dmplex for linear advection

[Praveen C]

Dear all

I am trying to write a simple first order upwind FVM to solve linear advection in 2d on triangular meshes using dmplex.

Based on some suggestions in recent discussions, I cooked up this simple example

https://github.com/cpraveen/cfdlab/tree/master/petsc/dmplex_convect2d

A gaussian profile must rotate around origin but in my code, there is some errors developing in just one or two cells, see picture.

I am looking for any suggestions/hints/tips on how to identify the problem. Since I have never written a dmplex code before, I am not able to ask a specific question.

Thanks
praveen

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190828/e124c939/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0000.png
Type: image/png
Size: 117214 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190828/e124c939/attachment-0001.png>