GROWTH TIME OF ACOUSTIC PERTURBATIONS IN ISENTROPICALLY UNSTABLE HEAT-RELEASING MEDIUM Текст научной статьи по специальности «Физика»

Scientific article

DOI: 10.18287/2541-7525-2022-28-1-2-113-119 Submitted: 06.09.2022

Revised: 12.10.2022 Accepted: 14.11.2022

D.S. Riashchikov

Lebedev Physical Institute, Samara, Russian Federation Samara National Research University, Samara, Russian Federation E-mail: ryashchikovd@gmail.com. ORCID: https://orcid.org/0000-0001-7143-2968

I.A. Pomelnikov

Samara National Research University, Samara, Russian Federation E-mail: vanidzepomelnikov@gmail.com. ORCID: https://orcid.org/0000-0001-7839-5784

N.E. Molevich

Lebedev Physical Institute, Samara, Russian Federation Samara National Research University, Samara, Russian Federation E-mail: nonna.molevich@mail.ru. ORCID: https://orcid.org/0000-0001-5950-5394

GROWTH TIME OF ACOUSTIC PERTURBATIONS IN ISENTROPICALLY UNSTABLE HEAT-RELEASING MEDIUM1

ABSTRACT

Isentropic instability is a type of thermal instability that leads to the growth of acoustic waves. As a result of wave growth in such media, autowave structures are formed, the parameters of which depend only on the properties of the medium and can be predicted both analytically and numerically. This study aims to answer the question of how quickly these structures can form in an isentropically unstable medium with parameters similar to Orion Bar. It is shown that the growth time depends on the characteristic size of the initial perturbation. The fastest growing structures take 3-6 thousand years to reach half their maximum amplitude. Further growth to the maximum value takes 15-20 thousand years.

Key words: instability; thermal instability; nonlinear waves; shock waves; autowaves; interstellar gas; photodissociation region; Orion Bar.

Citation. Riashchikov D.S., Pomelnikov I.A., Molevich N.E. Growth time of acoustic perturbations in isentropically unstable heat-releasing medium. Vestnik Samarskogo universiteta. Estestvennonauchnaia seriia = Vestnik of Samara University. Natural Science Series, 2022, vol. 28, no. 1-2, pp. 113-119. DOI: http://doi.org/10.18287/2541-7525-2022-28-1-2-113-119.

Information about the conflict of interests: authors and reviewers declare no conflict of interests.

© Riashchikov D.S., 2022

Dmitrii S. Riashchikov — Candidate of Physical and Mathematical Sciences, research associate of the Theoretical Department of Lebedev Physical Institute, 221, Novo-Sadovaya Street, Samara, 443011, Russian Federation; senior lecturer of the Department of Physics, Samara National Research University, 34, Moskovskoye shosse, Samara, 443086, Russian Federation. ©c Pomelnikov I.A., 2022

Ivan A. Pomelnikov — student of the Institute of IT and Cybernatics, Samara National Research University, 34, Moskovskoye shosse, Samara, 443086, Russian Federation. ©c Molevich N.E., 2022

Nonna E. Molevich — Doctor of Physical and Mathematical Sciences, chief researcher of the Theoretical Department of Lebedev Physical Institute, 221, Novo-Sadovaya Street, Samara, 443011, Russian Federation;

1The work was supported in part by the Ministry of Education and Science (projects FSSS-2020-0014, 0023-2019-0003).

professor of the Department of Physics, Samara National Research University, 34, Moskovskoye shosse, Samara, 443086, Russian Federation.

Introduction

Isentropic instability is a type of thermal instability that leads to the amplification of acoustic waves. It may occur in a medium with heating r(p, T) and cooling L(p,T) processes, which powers depend on the density and temperature. In a state of equilibrium po,To they compensate each other, so r(po,To) = L(po,To). Acoustic perturbations violate this equilibrium, and heat release may further amplify (instability) or suppress them.

Conditions for isentropic instability may exist in photodissociation regions [1-3]. The instability results in a periodic wave structure, which in the later stages of evolution is a sequence of shock waves with autowave properties [4-6]. Similar structures are found, for example, near RCW120 [7; 8] or Orion nebulae [9; 10].

