[Commits] [svn:einsteintoolkit] www/about/gallery/ (Rev. 940)
knarf at cct.lsu.edu
knarf at cct.lsu.edu
Sun Mar 17 02:33:49 CDT 2013
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> </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
+
More information about the Commits
mailing list