We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We carry out a linear stability analysis of a magnetized relativistic rotating cylindrical jet flow using the approximation of zero thermal pressure. We identify several modes of instability in the jet: Kelvin-Helmholtz, current driven and two kinds of centrifugal-buoyancy modes -- toroidal and poloidal. The Kelvin-Helmholtz mode is found at low magnetization and its growth rate depends very weakly on the pitch parameter of the background magnetic field and on rotation. The current driven mode is found at high magnetization, the values of its growth rate and the wavenumber, corresponding to the maximum growth, increase as we decrease the pitch parameter of the background magnetic field. This mode is stabilized by rotation, especially, at high magnetization. The centrifugal-buoyancy modes, arising due to rotation, tend also to be more stable when magnetization is increased. Overall, relativistic jet flows appear to be more stable with respect to their non-relativistic counterpart.

We studied dynamical balances in magnetorotational instability (MRI) turbulence with a net vertical field in the shearing box model of disks. Analyzing the turbulence dynamics in Fourier (${\bf k}$-)space, we identified three types of active modes that define turbulence characteristics. These modes have lengths similar to the box size, i.e., lie in the small wavenumber region in Fourier space labeled the vital area and are: (i) the channel mode - uniform in the disk plane with the smallest vertical wavenumber,(ii) the zonal flow mode - azimuthally and vertically uniform with the smallest radial wavenumber and (iii) the rest modes. The rest modes comprise those harmonics in the vital area whose energies reach more than $50 \%$ of the maximum spectral energy. The rest modes individually are not so significant compared to the channel and zonal flow modes, however, the combined action of their multitude is dominant over these two modes. These three mode types are governed by interplay of the linear and nonlinear processes, leading to their interdependent dynamics. The linear processes consist in disk flow nonmodality-modified classical MRI with a net vertical field. The main nonlinear process is transfer of modes over wavevector angles in Fourier space - the transverse cascade. The channel mode exhibits episodic bursts supplied by linear MRI growth, while the nonlinear processes mostly oppose this, draining the channel energy and redistributing it to the rest modes. As for the zonal flow, it does not have a linear source and is fed by nonlinear interactions of the rest modes.

Magnetorotational instability (MRI) is one of the fundamental processes in astrophysics, driving angular momentum transport and mass accretion in a wide variety of cosmic objects. Despite much theoretical/numerical and experimental efforts over the last decades, its saturation mechanism and amplitude, which sets the angular momentum transport rate, remains not well understood, especially in the limit of high resistivity, or small magnetic Prandtl numbers typical to interiors (dead zones) of protoplanetary disks, liquid cores of planets and liquid metals in laboratory. Using direct numerical simulations, in this paper we investigate the nonlinear development and saturation properties of the helical magnetorotational instability (HMRI) -- a relative of the standard MRI -- in a magnetized Taylor-Couette flow at very low magnetic Prandtl number (correspondingly at low magnetic Reynolds number) relevant to liquid metals. For simplicity, the ratio of azimuthal field to axial field is kept fixed. From the linear theory of HMRI, it is known that the Elsasser number, or interaction parameter determines its growth rate and plays a special role in the dynamics. We show that this parameter is also important in the nonlinear problem. By increasing its value, a sudden transition from weakly nonlinear, where the system is slightly above the linear stability threshold, to strongly nonlinear, or turbulent regime occurs. We calculate the azimuthal and axial energy spectra corresponding to these two regimes and show that they differ qualitatively. Remarkably, the nonlinear state remains in all cases nearly axisymmetric suggesting that HMRI turbulence is quasi two-dimensional in nature. Although the contribution of non-axisymmetric modes increases moderately with the Elsasser number, their total energy remains much smaller than that of the axisymmetric ones.

Helical and azimuthal magnetorotational instabilities operate in rotating magnetized flows with relatively steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear, which determine the threshold of modal growth of these instabilities, are continuously connected when some axial electrical current is allowed to pass through the rotating fluid. We investigate the nonmodal dynamics of these instabilities arising from the non-normality of shear flow in the local approximation, generalizing the results of the modal approach. It is demonstrated that moderate transient/nonmodal amplification of both types of magnetorotational instability occurs within the Liu limits, where the system is stable according to modal analysis. We show that for the helical magnetorotational instability this magnetohydrodynamic behavior is closely connected with the nonmodal growth of the underlying purely hydrodynamic problem.

