Search Archive | News and Event | Author Guideline |

Sorry.

You are not permitted to access the full text of articles.

If you have any questions about permissions,

please contact the Society.

죄송합니다.

회원님은 논문 이용 권한이 없습니다.

권한 관련 문의는 학회로 부탁 드립니다.

Journal of the Korean Astronomical Society -
Vol. 53 ,
No. 3

[ Article ] | |

Journal of the Korean Astronomical Society - Vol. 53, No. 2, pp.49-57 | |

Abbreviation: JKAS | |

ISSN: 1225-4614 (Print) 2288-890X (Online) | |

Print publication date 30 Apr 2020 | |

Received 10 Oct 2019 Accepted 18 Mar 2020 | |

DOI: https://doi.org/10.5303/JKAS.2020.53.2.49 | |

DECAY OF TURBULENCE IN FLUIDS WITH POLYTROPIC EQUATIONS OF STATE | |

JEONGHOON LIM ^{1} ; JUNGYEON CHO^{1}^{, 2}
| |

1Department of Astronomy and Space Science, Chungnam National University, 99, Daehak-ro, Yuseong-gu, Daejeon 34134, Republic of Korea (jhlim0918@o.cnu.ac.kr)(jcho@cnu.ac.kr) | |

2Korea Astronomy and Space Science Institute, 776, Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Republic of Korea | |

Correspondence to : J. Cho | |

Published under Creative Commons license CC BY-SA 4.0 | |

Funding Information ▼ |

Abstract

We present numerical simulations of decaying hydrodynamic turbulence initially driven by solenoidal (divergence-free) and compressive (curl-free) drivings. Most previous numerical studies for decaying turbulence assume an isothermal equation of state (EOS). Here we use a polytropic EOS, *P* ∝ *ρ*^{γ} , with polytropic exponent *γ* ranging from 0.7 to 5/3. We mainly aim at determining the effects of *γ* and driving schemes on the decay law of turbulence energy, *E* ∝ *t*^{−α}. We additionally study probability density function (PDF) of gas density and skewness of the distribution in polytropic turbulence driven by compressive driving. Our findings are as follows. First of all, we find that even if *γ* does not strongly change the decay law, the driving schemes weakly change the relation; in our all simulations, turbulence decays with *α* ≈ 1, but compressive driving yields smaller *α* than solenoidal driving at the same sonic Mach number. Second, we calculate compressive and solenoidal velocity components separately and compare their decay rates in turbulence initially driven by compressive driving. We find that the former decays much faster so that it ends up having a smaller fraction than the latter. Third, the density PDF of compressively driven turbulence with *γ* > 1 deviates from log-normal distribution: it has a power-law tail at low density as in the case of solenoidally driven turbulence. However, as it decays, the density PDF becomes approximately log-normal. We discuss why decay rates of compressive and solenoidal velocity components are different in compressively driven turbulence and astrophysical implication of our findings.

Keywords: ISM: general, hydrodynamics, turbulence |

1. INTRODUCTION

Supersonic turbulence in the interstellar medium (ISM) is a well-known phenomenon and plays an essential role in star formation processes (Larson 1981; Padoan & Nordlund 2002; Mac Low & Klessen 2004). Given that driving mechanisms of astrophysical turbulence are usually intermittent in both space and time, it is natural for turbulence to decay. Earlier studies showed that non-driven turbulence decays quickly in approximately one large-eddy turnover time (for hydrodynamic turbulence, see e.g., Lesieur 2008; for magnetohydrodynamic turbulence, see Mac Low et al. 1998; Stone et al. 1998), which is consistent with the fact that an energy cascade occurs within one large-scale eddy turnover time even in the case of strongly magnetized turbulence (Goldreich & Sridhar 1995).

It has been analytically suggested that turbulence energy decays with a power-law of the form *E* ∝ *t*^{−α} (see e.g., Chapter 7 of Lesieur 2008). Results from previous numerical studies of turbulence have converged to a value of *α* of approximately unity, and it does not strongly depend on the degree of magnetization and compressibility (Mac Low et al. 1998; Stone et al. 1998; Biskamp & Müller 1999; Ostriker et al. 2001; Cho et al. 2002).

