Source author record

Almog Yalinewich

Almog Yalinewich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Plug Disintegration in GRB Jet Eruption

In this work we consider the eruption of a tenuous relativistic hydrodynamic jet from a dense baryonic envelope. As the jet moves out and away, it carries along and continues to accelerate a layer of baryonic material which we refer to as the plug. We solve the relativistic equations of motion for the trajectory of the plug, and verify it using a relativistic hydrodynamic simulation. We show that under these conditions, the plug breaks up at a radius larger by a factor of a few from the radius of the envelope, due to the onset of the Rayleigh Taylor instability. After breakup the jet continues to accelerate to higher Lorentz factors while the plug fragments maintain a moderate Lorentz factor. The presence of slower moving ejecta can explain late time features of GRBs such as X ray flares without recourse to a long lived engine.

preprint2022arXiv

The Moving Mirror model for Fast Radio Bursts

Recent observations of coherent radiation from the Crab pulsar (Bij et al 2021) suggest the emission is driven by an ultra - relativistic ($γ\sim 10^4$), cold plasma flow. A relativistically expanding plasma shell can compress the ambient magnetic field, like a moving mirror, and thus produce coherent radiation whose wavelength is shorter than that of the ambient medium by $γ^2$. This mechanism has been studied in the past by Colgate and Noerdelinger (1971), in the context of radio loud supernova explosions. In this work we propose that a similar mechanism drives the coherent emission in fast radio bursts. The high Lorenz factors dramatically lower the implied energy and magnetic field requirements, allowing the spin down energy of regular (or even recycled), fast spinning pulsars, rather than slow spinning magnetars, to explain FRBs. We show that this model can explain the frequency and the time evolution of observed FRBs, as well as their duration, energetics and absence of panchromatic counterparts. We also predict that the peak frequency of sub pulses decline with observation time as $ω_{\rm obs} \propto t_{\rm obs}^{-1/2}$. Unfortunately, with current capabilities it is not possible to constrain the shape of the curve $ω_{\rm obs} \left(t_{\rm obs} \right)$. Finally, we find that a variation of this model can explain weaker radio transients, such as the one observed from a galactic magnetar. In this variant, the shock wave produces low frequency photons which are then Compton scattered to the GHz range.

preprint2020arXiv

Nekhoroshev Estimates for the Survival Time of Tightly Packed Planetary Systems

$N$-body simulations of non-resonant tightly-packed planetary systems have found that their survival time (i.e. time to first close encounter) grows exponentially with their interplanetary spacing and planetary masses. Although this result has important consequences for the assembly of planetary systems by giants collisions and their long-term evolution, this underlying exponential dependence is not understood from first principles, and previous attempts based on orbital diffusion have only yielded power-law scalings. We propose a different picture, where large deviations of the system from its initial conditions is due to few slowly developing high-order resonances. Thus, we show that the survival time of the system $T$ can be estimated using a heuristic motivated by Nekhoroshev's theorem, and obtain a formula for systems away from overlapping two-body mean-motion resonances as: $T/P=c_1 \frac{a}{Δa} \exp \left(c_2 \frac{Δa}{a} /μ^{1/4}\right)$, where $P$ is the average Keplerian period, $a$ is the average semi major axis, $Δa\ll a$ is the difference between the semi major axes of neighbouring planets, $μ$ is the planet to star mass ratio, and $c_1$ and $c_2$ are dimensionless constants. We show that this formula is in good agreement with numerical N-body experiments for $c_1=5 \cdot 10^{-4}$ and $c_2=8$.

preprint2020arXiv

Self Similar Adiabatic Strong Explosion in a Medium Gravitationally Free Falling to a Point Mass

We develop a generalisation to the classical Sedov Taylor explosion where the medium free falls to a point mass at the centre of the explosion. To verify our analytic results, we compare them to a suite of numerical simulations. We find that there exists a critical energy below which, instead of propagating outward the shock stalls and collapses under gravity. Furthermore, we find that the value of the critical energy threshold decreases when the adiabatic index increases and material is more evenly distributed within the shocked region. We apply this model to the problem of a shock bounce in core collapse supernova, in which the proto neutron star serves as the point mass. The relation between the threshold energy and the distribution of mass in the shock might help explain how turbulence prevents shock stalling and recession in a core collapse supernova explosion.