We investigate magnetohydrodynamic turbulence driven by the magnetorotational instability (MRI) in Keplerian disks with a nonzero net azimuthal magnetic field using shearing box simulations. As distinct from most previous studies, we analyze turbulence dynamics in Fourier (${\bf k}$-) space to understand its sustenance. The linear growth of MRI with azimuthal field has a transient character and is anisotropic in Fourier space, leading to anisotropy of nonlinear processes in Fourier space. As a result, the main nonlinear process appears to be a new type of angular redistribution of modes in Fourier space -- the \emph{nonlinear transverse cascade} -- rather than usual direct/inverse cascade. We demonstrate that the turbulence is sustained by interplay of the linear transient growth of MRI (which is the only energy supply for the turbulence) and the transverse cascade. These two processes operate at large length scales, comparable to box size and the corresponding small wavenumber area, called \emph{vital area} in Fourier space is crucial for the sustenance, while outside the vital area direct cascade dominates. The interplay of the linear and nonlinear processes in Fourier space is generally too intertwined for a vivid schematization. Nevertheless, we reveal the \emph{basic subcycle} of the sustenance that clearly shows synergy of these processes in the self-organization of the magnetized flow system. This synergy is quite robust and persists for the considered different aspect ratios of the simulation boxes. The spectral characteristics of the dynamical processes in these boxes are qualitatively similar, indicating the universality of the sustenance mechanism of the MRI-turbulence.

To understand the self-sustenance of subcritical turbulence in spectrally stable shear flows, we performed direct numerical simulations of homogeneous shear turbulence for different aspect ratios of the flow domain and analyzed the dynamical processes in Fourier space. There are no exponentially growing modes in such flows and the turbulence is energetically supported only by the linear growth of perturbation harmonics due to the shear flow non-normality. This non-normality-induced, or nonmodal growth is anisotropic in spectral space, which, in turn, leads to anisotropy of nonlinear processes in this space. As a result, a transverse (angular) redistribution of harmonics in Fourier space appears to be the main nonlinear process in these flows, rather than direct or inverse cascades. We refer to this type of nonlinear redistribution as the nonlinear transverse cascade. It is demonstrated that the turbulence is sustained by a subtle interplay between the linear nonmodal growth and the nonlinear transverse cascade that exemplifies a well-known bypass scenario of subcritical turbulence. These two basic processes mainly operate at large length scales, comparable to the domain size. Therefore, this central, small wave number area of Fourier space is crucial in the self-sustenance; we defined its size and labeled it as the vital area of turbulence. Outside the vital area, the nonmodal growth and the transverse cascade are of secondary importance. Although the cascades and the self-sustaining process of turbulence are qualitatively the same at different aspect ratios, the number of harmonics actively participating in this process varies, but always remains quite large. This implies that the self-sustenance of subcritical turbulence cannot be described by low-order models.

We perform a linear stability analysis of magnetized rotating cylindrical jet flows in the approximation of zero thermal pressure. We focus our analysis on the effect of rotation on the current driven mode and on the unstable modes introduced by rotation. We find that rotation has a stabilizing effect on the current driven mode only for rotation velocities of the order of the Alfv\'en velocity. Rotation introduces also a new unstable centrifugal buoyancy mode and the "cold" magnetorotational instability. The first mode is analogous to the Parker instability with the centrifugal force playing the role of effective gravity. The magnetorotational instability can be present, but only in a very limited region of the parameter space and is never dominant. The current driven mode is characterized by large wavelenghts and is dominant at small values of the rotational velocity, while the buoyancy mode becomes dominant as rotation is increased and is characterized by small wavelenghts.

The helical magnetorotational instability is known to work for resistive rotational flows with comparably steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear are continuously connected when some axial electrical current is allowed to flow through the rotating fluid. Using a local approximation we demonstrate that the magnetohydrodynamic behavior of this dissipation-induced instability is intimately connected with the nonmodal growth and the pseudospectrum of the underlying purely hydrodynamic problem.

