Об одном методе приближенного решения нелинейного уравнения теплопроводности

ЧЕБЫШЕВСКИЙ СБОРНИК

Том 17. Выпуск 4.

УДК 539.3 DOI 10.22405/2226-8383-2016-17-4-11-22

ОБ ОДНОМ МЕТОДЕ ПРИБЛИЖЕННОГО РЕШЕНИЯ НЕЛИНЕЙНОГО УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

В. Б. Беднова

Аннотация

В данной работе рассматривается одномерная нестационарная задача теплопроводности, моделирующая процесс быстрого локального нагрева образца балочного типа по боковой поверхности. При этом характер нагрева таков, что можно выделить единственное определенное направление, в котором распространяется тепло.

Температурные поля определяются приближенным методом, основанным на идее теплового фронта. Решение ищется в виде степенного ряда по координате с коэффициентами, зависящими от времени. Границы фронта распространения тепла как функции времени определяются из условия интегрального удовлетворения уравнению теплопроводности.

Рассматриваемые температурные поля возникают во многих технологических процессах, например, при лазерной обработке материалов, когда из-за больших градиентов температур могут возникать температурные напряжения, приводящие к микрорастрескиванию внутренних слоев или разрушению элементов конструкций. Аналитический вид решения задачи теплопроводности позволяет получить аналитические выражения для температурных напряжений и в дальнейшем облегчает анализ результатов.

В работе получены решения задач с граничными условиями первого и второго родов для двух монотонных и одной немонотонной зависимостей коэффициента теплопроводности от температуры. Проведено сравнение полученного решения нестационарной линейной задачи с точным и показана приемлемость метода для дальнейшего использования.

Ключевые слова: физическая нелинейность уравнения теплопроводности, приближенное решение уравнения теплопроводности, тепловой фронт.

Библиография: 16 названий.

ON A METHOD FOR APPROXIMATE SOLUTION NONLINEAR HEAT CONDUCTION EQUATION

V. B. Bednova (Moscow) Abstract

In this paper a one-dimensional non-stationary heat conduction problem, modeling the process of rapid local heating of the sample beam type on the lateral surface is considered. The character of heating is such that it is possible to allocate only a certain directionof the heat propagation.

Temperature fields are determined by an approximate method based on the idea of the thermal front. The solution is sought in the form of a power series in the coordinate with coefficients depending on time. The boundaries of the front heat distribution as a function of time are determined by the condition of the integral satisfaction of the heat conduction equation.

Considered temperature fields arise in many industrial processes, such as laser material processing, when due to large temperature gradients can arise thermal stresses, leading to microcracking inner layers or the destruction of structural elements. Analytical view of the heat conduction problem's solution allows to obtain analytical expressions for the thermal stresses and further facilitates the results analysis.

The paper presents the solution of problems with boundary conditions of the first and second kinds for two monotonic and one non-monotonic dependencies of the thermal conductivity coefficient on temperature. The approximate solution and the exact solution of the non-stationary linear problem are compared and shows the suitability of the method for future use.

Keywords: Physical non-linearity of the heat conduction equation, approximate solution of the heat conduction equation, thermal front.

Bibliography: 16 titles.

1. Введение

Известно, что классическое уравнение притока тепла при описании процессов теплопередачи механизмом теплопроводности дает не соответствующую действительности бесконечную скорость распространения возмущений. Кроме того, во многих прикладных задачах граничные условия невозможно задать с высокой степенью точности. Погрешность вполне может достигать 10 — 20%. Поэтому в реальных процессах актуальным может являться построение приближенных решений, удовлетворяющих некоторым интегральным энергетическим условиям.

В работе предлагается приближенный метод решения нелинейного уравнения теплопроводности с различными зависимостями коэффициента теплопроводности от температуры. Метод основывается на работах [1, 9, 14], развитых для линейного уравнения теплопроводности. В работе [15] с помощью приближенного решения уравнения теплопроводности показана возможность возникновения зон макроразрушения при локальном импульсном нагреве диска.

Рассматривается обезразмеренное уравнение теплопроводности

Тт = (К (Т )ТХ)Х,

