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.<div>
<br></div><div>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.<br>
<br><div class="gmail_quote">On Sun, May 27, 2012 at 9:30 AM, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Developers,<br><br>I want to process each solution iteration of SNES by limiter function (preserving TVD condition for fluid flow). How should I deal with this?<br><br>Thanks,<br>BehZad<br>
</blockquote></div><br></div>