diff --git a/Examples/Fluid_Kerr/params.txt b/Examples/Fluid_Kerr/params.txt index de6892a03..6ac48473a 100644 --- a/Examples/Fluid_Kerr/params.txt +++ b/Examples/Fluid_Kerr/params.txt @@ -8,21 +8,21 @@ verbosity = 0 # location / naming of output files # output_path = "" # Main path for all files. Must exist! -chk_prefix = PerfectFluid_ -plot_prefix = PerfectFluidp_ -# restart_file = PerfectFluid_000000.3d.hdf5 +chk_prefix = FluidKerrBH_ +plot_prefix = FluidKerrBHp_ +#restart_file = # HDF5files are written every dt = L/N*dt_multiplier*checkpoint_interval -checkpoint_interval = 500 +checkpoint_interval = 1 # set to 0 to turn off plot files (except at t=0 and t=stop_time) # set to -1 to never ever print plotfiles -plot_interval = 1 -num_plot_vars = 2 -plot_vars = rho eps +plot_interval = 0 +num_plot_vars = 0 +# plot_vars = # subpaths - specific directories for hdf5, pout, extraction data # (these are created at runtime) -hdf5_subpath = "hdf5" +hdf5_subpath = "/sakura/ptmp/dtraykova/fluid_test_Kerr_4lvl/hdf5" pout_subpath = "pout" data_subpath = "data" @@ -41,16 +41,17 @@ print_progress_only_to_rank_0 = 1 # Here we decouple the evolution so the scalar evolved on the # metric background without backreaction (this avoids the need # to solve the constaints) -G_Newton = 0.0 - -# Scalar field initial data -#scalar_amplitude = 0.1 -#scalar_width = 5.0 -#scalar_mass = 0.2 +G_Newton = 1.0 +# Fluid params +fluid_rho0 = 1.0 +fluid_width = 5. +fluid_delta = 0.2 + # Kerr BH data -kerr_mass = 0. -kerr_spin = 0.0 +kerr_mass = 1.0 +kerr_spin = 0.5 +# kerr_center = 64 64 64 # defauts to center ################################################# # Grid parameters @@ -62,28 +63,25 @@ kerr_spin = 0.0 # NB - if you have a non-cubic grid, you can specify 'N1' or 'N1_full', # 'N2' or 'N2_full' and 'N3' or 'N3_full' ( then dx_coarsest = L/N(max) ) # NB - the N values need to be multiples of the block_factor -N_full = 64 #128 -L_full = 128 #256 +N_full = 32 +L_full = 128 # Maximum number of times you can regrid above coarsest level -max_level = 5 # There are (max_level+1) grids, so min is zero +max_level = 4 # There are (max_level+1) grids, so min is zero # Frequency of regridding at each level and thresholds on the tagging # Need one for each level except the top one, ie max_level items # Generally you do not need to regrid frequently on every level -# in this example turn off regridding on all levels -# Level Regridding: 0 1 2 3 4 5 -regrid_interval = 0 0 0 0 0 0 -# regrid_threshold = 0.5 - -regrid_interval = 0 0 0 0 0 0 +# Level Regridding: 0 1 2 3 4 +#regrid_interval = 50 25 5 5 5 +#regrid_thresholds = 0.3 0.3 0.3 0.3 0.3 +regrid_interval = 0 0 0 0 # Max and min box sizes max_box_size = 16 min_box_size = 16 -tag_buffer_size = 0 # this example uses a fixed grid - +# tag_buffer_size = 3 # grid_buffer_size = 8 # fill_ratio = 0.7 # num_ghosts = 3 @@ -92,14 +90,13 @@ tag_buffer_size = 0 # this example uses a fixed grid ################################################# # Boundary Conditions parameters -#Periodic directions - 0 = false, 1 = true -isPeriodic = 1 1 1 +# Periodic directions - 0 = false, 1 = true +isPeriodic = 0 0 0 # if not periodic, then specify the boundary type # 0 = static, 1 = sommerfeld, 2 = reflective -# 3 = extrapolating, 4 = mixed # (see BoundaryConditions.hpp for details) -hi_boundary = 4 4 4 -lo_boundary = 2 2 2 +hi_boundary = 3 3 3 +lo_boundary = 3 3 2 # if reflective boundaries selected, must set # parity of all vars (in order given by UserVariables.hpp) @@ -111,8 +108,8 @@ vars_parity = 0 0 4 6 0 5 0 #chi and hij 0 0 4 6 0 5 0 #K and Aij 0 1 2 3 #Theta and Gamma 0 1 2 3 1 2 3 #lapse shift and B - 0 1 2 3 0 #D Sj and tau - 1 2 3 0 0 #vi rho and epsilon + 0 1 2 3 0 #D Sj and tau + 1 2 3 0 0 #vi rho and epsilon vars_parity_diagnostic = 0 1 2 3 #Ham and Mom # if sommerfeld boundaries selected, must select @@ -122,17 +119,17 @@ nonzero_asymptotic_vars = chi h11 h22 h33 lapse nonzero_asymptotic_values = 1.0 1.0 1.0 1.0 1.0 # if you are using extrapolating BC: -extrapolation_order = 0 #1 -num_extrapolating_vars = 2 -extrapolating_vars = phi Pi +extrapolation_order = 0 +num_extrapolating_vars = 10 +extrapolating_vars = D Sj1 Sj2 Sj3 tau vi1 vi2 vi3 rho eps ################################################# # Evolution parameters # dt will be dx*dt_multiplier on each grid level -dt_multiplier = 0.078125 #0.25 -stop_time = 5.0 -# max_steps = 4 +dt_multiplier = 0.1 +stop_time = 100.0 +#max_steps = 4 # Spatial derivative order (only affects CCZ4 RHS) max_spatial_derivative_order = 4 # can be 4 or 6 diff --git a/Examples/PerfectFluid_Kerr/test_Gauss_small/params.txt b/Examples/PerfectFluid_Kerr/test_Gauss_small/params.txt new file mode 100644 index 000000000..6627cb810 --- /dev/null +++ b/Examples/PerfectFluid_Kerr/test_Gauss_small/params.txt @@ -0,0 +1,162 @@ +# See the wiki page for an explanation of the params! +# https://github.com/GRChombo/GRChombo/wiki/Guide-to-parameters + +################################################# +# Filesystem parameters + +verbosity = 0 + +# location / naming of output files +# output_path = "" # Main path for all files. Must exist! +chk_prefix = FluidKerrBH_ +plot_prefix = FluidKerrBHp_ +#restart_file = /sakura/ptmp/dtraykova/Kerr/hdf5/KerrBH_000060.3d.hdf5 + +# HDF5files are written every dt = L/N*dt_multiplier*checkpoint_interval +checkpoint_interval = 1 +# set to 0 to turn off plot files (except at t=0 and t=stop_time) +# set to -1 to never ever print plotfiles +plot_interval = 0 +num_plot_vars = 0 +# plot_vars = + +# subpaths - specific directories for hdf5, pout, extraction data +# (these are created at runtime) +hdf5_subpath = "/sakura/ptmp/dtraykova/fluid_test_Kerr_4lvl/hdf5" +pout_subpath = "pout" +data_subpath = "data" + +# change the name of output files +# pout_prefix = "pout" +print_progress_only_to_rank_0 = 1 + +# ignore_checkpoint_name_mismatch = 0 +# write_plot_ghosts = 0 + +################################################# +# Initial Data parameters + +# Change the gravitational constant of the Universe! +# Default is 1.0, for standard geometric units +# Here we decouple the evolution so the scalar evolved on the +# metric background without backreaction (this avoids the need +# to solve the constaints) +G_Newton = 1.0 + +# Fluid params +fluid_rho0 = 1.0 +fluid_width = 5. +fluid_delta = 0.2 + +# Kerr BH data +kerr_mass = 1.0 +kerr_spin = 0.5 +# kerr_center = 64 64 64 # defauts to center + +################################################# +# Grid parameters + +# 'N' is the number of subdivisions in each direction of a cubic box +# 'L' is the length of the longest side of the box, dx_coarsest = L/N +# NB - If you use reflective BC and want to specify the subdivisions and side +# of the box were there are no symmetries, specify 'N_full' and 'L_full' instead +# NB - if you have a non-cubic grid, you can specify 'N1' or 'N1_full', +# 'N2' or 'N2_full' and 'N3' or 'N3_full' ( then dx_coarsest = L/N(max) ) +# NB - the N values need to be multiples of the block_factor +N_full = 32 +L_full = 128 + +# Maximum number of times you can regrid above coarsest level +max_level = 4 # There are (max_level+1) grids, so min is zero + +# Frequency of regridding at each level and thresholds on the tagging +# Need one for each level except the top one, ie max_level items +# Generally you do not need to regrid frequently on every level +# Level Regridding: 0 1 2 3 4 +#regrid_interval = 50 25 5 5 5 +#regrid_thresholds = 0.3 0.3 0.3 0.3 0.3 +regrid_interval = 0 0 0 0 + +# Max and min box sizes +max_box_size = 16 +min_box_size = 16 + +# tag_buffer_size = 3 +# grid_buffer_size = 8 +# fill_ratio = 0.7 +# num_ghosts = 3 +# center = 256.0 256.0 256.0 # defaults to center of the grid + +################################################# +# Boundary Conditions parameters + +# Periodic directions - 0 = false, 1 = true +isPeriodic = 0 0 0 +# if not periodic, then specify the boundary type +# 0 = static, 1 = sommerfeld, 2 = reflective +# (see BoundaryConditions.hpp for details) +hi_boundary = 3 3 3 +lo_boundary = 3 3 2 + +# if reflective boundaries selected, must set +# parity of all vars (in order given by UserVariables.hpp) +# 0 = even +# 1,2,3 = odd x, y, z +# 4,5,6 = odd xy, yz, xz +# 7 = odd xyz +vars_parity = 0 0 4 6 0 5 0 #chi and hij + 0 0 4 6 0 5 0 #K and Aij + 0 1 2 3 #Theta and Gamma + 0 1 2 3 1 2 3 #lapse shift and B + 0 1 2 3 0 #D Sj and tau + 1 2 3 0 0 #vi rho and epsilon +vars_parity_diagnostic = 0 1 2 3 #Ham and Mom + +# if sommerfeld boundaries selected, must select +# non zero asymptotic values +num_nonzero_asymptotic_vars = 5 +nonzero_asymptotic_vars = chi h11 h22 h33 lapse +nonzero_asymptotic_values = 1.0 1.0 1.0 1.0 1.0 + +# if you are using extrapolating BC: +extrapolation_order = 0 +num_extrapolating_vars = 10 +extrapolating_vars = D Sj1 Sj2 Sj3 tau vi1 vi2 vi3 rho eps + +################################################# +# Evolution parameters + +# dt will be dx*dt_multiplier on each grid level +dt_multiplier = 0.1 +stop_time = 100.0 +#max_steps = 4 + +# Spatial derivative order (only affects CCZ4 RHS) +max_spatial_derivative_order = 4 # can be 4 or 6 + +nan_check = 1 + +# Lapse evolution +lapse_advec_coeff = 1.0 +lapse_coeff = 2.0 +lapse_power = 1.0 + +# Shift evolution +shift_advec_coeff = 0.0 # Usually no advection for beta +shift_Gamma_coeff = 0.75 +eta = 1.0 # eta of gamma driver, should be of order ~1/M_ADM of spacetime + +# CCZ4 parameters +formulation = 1 # 1 for BSSN, 0 for CCZ4 +kappa1 = 0. +kappa2 = 0. +kappa3 = 0. +covariantZ4 = 1 # 0: keep kappa1; 1 [default]: replace kappa1 -> kappa1/lapse + +# coefficient for KO numerical dissipation +sigma = 0.3 + +# min_chi = 1.e-4 +# min_lapse = 1.e-4 + +#################################################