The fundamentals of the theory of thermal instabilities were thoroughly developed by Field [11]. His and many later studies have focused on the dispersion properties of gasdynamic perturbations in heat-releasing media. This allows us to judge the initial stage of perturbation evolution. Here, it is worth noting the work [12] in which a linear equation is obtained and analytically solved, which allows us to study the behavior of acoustic waves at this stage depending on the parameters of the initial perturbation.

However, linear approach works only for small amplitude waves. A study of the subsequent evolution of the acoustic waves at the nonlinear stage using the nonlinear equation [1; 5] for small amplitude waves and numerical simulations revealed the formation of autowave structures which parameters do not depend on the parameters of the initial perturbation. The parameters of the autowave structures were afterwards analytically evaluated in [6] without any restrictions on the amplitude.

The aforementioned studies allow us to answer the question of what structures can be observed under the known parameters of the medium and the heating and cooling functions, but do not answer the question of how quickly these structures can emerge. The importance of this issue stems from the fact that acoustic waves propagate during amplification, and the size of the medium in which their amplification can occur may be limited. Thus, the acoustic waves may not have time to reach the predicted amplitudes.

In this work, we numerically estimate the growth time of autowave structures that can emerge in photodissociation region Orion Bar using the model of heating and cooling functions from [3].

1. Time of the formation of autowave pulse

Isentropic instability, as mentioned above, leads to the formation of a periodic structure. Ahead of this structure the so-called autowave pulse propagates, which is the final stage of evolution of any small acoustic perturbation in an isentropically unstable medium. The parameters of this pulse depend only on form of heating and cooling functions of the medium from temperature and density. In this section, we estimate the time of formation of the autowave pulse in Orion Bar.

The dynamics of acoustic perturbations in such a medium can be described by the following system of equations:

ff + div (pv) = 0,

p dv = -VP' nn

P = m pT,

m ' '

where p,T,P are density, temperature, and pressure, v is speed vector, cv is heat capacity under constant volume, kB is Boltzmann constant; d/dt is substantial derivative, W(p,T) = L(p,T) — r(p, T) is generalized heat-loss function.

In the current work, we use the model of heat-loss function W(p, T) proposed by [3] for photodissociation regions of the interstellar medium. The following parameters of heat-loss function were used: FUV field Go = = 4 x 104, cooling line opacity tc = 0.5, the ratio of visual extinction to reddening rv = 5.5, C abundance per H nucleus in very small grains bC = 3 x 10~5 [3], the abundances of carbon £C = 1.2 x 10~4 and oxygen

= 2.56 x 10-4 [13].

Studies show that the maximum growth rate of acoustic waves is expected at temperatures about 1000 K [3] which also lies in the limit of observable temperatures in Orion Bar. So, we use this temperature To = = 1000^ as an equilibrium one. Then, using the equilibrium condition W (po, To) = W (no, To) = 0, one can find the corresponding equilibrium number density in Orion Bar as no = 2.26 x 105crn-3.

Since there is no analytical solution capable of describing the evolution of an arbitrary perturbation in a medium with strong dispersion and nonlinearity, we will investigate the growth time numerically using the Athena MHD code [14] for astrophysical simulations.

The initial condition for numerical simulations is Gaussian perturbation in form

P = Po

(1.2)

where a is the dimensionless amplitude of density perturbation, a is the characteristic size of perturbation, y is the adiabatic index which equals 5/3 in Orion Bar. In our simulations, we used initial amplitude a = 0.01.

Simulations show the splitting of initial perturbation into two waves propagating in the opposite directions. Then, from any single wave, the periodic wave structure with period determined by heat-loss function appear. Each wave in this structure grows forming several autowave pulses in front of the wave sequence and a series of smaller waves at an earlier stage of evolution behind them (Figure 1.1).

0.00

0.05

0.10

0.15

x, pc

Fig. 1.1. Numerical simulation of autowave pulse formation in isentropically unstable medium Рис. 1.1. Численное моделирование формирования автоволнового импульса в изоэнтропически

неустойчивой среде

