New pseudospectral code for the construction of initial data
New pseudospectral code for the construction of initial data
复制标题
DOI：
10.1103/physrevd.105.104027
复制
发表时间：
2021-09
影响因子：
5
通讯作者：
A. Rashti;F. Fabbri;B. Brügmann;Swami Vivekanandji Chaurasia;T. Dietrich;M. Ujevic;W. Tichy
中科院分区：
文献类型：
--
作者：
A. Rashti;F. Fabbri;B. Brügmann;Swami Vivekanandji Chaurasia;T. Dietrich;M. Ujevic;W. Tichy
引力系统动力学的数值研究，例如黑洞-中子星系统，需要满足物理和约束的初始数据。在本文中，我们提出了新开发的伪谱代码Elliptica，它是构造各种双引力和单引力系统初始数据的基础。所考虑的椭圆方程在时空流形的单个空间超曲面上求解。利用坐标映射，超表面被补丁覆盖，这些补丁的边界可以适应紧凑物体的表面。为了求解具有任意边界条件的椭圆型方程，Elliptica采用了直接求解器的Schur补域分解方法。在这个版本中，我们使用立方球坐标映射，并使用第一类切比雪夫多项式展开域。在这里，我们解释了椭圆的构建模块和黑洞-中子星双星系统的初始数据构建算法。我们执行收敛测试并改进数据以验证我们的结果。在我们的框架内，中子星可以在任意方向达到接近分裂的自旋值，而黑洞可以具有任意自旋，无量纲自旋量级为0。8.
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica , an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, Elliptica deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of Elliptica and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude ∼ 0 . 8.