<div dir="ltr">Sorry lan Hinder,<div><br></div><div>I miss-understood your words "constraints". I thought you are asking me whether I give the initial data some initial velocity so that it can surround the BH stably (if there is only one BH).</div><div><br></div><div>Now I know your meaning, and I am trying to solve them.</div><div><br></div><div>Thank you.</div><div><br></div><div>Yours,</div><div><br></div><div>Zhi-Chao.</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-08-18 22:55 GMT+08:00 赵志超 <span dir="ltr"><<a href="mailto:yanyuechuixue@gmail.com" target="_blank">yanyuechuixue@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thank you Ian Hinder.<div><br></div><div>At the first step, I just want to add a low density fluid to the BBH simulation (mimic the atmosphere around the BH) to see how it evolution (I remember I had read a paper which describe this kind simulation, and I want to repeat it.). I put the fluid static in the initial data.</div><div><br></div><div>I think I will give some constraints on the fluid static in the initial data in the future. Maybe it need to write a new thorn, I think.</div><div><br></div><div>Yours,</div><div>Zhi-Chao.</div><div><br></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">2017-08-18 22:43 GMT+08:00 Ian Hinder <span dir="ltr"><<a href="mailto:ian.hinder@aei.mpg.de" target="_blank">ian.hinder@aei.mpg.de</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><br><div><span><div>On 18 Aug 2017, at 10:05, 赵志超 <<a href="mailto:yanyuechuixue@gmail.com" target="_blank">yanyuechuixue@gmail.com</a>> wrote:</div><br class="m_-7362986380194826924m_-135720166575567396Apple-interchange-newline"><blockquote type="cite"><div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div dir="ltr">Thank you Ian Hinder.<div><br></div><div>Yes I want to add a low density fluid to a BBH simulation.</div><div><br></div><div>I review the source code of GRHydro_InitialData Thorn's GRHydro_OnlyAtmo.F90 . I found it give any position ( I mean every nx,ny,nz) a rho value which is a constant equal to "rho_abs_min" or "initial_rho_abs_min".</div><div><br></div><div>so I think it can give a low density fluid to a BBH simulation.</div><div><br></div><div>But it report "NAN" or "Con2Prim: count > GRHydro_countmax".</div></div></div></blockquote><div><br></div></span><div>OK, thanks for the clarification! </div><div><br></div><div>I'm not an expert on GRHydro, but maybe someone else can suggest something? </div><div><br></div><div>Are you planning to solve the constraints?</div><div><br></div></div><span><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>-- </div><div>Ian Hinder</div><div><a href="http://members.aei.mpg.de/ianhin" target="_blank">http://members.aei.mpg.de/ianh<wbr>in</a></div></div></div></div></div>
</div>
<br></span></div></blockquote></div><br><br clear="all"><div><br></div></div></div><span class="">-- <br><div class="m_-7362986380194826924gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>赵志超</div><div>中国科学院高能物理研究所</div><div>中国 北京</div><div>18401696963</div><div><a href="mailto:yanyuechuixue@gmail.com" target="_blank">yanyuechuixue@gmail.com</a></div></div></div>
</span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>赵志超</div><div>中国科学院高能物理研究所</div><div>中国 北京</div><div>18401696963</div><div><a href="mailto:yanyuechuixue@gmail.com" target="_blank">yanyuechuixue@gmail.com</a></div></div></div>
</div>