# Hydrodynamic Modeling of the Interaction of Winds within a Collapsing Turbulent Gas Cloud

Authors: Guillermo Arreaga-García; Julio Saucedo-Morales
Journal: Advances in Astronomy
 (2015)

---

## Abstract

By using the particle-based code Gadget2, we follow the evolution of a gas giant molecular cloud, in which a set of gas particles representing the wind are created by a Monte Carlo scheme and suddenly move outwards from the cloud’s center. The particles representing the gas cloud initially have
a velocity according to a turbulent spectrum built in a Fourier space of 643 grid elements. The level of turbulence and the temperature of the cloud are both adjusted so that a gravitational collapse of the cloud is initially induced. All the winds are activated in a very early stage of evolution of the cloud. We consider only two kinds of winds, namely, one with spherical symmetry and the second one of a bipolar collimated jet. In order to assess the dynamical change in the cloud due to interactions with the winds, we show isovelocity and isodensity plots for all our simulations. We also report on the accretion centers detected at the last simulation time available for each model.


---

## Body

## 1. Introduction

Stars are born in large gas structures made of molecular hydrogen. According to [1] a cloud is a gas structure with mass and size around 1000–10000M⊙ and 2–15 pc, respectively.The physical process by which the molecular gas is transformed from a gas structure into stars is mainly gravitational collapse, whose main effect on the cloud is a reduction in size with the corresponding increase in density at various points of the cloud. These small overdensities within the larger gas structure are defined as cores, which have a typical mass and size around0.5–5M⊙ and 0.03–0.2pc, respectively, [1].The basic idea ofprompt fragmentation is that during gravitational collapse a molecular cloud may spontaneously break into several cores. At some point during the transformation process from clouds to cores and then to stars, the densest gas structures settle down in more or less dynamically stable gas objects called protostars. If the cores do not undergo further subfragmentation upon collapsing to higher densities, they will become protostars.To achieve a better understanding of the huge scale changes due to this transformation process, we may recall that the number density of a typical star is around1024 molecules per cm-3. Meanwhile, the number density of a typical typical cloud structure ranges around 103 molecules per cm-3 and that of a typical protostar ranges around 1014 molecules per cm-3.Many other phenomena occurring in the interstellar medium can have a great influence on the evolution of the cloud, among others: (i) the highly ionized gas ejected by supernovae explosions; (ii) the bipolar collimated winds ejected by massive protostars; (iii) rapidly expandingHII regions which hit the slower and less dense gas structures.For instance, in [2] a set of SPH simulations were conducted to study star formation triggered by an expanding HII region within a spherical gas structure of uniform density that was not globally collapsing. The expanding shock plays the role of a snowplow, which sweeps out the gas of the surrounding medium. The density of the swept out gas increases as a consequence of this agglomeration process, and a gravitational collapse may then be locally initiated. Small gas overdensities can be formed in this way, which may achieve the protostar stage. Because these two processes are complementary in forming protostars, this star formation scenario is known asthe collect and collapse model. The one-dimensional model of freely expanding wind into an uniform ambient medium was analitically and numerically solved some time ago; see [3] and references there in.In this paper, we investigate the change in the dynamical configuration of a typical turbulent cloud when a wind of particles is outwardly ejected from the central region of the cloud. The most important differences between the present paper and [2] are the turbulent nature of our cloud and its induced global tendency to gravitational collapse. Turbulence makes a big change in the spatial distribution as the cloud becomes filamentary and flocculent. The tendency to collapse favors the condensation of gas towards the central region of the cloud. Thus, a gas wind flowing through this collapsing and turbulent gas cloud is a very difficult hydrodynamical phenomenon that can only be solved by numerical simulations. It is shown that the presence of winds (i) keeps the collapsing nature of the turbulent cloud, but they cause a delay in the runway collapse; (ii) changes the number and location of accretion centers detected in the densest central region of the cloud.Furthermore, in order to study the effects of protostellar outflows on the turbulence of a star forming region, in [4, 5], magneto-hydrodynamics simulations (MHD) were conducted with a mesh-based code which implements the adaptive mesh refinement technique (AMR). As we use here the fully parallelized Gadget2 code which implements the SPH technique, this will make it possible to carry out an interesting comparison between the two computational techniques.


## 2. The Physical System

In this section we briefly describe the physics of the cloud and the winds which will be considered in the following sections.


### 2.1. The Cloud

We consider here a typical spherical cloud with a radiusR0=2pc and mass M0=1219M⊙. Initially, it has a radially uniform density distribution with an average density given by ρ0=2.4×10-21 gr cm-3, which is equivalent to a number density n0≈600molecules cm-3 for molecular hydrogen with molecular mass μ=4.0×10-24gr/mol. The size and mass of this cloud are chosen here to be typical in statistical sense, in accordance with [1].The free fall timetff is defined as the time needed for an external particle to reach the center of the cloud when gravity is the only force pulling the particle. In this idealized gravitational collapse, we have(1)tff≈3π32Gρ0,where G is Newton’s universal gravitational constant. For our cloud, we have tff=1.3×106years.For a spherical cloud, the approximate total gravitational potential energy isEgrav≈-3/5GM02/R0. The average total thermal energy Etherm (kinetic plus potential interaction terms of the molecules) is Etherm≈3/2NkT=3/2M0c02, where k is the Boltzmann constant, T is the equilibrium temperature, N is the total number of molecules in the gas, and c0 is the speed of sound; see Section 3.5 for a more precise definition.The kinetic energyEkin can be estimated by M0vav2/2, where vav is the average translational velocity of the cloud. In order to have both energies of the same order of magnitude, Ekin≈Egrav, the gas elements of the cloud must attain average velocities within the range(2)vavc0≈3-3.5or vav≈1.6km/sec, for a speed of sound given by(3)c0=0.54km/secso that the corresponding temperature associated with the cloud is T≈15K.It is possible to define the crossing time by means of(4)tcr≈R0c0=3.56×106yrwhich sets a time scale for a sound wave to travel across the cloud. To make the crossing time comparable in magnitude to the free fall time of (1), the front wave must have velocities around vreq/c0≈2.6 or vreq≈1.45 km/sec, which are velocities a little bit slower than the ones estimated above; see (2). In this paper, we will treat propagation velocities of gas particles ranging around 2-3 Mach.


### 2.2. The Wind

