 first order $\alpha_{\rm SPH}$ term, while preventing particle interpenetration in the strong accretion shock. We have tested lowering the $\beta_{\rm SPH}$ value below the adopted value and found that it can cause numerical artefacts.

Following @artymowicz96, we calculate the time scale on which the disc undergoes significant viscous evolution, $\tau_{\nu}=Re\,\Omega^{-1}$, assuming the Reynolds number, $Re$, and angular frequency, $\Omega$, are given by: $$\label{eq:reynolds_eda1}
Re = \frac{(r/\langle h \rangle)^2}{\alpha_{\nu}}\\
\Omega = \sqrt{\frac{GM_{*}}{r^3}}\\$$ with $r$ the radial distance to the star and $\langle h \rangle$ the (resolution-dependent) average smoothing length in the disc at that radius. For a simulation of 16.000 disc particles, and assuming the corresponding $\alpha_{\nu}$, this gives us a viscous time scale, $\tau_{\nu}$, of roughly 10.000 years at 10 AU. We can also calculate the physical viscous time scale at 10 AU by replacing $\langle h \rangle$ with the scale height of the disc, $H=\frac{c_s}{\Omega}$, at that radius. This leads to a time scale of $3\times10^5$ years. In order to remain safely in the regime where the simulation is not dominated by viscous evolution of the disc, we evolve the whole set-up for 2500 years and start the inflow at $t=0$. The numerical viscosity in the disc changes during the simulations with Gadget2 and we will discuss the numerical effects in Sec. \[sec:convergence\_eda1\] and \[sec:consistency\_eda1\].

