 good conservation requires the corrected terms above.

Implicitly, we have been working in the frame moving with the feedback “source” (${\bf x}_{a}={\bf 0}$, ${\bf v}_{a} \equiv d{\bf x}_{a}/dt = {\bf 0}$), in which the source is statistically isotropic. However, in coupling the fluxes to surrounding gas elements, we also must account for the frame motion. Boosting back to the lab/simulation frame, the total ejecta velocity entering an element is of course $\Delta m_{b}^{-1}\,\Delta {\bf p}_{b} + {\bf v}_{a}$. This change of frame has no effect on the mass fluxes, but it does modify the momentum and energy fluxes: to be properly conservative, we must take: $$\begin{aligned}
\label{eqn:flux.mz.framecorr} \Delta m_{b}^{\prime} &\equiv \Delta m_{b}\ \ \, \ \ \  \Delta m_{Z,\,b}^{\prime} \equiv \Delta m_{Z,\,b} \\ 
\label{eqn:flux.p.framecorr} \Delta {\bf p}_{b}^{\prime} &\equiv \Delta {\bf p}_{b} + \Delta m_{b}\,{\bf v}_{a} \\ 
\label{eqn:flux.e.framecorr} \Delta E_{b}^{\prime} &\equiv \Delta {E}_{b} + \frac{1}{2\,\Delta m_{b}}\,\left( |\Delta {\bf p}_{b}^{\prime}|^{2} - |\Delta {\bf p}_{b}|^{2} \right)\end{aligned}$$ where the prime (e.g. “$ \Delta m_{b}^{\prime}$”) notation denotes the lab frame. Note that the extra momentum added to the neighbors ($\sum_{b}\,\Delta m_{b}\,{\bf v}_{a} = m_{\rm ej}\,{\bf v}_{a}$) is exactly the momentum lost by the feedback source $a$, by virtue of its losing $m_{\rm ej}$ in mass.[^9]