The dynamical characteristics of the wind strongly depend on its type of source. All stars eject winds of the first kind, which are driven by stellar radiation. For instance, in cool stars, like the ones observed in the AGB (asymtotic giant branch) of the Galaxy, the winds cause a mass loss in the range10-8–10-4M⊙/yr whereas the terminal wind velocities are around 10–45km/sec. In OB stars, the mass loss ranges over 10-6–10-4M⊙/yr and the terminal wind velocities can go up to thousands of km/sec.Supernovae dump around104 joules of thermal and kinetic energy into the interstellar medium. But there are several types of supernovas, so that the mass losses and terminal velocities are very different; see [6]. For example, for a supernova whose progenitor was a He star, the mass loss and terminal velocities are within the ranges 10-7–10-4M⊙/yr and 100–1000km/sec, respectively. When the progenitor was a RSG star, their values range over 10-5–10-4M⊙/yr and 10–40km/sec, respectively.It seems that all protostars eject highly collimated jets of gas during their formation process by gravitational collapse. The origin of these jets is still unclear, but it is very likely that the accretion disk and magnetic field around protostars play a crucial role in determining the velocities and degree of collimation of the jets. For the molecular winds associated with protostarts of Class 0 and Class 1, the characteristic velocities are around 20 km/sec. However, for optical jets of highly ionized gas, the typical jet velocities are a few hundred km/sec; see [7] and the references therein.


## 2.1. The Cloud

We consider here a typical spherical cloud with a radiusR0=2pc and mass M0=1219M⊙. Initially, it has a radially uniform density distribution with an average density given by ρ0=2.4×10-21 gr cm-3, which is equivalent to a number density n0≈600molecules cm-3 for molecular hydrogen with molecular mass μ=4.0×10-24gr/mol. The size and mass of this cloud are chosen here to be typical in statistical sense, in accordance with [1].The free fall timetff is defined as the time needed for an external particle to reach the center of the cloud when gravity is the only force pulling the particle. In this idealized gravitational collapse, we have(1)tff≈3π32Gρ0,where G is Newton’s universal gravitational constant. For our cloud, we have tff=1.3×106years.For a spherical cloud, the approximate total gravitational potential energy isEgrav≈-3/5GM02/R0. The average total thermal energy Etherm (kinetic plus potential interaction terms of the molecules) is Etherm≈3/2NkT=3/2M0c02, where k is the Boltzmann constant, T is the equilibrium temperature, N is the total number of molecules in the gas, and c0 is the speed of sound; see Section 3.5 for a more precise definition.The kinetic energyEkin can be estimated by M0vav2/2, where vav is the average translational velocity of the cloud. In order to have both energies of the same order of magnitude, Ekin≈Egrav, the gas elements of the cloud must attain average velocities within the range(2)vavc0≈3-3.5or vav≈1.6km/sec, for a speed of sound given by(3)c0=0.54km/secso that the corresponding temperature associated with the cloud is T≈15K.It is possible to define the crossing time by means of(4)tcr≈R0c0=3.56×106yrwhich sets a time scale for a sound wave to travel across the cloud. To make the crossing time comparable in magnitude to the free fall time of (1), the front wave must have velocities around vreq/c0≈2.6 or vreq≈1.45 km/sec, which are velocities a little bit slower than the ones estimated above; see (2). In this paper, we will treat propagation velocities of gas particles ranging around 2-3 Mach.


## 2.2. The Wind

The dynamical characteristics of the wind strongly depend on its type of source. All stars eject winds of the first kind, which are driven by stellar radiation. For instance, in cool stars, like the ones observed in the AGB (asymtotic giant branch) of the Galaxy, the winds cause a mass loss in the range10-8–10-4M⊙/yr whereas the terminal wind velocities are around 10–45km/sec. In OB stars, the mass loss ranges over 10-6–10-4M⊙/yr and the terminal wind velocities can go up to thousands of km/sec.Supernovae dump around104 joules of thermal and kinetic energy into the interstellar medium. But there are several types of supernovas, so that the mass losses and terminal velocities are very different; see [6]. For example, for a supernova whose progenitor was a He star, the mass loss and terminal velocities are within the ranges 10-7–10-4M⊙/yr and 100–1000km/sec, respectively. When the progenitor was a RSG star, their values range over 10-5–10-4M⊙/yr and 10–40km/sec, respectively.It seems that all protostars eject highly collimated jets of gas during their formation process by gravitational collapse. The origin of these jets is still unclear, but it is very likely that the accretion disk and magnetic field around protostars play a crucial role in determining the velocities and degree of collimation of the jets. For the molecular winds associated with protostarts of Class 0 and Class 1, the characteristic velocities are around 20 km/sec. However, for optical jets of highly ionized gas, the typical jet velocities are a few hundred km/sec; see [7] and the references therein.


## 3. The Computational Methods and the Models

In this section we briefly describe the way we set up the physical system outlined in Section2 in computational terms as well as the models we consider in this paper.


### 3.1. The Initial Configuration of Particles

We setN=10 million SPH particles for representing the gas cloud. By means of a rectangular mesh we make the partition of the simulation volume in small elements each with a volume ΔxΔyΔz; at the center of each volume we place a particle (the ith, say), with a mass determined by its location according to the density profile being considered, that is, mi=ρ(xi,yi,zi)∗ΔxΔyΔz with i=1,…,N. Next, we displace each particle from its location by a distance on the order of Δx/4.0 in a random spatial direction.As stated earlier, in this paper we only consider a uniform density cloud, for whichρ(xi,yi,zi)≡ρ0, for all the simulations (see Section 2.1). Therefore, all the cloud particles have the same mass.


### 3.2. The Initial Turbulent Velocity of Particles

