Parameters

Paramters adopted in this simulation aregiven below. For simplicity, the disc is assumed to be isothermal and Shakura-Sunyaev's viscosity prescription is adopted. The origin is at the centre of mass of the system, the semimajor axis is on the $x$-axis, and the orbital plane is on the $xy$-plane. At $t=0$ the neutron star is at apastron (Phase $=$ 0.5).

Stellar parameters
Mass of the Be star $M_{*} =$ $18\,M\odot$	Radius of the Be star $R_{*} =$ $8\,R\odot$
Disk temperature $T_{\rm d} = 13000\,{\rm K}$	Disk thickness $H(R_*)/R_* = 0.02$
Sound speed $c_{\rm s} = 13.4\,{\rm km\,s}^{-1}$	Break-up velocity $V_{\rm K}(R_*) = 655\,{\rm km\,s}^{-1}$
Orbital parameters
Orbital period $P_{\rm orb} = 24.3\,{\rm d}$ ( $= 0.067\,{\rm yr}$)	Orbital eccentricity $e = 0.34$
Mass ratio $q = 0.078$ ( $M_x = 1.41\,M\odot$)	Semimajor axis $a = 11.85\,R_*$ $= 6.6\cdot 10^{12}\,{\rm cm}$
Initial tilt angles: $0^\circ$ about $x$-axis & $0^\circ$ about $y$-axis
Simulation parameters
Unit of time $=$ $24.3\,{\rm d}$ ( $= 0.067\,{\rm yr}$)	Unit of length $=$ $11.85\,R_*$ $= 6.6\cdot 10^{12}\,{\rm cm}$
$\alpha_{\rm SS} = 1$ & $\beta_{\rm SPH} = 0$
$\dot{N}_{\rm injected} = 1.57\cdot 10^{4}$ ($t=0.03$) $\to$ $1.57\cdot 10^{4}$ ($t = 30.02$) ( $r_{\rm injected} = 1.0201R_{*}$)
$N_{\rm SPH}\mbox{ (active)} = 100$ ($t=0$) $\to$ $16250$ ($t = 30.02$)