We investigate the feedback of the stellar jets on the surrounding interstellar gas based on 2D and 3D simulations applying HD and MHD module of PLUTO 4.2 code. The main question we address is whether the stellar jet can be considered as a turbulence driver into the interstellar gas. In addition, we investigate the most effective circumstances in which the driven turbulence is larger and can survive for a longer time scale in the ambient gas. We present a case study of different parameters runs including the jet Mach number, the initial jet velocity field and the background magnetic field geometries and the interacting jets. Also, we study the environmental effects on the jet-gas interaction by considering the non-homogeneous surrounding gas containing the clumps in the model setup. Among different setups, we find that for (1) a higher jet Mach number, (2) a rotating jet, 16 pages(3) a jet propagating in a magnetized environment, (4)a jet propagating in a non-homogeneous environment, and (5) the interacting jets more fluctuations and random motions are produced in the entrained gas which can survive for a longer time scale. In addition, we perform the 3D simulations of jet-ambient gas interaction and we find that the amount of (subsonic-supersonic) fluctuations increases compared to the axisymmetric run and the entrained gas gains higher velocities in a 3D run. In total, we confirm the previous finding that the stellar jets can transfer the turbulence on neighboring regions and are not sufficient drivers of the large-scale supersonic turbulence in molecular clouds.

We present fully three-dimensional magnetohydrodynamic jet launching simulations of a jet source orbiting in a binary system. We consider a time-dependent binary gravitational potential, thus all tidal forces that are experienced in the non-inertial frame of the jet-launching primary. We investigate systems with different binary separation, different mass ratio, and different inclination between the disk plane and the orbital plane. The simulations run over a substantial fraction of the binary orbital period. All simulations show similar local and global non-axisymmetric effects such as local instabilities in the disk and the jet, or global features such as disk spiral arms and warps, or a global re-alignment of the inflow-outflow structure. The disk accretion rate is higher than for axisymmetric simulations, most probably due to the enhanced angular momentum transport by spiral waves. The disk outflow leaves the Roche lobe of the primary and becomes disturbed by tidal effects. While a disk-orbit inclination of $10\degr$ still allows for a persistent outflow, an inclination of $30\degr$ does not, suggesting a critical angle in between. For moderate inclination we find indication for jet precession such that the jet axis starts to follow a circular pattern with an opening cone of $\simeq 8\degr$. Simulations with different mass ratio indicate a change of time scales for the tidal forces to affect the disk-jet system. A large mass ratio (a massive secondary) leads to stronger spiral arms, a higher (average) accretion and a more pronounced jet-counter jet asymmetry.

We present results of the first ever three-dimensional (3D) magnetohydrodynamic (MHD) simulations of the accretion-ejection structure. We investigate the 3D evolution of jets launched symmetrically from single stars but also jets from warped disks in binary systems. We have applied various model setups and tested them by simulating a stable and bipolar symmetric 3D structure from a single star-disk-jet system. Our reference simulation maintains a good axial symmetry and also a bipolar symmetry for more than 600 rotations of the inner disk confirming the quality of our model setup. We have then implemented a 3D gravitational potential (Roche potential) due to a companion star and run a variety of simulations with different binary separations and mass ratios. These simulations show typical 3D deviations from axial symmetry, such as jet inclination outside the Roche lobe or spiral arms forming in the accretion disk. In order to find indication for precession effects, we have also run an exemplary parameter setup, essentially governed by a small binary separation of only $\simeq{200}$ inner disk radii. This simulation shows strong indication that we observe the onset of a jet precession caused by the wobbling of the jet launching disk. We estimate the opening angle of the precession cone defined by the lateral motion of the the jet axis of about 4 degree after about 5000 dynamical time steps.