Even if the consensus that turbulence quickly decays has been numerically established for the last two decades, the previous numerical results rely heavily on the isothermal condition. However, as long as various density and temperature phases in the ISM (Ferrière 2001) are concerned, the use of a polytropic equation of state (EOS)

(1) |

where *P* is the pressure, *ρ* is the density, and both *K* and *γ* are constants, is a valid approach (see Vazquez-Semadeni et al. 1996 and references therein). The polytropic EOS has been used for many astrophysical problems, such as complex chemical processes (Spaans & Silk 2000; Glover & Mac Low 2007a) or turbulence (Scalo et al. 1998; Li et al. 2003; Glover & Mac Low 2007b; Federrath & Banerjee 2015).

Besides a variety of density and temperature phases, a wide range of driving agents of turbulence also characterizes interstellar turbulence (see Federrath et al. 2017 for a review). Based on its compressibility, we may consider two extreme types of driving: solenoidal (divergence-free) and compressive (curl-free). Until recently, solenoidal driving had been mainly used for turbulence studies. However, Federrath et al. (2010) showed that compressive driving and solenoidal driving can have different statistics. For example, they showed that “the former yields stronger compression at the same RMS Mach number than the latter, resulting in a three times larger standard deviation of volumetric and column density probability distributions.” To the best of our knowledge, scaling relations of decaying polytropic turbulence initially driven by compressive driving have not been studied yet.

The main goal of this paper is to examine whether the decay exponent *α* depends on the value of the polytropic exponent *γ*. Here we concentrate on decay of polytropic turbulence driven by either solenoidal or compressive driving in both transonic and supersonic regimes. Hence, we expect to demonstrate how polytropic EOS and types of driving affect decaying turbulence. In addition, we also investigate the probability density function (PDF) of gas density and skewness of the PDF in decaying polytropic turbulence initially driven by compressive driving.

The paper is organized as follows. We explain our motivation and numerical method in Section 2, and present the results from our numerical simulations in Section 3. We discuss our finding and its astrophysical implication and give a summary in Section 4.

2. MOTIVATION AND NUMERICAL METHOD

2.1. Motivation

As we described earlier, decay of solenoidally driven isothermal turbulence follows *E* ∝ *t*^{−α} with *α* ≈ 1. The type of driving or *γ* may affect this scaling relation.

First, compressive driving yields more compressions at the same Mach number. Therefore, while decaying, compressed regions could generate additional kinetic energy via expansion, which could affect the rate of turbulence decay.

Second, regarding the effects of *γ* on decaying turbulence, only limited regions of the parameter space have been studied. Mac Low et al. (1998) found that supersonic turbulence with *γ* = 1.4 decays with *α* ~ 1.2. For isothermal cases (i.e., *γ* = 1), they found that *α* is nearly unity. This suggests that the scaling exponent *α* in *E* ∝ *t*^{−α} only weakly depends on *γ* as assumed by Davidovits & Fisch (2017). However, Mac Low et al. (1998) used a random initial velocity perturbation which follows a power-law, and a constant initial density. Therefore, it is necessary to test the decay law using initial velocity and density data cubes from actual turbulence simulations with both soft EOS (i.e., polytropic *γ* < 1) and stiff EOS (i.e., polytropic *γ* > 1).

Third, the effects of *γ* and the type of driving on density PDF of turbulence have been addressed in several previous studies. For example, Federrath et al. (2010) showed that solenoidal driving and compressive driving can produce different statistics of isothermal turbulence. In addition, Federrath & Banerjee (2015) found that the density PDF of solenoidally driven turbulence with polytropic *γ* = 5/3 has a power-law tail at low density, which is not observed in isothermal turbulence. However, earlier studies have not addressed turbulence with polytropic EOS and compressive driving. In this paper, we use both compressive driving and polytropic EOS to investigate density PDF and skewness of driven and decaying turbulence.

2.2. Numerical Method

2.2.1. Numerical Code

We use an Essentially Non-Oscillatory (ENO) scheme (see Cho & Lazarian 2002) to solve the ideal hydrodynamic equations in a periodic box of size 2*π* (≡ *L*_{0}):

(2) |

(3) |

where f is a driving force, *P* is the pressure (see Section 2.2.2), ρ is the density, and v is the velocity. The density and velocity are set to 1 and zero at *t* = 0 to assume a static medium with a constant density at the beginning.

