We conducted three-dimensional hydrodynamical simulations to study the interaction of two late opposite jets with the ejecta of a core collapse supernova (CCSN), and study the bipolar structure that results from this interaction as the jets inflate hot-low-density bubbles. The newly born central object, a neutron star (NS; or a black hole), launches these jets at about 50 to 100 days after explosion. The bubbles cross the photosphere in the polar directions at much earlier times than the regions at the same radii near the equatorial plane. The hot bubbles releases more radiation and the photosphere recedes more rapidly in the tenuous bubble. Our results strengthen earlier claims that were based on toy models that such an interaction might lead to a late peak in the light curve, and that an equatorial observer might see a rapid drop in the light curve. Our results have implications to much earlier jets that explode the star, either jets that the newly born NS launches in a CCSN, or jets that a NS companion that merges with the core of a massive star launches in a common envelope jets supernova (CEJSN) event. Our results add indirect support to the CEJSN scenario for fast blue optical transients, e.g., AT2018cow, ZTF18abvkwla, and CSS161010.

The accretion model for the nineteenth century Great Eruption (GE) of $\eta$ Carinae suggests that mass outflowing from the primary was accreted onto the secondary, and the gravitational energy of that mass accounts for the increase in luminosity and most of the kinetic energy of the ejecta. It further argues that the accretion was accompanied by the ejection of two jets that shaped the bipolar Homunculus nebula. Observations of echos from the GE found emission lines with broad wings suggesting some of the mass in equatorial directions reached more than $10\,000 \,\rm{km\,s^{-1}}$. We run hydrodynamic simulations following periastron passage during the GE, launching jets from the secondary as it accreted gas erupted from the primary. We then follow the interaction of the polar jets with the surrounding primary wind, as they accelerate part of the flow to velocities $ > 10\,000 \,\rm{km\,s^{-1}}$ and deflect it towards lower latitudes. We find that the amount of mass that reached these high velocities during the GE is $M_h \approx 0.02 \,\rm{M_{\odot}}$. This value reaches maximum and then decreases with time. Our simulations agree with previous results of the accretion model from which we estimate $M_h$ taking into account the energy budget of the GE. The accretion model can explain the observations of high velocity gas in light echos with the known two stars, and a triple star system is not required.

We conduct three-dimensional hydrodynamical simulations of two opposite jets that interact with a spherical slow wind that includes a denser shell embedded within it, and obtain a bipolar nebula where each of the two lobes is composed of two connected bubbles and Rayleigh-Taylor instability tongues that protrude from the outer bubble and form the `columns crown'. The jets are launched for a short time of 17 years and inflate a bipolar nebula inside a slow wind. When the bipolar structure encounters the dense shell, the interaction causes each of the two lobes to split to two connected bubbles. The interaction is prone to Rayleigh-Taylor instabilities that form tongues that protrude as columns from the outer bubble. The bases of the columns form a ring on the surface of the outer bubble, and the structure resemble a crown that we term the columns crown. This structure resembles, but is not identical to, the many filaments that protrude from the lobes of the bipolar planetary nebula Menzel~3. We discuss our results in comparison to the structure of Menzel~3 and the ways by which the discrepancies can be reconciled, and possibly turn our failure to reproduce the exact structure of Menzel~3 to a success with jets-shell interaction simulations that include more ingredients.

We conduct three-dimensional hydrodynamical simulations of two opposite jets launched from a binary stellar system into a previously ejected shell and show that the interaction can form barrel-like and H-like shapes in the descendant nebula. Such features are observed in planetary nebulae and supernova remnants. Under our assumption the dense shell is formed by a short instability phase of the giant star as it interacts with a stellar companion, and the jets are then launched by the companion as it accretes mass through an accretion disk from the giant star. We find that the H-shaped and barrel-shaped morphological features that the jets form evolve with time, and that there are complicated flow patterns, such as vortices, instabilities, and caps moving ahead along the symmetry axis. We compare our numerical results with images of 12 planetary nebulae, and show that jet-shell interaction that we simulate can account for the barrel-like or H-like morphologies that are observed in these PNe.

We conduct three-dimensional hydrodynamical simulations of two opposite jets launched obliquely to the orbital plane around an asymptotic giant branch (AGB) star and within its dense wind, and demonstrate the formation of a `messy' planetary nebula (PN), namely, a PN lacking any type of symmetry (highly irregular). In building the initial conditions we assume that a tight binary system orbits the AGB star, and that the orbital plane of the tight binary system is inclined to the orbital plane of the binary system and the AGB star (the triple system plane). We further assume that the accreted mass on to the tight binary system forms an accretion disk around one of the stars, and that the plane of the disk is tilted to the orbital plane of the triple system. The highly asymmetrical and filamentary structure that we obtain support the notion that messy PNe might be shaped by triple stellar systems.

