 a shell. The superscript denotes the time. For instance, $\rho_{j+1/2}^{n+1}$ is the zone-centred density of shell $j+1$ at time $t^{n+1}$, and $r_{j}^{n}$ is the zone-edge-centred radius of shell $j$ at time $t^{n}$.

The Gas Dynamical Conservation Equations
----------------------------------------

Hydrodynamic conservation equations for the baryonic component (eqs. \[\[eq:realhydro\_mass\]\] - \[\[eq:realhydro\_energy\]\]) are solved following the finite-difference scheme by @1995ApJ...442..480T. We first update the velocity and position using the so-called “leap-frog” scheme, so that the velocity and the position are staggered in time: $$v_{j}^{n+1/2} = v_{j}^{n-1/2} -\left[
  4\pi(r_{j}^{n})^{2}\frac{p_{j+1/2}^{n} -p_{j-1/2}^{n}}{dm_{j}}
  +\frac{m_{j}^{n}}{(r_{j}^{n})^{2}}\right] dt^{n},
\label{eq:leapfrog1}$$ and $$r_{j}^{n+1} = r_{j}^{n} + v_{j}^{n+1/2}dt^{n+1/2},
\label{eq:leapfrog2}$$ which are second-order accurate. As the mass of each shell is conserved for such a Lagrangian scheme, density is updated following $$\rho_{j+1/2}^{n+1} = \frac{dm_{j+1/2}}{(4/3) \pi[(r_{j+1}^{n+1})^{3}
    -(r_{j}^{n+1})^{3}]}.$$ In these equations, $$dt^{n}=\frac{1}{2}(dt^{n-1/2}+dt^{n+1/2}),$$ and $$dm_{j}=\frac{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ We then advance the energy by $$\begin{aligned}
e_{i+1/2}^{n+1} &=& e_{i+1/2}^{n} -
p_{i+1/
a shell. The superscript denotes the time. For instance, $ \rho_{j+1/2}^{n+1}$ is the zone - center concentration of shell $ j+1 $ at time $ t^{n+1}$, and $ r_{j}^{n}$ is the zone - edge - center spoke of shell $ j$ at time $ t^{n}$. 

 The Gas Dynamical Conservation Equations 
 ---------------------------------------- 

 Hydrodynamic conservation equations for the baryonic component (equivalent. \[\[eq: realhydro\_mass\]\ ] - \[\[eq: realhydro\_energy\]\ ]) are resolve following the finite - difference scheme by @1995ApJ... 442.. 480T. We first update the speed and position using the so - called “ jump - frog ” scheme, so that the velocity and the placement are staggered in time: $ $ v_{j}^{n+1/2 } = v_{j}^{n-1/2 } -\left [ 
   4\pi(r_{j}^{n})^{2}\frac{p_{j+1/2}^{n } -p_{j-1/2}^{n}}{dm_{j } } 
   + \frac{m_{j}^{n}}{(r_{j}^{n})^{2}}\right ] dt^{n }, 
 \label{eq: leapfrog1}$$ and $ $ r_{j}^{n+1 } = r_{j}^{n } + v_{j}^{n+1/2}dt^{n+1/2 }, 
 \label{eq: leapfrog2}$$ which are second - ordering accurate. As the mass of each shell is conserved for such a Lagrangian outline, density is updated following $ $ \rho_{j+1/2}^{n+1 } = \frac{dm_{j+1/2}}{(4/3) \pi[(r_{j+1}^{n+1})^{3 } 
     -(r_{j}^{n+1})^{3}]}.$$ In these equations, $ $ dt^{n}=\frac{1}{2}(dt^{n-1/2}+dt^{n+1/2}),$$ and $ $ dm_{j}=\frac{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ We then boost the energy by $ $ \begin{aligned } 
 e_{i+1/2}^{n+1 } & = & e_{i+1/2}^{n } - 
 p_{i+1/
 a dhell. The superscript dekotes the time. For instence, $\rhk_{j+1/2}^{n+1}$ is tfe zone-centred density of shxll $h+1$ at upme $t^{n+1}$, and $r_{j}^{n}$ is the zone-edge-bentred rqdiuw of shell $o$ at timc $t^{n}$.

Tgc Gas Vynamical Consetvation Equadions
----------------------------------------

Hydrodynakiz eonservation equations for the baryogic comlojent (eqs. \[\[eq:reajhydgo\_iass\]\] - \[\[eq:realhydro\_energy\]\]) are solved folmowing uhe finite-differenve scheme by @1995ApJ...442..480T. We first updwte the velocity ajd position usigt the so-calldd “leap-frog” scheme, so jhat the velocity and the positiun arz staggered ib tlke: $$v_{j}^{n+1/2} = v_{j}^{n-1/2} -\left[
  4\pi(r_{j}^{n})^{2}\frac{p_{m+1/2}^{m} -p_{j-1/2}^{n}}{dk_{j}}
  +\frav{m_{j}^{n}}{(r_{j}^{n})^{2}}\right] db^{n},
\labxl{eq:oeapfrog1}$$ and $$r_{j}^{n+1} = r_{j}^{n} + v_{j}^{n+1/2}dt^{n+1/2},
\label{eq:leakfrog2}$$ whicv cre second-order accurqtw. As jhe mdss ud exch sielm is clnsxrved for shch a Lagrabgian scheme, densitu yw updated folmowing $$\rro_{j+1/2}^{n+1} = \frac{dm_{j+1/2}}{(4/3) \pi[(r_{j+1}^{n+1})^{3}
    -(r_{j}^{n+1})^{3}]}.$$ In these equdtikns, $$dt^{n}=\frac{1}{2}(dt^{n-1/2}+dt^{n+1/2}),$$ and $$dn_{j}=\frac{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ We then wdvance tre energy by $$\begin{aligned}
e_{i+1/2}^{n+1} &=& e_{i+1/2}^{n} -
p_{i+1/
a shell. The superscript denotes the time. $\rho_{j+1/2}^{n+1}$ the zone-centred of shell $j+1$ is zone-edge-centred radius of $j$ at time The Gas Dynamical Conservation Equations ---------------------------------------- conservation equations for the baryonic component (eqs. \[\[eq:realhydro\_mass\]\] - \[\[eq:realhydro\_energy\]\]) are solved following finite-difference scheme by @1995ApJ...442..480T. We first update the velocity and position using the “leap-frog” so the and the position are staggered in time: $$v_{j}^{n+1/2} = v_{j}^{n-1/2} -\left[ 4\pi(r_{j}^{n})^{2}\frac{p_{j+1/2}^{n} -p_{j-1/2}^{n}}{dm_{j}} +\frac{m_{j}^{n}}{(r_{j}^{n})^{2}}\right] dt^{n}, \label{eq:leapfrog1}$$ $$r_{j}^{n+1} = r_{j}^{n} + v_{j}^{n+1/2}dt^{n+1/2}, \label{eq:leapfrog2}$$ which are accurate. As the mass each shell is conserved for a scheme, density updated $$\rho_{j+1/2}^{n+1} \frac{dm_{j+1/2}}{(4/3) \pi[(r_{j+1}^{n+1})^{3} -(r_{j}^{n+1})^{3}]}.$$ these equations, $$dt^{n}=\frac{1}{2}(dt^{n-1/2}+dt^{n+1/2}),$$ and $$dm_{j}=\frac{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ We then advance the energy by $$\begin{aligned} e_{i+1/2}^{n+1} &=& e_{i+1/2}^{n} - p_{i+1/
 a shell. The superscript denotEs the time. FOr insTanCe, $\rHo_{J+1/2}^{n+1}$ is The zOne-centred densITy of Shell $j+1$ at time $t^{n+1}$, and $r_{j}^{n}$ is The zoNe-EDge-cENtRed raDius of sHElL $J$ At tImE $t^{N}$.

ThE GAS DYnamiCal conservAtion EquatIonS
----------------------------------------

HYdrodynamic cONsErvation eqUatIons for the baRyoNic comPoNenT (Eqs. \[\[eq:ReaLhydrO\_mass\]\] - \[\[eQ:RealhyDro\_energy\]\]) ArE Solved FOllowinG THe FiniTe-difference schemE By @1995aPJ...442..480T. We first updaTe the vElOCiTY And PosItion using ThE so-caLLed “leap-FRoG” SCHemE, So that the veloCity and the pOSitIon are StAggERed in tIme: $$v_{j}^{N+1/2} = v_{J}^{N-1/2} -\leFt[
  4\pi(r_{j}^{n})^{2}\frac{P_{j+1/2}^{n} -p_{J-1/2}^{n}}{dm_{j}}
  +\frac{M_{j}^{n}}{(r_{j}^{n})^{2}}\RIght] dt^{n},
\LAbel{eq:lEapfroG1}$$ anD $$r_{j}^{N+1} = r_{j}^{n} + V_{J}^{n+1/2}Dt^{N+1/2},
\laBeL{Eq:lEApFroG2}$$ WhiCh are secOnD-oRder aCcurATE. aS the MasS of eAch shEll is conserveD foR sucH A LaGrangIan scHeme, DeNsity Is updaTed foLlOwing $$\rho_{j+1/2}^{n+1} = \frac{dM_{j+1/2}}{(4/3) \pi[(R_{j+1}^{n+1})^{3}
    -(r_{j}^{n+1})^{3}]}.$$ In tHesE eQuaTiOns, $$dt^{N}=\Frac{1}{2}(dt^{N-1/2}+dt^{N+1/2}),$$ anD $$dm_{j}=\fraC{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ wE thEn ADVAnCe the energy by $$\begin{AlIGNeD}
e_{i+1/2}^{n+1} &=& e_{i+1/2}^{n} -
p_{I+1/
 a shell. The superscriptdenotes th e tim e.For  i nsta nce,  $\rho_{j+1/2} ^ {n+1 }$ is the zone-centred  dens it y  ofs he ll $j +1$ att im e $t^ {n +1 }$,  a n d$r_{j }^{ n}$ isthe zone-e dge -c entred radiu s  o f shell $j $ a t time $t^{n }$. 

TheGa s D y namic alConse rvatio n  Equat ions
---- -- - ------ - ------- - - -- ---- -----------

Hydr o dy n amic conservat ion eq ua t io n s  fo r t he baryoni ccompo n ent (eq s .\ [ \ [eq : realhydro\_ma ss\]\] - \[ \ [eq :realh yd ro\ _ energy \]\])  a r e s olved follo wing  the fini te-dif f erences cheme b y @199 5Ap J.. .442 . .4 80 T.We fir s tupd a tethe velo ci ty  andposi t i o n  usi ngtheso-ca lled “leap-fr og”  sch e me,  so t hat t he v el ocity  and t he po si tion are stagge redin time:$$v _{ j}^ {n +1/2} = v_{j }^{ n-1 /2} -\l eft[
 4 \pi (r _ { j }^ {n})^{2}\frac{p_{j +1 / 2 }^ {n} -p_{ j-1/2} ^ {n }} { dm_{j}} +\f rac{ m _ {j}^{ n}}{ ( r_ {j}^{n}) ^{2}}\ r ig ht ] dt^{n }, 
\labe l{ eq: lea pfrog 1 }$$and $$ r_{j}^{n +1} = r_{j}^{n} + v_ { j}^{n+1/2}dt^ { n+ 1 / 2} , 
\la bel {eq:leapfro g2}$ $  whi ch a r esec o nd-or der a cc u ra t e. As the mass of e ac h shel l isconserved for  such a La g r a ngian sc heme ,  d e nsity is updat ed fo llowing $$ \ rho_{j+1 /2}^{ n+1} = \ frac{dm_{ j + 1/2}}{(4 /3)  \p i[( r_{ j + 1} ^{n+1})^{3}
 -(r _{ j}^{n+1 })^ {3}]}.$ $ I n t hes e e qu ations, $ $dt^{n}= \f ra c{ 1} {2} (dt^{ n -1/2}+dt ^{ n+1 /2 }), $$ an d  $$dm_ {j}=\ frac {1 }{ 2 }(d m_{j-1/ 2 }+ d m _{j+ 1/ 2} ).$$  We  t hen a dvan c e t he ener gy by $$\ beg i n{al ig ne d}
e_{i +1/2}^{n+1} & =&  e_{i+1/2} ^{ n}-
p_{i + 1 /
 a_shell. The_superscript denotes the time._For instance,_$\rho_{j+1/2}^{n+1}$_is the_zone-centred_density of shell_$j+1$ at time_$t^{n+1}$, and $r_{j}^{n}$ is_the zone-edge-centred radius_of_shell $j$ at time $t^{n}$.

The Gas Dynamical Conservation Equations
----------------------------------------

Hydrodynamic conservation equations for the baryonic_component_(eqs. \[\[eq:realhydro\_mass\]\]_-_\[\[eq:realhydro\_energy\]\])_are solved following the finite-difference_scheme by @1995ApJ...442..480T. We first_update the_velocity and position using the so-called “leap-frog” scheme,_so_that the velocity_and the position are staggered in time: $$v_{j}^{n+1/2} =_v_{j}^{n-1/2} -\left[
  4\pi(r_{j}^{n})^{2}\frac{p_{j+1/2}^{n} -p_{j-1/2}^{n}}{dm_{j}}
 _+\frac{m_{j}^{n}}{(r_{j}^{n})^{2}}\right] dt^{n},
\label{eq:leapfrog1}$$ and_$$r_{j}^{n+1}_=_r_{j}^{n} + v_{j}^{n+1/2}dt^{n+1/2},
\label{eq:leapfrog2}$$ which_are second-order accurate. As the mass_of each shell is conserved for_such a Lagrangian scheme, density is updated_following $$\rho_{j+1/2}^{n+1} = \frac{dm_{j+1/2}}{(4/3) \pi[(r_{j+1}^{n+1})^{3}
 _  -(r_{j}^{n+1})^{3}]}.$$ In these_equations, $$dt^{n}=\frac{1}{2}(dt^{n-1/2}+dt^{n+1/2}),$$_and $$dm_{j}=\frac{1}{2}(dm_{j-1/2}+dm_{j+1/2}).$$ We then advance_the energy by_$$\begin{aligned}
e_{i+1/2}^{n+1} &=&_e_{i+1/2}^{n} -
p_{i+1/