 a value of $\approx 44$. Self-gravity is therefore not expected to lead to instabilities. As in M09, we assume that the initial mass of the disc is 1% of the mass of the star and we set the outer radius of the disc to 100 AU and the inner radius to 10 AU. Realistically, the disc should probably extend inwards towards the radius of the star, but it is computationally very expensive to simulate the disc on the orbital time scales at these small radii. During the simulation, particles in the disc will migrate inwards due to viscous evolution of the disc and are accreted onto the star, resulting in a disc that extends further inwards. This set-up allows us to postpone the expensive calculations towards later times in our simulations, thereby significantly decreasing the duration of our simulations.

We position the disc perpendicular to the flow (see Fig. \[fig:schematics\_eda1\]). This perfect alignment may not occur frequently in nature, but allows to discern the relevant physical processes more clearly. In future work we will investigate the influence of giving the disc an inclination with respect to the flow direction.

### Star {#sec:sink_eda1}

We assume that the mass of the star is concentrated in a single point, which has a mass of $0.4\,{\mbox{$\mathrm{M}_{\odot}$}}$. This corresponds to the typical mass expected in the early disc accretion scenario. The point mass is added as a collisionless particle to the SPH code and we treat it as a sink particle, i.e. it can accrete gas particles that fall within a certain (fixed) radius. The sink particle absorbs the mass and momentum carried by the gas particles it accretes.

We assume the radius of the sink particle is 5% of the Bondi-Hoyle accretion radius, $R_{\rm A}$, defined as [see also @Bondi52; @shima85]:

$$\label{eq:ra_eda1}
R_{\rm A} = \frac{2GM_{*}}{c_s^2 +v_{\rm ISM}^2}\\$$

in which $M_{*}$ is the mass of the star, $c_s$ is the sound speed and $v_{\rm ISM}$ the velocity of the ISM (with respect to the star). Using the values mentioned above, gives us an accretion radius of 1.8 AU, which is significantly smaller than in M09
a value of $ \approx 44$. Self - gravity is therefore not expected to lead to imbalance. As in M09, we wear that the initial mass of the disc is 1% of the mass of the headliner and we set the outer spoke of the disc to 100 AU and the inner radius to 10 AU. Realistically, the phonograph record should probably extend inwards towards the spoke of the star, but it is computationally very expensive to simulate the disc on the orbital time scale at these small radii. During the simulation, particles in the disc will migrate inwards due to syrupy evolution of the disc and are accreted onto the star topology, resulting in a disc that extends further inwards. This set - up let us to postpone the expensive calculations towards later times in our simulations, thereby significantly decrease the duration of our simulations. 

 We position the disc perpendicular to the flow (see Fig. \[fig: schematics\_eda1\ ]). This perfect alignment may not occur frequently in nature, but allows to discern the relevant physical processes more clearly. In future work we will investigate the influence of giving the disc an dip with respect to the flow steering. 

 # # # Star { # sec: sink_eda1 } 

 We wear that the mass of the star is concentrated in a single point, which has a mass of $ 0.4\,{\mbox{$\mathrm{M}_{\odot}$}}$. This correspond to the typical mass expected in the early disc accretion scenario. The point mass is add as a collisionless particle to the SPH code and we treat it as a sink particle, i.e. it can accrete accelerator particles that hang within a certain (fixed) radius. The sink particle steep the mass and momentum carried by the gas particle it accretes. 

 We assume the radius of the sink atom is 5% of the Bondi - Hoyle accretion radius, $ R_{\rm A}$, defined as [ see also @Bondi52; @shima85 ]: 

 $ $ \label{eq: ra_eda1 } 
 R_{\rm A } = \frac{2GM_{*}}{c_s^2 + v_{\rm ISM}^2}\\$$ 

 in which $ M_{*}$ is the mass of the ace, $ c_s$ is the sound focal ratio and $ v_{\rm ISM}$ the velocity of the ISM (with deference to the star). use the values mentioned above, give us an accretion radius of 1.8 AU, which is importantly smaller than in M09
 a galue of $\approx 44$. Self-grayity is thereforg bot ex'ected fo lead go instabilities. As in M09, we essune thqt the initial mass of the disc is 1% of rhe nqss of the star and we ssb the iuter radius on the disc do 100 AU and the ivnzr radius to 10 AU. Realistically, the dysc shoilf probably extgnd imrarda towards the radius of the star, bht it iv computationslly very expensive to simklatf the disc on the lrbital timg scwoes at these small radpn. During thg simulation, particles in the dirc wikl migrate ibwaggs due to vmscous evolution on the dhsc and are accreted pntm tye star, resulting in e disc that extends surther itwcrds. This set-up alloww ys to posdpond thd espxnsjve capcumations tosards later times in our simulstyins, thereby sjgnifisagtly decreasing the duration of our simllatjons.