To generate the turbulent velocity spectrum for the cloud particleswe follow a procedure based on the papers [8, 9]. We set a second mesh Nx=Ny=Nz=64 with the size of each element given by δx=R0/Nx, δy=R0/Ny, and δz=R0/Nz. In Fourier space, the partition is δKx=1.0/Nx×δx, δKy=1.0/Ny×δy, and δKz=1.0/Nz×δz. Each Fourier mode has the components Kx=iKxδKx, Ky=iKyδKy, and Kz=iKzδKz, where the indices iKx, iKy, and iKz take values in [-Nx/2,Nx/2], [-Ny/2,Ny/2], and [-Nz/2,Nz/2], respectively. The wave number magnitude is K=Kx2+Ky2+Kz2, and so Kmin=0 and Kmax=3Nx/2R0. The Fourier wave can be equally described by a wave length λ=2π/K, and then we see that K≈1/R0 and λ≈R0.The components of the particle velocity are(5)v→=Σ-Kmax⁡Kmax⁡K-n-2/2×KzCKysin⁡K→·r→+ΦKy-KyCKzsin⁡K→·r→+ΦKzforvx,-KxCKzsin⁡K→·r→+ΦKz+KzCKxsin⁡K→·r→+ΦKxforvy,-KxCKysin⁡K→·r→+ΦKy+KyCKxsin⁡K→·r→+ΦKxforvz,where the spectral index n was fixed at n=-1 and thus we have v2≈K-3. The vector C→K whose components are denoted by CKx,CKy,CKz takes values obeying a Rayleigh distribution. The wave phase vector, Φ→K, given by ΦKx,ΦKy,ΦKz takes random values on the interval [0,2π]. The components of the vector C→ are calculated by means of C=σ×-2.0×log1.0-u, where u is a random number in (0,1). σ is a fixed parameter with value σ=1.0.


### 3.3. The Set-Up of the Particle Wind

We generate the initial wind particles by a traditional Monte Carlo scheme, in which particles are initially located randomly in the central volume of the cloud. Letu, v, and w be random uniform variables taking real values within the interval [0,1]. According to the fundamental probability conservation law for a system with spherical symmetry, we must have 4πr2ρrdrsinθdθdϕ=dudvdw. By means of integration we obtain that the spherical coordinates (r,θ,ϕ) of the particles are related to the uniform random variables by the following equations:(6)u=4πMg∫0rρgr′r′2dr′,v=1-cos⁡θ2=12∫0θsin⁡θ′dθ′,w=ϕ2π=12π∫0ϕdϕ′,where Mg is the total mass contained in the central cloud region of radius Rg and ρg(r) its mass density. It must be noticed that the first relation should be integrated numerically to obtain r once u has taken an allowed uniform random value.In order to assign a velocity vector to each wind particle, let us consider the equation of mass conservation for a set of particles moving radially outwards; that is,(7)M˙=4πr2ρrvr.In order to determine the wind velocities according to (7), where the radius r ranges in the interval r=(0,Rg), the mass loss M˙ and the mass density ρg will be considered as free parameters in the simulations; see Section 3.8. As the velocity magnitude diverges for particles around r≈0, we must then set up a cut velocity value, vcut. Therefore, the velocity vector for the wind particle is given by v→=ve^r, where v is the speed of the particle obtained by (7) and e^r is the unit vector pointing radially outwards.Of course, there are other possibilities, which will be considered elsewhere: one is to fix the radial densityρ(r) and/or the velocity profile v(r) in order to obtain the mass loss M˙ as a result. For example, for modeling an expanding HII region, the authors of [2] proposed a complicated velocity function, which gives a constant expansion velocity at the last stages of time evolution, so that the average velocity of the shocked shell considered by [2] is vcc/c0≈5.6 or vcc≈3.7 km/sec.


### 3.4. Resolution and Thermodynamical Considerations

Following [10, 11], in order to avoid artificial fragmentation, the SPH code must fulfil certain resolution criteria, imposed on the Jeans wavelength λJ, which is given by(8)λJ=πc2Gρ,where c is the instantaneous speed of sound and ρ is the local density. To obtain a more useful form for a particle based code, the Jeans wavelength λJ is transformed into the Jeans mass given by(9)MJ≡43πρλJ23=π5/26c3G3ρ.In this paper, the values of the density and speed of sound are updated according to the following equation of state:(10)p=c02ρ1+ρρcritγ-1,as proposed by [12], where γ≡5/3, and for the critical density we assume the value ρcrit=5.0×10-14gcm-3.For the turbulent cloud under consideration, we havemr≈MJ/(2Nneigh)≈7.47×1033g, where we take Nneigh=40.In this paper, the mass of an SPH particle ismp=1.98×1029 g, so that mp/mr=2.5×10-4 and therefore the Jeans resolution requirement is very easily satisfied.


### 3.5. Initial Energies

Following [13], the dynamical properties of the initial distribution of gas are usually characterized by α, the ratio of the thermal energy to the gravitational energy, and β, the ratio of the rotational energy to the gravitational energy.In a particle based code, we approximate the thermal energy of the cloud by calculating the sum over all theN particles described in Section 3.1; that is,(11)Etherm=∑i=1N32Piρmiρi,where Pi and ρi are the pressure and density associated with particle i by means of the equation of state given in (10). In a similar way, the approximate potential energy is(12)Epot=∑i=1N12miΦi,where Φi is the gravitational potential of particle i.For all the simulations considered in this paper, the values of the speed of sound c0 (see (3)) and the level of turbulence are chosen so that the energy ratios have the numerical values(13)α≡EthermEgrav=0.3,β≡EkinEgrav=1.0.


### 3.6. The Evolution Code

We carry out the time evolution of the initial distribution of particles with the fully parallelGadget2 code, which is described in detail by [14]. Gadget2 is based on the tree-PM method for computing the gravitational forces and on the standard SPH method for solving the Euler equations of hydrodynamics. Gadget2 incorporates the following standard features: (i) each particle i has its own smoothing length hi; (ii) the particles are also allowed to have individual gravitational softening lengths ϵi, whose values are adjusted so that for every time step ϵihi is of order unity. Gadget2 fixes the value of ϵi for each time-step using the minimum value of the smoothing length of all particles; that is, if hmin=min(hi) for i=1,2,…N, then ϵi=hmin.TheGadget2 code has an implementation of a Monaghan-Balsara form for the artificial viscosity; see [15, 16]. The strength of the viscosity is regulated by setting the parameters αν=0.75 and βν=3/2×αν; see (14) in [14]. Here we fix the Courant factor to 0.1.


### 3.7. The Detection of Accretion Centers

Let us now describe the code implemented to detect accretion centers in the models. Any gas particle with density higher thanρs is a candidate to be an accretion center. We localize all the candidate particles for a given snapshot at time t.We then test the separation between candidate particles: if there is one candidate with no other candidates closer than10×racc, then this particle is identified as an accretion center at time t. We define racc as the neighbor radius for an accretion center, given by racc=15×hmin, where hmin is the minimum value of the smoothing length of all particles; that is, if hmin=min(hi) with hi being the smoothing length of each particle i, for i=1,2,…N in this way racc determines a set of particles which are within the sphere having this radius and whose center is the accretion center itself. All those particles will contribute to the mass and momentum of the accretion center.We carry out the detection of accretion centers for all the models, using a fixed value ofρs given by ρs=100.0×ρ0, as was done by [5].The code described here is a first step towards the full sink-particle technique implementation already in use worldwide. Our code is still incomplete as we have not yet included all the particle tests suggested in [17, 18].


