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 present a "multipatch" infrastructure for numerical simulation of fluid problems in which sub-regions require different gridscales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client-router-server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. This infrastructure may be used with a wide variety of fluid dynamics codes; the only requirement is that their primary dependent variables be the same in all patches, e.g., fluid mass density, internal energy density, and velocity. Its structure can accommodate either Newtonian or relativistic dynamics. The overhead imposed by this system is both problem- and computer cluster architecture-dependent. Compared to a conventional simulation using the same number of cells and processors, the increase in runtime can be anywhere from negligible to a factor of a few; however, one of the infrastructure's advantages is that it can lead to a very large reduction in the total number of zone-updates.

A tidal disruption event (TDE) ensues when a star passes too close to the supermassive black hole (SMBH) in a galactic center and is ripped apart by the tidal field of the SMBH. The gaseous debris produced in a TDE can power a bright electromagnetic flare as it is accreted by the SMBH; so far, several dozen TDE candidates have been observed. For SMBHs with masses above $\sim 10^7 M_\odot$, the tidal disruption of solar-type stars occurs within ten gravitational radii of the SMBH, implying that general relativity (GR) is needed to describe gravity. Three promising signatures of GR in TDEs are: (1) a super-exponential cutoff in the volumetric TDE rate for SMBH masses above $\sim 10^8 M_\odot$ due to direct capture of tidal debris by the event horizon, (2) delays in accretion disk formation (and a consequent alteration of the early-time light curve) caused by the effects of relativistic precession on stream circularization, and (3) quasi-periodic modulation of X-ray emission due to global precession of misaligned accretion disks and the jets they launch. We review theoretical models and simulations of TDEs in Newtonian gravity, then describe how relativistic modifications give rise to these proposed observational signatures, as well as more speculative effects of GR. We conclude with a brief summary of TDE observations and the extent to which they show indications of these predicted relativistic signatures.

ASASSN-14li is a recently-discovered tidal disruption event with an exceptionally rich data-set: spectra and lightcurves in soft X-rays, UV, optical, and radio. To understand its emission properties in all these bands, we have extended our model for post-tidal disruption accretion and photon production to estimate both soft X-ray radiation produced by the prompt accretion phase and synchrotron emission associated with the bow shock driven through an external medium by the unbound tidal debris, as well as optical and UV light. We find that fiducial values of the stellar mass ($1 M_\odot$) and black hole mass ($10^{6.5} M_{\odot}$) yield: quantitative agreement with the optical/UV luminosity, lightcurve, and color temperature; approximate agreement with the somewhat uncertain soft X-ray spectrum and lightcurve; and quantitative agreement with the radio luminosity, spectrum and lightcurve. Equipartition analysis of the radio data implies that the radio-emitting region expands with a constant speed, and its magnitude is comparable to the speed expected for the unbound stellar ejecta. Both facts provide strong support to our model. We find that the disruption event took place in mid-September 2014. Two independent parameters, the magnitude and logarithmic radial gradient of the ambient gas density near the black hole, must be fit to the data to explain the radio emission; their inferred values are comparable to those found near both Sgr A* and the TDE candidate Swift J1644.

A tidal disruption event (TDE) takes place when a star passes near enough to a massive black hole to be disrupted. About half the star's matter is given elliptical trajectories with large apocenter distances, the other half is unbound. To "circularize", i.e., to form an accretion flow, the bound matter must lose a significant amount of energy, with the actual amount depending on the characteristic scale of the flow measured in units of the black hole's gravitational radius ($\sim 10^{51} (R/1000R_g)^{-1}$~erg). Recent numerical simulations \citep{Shiokawa+2015} have revealed that the circularization scale is close to the scale of the most-bound initial orbits, $\sim 10^3 M_{BH,6.5}^{-2/3} R_g \sim 10^{15} M_{BH,6.5}^{1/3}$~cm from the black hole, and the corresponding circularization energy dissipation rate is $\sim 10^{44} M_{BH,6.5}^{-1/6}$~erg/s. We suggest that the energy liberated during circularization, rather then energy liberated by accretion onto the black hole, powers the observed optical TDE candidates. The observed rise times, luminosities, temperatures, emission radii, and line widths seen in these TDEs \citep[e.g.][]{Arcavi+2014} are all more readily explained in terms of heating associated with circularization than in terms of accretion.

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