These fluxes are simply added to each neighbor in a fully-conservative manner: $$\begin{aligned}
\label{eqn:flux.m.coupling} m_{b}^{\rm new} &= m_{b} + \Delta m_{b}^{\prime} \\ 
\label{eqn:flux.z.coupling
good conservation requires the corrected terms above. 

 Implicitly, we have been work in the frame of reference moving with the feedback “ source ” ($ { \bf x}_{a}={\bf 0}$, $ { \bf v}_{a } \equiv d{\bf x}_{a}/dt = { \bf 0}$), in which the source is statistically isotropic. However, in couple the fluxes to surrounding natural gas element, we also must report for the frame of reference motion. boost back to the lab / simulation frame, the total ejecta speed entering an element is of course $ \Delta m_{b}^{-1}\,\Delta { \bf p}_{b } + { \bf v}_{a}$. This change of frame has no effect on the batch fluxes, but it does modify the momentum and energy flux: to be properly conservative, we must take: $ $ \begin{aligned } 
 \label{eqn: flux.mz.framecorr } \Delta m_{b}^{\prime } & \equiv \Delta m_{b}\ \ \, \ \ \   \Delta m_{Z,\,b}^{\prime } \equiv \Delta m_{Z,\,b } \\ 
 \label{eqn: flux.p.framecorr } \Delta { \bf p}_{b}^{\prime } & \equiv \Delta { \bf p}_{b } + \Delta m_{b}\,{\bf v}_{a } \\ 
 \label{eqn: flux.e.framecorr } \Delta E_{b}^{\prime } & \equiv \Delta { E}_{b } + \frac{1}{2\,\Delta m_{b}}\,\left (|\Delta { \bf p}_{b}^{\prime}|^{2 } - |\Delta { \bf p}_{b}|^{2 } \right)\end{aligned}$$ where the prime (for example   “ $ \Delta m_{b}^{\prime}$ ”) notation denote the lab frame. Note that the extra momentum add to the neighbors ($ \sum_{b}\,\Delta m_{b}\,{\bf v}_{a } = m_{\rm ej}\,{\bf v}_{a}$) is exactly the momentum lost by the feedback source $ a$, by virtue of its lose $ m_{\rm ej}$ in mass.[^9 ] 

 These fluxes are simply added to each neighbor in a fully - conservative manner: $ $ \begin{aligned } 
 \label{eqn: flux.m.coupling } m_{b}^{\rm new } & = m_{b } + \Delta m_{b}^{\prime } \\ 
 \label{eqn: flux.z.coupling
 gold conservation requires the corrected jeems abmve.

Impmicitly, de have been working in the hramw movung with the feedback “rource” (${\bf x}_{a}={\bf 0}$, ${\bd v}_{a} \wquiv d{\bf e}_{z}/dt = {\bf 0}$), in wglch tke source is stajistically ivotropic. Howevar, iu coupling the fluxes to surrounding gas elrmfnts, we also most abcjunt for the frame motion. Boosting baci to tht lab/simulation frsme, the total ejecta veloclty fntering an elemenh is of coutae $\Qwlta m_{b}^{-1}\,\Delta {\bf p}_{b} + {\bf v}_{a}$. This chahge of frame has no effect on tfe maxs fluxes, vur ih does modifb the iomentum and energy xluxes: yo be properly coisercative, we must take: $$\bxgin{aligned}
\label{eqn:fjux.mz.framackrr} \Delta m_{b}^{\prime} &\ewuiv \Geltd m_{b}\ \ \, \ \ \  \Sektz m_{Z,\,b}^{\pgimx} \equiv \Delfa m_{Z,\,b} \\ 
\labwl{eqn:flux.p.framecorr} \Dqota {\bf p}_{b}^{\prims} &\equid \Qelta {\bf p}_{b} + \Delta m_{b}\,{\bf v}_{a} \\ 
\label{eqn:fluq.e.frzmecorr} \Delta E_{b}^{\prime} &\ewuiv \Delta {E}_{b} + \frac{1}{2\,\Depta m_{b}}\,\lefe( |\Delta {\bf p}_{b}^{\prime}|^{2} - |\Delta {\bf p}_{b}|^{2} \right)\end{aligned}$$ wvere vhd pxlid (e.h. “$ \Delta m_{b}^{\prime}$”) notation denotes the lab framq. Nptv that the extra iomentum adcef yj the neighbots ($\sum_{b}\,\Bslfa m_{b}\,{\bf v}_{a} = m_{\rm ej}\,{\hf v}_{a}$) if exaxtly the iomemtum lost by the feedback siurce $a$, by vprtuw of its losing $m_{\rl ej}$ in masr.[^9]

Thexe flixes are simply added tu eadh neighbor in a fulmh-conservative mavneg: $$\becin{alignta}
\label{eqn:flux.m.cou[ling} m_{b}^{\rn nec} &= m_{b} + \Ddlta m_{b}^{\priie} \\ 
\label{een:fluw.s.coupling
good conservation requires the corrected terms above. have working in frame moving with 0}$, v}_{a} \equiv d{\bf = {\bf 0}$), which the source is statistically isotropic. in coupling the fluxes to surrounding gas elements, we also must account for frame motion. Boosting back to the lab/simulation frame, the total ejecta velocity entering element of $\Delta {\bf p}_{b} + {\bf v}_{a}$. This change of frame has no effect on the mass fluxes, it does modify the momentum and energy fluxes: be properly conservative, we take: $$\begin{aligned} \label{eqn:flux.mz.framecorr} \Delta m_{b}^{\prime} \Delta \ \, \ \Delta \equiv \Delta m_{Z,\,b} \label{eqn:flux.p.framecorr} \Delta {\bf p}_{b}^{\prime} &\equiv \Delta {\bf p}_{b} + \Delta m_{b}\,{\bf v}_{a} \\ \label{eqn:flux.e.framecorr} \Delta E_{b}^{\prime} &\equiv {E}_{b} + |\Delta {\bf - {\bf \right)\end{aligned}$$ where the “$ \Delta m_{b}^{\prime}$”) notation denotes the that the extra momentum added to the neighbors m_{b}\,{\bf v}_{a} m_{\rm ej}\,{\bf v}_{a}$) is exactly the lost by the feedback source $a$, by virtue its losing $m_{\rm ej}$ in mass.[^9] These fluxes are simply added to each neighbor in manner: $$\begin{aligned} \label{eqn:flux.m.coupling} m_{b}^{\rm &= m_{b} + m_{b}^{\prime} \label{eqn:flux.z.coupling
 good conservation requires tHe correcteD termS abOve.

imPlicItly, We have been workINg in The frame moving with the fEedbaCk “SOurcE” (${\Bf X}_{a}={\bf 0}$, ${\bF v}_{a} \equiV D{\bF X}_{A}/dt = {\Bf 0}$), In WhiCh THe SourcE is StatistIcally isotRopIc. however, in couPLiNg the fluxeS to Surrounding gAs eLementS, wE alSO must AccOunt fOr the fRAme motIon. BoostiNg BAck to tHE lab/simULAtIon fRame, the total ejectA VeLOcity entering aN elemeNt IS oF COurSe $\DElta m_{b}^{-1}\,\DeltA {\bF p}_{b} + {\bf V}_{A}$. This chANgE OF FraME has no effect oN the mass fluXEs, bUt it doEs ModIFy the mOmentUm ANd eNergy fluxes: To be Properly cOnservATive, we mUSt take: $$\bEgin{alIgnEd}
\lAbel{EQn:FlUx.mZ.fRAmeCOrR} \DeLTa m_{B}^{\prime} &\eqUiV \DElta m_{B}\ \ \, \ \ \  \DelTA M_{z,\,B}^{\priMe} \eQuiv \delta M_{Z,\,b} \\ 
\label{eqn:flUx.p.FramECorR} \DeltA {\bf p}_{b}^{\PrimE} &\eQuiv \DElta {\bf P}_{b} + \DelTa M_{b}\,{\bf v}_{a} \\ 
\label{eqn:fLux.e.Framecorr} \delTa e_{b}^{\pRiMe} &\equIV \Delta {e}_{b} + \fRac{1}{2\,\delta m_{b}}\,\Left( |\DelTA {\bf P}_{b}^{\PRIMe}|^{2} - |\delta {\bf p}_{b}|^{2} \right)\end{aLiGNEd}$$ Where the Prime (e.G. “$ \deLtA M_{b}^{\prime}$”) nOtAtiOn deNOTes thE lab FRaMe. Note thAt the eXTrA mOmentum AdDed to tHe NeiGhbOrs ($\suM_{B}\,\DelTa m_{b}\,{\bf V}_{a} = m_{\rm ej}\,{\bF v}_{a}$) is EXactly the momenTUm lost by the feEDbACK sOUrce $A$, by Virtue of its LosiNG $m_{\rm Ej}$ in MAsS.[^9]

ThESe fluXes arE sIMpLY added to each neighboR iN a fullY-consErvative manneR: $$\begin{aligNED}
\Label{eqn:Flux.M.CoUPling} m_{b}^{\rm new} &= m_{b} + \delta M_{b}^{\prime} \\ 
\labEL{eqn:flux.Z.coupLing
 good conservation require s the corr ected  te rms  a bove .

I mplicitly, weh avebeen working in the fr ame m ov i ng w i th  thefeedbac k  “ s o urc e”  ( ${\ bf x} _{a}= {\b f 0}$,${\bf v}_{ a}\e quiv d{\bf x } _{ a}/dt = {\ bf0}$), in whi chthe so ur cei s sta tis tical ly iso t ropic.  However,  i n  coupl i ng thef l ux es t o surrounding gas el e ments, we also  mustac c ou n t  fo r t he frame m ot ion.B oosting ba c k tot he lab/simula tion frame, the  total  e jec t a velo cityen t eri ng an eleme nt i s of cour se $\D e lta m_{ b }^{-1}\ ,\Delt a { \bf  p}_ { b}  +  {\ bf v}_ { a} $.T his  changeof  f ramehasn o e ffec t o n th e mas s fluxes, but  it  doe s  mo difythe m omen tu m and  energ y flu xe s: to be proper ly c onservati ve,  w e m us t tak e : $$\b egi n{a ligned} 
\label { eqn :f l u x .m z.framecorr} \Delt am _ {b }^{\prim e} &\e q ui v\ Delta m_ {b }\\ \, \  \ \ \De l ta  m_{Z,\, b}^{\p r im e}  \equiv  \ Deltam_ {Z, \,b } \\\lab el{eqn :flux.p. frame c orr} \Delta {\ b f p}_{b}^{\pr i me } &\ e quiv  \D elta {\bf p }_{b }  + \ Delt a  m _{b } \,{\b f v}_ {a }  \ \  
\label{eqn:flux.e .f rameco rr} \ Delta E_{b}^{ \prime} &\ e q u iv \Delt a {E } _{ b } + \frac{1}{2 \,\De lta m_{b}} \ ,\left(|\Del ta {\bfp}_{b}^{\ p r ime}|^{2 } -  |\ Del ta{ \ bf  p}_{b}|^{2}\ r ight )\ end{ali gne d}$$ wh ere  th e p rim e(e.g. “$\Delta m _{ b} ^{ \p rim e}$”) notation  d eno te s t he la b  frame . Not e th at  t h e e xtra mo m en t u m ad de dto t hene ighbo rs ( $ \su m_{b}\, \Delta m_ {b} \ ,{\b fv} _{a} =m_{\rm ej}\,{ \b f v}_{a}$)  i s e xactly t he momen tum lost by the feedbac k  source  $a $, by  vir tue of it s l osing$m_ { \rm ej }$ inmass. [^ 9]T hesef l ux esar e simply a d d edto ea ch  nei ghbor i n a fully-conserva t ive  manner: $$\b egi n{al i g ne d}\ la b el{ eq n :fl u x .m.coupling} m_ {b}^{\rm n ew }  & = m_{b} +\ Del ta  m_{b}^ {\prime } \\\label{ eqn:flux. z.couplin g
 good_conservation requires_the corrected terms above.

Implicitly,_we have_been_working in_the_frame moving with_the feedback “source”_(${\bf x}_{a}={\bf 0}$, ${\bf_v}_{a} \equiv d{\bf_x}_{a}/dt_= {\bf 0}$), in which the source is statistically isotropic. However, in coupling the_fluxes_to surrounding_gas_elements,_we also must account for_the frame motion. Boosting back_to the_lab/simulation frame, the total ejecta velocity entering an_element_is of course_$\Delta m_{b}^{-1}\,\Delta {\bf p}_{b} + {\bf v}_{a}$. This change_of frame has no effect on_the mass fluxes,_but_it_does modify the momentum_and energy fluxes: to be properly_conservative, we must take: $$\begin{aligned}
\label{eqn:flux.mz.framecorr} \Delta_m_{b}^{\prime} &\equiv \Delta m_{b}\ \ \, \_\ \  \Delta m_{Z,\,b}^{\prime} \equiv_\Delta m_{Z,\,b} \\ 
\label{eqn:flux.p.framecorr} \Delta_{\bf p}_{b}^{\prime}_&\equiv \Delta {\bf p}_{b} +_\Delta m_{b}\,{\bf v}_{a}_\\ 
\label{eqn:flux.e.framecorr}_\Delta E_{b}^{\prime} &\equiv_\Delta {E}_{b} + \frac{1}{2\,\Delta m_{b}}\,\left( |\Delta_{\bf p}_{b}^{\prime}|^{2} -_|\Delta {\bf p}_{b}|^{2} \right)\end{aligned}$$ where the_prime_(e.g. “$ \Delta m_{b}^{\prime}$”)_notation_denotes_the lab_frame. Note that_the_extra momentum_added_to the neighbors ($\sum_{b}\,\Delta m_{b}\,{\bf v}_{a}_=_m_{\rm ej}\,{\bf v}_{a}$) is exactly the momentum_lost by the feedback_source_$a$, by virtue of_its losing $m_{\rm ej}$ in_mass.[^9]

These fluxes are simply added to_each neighbor_in a_fully-conservative manner: $$\begin{aligned}
\label{eqn:flux.m.coupling} m_{b}^{\rm new} &= m_{b} + \Delta m_{b}^{\prime} \\_
\label{eqn:flux.z.coupling