Numerical simulations have historically played a major role in understanding the hydrodynamics of the tidal disruption process. Given the complexity of the geometry of the system, the challenges posed by the problem have indeed stimulated much work on the numerical side. Smoothed Particles Hydrodynamics methods, for example, have seen their very first applications in the context of tidal disruption and still play a major role to this day. Likewise, initial attempts at simulating the evolution of the disrupted star with the so-called affine method have been historically very useful. In this Chapter, we provide an overview of the numerical techniques used in the field and of their limitations, and summarize the work that has been done to simulate numerically the tidal disruption process.

We propose a semi-analytical model for the self-intersection of the fallback stream in tidal disruption events (TDEs). When the initial periapsis is less than about 15 gravitational radii, a large fraction of the shocked gas is unbound in the form of a collision-induced outflow (CIO). This is because large apsidal precession causes the stream to self-intersect near the local escape speed at radius much below the apocenter. The rest of the fallback gas is left in more tightly bound orbits and quickly joins the accretion flow. We propose that the CIO is responsible for reprocessing the hard emission from the accretion flow into the optical band. This picture naturally explains the large photospheric radius (or low blackbody temperature) and typical widths of the H and/or He emission lines seen in optical TDEs. We predict the CIO-reprocessed spectrum in the infrared to be L_\nu ~ \nu^{~0.5}, shallower than a blackbody. The partial sky coverage of the CIO also provides a unification of the diverse X-ray behaviors of optical TDEs. According to this picture, optical surveys filter out a large fraction of TDEs with low-mass blackholes due to lack of a reprocessing layer, and the volumetric rate of optical TDEs is nearly flat wrt. the blackhole mass for M < 10^7 solar masses. This filtering causes the optical TDE rate to be lower than the total rate by a factor of ~10 or more. When the CIO is decelerated by the ambient medium, radio emission at the level of that in ASASSN-14li may be produced, but the timescales and peak luminosities can be highly diverse. Finally, our method paves the way for global simulations of the disk formation process by injecting gas at the intersection point according to the prescribed velocity and density profiles.

A star coming too close to a supermassive black hole gets disrupted by the tidal force of the compact object in a tidal disruption event, or TDE. Following this encounter, the debris evolves into an elongated stream, half of which coming back to pericenter. Relativistic apsidal precession then leads to a self-crossing shock that initiates the formation of an accretion disc. We perform the first simulation of this process considering a realistic stellar trajectory and black hole mass, which has so far eluded investigations for computational reasons. This numerical issue is alleviated by using as initial conditions the outflow launched by the self-crossing shock according the local simulation of Lu & Bonnerot (2019). We find that the gas leaving the intersection point experiences numerous secondary shocks that result in the rapid formation of a thick and marginally-bound disc. The mass distribution features two overdensities identified as spiral shocks that drive slow gas inflow along the mid-plane. Inward motion primarily takes place along the funnels of the newly-formed torus, from which a fraction of the matter can get accreted. Further out, the gas moves outward forming an extended envelope completely surrounding the accretion flow. Secondary shocks heat the debris at a rate of a few times $10^{44} \, \rm erg \, s^{-1}$ with a large fraction likely participating to the bolometric luminosity. These results pave the way towards a complete understanding of the early radiation from TDEs that progressively becomes accessible from observations.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

The rate of tidal disruption events (TDEs) can vary by orders of magnitude depending on the environment and the mechanism that launches the stars towards the black hole's vicinity. For the largest rates, two disruptions can take place shortly one after the other in a double TDE. In this case, the two debris streams may collide with each other before falling back to the black hole resulting in an electromagnetic emission that is absent from single TDEs. We analytically evaluate the conditions for this streams collision to occur. It requires that the difference in pericenter location between the two disruptions makes up for the time delay between them. In addition, the width of the streams must compensate for the vertical offset induced by the inclination of their orbital planes. If the double TDE happens following the tidal separation of a binary, we find that the streams can collide with a probability as high as $44\%$. We validate our analytical conditions for streams collision through hydrodynamical simulations and find that the associated shocks heat the gas significantly. If photons are able to rapidly escape, a burst of radiation ensues lasting a few days with a luminosity $\sim 10^{43}\, \rm erg\, s^{-1}$, most likely in the optical band. This signal represents a precursor to the main flare of TDEs that could in particular be exploited to determine the efficiency of disc formation from the stellar debris.

We demonstrate that the compact, thick disc formed in a tidal disruption event may be unstable to non-axisymmetric perturbations in the form of the Papaloizou-Pringle instability. We show this can lead to rapid redistribution of angular momentum that can be parameterised in terms of an effective Shakura-Sunyaev $\alpha$ parameter. For remnants that have initially weak magnetic fields, this may be responsible for driving mass accretion prior to the onset of the magneto-rotational instability. For tidal disruptions around a $10^6$ M$_{\odot}$ black hole, the measured accretion rate is super-Eddington but is not sustainable over many orbits. We thus identify a method by which the torus formed in tidal disruption event may be significantly accreted before the magneto-rotational instability is established.