preprint2016arXiv

Grid Noise in Moving Mesh Codes: Fixing The Volume Inconsistency Problem

Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to attain second order convergence. The fix is based on the understanding that the volume exchanged between cells should be consistent with the hydrodynamical calculations. We benchmark our fix with three test problems and show that it can significantly improve the computational accuracy. We also examine the effect of initial mesh initialization and present an improved model for the Green-Gauss based gradient estimator.

preprint2015arXiv

Asymptotic Steady State Solution to a Bow Shock with an Infinite Mach Number

The problem of a cold gas flowing past a stationary object is considered. It is shown that at large distances from the obstacle the shock front forms a parabolic solid of revolution. The interior of the shock front is obtained by solution of the hydrodynamic equations in parabolic coordinates. The results are verified with a hydrodynamic simulation. The drag force and expected spectra are calculated for such shock, both in case of an optically thin and thick media. Finally, relations to astrophysical bow shocks and other analytic works on oblique shocks are discussed.

preprint2014arXiv

Angular Density Perturbations to Filled Type I Strong Explosions

In this paper we extend the Sedov - Taylor - Von Neumann model for a strong explosion to account for small angular and radial variations in the density. We assume that the density profile is given by $ρ\left(r,θ,ϕ\right)=kr^{-ω}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(θ,ϕ\right)\right)$, where $\varepsilon\ll1$ and $ω\le\frac{7-γ}{γ+1}$. In order to verify our results we compare them to analytical approximations and full hydrodynamic simulations. We demonstrate how this method can be used to describe arbitrary (not just self similar) angular perturbations. This work complements our previous analysis on radial, spherically symmetric perturbations, and allows one to calculate the response of an explosion to arbitrary perturbations in the upstream density. Together, they settle an age old controversy about the inner boundary conditions

preprint2014arXiv

Atmospheric Mass Loss During Planet Formation: The Importance of Planetesimal Impacts

We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($γ=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $χ_{loss}=(1.71 v_g/v_{esc})^{4.9}$, where $v_{esc}$ is the escape velocity from the target. Yet, the global atmospheric mass loss is a weaker function of the impactor velocity $v_{Imp}$ and mass $m_{Imp}$ and given by $X_{loss} ~ 0.4x+1.4x^2-0.8x^3$ (isothermal atmosphere) and $X_{loss} ~ 0.4x+1.8x^2-1.2x^3$ (adiabatic atmosphere), where $x=(v_{Imp}m/v_{esc}M)$. Atmospheric mass loss due to planetesimal impacts proceeds in two different regimes: 1) Large enough impactors $m > \sqrt{2} ρ_0 (πh R)^{3/2}$ (25~km for the current Earth), are able to eject all the atmosphere above the tangent plane of the impact site, which is $h/2R$ of the whole atmosphere, where $h$, $R$ and $ρ_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 πρ_0 h^3$ (1~km for the current Earth) are only able to eject a fraction of the atmospheric mass above the tangent plane. We find that the most efficient impactors (per unit impactor mass) for atmospheric loss are planetesimals just above that lower limit and that the current atmosphere of the Earth could have resulted from an equilibrium between atmospheric erosion and volatile delivery to the atmosphere from planetesimals. We conclude that planetesimal impacts are likely to have played a major role in atmospheric mass loss over the formation history of the terrestrial planets. (Abridged)

preprint2014arXiv

Balancing the load: A Voronoi based scheme for parallel computations

The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs. We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a 2D moving mesh hydrodynamical code, but it can be extended trivially to other codes in 2D or 3D. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes as well as other applications.

preprint2014arXiv

Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh

We present here RICH, a state of the art 2D hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time advancement scheme as well as a novel parallelization scheme based on Voronoi tessellation. Using our code we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time advancement scheme is more accurate and robust than AREPO's, when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Llyod iterations) and it effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way, and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. Moreover, we show that Voronoi based moving mesh schemes suffer from an error, that is resolution independent, due to inconsistencies between the flux calculation and change in the area of a cell. Our code is publicly available as open source and designed in an object oriented, user friendly way that facilitates incorporation of new algorithms and physical processes.

preprint2011arXiv

Slightly Two or Three Dimensional Self-Similar Solutions

Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.