2.2.2. Simulations

^{a} Driving schemes – either solenoidal or compressive driving.^{b} Polytropic exponent.^{c} The driving wavenumber at which the energy injection rate peaks.^{d} The sonic Mach number which is defined in Equation (5).^{e} Average of power law exponent for *v*^{2} ∝ *t*^{−α1} .^{f} Average of power law exponent for *σ*_{ρ/ρ0} ∝ *t*^{−α2} .^{g} Time interval in units of large-eddy turnover time (*t*_{ed}) for averaging *α*_{1} and *α*_{2}.

We use 512^{3} grid points in our periodic computational box. The peak of energy injection occurs at *k* ~ *k _{f}* , where

In both drivings, the driving vectors continuously change with a correlation time comparable to the large-eddy turnover time. We also adopt a polytropic EOS

(4) |

where *P* is the normalized pressure, and *P*_{0}, *c*_{s0}, and *ρ*_{0} are the initial pressure, sound speed, and density, respectively. The sonic Mach number *M _{s}* is defined by

(5) |

where *v*_{rms} is the rms velocity. We vary *γ* and the sonic Mach number *M _{s}* to consider both soft and stiff EOS in transonic and supersonic regimes.

Table 1 lists our simulation models. We use the notation XMSY-*γ*Z, where X = S or C refers to solenoidal (S) or compressive (C) driving; Y = 1, 3, or 5 refers to the sonic Mach number *M _{s}*; Z = 0.7, 1.0, 1.5, or 5/3 refers to the value of

Table 1

Simulation runs and parameter values

Simulation runs and parameter values

Run | Driving^{a} |
γ^{b} |
Resolution | k^{c}_{f} |
M^{d}_{s} |
α^{e}_{1} |
α^{f}_{2} |
(t_{1}, t_{2})^{g} |
---|---|---|---|---|---|---|---|---|

CMS5-γ1.0 |
Compressive | 1.0 | 512^{3} |
8.0 | ~5 | 0.767 | 0.513 | (5,10) |

We did not include a stiff EOS simulation (i.e., *γ* > 1) for compressively driven turbulence with *M*_{s} ≿ 3 and for solenoidally driven turbulence with *M*_{s} ≿ 5. This is because the simulations crashed before turbulence saturates due to the high *M*_{s} and large *γ*, which can lead to extremely low density values (see the PDFs in Figures 8 and 9 for *γ* > 1). Therefore, we could not run CMS3-*γ*5/3, CMS5-*γ*5/3, or SMS5-*γ*5/3. When *γ* = 5/3 simulations are unstable, we use *γ* = 1.5 instead for SMS5 and CMS3. Using different *γ* for the stiff EOS would not be appropriate for a comparison. Nevertheless, we think that the difference between *γ* = 1.5 and 5/3 would not significantly affect the decay law. Even *γ* = 1.5 is unstable for CMS5. Therefore, we do not have data for CMS5 with a stiff EOS.

3. RESULTS

3.1. Decay of Hydrodynamic Turbulence with Polytropic EOS

3.1.1. Decay of Turbulence Driven by Solenoidal Driving

Figure 1.