We present MHD simulations exploring the launching, acceleration and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an $\alpha^2\Omega$ mean-field dynamo. The disk magnetization remains on a rather low level, that helps to evolve the simulations for $T > 10,000$ dynamical time steps on a domain extending 1500 inner disk radii. We find a magnetic field of the inner disk similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component that plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and {\em large-scale jet knots} by a {\em disk-intrinsic dynamo} that is time-dependent and regenerates the jet-launching magnetic field.

We present results of magnetohydrodynamic (MHD) simulations of jet formation and propagation, discussing a variety of astrophysical setups. In the first approach we consider simulations of relativistic MHD jet formation, considering jets launched from the surface of a Keplerian disk, demonstrating numerically - for the first time - the self-collimating ability of relativistic MHD jets. We obtain Lorentz factors up to about 10 while acquiring a high degree of collimation of about 1 degree. We then present synchrotron maps calculated from the intrinsic jet structure derived from the MHD jet formation simulation. We finally present (non-relativistic) MHD simulations of jet lauching, treating the transition between accretion and ejection. These setups include a physical magnetic diffusivity which is essential for loading the accretion material onto the outflow. We find relatively high mass fluxes in the outflow, of the order of 20-40 % of the accretion rate.

We investigate the jet launching process from accretion disks extending our recent study (paper I) to a truly bipolar setup. We perform axisymmetric MHD simulations of the disk-jet interaction on a computational domain covering both hemispheres, in particular addressing the question of an intrinsically asymmetric origin of jet / counter jet systems. Treating both hemispheres simultaneously, we overcome the equatorial plane symmetry boundary condition used in most previous studies which naturally fosters a symmetric evolution. For the magnetic diffusivity prescription we apply an alpha-parametrisation, considering both, globally models of diffusivity, and local models. We first approve the quality of our numerical setup by generating perfectly symmetric jets, lasting over a 1000s of dynamical time scales. We then disturb the hemispheric symmetry in the disk, and investigate the subsequent evolution of the outflow. The evolution first leads to a substantial disk warping with electric currents intersecting the equatorial plane. We investigate two models, i) a disk with (initially) different thermal scale height in both hemispheres, and ii) a symmetric disk into which a local disturbance is injected in one hemisphere. In both cases the disk asymmetry results in asymmetric outflows with mass fluxes differing by 10-20%. We find up to 30% difference in mass flux between jet and counter jet for this setup, lasting over 1000s of dynamical time scales (i.e. lasting for the whole simulation). In summary, our results suggest that the jet asymmetries in protostellar and extragalactic jets can indeed be generated intrinsically and maintained over long time by disk asymmetries and the standard jet launching mechanism.

We investigate the launching of jets and outflows from magnetically diffusive accretion disks. Using the PLUTO code we solve the time-dependent resistive MHD equations taking into account the disk and jet evolution simultaneously. The main question we address is which kind of disks do launch jets and which kind of disks do not? In particular, we study how the magnitude and distribution of the (turbulent) magnetic diffusivity affect mass loading and jet acceleration. We have applied a turbulent magnetic diffusivity based on \alpha-prescription, but have also investigate examples where the scale height of diffusivity is larger than that of the disk gas pressure. We further investigate how the ejection efficiency is governed by the magnetic field strength. Our simulations last for up to 5000 dynamical time scales corresponding to 900 orbital periods of the inner disk. As a general result we observe a continuous and robust outflow launched from the inner part of the disk, expanding into a collimated jet of super fast magneto-sonic speed. For long time scales the disk internal dynamics changes, as due to outflow ejection and disk accretion the disk mass decreases. For magneto-centrifugally driven jets we find that for i) less diffusive disks, ii) a stronger magnetic field, iii) a low poloidal diffusivity, or a iv) lower numerical diffusivity (resolution), the mass loading of the outflow is increased - resulting in more powerful jets with high mass flux. For weak magnetization the (weak) outflow is driven by the magnetic pressure gradient. We further investigate the jet asymptotic velocity and the jet rotational velocity in respect of the different launching scenarios. We find a lower degree of jet collimation than previous studies, most probably due to our revised outflow boundary condition.