When a star gets tidally disrupted by a supermassive black hole, its magnetic field is expected to pervade its debris. In this paper, we study this process via smoothed particle magnetohydrodynamical simulations of the disruption and early debris evolution including the stellar magnetic field. As the gas stretches into a stream, we show that the magnetic field evolution is strongly dependent on its orientation with respect to the stretching direction. In particular, an alignment of the field lines with the direction of stretching induces an increase of the magnetic energy. For disruptions happening well within the tidal radius, the star compression causes the magnetic field strength to sharply increase by an order of magnitude at the time of pericentre passage. If the disruption is partial, we find evidence for a dynamo process occurring inside the surviving core due to the formation of vortices. This causes an amplification of the magnetic field strength by a factor of $\sim 10$. However, this value represents a lower limit since it increases with numerical resolution. For an initial field strength of 1 G, the magnetic field never becomes dynamically important. Instead, the disruption of a star with a strong 1 MG magnetic field produces a debris stream within which magnetic pressure becomes similar to gas pressure a few tens of hours after disruption. If the remnant of one or multiple partial disruptions is eventually fully disrupted, its magnetic field could be large enough to magnetically power the relativistic jet detected from Swift J1644+57. Magnetized streams could also be significantly thickened by magnetic pressure when it overcomes the confining effect of self-gravity.

Most massive galaxies are thought to contain a supermassive black hole in their centre surrounded by a tenuous gas environment, leading to no significant emission. In these quiescent galaxies, tidal disruption events represent a powerful detection method for the central black hole. Following the disruption, the stellar debris evolves into an elongated gas stream, which partly falls back towards the disruption site and accretes onto the black hole producing a luminous flare. Using an analytical treatment, we investigate the interaction between the debris stream and the gas environment of quiescent galaxies. Although we find dynamical effects to be negligible, we demonstrate that Kelvin-Helmholtz instability can lead to the dissolution of the stream into the ambient medium before it reaches the black hole, likely dimming the associated flare. This result is robust against the presence of a typical stellar magnetic field and fast cooling within the stream. Furthermore, we find this effect to be enhanced for disruptions involving more massive black holes and/or giant stars. Consequently, although disruptions of evolved stars have been proposed as a useful probe of black holes with masses $\gtrsim 10^8 \, M_{\odot}$, we argue that the associated flares are likely less luminous than expected.

A large number of tidal disruption event (TDE) candidates have been observed recently, often differing in their observational features. Two classes appear to stand out: X-ray and optical TDEs, the latter featuring lower effective temperatures and luminosities. These differences can be explained if the radiation detected from the two categories of events originates from different locations. In practice, this location is set by the evolution of the debris stream around the black hole and by the energy dissipation associated with it. In this paper, we build an analytical model for the stream evolution, whose dynamics is determined by both magnetic stresses and shocks. Without magnetic stresses, the stream always circularizes. The ratio of the circularization timescale to the initial stream period is $t_{\rm ev}/t_{\rm min} = 8.3 (M_{\rm h}/10^6 M_{\odot})^{-5/3} \beta^{-3}$, where $M_{\rm h}$ is the black hole mass and $\beta$ is the penetration factor. If magnetic stresses are strong, they can lead to the stream ballistic accretion. The boundary between circularization and ballistic accretion corresponds to a critical magnetic stresses efficiency $v_{\rm A}/v_{\rm c} \approx 10^{-1}$, largely independent of $M_{\rm h}$ and $\beta$. However, the main effect of magnetic stresses is to accelerate the stream evolution by strengthening self-crossing shocks. Ballistic accretion therefore necessarily occurs on the stream dynamical timescale. The shock luminosity associated to energy dissipation is sub-Eddington but decays as $t^{-5/3}$ only for a slow stream evolution. Finally, we find that the stream thickness rapidly increases if the stream is unable to cool completely efficiently. A likely outcome is its fast evolution into a thick torus, or even an envelope completely surrounding the black hole.

The potential of tidal disruption of stars to probe otherwise quiescent supermassive black holes cannot be exploited, if their dynamics is not fully understood. So far, the observational appearance of these events has been derived from analytical extrapolations of the debris dynamical properties just after disruption. By means of hydrodynamical simulations, we investigate the subsequent fallback of the stream of debris towards the black hole for stars already bound to the black hole on eccentric orbits. We demonstrate that the debris circularize due to relativistic apsidal precession which causes the stream to self-cross. The circularization timescale varies between 1 and 10 times the period of the star, being shorter for more eccentric and/or deeper encounters. This self-crossing leads to the formation of shocks that increase the thermal energy of the debris. If this thermal energy is efficiently radiated away, the debris settle in a narrow ring at the circularization radius with shock-induced luminosities of $\sim 10-10^3 \, L_{\rm Edd}$. If instead cooling is impeded, the debris form an extended torus located between the circularization radius and the semi-major axis of the star with heating rates $\sim 1-10^2 \, L_{\rm Edd}$. Extrapolating our results to parabolic orbits, we infer that circularization would occur via the same mechanism in $\sim 1$ period of the most bound debris for deeply penetrating encounters to $\sim 10$ for grazing ones. We also anticipate the same effect of the cooling efficiency on the structure of the disc with associated luminosities of $\sim 1-10 \, L_{\rm Edd}$ and heating rates of $\sim 0.1-1 \, L_{\rm Edd}$. In the latter case of inefficient cooling, we deduce a viscous timescale generally shorter than the circularization timescale. This suggests an accretion rate through the disc tracing the fallback rate, if viscosity starts acting promptly.