<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Andrea,<div class=""><br class=""></div><div class="">I have used the 2.5D UHD approach and it seems to work for me. Behavior is as expected, though I guess I have not done an extremely thorough job of testing the angular momentum conservation. Do you have a simple test problem that could be used?</div><div class=""><br class=""></div><div class="">Sean</div><div class=""><br class=""></div><div class=""><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">-----------------------------------------------------------<br class="">Sean M. Couch<br class="">Theoretical Astrophysics Including Relativity<br class="">California Institute of Technology<br class="">Mail Code 350-17<br class="">Pasadena, CA 91125<br class="">(626) 395-4282<br class=""><a href="http://www.tapir.caltech.edu/~smc" class="">www.tapir.caltech.edu/~smc</a></div></div></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On Apr 13, 2015, at 12:42 PM, Andrea Derdzinski <<a href="mailto:aderdzinski@astro.columbia.edu" class="">aderdzinski@astro.columbia.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class="">Hi everyone,<br class=""><br class=""></div>I'm trying to simulate an isothermal disk with in 2.5D cylindrical geometry (r,z) around a point mass (with the computational domain outside of the origin). To my understanding I need to use the unsplit Hydro solver in order to conserve angular momentum (<a href="http://flash.uchicago.edu/pipermail/flash-users/2012-September/001163.html" class="">http://flash.uchicago.edu/pipermail/flash-users/2012-September/001163.html</a>). <br class=""><br class=""></div>There is no viscosity, so this is essentially a rotating fluid. Initially I want the fluid to be steady, with no radial velocity (v_r = 0) and no vertical velocity (v_z = 0), only a rotational phi velocity. I solve for the initial conditions analytically to ensure that there is both radial and vertical support, but at the first time step there seems to be a lack of vertical pressure support and the gas develops a negative vertical velocity. <br class=""><br class=""></div>Analytically the gravitational force is balanced by the pressure gradient force, but this doesn't seem to be the case in the simulation. The issue does not go away with higher resolution. I've checked that the isothermal condition is enforced in the EOS unit, and that the gravitational acceleration is implemented correctly. <br class=""><br class=""></div>Has anyone had success running a simulation in 2.5D cylindrical geometry using the unsplit solver? I'd appreciate hearing about any issues other users have ran into or any advice.<br class=""><br class=""></div>The unsplit setup:<br class=""></div>order = 3<br class=""></div>slopeLimiter = "mc"<br class=""></div>RiemannSolver = "HLLC"<br class=""></div>conserveAngmom = .true.<br class=""></div>use_gravhalfupdate = .true.<br class=""></div>fluxcorrect = .true.<br class=""><div class=""><div class=""><div class=""><br class=""><div class=""><div class=""><div class=""><div class=""> <br class=""><br class=""></div><div class="">Thanks,<br class=""></div><div class="">Andrea<br class=""></div><div class=""><br class=""><div class=""><br class=""><br class=""></div></div></div></div></div></div></div></div></div>
</div></blockquote></div><br class=""></div></body></html>