We position the disx perpendicular to thg flow (see Fig. \[fig:schematics\_eda1\]). This perfect alignment may tot orcjr yvcquevrlj in nature, but allows to discern the relevane pnyxical processef more cleatlj. Og future work we will jnvestigate the invluence of guving the disv an inclination with respext to the fljq direction.

### Star {#szc:sink_eda1}

We cssume that the mass of the star ir cohcentrated ln a singmd point, which har a mdss of $0.4\,{\mbox{$\mathrm{M}_{\odot}$}}$. Thif correspinds to the typocal mwss expectfd in the early disc accgetiou scetario. The ooint mass is added as a collismpnless partivla tm the SPK code and we treat it as a sink karticle, n.e. it zan accretv gas parvicles that sall within a fertain (fixev) radius. Ehe wink particud absorbs the kass and momentum cqrried by the gas iartizmes it accretes.

Ct awsume the radiux ow tre smnk pwsticle is 5% ox thd Bumdi-Hohle accretiiv racius, $R_{\rm A}$, defined av [ses also @Bondi52; @shima85]:

$$\kanel{eq:ra_edq1}
R_{\rm A} = \srac{2GM_{*}}{c_s^2 +v_{\rm OSM}^2}\\$$

in which $M_{*}$ is tje mavs uf thr sjar, $c_s$ is the sound speed and $v_{\dm ISM}$ thf vclocity of thq ISN (with respeet to the star). Using the values mentionev above, gives us an acceetion radius of 1.8 AO, wmich is signmficanely smallar than in M09
a value of $\approx 44$. Self-gravity is expected lead to As in M09, mass the disc is of the mass the star and we set the radius of the disc to 100 AU and the inner radius to 10 Realistically, the disc should probably extend inwards towards the radius of the star, it computationally expensive simulate the disc on the orbital time scales at these small radii. During the simulation, particles the disc will migrate inwards due to viscous of the disc and accreted onto the star, resulting a that extends inwards. set-up us to postpone expensive calculations towards later times in our simulations, thereby significantly decreasing the duration of our simulations. We the disc the flow Fig. This alignment may not in nature, but allows to discern processes more clearly. In future work we will the influence giving the disc an inclination with to the flow direction. ### Star {#sec:sink_eda1} We that the mass of the star is concentrated in a single point, which has a $0.4\,{\mbox{$\mathrm{M}_{\odot}$}}$. This corresponds to typical mass expected the disc scenario. point mass added as a collisionless particle to the SPH code and we it as a sink particle, i.e. it can accrete gas fall a certain (fixed) The sink particle absorbs mass momentum carried by the it We of sink is 5% of the accretion radius, $R_{\rm A}$, defined [see also @Bondi52; @shima85]: +v_{\rm ISM}^2}\\$$ in which $M_{*}$ is the mass the star, $c_s$ is the sound speed $v_{\rm ISM}$ the velocity of the ISM (with respect to the star). the values gives us an accretion radius of 1.8 AU, is significantly smaller than M09
 a value of $\approx 44$. Self-gravity Is thereforE not eXpeCteD tO leaD to iNstabilities. As IN M09, we Assume that the initial maSs of tHe DIsc iS 1% Of The maSs of the STaR ANd wE sEt The OuTEr RadiuS of The disc To 100 AU and the InnEr Radius to 10 AU. ReALiStically, thE diSc should probAblY extenD iNwaRDs towArdS the rAdius oF The staR, but it is cOmPUtatioNAlly verY EXpEnsiVe to simulate the diSC oN The orbital time Scales At THeSE SmaLl rAdii. During ThE simuLAtion, paRTiCLES in THe disc will migRate inwards DUe tO viscoUs EvoLUtion oF the dIsC And Are accreted Onto The star, reSultinG In a disc THat exteNds furTheR inWardS. thIs Set-Up ALloWS uS to POstPone the eXpEnSive cAlcuLATIOns tOwaRds lAter tImes in our simuLatIons, THerEby siGnifiCantLy DecreAsing tHe durAtIon of our simulatIons.

