% Example compile-time options: % % BOX_SPATIAL_DIMENSION=2 % SELFGRAVITY_OFF % EOS_TILLOTSON % EOS_ELASTIC % KERNEL_FUNCTION=6 % % (The choice of KERNEL_FUNCTION=6 is optional here. Experiment with different kernels. % But be sure to set DesNumNgb accordingly; for e.g. KERNEL_FUNCTION=5 or 7, use 40, =3, use 20, =4 or 6 use 32) % % Experiment with both HYDRO_MESHLESS_FINITE_MASS and HYDRO_DENSITY_SPH. On some elastic problems, SPH can % give improved problems, though if you have fluid phases, MFM tends to give better behavior. % InitCondFile ring_collision_ics % choose IC and resolution OutputDir output TimeMax 300 TimeBetSnapshot 5 MaxSizeTimestep 0.1 DesNumNgb 32 % make sure to set appropriate for kernel %--- parameters here make EOS: P=cs*cs*(rho-rho_0), where cs=rho_0=1 in code units ---% Tillotson_EOS_params_a 0 Tillotson_EOS_params_b 0 Tillotson_EOS_params_u_0 1e+20 Tillotson_EOS_params_rho_0 1 Tillotson_EOS_params_A 1 Tillotson_EOS_params_B 0 Tillotson_EOS_params_u_s 1e+20 Tillotson_EOS_params_u_s_prime 1e+20 Tillotson_EOS_params_alpha 0 Tillotson_EOS_params_beta 0 Tillotson_EOS_params_mu 0.22 Tillotson_EOS_params_Y0 1e+20 %-- optional [if having trouble on initial iteration to correct kernel sizes], can prompt with %Softening_Type0 0.001 % this is our elastic material