Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user3D simulations of disc-winds extending radially self-similar MHD models
111
0
0.0
(
0
)
Authors
Matthias Stute
Ask ChatGPT about the research
No Arabic abstract
Disc-winds originating from the inner parts of accretion discs are considered as the basic component of magnetically collimated outflows. The only available analytical MHD solutions to describe disc-driven jets are those characterized by the symmetry of radial self-similarity. However, radially self-similar MHD jet models, in general, have three geometrical shortcomings, (i) a singularity at the jet axis, (ii) the necessary assumption of axisymmetry, and (iii) the non-existence of an intrinsic radial scale, i.e. the jets formally extend to radial infinity. Hence, numerical simulations are necessary to extend the analytical solutions towards the axis, by solving the full three-dimensional equations of MHD and impose a termination radius at finite radial distance. We focus here on studying the effects of relaxing the (ii) assumption of axisymmetry, i.e. of performing full 3D numerical simulations of a disc-wind crossing all magnetohydrodynamic critical surfaces. We compare the results of these runs with previous axisymmetric 2.5D simulations. The structure of the flow in all simulations shows strong similarities. The 3D runs reach a steady state and stay close to axisymmetry for most of the physical quantities, except for the poloidal magnetic field and the toroidal velocity which slightly deviate from axisymmetry. The latter quantities show signs of instabilities, which, however, are confined to the region inside the fast magnetosonic separatrix surface. The forces present in the flow, both of collimating and accelerating nature, are in good agreement in both the 2.5D and the 3D runs. We conclude that the analytical solution behaves well also after relaxing the basic assumption of axisymmetry.
rate research
Read More
We derive a self-similar description for the 2D streamline topology and flow structure of an axi-symmetric, thermally driven wind originating from a disc in which the density is a power law function of radius. Our scale-free solution is strictly only valid in the absence of gravity or centrifugal support; comparison with 2D hydrodynamic simulations of winds from Keplerian discs however demonstrates that the scale-free solution is a good approximation also in the outer regions of such discs, and can provide a reasonable description even for launch radii well within the gravitational radius of the flow. Although other authors have considered the flow properties along streamlines whose geometry has been specified in advance, this is the first isothermal calculation in which the flow geometry and variation of flow variables along streamlines is determined self-consistently. It is found that the flow trajectory is very sensitive to the power-law index of radial density variation in the disc: the steeper the density gradient, the stronger is the curvature of streamlines close to the flow base that is required in order to maintain momentum balance perpendicular to the flow. Steeper disc density profiles are also associated with more rapid acceleration, and a faster fall-off of density, with height above the disc plane. The derivation of a set of simple governing equations for the flow structure of thermal winds from the outer regions of power law discs offers the possibility of deriving flow observables without having to resort to hydrodynamical simulation.
Thermal disc winds occur in many contexts and may be particularly important to the secular evolution and dispersal of protoplanetary discs heated by high energy radiation from their central star. In this paper we generalise previous models of self-similar thermal winds - which have self-consistent morphology and variation of flow variables - to the case of launch from an elevated base and to non-isothermal conditions. These solutions are well-reproduced by hydrodynamic simulations, in which, as in the case of isothermal winds launched from the mid-plane, we find winds launch at the maximum Mach number for which the streamline solutions extend to infinity without encountering a singularity. We explain this behaviour based on the fact that lower Mach number solutions do not fill the spatial domain. We also show that hydrodynamic simulations reflect the corresponding self-similar models across a range of conditions appropriate to photoevaporating protoplanetary discs, even when gravity, centrifugal forces, or changes in the density gradient mean the problem is not inherently scale free. Of all the parameters varied, the elevation of the wind base affected the launch velocity and flow morphology most strongly, with temperature gradients causing only minor differences. We explore how launching from an elevated base affects Ne II line profiles from winds, finding it increases (reduces) the full width at half maximum (FWHM) of the line at low (high) inclination to the line of sight compared with models launched from the disc mid-plane and thus weakens the dependence of the FWHM on inclination.
Coronal rain consists of cool and dense plasma condensations formed in coronal loops as a result of thermal instability. Previous numerical simulations of thermal instability and coronal rain formation have relied on artificially adding a coronal heating term to the energy equation. To reproduce large-scale characteristics of the corona, using more realistic coronal heating prescription is necessary. We analyse coronal rain formation and evolution in a 3-dimensional radiative magnetohydrodynamic simulation spanning from convection zone to corona which is self-consistently heated by magnetic field braiding as a result of convective motions. We investigate the spatial and temporal evolution of energy dissipation along coronal loops which become thermally unstable. Ohmic dissipation in the model leads to the heating events capable of inducing sufficient chromospheric evaporation into the loop to trigger thermal instability and condensation formation. The cooling of the thermally unstable plasma occurs on timescales comparable to the duration of the individual impulsive heating events. The impulsive heating has sufficient duration to trigger thermal instability in the loop but does not last long enough to lead to coronal rain limit cycles. We show that condensations can either survive and fall into the chromosphere or be destroyed by strong bursts of Joule heating associated with a magnetic reconnection events. In addition, we find that condensations can also form along open magnetic field lines.
We study the magnetic field and current structure associated with a coronal loop. Through this we investigate to what extent the assumptions of a force-free magnetic field break down and where they might be justified. We analyse a 3D MHD model of the solar corona in an emerging active region with the focus on the structure of the forming coronal loops. The lower boundary of this simulation is taken from a model of an emerging active region. As a consequence of the emerging magnetic flux and the horizontal motions at the surface a coronal loop forms self-consistently. We investigate the current density along magnetic field inside (and outside) this loop and study the magnetic and plasma properties in and around it. We find that the total current along the loop changes its sign from being antiparallel to parallel to the magnetic field. This is caused by the inclination of the loop together with the footpoint motion. Around the loop the currents form a complex non-force-free helical structure. This is directly related to a bipolar current structure at the loop footpoints at the base of the corona and a local reduction of the background magnetic field (i.e. outside the loop) caused by the plasma flow into and along the loop. The locally reduced magnetic pressure in the loop allows the loop to sustain a higher density, which is crucial for the emission in extreme UV. The acting of the flow on the magnetic field hosting the loop turns out to be also responsible for the observed squashing of the loop. The complex magnetic field and current system surrounding it can be modeled only in 3D MHD models where the magnetic field has to balance the plasma pressure. A 1D coronal loop model or a force-free extrapolation can not capture the current system and the complex interaction of the plasma and the magnetic field in the coronal loop, despite the fact that the loop is under low-$beta$ conditions.
A simple 1D dynamical model of thermally driven disc winds is proposed, based on the results of recent, 2.5D axi-symmetric simulations. Our formulation of the disc wind problem is in the spirit of the original Parker (1958) and Bondi (1952) problems, namely we assume an elementary flow configuration consisting of an outflow following pre-defined trajectories in the presence of a central gravitating point mass. Viscosity and heat conduction are neglected. We consider two different streamline geometries, both comprised of straight lines in the (x,z)-plane: (i) streamlines that converge to a geometric point located at (x,z)=(0,-d) and (ii) streamlines that emerge at a constant inclination angle from the disc midplane (the x-axis, as we consider geometrically thin accretion discs). The former geometry is commonly used in kinematic models to compute synthetic spectra, while the latter, which exhibits self-similarity, is likely unused for this purpose, although it easily can be with existing kinematic models. We make the case that it should be, i.e. that geometry (ii) leads to transonic wind solutions with substantially different properties owing to its lack of streamline divergence. Both geometries can be used to complement recent efforts to estimate photoevaporative mass loss rates from protoplanetary discs. Pertinent to understanding our disc wind results, which are also applicable to X-ray binaries and active galactic nuclei, is a focused discussion on lesser known properties of classic Parker wind solutions. We find that the parameter space corresponding to decelerating Parker wind solutions is made larger due to rotation and leads instead to disc wind solutions that always accelerate after the bulk velocity is slowed to a minimum value. Surprisingly, Keplerian rotation may allow for two different transonic wind solutions for the same physical conditions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