We show that a fast wind that expands into a bipolar nebula composed of two opposite jet-inflated bubbles can form a pair of bipolar rings around giant stars. Our model assumes three mass loss episodes: a spherical slow and dense shell, two opposite jets, and a spherical fast wind. We use the FLASH hydrodynamical code in three-dimensions to simulate the flow, and obtain the structure of the nebula. We assume that the jets are launched from an accretion disk around a stellar companion to the giant star. The accretion disk is assumed to be formed when the primary giant star and the secondary star suffer a strong interaction accompanied by a rapid mass transfer process from the primary to the secondary star, mainly a main sequence star. Later in the evolution the primary star is assumed to shrink and blow a fast tenuous wind that interacts with the dense gas on the surface of the bipolar structure. We assume that the dense mass loss episode before the jets are launched is spherically symmetric. Our results might be applicable to some planetary nebulae, and further emphasize the large variety of morphological features that can be formed by jets. But we could not reproduce some of the properties of the outer rings of SN1987A. It seems that some objects, like SN1987A, require a pre-jets mass loss episode with a mass concentration at mid-latitudes.

We suggest that clumpy-dense outflowing equatorial rings around evolved giant stars, such as in supernova 1987A and the Necklace planetary nebula, are formed by bipolar jets that compress gas toward the equatorial plane. The jets are launched from an accretion disk around a stellar companion. Using the FLASH hydrodynamics numerical code we perform 3D numerical simulations, and show that bipolar jets expanding into a dense spherical shell can compress gas toward the equatorial plane and lead to the formation of an expanding equatorial ring. Rayleigh-Taylor instabilities in the interaction region break the ring to clumps. Under the assumption that the same ring-formation mechanism operates in massive stars and in planetary nebulae, we find this mechanism to be more promising for ring formation than mass loss through the second Lagrangian point. The jets account also for the presence of a bipolar nebula accompanying many of the rings.

We simulate the formation of bipolar planetary nebulae (PNe) through very short impulsive mass ejection events from binary systems, where the asymptotic giant branch (AGB) star ejects a mass shell that is accelerated by jets launched from a compact companion. The acceleration process takes place at very short distances from the binary system, such that the photon-diffusion time is long enough to prevent rapid cooling of the shocked jets' material. When the shocked jets' gas density is lower than the shell density the flow becomes Rayleigh-Taylor unstable and dense clumps are formed in the flow. At later times a PN with clumpy lobes that have a linear distance-velocity relation will be observed. This process might account for the formation of bipolar PNe with clumpy lobes, such as NGC 6302. The energy radiated during the months to years duration of such an event will appear as an intermediate-luminosity optical transient (ILOT)

We perform three-dimensional gas-dynamical simulations and show that the asymmetric morphology of the blue and red-shifted components of the outflow at hundreds of astronomical units (AU) from the massive binary system eta Carinae can be accounted for from the collision of the free primary stellar wind with the slowly expanding dense equatorial gas. Owing to the very complicated structure of the century-old equatorial ejecta, that is not fully spatially resolved by observations, we limit ourselves to modelling the equatorial dense gas by one or two dense spherical clouds. Because of that we reproduce the general qualitative properties of the velocity maps, but not the fine details. The fine details of the velocity maps can be matched by simply structuring the dense ejecta in an appropriate way. The blue and red-shifted components are formed in the post-shock flow of the primary wind, on the two sides of the equatorial plane, respectively. The fast wind from the secondary star plays no role in our model, as for most of the orbital period in our model the primary star is closer to us. The dense clouds are observed to be closer to us than the binary system is, and so in our model the primary star faces the dense equatorial ejecta for the majority of the orbital period.