### 3.8. The Models

The initial configuration of particles has many parameters and their values determine the models to consider as follows. To sum up the meaning of these parameters in order of importance, we make a list in Table1, without any reference to some specific physical system. In the first line, tstart determines the time when the wind appears within the turbulent cloud. In the second and third lines, M˙ and ρg determine the velocity by (7). Next, the radius Rg determines the initial spatial extension of the wind particles. Ng (Mg) is the total number (the total mass) of wind particles within Rg. vcut is the peak wind velocity allowed at the initial simulation time to avoid divergences in the velocity given by (7). Finally, we need to define two angles in order to characterize the collimated jets: the polar angular aperture θg and the azimuthal angular aperture ϕg, centered around the Y-axis with ϕ=π/2, and the equator with θ=π/2, respectively.Table 1
The parameters of the initial configuration code.
Parameter
Name or concept
t
start
The time at which the wind starts
M
˙
The ejected mass rate
ρ
g
The average density of wind particles
R
g
The outer boundary wind radius
N
g or Mg
The total number (or total mass) of the wind particles
v
cut
The cutting velocity
θ
g
The polar angular aperture
ϕ
g
The azimuthal angular apertureIn this paper, we fixtstart/tff=0.05 and vcut=100km/sec for all the simulations while θg=10° and θg=10° for the collimated jet simulations. The choice of a wind density value ρg must not be so different to the cloud uniform value ρ0, to avoid discontinuity problems between the gas cloud and the gas wind when Gadget2 calculates the tree. Besides, it should be noticed that values of Rg have been chosen to represent a core (with typical values of size and mass) within the cloud.In Table2 we present the values used to define the models calculated below as well as some results described in the coming sections. In this Table and throughout the paper, we use the following labels: the label “T” to indicate the turbulent gas cloud when the winds appear within the cloud, the label “TW” to refer to the spherically symmetric case, and the label “TWC” for the bipolar collimated jet.Table 2
Parameters and some outcomes for the models considered.
Model
M
˙  [M⨀/yr]
R
g
/
R
0
N
g
ρ
g
/
ρ
0
M
g
[
M
⨀
]
M
∞
[
M
⨀
]
v
∞
/
c
s
t
max
⁡
/
t
ff
ρ
max
⁡ [gr/cm3]
T
—
—
10,000,000
1
1219
333
10
0.95
2.75
×
1
0
-
12
TW1
1.0
×
1
0
-
4
0.1
56336
1
6.84
359
11
1.06
1.37
×
1
0
-
16
TW2
1.0
×
1
0
-
7
0.05
187769
0.01
0.33
333
10
0.96
9.98
×
1
0
-
16
TWC1
1.0
×
1
0
-
3
0.05
37170
1
4.51
422
72
1.11
2.67
×
1
0
-
14
TWC2
1.0
×
1
0
-
7
0.05
82320
0.01
0.21
333
10
0.96
1.45
×
1
0
-
15
TWC3
1.0
×
1
0
-
4
0.1
119425
1
14
358
8
1.0
1.84
×
1
0
-
16The first six columns of Table2 were already defined in the preceding paragraphs. In the seventh column of Table 2, we show the total mass M∞ swept out by the winds outside R0 by the end of the simulation. In column 8 of Table 2, we show the velocities attained v∞ by those particles located far outside the initial cloud. In columns 9 and 10 of Table 2, we show the maximum evolution time and peak density reached in the simulations.


## 3.1. The Initial Configuration of Particles

We setN=10 million SPH particles for representing the gas cloud. By means of a rectangular mesh we make the partition of the simulation volume in small elements each with a volume ΔxΔyΔz; at the center of each volume we place a particle (the ith, say), with a mass determined by its location according to the density profile being considered, that is, mi=ρ(xi,yi,zi)∗ΔxΔyΔz with i=1,…,N. Next, we displace each particle from its location by a distance on the order of Δx/4.0 in a random spatial direction.As stated earlier, in this paper we only consider a uniform density cloud, for whichρ(xi,yi,zi)≡ρ0, for all the simulations (see Section 2.1). Therefore, all the cloud particles have the same mass.


## 3.2. The Initial Turbulent Velocity of Particles

To generate the turbulent velocity spectrum for the cloud particleswe follow a procedure based on the papers [8, 9]. We set a second mesh Nx=Ny=Nz=64 with the size of each element given by δx=R0/Nx, δy=R0/Ny, and δz=R0/Nz. In Fourier space, the partition is δKx=1.0/Nx×δx, δKy=1.0/Ny×δy, and δKz=1.0/Nz×δz. Each Fourier mode has the components Kx=iKxδKx, Ky=iKyδKy, and Kz=iKzδKz, where the indices iKx, iKy, and iKz take values in [-Nx/2,Nx/2], [-Ny/2,Ny/2], and [-Nz/2,Nz/2], respectively. The wave number magnitude is K=Kx2+Ky2+Kz2, and so Kmin=0 and Kmax=3Nx/2R0. The Fourier wave can be equally described by a wave length λ=2π/K, and then we see that K≈1/R0 and λ≈R0.The components of the particle velocity are(5)v→=Σ-Kmax⁡Kmax⁡K-n-2/2×KzCKysin⁡K→·r→+ΦKy-KyCKzsin⁡K→·r→+ΦKzforvx,-KxCKzsin⁡K→·r→+ΦKz+KzCKxsin⁡K→·r→+ΦKxforvy,-KxCKysin⁡K→·r→+ΦKy+KyCKxsin⁡K→·r→+ΦKxforvz,where the spectral index n was fixed at n=-1 and thus we have v2≈K-3. The vector C→K whose components are denoted by CKx,CKy,CKz takes values obeying a Rayleigh distribution. The wave phase vector, Φ→K, given by ΦKx,ΦKy,ΦKz takes random values on the interval [0,2π]. The components of the vector C→ are calculated by means of C=σ×-2.0×log1.0-u, where u is a random number in (0,1). σ is a fixed parameter with value σ=1.0.


## 3.3. The Set-Up of the Particle Wind

