User: knarf Date: 2013/03/17 02:33 AM

Added: /about/gallery/ tov.par

Modified: /about/gallery/ index.php

Log: add tov example

File Changes:

Directory: /about/gallery/ ==========================

File [modified]: index.php Delta lines: +15 -9 =================================================================== --- about/gallery/index.php 2013-03-01 18:49:24 UTC (rev 939) +++ about/gallery/index.php 2013-03-17 07:33:48 UTC (rev 940) @@ -16,20 +16,26 @@ <ul> <li> Single black hole - <p>This simulation shows how to evolve a stable, single black hole. - Nothing much happens, apart from ...</p> - <table> - <tr><th>Parameter file</th><td>&nbsp;</td></tr> - <tr><th>Thornlist</th><td><a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2012_11/einsteintoolkit.th">ET Ørsted release</a></td></tr> - <tr><th>approx. memory</th><td>X GB</td></tr> - <tr><th>approx. runtime</th><td>X SU</td></tr> - </table> - </li> <li> Binary black hole inspiral </li> <li> Single, stable neutron star + <p>This simulation shows how to evolve a stable, single neutron star. It uses a TOV solution + as initial data and follows the initial (numerical) perturbations, writing the maximum + density to disk. What can be seen is that the oscillates for a few cycles (only very + slitely because the solution is stationary and the induced numerical perturbation is small), + before numerical errors due to the lack of resolution spoil the results. The size of this + simulation is intentionally kept as small as possible, thus the resolution is about as low + as youcan go without loosing any confidence in at least quantitatively correct results. + It is intended to be used as quick examples in lectures, workshops or summer schools where + both time and computational power are typically scarce.</p> + <table> + <tr><th>Parameter file</th><td><a href="tov.par">tov.par</a></td></tr> + <tr><th>Thornlist</th><td><a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2012_11/einsteintoolkit.th">ET Ørsted release</a></td></tr> + <tr><th>approx. memory</th><td>0.25 GB</td></tr> + <tr><th>approx. runtime</th><td>2.5 SU</td></tr> + </table> </li> <li> Single, collapsing neutron star