The mechanism through which meter-sized boulders grow to km-sized planetesimals in protoplanetary discs is a subject of active research, since it is critical for planet formation. To avoid spiralling into the protostar due to aerodynamic drag, objects must rapidly grow from cm-sized pebbles, which are tightly coupled to the gas, to large boulders of 1-100m in diameter. It is already well known that over-densities in the gaseous component of the disc provide potential sites for the collection of solids, and that significant density structures in the gaseous component of the disc (e.g., spiral density waves) can trap solids efficiently enough for the solid component of the disc to undergo further gravitational collapse due to their own self-gravity. In this work, we employ the PENCIL CODE to conduct local shearing sheet simulations of massive self-gravitating protoplanetary discs, to study the effect of anticyclonic transient vortices, or eddies, on the evolution of solids in these discs. We find that these types of structures are extremely efficient at concentrating small and intermediate-sized dust particles with friction times comparable to, or less than, the local orbital period of the disc. This can lead to significant over-densities in the solid component of the disc, with density enhancements comparable to, and even higher, than those within spiral density waves; increasing the rate of gravitational collapse of solids into bound structures.

We find and investigate via numerical simulations self-sustained two-dimensional turbulence in a magnetohydrodynamic flow with a maximally simple configuration: plane, noninflectional (with a constant shear of velocity) and threaded by a parallel uniform background magnetic field. This flow is spectrally stable, so the turbulence is subcritical by nature and hence it can be energetically supported just by transient growth mechanism due to shear flow nonnormality. This mechanism appears to be essentially anisotropic in spectral (wavenumber) plane and operates mainly for spatial Fourier harmonics with streamwise wavenumbers less than a ratio of flow shear to the Alfv\'{e}n speed, $k_y < S/u_A$ (i.e., the Alfv\'{e}n frequency is lower than the shear rate). We focused on the analysis of the character of nonlinear processes and underlying self-sustaining scheme of the turbulence, i.e., on the interplay between linear transient growth and nonlinear processes, in spectral plane. Our study, being concerned with a new type of the energy-injecting process for turbulence -- the transient growth, represents an alternative to the main trends of MHD turbulence research. We find similarity of the nonlinear dynamics to the related dynamics in hydrodynamic flows -- to the \emph{bypass} concept of subcritical turbulence. The essence of the analyzed nonlinear MHD processes appears to be a transverse redistribution of kinetic and magnetic spectral energies in wavenumber plane [as occurs in the related hydrodynamic flow, see Horton et al., Phys. Rev. E {\bf 81}, 066304 (2010)] and differs fundamentally from the existing concepts of (anisotropic direct and inverse) cascade processes in MHD shear flows.

We study particle dynamics in self-gravitating gaseous discs with a simple cooling law prescription via two-dimensional simulations in the shearing sheet approximation. It is well known that structures arising in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results have shown that spiral density waves can be highly efficient at collecting dust particles, creating significant local over-densities of particles. The degree of such concentrations has been shown to be dependent on two parameters: the size of the dust particles and the rate of gas cooling. We expand on these findings, including the self-gravity of dust particles, to see how these particle over-densities evolve. We use the PENCIL CODE to solve the local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas through an aerodynamic drag force. We find that the enhancements in the surface density of particles in spiral density wave crests can reach levels high enough to allow the solid component of the disc to collapse under its own self-gravity. This produces many gravitationally bound collections of particles within the spiral structure. The total mass contained in bound structures appears nearly independent of the cooling time, suggesting that the formation of planetesimals through dust particle trapping by self-gravitating density waves may be possible at a larger range of radii within a disc than previously thought. So, density waves due to gravitational instabilities in the early stages of star formation may provide excellent sites for the rapid formation of many large, planetesimal-sized objects.