Using numerical simulation, we determine the amplitude of the autowave pulse as a function of time. Since the amplitude of the autowave pulses in front of wave sequence is greater than the amplitude of waves behind them, we take as the pulse amplitude the maximum value of the amplitude of the waves in the sequence.

When conducting numerical simulation, the important issue is to investigate the influence of grid step on the time of the formation of autowave pulses (Figure 1.2). The need for this is due to the fact that the numerical scheme introduces diffusion, much larger than that observed in real media. Decreasing the grid step reduces this effect, but it is impossible to reduce the numerical diffusion to an order of magnitude observed in real media due to limited computational capabilities.

One can see from Figure 1.2 that the finer the grid, the faster the waves grow. Let us note that three plots with a coarser grid have a long interval with almost unchanged amplitude. We suppose this is due to the dispersion properties of the medium and the spectrum of the initial perturbation.

Heating and cooling processes acts at characteristic time tq and its corresponding characteristic length lq [5]:

CW Y

TQ

2n-

^ , Lq = VхквTo/mTQ,

(1.3)

Wto\ 1 - (PoWpo) / (T0Wto), where WT0 = (dW/dT)T=Tg p=pg, Wp0 = (dW/dp)T=Tg p=pg. Estimations of these parameters for the Orion Bar give tq = 875 years and Lq = 2.6 x 10~3 pc.

In a heat-releasing medium, the high-frequency harmonics of the initial perturbation (vtq ^ 1, where v is the frequency of the wave) have the largest increment, while the low-frequency harmonics (vtq ^ 1) have a relatively small one. Scheme diffusion significantly slows down growth of high-frequency harmonics. And if an amplitude of high-frequency harmonics in the initial perturbation is small, then for their growth it is required sufficiently long time, and their contribution to amplitude of structure as a whole will be imperceptible for this time. This effect is clearly seen in Figure 1.3 where the amplitude of perturbations with a large initial

1.2-

10000 20000 30000

г, years

40000

50000

Fig. 1.2. Influence of the grid step on the time of the autowave pulse formation. Characteristic size of initial

perturbation a = 1.0Lq

Рис. 1.2. Влияние шага по координате в численном моделировании на время формирования автоволнового импульса. Характерный размер начального возмущения a = 1.0Lq

characteristic size begins to increase later than that of perturbations with a smaller characteristic size, with the same numerical grid parameters for all perturbations.

«1 e 1-2

20000

40000 r, years

60000

80000

Fig. 1.3. Influence of the characteristic size of initial perturbation on the time of the autowave formation. Grid

step Ax = 0.015Lq

Рис. 1.3. Влияние характерного размера начального возмущения на время формирования автоволнового

импульса. Шаг сетки по координате Ax = 0.015Lq

Thus, we believe that as the grid size is further reduced, the growth rate of acoustic waves until the amplitude reaches half of the maximum value will tend to that predicted by Field's theory for inviscid heat-releasing medium [11]. Characteristic time of isentropic instability tinst (the time at which the amplitude of the waves increases by a factor of e) of the high-frequency acoustic waves is determined by expression

¿inst

2T0lCv

ToWtо - poWpo - ToY

(1.4)

The estimations of instability time for the chosen parameters of the mediums give tinst = 1065 years. Thus, the growth from the half of the amplitude of initial perturbation (since the waves propagate in two opposite directions) to a half of maximum amplitude may take about 4000 years, which agrees well with the results of numerical simulations on the finest grid.

However, the subsequent growth of the waves from half to maximum amplitude takes about 15-20 thousand years with little dependence on grid size and the characteristic size of initial perturbation.

Conclusion

Growth time of acoustic perturbations for parameters of isentropically unstable photodissociation region Orion Bar is estimated. The growth of acoustic waves can be divided into 3 stages. At the first stage, the high-frequency components of the initial perturbation grow without a significant change in the amplitude of the wave packet as a whole. The duration of this stage depends substantially on the characteristic size of the initial perturbation and can occupy fractions of the characteristic instability time tinst for high-frequency perturbations, and take tens of tinst for low-frequency ones, which corresponds to several tens of thousands of years for Orion Bar. In the second stage, there is an explosive growth of the wave amplitude up to half of the maximum value. This process takes several tinst or about 4000 years for Orion Bar for initial perturbation with an amplitude of 0.01 of the equilibrium concentration. In the third stage, there is a smooth increase in the amplitude of the wave up to the maximum value, which takes about 15-20 thousand years.

