[petsc-users] Adding a processing on the solution

Jed Brown jedbrown at mcs.anl.gov
Sun May 27 10:09:25 CDT 2012

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>

More information about the petsc-users mailing list