We study the similarities of jet-medium interactions in several quite different astrophysical systems using 2D and 3D hydrodynamical numerical simulations, and find many similarities. The systems include cooling flow (CF) clusters of galaxies, core-collapse supernovae (CCSNe), planetary nebulae (PNe), and common envelope (CE) evolution. The similarities include hot bubbles inflated by jets in a bipolar structure, vortices on the sides of the jets, vortices inside the inflated bubbles, fragmentation of bubbles to two and more bubbles, and buoyancy of bubbles. The activity in many cases is regulated by a negative feedback mechanism. Namely, higher accretion rate leads to stronger jet activity that in turn suppresses the accretion process. After the jets power decreases the accretion resumes, and the cycle restarts. In the case of CF in galaxies and clusters of galaxies we also study the accretion process, which is most likely by cold clumps, i.e., the cold feedback mechanism. In CF clusters we find that heating of the intra-cluster medium (ICM) is done by mixing hot shocked jet gas with the ICM, and not by shocks. Our results strengthen the jet feedback mechanism (JFM) as a common process in many astrophysical objects.

We perform 3D hydrodynamical numerical simulations of the winds interaction process in the massive binary system $\eta$ Carinae, and find the secondary star to accrete mass from the dense primary wind close to periastron passage. This accretion is thought to result in the spectroscopic event and X-ray minimum observed in the system every revolution. In this study we limit ourselves to explore the role of clumps in the primary wind in triggering the accretion process. We include the gravity of the secondary star and the orbital motion starting 19 days (90 degrees) before periastron passage. The accretion process is triggered by dense clumps that cannot be decelerated by the ram pressure of the secondary wind. The dense clumps are formed by instabilities in the thin dense shell formed by the shocked primary wind gas. We explore the role of the numerical viscosity and some physical parameters on the initiation of the accretion process, and explain the unique properties of $\eta$ Car that allow for the periastron accretion process to occur. The accretion starts about a week before periastron passage, as is required to explain the several weeks long X-ray minimum.

We use recent observations of the He I $\lambda10830 \AA$ absorption line and 3D hydrodynamical numerical simulations of the winds collision, to strengthen the case for an orientation of the semimajor axis of the massive binary system Eta Carinae where the secondary star is toward us at periastron passage. Those observations show that the fast blue absorption component exists for only several weeks prior to the periastron passage. We show that the transient nature of the fast blue absorption component supports a geometry where the fast secondary wind, both pre and post-shock material, passes in front of the primary star near periastron passage.

We conduct 3D numerical simulations of the winds collision process in the massive binary system Eta Carinae, and conclude that accretion occurs during periastron passage. We include radiative cooling of the two winds, one from each star, and the gravity of the secondary and primary stars. Our new numerical finding is that at an orbital separation of ~3-4AU, about three weeks before periastron passage, accretion of dense primary wind gas onto the secondary star begins. To isolate the basic role of the secondary stellar gravity, we neglect the orbital motion and the acceleration zone of the primary wind. Including these effects will strengthen even more our conclusion that accretion near periastron passage of Eta Car is inevitable. Accretion of the primary wind gas onto the secondary star for several weeks near periastron passage accounts for the otherwise puzzling behavior of the binary system near periastron passage.

We conduct numerical simulations of axisymmetrical jets expanding into a spherical AGB slow wind. The three-dimensional flow is simulated with an axially symmetric numerical code. We concentrate on jets that are active for a relatively short time. Our results strengthen other studies that show that jets can account for many morphological features observed in planetary nebulae (PNs). Our main results are as follows. (1) With a single jet's launching episode we can reproduce a lobe structure having a `front-lobe', i.e., a small bulge on the front of the main lobe, such as that in the PN Mz~3. (2) In some runs dense clumps are formed along the symmetry axis, such as those observed in the pre-PN M1-92. (3) The mass loss history of the slow wind has a profound influence on the PN structure. (4) A dense expanding torus (ring; disk) is formed in most of our runs. The torus is formed from the inflated lobes, and not from a separate equatorial mass loss episode. (5) The torus and lobes are formed at the same time and from the same mass loss rate episode. However, when the slow wind density is steep enough, the ratio of the distance divided by the radial velocity is larger for regions closer to the equatorial plane than for regions closer to the symmetry axis. (6) With the short jet-active phase a linear relation between distance and expansion velocity is obtained in many cases. (7) Regions at the front of the lobe are moving sufficiently fast to excite some visible emission lines.

