Full implementation of AGN feedback for the cool-core destruction project

doc/wiki/Runtime-Parameters-related/Runtime-Parameters:-All.md

@@ -398,6 +398,10 @@ For variables with `Default/Min/Max` labeled as `Depend`, click the parameter na
 | SPEC_INT_VORTEX_THRESHOLD                                                                            |             0.1 |             0.0 |            None | vortex detection threshold for SPEC_INT_XY_INSTEAD_DEPHA [0.1] ##ELBDM & SUPPORT_SPECTRAL_INT ONLY## |
 | SPEC_INT_XY_INSTEAD_DEPHA                                                                            |               1 |            None |            None | interpolate x and y (real and imaginary parts in current implementation) around vortices instead of density and phase for the spectral interpolation, which has the advantage of being well-defined across vortices [1] ##ELBDM & SUPPORT_SPECTRAL_INT ONLY## |
 | SRC_DELEPTONIZATION                                                                                  |               0 |            None |            None | deleptonization (for simulations of stellar core collapse) [0] ##HYDRO ONLY## |
+| SRC_EC_DTCOEF                                                                                        |             0.5 | 2.22507386e-308 |            None | the coefficient of the cooling time step [0.5] |
+| SRC_EC_SUBCYCLING                                                                                    |               0 |            None |            None | perform subcycling when the cooling time step is small (0=off, 1=on) [0] |
+| SRC_EC_TEF_N                                                                                         |            1501 |               1 |            None | number of points for lambda(T) sampling in LOG [1501] |
+| SRC_EXACTCOOLING                                                                                     |               0 |            None |            None | exact cooling scheme from Gaspari (2009) [0] ##HYDRO ONLY## |
 | SRC_GPU_NPGROUP                                                                                      |              -1 |            None |            None | number of patch groups sent into the CPU/GPU source-term solver (<=0=auto) [-1] |
 | SRC_USER                                                                                             |               0 |            None |            None | user-defined source terms -> edit "Src_User.cpp" [0] |
 

example/test_problem/Hydro/ClusterMerger/Input_Options/case1/Input__Flag_Lohner

@@ -0,0 +1,13 @@
+# Level    Threshold_Refine  Threshold_Derefine              Filter              Soften          MinDensity
+    0                  0.80                0.80                0.01                0.00             1.50e-2
+    1                  0.80                0.80                0.01                0.00             1.50e-2
+    2                  0.80                0.80                0.01                0.00             1.50e-2
+    3                  0.80                0.80                0.01                0.00             1.50e-2
+    4                  8000                8000                10                  0.00             1.50e-2
+    5                  8000                8000                10                  0.00             1.50e-2
+    6                  8000                8000                10                  0.00             1.50e-2
+    7                  8000                8000                10                  0.00             1.50e-2
+    8                  8000                8000                10                  0.00             1.50e-2
+    9                  8000                8000                10                  0.00             1.50e-2
+   10                  8000                8000                10                  0.00             1.50e-2
+   11                  8000                8000                10                  0.00             1.50e-2

example/test_problem/Hydro/ClusterMerger/Input_Options/case1/Input__Flag_NParPatch

@@ -0,0 +1,13 @@
+# Level   Number of Particles per Patch
+      0                             600
+      1                             600
+      2                             600
+      3                             600
+      4                         6000000
+      5                         6000000
+      6                         6000000
+      7                         6000000
+      8                         6000000
+      9                         6000000
+     10                         6000000
+     11                         6000000

example/test_problem/Hydro/ClusterMerger/Input_Options/case1/Input__TestProb