где Т, т, х — соответственно безразмерные температура, время, координата, К(Т) — безразмерный коэффициент теплопроводности. Обезразмеривание задачи теплопроводности проведено в работе [8], зависимость размерной температуры от безразмерной имеет вид Т = > , где _ температура среды, вт — температура плавления материала. Во всех задачах начальная (лабораторная) температура Т|¿=о равна нулю. Также при решении используется понятие границы теплового фронта 1(т) [14, 15, 16], которую можно интерпретировать как линию уровня пренебрежимо малой температуры, функция 1(т) определяется из условия интегрального удовлетворения уравнению теплопроводности.

2. Построение приближенных решений линейной задачи теплопроводности

Рассмотрим линейную задачу (К(Т) = 1), моделирующую процесс нагрева балки по границе ж = 0 с постоянной по времени температурой То = 1, причем на границе теплового фронта 1(т) температура Т(т, х) и поток Тх(т,х) считаются равными нулю:

Т = Т

±Т - ±хх

Т\х=1(т) = 0 Тх\х=1(т) =

Перейдем от переменных (т, ж) к переменным (t, z) следующим образом:

г = ж - 1

£ = т.

Преобразования производных в этом случае:

д д д2 д2 д д , д дх дг дх2 дг21 дт сЯ дг'

В новых переменных задача имеет вид

Тг - V(г)Т2 = тхх

т\*=-т = 1

т\,=0 = 0

Т.г и =

Точные решения уравнения теплопроводности хорошо аппроксимируются степенными функциями или многочленами, содержащими три-четыре слагаемых [5]. Поэтому приближенные решения будем искать в виде степенных рядов по координате х, ограниченных сначала до трех, потом — до четырех слагаемых. Сначала возьмем

т ^ 2) = ( ио(Ь) + щ(г)г + и2&) ^, если -1(Ь) <г< 0; ^

0

Из граничных условий задачи (2) при учете (3) получаем ио = 0; и\ = 0; и2 = откуда 2

Т(г, г) = тщЩ>и -Щ) < г < 0.

Для нахождения 1(Ь) подставим данное решение в уравнение из (2), проинтегрируем полученное выражение в пределах от -1(Ь) до 0 и приравняем к нулю:

о