We investigate linear dynamics of non-axisymmetric perturbations in incompressible, vertically stratified Keplerian discs with a weak vertical magnetic field in the shearing box approximation. Perturbations are decomposed into shearing waves whose evolution is followed via numerical integration of the linearized ideal MHD equations. There are two basic modes in the system -- inertia-gravity waves and magnetic mode, which displays the magnetorotational instability (MRI). As distinct from previous studies, we introduce `eigenvariables' characterizing each (counter-propagating) component of the inertia-gravity and magnetic modes, which are governed by a set of four first order coupled ordinary differential equations. This allowed us to identify a new process of linear coupling of the two above non-axisymmetric modes due to the disc's differential rotation. We did a comparative analysis of the dynamics of non-axisymmetric and axisymmetric magnetic mode perturbations. It is shown that the growth of optimal and close-to-optimal non-axisymmetric harmonics of this mode, having transient nature, can prevail over the exponential growth of axisymmetric ones (i.e., over the axisymmetric MRI) during dynamical time. A possible implication of this result for axisymmetric channel solutions is discussed. Specifically, the formation of the channel may be affected/impeded by non-axisymmetric modes already at the early linear stage leading to its untimely disruption -- the outcome strongly depends on the amplitude and spectrum of initial perturbation. So, this competition may result in an uncertainty in the magnetic mode's non-linear dynamics. It is also shown that a maximum non-axisymmetric growth is at vertical wavelengths close to the scale-height for which compressibility effects are important. This indirectly suggests that compressibility plays a role in the dynamics of the non-axisymmetric MRI.

We perform a linear analysis of the stability of a magnetized relativistic non-rotating cylindrical flow in the aproximation of zero thermal pressure, considering only the m = 1 mode. We find that there are two modes of instability: Kelvin-Helmholtz and current driven. The Kelvin-Helmholtz mode is found at low magnetizations and its growth rate depends very weakly on the pitch parameter. The current driven modes are found at high magnetizations and the value of the growth rate and the wavenumber of the maximum increase as we decrease the pitch parameter. In the relativistic regime the current driven mode is splitted in two branches, the branch at high wavenumbers is characterized by the eigenfunction concentrated in the jet core, the branch at low wavenumbers is instead characterized by the eigenfunction that extends outside the jet velocity shear region.

We study particle dynamics in local two-dimensional simulations of self-gravitating accretion discs with a simple cooling law. It is well known that the structure which arises in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results using global simulations indicate that spiral density waves are highly efficient at collecting dust particles, creating significant local over-densities which may be able to undergo gravitational collapse. We expand on these findings, using a range of cooling times to mimic the conditions at a large range of radii within the disc. Here we use the Pencil Code to solve the 2D local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas solely through aerodynamic drag force. We find that spiral density waves can create significant enhancements in the surface density of solids, equivalent to 1-10cm sized particles in a disc following the profiles of Clarke (2009) around a solar mass star, causing it to reach concentrations several orders of magnitude larger than the particles mean surface density. We also study the velocity dispersion of the particles, finding that the spiral structure can result in the particle velocities becoming highly ordered, having a narrow velocity dispersion. This implies low relative velocities between particles, which in turn suggests that collisions are typically low energy, lessening the likelihood of grain destruction. Both these findings suggest that the density waves that arise due to gravitational instabilities in the early stages of star formation provide excellent sites for the formation of large, planetesimal-sized objects.

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

Motivated by the recent results of \citet{Lesur_Ogilvie10} on the transport properties of incompressible convection in protoplanetary discs, in this paper we study the role of compressibility and hence of another basic mode -- spiral density waves -- in convective instability in discs. We analyse the linear dynamics of non-axisymmetric convection and spiral density waves in a Keplerian disc with superadiabatic vertical stratification using the local shearing box approach. It is demonstrated that the shear associated with Keplerian differential rotation introduces a novel phenomenon, it causes these two perturbation modes to become coupled: during evolution the convective mode generates (trailing) spiral density waves and can therefore be regarded as a new source of spiral density waves in discs. The wave generation process studied here owes its existence solely to shear of the disc's differential rotation, and is a special manifestation of a more general linear mode coupling phenomena universally taking place in flows with an inhomogeneous velocity profile. We quantify the efficiency of spiral density wave generation by convection as a function of azimuthal and vertical wavenumbers of these modes and find that it is maximal and most powerful when both these length-scales are comparable to the disc scale height. We also show that unlike the convective mode, which tends to transport angular momentum inwards in the linear regime, the spiral density waves transport angular momentum outwards. Based on these findings, we suggest that in the non-linear regime spiral density waves generated by convection may play a role in enhancing the transport of angular momentum due the convective mode alone, which is actually being changed to outward by non-linearity, as indicated by above-mentioned recent developments.

We perform linear analysis of axisymmetric vertical normal modes in stratified compressible self-gravitating polytropic discs in the shearing box approximation. We study specific dynamics for subadiabatic, adiabatic and superadiabatic vertical stratifications. In the absence of self-gravity, four well-known principal modes can be identified in a stratified disc: acoustic p-, surface gravity f-, buoyancy g- and inertial r-modes. After characterizing modes in the non-self-gravitating case, we include self-gravity and investigate how it modifies the properties of these modes. We find that self-gravity, to a certain degree, reduces their frequencies and changes the structure of the dispersion curves and eigenfunctions at radial wavelengths comparable to the disc height. Its influence on the basic branch of the r-mode, in the case of subadiabatic and adiabatic stratifications, and on the basic branch of the g-mode, in the case of superadiabatic stratification (which in addition exhibits convective instability), does appear to be strongest. Reducing the three-dimensional Toomre's parameter Q_{3D} results in the latter modes becoming unstable due to self-gravity, so that they determine the onset criterion and nature of gravitational instability of a stratified disc. By contrast, the p-, f- and convectively stable g-modes, although their corresponding \omega^2 are reduced by self-gravity, never become unstable however small the value of Q_{3D}. This is a consequence of the three-dimensionality of the disc. The eigenfunctions corresponding to the gravitationally unstable modes are intrinsically 3D. We also contrast the more exact instability criterion based on our 3D model with that of density waves in 2D (razor-thin) discs. Based on these findings, we comment on the origin of surface distortions seen in numerical simulations of self-gravitating discs.

(abridged) Vortices are believed to play a role in the formation of km-sized planetesimals. However, vortex dynamics is commonly studied in non-self-gravitating discs. The main goal here is to examine the effects of disc self-gravity on vortex dynamics. For this purpose, we employ the 2D self-gravitating shearing sheet approximation. A simple cooling law with a constant cooling time is adopted, such that the disc settles down into a quasi-steady gravitoturbulent state. In this state, vortices appear as transient structures undergoing recurring phases of formation, growth to sizes comparable to a local Jeans scale and eventual shearing and destruction due to the combined effects of self-gravity and background Keplerian shear. Each phase typically lasts about 2 orbital periods or less. As a result, in self-gravitating discs the overall dynamical picture of vortex evolution is irregular consisting of many transient vortices at different evolutionary stages and, therefore, with various sizes up to the local Jeans scale. Vortices generate density waves during evolution, which turn into shocks. Therefore, the dynamics of density waves and vortices are coupled implying that, in general, one should consider both vortex and spiral density wave modes in order to get a proper understanding of self-gravitating disc dynamics. Our results suggest that given such an irregular and rapidly varying character of vortex evolution in self-gravitating discs, it may be difficult for such vortices to effectively trap dust particles. Further study of the behaviour of dust particles embedded in a self-gravitating gaseous disc is, however, required to strengthen this conclusion.

Vortices are believed to greatly help the formation of km sized planetesimals by collecting dust particles in their centers. However, vortex dynamics is commonly studied in non-self-gravitating disks. The main goal here is to examine the effects of disk self-gravity on the vortex dynamics via numerical simulations. In the self-gravitating case, when quasi-steady gravitoturbulent state is reached, vortices appear as transient structures undergoing recurring phases of formation, growth to sizes comparable to a local Jeans scale, and eventual shearing and destruction due to gravitational instability. Each phase lasts over 2-3 orbital periods. Vortices and density waves appear to be coupled implying that, in general, one should consider both vortex and density wave modes for a proper understanding of self-gravitating disk dynamics. Our results imply that given such an irregular and rapidly changing, transient character of vortex evolution in self-gravitating disks it may be difficult for such vortices to effectively trap dust particles in their centers that is a necessary process towards planet formation.

Flow nonnormality induced linear transient phenomena in thin self-gravitating astrophysical discs are studied in the shearing sheet approximation. The considered system includes two modes of perturbations: vortex and (spiral density) wave. It is shown that self-gravity considerably alters the vortex mode dynamics -- its transient (swing) growth may be several orders of magnitude stronger than in the non-self-gravitating case and 2-3 times larger than the transient growth of the wave mode. Based on this finding, we comment on the role of vortex mode perturbations in a gravitoturbulent state. Also described is the linear coupling of the perturbation modes, caused by the differential character of disc rotation. The coupling is asymmetric -- vortex mode perturbations are able to excite wave mode ones, but not vice versa. This asymmetric coupling lends additional significance to the vortex mode as a participant in spiral density waves and shocks manifestations in astrophysical discs.