@@ -0,0 +1,61 @@
+# problem-specific runtime parameters
+Merger_Coll_NumHalos    2                        # number of halos
+
+# parameters for cluster 1
+Merger_File_Prof1       profile_CCD_C1.h5        # profile table of cluster 1
+Merger_File_Par1        1to3_b0.0_CCD_C1.h5      # particle file of cluster 1
+Merger_Coll_PosX1       5740.67470419            # X-center   of cluster 1 in kpc
+Merger_Coll_PosY1       7500.0                   # Y-center   of cluster 1 in kpc
+Merger_Coll_VelX1       600.0                    # X-velocity of cluster 1 in km/s
+Merger_Coll_VelY1       0.0                      # Y-velocity of cluster 1 in km/s
+Merger_Coll_IsGas1      1                        # If cluster 1 should have gas and not have DM only
+Bondi_MassBH1           7.11601760524e10         # black hole mass (in Msun)
+Jet_HalfHeight1         2.0                      # jet 1: half height of the cylinder-shape jet source (in kpc)
+Jet_Radius1             2.5                      # jet 1: radius of the cylinder-shape jet source (in kpc)
+
+# parameters for cluster 2
+Merger_File_Prof2       profile_CCD_C1.h5        # profile table of cluster 2
+Merger_File_Par2        1to3_b0.0_CCD_C1.h5      # particle file of cluster 2
+Merger_Coll_PosX2       9259.32529581            # X-center   of cluster 2 in kpc
+Merger_Coll_PosY2       7500.0                   # Y-center   of cluster 2 in kpc
+Merger_Coll_VelX2       -600.0                   # X-velocity of cluster 2 in km/s
+Merger_Coll_VelY2       0.0                      # Y-velocity of cluster 2 in km/s
+Merger_Coll_IsGas2      1                        # If cluster 2 should have gas and not have DM only
+Bondi_MassBH2           7.11601760524e10         # black hole mass (in Msun)
+Jet_HalfHeight2         2.0                      # jet 2: half height of the cylinder-shape jet source (in kpc)
+Jet_Radius2             2.5                      # jet 2: radius of the cylinder-shape jet source (in kpc)
+
+# parameters for cluster 3 (disabled for this example)
+Merger_File_Prof3       none                     # profile table of cluster 3
+Merger_File_Par3        none                     # particle file of cluster 3
+Merger_Coll_PosX3       0                        # X-center   of cluster 3 in kpc
+Merger_Coll_PosY3       0.0                      # Y-center   of cluster 3 in kpc
+Merger_Coll_VelX3       0.0                      # X-velocity of cluster 3 in km/s
+Merger_Coll_VelY3       0.0                      # Y-velocity of cluster 3 in km/s
+Merger_Coll_IsGas3      0                        # If cluster 3 should have gas and not have DM only
+Bondi_MassBH3           3.4e8                    # black hole mass (in Msun)
+Jet_HalfHeight3         0.70                     # jet 3: half height of the cylinder-shape jet source (in kpc)
+Jet_Radius3             0.35                     # jet 3: radius of the cylinder-shape jet source (in kpc)
+
+# parameters of AGN jet feedback
+AGN_feedback            1                        # turn on/off (1/0) AGN feedback
+Accretion_Mode          3                        # 1: hot mode; 2: cold mode; 3: combine (hot + cold)
+eta                     1.0                      # mass loading factor in jet feedback
+eps_f                   0.01                     # the radiative efficiency in jet feedback
+eps_m                   0.0                      # the fraction of total energy that goes into the thermal
+                                                 # energy in jet feedback
+R_acc                   4.0                      # accretion radius: compute the accretion rate (in kpc)
+R_dep                   4.0                      # radius to deplete the accreted gas (in kpc)
+JetDirection_case       3                        # Methods for choosing the jet direction:
+                                                 #    1: Fixed at x-axis
+                                                 #    2: Import from table (generate JetDirection_file)
+                                                 #    3: Align with angular momentum
+JetDirection_file       JetDirection.txt         # File containing the jet direction as a function of time
+
+# other parameters
+Merger_Coll_UseMetals   0                        # look for a metal field in the profile files [1]
+Merger_Coll_LabelCenter 1                        # label the particle closest to the center of each cluster [1]
+AdjustBHPos             1                        # (true/false) --> Adjust the BH position
+AdjustBHVel             1                        # (true/false) --> Adjust the BH velocity
+AdjustPeriod            0.001                    # the time interval of adjustment (in Myr)
+fixBH                   0                        # fix the BH at the simulation box center and set its velocity to be zero (1 cluster only)

example/test_problem/Hydro/ClusterMerger/Input_Options/case2/Input__Flag_Lohner

@@ -0,0 +1,13 @@
+# Level    Threshold_Refine  Threshold_Derefine              Filter              Soften          MinDensity
+    0                  0.80                0.80                0.01                0.00             1.50e-2
+    1                  0.80                0.80                0.01                0.00             1.50e-2
+    2                  0.80                0.80                0.01                0.00             1.50e-2
+    3                  0.80                0.80                0.01                0.00             1.50e-2
+    4                  8000                8000                10                  0.00             1.50e-2
+    5                  8000                8000                10                  0.00             1.50e-2
+    6                  8000                8000                10                  0.00             1.50e-2
+    7                  8000                8000                10                  0.00             1.50e-2
+    8                  8000                8000                10                  0.00             1.50e-2
+    9                  8000                8000                10                  0.00             1.50e-2
+   10                  8000                8000                10                  0.00             1.50e-2
+   11                  8000                8000                10                  0.00             1.50e-2

example/test_problem/Hydro/ClusterMerger/Input_Options/case2/Input__Flag_NParPatch

@@ -0,0 +1,13 @@
+# Level   Number of Particles per Patch
+      0                             600
+      1                             600
+      2                             600
+      3                             600
+      4                         6000000
+      5                         6000000
+      6                         6000000
+      7                         6000000
+      8                         6000000
+      9                         6000000
+     10                         6000000
+     11                         6000000