We generate the initial wind particles by a traditional Monte Carlo scheme, in which particles are initially located randomly in the central volume of the cloud. Letu, v, and w be random uniform variables taking real values within the interval [0,1]. According to the fundamental probability conservation law for a system with spherical symmetry, we must have 4πr2ρrdrsinθdθdϕ=dudvdw. By means of integration we obtain that the spherical coordinates (r,θ,ϕ) of the particles are related to the uniform random variables by the following equations:(6)u=4πMg∫0rρgr′r′2dr′,v=1-cos⁡θ2=12∫0θsin⁡θ′dθ′,w=ϕ2π=12π∫0ϕdϕ′,where Mg is the total mass contained in the central cloud region of radius Rg and ρg(r) its mass density. It must be noticed that the first relation should be integrated numerically to obtain r once u has taken an allowed uniform random value.In order to assign a velocity vector to each wind particle, let us consider the equation of mass conservation for a set of particles moving radially outwards; that is,(7)M˙=4πr2ρrvr.In order to determine the wind velocities according to (7), where the radius r ranges in the interval r=(0,Rg), the mass loss M˙ and the mass density ρg will be considered as free parameters in the simulations; see Section 3.8. As the velocity magnitude diverges for particles around r≈0, we must then set up a cut velocity value, vcut. Therefore, the velocity vector for the wind particle is given by v→=ve^r, where v is the speed of the particle obtained by (7) and e^r is the unit vector pointing radially outwards.Of course, there are other possibilities, which will be considered elsewhere: one is to fix the radial densityρ(r) and/or the velocity profile v(r) in order to obtain the mass loss M˙ as a result. For example, for modeling an expanding HII region, the authors of [2] proposed a complicated velocity function, which gives a constant expansion velocity at the last stages of time evolution, so that the average velocity of the shocked shell considered by [2] is vcc/c0≈5.6 or vcc≈3.7 km/sec.


## 3.4. Resolution and Thermodynamical Considerations

Following [10, 11], in order to avoid artificial fragmentation, the SPH code must fulfil certain resolution criteria, imposed on the Jeans wavelength λJ, which is given by(8)λJ=πc2Gρ,where c is the instantaneous speed of sound and ρ is the local density. To obtain a more useful form for a particle based code, the Jeans wavelength λJ is transformed into the Jeans mass given by(9)MJ≡43πρλJ23=π5/26c3G3ρ.In this paper, the values of the density and speed of sound are updated according to the following equation of state:(10)p=c02ρ1+ρρcritγ-1,as proposed by [12], where γ≡5/3, and for the critical density we assume the value ρcrit=5.0×10-14gcm-3.For the turbulent cloud under consideration, we havemr≈MJ/(2Nneigh)≈7.47×1033g, where we take Nneigh=40.In this paper, the mass of an SPH particle ismp=1.98×1029 g, so that mp/mr=2.5×10-4 and therefore the Jeans resolution requirement is very easily satisfied.


## 3.5. Initial Energies

Following [13], the dynamical properties of the initial distribution of gas are usually characterized by α, the ratio of the thermal energy to the gravitational energy, and β, the ratio of the rotational energy to the gravitational energy.In a particle based code, we approximate the thermal energy of the cloud by calculating the sum over all theN particles described in Section 3.1; that is,(11)Etherm=∑i=1N32Piρmiρi,where Pi and ρi are the pressure and density associated with particle i by means of the equation of state given in (10). In a similar way, the approximate potential energy is(12)Epot=∑i=1N12miΦi,where Φi is the gravitational potential of particle i.For all the simulations considered in this paper, the values of the speed of sound c0 (see (3)) and the level of turbulence are chosen so that the energy ratios have the numerical values(13)α≡EthermEgrav=0.3,β≡EkinEgrav=1.0.


## 3.6. The Evolution Code

We carry out the time evolution of the initial distribution of particles with the fully parallelGadget2 code, which is described in detail by [14]. Gadget2 is based on the tree-PM method for computing the gravitational forces and on the standard SPH method for solving the Euler equations of hydrodynamics. Gadget2 incorporates the following standard features: (i) each particle i has its own smoothing length hi; (ii) the particles are also allowed to have individual gravitational softening lengths ϵi, whose values are adjusted so that for every time step ϵihi is of order unity. Gadget2 fixes the value of ϵi for each time-step using the minimum value of the smoothing length of all particles; that is, if hmin=min(hi) for i=1,2,…N, then ϵi=hmin.TheGadget2 code has an implementation of a Monaghan-Balsara form for the artificial viscosity; see [15, 16]. The strength of the viscosity is regulated by setting the parameters αν=0.75 and βν=3/2×αν; see (14) in [14]. Here we fix the Courant factor to 0.1.


## 3.7. The Detection of Accretion Centers

Let us now describe the code implemented to detect accretion centers in the models. Any gas particle with density higher thanρs is a candidate to be an accretion center. We localize all the candidate particles for a given snapshot at time t.We then test the separation between candidate particles: if there is one candidate with no other candidates closer than10×racc, then this particle is identified as an accretion center at time t. We define racc as the neighbor radius for an accretion center, given by racc=15×hmin, where hmin is the minimum value of the smoothing length of all particles; that is, if hmin=min(hi) with hi being the smoothing length of each particle i, for i=1,2,…N in this way racc determines a set of particles which are within the sphere having this radius and whose center is the accretion center itself. All those particles will contribute to the mass and momentum of the accretion center.We carry out the detection of accretion centers for all the models, using a fixed value ofρs given by ρs=100.0×ρ0, as was done by [5].The code described here is a first step towards the full sink-particle technique implementation already in use worldwide. Our code is still incomplete as we have not yet included all the particle tests suggested in [17, 18].


## 3.8. The Models