Distinguishing the disc from the ISM {#sec:disc_recognition_eda1}
------------------------------------

\

\

To determine how much ISM has been swept-up by the disc at any moment in time, we need to differentiate between the disc and the ISM flow. This is not straightforward since there is no clear separation between the outer edge of the disc and the flow. In order to distinguish between the disc and the ISM flow, we use a clump finding algorithm to identify different groups in the parameter space of the angular speed ($v_{\theta}$), the density ($\rho$) and the speed along the axis of the cylinder($v_
first order $ \alpha_{\rm SPH}$ term, while preventing atom interpenetration in the impregnable accretion shock. We have tested turn down the $ \beta_{\rm SPH}$ value below the adopt value and found that it can cause numeral artifact. 

 Following @artymowicz96, we calculate the time scale on which the disc undergo significant viscous evolution, $ \tau_{\nu}=Re\,\Omega^{-1}$, assuming the Reynolds numeral, $ Re$, and angular frequency, $ \Omega$, are given by: $ $ \label{eq: reynolds_eda1 } 
 Re = \frac{(r/\langle h \rangle)^2}{\alpha_{\nu}}\\ 
 \Omega = \sqrt{\frac{GM_{*}}{r^3}}\\$$ with $ r$ the radial distance to the star and $ \langle h \rangle$ the (resolution - dependent) modal smoothing length in the disc at that radius. For a pretense of 16.000 disc particles, and assuming the corresponding $ \alpha_{\nu}$, this pass us a viscous time scale, $ \tau_{\nu}$, of roughly 10.000 years at 10 AU. We can besides calculate the physical viscous time scale at 10 AU by replacing $ \langle planck's constant \rangle$ with the scale height of the disc, $ H=\frac{c_s}{\Omega}$, at that radius. This leads to a time scale of $ 3\times10 ^ 5 $ years. In order to remain safely in the government where the simulation is not dominated by gluey development of the disk, we evolve the whole hardening - up for 2500 years and start the inflow at $ t=0$. The numerical viscosity in the disc change during the simulations with Gadget2 and we will discuss the numerical effects in Sec. \[sec: convergence\_eda1\ ] and \[sec: consistency\_eda1\ ]. 

 Distinguishing the disc from the ISM { # sec: disc_recognition_eda1 } 
 ------------------------------------ 

 \ 

 \ 

 To decide how much ISM has been swept - up by the disc at any moment in time, we necessitate to differentiate between the phonograph record and the ISM flow. This is not straightforward since there is no clear separation between the out edge of the disc and the flow. In decree to distinguish between the disc and the ISM flow, we use a ball finding algorithm to identify unlike groups in the parameter space of the angular speed ($ v_{\theta}$), the concentration ($ \rho$) and the speed along the axis of the cylinder($v _
 figst order $\alpha_{\rm SPH}$ tevm, while preventnbg parvicle ihterpenegration in the strong accretmon whock. We have tested lowerivg the $\bena_{\rm SPH}$ calut below the adoptxs value and fkmnd tkav it can cause kumerical astefacts.

Followhne @crtymowicz96, we calculate the time scaje on wnifh the disc unqergpqs sjgnificant viscous evolution, $\tau_{\nu}=De\,\Omega^{-1}$, assuming the Reynolds number, $Re$, and anhulag frequency, $\Omega$, wre given bi: $$\lafwl{eq:reynolds_dda1}
Re = \frab{(x/\langle h \rzngle)^2}{\alpha_{\nu}}\\
\Omega = \sqrt{\frac{GM_{*}}{r^3}}\\$$ dith $x$ the radiao eishdnce to the star and $\langle m \rangla$ the (rrsolution-depencenv) avwrage smoothing lengti in the disc at thaj radius. Fmr a simulation of 16.000 dusc pdrtiwles, qnd asaukihg the cocresponding $\alpha_{\nu}$, thus gives us a viscois nome scale, $\tah_{\nu}$, of rjughly 10.000 years at 10 AU. We can also calculdte the physical viscous tume scale at 10 AU by rgplacing $\lwngle h \rangle$ with the scale height of the disc, $V=\frac{r_s}{\Umeyq}$, at gyah radius. This leads to a time scale of $3\times10^5$ rsats. In order to rcmain safely in thr gebyme where the simulcfikn is not dominatef by vifcous evolutiog of the disc, we evolve the whooe set-up for 2500 years and start tke inflow at $t=0$. Tne nukerical viscosity in thz disc changes duging the akmulations with Eadnet2 and we will discuss the ntmerical xffeccs in Sez. \[sev:convewgence\_eda1\] wnd \[scw:consistency\_eda1\].

Dishingunshinc the disc from the ISM {#sec:disc_recognitioi_xda1}
------------------------------------

\

\

To determime hof much IFM hax been swept-u[ by the disc ct any mjment in time, wv need to differentiwte between tvg disc and thx ISM flor. Thus iw not sgfaightforward xince thege is bo clear separatiok betdsen the outer ebye of the disc anc tfe sllw. Ig order to divtineuirn betdeen tht dlsc and the ISM flow, we use a cmump finding algorotmm to idebtify disferent groupx in the parameter spacx of tie angilat speed ($v_{\theta}$), the density ($\rho$) znd the soeeq along the wxis of the cyliuder($v_
first order $\alpha_{\rm SPH}$ term, while preventing in strong accretion We have tested below adopted value and that it can numerical artefacts. Following @artymowicz96, we calculate time scale on which the disc undergoes significant viscous evolution, $\tau_{\nu}=Re\,\Omega^{-1}$, assuming the number, $Re$, and angular frequency, $\Omega$, are given by: $$\label{eq:reynolds_eda1} Re = \frac{(r/\langle \rangle)^2}{\alpha_{\nu}}\\ = with the radial distance to the star and $\langle h \rangle$ the (resolution-dependent) average smoothing length in disc at that radius. For a simulation of disc particles, and assuming corresponding $\alpha_{\nu}$, this gives us viscous scale, $\tau_{\nu}$, roughly years 10 AU. We also calculate the physical viscous time scale at 10 AU by replacing $\langle h \rangle$ with the height of $H=\frac{c_s}{\Omega}$, at radius. leads a time scale years. In order to remain safely where the simulation is not dominated by viscous of the we evolve the whole set-up for years and start the inflow at $t=0$. The viscosity in the disc changes during the simulations with Gadget2 and we will discuss the in Sec. \[sec:convergence\_eda1\] and Distinguishing the disc the {#sec:disc_recognition_eda1} \ To determine much ISM has been swept-up by the disc at any moment time, we need to differentiate between the disc and the This not straightforward since is no clear separation the edge of the disc flow. order the and ISM flow, we use clump finding algorithm to identify groups in the parameter ($v_{\theta}$), the density ($\rho$) and the speed along axis of the cylinder($v_
 first order $\alpha_{\rm SPH}$ term, wHile prevenTing pArtIclE iNterPeneTration in the stROng aCcretion shock. We have tesTed loWeRIng tHE $\bEta_{\rm sPH}$ valuE BeLOW thE aDoPteD vALuE and fOunD that it Can cause nuMerIcAl artefacts.

FOLlOwing @artymOwiCz96, we calculatE thE time sCaLe oN Which The Disc uNdergoES signiFicant visCoUS evoluTIon, $\tau_{\nU}=rE\,\OMega^{-1}$, Assuming the ReynolDS nUMber, $Re$, and angulAr freqUeNCy, $\oMEga$, Are Given by: $$\labEl{Eq:reyNOlds_eda1}
rE = \fRAC{(R/\laNGle h \rangle)^2}{\alpHa_{\nu}}\\
\Omega = \sqRT{\frAc{GM_{*}}{r^3}}\\$$ wItH $r$ tHE radiaL distAnCE to The star and $\lAnglE h \rangle$ tHe (resoLUtion-dePEndent) aVerage SmoOthIng lENgTh In tHe DIsc AT tHat RAdiUs. For a siMuLaTion oF 16.000 disC PARTiclEs, aNd asSuminG the corresponDinG $\alpHA_{\nu}$, This gIves uS a viScOus tiMe scalE, $\tau_{\nU}$, oF roughly 10.000 years at 10 aU. We Can also caLcuLaTe tHe PhysiCAl viscOus TimE scale aT 10 AU by rePLacInG $\LANgLe h \rangle$ with the scAlE HEiGht of the Disc, $H=\fRAc{C_s}{\oMega}$, at thAt RadIus. THIS leadS to a TImE scale of $3\Times10^5$ yEArS. IN order tO rEmain sAfEly In tHe regIMe whEre the SimulatiOn is nOT dominated by viSCous evolution OF tHE DiSC, we eVolVe the whole sEt-up FOr 2500 yeArs aND sTarT The inFlow aT $t=0$. tHe NUmerical viscosity in ThE disc cHangeS during the simUlations wiTH gAdget2 and We wiLL dIScuss the numeriCal efFects in Sec. \[SEc:converGence\_Eda1\] and \[seC:consisteNCY\_eda1\].

DistIngUisHinG thE DIsC from the ISM {#seC:DIsc_rEcOgnitioN_edA1}
------------------------------------

\

\

To deteRmiNe hOw mUch iSm has been sWept-up by ThE dIsC aT anY momeNT in time, wE nEed To DifFerenTIate beTween The dIsC aND thE ISM floW. thIS Is noT sTrAighTfoRwArd siNce tHEre Is no cleAr separatIon BEtweEn ThE outer eDge of the disc aNd The flow. In oRdEr tO distiNGUish betwEen the disc and the ISM flow, WE use a clUmp FindiNg alGorithm to IdeNtify dIffERent grOups in The paRaMetER Space OF ThE anGuLar speed ($v_{\tHETa}$), tHe denSiTy ($\rhO$) and the Speed along the axis oF The Cylinder($v_
 first order $\alpha_{\rmSPH}$ term , whi lepre ve ntin g pa rticle interpe n etra tion in the strong acc retio ns hock .  W e hav e teste d  l o w eri ng  t he$\ b et a_{\r m S PH}$ va lue belowthe  a dopted value an d found th atit can cause  nu merica lart e facts .
Follo wing @ a rtymow icz96, we  c a lculat e  the ti m e  s cale  on which the dis c  u n dergoes signif icantvi s co u s  ev olu tion, $\ta u_ {\nu} = Re\,\Om e ga ^ { - 1}$ ,  assuming the  Reynolds n u mbe r, $Re $,  an d  angul ar fr eq u enc y, $\Omega$ , ar e given b y: $$\ l abel{eq : reynold s_eda1 }
R e =  \fr a c{ (r /\l an g leh  \ ran g le) ^2}{\alp ha _{ \nu}} \\
\ O m e g a =\sq rt{\ frac{ GM_{*}}{r^3}} \\$ $ wi t h $ r$ th e rad ialdi stanc e to t he st ar  and $\langle h  \ra ngle$ the  (r es olu ti on-de p endent ) a ver age smo othingl eng th i n  t he disc at that ra di u s .For a si mulati o nof 16.000 d is c p arti c l es, a nd a s su ming the  corre s po nd ing $\a lp ha_{\n u} $,thi s giv e s us  a vis cous tim e sca l e, $\tau_{\nu} $ , of roughly1 0. 0 0 0y ears  at  10 AU. Wecana lsocalc u la tet he ph ysica lv is c ous time scale at 1 0AU byrepla cing $\langle  h \rangle $ w ith thescal e  h e ight of the di sc, $ H=\frac{c_ s }{\Omega }$, a t that r adius. Th i s  leads t o a  ti mesca l e  o f $3\times10^ 5 $  yea rs . In or der  to rem ain  sa fel y i nthe regim e whereth esi mu lat ion i s  not dom in ate dbyvisco u s evol ution  ofth ed isc , we ev o lv e thewh ol e se t-u pfor 2 500y ear s and s tart theinf l ow a t$t =0$. Th e numerical v is cosity inth e d isc ch a n ges duri ng the simulations with Gadget2  an d wewill  discussthe  numer ica l  effec ts inSec.\[ sec : c onver g e nc e\_ ed a1\] and \ [ s ec: consi st ency \_eda1\ ].

Distinguishing the  disc from th e I SM { # s ec :di s c_ r eco gn i tio n _ eda1}
--------- ---------- -- - -- ---------- - -
\
\

Todetermi ne ho w  much I SM has be en swept- up  byt h e d isc at any  momentin time,w e nee d  t o dif fer entiat ebet weenthe di s c a nd th e ISMfl ow. Th is is  n ot strai ghtforward since thereis noclear  se parationbet w een  the oute r ed ge of thedis c a nd th e f l ow. I n or d er  to disti ngui s h between th e d i s cand the ISM f l ow,  we u sea  clump  fin ding algorithm to identify diffe rent g rou psi n th eparameter spac e o ft h e angula rspeed ($v_{ \theta}$ ), the d ensity  ($\rh o$) and t he speedalon g t he axis o f t he cylinde r( $v _ 
 first_order $\alpha_{\rm_SPH}$ term, while preventing_particle interpenetration_in_the strong_accretion_shock. We have_tested lowering the_$\beta_{\rm SPH}$ value below_the adopted value_and_found that it can cause numerical artefacts.

Following @artymowicz96, we calculate the time scale on_which_the disc_undergoes_significant_viscous evolution, $\tau_{\nu}=Re\,\Omega^{-1}$, assuming the_Reynolds number, $Re$, and angular_frequency, $\Omega$,_are given by: $$\label{eq:reynolds_eda1}
Re = \frac{(r/\langle h \rangle)^2}{\alpha_{\nu}}\\
\Omega_=_\sqrt{\frac{GM_{*}}{r^3}}\\$$ with $r$_the radial distance to the star and $\langle h_\rangle$ the (resolution-dependent) average smoothing length_in the disc_at_that_radius. For a simulation_of 16.000 disc particles, and assuming_the corresponding $\alpha_{\nu}$, this gives us_a viscous time scale, $\tau_{\nu}$, of roughly_10.000 years at 10 AU. We_can also calculate the physical_viscous time_scale at 10 AU by_replacing $\langle h_\rangle$ with_the scale height_of the disc, $H=\frac{c_s}{\Omega}$, at that_radius. This leads_to a time scale of $3\times10^5$_years._In order to_remain_safely_in the_regime where the_simulation_is not_dominated_by viscous evolution of the disc,_we_evolve the whole set-up for 2500 years_and start the inflow_at_$t=0$. The numerical viscosity_in the disc changes during_the simulations with Gadget2 and we_will discuss_the numerical_effects in Sec. \[sec:convergence\_eda1\] and \[sec:consistency\_eda1\].

Distinguishing the disc from the ISM_{#sec:disc_recognition_eda1}
------------------------------------

\

\

To determine how much ISM has_been swept-up by the_disc at_any_moment in time,_we_need to_differentiate between the disc and the ISM_flow. This_is not straightforward since there is_no clear separation between_the_outer edge of the disc and_the flow. In order to distinguish_between the disc and the_ISM_flow,_we use a clump finding_algorithm to identify different groups in_the parameter space_of the angular speed ($v_{\theta}$), the density_($\rho$)_and the speed along the axis_of_the cylinder($v_