From 4fbfc1c5921bd472fb989e87eaf05dc99a2d109a Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sat, 14 Nov 2009 22:59:25 -0600 Subject: Add description for converting BSSN_MoL parameter files to McLachlan --- doc/PARAMETERS | 100 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 100 insertions(+) create mode 100644 doc/PARAMETERS (limited to 'doc') diff --git a/doc/PARAMETERS b/doc/PARAMETERS new file mode 100644 index 0000000..14e4369 --- /dev/null +++ b/doc/PARAMETERS @@ -0,0 +1,100 @@ +Comparison between BSSN_MoL (ADM_BSSN) and McLachlan (ML_BSSN) parameters + +2009-11-14: Written by Erik Schnetter, not yet trustworthy +XXXX-XX-XX: Cross-checked by NNN + + + +Formulation: + +ADMBase::evolution_method = "ADM_BSSN" + ADMBase::evolution_method = "ML_BSSN" + +ADM_BSSN::RfromH = yes + [always for ML_BSSN] + +ADM_BSSN::bound = "radiative" + ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" + ML_BSSN::my_rhs_boundary_condition = "NewRad" + # These are necessary + ML_BSSN::ML_log_confac_bound = "none" + ML_BSSN::ML_metric_bound = "none" + ML_BSSN::ML_Gamma_bound = "none" + ML_BSSN::ML_trace_curv_bound = "none" + ML_BSSN::ML_curv_bound = "none" + ML_BSSN::ML_lapse_bound = "none" + ML_BSSN::ML_dtlapse_bound = "none" + ML_BSSN::ML_shift_bound = "none" + ML_BSSN::ML_dtshift_bound = "none" + +Boundary::radpower + Boundary::radpower + + + +Lapse: + +ADMBase::lapse_evolution_method = "1+log" + ADMBase::lapse_evolution_method = "ML_BSSN" + ADMBase::dtlapse_evolution_method = "ML_BSSN" + ML_BSSN::harmonicN = 1 # 1+log + +ADM_BSSN::lapsesource = "straight" + ML_BSSN::AlphaDriver = 0 +ADM_BSSN::lapsesource = "modified" + ML_BSSN::AlphaDriver = 1 + +ADM_BSSN::harmonic_f + ML_BSSN::harmonicF + +ADM_BSSN::force_lapse_positive = yes + [always for ML_BSSN, but MinimumLapse can be negative] +ADM_BSSN::force_lapse_minimum + ML_BSSN::MinimumLapse + +ADM_BSSN::lapse_advection_coeff + ML_BSSN::LapseAdvectionCoeff + [should be 0 or 1] + + + +Shift: + +ADMBase::shift_evolution_method = "gamma0" + ADMBase::shift_evolution_method = "ML_BSSN" + ADMBase::dtshift_evolution_method = "ML_BSSN" + ML_BSSN::ShiftGammaCoeff = 0 +ADMBase::shift_evolution_method = "gamma2" + ADMBase::shift_evolution_method = "ML_BSSN" + ADMBase::dtshift_evolution_method = "ML_BSSN" + ML_BSSN::ShiftGammaCoeff = 1 + +ADM_BSSN::ShiftGammaCoeff + ML_BSSN::ShiftGammaCoeff + +ADM_BSSN::ShiftAlpPower = 0 + [always for ML_BSSN] + [There is a parameter ML_BSSN::ShiftAlphaPower, but it is unused] + +ADM_BSSN::GammaDriver + [has no effect in ADM_BSSN] + +ADM_BSSN::BetaDriver + ML_BSSN::BetaDriver + +ADM_BSSN::gamma_driver_advection_coeff + ML_BSSN::ShiftAdvectionCoeff + [should be 0 or 1] + +ADM_BSSN::ApplyShiftBoundary = no + [not possible with ML_BSSN; shift boundary is always applied] + + + +Discretisation: + +ADM_BSSN::timelevels + ML_BSSN::timelevels + +ADM_BSSN::advection = "upwind4" + [always for ML_BSSN when 4th order is used] -- cgit v1.2.3