/ (-%

'(г) 2 21'(г) 2 \ 21'(г) 1(г)3 21'(г) 12(г) 2

. Г13 " " " '

-1

,2 - г - - аг = + - т^т= 0-

13

Получим дифференциальное уравнение на 1(1): II' = Ь, откуца ^ = Ы + с. Так как 1(0) = 0, то с = 0 и 1(г) = 2^31.

2

В итоге, когда ряд ограничен тремя членами, в прогретой зоне получаем Т(I, г) = или

г{т,х) = к-ур2.

Теперь будем искать решение в виде

{2 3

ио(^ + иг(г)г + щ^)\ + и3(г), если -1(г) < г < 0 ^

0

В данном случае у нас четыре неизвестных функции щ, г = 0,1, 2, 3. При тех же граничных условиях получаем

I2 ш 13а)

ио = 0; ил = 0; и,2—--= 1.

2 Ь

Четвертое условие получим, продифференцировав уравнение теплопроводности задачи (2) по г и положив г — 0 [2]:

откуда при учете (4) и (5) —21'{Ъ)и2 — 6и3.

Таким образом и0 — 0, и\ — 0 и2 — р^+и'у и3 — — ¡2(3+ и решение (2) при —I(¿) < г < 0

имеет вид

т (' *) — 12(3 + 11') ** 12(3 + 11')г 3 Аналогично предыдущему случаю получаем выражение для подынтегральной функции

_ 6___31' 3121'3 — 18 I' — 31211" 2 _ 3 II'' — 311'3 — 61 '2 3

( ) — —12(3 + 11') — 12(3 + 11')* + 13(3 + 11')2 * I3(3 + 11')2 23

интегрируя которое по г от —I(í) до 0 и приравнивая к 0 получаем уравнение на 1(Ь):

13113 — 131 '' — 30 II' — 72 — 0.

Данное дифференциальное уравнение можно решить, сначала понизив порядок дифференцирования с помощью замены I' — р(1), I" — р'(1)р(1). Получим дифференциальное уравнение 13рр' — 13р3 + 30 1р + 72 — 0. Решение этого дифференциального уравнения получаем в неявном виде: 1пI — с\ — п(р1 + 3) — 151п(р1 — 6) + 511п(р1 + 4) — 0. Данный вид позволяет получить достаточно точную зависимость функции теплового фронта от температуры.

Также для нахождения /(¿), удовлетворяющей (6), можно предположить, что 1(1 ) — кл/1. Тогда, подставив 1(1) в (6), получим уравнение на к: к6 + 2к4 — 120к2 — 576 — 0, откуда к — 2\/3 и 1(Ъ) — 2^3. Тогда

11

3

т(г, г) — — х2--з-г

36 36^312

или, перейдя к переменным

Т

36т 36л/3т 2

Проведем сравнение полученных приближенных решений с точным решением аналогичной задачи для полуограниченного твердого тела [8]. Данное сравнение корректно, потому что нас интересуют промежутки времени, когда тело еще не прогрелось полностью и температура на второй границе практически не изменилась. В этом случае применимо решение для полуограниченного твердого тела:

х

2 —

т

о

Сравнение приближенных решений с точным показано на рис. 1. Из графика видно, что приближенные и точное решения отличаются незначительно.

Рис. 1: Безразмерная температура в моменты безразмерного времени т = у^д — 1, т = 10 — 2. Сплошная .линия точное решение (8), короткий пунктир приближенное решение (7), длинный пунктир приближенное решение (6).

3. Построение приближенных решений задач теплопроводности с коэффициентом К (Т) = Т

1. В работе [2| доказывается существование и единственность аналитического решения задачи

( Тт = (К{Т)ТХ)Х

\ Т|Х=1(т) = 0,

описывающего распространение тепловой волны по холодному фону. Коэффициент теплопроводности берется в виде К(Т) = аТа, а и а — константы.

Решение задачи (9) строится в виде бесконечного ряда по координате. Найдем первые четыре члена ряда по методу, предложенному в [2] в случае К(Т) = Т при а = 1 и а = 1.

Задача теплопроводности (9) после замены переменных (1) будет иметь вид

( Тг - 1'Тг = ТТХХ + Т2 (

\ Т |,=0 = 0

Температуру ищем в виде ряда (4). Из граничного условия следует, что ио = 0. Далее, уравнение задачи (10) при г = 0 и учете найденного ио = 0 переходит в соотношение и,\ = —V. Для нахождения щ н и3 уравнение задачи (10) дифференцируется соответственно один и два раза по г, полагается г = 0 и, при учете уже пай денных щ, получаем соотнош ения и'у = = (3и\ + V)«2 и и'2 — I1 и3 = 4и\и3 + 3и22. После преобразований получаем

I" Ы"2 — 21'1Ш

иг = —1 ,и2 = -,щ = Шз .

Таким образом функция температуры по методу, описанному в работе [2] в зоне —I (t) < z < 0 имеет вид

^ . „ I" z2 I''2 — 21'I'''z3

Т

решенной выше. Если оставить только первые два слагаемых в (11), то есть взять

I''

Т(t, z) = —I'z +—— -— v ' ! 21' 2

и перейти к координатам (т,х), получим

2 2

Т(т,х)=11 + — —1 х — —х + —х •

Проинтегрировав выражение L(T) = Тт —Т2 — ТТХХ от 0 до ¿(т) по координате х и приравняв результат к нулю, получим дифференциальное уравнение 214"' — Ы"2 = 0. Если в него подставить функцию I = кта, получим а = | и 1(т) = кт1/3, откуда

к2 х2 ( ' ) 6т1/3 6т

Такое решение может описать ситуацию с изначально неравномерно прогретым телом, когда процесс состоит только в выравнивании температуры по всему образцу.

2. Теперь рассмотрим задачу о нагреве с заданным постоянным тепловым потоком qoo на нагреваемой поверхности с предположением, что на границе теплового фронта температура равна нулю:

Т = Т2 + ТТ

Тт - Т Х + Т ТХХ

Т\х=1(т) = 0

—ТТх\Х=О = qoo• Перейдя к координатам (t, z), получим

Тг — I 'Т2 = ТТгг +Т2

Т\z=o = 0

—ТТх \z=-l = Qoo

Температуру будем искать в виде ряда (3). В качестве третьего условия берем выражение Т — I'Т\'z = TTzz + Т2 при z = 0. Получим:

и0 = 0 (ml — U2l~2 )(ui —U2I) = qoo

U\ = 0 U\ = —I'

При U\ = 0 получим Т(t' z) = z2, или Т(т, х) = \JЩг (1 — Х)2-

При U\ = —l' получаем

rr,. \ ' , —3 IV + Vl2l'2 + Sqool Z2

(t' Z^ = —U +-22-~2 •

Уравнение на функцию теплового фронта можно получить из условия

qoot = J Т dz, -i

откуда

8 14l12 - 72qootl2l' = 8 qQQl3 - 144qoot2.

Второй вариант для нахождения l(t) — интегрально удовлетворить уравнению теплопроводности

o

J (Т - l'Tz - Т2 - TTZZ) dz = 0.

-l

В этом случае получаем

131Ч" + 3/2l Wl2l'2 + 8 qool + 2 12l'3 + 6 11,2^l2l12 + 8 q ool + 12 qool I1--12 q ooVl2l'2 + 8 q ool = 0.

4. Построение приближенных решений задач теплопроводности с коэффициентом К (Т) = 1 - Т

1. Рассматривается задача с коэффициентом теплопроводности К(Т) = 1 - Т. На нагреваемой поверхности температура считается равной единице, на границе теплового фронта — нулю:

Тт - Тх + Тхх ТТХХ

Т\х=1(г) = 0

Пх^ = 1

Сделаем замены переменных и = 1 - Т и (1), получим

Ut - I'(t)uz = uuzz + и2

и

u\z=-l(t) =

Ищем решение в виде

u(z, t) =

{uo(t) + u,\(t)z + u2(t)z2, если -l(t) < z < 0; 02

2

иначе,

где необходимо определить три неизвестных функции времени. В качестве третьего условия будет уравнение системы (13), в котором положим г = 0. Из граничных условий получаем

ио = 1 _

= л/(2+И' )2 +8-11'-2 и'1 = 21

л/(2+Й7)2+8-II'

и2 = --2Г2-.

o

Таким образом

т , г) = -11'-2г — г и- (^ <г< 0

[О иначе.

2. Теперь рассмотрим случай, когда на поверхности нагрева задан постоянный начальный поток доо. Предполагается, что на границе теплового фронта температура и тепловой поток равны нулю. После замены переменных аналогично (12) получим

и — I '{Ъ)их = иигг + и2 ииг \г=-1 = доо иг\г=о = О и| г=0 = 1.

При данных граничных условиях получаем ио = 1, и\ = 0, и2 = Функция температуры после возврата к переменным (Т, т, х) будет иметь вид

т (тх) = { 1-^2-2Я°о1 (х — 1)2, если0 <х< 1(т)-, О

Функция теплового фронта ищется из условия

1

I (Тт — т2 + (1 — т) Тхх) йх = О,

откуда

Решение дифференциального уравнения (14) имеет вид

13 — 6т I + 18 доот2 = 0.

о

5. Построение приближенного решения нелинейного уравнения теплопроводности с коэффициентом К (Т) = 2Т

1+Т2

1. Рассмотрим задачу

{

Тт = (К (Т)ТХ)Х

Т (I ,х)\х=Кт) = 0

с условием К(Т)\х=1(т) = 0.

После замены переменных (1) получаем

{

Строим решение в виде ряда

Т = 1'(1)ТХ + К' (т)Т2 + К (Т)ТХ,

Т(I, г)\2=о = 0.

Т (I, г) = ¡То+Т1(1)г + Т2(1) ^, если —I (г) <г< 0 ^

1 0 иначе,

где по условию То(Ь) = 0, Т0(Ь) = 0. В уравнение в системе (16) подставляем г = 0. Получаем

1'(г)Т1 + К' (Т)Т2 = 0,

откуда

* Т т 1'

Т  = - КЩ, Т  =-КЩТ-)-.

Дифференцируем уравнение в системе (16) по г:

Тг = I 'Тхх + 2К'ТХТХХ + К'Т + КТХХХ + К'ТХ Тхх. Подставляем х = 0 ш (17):

Т[ = 1'Т2 + 3К'Т\Т2 + К''Т—.

Подставляем в (19) Т' (¿) и Т\ (1) из (18). Получаем

^ К'21'' - К'' I'3 - К'К'П'

Т2 = -

21'К'3

и в итоге при -I(¿) < г < 0

^ , 1'(г) К'21''-К''I'3 -К'К'П' г2 , ,

Выражение (20) дает общий вид зависимости температуры от времени и координаты для коэффициента теплопроводности, удовлетворяющего условию К(Т)1Х=1(Т) = 0.

2. Теперь решим задачу теплопроводности (15) с коэффициентом К(Т) = 2, опираясь на (20). Поскольку К'(Т) = 2 5 система (16) примет вид

= 1,($)Тг + ((+у2)2 Тг + 1++Т2

Т^, г)12=о = 0.

В этом случае получаем Т\ = -,Т2 = то есть Т = - ^ + \2

=1

уравнение

Т = 1 то есть выполнено соотношение Т — Ч7 + ^¡т — 1, откуда получаем дифференциальное

411'2 - 8I' + 121'' = 0.

3 в2 в2 82\/3 (2 1 \

16!2 - + + ^ а21 + *2) + — «**{ ^-73з)+с1 =0

где в = (3С2)3 , а константы С, С связаны зависимостью из условия 1(0) = 0. Также можно найти вид функции границы теплового фронта, подставив в (21) 1(1) = кл/1. Тогда к = — и

Т , г) = 1- 7Шг - ш ^ если- — <*< 0 0

6. Заключение

В настоящей работе предложены приближенные решения одномерного уравнения теплопроводности с граничными условиями первого и второго родов. В нелинейном случае решения получены для двух монотонных зависимостей коэффициента теплопроводности от температуры и одной немонотонной. В линейном случае проведено сравнение построенного решения с точным, из которого видно, что точность приближенного метода не превышает точности задания физических констант и граничных условий в прикладных задачах. Аналитический вид полученных решений облегчает анализ результатов и позволяет получить аналитические выражения для температурных напряжений.

Автор приносит благодарность М. В. Юмашеву за постановку задачи и ценные обсуждения.

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. Баренблатт Г. И. О некоторых приближенных методах в теории одномерной неустановившейся фильтрации жидкости при упругом режиме. Изв.АН СССР, ОТН. 1954. №9. С. 35 19.

2. Баутин С. П. Аналитическая тепловая волна. - М.: Физматлит, 2003. - 88 с.

3. Бахарев М. С., Миркин Л. И., Шестериков С. А., Юмашева М. А. Структура и прочность материалов при лазерных воздействиях. - М.: Изд-во Моск. ун-та, 1988. - 224 с.

4. Беднова В. Б. Приближенный метод определения температурного поля при быстром локальном нагреве образца // Труды конференции-конкурса молодых ученых НИИ механики МГУ имени М. В. Ломоносова. 2013. С. 73-76.

5. Боли Б., Уэйнер Дж. Теория температурных напряжений. - М.: Мир, 1964. - 517 с.

6. Галактионов В. А., Курдюмов С. П., Михайлов А. П., Самарский А. А. Локализация тепла в нелинейных средах // Дифференциальные уравнения. Октябрь 1981. Т. XVII. №10. С. 1826-1841.

7. Калашников А. С. Об уравнениях типа нестационарной фильтрации с бесконечной скоростью распространения возмущений // Вестник Московского университета. 1972. №6. С. 45-49.

8. Карслоу Г., Егер Д. Теплопроводность твердых тел. - М.: Наука, 1964. - 488 с.

9. Локощенко А. М. Ползучесть и длительная прочность металлов в агрессивных средах. -М.: Изд-во Моск. ун-та, 2000. - 178 с.

10. Лыков А. В. Теория теплопроводности. - М.: Высшая школа, 1967. - 599 с.

11. Олейник О. А., Калашников А. С., Чжоу Юй-Линь Задача Коши и краевые задачи для уравнений типа нестационарной фильтрации // Известия АН СССР. Серия математическая. 1958. Т. 22. С. 667-704.

12. Паркус Г. Неустановившиеся температурные напряжения. - М.: Государственное изд-во физико-математической литературы, 1963. - 252с.

13. Тимошенко С. П., Гудьер Дж. Теория упругости. - М.: Наука, 1975. - 576 с.

14. Шестериков С. А., Юмашева М. А. Приближенный метод оценки нестационарных температурных полей // Институт механики МГУ. Научные труды. Деформирование и разрушение твердых тел. Вып. 23. - М.: Изд-во МГУ, 1973. С. 15—20.

15. Юмашев М. В., Беднова В. Б., Вергазов М. \!.. Юмашева М. А. Разрушение хрупких материалов в условиях локального воздействия на поверхность энергетическим потоком // Машиностроение и инженерное образование. 2014. №4. С. 52^58.

16. Юмашев М. В., Юмашева М. А., Краснова П. А. Моделирование процесса нагрева тела при интенсивном тепловом воздействии на поверхность // Вестник Московского университета. 2010. №4. С. 44^54.

REFERENCES

1. Barenblatt, G. I. 1954, "On some approximate methods in the theory of one-dimensional non-stationary filtration in the elastic drive regime", Izv. AN SSSR, OTN, no. 9, pp. 35 19.

2. Bautin, S. P. 2003, "Analiticheskava teplovava volna", [Analytical heat wave], Fizmatlit, Moscow, 88p.

3. Baharev, M. S., Mirkin, L. I., Shesterikov, S. A. к Yumasheva M. A. 1988, "Struktura i prochnost' materialov pri lazernvh vozdejstvivah", [Structure and strength of materials with the laser action], Izd-vo Mosk. un-ta, Moscow 224p.

4. Bednova, V. B. 2013, "An approximate method for determining the temperature field with rapid local heating of the sample", Trudy konferencii-konkursa molodyh uchenyh N11 mekhaniki MGU imeni M. V.Lomonosova, pp. 73-76.

5. Bolev, В. A. к Weiner J. H. 1964, "Teoriva temperaturnvh naprvazhenij", [Theory of Thermal Stresses], Mir, Moscow, 517p.

6. Galaktionov, V. A., Kurdvumov, S. P., Mihajlov, A. P. к Samarskij, A. A. 1981, "Heat localization in nonlinear media", Differencial'nye uravneniya, Vol. XVII, no. 10, pp. 1826-1841.

7. Kalashnikov, A. S. 1972, "On the equations of non-stationary filtration type with infinite speed of propagation of disturbances", Vestnik Moskovskogo universiteta, no. 6, pp. 45-49.

8. Carslaw, H.S. к Jaeger, J. C. 1964, "Teploprovodnost' tverdvh tel", [Conduction of Heat in Solids], Nauka, Moscow, 488p.

9. Lokoshchenko, A. M. 2000, "Polzuchest' i dlitel'nava prochnost' metallov v agressivnvh sredah", [The creep and the long-term strength of metals in the aggressive media], Izd-vo Mosk. un-ta, Moscow, 178p.

10. Lvkov, A. V. 1967, "Teoriva teploprovodnosti", [Theory of heat conduction], Vvsshava shkola, Moscow, 599p.

11. Olejnik, O. A., Kalashnikov, A. S. к Chzhou Yuj-Lin' 1958, "The Cauchv problem and boundary value problems for equations of non-stationary filtration", Izv. AN SSSR, Seriya matematicheskaya, Vol. 22, pp. 667-704.

12. Parcus, H. 1963, "Neustanovivshiesva temperaturnve naprvazheniva", [Non-stationary thermal stresses], Gosudarstvennoe izd-vo fiziko-matematicheskoj literaturv, Moscow, 252p.

13. Timoshenko, S. Р. & Goodier, J. 1975, "Teoriva uprugosti", [Theory of Elasticity], Nauka, Moscow, 576p.

14. Shesterikov, S. A. k, Yumasheva, M. A. 1973, "An approximate method for estimating non-stationary temperature fields", Institut mekhaniki MGU. Nauchnye trudy. Deformirovanie i razrushenie tverdyh tel. no. 23, Izd-vo Mosk. un-ta, Moscow, pp. 15 21).

15. Yumashev, M. V., Bednova, V. В., Vergazov, M. M. k, Yumasheva, M. A. 2014, "The destruction of brittle materials under local effect on the surface energy flow", Mashinostroenie i inzhenernoe obrazovanie, no. 4, pp. 52^58.

16. Yumashev, M. V., Yumasheva, M. A. k, Krasnova, P. A. 2010, "Modeling of the body heating process with intensive exposure to heat to the surface", Vestnik Moskovskogo universiteta, no. 4, pp. 44^54.

Московский государственный университет имени М. В. Ломоносова

Получено 25.07.2016 г.

Принято в печать 13.12.2016 г.