<div class="gmail_quote">On Wed, Dec 14, 2011 at 19:31, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">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 id=":46f">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>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><br></div></div>