Sorry lan Hinder,
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).
Now I know your meaning, and I am trying to solve them.
Thank you.
Yours,
Zhi-Chao.
2017-08-18 22:55 GMT+08:00 赵志超 yanyuechuixue@gmail.com:
Thank you Ian Hinder.
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.
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.
Yours, Zhi-Chao.
2017-08-18 22:43 GMT+08:00 Ian Hinder ian.hinder@aei.mpg.de:
On 18 Aug 2017, at 10:05, 赵志超 yanyuechuixue@gmail.com wrote:
Thank you Ian Hinder.
Yes I want to add a low density fluid to a BBH simulation.
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".
so I think it can give a low density fluid to a BBH simulation.
But it report "NAN" or "Con2Prim: count > GRHydro_countmax".
OK, thanks for the clarification!
I'm not an expert on GRHydro, but maybe someone else can suggest something?
Are you planning to solve the constraints?
-- Ian Hinder http://members.aei.mpg.de/ianhin
-- 赵志超 中国科学院高能物理研究所 中国 北京 18401696963 yanyuechuixue@gmail.com