The initial configuration of particles has many parameters and their values determine the models to consider as follows. To sum up the meaning of these parameters in order of importance, we make a list in Table1, without any reference to some specific physical system. In the first line, tstart determines the time when the wind appears within the turbulent cloud. In the second and third lines, M˙ and ρg determine the velocity by (7). Next, the radius Rg determines the initial spatial extension of the wind particles. Ng (Mg) is the total number (the total mass) of wind particles within Rg. vcut is the peak wind velocity allowed at the initial simulation time to avoid divergences in the velocity given by (7). Finally, we need to define two angles in order to characterize the collimated jets: the polar angular aperture θg and the azimuthal angular aperture ϕg, centered around the Y-axis with ϕ=π/2, and the equator with θ=π/2, respectively.Table 1
The parameters of the initial configuration code.
Parameter
Name or concept
t
start
The time at which the wind starts
M
˙
The ejected mass rate
ρ
g
The average density of wind particles
R
g
The outer boundary wind radius
N
g or Mg
The total number (or total mass) of the wind particles
v
cut
The cutting velocity
θ
g
The polar angular aperture
ϕ
g
The azimuthal angular apertureIn this paper, we fixtstart/tff=0.05 and vcut=100km/sec for all the simulations while θg=10° and θg=10° for the collimated jet simulations. The choice of a wind density value ρg must not be so different to the cloud uniform value ρ0, to avoid discontinuity problems between the gas cloud and the gas wind when Gadget2 calculates the tree. Besides, it should be noticed that values of Rg have been chosen to represent a core (with typical values of size and mass) within the cloud.In Table2 we present the values used to define the models calculated below as well as some results described in the coming sections. In this Table and throughout the paper, we use the following labels: the label “T” to indicate the turbulent gas cloud when the winds appear within the cloud, the label “TW” to refer to the spherically symmetric case, and the label “TWC” for the bipolar collimated jet.Table 2
Parameters and some outcomes for the models considered.
Model
M
˙  [M⨀/yr]
R
g
/
R
0
N
g
ρ
g
/
ρ
0
M
g
[
M
⨀
]
M
∞
[
M
⨀
]
v
∞
/
c
s
t
max
⁡
/
t
ff
ρ
max
⁡ [gr/cm3]
T
—
—
10,000,000
1
1219
333
10
0.95
2.75
×
1
0
-
12
TW1
1.0
×
1
0
-
4
0.1
56336
1
6.84
359
11
1.06
1.37
×
1
0
-
16
TW2
1.0
×
1
0
-
7
0.05
187769
0.01
0.33
333
10
0.96
9.98
×
1
0
-
16
TWC1
1.0
×
1
0
-
3
0.05
37170
1
4.51
422
72
1.11
2.67
×
1
0
-
14
TWC2
1.0
×
1
0
-
7
0.05
82320
0.01
0.21
333
10
0.96
1.45
×
1
0
-
15
TWC3
1.0
×
1
0
-
4
0.1
119425
1
14
358
8
1.0
1.84
×
1
0
-
16The first six columns of Table2 were already defined in the preceding paragraphs. In the seventh column of Table 2, we show the total mass M∞ swept out by the winds outside R0 by the end of the simulation. In column 8 of Table 2, we show the velocities attained v∞ by those particles located far outside the initial cloud. In columns 9 and 10 of Table 2, we show the maximum evolution time and peak density reached in the simulations.


## 4. Results

To present the results of our simulations, we consider a slice of particles around the equatorial plane of the cloud; with these particles (around 10000) we make plots containing two complementary panels: one to show colored regions of isodensity and the other one to show the velocity field of the particles. Later, we present plots of the velocity distributions and the radial velocity profile, for which we use all the simulation particles. Finally, in order to make a comparison of the changes produced in the dynamical state of the central region of the cloud, we present2d and 3d plots with the accretion centers detected by locating the densest particles of each simulation.


### 4.1. The Evolution of the Turbulent Cloud

Because of the random nature of the initial velocity distribution, a huge number of particle collisions occur across the cloud. In the first two panels of Figure1 one can already notice one of the main characteristics of turbulence: the appearance of a filamentary and flocculent structure across the cloud, a spatial structure similar to “Swiss cheese.”Figure 1
Iso-density plots for a thin slice of matter parallel toXY plane illustrating the evolution of the turbulent cloud, model T.The particle collisions dissipate very quickly much of the kinetic energy of the cloud, so that the cloud goes to a physical state similar to a free fall collapse; then the gravitational force pulls the cloud mass to a global collapse towards its central region, as can be seen in the last panels of Figure1.As one can see in the top panel of Figure2, very few particles can attain velocities much higher than those average velocities provided by turbulence alone. But almost all the simulation particles have velocities v/cs<4. In the middle panel of Figure 2 we show the radial profile of the velocity, which clearly indicates that almost all particles have the same average velocity across the cloud, which is around 2 Mach. At the time t/tff=0.05, the fully flocculent cloud reaches a the peak in its density curve, as we show in the bottom panel of Figure 2.Figure 2
For the initial snapshot of each simulation we show (a) the velocity distribution, (b) the velocity radial profile, and (c) the peak density time evolution.
(a)
(b)
(c)We emphasize now a very important fact occurring at the outer regions of the turbulent cloud. Near the cloud boundary, the interpolating kernel of the Gadget2 code introduces small errors in the calculation of the particle density that makes the boundary density lower than it should be. So to speak, the mass of each SPH particle weighs less for particles located near the boundary. There are methods for correcting this density deficiency problem; see, for example, [19, 20]. Unfortunately, such corrections have not been incorporated so far into the public Gadget2 code. Furthermore, the turbulent cloud is not in hydrodynamic equilibrium and there is no external pressure acting upon the cloud, and then the outermost particles have nonequilibrated thermal pressure. All these features act together in such a way that the outer cloud particles expand outwards.So, we have to keep in mind this expansion effect for the problem at hand, as we shall quantify the mass of the cloud swept out by the winds; see, for instance, Figure3 and the next section.Figure 3
For the last snapshot available in each simulation, we plot (a) the velocity distribution, (b) the velocity radial profile, and (c) a zoom-in of the middle panel.
(a)
(b)
(c)


### 4.2. The Effects of Winds in the Evolution of the Turbulent Cloud