This work was supported in part by the Ministry of Education and Science (projects FSSS-2020-0014, 0023-2019-0003).

References

[1] Molevich N.E., Zavershinsky D.I., Galimov R.N., Makaryan V.G. Travelling self-sustained structures in interstellar clouds with the isentropic instability. Astrophysics and Space Science, 2011, vol. 334, no. 1, pp. 35-44. DOI: http://doi.org/10.1007/s10509-011-0683-0.

[2] Krasnobaev K.V., Tagirova R.R., Arafailov S.I., Kotova G.Y. Evolution and saturation of Autowaves in photodissociation regions. Astronomy Letters, 2016, vol. 42, no. 7, pp. 460-473. DOI: http://doi.org/10.1134/S1063773716070057.

[3] Krasnobaev K.V., Tagirova R.R. Isentropic thermal instability in atomic surface layers of photodissociation regions. Monthly Notices of the Royal Astronomical Society, 2017, vol. 469, issue 2, pp. 1403-1413. DOI: http://dx.doi.org/10.1093/mnras/stx884.

[4] Zavershinskii D.I., Kolotkov D.Y., Nakariakov V.M., Molevich N.E., Ryashchikov D.S. Formation of quasi-periodic slow magnetoacoustic wave trains by the heating/cooling misbalance. Physics of Plasmas, 2019, vol. 26, issue 8, p. 82113. DOI: http://doi.org/10.1063/1.5115224.

[5] Zavershinskii D.I., Molevich N.E., Riashchikov D.S., Belov S.A. Nonlinear magnetoacoustic waves in plasma with isentropic thermal instability. Physical Review E, 2020, vol. 101, issue 4, p. 43204. DOI: http://doi.org/10.1103/PhysRevE.101.043204.

[6] Molevich N., Riashchikov D. Shock wave structures in an isentropically unstable heat-releasing gas. Physics of Fluids, 2021, vol. 33, issue 7, p. 076110. DOI: http://doi.org/10.1063/5.0053394.

[7] Zavagno A., Pomares G., Deharveng L., Hosokawa T., Russeil D., Caplan J. Triggered star formation on the borders of the Galactic H ii region RCW 120. Astronomy & Astrophysics, 2007, vol. 472, issue 3, pp. 835-846. DOI: http://dx.doi.org/10.1051/0004-6361:20077474.

[8] Deharveng L., Zavagno A., Schuller F. Caplan J., Pomares G., C. De Breuck. Star formation around RCW 120, the perfect bubble. Astronomy & Astrophysics, 2009, vol. 496, issue 1, pp. 177-190. DOI: http://dx.doi.org/10.1051/0004-6361/200811337.

[9] Goicoechea J.R. [et al.] Compression and ablation of the photo-irradiated molecular cloud the Orion Bar. Nature, 2016, vol. 537, issue 7619, pp. 207-209. DOI: http://doi.org/10.1038/nature18957.

[10] Berne O., Marcelino N., Cernicharo J. Waves on the surface of the Orion molecular cloud. Nature, 2010, vol. 466, issue 7309, pp. 947-949. DOI: http://doi.org/10.1038/nature09289

[11] Field G.B. Thermal Instability. The Astrophysical Journal, 1965, Vol. 142, pp. 531-567.

[12] Zavershinskii D., Kolotkov D., Riashchikov D., Molevich N. Mixed Properties of Slow Magnetoacoustic and Entropy Waves in a Plasma with Heating/Cooling Misbalance. Solar Physics, 2021, vol. 296, issue 6, p. 96. DOI: http://doi.org/10.1007/s11207-021-01841-1.

[13] Wakelam V., Herbst E. Polycyclic Aromatic Hydrocarbons in Dense Cloud Chemistry. The Astrophysical Journal, 2008, vol. 680, issue 1, pp. 371-383. DOI: http://dx.doi.org/10.1086/587734