Motivated by an improved multi-wavelength observational coverage of the transient sky, we investigate the importance of relativistic effects in disruptions of stars by non-spinning black holes (BHs). This paper focuses on calculating the ballistic rate of return of debris to the black hole as this rate is commonly assumed to be proportional to the light curve of the event. We simulate the disruption of a low mass main sequence star by BHs of varying masses ($10^5,10^6,10^7 M_\odot$) and of a white dwarf by a $10^5 M_\odot$ BH. Based on the orbital energy as well as angular momentum of the debris, we infer the orbital distribution and estimate the return rate of the debris following the disruption. We find two signatures of relativistic disruptions: a gradual rise as well as a delayed peak in the return rate curves relative to their Newtonian analogs. Assuming that the return rates are proportional to the light curves, we find that relativistic effects are in principle measurable given the cadence and sensitivity of the current transient sky surveys. Accordingly, using a simple model of a relativistic encounter with a Newtonian parametric fit of the peak leads to an overestimate in the BH mass by a factor of $\sim {\rm few}\times0.1$ and $\sim {\rm few}$ in the case of the main sequence star and white dwarf tidal disruptions, respectively.

The development of a new generation of theoretical models for tidal disruptions is timely, as increasingly diverse events are being captured in surveys of the transient sky. Recently, Gezari et al. reported a discovery of a new class of tidal disruption events: the disruption of a helium-rich stellar core, thought to be a remnant of a red giant (RG) star. Motivated by this discovery and in anticipation of others, we consider tidal interaction of an RG star with a supermassive black hole (SMBH) which leads to the stripping of the stellar envelope and subsequent inspiral of the compact core toward the black hole. Once the stellar envelope is removed the inspiral of the core is driven by tidal heating as well as the emission of gravitational radiation until the core either falls into the SMBH or is tidally disrupted. In the case of tidal disruption candidate PS1-10jh we find that there is a set of orbital solutions at high eccentricities in which the tidally stripped hydrogen envelope is accreted by the SMBH before the helium core is disrupted. This places the RG core in a portion of parameter space where strong tidal heating can lift the degeneracy of the compact remnant and disrupt it before it reaches the tidal radius. We consider how this sequence of events explains the puzzling absence of the hydrogen emission lines from the spectrum of PS1-10jh and gives rise to its other observational features.

We consider tidal encounters between a white dwarf and an intermediate mass black hole. Both weak encounters and those at the threshold of disruption are modeled. The numerical code combines mesh-based hydrodynamics, a spectral method solution of the self-gravity, and a general relativistic Fermi normal coordinate system that follows the star and debris. Fermi normal coordinates provide an expansion of the black hole tidal field that includes quadrupole and higher multipole moments and relativistic corrections. We compute the mass loss from the white dwarf that occurs in weak tidal encounters. Secondly, we compute carefully the energy deposition onto the star, examining the effects of nonradial and radial mode excitation, surface layer heating, mass loss, and relativistic orbital motion. We find evidence of a slight relativistic suppression in tidal energy transfer. Tidal energy deposition is compared to orbital energy loss due to gravitational bremsstrahlung and the combined losses are used to estimate tidal capture orbits. Heating and partial mass stripping will lead to an expansion of the white dwarf, making it easier for the star to be tidally disrupted on the next passage. Finally, we examine angular momentum deposition. By including the octupole tide, we are able for the first time to calculate deflection of the center of mass of the star and debris. With this observed deflection, and taking into account orbital relativistic effects, we compute directly the change in orbital angular momentum and show its balance with computed spin angular momentum deposition.