File [added]: tov.par Delta lines: +291 -0 =================================================================== --- about/gallery/tov.par (rev 0) +++ about/gallery/tov.par 2013-03-17 07:33:48 UTC (rev 940) @@ -0,0 +1,291 @@ +# Example parameter file for a static TOV star. Everything is evolved, but +# because this is a solution to the GR and hydro equations, nothing changes +# much. What can be seen is the initial perturbation (due to numerical errors) +# ringing down (look at the density maximum), and later numerical errors +# governing the solution. Try higher resolutions to decrease this error. + +# Some basic stuff +ActiveThorns = "Time MoL" +ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal" +ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros" +ActiveThorns = "IOUtil" +ActiveThorns = "Formaline" +ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce aeilocalinterp LoopControl" +ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp" +ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic" + +# Finalize +Cactus::terminate = "time" +Cactus::cctk_final_time = 800 # devide by ~203 to get ms + +# Termination Trigger +ActiveThorns = "TerminationTrigger" +TerminationTrigger::max_walltime = 24 # hours +TerminationTrigger::on_remaining_walltime = 15 # minutes +TerminationTrigger::check_file_every = 512 +TerminationTrigger::termination_file = "TerminationTrigger.txt" +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" + +# grid parameters +Carpet::domain_from_coordbase = "yes" +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" +CoordBase::xmin = 0.0 +CoordBase::ymin = 0.0 +CoordBase::zmin = 0.0 +CoordBase::xmax = 48.0 +CoordBase::ymax = 48.0 +CoordBase::zmax = 48.0 +# Change these parameters to change resolution. The ?max settings above +# have to be multiples of these. 'dx' is the size of one cell in x-direction. +# Making this smaller means using higher resolution, because more points will +# be used to cover the same space. +CoordBase::dx = 4.0 +CoordBase::dy = 4.0 +CoordBase::dz = 4.0 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::num_levels_1 = 3 +CarpetRegrid2::radius_1[1] = 24.0 +CarpetRegrid2::radius_1[2] = 12.0 + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + + +ActiveThorns = "ReflectionSymmetry" + +ReflectionSymmetry::reflection_x = "yes" +ReflectionSymmetry::reflection_y = "yes" +ReflectionSymmetry::reflection_z = "yes" +ReflectionSymmetry::avoid_origin_x = "no" +ReflectionSymmetry::avoid_origin_y = "no" +ReflectionSymmetry::avoid_origin_z = "no" + +# storage and coupling +TmunuBase::stress_energy_storage = yes +TmunuBase::stress_energy_at_RHS = yes +TmunuBase::timelevels = 1 +TmunuBase::prolongation_type = none +# This is still necessary for the ET_2012_11 release, but is the default +# after that. +TmunuBase::support_old_CalcTmunu_mechanism = "no" + + +HydroBase::timelevels = 3 + +ADMMacros::spatial_order = 4 + +SpaceMask::use_mask = "yes" + +Carpet::enable_all_storage = no +Carpet::use_buffer_zones = "yes" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" + +Carpet::init_3_timelevels = no +Carpet::init_fill_timelevels = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +# system specific Carpet paramters +Carpet::max_refinement_levels = 10 +driver::ghost_size = 3 +Carpet::prolongation_order_space = 3 +Carpet::prolongation_order_time = 2 + +# Time integration +time::dtfac = 0.25 + +MoL::ODE_Method = "rk4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +# check all physical variables for NaNs +# This can save you computing time, so it's not a bad idea to do this +# once in a whioe. +ActiveThorns = "NaNChecker" +NaNChecker::check_every = 16384 +NaNChecker::action_if_found = "terminate" #"terminate", "just warn", "abort" +NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" + +# Hydro paramters + +ActiveThorns = "EOS_Omni GRHydro" + +HydroBase::evolution_method = "GRHydro" + +GRHydro::riemann_solver = "HLLE" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" +GRHydro::recon_method = "ppm" +GRHydro::GRHydro_stencil = 3 +GRHydro::bound = "none" +GRHydro::rho_abs_min = 1.e-10 +GRHydro::GRHydro_atmo_tolerance = 1.e-3 + +# Curvature evolution parameters + +ActiveThorns = "GenericFD NewRad" +ActiveThorns = "ML_BSSN ML_BSSN_Helper" +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::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftBCoeff = 1 +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 2.66 +ML_BSSN::LapseAdvectionCoeff = 0.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" + +# Some dissipation to get rid of high-frequency noise +ActiveThorns = "SphericalSurface Dissipation" +Dissipation::verbose = "no" +Dissipation::epsdis = 0.01 +Dissipation::vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_curv + ML_BSSN::ML_trace_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift +" + + +# init parameters +InitBase::initial_data_setup_method = "init_some_levels" + +# Use TOV as initial data +ActiveThorns = "TOVSolver" + +HydroBase::initial_hydro = "tov" +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +# Parameters for initial star +TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::TOV_Gamma[0] = 2.75 +TOVSolver::TOV_K[0] = 30000 + +# Set equation of state for evolution +EOS_Omni::poly_gamma = 2.75 +EOS_Omni::poly_gamma_ini = 2.75 +EOS_Omni::poly_k = 30000 +EOS_Omni::gl_gamma = 2.75 +EOS_Omni::gl_k = 30000 + +# I/O + +cactus::cctk_timer_output = "full" + +# Use (create if necessary) an output directory named like the +# parameter file (minus the .par) +IO::out_dir = ${parfile} + +# Write one file overall per output (variable/group) +# In production runs, comment this or set to "proc" to get one file +# per MPI process +IO::out_mode = "onefile" + +# Some screen output +IOBasic::outInfo_every = 16 +IOBasic::outInfo_vars = "Carpet::physical_time_per_hour HydroBase::rho{reductions='maximum'}" + +# Scalar output +IOScalar::outScalar_every = 512 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "norm1 norm2 norm_inf sum maximum minimum" +IOScalar::outScalar_vars = " + HydroBase::rho{reductions='maximum'} + HydroBase::press{reductions='maximum'} + HydroBase::eps{reductions='minimum maximum'} + HydroBase::vel{reductions='minimum maximum'} + HydroBase::w_lorentz{reductions='minimum maximum'} + ADMBase::lapse{reductions='minimum maximum'} + ADMBase::shift{reductions='minimum maximum'} + ML_BSSN::ML_Ham{reductions='norm1 norm2 maximum minimum norm_inf'} + ML_BSSN::ML_mom{reductions='norm1 norm2 maximum minimum norm_inf'} + GRHydro::dens{reductions='minimum maximum sum'} + Carpet::timing{reductions='average'} +" + +# 1D ASCII output. Disable for production runs! +IOASCII::out1D_every = 512 +IOASCII::one_file_per_group = yes +IOASCII::output_symmetry_points = no +IOASCII::out1D_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ML_BSSN::ML_Ham + ML_BSSN::ML_mom +" + +# 2D HDF5 output +CarpetIOHDF5::output_buffer_points = "no" + +CarpetIOHDF5::out2D_every = 32768 +CarpetIOHDF5::out2D_vars = " + HydroBase::rho + HydroBase::eps + HydroBase::vel + HydroBase::w_lorentz + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + " + +# Checkpointing options +IOHDF5::checkpoint = "yes" +IO::checkpoint_dir = $parfile +IO::recover_dir = $parfile +IO::recover = "autoprobe" +IO::checkpoint_ID = "yes" +IO::checkpoint_every = 1048576 +IO::checkpoint_keep = 3 +IO::checkpoint_on_terminate = "yes" + +# Enable to get detailed timing information +#ActiveThorns = "TimerReport" +#TimerReport::out_every = 1024 +#TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::output_all_timers = "yes" + +# Enable for profiling +#Carpet::output_timers_every = 1024 +