We conduct 2D numerical simulations of jets expanding into the slow wind of asymptotic giant branch stars. We show that the post-shock jets' material can explain the observed extended X-ray emission from some planetary nebulae (PNs). Such jets are thought to shape many PNs, and therefore it is expected that this process will contribute to the X-ray emission from some PNs. In other PNs (not simulated in this work) the source of the extended X-ray emission is the shocked spherical wind blown by the central star. In a small fraction of PNs both sources might contribute, and a two-temperatures gas will fit better the X-ray properties than a one-temperature gas. A spacial separation between these two components is expected.

We conduct 2D axisymmetrical hydrodynamical simulations to investigate the interaction of a collimated fast wind (CFW; wide jets) with a spherical AGB wind. The code includes radiative cooling. We find that the shape of the planetary nebula (PN) is sensitive to the exact mass loss history of the AGB wind, and the opening angle of the CFW. Some typical PN morphologies are obtained, but many other observed morphologies seem to require more ingredients than what we assume in our present simulations, e.g., equatorial AGB wind, and ionization and fast wind during the PN phase. The hot bipolar bubble formed by the jets is an X-ray source.

We propose that the large, radius of ~1000 AU, circumbinary rotating disks observed around some post-asymptotic giant branch (post-AGB) binary stars are formed from slow AGB wind material that is pushed back to the center of the nebula by wide jets. We perform 2D-axisymmetrical numerical simulations of fast and wide jets that interact with the previously ejected slow AGB wind. In each system there are two oppositely launched jets, but we use the symmetry of the problem and simulate only one jet. A large circularization-flow (vortex) is formed to the side of the jet which together with the thermal pressure of the shocked jet material accelerate cold slow-wind gas back to the center from distances of ~1000-10000 AU. We find for the parameters we use that up to 0.001 Mo is back-flowing to the center. We conjecture that the orbital angular momentum of the disk material results from the non-axisymmetric structure of jets launched by an orbiting companion. This conjecture will have to be tested with 3D numerical codes.

We calculate the X-ray emission from both constant and time evolving shocked fast winds blown by the central stars of planetary nebulae (PNs) and compare with observations. Using spherically symmetric numerical simulations with radiative cooling, we calculate the flow structure, and the X-ray temperature and luminosity of the hot bubble formed by the shocked fast wind. We find that a constant fast wind gives results that are very close to those obtained from the self-similar solution. We show that in order for a fast shocked wind to explain the observed X-ray properties of PNs, rapid evolution of the wind is essential. More specifically, the mass loss rate of the fast wind should be high early on when the speed is ~300-700 km/s, and then it needs to drop drastically by the time the PN age reaches ~1000 yr. This implies that the central star has a very short pre-PN (post-AGB) phase.

We calculate the X-ray emission from the shocked fast wind blown by the central stars of planetary nebulae (PNs) and compare with observations. Using spherically symmetric self similar solutions, we calculate the flow structure and X-ray temperature for a fast wind slamming into a previously ejected slow wind. We find that the observed X-ray emission of six PNs can be accounted for by shocked wind segments that were expelled during the early PN phase, if the fast wind speed is moderate, v_2~400-600 km/sec, and the mass loss rate is a few times 10^{-7} Mo/year. We find, as proposed previously, that the morphology of the X-ray emission is in the form of a narrow ring inner to the optical bright part of the nebula. The bipolar X-ray morphology of several observed PNs, which indicates an important role of jets rather than a spherical fast wind, cannot be explained by the flow studied here.

We calculate the X-ray luminosity and light curve for the stellar binary system Eta Carinae for the entire orbital period of 5.54 years. By using a new approach we find, as suggested before, that the collision of the winds blown by the two stars can explain the X-ray emission and temporal behavior. Most X-ray emission in the $2-10 \kev$ band results from the shocked secondary stellar wind. The observed rise in X-ray luminosity just before minimum is due to increase in density and subsequent decrease in radiative cooling time of the shocked fast secondary wind. Absorption, particularly of the soft X-rays from the primary wind, increases as the system approaches periastron and the shocks are produced deep inside the primary wind. However, absorption can not account for the drastic X-ray minimum. The 70 day minimum is assumed to result from the collapse of the collision region of the two winds onto the secondary star. This process is assumed to shut down the secondary wind, hence the main X-ray source. We show that this assumption provides a phenomenological description of the X-ray behavior around the minimum.