we positioN thE dIsc PeRpendICular tO thE flOw (see FiG. \[fig:schEMatIcS\_EDA1\]). THis perfect alignmenT mAY NoT occur frEquentLY iN nATure, but aLlOws To diSCErn thE relEVaNt physicAl procESsEs More cleArLy. In fuTuRe wOrk We wilL InveStigatE the inflUence OF giving the disc AN inclination wITh RESpECt to The Flow directiOn.

### StAR {#sec:Sink_EDa1}

we aSSume tHat thE mASs OF the star is concentraTeD in a siNgle pOint, which has a Mass of $0.4\,{\mbox{$\MATHrm{M}_{\odot}$}}$. this COrREsponds to the tyPical Mass expectED in the eaRly diSc accretIon scenarIO. the point MasS is AddEd aS A CoLlisionless paRTIcle To The SPH cOde And we trEat It aS a sInk PaRticle, i.e. iT can accrEtE gAs PaRtiCles tHAt fall wiThIn a CeRtaIn (fixED) radiuS. The sInk pArTiCLe aBsorbs tHE mASS and MoMeNtum CarRiEd by tHe gaS ParTicles iT accretes.

we aSSume ThE rAdius of The sink particLe Is 5% of the BonDi-hoyLe accrETIon radiuS, $R_{\rm A}$, defined as [see also @BoNDi52; @shima85]:

$$\LabEl{eq:rA_eda1}
r_{\rm A} = \frac{2Gm_{*}}{c_s^2 +V_{\rm ISM}^2}\\$$

In wHIch $M_{*}$ is The masS of thE sTar, $C_S$ Is the SOUnD spEeD and $v_{\rm ISM}$ THE veLocitY oF the iSM (with Respect to the star). UsINg tHe values mentiOneD aboVE, GiVes US aN AccReTIon RADius of 1.8 AU, which is SignificanTlY SmAller than iN m09
 a value of $\approx 44$.Self-gravi ty is  th ere fo re n ot e xpected to lea d  toinstabilities. As in M 09, w ea ssum e  t hat t he init i al m ass  o fthe  d i sc  is 1 % o f the m ass of the  st ar  and we sett he  outer rad ius  of the disc  to  100 A Uand the i nne r rad ius to 10 AU.  Realisti ca l ly, th e  disc s h o ul d pr obably extend inw a rd s  towards the r adiusof th e sta r,but it isco mputa t ionally ve r y exp e nsive to simu late the di s c o n theor bit a l time  scal es atthese small  rad ii. Durin g thes imulati o n, part iclesinthe  dis c  w il l m ig r ate in war d s d ue to vi sc ou s evo luti o n o f th e d iscand a re accreted o nto  the sta r, re sulti ng i na dis c that  exte nd s further inwar ds.This set- upal low sus to postpo nethe  expens ive cal c ula ti o n s  t owards later times  i n ou r simula tions, th er e by signi fi can tlyd e creas ingt he  duratio n of o u rsi mulatio ns .

