[Users] Excision with McLachlan
Jonah Miller
jonah.maxwell.miller at gmail.com
Mon Mar 27 09:59:51 CDT 2017
Hi Miguel,
I'm not sure about excising dynamically during evolution--maybe someone
else has insight there. However, excision is definitely possible. I have
been using the sixpatches Llama thorn to perform excision for a single
black hole. I attach a parameter file to excise a Kerr-Schild black hole
that I based my simulations on.
Best,
Jonah
On 2017-03-27 05:01 AM, Miguel Zilhão wrote:
> hi all,
>
> i was wondering whether one could straightforwardly use the ET to perform black hole
> evolutions using excision instead of punctures.
>
> specifically, i had in mind a collapsing configuration where i'd like to excise the inner
> domain once a black hole horizon is formed... is this possible?
> i had a look around and noticed the LegoExcision thorn, but i'm not sure whether it's
> actively maintained, nor whether it can interface with McLachlan...
>
> thanks,
> Miguel
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users
-------------- next part --------------
ActiveThorns = "
ADMBase
ADMCoupling
ADMMacros
AEILocalInterp
AHFinderDirect
Boundary
Carpet
CarpetIOASCII
CarpetIOBasic
CarpetIOHDF5
CarpetIOScalar
CarpetInterp
CarpetInterp2
CarpetLib
CarpetReduce
CartGrid3D
CoordBase
CoordGauge
Coordinates
Formaline
GenericFD
IOUtil
InitBase
Interpolate2
KerrSchild
ML_BSSN
ML_BSSN_Helper
ML_LocalInterp
MoL
NaNChecker
NewRad
QuasiLocalMeasures
SpaceMask
SphericalSurface
StaticConformal
SymBase
SystemStatistics
Time
TimerReport
TmunuBase
"
#############################################################
# Grid
#############################################################
CartGrid3D::type = "multipatch"
Coordinates::coordinate_system = "Thornburg04nc"
Coordinates::h_radial = 1.0/(5*@RHO@) # (order+1)*n cells
Coordinates::sphere_inner_radius = 1.0
Coordinates::sphere_outer_radius = 10.0 # 100
Coordinates::n_angular = 5*@RHO at -1 # (order+1)*n-1
Coordinates::stagger_patch_boundaries = yes
Coordinates::stagger_outer_boundaries = yes
Interpolate2::interpolator_order = 4
Interpolate2::shift_edges = no
Coordinates::store_jacobian = yes # I need to
Coordinates::store_inverse_jacobian = no # Check for runtime errors
Coordinates::store_jacobian_derivative = no
Coordinates::store_volume_form = no
GenericFD::jacobian_group = "Coordinates::jacobian"
GenericFD::jacobian_derivative_group = "Coordinates::jacobian2"
# OVERLAP
Coordinates::additional_overlap_size = 3
# BOUNDARY ZONES (needs to be set to three) (also in coordinates thorne)
Coordinates::patch_boundary_size = 3
Coordinates::outer_boundary_size = 3
#############################################################
# Carpet
#############################################################
Carpet::domain_from_multipatch = yes
Carpet::ghost_size = 3
# not for FD
#Carpet::granularity = 5 # (order+1)*m
#Carpet::granularity_boundary = 1
Carpet::poison_new_timelevels = no # does not work with persistent RHS
#############################################################
# Time integration
#############################################################
#Cactus::terminate = "time"
#Cactus::cctk_final_time = @FINAL@
Cactus::cctk_itlast = 1
Time::timestep_method = "given"
#Time::timestep = 1E-3
Time::timestep = 0.05/@RHO@
#MethodOfLines::ODE_method = "Euler"
MethodOfLines::ODE_method = "RK4"
#MethodOfLines::MoL_Intermediate_Steps = 1
MethodOfLines::MoL_Intermediate_Steps = 4
MethodOfLines::MoL_Num_Scratch_Levels = 1
MethodOfLines::init_RHS_zero = yes
MethodOfLines::MoL_NaN_Check = no
#MethodOfLines::MoL_NaN_Check = yes
#############################################################
# ML_BSSN
#############################################################
ADMBase::initial_data = "KerrSchild"
ADMBase::initial_lapse = "KerrSchild"
ADMBase::initial_shift = "KerrSchild"
ADMBase::initial_dtlapse = "KerrSchild"
ADMBase::initial_dtshift = "KerrSchild"
KerrSchild::M = 1.0
KerrSchild::a = 0.0
KerrSchild::exact_method = "KerrSchild"
ADMBase::evolution_method = "ML_BSSN"
ADMBase::lapse_evolution_method = "ML_BSSN"
ADMBase::shift_evolution_method = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"
ML_BSSN::timelevels = 3
ML_BSSN::other_timelevels = 3
ML_BSSN::rhs_timelevels = 1
#ML_BSSN::rhs_scheduling = "persistent"
ML_BSSN::conformalMethod = 1 # W
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::harmonicF = 2.0 # 1+log
ML_BSSN::shiftFormulation = 0 # Gamma driver
ML_BSSN::shiftGammaCoeff = 0.75
ML_BSSN::alphaDriver = 1.0
ML_BSSN::betaDriver = 1.0
ML_BSSN::evolveA = 1
ML_BSSN::evolveB = 1
ML_BSSN::advectLapse = 0
ML_BSSN::advectShift = 0
ML_BSSN::fixAdvectionTerms = 1
ML_BSSN::minimumLapse = 1.0e-10
ML_BSSN::confromalMethod = 1 # W
ML_BSSN::fixAdcectionTerms = 1
ML_BSSN::initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::rhs_boundary_condition = "NewRad"
Boundary::radpower = 2
#############################################################
# Analysis
#############################################################
AHFinderDirect::find_every = 20*@RHO@
AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation"
AHFinderDirect::N_horizons = 1
AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0
AHFinderDirect::reset_horizon_after_not_finding[1] = no
AHFinderDirect::verbose_level = "physics details"
#AHFinderDirect::verbose_level = "algorithm highlights"
AHFinderDirect::output_h_every = 1
AHFinderDirect::output_OpenDX_control_files = false
#############################################################
# NaNChecker
#############################################################
# Debugging stuff. Comment out for production runs.
NaNChecker::check_every = 1
NaNChecker::check_after = 0
NaNChecker::check_vars = "
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
"
NaNChecker::action_if_found = "just warn"
NaNChecker::verbose = "all"
#############################################################
# Output
#############################################################
IO::out_dir = "$parfile"
CarpetIOBasic::outInfo_every = 20*@RHO@
CarpetIOBasic::outInfo_reductions = "sum norm2"
CarpetIOBasic::outInfo_vars = "
Carpet::physical_time_per_hour
SystemStatistics::maxrss_mb
SystemStatistics::swap_used_mb
ML_BSSN::alpha
ML_BSSN::H
"
IOScalar::one_file_per_group = yes
IOScalar::out_precision = 17
IOScalar::outScalar_every = 20*@RHO@
IOScalar::outScalar_vars = "
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOASCII::one_file_per_group = yes
IOASCII::out_precision = 17
IOASCII::out1D_every = 20*@RHO@
IOASCII::out1D_xline_z = 3.5
IOASCII::out1D_yline_z = 3.5
IOASCII::out1D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOHDF5::one_file_per_group = yes
IOHDF5::out1D_every = 0 #1*20*@RHO@
IOHDF5::out1D_xline_z = 3.5
IOHDF5::out1D_yline_z = 3.5
IOHDF5::out1D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOHDF5::out2D_every = 0*20*@RHO@
IOHDF5::out2D_xyplane_z = 3.5
IOHDF5::out2D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOASCII::out2D_every = 0*20*@RHO@
IOASCII::out2D_xyplane_z = 3.5
IOASCII::out2D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOHDF5::out3D_every = 0*20*@RHO@
IOHDF5::out3D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
IOASCII::out3D_every = 0 #0*20*@RHO@
IOASCII::out3D_vars = "
grid::coordinates
CarpetReduce::weight
Coordinates::jacobian
Interpolate::source_patch
ML_BSSN::ML_confac
ML_BSSN::ML_metric
ML_BSSN::ML_Gamma
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
ML_BSSN::ML_confacrhs
ML_BSSN::ML_metricrhs
ML_BSSN::ML_Gammarhs
ML_BSSN::ML_trace_curvrhs
ML_BSSN::ML_curvrhs
ML_BSSN::ML_lapserhs
ML_BSSN::ML_shiftrhs
ML_BSSN::ML_dtlapserhs
ML_BSSN::ML_dtshiftrhs
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
ML_BSSN::ML_cons_detg
ML_BSSN::ML_cons_Gamma
ML_BSSN::ML_cons_traceA
ML_BSSN::At11rhs
ML_BSSN::At12rhs
ML_BSSN::At13rhs
ML_BSSN::At22rhs
ML_BSSN::At23rhs
ML_BSSN::At33rhs
ML_BSSN::Arhs
ML_BSSN::B1rhs
ML_BSSN::B2rhs
ML_BSSN::B3rhs
ML_BSSN::Xt1rhs
ML_BSSN::Xt2rhs
ML_BSSN::Xt3rhs
ML_BSSN::gt11rhs
ML_BSSN::gt12rhs
ML_BSSN::gt13rhs
ML_BSSN::gt22rhs
ML_BSSN::gt23rhs
ML_BSSN::gt33rhs
ML_BSSN::trKrhs
ADMBase::metric
ADMBase::curv
ADMBase::lapse
ADMBase::shift
ADMBase::dtlapse
ADMBase::dtshift
KerrSchild::metric_exact
KerrSchild::curv_exact
KerrSchild::lapse_exact
KerrSchild::shift_exact
KerrSchild::dtlapse_exact
KerrSchild::dtshift_exact
"
# The timer output
Cactus::cctk_timer_output = "full"
TimerReport::out_every = 10
TimerReport::out_filename = "timer.asc"
TimerReport::output_all_timers_readable = yes
Carpet::output_timers_every = 10
CarpetLib::print_timestats_every = 10
More information about the Users
mailing list