Let us again consider the top and middle panels of Figure2, which show that the wind particles near the cloud’s center can attain velocities much higher than v/cs≈4. It must be kept in mind that these curves are smoothed; thus they appear to steep even for r much greater than Rg.Because the wind particles are located at the center of the cloud by a Monte Carlo method, the particle density there significantly increases at timet/tff=0.05, when the winds are activated, as can be seen in the bottom panel of Figure 2. It should also be noticed by looking at this panel that the global collapse of the cloud does not change even when the winds particles are present, as the peak density curve of each run goes to higher values by the end of the simulation time.We emphasize that for all the wind models the gravitational collapse of the cloud is slightly slower as their peak the curves begin to slope upwards a little later the middle up to the final stages of their evolution time. The reason to claim this delay is that the peak density values reached by the wind models are always lower than the turbulent cloud compared at a similar evolution time. For instance, let us mention that the turbulent cloud alone needs0.95tff to achieve a central peak density of 2.75×10-12gr/cm3 while the models TW1 and TWC3 have peak densities several orders of magnitude smaller at a evolution time slightly greather than tff. The model TWC1 needs 1.11tff to reach a peak density of 2.67×10-14gr/cm3, as can be seen in the last two columns of Table 2.The isodensity and velocity plots for all the simulations are shown in Figures4, 5, 6, 7, and 8. It is recommended to make zooming in the panels of these figures in order to notice in all the model the occurrence of a common transformation process: from a very disordered velocity distribution to a very clear tendency to condensate toward the central cloud region. There are still some differences in the final morphology and in the accretion centers detected in the central region of the cloud.Figure 4
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TW1.”Figure 5
Isodensity and velocity plotsfor a thin slice of matter parallel to XY plane for the model “TW2.”Figure 6
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TWC1.”Figure 7
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TWC2.”Figure 8
Isodensity and velocity distribution plots for a thin slice of matter parallel toXY plane for the model “TWC3.”Let us first consider the models with spherically symmetric winds, TW1 and TW2. The wind particles for model TW1 push the cloud particles outwards. Thus, by the timet/tff=0.36, we see a void created in the central region of the cloud. However, gravity and viscosity act together in such a way that the particles quickly fill the void, as can be seen in the last panel of Figure 4. Instead, the wind particles for model TW2 diffuse through the cloud without creating a void. As there is no much more mass swept out than that of the turbulent cloud alone, it is possible that all the wind particles reach the outer cloud boundary, but they do not push gas particles outwards.Let us now consider the bipolar collimated jet models TWC1 and TWC2 illustrated in Figures6 and 7. Essentially, we see either collision (for model TWC1) or diffusion (for model TWC2) of the wind particles taking place in the directions in which the collimated jets are initially emitted.For model TWC1, the lack of wind particles in the perpendicular directions of the emitted jet produces a slightly elongated void, as can be seen in the second panel of Figure6. As expected, the cloud mass located in the perpendicular directions to the emitted jets is strongly attracted by the gravitational force of cloud mass located in the other side. When the central cloud region is filled again with mass, it adopts a slightly elongated condensation along the direction of the emitted jets. With the purpose of emphasizing this effect, we include another collimated model in which we increase the initial mass, the model TWC3, which is illustrated in Figure 8.We emphasize that the fraction of gas swept out by the winds is quite significant: from nothing to dozens of solar masses moved even far beyond the cloud radiusR0, as can be seen in the seventh column of Table 2. We also report in the eighth column the expansion velocities reached during the time we follow the simulations. The most evolved simulation we have is for model TWC1. This fact alone explains why it has the highest values reported in Table 2 for the the masses M∞ and v∞.As expected, all the wind models show a strong tendency to form accretion centers in the central region of the cloud. It must be noticed that the great differences appear when we compare the number and location of the accretion centers obtained for the turbulent cloud and for the wind models, as can be seen in Figures9 and 10. In these plots, it is possible to see the spatial distribution of the accretion centers: every point represents a SPH particle, and then we show all the particles captured by the formed accretion centers.Figure 9
For the last snapshot available for each simulation, we show the accretion centers detected with the densityρs=100ρ0 and projected on the XY plane.Figure 10
For the last snapshot available for each simulation, we show the accretion centers detected with the densityρs=100ρ0.
(a)
(b)
(c)
(d)
(e)
(f)


## 4.1. The Evolution of the Turbulent Cloud

Because of the random nature of the initial velocity distribution, a huge number of particle collisions occur across the cloud. In the first two panels of Figure1 one can already notice one of the main characteristics of turbulence: the appearance of a filamentary and flocculent structure across the cloud, a spatial structure similar to “Swiss cheese.”Figure 1
Iso-density plots for a thin slice of matter parallel toXY plane illustrating the evolution of the turbulent cloud, model T.The particle collisions dissipate very quickly much of the kinetic energy of the cloud, so that the cloud goes to a physical state similar to a free fall collapse; then the gravitational force pulls the cloud mass to a global collapse towards its central region, as can be seen in the last panels of Figure1.As one can see in the top panel of Figure2, very few particles can attain velocities much higher than those average velocities provided by turbulence alone. But almost all the simulation particles have velocities v/cs<4. In the middle panel of Figure 2 we show the radial profile of the velocity, which clearly indicates that almost all particles have the same average velocity across the cloud, which is around 2 Mach. At the time t/tff=0.05, the fully flocculent cloud reaches a the peak in its density curve, as we show in the bottom panel of Figure 2.Figure 2
For the initial snapshot of each simulation we show (a) the velocity distribution, (b) the velocity radial profile, and (c) the peak density time evolution.
(a)
(b)
(c)We emphasize now a very important fact occurring at the outer regions of the turbulent cloud. Near the cloud boundary, the interpolating kernel of the Gadget2 code introduces small errors in the calculation of the particle density that makes the boundary density lower than it should be. So to speak, the mass of each SPH particle weighs less for particles located near the boundary. There are methods for correcting this density deficiency problem; see, for example, [19, 20]. Unfortunately, such corrections have not been incorporated so far into the public Gadget2 code. Furthermore, the turbulent cloud is not in hydrodynamic equilibrium and there is no external pressure acting upon the cloud, and then the outermost particles have nonequilibrated thermal pressure. All these features act together in such a way that the outer cloud particles expand outwards.So, we have to keep in mind this expansion effect for the problem at hand, as we shall quantify the mass of the cloud swept out by the winds; see, for instance, Figure3 and the next section.Figure 3
For the last snapshot available in each simulation, we plot (a) the velocity distribution, (b) the velocity radial profile, and (c) a zoom-in of the middle panel.
(a)
(b)
(c)


## 4.2. The Effects of Winds in the Evolution of the Turbulent Cloud