Time evolution of <*v*^{2}> in decaying turbulence initially driven by solenoidal driving with *γ* = 0.7 (blue), 1.0 (red), 1.5 (cyan), or 5/3 (green), where <· · ·> indicates spatial averaging. *Left:* *M*_{s} ~ 1. *Center:* *M*_{s} ~ 3. *Right:* *M*_{s} ~ 5. Turbulence starts decaying at *t*_{code}/*t*_{ed} = 0. We normalize <*v*^{2}> by υ 0 2 which is the value at *t*_{code}/*t*_{ed} = 0. The black dotted lines in all three panels are reference lines for different power-law exponents.

Figure 2.

Same as Figure 1, but for the decay of the standard deviation of the density fluctuation*σ*_{ρ/ρ0} . We normalize *σ*_{ρ/ρ0} by (*σ*_{ρ/ρ0})_{0} which is the value at *t*_{code}/*t*_{ed} = 0.

Here, we consider decaying polytropic turbulence initially driven by solenoidal driving. Figures 1 and 2 show the decay of <*v*^{2}>^{1} and the standard deviation of density fluctuations *σ*_{ρ/ρ0} , respectively, where <· · ·> denoted spatial averaging. From left to right, the sonic Mach number *M _{s}* is ~ 1, ~ 3, and ~ 5, respectively. Blue, red, cyan, and green curves in each panel correspond to

Figure 1.

Time evolution of <

Figure 2.

Same as Figure 1, but for the decay of the standard deviation of the density fluctuation

First of all, we clearly see that the decay of <*v*^{2}> follows a power-law of the form <*v*^{2}> ∝ *t*^{−α}, and *α* is almost the same at the same *M _{s}* regardless of the value of

3.1.2. Decay of Turbulence Driven by Compressive Driving

Figure 3.

Same as Figure 1, but for compressively driven turbulence. Note that unlike Figure 1, <*v*^{2}> slightly increases along the decay as indicated by a black arrow in each panel; this is most apparent in the case of *M*_{s} ~ 1 (left panel).

We now deal with the decay of polytropic turbulence initially driven by compressive driving. Figures 3 and 4 show the decay of <*v*^{2}> and *σ*_{ρ/ρ0} , respectively. As in the case of solenoidal turbulence,^{2} we use different values of *M _{s}* (from left to right panels) and

Figure 3.

Same as Figure 1, but for compressively driven turbulence. Note that unlike Figure 1, <

Similar to the result from solenoidal turbulence, both <*v*^{2}> and *σ*_{ρ/ρ0} in compressively driven turbulence exhibit power-law decay, and *γ* hardly affects the decay rate (see Table 1 for measurements of *α*). According to Figure 3, the power-law exponent *α* is ~1.0 for *M _{s}* ~ 1, and ~0.8 for

Unlike in the case of solenoidal turbulence, *γ* slightly affects the decay of compressively driven turbulence. First, <*v*^{2}> in compressively driven turbulence shows bump-like features (indicated by the black arrow in each panel of Figure 3). This slight increase of <*v*^{2}> is most pronounced in the case *M _{s}* ~ 1. Second, dip-like features (indicated by the black arrow in each panel of Figure 4) are visible in the evolution of

3.2. Decay of Solenoidal and Compressive Velocity Components in Turbulence Driven by Compressive and Solenoidal Drivings

Figure 5.

Decay of the ratio of compressive (upper panels) and solenoidal (bottom panels) velocity components in compressively driven turbulence.*Left panels:* *M*_{s} ~ 1. *Center panels:* *M*_{s} ~ 3. *Right panels:* *M*_{s} ~ 5. Blue, red, cyan, and green curves denote a *γ* of 0.7, 1.0, 1.5, and 5/3, respectively.

Figure 6.

Decay of solenoidal and compressive velocity components in compressively driven turbulence. Only the case of isothermal turbulence (i.e.,*γ* = 1.0) is presented here. Left: *M*_{s} ~ 1. *Center:* *M*_{s} ~ 3. *Right:* *M*_{s} ~ 5. Black, cyan, and magenta curves indicate total, solenoidal, and compressive velocity components, respectively. The dotted lines in different colors are reference lines for different power-law forms. Note that the compressive velocity component decays much faster than the solenoidal one.

Figure 7.

Same as Figure 5 but for solenoidally driven turbulence and showing the compressive ratio only.

Here, we deal with the differences in the decay of solenoidal and compressive modes. We decompose 3D velocity fields of both solenoidally and compressively driven turbulences into solenoidal and compressive components.^{3} Figure 5 illustrates the decay of the compressive ratio _{sol} and v_{comp} are the solenoidal and compressive velocity components, respectively. Blue, red, cyan, and green curves indicate values of *γ* of 0.7, 1.0, 1.5, and 5/3, respectively.

Figure 5.

Decay of the ratio of compressive (upper panels) and solenoidal (bottom panels) velocity components in compressively driven turbulence.

First, Figure 5 shows that the compressive ratio decreases as turbulence decays. When *M _{s}* ~ 1 or ~ 3, the smaller

Figure 6.

Decay of solenoidal and compressive velocity components in compressively driven turbulence. Only the case of isothermal turbulence (i.e.,

While the compressive ratio in compressively driven turbulence becomes less than the solenoidal ratio as the turbulence decays, Figure 7 shows that the compressive ratio in solenoidal turbulence does not. In the figure, blue, red, cyan, and green curves indicate *γ* of 0.7, 1.0, 1.5, and 5/3, respectively. As we can see from Figure 7, the compressive ratio in solenoidal turbulence is always less than 0.5 regardless of both *M _{s}* and polyrtropic

Figure 7.

Same as Figure 5 but for solenoidally driven turbulence and showing the compressive ratio only.

3.3. Density PDF and Skewness

Figure 8.

PDF of the logarithmic density*s* = ln(*ρ*/*ρ*_{0}) of decaying polytropic turbulence initially driven by compressive driving with *M*_{s} ~ 1. *Left:* polytropic *γ* = 0.7. *Center:* *γ* = 1.0. *Right:* *γ* = 5/3. Note that the *p*(*s*) axis is logarithmic. Black, purple, and orange curves represent different times along the decay. Dotted lines are log-normal fitting line.

