TVD methods need to modify the flux function. If you just do an unlimited flux function and then project after each step, you lose local conservation, which if you look in the dictionary, is the definition of throwing the baby out with the bath water.
So put the limiter in the residual evaluation (as slope reconstruction or flux limiting; as you prefer), but try to use a limiter that has a continuous first derivative, otherwise Newton may converge slowly or even stagnate before reaching a tight tolerance. On Sun, May 27, 2012 at 9:30 AM, behzad baghapour < behzad.baghapour at gmail.com> wrote: > Dear Developers, > > I want to process each solution iteration of SNES by limiter function > (preserving TVD condition for fluid flow). How should I deal with this? > > Thanks, > BehZad > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120527/0cf40923/attachment.html>