Let us again consider the top and middle panels of Figure2, which show that the wind particles near the cloud’s center can attain velocities much higher than v/cs≈4. It must be kept in mind that these curves are smoothed; thus they appear to steep even for r much greater than Rg.Because the wind particles are located at the center of the cloud by a Monte Carlo method, the particle density there significantly increases at timet/tff=0.05, when the winds are activated, as can be seen in the bottom panel of Figure 2. It should also be noticed by looking at this panel that the global collapse of the cloud does not change even when the winds particles are present, as the peak density curve of each run goes to higher values by the end of the simulation time.We emphasize that for all the wind models the gravitational collapse of the cloud is slightly slower as their peak the curves begin to slope upwards a little later the middle up to the final stages of their evolution time. The reason to claim this delay is that the peak density values reached by the wind models are always lower than the turbulent cloud compared at a similar evolution time. For instance, let us mention that the turbulent cloud alone needs0.95tff to achieve a central peak density of 2.75×10-12gr/cm3 while the models TW1 and TWC3 have peak densities several orders of magnitude smaller at a evolution time slightly greather than tff. The model TWC1 needs 1.11tff to reach a peak density of 2.67×10-14gr/cm3, as can be seen in the last two columns of Table 2.The isodensity and velocity plots for all the simulations are shown in Figures4, 5, 6, 7, and 8. It is recommended to make zooming in the panels of these figures in order to notice in all the model the occurrence of a common transformation process: from a very disordered velocity distribution to a very clear tendency to condensate toward the central cloud region. There are still some differences in the final morphology and in the accretion centers detected in the central region of the cloud.Figure 4
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TW1.”Figure 5
Isodensity and velocity plotsfor a thin slice of matter parallel to XY plane for the model “TW2.”Figure 6
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TWC1.”Figure 7
Isodensity and velocity plots for a thin slice of matter parallel toXY plane for the model “TWC2.”Figure 8
Isodensity and velocity distribution plots for a thin slice of matter parallel toXY plane for the model “TWC3.”Let us first consider the models with spherically symmetric winds, TW1 and TW2. The wind particles for model TW1 push the cloud particles outwards. Thus, by the timet/tff=0.36, we see a void created in the central region of the cloud. However, gravity and viscosity act together in such a way that the particles quickly fill the void, as can be seen in the last panel of Figure 4. Instead, the wind particles for model TW2 diffuse through the cloud without creating a void. As there is no much more mass swept out than that of the turbulent cloud alone, it is possible that all the wind particles reach the outer cloud boundary, but they do not push gas particles outwards.Let us now consider the bipolar collimated jet models TWC1 and TWC2 illustrated in Figures6 and 7. Essentially, we see either collision (for model TWC1) or diffusion (for model TWC2) of the wind particles taking place in the directions in which the collimated jets are initially emitted.For model TWC1, the lack of wind particles in the perpendicular directions of the emitted jet produces a slightly elongated void, as can be seen in the second panel of Figure6. As expected, the cloud mass located in the perpendicular directions to the emitted jets is strongly attracted by the gravitational force of cloud mass located in the other side. When the central cloud region is filled again with mass, it adopts a slightly elongated condensation along the direction of the emitted jets. With the purpose of emphasizing this effect, we include another collimated model in which we increase the initial mass, the model TWC3, which is illustrated in Figure 8.We emphasize that the fraction of gas swept out by the winds is quite significant: from nothing to dozens of solar masses moved even far beyond the cloud radiusR0, as can be seen in the seventh column of Table 2. We also report in the eighth column the expansion velocities reached during the time we follow the simulations. The most evolved simulation we have is for model TWC1. This fact alone explains why it has the highest values reported in Table 2 for the the masses M∞ and v∞.As expected, all the wind models show a strong tendency to form accretion centers in the central region of the cloud. It must be noticed that the great differences appear when we compare the number and location of the accretion centers obtained for the turbulent cloud and for the wind models, as can be seen in Figures9 and 10. In these plots, it is possible to see the spatial distribution of the accretion centers: every point represents a SPH particle, and then we show all the particles captured by the formed accretion centers.Figure 9
For the last snapshot available for each simulation, we show the accretion centers detected with the densityρs=100ρ0 and projected on the XY plane.Figure 10
For the last snapshot available for each simulation, we show the accretion centers detected with the densityρs=100ρ0.
(a)
(b)
(c)
(d)
(e)
(f)


## 5. Discussion

In this paper we carry out a fully 3D set of numerical hydrodynamical simulations within the framework of theSPH technique, aimed at following the interaction of winds with a turbulent cloud. We comment now about three technical limitations we ougth to keep in mind when looking at the results reported here.Firstly, the Monte Carlo method can generate couples of very close particles. These couples should be removed because the Gadget2 code gets into troubles when the construction of the tree is iniciatted. We romove close particles by running a code that searches for those particles which are closer than an allowed minimun distance. Thus more particles than really needed are removed by this procedure. The net effect is that there is no control on the total number of wind particles to be evolved in a simulation. This explains why the number of particles entering in the wind models are so different, as can be seen in Table2.Secondly, the Gadget2 code used here was written in the so-calleddensity-entropy formulation, for which the adaptive smoothing length of each particle hi is determined by creating a kernel volume containing a constant mass, in which a number of neighbour particles is kept roughly constant for each particle i. For this reason, we can not have a wind density so different from the cloud density as the smoothing lengths hi may become too large for those less clustered wind particles; see [21].Thirdly, it must be taken into consideration that (i) we have used a low density thresholdρs for a particle to become an accretion center and (ii) there is no restriction on the minimum number of particles that must be captured for a particle to really become an accretion center: so we have seen that there are a large number of accretion centers having just a very few captured particles. For these two reasons we must expect overproduction of accretion centers in some of our simulations. However, as we have applied the same code to all the models, a comparison is possible. Finally, we recall that accretion centers are made of gas particles which can collide or merge. In view of this, the future of the accretion centers is still unclear and requires taking the calculations further in time.


## 6. Concluding Remarks

A star formation scenario based only on the collapse of turbulent gas structures provides a very highly efficient transformation of the gas into protostars. However, this point is in contradiction with observations.Because of this, another scenario must be considered, or at least a theoretical complement to the turbulent model is needed. As we have shown in this paper, winds must be considered an additional ingredient to complement the turbulent model with the hope that these models TW and TWC can alleviate the problem mentioned above as they make a delay in the runaway collapse of the cloud.We have obtained a physical system showing simultaneously in-fall and out-flow motions of particles. We emphasize that these simultaneous motions are observed in clusters like NGC 1333 and NGC 2264; see [22, 23], respectively.Another interesting objet is IRAS 16547-4247. It has a very dense single central region surrounded by a cloud of about1300M⊙. A couple of wind jets have been detected recently in this star forming region and they are belived to be generated by a massive protostar; see [24]. The model studied here can be modified in order to take into account the new findings.This paper must be considered a first step towards a more complete study of winds within a collapsing turbulent cloud. A forthcoming publication will treat a different mass and density for the wind particles and the different directions of the emitted jets to model system like the IRAS source 16547-4247.



---

*Source: 196304-2015-05-20.xml*