Slope Limiters in FEniCS

Question

Hello all,

I have a general question - is there any literature/examples on implementing slope limiters in FEniCS. I looked for this, but couldn't find any. I'd like to use slope limiters to solve a hyperbolic pde (1st order transport equation) to address the oscillations that I'm currently getting. Any pointers would help. I could also provide a detailed example of what I'm looking to do.

Thanks,

Comments

The library timestepping has an example with a Slope limiter: https://bitbucket.org/dolfin-adjoint/dolfin-adjoint/src/ae7a40f983e930779e61b71226c9c86405cb6ae6/timestepping/tests/long/advection_2d_limiter_rk3?at=master&fileviewer=file-view-default

They use a embedded C++ code to postprocess the solution and apply the slope limiter. It's the only application I have seen in FEniCS. I tested it and it doesn't work very well, I believe the implementation is not complete. Although this can give you an idea about how to implement a particular slope limiter, i.e.: write C++ code to postprocess the solution.

---

Thank you for replying - This is exactly what I was looking forward - it'll get me started, and hopefully, I can make it work. If you come across any other examples, please share them. Thank you again :)