[14] Sone J.M., Gardiner T.A., Teben P., Hawley J.F., Simon J.B. Athena: a new code for astrophysical MHD. The Astrophysical Journal Supplement Series, 2008, vol. 178, issue 1, p. 137. DOI: http://dx.doi.org/10.1086/588755.

Научная статья DOI: 10.18287/2541-7525-2028-28-1-2-113-119

УДК 534.222.2; 524.5 Дата: поступления статьи: 06.09.2022

после рецензирования: 12.10.2022 принятия статьи: 14.11.2022

Д.С. Рящиков

СФ ФИАН, г. Самара, Российская Федерация; Самарский национальный исследовательский университет имени академика С.П. Королева, г. Самара, Российская Федерация E-mail: ryashchikovd@gmail.com. ORCID: https://orcid.org/0000-0001-7143-2968

И.А., Помельников Самарский национальный исследовательский университет имени академика С.П. Королева, г. Самара, Российская Федерация E-mail: vanidzepomelnikov@gmail.com. ORCID: https://orcid.org/0000-0001-7839-5784

Н.Е. Молевич СФ ФИАН, г. Самара, Российская Федерация; Самарский национальный исследовательский университет имени академика С.П. Королева, г. Самара, Российская Федерация E-mail: nonna.molevich@mail.ru. ORCID: https://orcid.org/0000-0001-5950-5394

ВРЕМЯ РОСТА АКУСТИЧЕСКИХ ВОЗМУЩЕНИЙ В ИЗОЭНТРОПИЧЕСКИ НЕУСТОЙЧИВОЙ ТЕПЛОВЫДЕЛЯЮЩЕЙ СРЕДЕ2

АННОТАЦИЯ

Изоэнтропическая неустойчивость является одним из типов тепловой неустойчивости, которая приводит к росту акустических волн. В результате их роста в таких средах образуются автоволновые структуры, параметры которых зависят только от свойств среды и могут быть предсказаны как аналитически, так и численно. Целью данного исследования является определения времени формирования автоволновых структур в изоэнтропической неустойчивой среде с параметрами, характерными для области фотодиссоциации Орион Бар. Показано, что время роста зависит от характерного размера начального возмущения. Наиболее быстро растущие структуры достигают половины от максимальной амплитуды за 3-6 тысяч лет. Дальнейший рост до максимального значения занимает 15-20 тысяч лет.

Ключевые слова: неустойчивость; тепловая неустойчивость; нелинейные волны; ударные волны; автоволны; межзвездный газ; области фотодиссоциации; Орион Бар.

Цитирование. Riashchikov D.S., Pomelnikov I.A., Molevich N.E. Growth time of acoustic perturbations in isentropically unstable heat-releasing medium // Вестник Самарского университета. Естественнонаучная серия. 2022. Т. 28, № 1-2. С. 113-119. DOI: http://doi.org/10.18287/2541-7525-2022-28-1-2-113-119.

Информация о конфликте интересов: авторы и рецензенты заявляют об отсутствии конфликта интересов.

© Рящиков Д.С., 2022

Дмитрий Сергеевич Рящиков — кандидат физико-математических наук, научный сотрудник, теоретический сектор, Самарский филиал Физического института им. П.Н. Лебедева РАН, 443011, Российская Федерация, г. Самара, ул. Ново-Садовая, 221.

Старший преподаватель кафедры физики, Самарский национальный исследовательский университет имени академика С.П. Королева, 443086, Российская Федерация, г. Самара, Московское шоссе, 34.

©c Помельников И.А., 2022 Иван Александрович Помельников — студент института информатики и кибернетики, Самарский национальный исследовательский университет имени академика С.П. Королева, 443086, Российская Федерация, г. Самара, Московское шоссе, 34.

2Работа выполнена при поддержке Минобрнауки РФ в рамках государственного задания вузам и научным организациям (проекты Р888-2020-0014, 0023-2019-0003).

© Молевич Н.Е., 2022

Нонна Евгеньевна Молевич — доктор физико-математических наук, главный научный сотрудник, и.о. зав. теоретическим сектором, Самарский филиал Физического института им. П.Н. Лебедева РАН, 443011, Российская Федерация, г. Самара, ул. Ново-Садовая, 221.

