<br><br><div class="gmail_quote">On Wed, Dec 14, 2011 at 9:35 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div class="gmail_quote"><div class="im">On Wed, Dec 14, 2011 at 19:31, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div>If you are truly near 0 in temperature or pressure and differencing can push you over the limit then I submit it is perfectly reasonable to solve them as VI.<br></div></blockquote><div><br></div></div><div>Note that this tends to grow the algebraic systems quite a bit. Usually such problems are solved in conservative variables (density, momentum, energy), but the more troublesome constraints (pressure, temperature) are nonlinear functions of the conservative variables. As currently implemented, you would have to solve (density, momentum, energy, pressure, temperature) in order to impose those constraints with SNESVI. I wonder if there is a way to eliminate these constraints earlier so that the systems don't grow so much.</div>

</div></blockquote><div>The trouble is that the constraints can get "eliminated" only when they become active. </div><div><br></div></div>