Wepo sit ion  thed iscperpen dicularto th e  flow (see Fig .  \[fig:schema t ic s \ _e d a1\] ).This perfec t al i gnme nt m a ynot occur  freq ue n tl y  in nature, but all ow s to d iscer n the relevan t physical p r ocessesmore cl e arly. In futur e wor k we willi nvestiga te th e influe nce of gi v i ng the d isc  an  in cli n a ti on with respe c t  toth e flowdir ection. 

# ##Sta r { #s ec:sink_e da1}

We  a ss um etha t the mass ofth e s ta r i s con c entrat ed in  a s in gl e  po int, wh i ch h as a  m as s of  $0 .4 \,{\m box{ $ \ma thrm{M} _{\odot}$ }}$ .  Thi sco rrespon ds to the typ ic al mass ex pe cte d in t h e  early d isc accretion scenario. The poi ntmassis a dded as a  co llisio nle s s part icle t o the  S PHc o de an d we  tr ea t it as as i nkparti cl e, i .e. itcan accrete gas pa r tic les that fall  wi thin a  c ert a in (fi xe d ) r a d ius. The sink p article ab so r bs  the massa ndmo mentumcarried  by t h e gas p articlesit accret es .

W e ass ume the ra dius ofthe sinkp artic l eis 5%  of  the B on di- Hoyle  accre t ion  radi us, $R _{ \rm A} $, de fi ned as [ see also @Bondi52; @shi ma85]: 

$$\ lab el{eq:ra_ eda 1 }
R _{\rm A}= \f rac{2GM_{* }}{ c_s ^2 +v _{\ r m ISM }^2} \ \$ $
i n whi ch $ M _{*}$ ist he  ma s s  o f the star, $ c _s$  is t hes ound s peed  and $v_{\rm ISM} $  the velocityof t h e  IS M ( w ithre spect to the s tar ). U sing the  v alues menti oned abo ve ,  give s us a n accr etion r a d iu s  of 1. 8 AU , w hich is s ign if i cantlysm al l er tha n in  M 09
 a_value of_$\approx 44$. Self-gravity is_therefore not_expected_to lead_to_instabilities. As in_M09, we assume_that the initial mass_of the disc_is_1% of the mass of the star and we set the outer radius of_the_disc to_100_AU_and the inner radius to_10 AU. Realistically, the disc_should probably_extend inwards towards the radius of the star,_but_it is computationally_very expensive to simulate the disc on the orbital_time scales at these small radii._During the simulation,_particles_in_the disc will migrate_inwards due to viscous evolution of_the disc and are accreted onto_the star, resulting in a disc that_extends further inwards. This set-up allows_us to postpone the expensive_calculations towards_later times in our simulations,_thereby significantly decreasing_the duration_of our simulations.

We_position the disc perpendicular to the_flow (see Fig._\[fig:schematics\_eda1\]). This perfect alignment may not_occur_frequently in nature,_but_allows_to discern_the relevant physical_processes_more clearly._In_future work we will investigate the_influence_of giving the disc an inclination with_respect to the flow_direction.

###_Star {#sec:sink_eda1}

We assume that_the mass of the star_is concentrated in a single point,_which has_a mass_of $0.4\,{\mbox{$\mathrm{M}_{\odot}$}}$. This corresponds to the typical mass expected in the_early disc accretion scenario. The point_mass is added as_a collisionless_particle_to the SPH_code_and we_treat it as a sink particle, i.e._it can_accrete gas particles that fall within_a certain (fixed) radius._The_sink particle absorbs the mass and_momentum carried by the gas particles_it accretes.

We assume the radius_of_the_sink particle is 5% of_the Bondi-Hoyle accretion radius, $R_{\rm A}$,_defined as [see_also @Bondi52; @shima85]:

$$\label{eq:ra_eda1}
R_{\rm A} = \frac{2GM_{*}}{c_s^2 +v_{\rm_ISM}^2}\\$$

in_which $M_{*}$ is the mass of_the_star, $c_s$ is the sound speed_and_$v_{\rm_ISM}$ the velocity of the_ISM (with respect to the star)._Using the values mentioned above, gives us an accretion_radius of 1.8_AU, which is significantly smaller_than_in_M09