Профессор кафедры физики, Самарский национальный исследовательский университет имени академика С.П. Королева, 443086, Российская Федерация, г. Самара, Московское шоссе, 34.

Литература

[1] Molevich N.E., Zavershinsky D.I., Galimov R.N., Makaryan V.G. Travelling self-sustained structures in interstellar clouds with the isentropic instability // Astrophysics and Space Science, 2011, Vol. 334, issue 1, pp. 35-44. DOI: http://doi.org/10.1007/s10509-011-0683-0.

[2] Krasnobaev K.V., Tagirova R.R., Arafailov S.I., Kotova G.Y. Evolution and saturation of Autowaves in photodissociation regions // Astronomy Letters, 2016, vol. 42, issue 7, pp. 460-473. DOI: http://doi.org/10.1134/S1063773716070057.

[3] Krasnobaev K.V., Tagirova R.R. Isentropic thermal instability in atomic surface layers of photodissociation regions // Monthly Notices of the Royal Astronomical Society, 2017, vol. 469, issue 2, pp. 1403-1413. DOI: http://dx.doi.org/10.1093/mnras/stx884.

[4] Zavershinskii D.I., Kolotkov D.Y., Nakariakov V.M., Molevich N.E., Ryashchikov D.S. Formation of quasi-periodic slow magnetoacoustic wave trains by the heating/cooling misbalance // Physics of Plasmas, 2019, vol. 26, no. 8, p. 82113. DOI: http://doi.org/10.1063/1.5115224.

[5] Zavershinskii D.I., Molevich N.E., Riashchikov D.S., Belov S.A. Nonlinear magnetoacoustic waves in plasma with isentropic thermal instability // Physical Review E, 2020, vol. 101, issue 4, p. 43204. DOI: http://doi.org/10.1103/PhysRevE.101.043204.

[6] Molevich N., Riashchikov D. Shock wave structures in an isentropically unstable heat-releasing gas // Physics of Fluids, 2021, vol. 33, issue 7, p. 076110. DOI: https://doi.org/10.1063/5.0053394.

[7] Zavagno A., Pomares G., Deharveng L., Hosokawa T., Russeil D., Caplan J. Triggered star formation on the borders of the Galactic H ii region RCW 120 // Astronomy & Astrophysics, 2007, vol. 472, issue 3, pp. 835-846. DOI: http://dx.doi.org/10.1051/0004-6361:20077474.

[8] Deharveng L., Zavagno A., Schuller F. Caplan J., Pomares G., C. De Breuck. Star formation around RCW 120, the perfect bubble // Astronomy & Astrophysics, 2009, vol. 496, issue 1, pp. 177-190. DOI: http://dx.doi.org/10.1051/0004-6361/200811337

[9] Goicoechea J.R. [et al.] Compression and ablation of the photo-irradiated molecular cloud the Orion Bar // Nature, 2016, vol. 537, issue 7619, pp. 207-209. DOI: http://doi.org/10.1038/nature18957.

[10] Berne O., Marcelino N., Cernicharo J. Waves on the surface of the Orion molecular cloud // Nature, 2010, vol. 466, no. 7309, pp. 947-949. DOI: http://doi.org/10.1038/nature09289.

[11] Field G.B. Thermal Instability // The Astrophysical Journal, 1965, vol. 142, pp. 531-567.

[12] Zavershinskii D., Kolotkov D., Riashchikov D., Molevich N. Mixed Properties of Slow Magnetoacoustic and Entropy Waves in a Plasma with Heating/Cooling Misbalance // Solar Physics, 2021, Vol. 296, no. 6, p. 96. DOI: http://doi.org/10.1007/s11207-021-01841-1.

[13] Wakelam V., Herbst E. Polycyclic Aromatic Hydrocarbons in Dense Cloud Chemistry // The Astrophysical Journal, 2008, vol. 680, no. 1, pp. 371-383. DOI: http://dx.doi.org/10.1086/587734.

[14] Sone J.M., Gardiner T.A., Teben P., Hawley J.F., Simon J.B. Athena: a new code for astrophysical MHD // The Astrophysical Journal Supplement Series, 2008, vol. 178, no. 1, p. 137. DOI: https://doi.org/10.1086/588755.