Figure 10.

Time evolution of the skewness of the density PDFs shown in Figure 8 (top panel) and Figure 9 (bottom panel). Blue, red, cyan, and green curves correspond to*γ* of 0.7, 1.0, 1.5, and 5/3, respectively.

In the following, we investigate the density PDF and its skewness for the case of compressively driven turbulence with polytropic EOS in driven and decay regimes. We define the skewness of the density PDF as

(6) |

where *N* is the total number of data points, *s* ≡ ln(*ρ*/*ρ _{0}*) is the natural logarithm of the density fluctuations, and <· · ·> denotes the spatial average. Skewness measures the asymmetry of a probability distribution. When the distribution is left-skewed (right-skewed), the skewness takes a negative (positive) value.

Figures 8 and 9 show the density PDFs of polytropic turbulence with *M _{s}* ~ 1 and ~ 3, respectively. Each colored solid curve corresponds to the density PDF at a different time along the decay. We carry out log-normal fitting for the PDFs using the relation

(7) |

Figure 8.

PDF of the logarithmic density

In the figures, fitting lines are indicated by dotted lines.

We first consider driven turbulence (black solid line in Figures 8 and 9). Compressively driven turbulence yields density PDFs which are not exactly log-normal even in the case *γ* = 1. The density PDF for *γ* = 0.7 is slightly right-skewed, and that for *γ* = 5/3 is strongly left-skewed. The latter has a pronounced power-law tail at low density as shown in Figure 8c. Figure 9c shows that the PDF of CMS3-*γ*1.5, which is for *M _{s}* ~ 3, deviates more strongly from the log-normal form at low density. However, it is not clear whether the low density tail follows a power-law.

We can interpret the behavior of density PDF with *γ* stated above as follows. When *γ* < 1 (*γ* > 1), compression decreases (increases) internal temperature and thus local sound speed, which allows (hinders) high density regions to develop easily (Federrath & Banerjee 2015). Since compression can commonly occur in our simulations of transonic/supersonic turbulence with compressive driving, the density PDFs with *γ* < 1 (*γ* > 1) become right-skewed (left-skewed). In addition, although density PDFs with *γ* = 0.7 are *slightly* right-skewed, we expect that they become more right-skewed as *γ* decreases further. In fact, Scalo et al. (1998) showed a clear power-law tail toward high density regions in their one-dimensional simulation with *γ* = 0.3 and *M _{s}* ~ 3.

Note that our results are for compressively driven turbulence with polytropic EOS. Earlier studies are available for solenoidal turbulence with polytropic EOS and *M _{s}* ~ 10 (Federrath & Banerjee 2015), and for compressively driven turbulence with

Figure 10 shows the time evolution of the skewness of the density PDF for compressively driven turbulence. The horizontal axis denotes the elapsed decay time normalized by the large-eddy turnover time. The top and bottom panels of Figure 10 show the skewness for *M _{s}* ~ 1 and ~ 3, respectively. At

Figure 10.

Time evolution of the skewness of the density PDFs shown in Figure 8 (top panel) and Figure 9 (bottom panel). Blue, red, cyan, and green curves correspond to

Run | t^{a}_{1} |
t^{b}_{2} |
t^{c}_{3} |
---|---|---|---|

CMS1-γ0.7CMS1- γ1.0CMS1- γ5/3CMS3- γ0.7CMS3- γ1.0CMS3- γ1.5 |
0.260 −0.335 −1.576 0.204 −0.301 −0.784 |
0.084 0.290 0.168 0.197 −0.149 −0.382 |
0.294 0.260 0.260 0.093 0.070 −0.013 |

Next, we consider the decay regime. As can be seen from the purple and orange solid curves in Figures 8 and 9, as turbulence decays, the density PDFs become narrow and approach log-normal forms in all cases. We can confirm this trend in Figure 10. The skewness for *M _{s}* ~ 1 and ~ 3 approaches, and fluctuates around, zero as turbulence decays, which is consistent with the temporal change of the PDFs shown in Figures 8 and 9.

4. DISCUSSION AND SUMMARY

The purpose of this study is to investigate the effects of the equation of state (i.e., the value of the polytropic index *γ*) and driving schemes (i.e., solenoidal and compressive drivings) on decaying turbulence and its statistics. In this paper, we demonstrated that the scaling relation for the decay law (<*v*^{2}> ∝ *t*^{−α}) does not show a strong dependence on *γ* and the driving schemes. Throughout our simulations, <*v*^{2}> decays with 0.8 ≾ *α* ≾ 1.2. The range is nearly same as the one found by Mac Low et al. (1998) (0.85 < *α* < 1.2).

For *γ* > 1 cases, *α* ranges from 1.0 to 1.2 in solenoidal turbulence and from 0.8 to 1.0 in compressively driven turbulence, with the largest value of *α* being obtained for the case *M _{s}* ~ 1 in both driving schemes. This result confirms the assumption of Davidovits & Fisch (2017) that for

Even though no relationship between polytropic index *γ* and scaling of the decay law (<*v*^{2}> ∝ *t*^{−α}) is found in our study, there are several noticeable characteristics in case of compressively driven turbulence. First, a slight increase of <*v*^{2}> and an associated decrease of *σ*_{ρ/ρ0} are found in Figures 3 and 4, respectively. As discussed earlier, those effects can be interpreted as being due to additional energy released from compressed regions via expansion. Second, in the case *M _{s}* ~ 1, the effect is most pronounced. This may be due to the pressure being strong compared to that of supersonic cases, resulting in stronger expansion. Third, for the same

More interestingly, the compressive velocity component decays faster than the solenoidal one in the case of turbulence initially driven by compressive driving as shown in Section 3.2. We can interpret this as follows. When turbulence initially driven by compressive driving decays, the energy of the compressive component is dissipated through both turbulent cascades and dissipation at shocks, and a fraction of the energy is converted into solenoidal energy. For the solenoidal component, only turbulent cascades can dissipate the energy. As it has fewer channels for energy dissipation available than the compressive component, the solenoidal velocity component in turbulence initially driven by compressive driving decays more slowly. However, a detailed analysis addressing, e.g., what fraction of compressive kinetic energy changes into solenoidal kinetic energy, is beyond the scope of this paper; further studies are required to understand this issue quantitatively.

Our study uses polytropic EOS to quantify decaying astrophysical turbulence. A fluid with polytropic EOS (*P* = *Kρ*^{γ}) is different from an ideal gas with adiabatic index *γ* > 1 as follows. First, we do not need to solve the energy conservation equation for the former, while we have to for the latter. Note that the EOS for the latter becomes *P* = (*γ* − 1)*E*, where *E* is the internal energy density. Second, they provide different descriptions of shocks. The polytropic EOS is only applicable to isentropic flows; the entropy is constant everywhere, even in shocked regions. The entropy is not constant for an ideal gas. Although entropy is conserved in regions with adiabatic expansion/compression, it increases when a parcel of gas passes through a shock front in an ideal gas. Therefore, a fluid with polytropic EOS and an ideal gas with adiabatic index *γ* have different descriptions of shocks and their shock jump conditions are different (see Bisnovatyi-Kogan & Moiseenko 2016 for shock jump conditions for adiabatic and isentropic flows).

Our results have astrophysical implications. The polytropic exponent *γ* is useful to describe a variety of components in the ISM. A polytropic EOS with *γ* ⋍ 0.8 can represent the density range of 10*cm*^{−3} ≤ *n* ≤ 10^{4}*cm*^{−3} , where *n* is the hydrogen number density (Glover & Mac Low 2007b); this corresponds to the density range for giant molecular clouds (Ferrière 2001). Furthermore, the EOS with *γ* ~ 1.4 describes the conditions in the centers of protostellar cores, with a corresponding density range 10^{12}*cm*^{−3} ≤ *n* ≤ 10^{17}*cm*^{−3} (Masunaga & Inutsuka 2000). Accordingly, our results suggest that once driving of turbulence ceases to act, turbulence quickly decays with a dynamical timescale that is characteristic of a certain system irrespective of its spatial scale. Moreover, even if turbulence is initially driven by compressive driving like supernova explosions, solenoidal motions will dominate as the turbulence decays due to the much faster decay of compressive motions.

In summary, we have studied the influence of polytropic EOS and driving schemes on decaying turbulence and its statistics and found the following results.

- 1. There is no significant correlation between the scaling of the decay law (
*E*∝*t*^{−α}) and polytropic exponent*γ*in the case of solenodially driven turbulence. - 2. Driving schemes have a non-negligible effect on the decay rate of turbulence: the power-law index
*α*for turbulence initially driven by compressive driving is smaller than that for turbulence initially driven by solenoidal driving. - 3. There is no significant effect of
*γ*on the decay rate of velocity in compressively driven turbulence. - 4. The polytropic exponent
*γ*has a small effect on the density fluctuations in compressively driven turbulence: the smaller*γ*is, the faster the standard deviation of density fluctuations of the turbulence decays. - 5. When considering the decay of solenoidal and compressive velocity components in compressively driven turbulence separately, the energy of the compressive velocity component decays much faster.
- 6. Regarding the statistics of compressively driven turbulence, the density PDF deviates from a log-normal distribution, especially for
*γ*> 1. In addition, we have found that the skewness of the density PDF of the turbulence becomes zero as it decays.

Notes
^{1}We are interested in the decay of <(*δv*)^{2}>, where *δv* is the magnitude of the velocity fluctuations. If <v> is non-zero, <(*δv*)^{2}> is different from <*v*^{2}>. That is, <*v*^{2}> = <(*δv*)^{2}>+<*v*_{x}>^{2} +<*v*_{y}>^{2} +<*v*_{z}>^{2} , where <*v*_{i}> (*i* = *x*, *y*, *z*) is the mean velocity component over the whole computational box. In Runs SMS1-*γ*1.0, SMS3-*γ*1.0, and SMS5-*γ*1.0, <*v*_{x}>, <*v*_{y}>, and <*v*_{z}> are less than ~0.01 at the moment turbulence starts decaying. Therefore, we can ignore the effect of non-zero <*v*_{i}>. We assume that this holds true also for the other simulations presented in this paper.^{2}By “solenoidal turbulence” we mean turbulence initially driven by solenoidal driving.^{3}To be specific, we carry out Fourier transforms of total velocity fields first. Then, we take only the velocity components perpendicular to k, which correspond to solenoidal velocity component in Fourier space, and carry out inverse Fourier transforms to produce solenoidal velocity fields in real space. Finally, we subtract the solenoidal velocity components from the total velocities in real space to yield the compressive velocity components.

Acknowledgments

This paper has been expanded from a chapter of Jeonghoon Lims MSc thesis. This work is supported by the National R&D Program through the National Research Foundation of Korea Grants funded by the Korean Government (NRF-2016R1A5A1013277 and NRF-2016R1D1A1B02015014).

References

1. |
Biskamp, D., & Müller, W.-C. 1999, Decay Laws for Three-Dimensional Magnetohydrodynamic Turbulence, Phys. Rev. Lett., 83, 2195 |

2. |
Bisnovatyi-Kogan, G. S. & Moiseenko, S. G. 2016, Isentropic “Shock Waves” in Numerical Simulations of Astrophysical Problems, Astrophys., 59, 1 |

3. |
Cho, J. & Lazarian, A. 2002, Compressible Sub-Alfvénic MHD Turbulence in Low-β Plasmas, Phys. Rev. Lett., 88, 245001 |

4. |
Cho, J., Lazarian, A., & Vishniac, E. T. 2002, Simulations of Magnetohydrodynamic Turbulence in a Strongly Magnetized Medium, ApJ, 564, 291 |

5. |
Davidovits, S. & Fisch, N. J. 2017, A Lower Bound on Adiabatic Heating of Compressed Turbulence for Simulation and Model Validation, ApJ, 838, 118 |

6. |
Federrath, C., Roman-Duval, J., Klessen, R. S., et al. 2010, Comparing the Statistics of Interstellar Turbulence in Simulations and Observations. Solenoidal versus Compressive Turbulence Forcing, A&A, 512, A81 |

7. |
Federrath, C. & Banerjee, S. 2015, The Density Structure and Star Formation Rate of Non-isothermal Polytropic Turbulence, MNRAS, 448, 3297 |

8. |
Federrath, C., Rathborne, J. M., Longmore, S. N., et al. 2017, The Link between Solenoidal Turbulence and Slow Star Formation in G0.253+0.016, Proc. IAU Symp. 322, 123 |

9. |
Ferrière, K. M. 2001, The Interstellar Environment of our Galaxy, Rev. Mod. Phys., 73, 1031 |

10. |
Glover, S. C. O. & Mac Low, M.-M. 2007, Simulating the Formation of Molecular Clouds. I. Slow Formation by Gravitational Collapse from Static Initial Conditions, ApJS, 169, 239 |

11. |
Glover, S. C. O. & Mac Low, M.-M. 2007, Simulating the Formation of Molecular Clouds. II. Rapid Formation from Turbulent Initial Conditions, ApJ, 659, 1317 |

12. |
Goldreich, P. & Sridhar, S. 1995, Toward a Theory of Interstellar Turbulence. II. Strong Alfvenic Turbulence, ApJ, 438, 763 |

13. |
Larson, R. B. 1981, Turbulence and Star Formation in Molecular Clouds, MNRAS, 194, 809 |

14. |
Lesieur, M. 2008, Turbulence in Fluids (Dordrecht: Springer) |

15. |
Li, Y., Klessen, R. S., & Mac Low, M.-M. 2003, The Formation of Stellar Clusters in Turbulent Molecular Clouds: Effects of the Equation of State, ApJ, 592, 975 |

16. |
Mac Low, M.-M., Klessen, R. S., Burkert, A., & Smith, M. D. 1998, Kinetic Energy Decay Rates of Supersonic and Super-Alfvénic Turbulence in Star-Forming Clouds, Phys. Rev. Lett., 80, 2754 |

17. |
Mac Low, M.-M. & Klessen, R. S. 2004, Control of Star Formation by Supersonic Turbulence, Rev. Mod. Phys., 76, 125 |

18. |
Masunaga, H. & Inutsuka, S. 2000, A Radiation Hydrodynamical Model for Protostellar Collapse. II. The Second Collapse and the Birth of a Protostar, ApJ, 531, 350 |

19. |
Ostriker, E. C., Stone, J. M., & Gammie, C. F. 2001, Density, Velocity, and Magnetic Field Structure in Turbulent Molecular Cloud Models, ApJ, 546, 980 |

20. |
Padoan, P., Nordlund, A., & Jones, B. J. T. 1997, The Universality of the Stellar Initial Mass Function, MNRAS, 288, 145 |

21. |
Padoan, P. & Nordlund, A. 2002, The Stellar Initial Mass Function from Turbulent Fragmentation, ApJ, 576, 870 |

22. |
Passot, T. & Vázquez-Semadeni, E. 1998, Density Probability Distribution in One-dimensional Polytropic Gas Dynamics, Phys. Rev. E, 58, 4501 |

23. |
Scalo, J., Vázquez-Semadeni, E., Chappell, D., & Passot, T. 1998, On the Probability Density Function of Galactic Gas. I. Numerical Simulations and the Significance of the Polytropic Index, ApJ, 504, 835 |

24. |
Spaans, M. & Silk, J. 2000, The Polytropic Equation of State of Interstellar Gas Clouds, ApJ, 538, 115 |

25. |
Stone, J. M., Ostriker, E. C., & Gammie, C. F. 1998, Dissipation in Compressible Magnetohydrodynamic Turbulence, ApJL, 508, L99 |

26. |
Vazquez-Semadeni, E., Passot, T., & Pouquet, A. 1996, Influence of Cooling-induced Compressibility on the Structure of Turbulent Flows and Gravitational Collapse, ApJ, 473, 881 |