Skip to content

A
AMSSNCKU
Commit 58db5f22 authored by Chen-Kai Qiao's avatar Chen-Kai Qiao
Browse files
Options

Replace AMSS_NCKU_Input.py

parent 15356f88
Pipeline #203 failed with stages
      Show whitespace changes
      Inline Side-by-side
      Showing with 200 additions and 213 deletions
      AMSS_NCKU_Input.py
      View file @ 58db5f22
      #################################################
      ##
      ## 这个文件包含了数值相对论程序 AMSS-NCKU 所需要的输入
      ## 这个文件包含了数值相对论所需要的输入
      ## 小曲
      ## 2024/03/19 --- 2024/12/02
      ##
      ......@@ -13,9 +13,9 @@ import numpy ## 导入 numpy 包
      ## 设置程序运行目录和计算资源
      File_directionary = "xiaoqu_Results_GW150914_test" ## 程序运行目录
      File_directionary = "xiaoqu_Results_GW150914_test1" ## 程序运行目录
      Output_directionary = "output_file" ## 存放二进制数据的子目录
      MPI_processes = 32 ## 想要调用的进程数目
      MPI_processes = 64 ## 想要调用的进程数目
      GPU_Calculation = "no" ## 是否开启 GPU 计算，可选 yes 或 no
      CPU_Part = 0.5
      ......@@ -31,7 +31,7 @@ GPU_Part = 0.5
      Symmetry = "equatorial-symmetry" ## 系统对称性，可选 equatorial-symmetry、no-symmetry
      Equation_Class = "BSSN" ## 设置方程形式，可选 BSSN、Z4C、BSSN-EScalar、BSSN-EM
      ## 注意：GPU 计算仅支持 BSSN
      Initial_Data_Method = "Ansorg-TwoPuncture" ## 设置初值的方法，可选 Ansorg-TwoPuncture、Lousto-Analytical、KerrSchild-Analytical
      Initial_Data_Method = "Ansorg-TwoPuncture" ## 设置求解数值相对论初值的方法，可选 Ansorg-TwoPuncture、Lousto-Analytical、KerrSchild-Analytical
      Time_Evolution_Method = "runge-kutta-45" ## 时间演化方法，可选 runge-kutta-45
      Finite_Diffenence_Method = "6th-order" ## 有限差分方法，可选 2nd-order、4th-order、6th-order、8th-order
      ......@@ -43,7 +43,7 @@ Finite_Diffenence_Method = "6th-order" ## 有限差分方法，可
      ## 设置时间演化信息
      Start_Evolution_Time = 0.0 ## 起始演化时间
      Final_Evolution_Time = 1500.0 ## 最终演化时间
      Final_Evolution_Time = 1000.0 ## 最终演化时间
      Check_Time = 1000.0
      Dump_Time = 50.0 ## 每隔一定时间间隔储存数据
      D2_Dump_Time = 300.0
      ......@@ -74,7 +74,7 @@ grid_minvalue = - grid_maxvalue ## 设置最外层格点的
      static_grid_number = 96 ## 设置固定格点每一层每一维数的格点数目
      moving_grid_number = 48 ## 设置可移动格点每一层每一维数的格点数目
      shell_grid_number = [32, 32, 80] ## 设置最外层球状网格（shell patch）的格点数目
      shell_grid_number = [32, 32, 100] ## 设置最外层球状网格（shell patch）的格点数目
      ## 以 phi、theta、r 的顺序给定
      devide_factor = 2.0 ## 设置相邻两层网格分辨率的比例（不要轻易改变）
      static_grid_type = 'Linear' ## 设置固定格点的类型，可选 'Linear'
      ......@@ -97,22 +97,25 @@ dimensionless_spin_BH = numpy.zeros( (puncture_number, 3) ) ## 初始化每个
      momentum_BH = numpy.zeros( (puncture_number, 3) ) ## 初始化每个黑洞的动量
      ## angular_momentum_BH = numpy.zeros( (puncture_number, 3) ) ## 初始化每个黑洞的自旋角动量
      puncture_data_set = "TwoPuncture" ## 设置 puncture data 的方式，可选 Manually 和 TwoPuncture
      puncture_data_set = "Manually" ## 设置双星轨道坐标的方式，可选 Manually 和 TwoPuncture
      #---------------------------------------------
      ## 如果设置 puncture data 的方式选为 TwoPuncture，只需要给定黑洞参数，偏心率，距离
      ## 这一步要与初值中的 Ansorg-TwoPuncture 配合使用，否则会报错
      ## 如果设置双星初始轨道坐标的方式选为 TwoPuncture，只需要给定黑洞参数，偏心率，距离即可
      ## 这一步要与初值求解中的 Ansorg-TwoPuncture 配合使用，否则可能会报错
      ## （也就是说，用 Ansorg-TwoPuncture 求解初值，轨道坐标设置可以设置 Manually 和 TwoPuncture，
      ## 但双星轨道坐标如果设置为 Manually 而不是 TwoPuncture，最好使用其它方法求解数值初值，否则可能会使 TwoPuncture 程序无法正确读入输入而报错）
      Distance = 10.0
      e0 = 0.0
      ## 设置每个黑洞的参数 (M Q a*)
      ## 质量 电荷 无量纲自旋
      ## 注意，如果设置 puncture data 的方式选为 TwoPuncture，第一个黑洞必须为质量较大的那个
      ## 设置每个黑洞的参数 (M Q* a*)
      ## 质量 无量纲电荷 无量纲自旋
      parameter_BH[0] = [ 36.0, 0.0, 0.31 ]
      parameter_BH[1] = [ 29.0, 0.0, -0.46 ]
      # parameter_BH[2] = [ 0.3247293, 0.0, 0.0 ] # 多黑洞手动补加
      # parameter_BH[2] = [ 0.3, 0.0, 0.0 ] # 多黑洞手动补加
      ## 注意，如果求解数值相对论初值的方法选为 Ansorg-TwoPuncture ，第一个黑洞必须为质量较大的那个，且则黑洞总质量会自动 rescale 为 M=1 （其它情况下必须手动 rescale）
      ## 设置每个黑洞的无量纲自旋
      ## 无对称性时 ，需要手动给 3 个方向的自旋角动量
      ......@@ -120,7 +123,8 @@ dimensionless_spin_BH[0] = [ 0.0, 0.0, 0.0 ]
      dimensionless_spin_BH[1] = [ 0.0, 0.0, 0.0 ]
      # dimensionless_spin_BH[2] = [ 0.0, 0.0, 0.0 ] # 多黑洞手动补加
      ## 注意，如果设置 puncture data 的方式选为 TwoPuncture，则自动调整将较大质量黑洞放在 y 轴正向，将较小质量黑洞放在 y 轴负向
      ## 注意，如果设置双星初始轨道坐标的方式选为 TwoPuncture，则程序自动调整将较大质量黑洞放在 y 轴正向，将较小质量黑洞放在 y 轴负向
      ## 如果设置双星初始轨道坐标的方式选为 Manually，且则需要手动调整到 y 轴方向
      ## use Brugmann's convention
      ## -----0-----> y
      ## - +
      ......@@ -130,14 +134,14 @@ dimensionless_spin_BH[1] = [ 0.0, 0.0, 0.0 ]
      ## 如果设置 puncture data 的方式选为 Manually，需要手动给定所有黑洞参数
      ## 设置每个黑洞的初始位置
      position_BH[0] = [ 0.0, +3.6321068, 0.0 ]
      position_BH[1] = [ 0.0, -3.6321068, 0.0 ]
      # position_BH[2] = [ 0.0, 0.0, 14.5284272 ] # 多黑洞手动补加
      position_BH[0] = [ 0.0, +4.4615385, 0.0 ]
      position_BH[1] = [ 0.0, -5.5384615, 0.0 ]
      # position_BH[2] = [ 0.0, 0.0, 0.0 ] # 多黑洞手动补加
      ## 设置每个黑洞的动量信息
      momentum_BH[0] = [ -0.0613136, 0.0, 0.0 ]
      momentum_BH[1] = [ +0.0613136, 0.0, 0.0 ]
      # momentum_BH[2] = [ 0.0, 0.0, -0.0668610 ] # 多黑洞手动补加
      momentum_BH[0] = [ -0.0953, -0.000845, 0.0 ]
      momentum_BH[1] = [ +0.0953, +0.000845, 0.0 ]
      # momentum_BH[2] = [ 0.0, 0.0, 0.0 ] # 多黑洞手动补加
      #################################################
      ......@@ -181,33 +185,16 @@ FR_phi0 = 0.00005
      ## boundary_choice = "BAM-choice" ## 索莫菲边界条件设定，可选 "BAM-choice" 和 "Shibata-choice"
      ## 目前的版本定死为 "BAM-choice"
      ## gauge_choice = # Gauge condition type
      ## gauge_choice = # 规范条件选取
      # 0: B^i gauge
      # 1: David's puncture gauge
      # 2: MB B^i gauge
      # 3: RIT B^i gauge
      # 4: MB beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
      # 5: RIT beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
      # 4: MB beta gauge
      # 5: RIT beta gauge
      # 6: MGB1 B^i gauge
      # 7: MGB2 B^i gauge
      ## 目前的版本定死为 0 或 1（根据方程不同选择）
      ## tetrad-type = 2 # tetradtype 0
      # v^a = (x,y,z)
      # orthonormal order: v,u,w
      # m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of PRD 75, 124018(2007)
      # tetradtype 1
      # orthonormal order: w,u,v
      # m = (theta + i phi)/sqrt(2) following Sperhake, Eq.(3.2) of PRD 85, 124062(2012)
      # tetradtype 2
      # v_a = (x,y,z)
      # orthonormal order: v,u,w
      # m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of PRD 75, 124018(2007)
      #################################################
      Markdown is supported
      0% or
      You are about to add 0 people to the discussion. Proceed with caution.
      Finish editing this message first!
      Please register or to comment