Интегрируемые модели динамики сжимаемой среды в собственном поле тяготения. Метод подстановок Коула – Хопфа Текст научной статьи по специальности «Математика»

УДК 532.5, 52-43, 514.86

В. М. Журавлев, Д. А. Зиновьев

ИНТЕГРИРУЕМЫЕ МОДЕЛИ ДИНАМИКИ СЖИМАЕМОЙ СРЕДЫ В СОБСТВЕННОМ ПОЛЕ ТЯГОТЕНИЯ. МЕТОД ПОДСТАНОВОК КОУЛА - ХОПФА1

Аннотация. На основе модификации метода обобщенных подстановок Коула -Хопфа строятся точные решения самосогласованной задачи динамики сжимаемой среды в собственном поле тяготения. Рассматриваются различные модели одномерного течения идеальной и вязкой самогравитирующей среды для общего вида начальных распределений ее плотности и скорости. Вычисляются интегралы движения и обобщенные подстановки Коула - Хопфа, позволяющие строить точные решения. Метод построения точных решений обобщается для идеальной среды в собственном поле тяготения на двумерные и трехмерные течения, обладающие цилиндрической и сферической симметриями. На основе построенных решений анализируются процессы формирования плотных объектов из пылевого облака с заданным начальным распределением плотности пыли вследствие неустойчивости Релея - Джинса для плоской и сферической симметрий. Вычисляется характерное время до образования сингулярности.

Ключевые слова: нелинейные уравнения, самогравитирующие среды, точно интегрируемые модели.

Abstract. The exact solutions of the self-coordinated problem of dynamics of the compressible fluid in self-gravitational field are developed by the authors on the basis of modifications to the Cole - Hopf generalized substitutions method. The article considers various models of a one-dimensional flow of the ideal and viscous self-gravitating fluid for the general view of initial distributions of its density and velosity. The authors calculate the integrals of movement and generalized Cole -Hopf substitutions, allowing to construct exact decisions. The method of construction of exact solutions is generalized for the ideal fluid in self-gravitational field on the two-dimensional and three-dimensional flows possessing cylindrical and spherical simmetry. The article analyzes the processes of formation of dense objects from a dust cloud with the set initial distribution of density of dust owing to Jeanse’s instability for flat and spherical simmetry on the basis of the constructed decisions. The researchers also calculate the characteristic time before singularity formation.

Key words: non-linear equations, self-gravitating fluids, exactly integrated models.

Введение

Исследования уравнения Бюргерса [1] одномерного течения вязкой жидкости и конструктивное доказательство его полной интегрируемости с помощью подстановки, открытой Коулом и Хопфом [2, 3], можно считать важным результатом в гидродинамике и математической физике середины

1 Настоящая работа выполнена в рамках федеральных целевых программ «Исследования и разработки по приоритетным направлениям развития научно технологического комплекса России на 2007-2012 годы» и «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы», а также работ в рамках государственного задания Минобрнауки России № 2.1894.2011 и при частичной поддержке РФФИ (проекты 11-01-00747-а, 12-01-00660-а).

ХХ в. Этот результат позволил разобраться с тем, как эволюционируют волны в вязкой среде с учетом нелинейности и понять, во что превращаются разрывные течения идеальной жидкости при учете вязкости [4]. Кроме этого, открытие подстановки Коула - Хопфа подтолкнуло (возможно не сразу [4]) к поискам методов точного интегрирования нелинейных уравнений в частных производных, что, в конце концов, вылилось в создание метода обратной задачи (МОЗ) в различных его формулировках (см. [5-7] и библиографию там) и целого ряда других методов в этой области.

К настоящему времени известно сравнительно небольшое количество точных решений уравнений Эйлера и Навье - Стокса [8, 9] для течений сжимаемой жидкости в двумерном и трехмерном случае, которые могут служить базой для решения конкретных прикладных задач. Поэтому при решении большинства прикладных задач гидродинамики приоритет отдается численным моделям динамики жидкости. Однако численные методы часто не дают возможности понять внутреннюю природу получаемых решений, поскольку процесс появления сложного решения оказывается скрыт от исследователя ненаблюдаемым потоком огромного числа простых математических операций. В этом потоке теряется физический смысл того, что происходит реально в системе, что, по всей видимости, не приближает нас к цельному пониманию гидродинамических процессов, особенно в случае течений сжимаемой жидкости. В связи с этим поиски методов построения точных решений уравнений Эйлера и На-вье - Стокса, вне зависимости от развития методов численного их анализа, будут оставаться одним из важных направлений исследования в гидромеханике и математической физике в целом (см., например, [10, 11] и библиографию там).

Часто точные решения гидродинамических уравнений строятся на основе параметризации некоторого класса течений, исходя из каких-либо сим-метрийных условий или аналитических свойств функций, входящих в параметризацию. Такой подход особенно эффективен в случае двумерных течений несжимаемой жидкости, для которых можно ввести функцию тока [8, 9] и использовать в явном виде дифференциальные законы сохранения. В трехмерном случае методы параметризации иные, но по сути сводятся к аналогичным методам (например, [10, 11]).

Для задач динамики идеальной сжимаемой среды одним из основных методов получения точной информации о решениях, в том числе и методов построения самих точных решений, является метод годографа и некоторые другие аналогичные методы (см. [12]). Однако эти методы позволяют получать решения в специфической форме зависимости координат и времени от параметров потока жидкости, что затрудняет их интерпретацию и построение решений начальных и краевых задач.

В качестве альтернативного подхода в работах [13-15] был предложен метод обобщенных подстановок Коула - Хопфа. В отличие от метода преобразования годографа метод, основанный на обобщенных подстановках Коула -Хопфа, позволяет строить решения либо в явном виде, либо в форме интеграла движения, задающего решение в неявной форме, но приспособленной для задания начальных условий для параметров среды.

Обобщение этого подхода было использовано в работах [13-15] для построения точных решений уравнений Эйлера и Навье - Стокса сжимаемой

жидкости для одномерных течений. А в работе [16] в частном случае квази-потенциальных течений аналогичный результат был получен для двумерных течений идеальной жидкости. В работе [17] был предложен еще один вариант построения аналитических решений уравнений Навье - Стокса и Эйлера на основе специальной подстановки Коула - Хопфа. Эти результаты для одномерных течений интерпретируются с несколько иных позиций, что позволяет получить еще более обобщенный результат, чем в [14, 15]. Теоретические основы приложения подхода, развитого в [17], к задачам одномерной динамики были частично исследованы в работе [18]. В настоящей работе эти результаты обобщаются и применяется к задаче построения аналитических решений течений самогравитирующей жидкости как для идеальной, так и для вязкой жидкости. На их основе анализируются решения, описывающие формирование структур в пылевых облаках вследствие неустойчивости Релея - Джинса для различных начальных конфигураций. Построение таких решений имеет важное значение для ряда задач астрофизики и космологии (см. [19] и библиографию там).

1. Обобщенные подстановки для уравнений сжимаемой жидкости

Из работы [17] следует, что параметризация одномерных течений может быть построена следующим образом. Рассмотрим представление для скорости и (х, t) потока в следующем виде:

А

и (х, 0 = , (1)

0 х

где 0 = 0(х, t) - некоторая вспомогательная функция. Эквивалентная запись этого соотношения имеет вид уравнения переноса изолиний функции 0(х, t):

01 + и (х, t )0 х = 0. (2)

Из (2) следует несколько простых тождеств, выполняющихся для

любой дифференцируемой функции 0(х, t). Первое тождество имеет

следующий вид:

Г д 1

^ (0) = ^ '(0)(0г + и0 х) = 0. (3)

д 1

— + и( х, t )— дt дх

Еще одно тождество можно получить, дифференцируя его по х:

^0 х +1Г (и( ^t )0 х ) = 0. дt дх

Это соотношение принимает форму уравнения неразрывности, если в качестве плотности жидкости р( х, t) рассматривать величину

р( х, t ) = 0Х. (4)

В этом случае

'Ч

Pt (ир) = 0. (5)

дх

Представим (5) в следующем виде:

д 1п р д 1п р

------ъ и------- их, (6)

дt дх х

и продифференцируем его по х. Полученное соотношение можно записать так:

д д ^ д 1п(р) 1 д

— + и (х, О— I—-----=-----—(рих ).

дt дх) дх р дх

Это соотношение также является тождеством, выполняющимся для произвольной функции 0(х, t).

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

и( х, 0 = Е (0) -у^^ПМ (7)

дх

где Е(0) - произвольная функция 0. Это дополнительное требование эквивалентно уравнению для функции 0 следующего вида:

0 + Е (0)0 х = У0 хх, (8)

которое получается из (7) подстановкой в него определения (1). Если 0 удовлетворяет уравнению (8), то функция и(х,^ удовлетворяет уравнению Бюргерса для сжимаемой жидкости:

1 д

иг + иих = V-—(рих), (9)

р дх

причем плотность р = 0 х удовлетворяет уравнению неразрывности (5).

Следовательно, функции и(х,t) и р(х,t) при условии, что функция 0 удовлетворяет уравнению (8), представляют точные решения задачи об одномерном течении вязкой сжимаемой жидкости при нулевых внешних силах и градиенте давления при условии сохранения ее массы.

Решения системы уравнений (5) и (9) строятся в виде подстановок (7) (или (1)) для скорости и (4) для плотности, в которых функция 0( х, t) - любое точное решение уравнения (8) при заданной функции Е(0). Наличие в схеме построения решений произвольной функции ставит вопрос о том, каким образом данная функция связана с начальными и граничными условиями. Анализ самих подстановок для скорости и плотности показывает, что вид функции Е(0) и сами полученные решения могут быть связаны с начальными условиями произвольного вида в распределении плотности и скорости жидкости в момент времени t = 0. Из (7) и (4) имеем:

р(х,0) = д0д'х), (10)

дх

и (х, 0) = и0( х) = Е (00( х)) -Уд 1п 00( х), (11)

дх

где 00 (х) = 0(х, 0) - начальное условие для функции 0(х, г) .

Явный вид функции Е(0), как видно, определяется начальным распределением скорости течения и плотности.

2. Одномерные течения самогравитирующей жидкости

Рассмотрим теперь одномерное течение вязкой жидкости в собственном поле тяготения без градиента давления. Уравнения таких течений должны удовлетворять трем уравнениям:

щ+иих = У (рих) -ф х; рг + (ри) х = 0; (12)

р дх

Ф хх =4^^р. (13)

Здесь ф - потенциал гравитационного поля; G - гравитационная

постоянная. Воспользуемся представлением для скорости (1) и плотности (4). В результате уравнение (13) примет следующий вид:

Ф хх =4^0 х.

Интегрируя это уравнение по х, находим

Фх = 4я£0 + go(t). (14)

Подставляя это выражение в (12) и умножая его на р, получаем

э

р(и + иих ) = V—(рих) -рф х. (15)

дх

Для того чтобы получить это уравнение из тождеств (3) и (6), потребуем выполнения соотношения

и = Е(0) -Уд 1п(р) - г4пО0 -у0(Г), (16)

дх

где

Мг ) = | go(г )^г.

Используя, тождество

^ д- + и д- ^ (г 4пО0 + У0 (г)) = 4я^0 + с(Г),

получаем, что функция и(х, г) (16) удовлетворяет уравнению (15) при условии, что функция 0 удовлетворяет уравнению

0Г + Е (0)0 х = У0 хх + 4пвг 00 х + У0 (Г )0 х (17)

с произвольной функцией Е(0).

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

Полученное представление решений исходной системы уравнений является полным в том смысле, что с его помощью есть возможность в принципе строить решения произвольной начальной задачи, которые определяются начальными условиями в виде двух произвольных функций р0 (х) = р(х, 0) и и0(х) = и (х,0) - начальными распределениями плотности и скорости в момент времени г = 0. Процедура вычисления функции Е(0) остается аналогичной той, которая была описана в предыдущем разделе и сводится к системе алгебраческих уравнений:

00(х) = 0(х,0) = |р(х,0)ёх, Е(00(х)) = и(х,0) + уд 1п(рхх,0)) + ^(0). (18)

В силу неотрицательности плотности вывод об однозначности вычисления Е (0) здесь также остается в силе.

3. Сферические и цилиндрические волны в самогравитирующей жидкости

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

сжимаемой жидкости в собственном поле тяготения можно записать

в следующем виде:

щ + ииг = -Фг; рг + -1 (гпри) = 0; (19)

гп дг

-1 ^ (гпФг ) = 4^р. (20)

гп дг

Здесь п = 1 для цилиндрической симметрии и п = 2 - для сферической, а г - соответствующая радиальная координата. Введем вспомогательную функцию р = гпр . В этом случае уравнения (19)-(20) можно переписать таким образом:

э

и{ + ииг = -фг; рг +—(ри) = 0; (21)

дг

'\

—(гпфг ) = 4^р. (22)

дг

Используя исходную параметризацию скорости и(г,г) = -0г / 0г (с заменой х на г) и вспомогательной функции р в форме р = 0г, уравнение

Пуассона (22) приведем к следующему виду:

Фг =~п 4^0 + go(Г), (23)

г

здесь go (Г) - произвольная функция времени, как и в одномерном случае, связанная с ускорением свободного падения на границе области, например

при г . В этом случае уравнение неразрывности выполняется тождественно, а уравнение для радиальной скорости примет такой вид:

4кG

щ + ииг =------— 0- go(t). (24)

г

По аналогии с одномерной моделью на полупрямой введем следующее условие для скорости радиального потока и (г, Г):

и = / (г )0 (0) + ^(Г), (25)

где ^(г), /(г) и Q(0) - некоторые пока произвольные функции. Используя основные тождества, получаем

иг + ииг = /'(г ^(0)и = /'(г) / (г )Q2 (0) +1>0.

Сравнивая это соотношение с (23), получаем:

Q 2(0) = 4^0, / (г) / '(г ) = -г -п, % = go(Г),

или

Q(0) = ^/4nGy|Q, ^(г) = |go(г)Л /(г)

= +

л/2^/ г0 1 + г 1, п = 2, л/2^1п г0 - 1п г, п = 1.

Здесь г0 - постоянная интегрирования.

Отсюда следует, что сферическая волна в сжимаемой самогравитиру-ющей среде без градиента давления описывается соотношениями

р = -20г, и = -0г / 0г = +>/8л^л/г0-1 + г-1 л/б + у0(Г) (26)

г2

при условии, что функция 0 удовлетворяет уравнению

0Г + МО ^л/^^/г0-1 + г-17б )0г = 0. (27)

Соответственно цилиндрическая волна в сжимаемой самогравитирую-щей системе описывается соотношениями

р = 10г, и = -0г / 0г = +л/8nG^JQJ/0 - 1п г + ^(г) (28)

г

с функцией 0 , удовлетворяющей уравнению

0г + (^(г) ++л/8л^^1п г0 - 1п )0г =0. (29)

Построенные модели для цилиндрического и сферического случая отличаются друг от друга важной особенностью. Если в цилиндрическом случае построенная модель, как и в одномерном плоском течении, дает описание течений на ограниченном интервале значений радиальной координаты [0, ^], то аналогичная модель для сферического случая при г0 > 0 дает описание на всем пространстве. Конечность интервала

вещественности в цилиндрической модели означает наличие цилиндрической стенки, на которой скорость потока обращается в ноль (при ^=0). Это ограничивает возможность применять модель с п = 1 для астрофизических задач. Однако модель с п = 2 свободна от этого недостатка. Как и в плоском случае, динамика среды, движущейся к началу координат и от него, описывается разными моделями, при этом начальное распределение скорости в этих моделях непроизвольное и связано с начальным распределением плотности среды.

Распределение плотности р0 (г) = р(г, 0) задает начальное распределение функции 00(г) = 0(г,0). Эта функция должна быть всюду неотрицательной и монотонно неубывающей с ростом г . Кроме этого, при выборе начального распределения р следует учитывать, что в отсутствие сингулярности плотности в начале координат необходимо требование Нт 0г = гп .

г ^0

4. Точные решения одномерных уравнений течений идеальной жидкости

В случае отсутствия вязкости система уравнений (12)-(13) допускает интеграл движения, позволяющий строить решения в неявной форме для произвольных начальных распределений плотности и скорости. Этот интеграл движения имеет следующий общий вид:

Н (0, х - X (0, г)) = 0, (30)

где Н = Н (0, ^) - произвольная дифференцируемая функция двух аргументов (0 и ^ = х - X(0,Г)). Функция X(0,Г) для рассматриваемого случая имеет следующий вид:

X (0, Г) = Е (0)г - 2^г 20 -1 у0 (Г )Л. (31)

Прямой проверкой убеждаемся, что функция 0 , удовлетворяющая алгебраическому уравнению (30), также удовлетворяет уравнению

0Г 0х =0, (32)

которое эквивалентно (17) при V = 0

0г + (Е (0) - 4^Г 0- у0(Г ))0 х = 0 (33)

для X(0, Г), определенного в (31).

Уравнение (33) (как и общее уравнение (32)) является уравнением простой волны и имеет многозначные решения, указывающие на обрушение волн и образование ударных волн в среде, динамику которой оно описывает. Как и в случае вязкой среды, наличие двух произвольных функций Н (0, ^) и Е (0) позволяет формально строить решения задачи о динамике сжимаемой среды с нулевым давлением для произвольных начальных распределений ее плотности и скорости. Процедура вычисления Е (0) остается прежней, а вид функции Н (0, ^) вычисляется из начального распределения плотности.

5. Джинсовская неустойчивость в одномерном случае

Полученные решения позволяют на точных решениях рассмотреть динамику локальных возмущений плотности под действием собственного поля тяготения в одномерном случае. Эта задача о джинсовской неустойчивости самогравитирующей среды важна для описания динамики крупномасштабной структуры Вселенной, как это уже объяснялось во введении. При анализе формирования крупномасштабной структуры предполагают, что точками среды являются галактики. При этом давление в такой среде равно нулю, что соответствует уравнению состояния пыли р = 0. Кроме этого, вязким трением в такой среде можно пренебречь. В результате динамика таких систем описывается в точности системой уравнений (12)-(13) с V = 0. Решения, полученные в предыдущем параграфе, позволяют строить точные решения этой задачи при заданных начальных условиях распределения плотности и скорости потока. Начальные условия при V = 0 задаются соотношениями

00 (х)

р( х,0) = ^А и( х,0) = Е (00( х)) - с(0), (34)

ох

где

00( х) = 0( х,0).

Для получения физически понятных решений необходимо предварительно выяснить то, какие распределения плотности, скорости и гравитационного потенциала допустимы с физической точки зрения. Основным требованием, ограничивающим выбор функции 00(х), является неотрицательность плотности жидкости во всех точках пространства: р> 0. В этом случае имеем

0(х, Г) = ©0 (Г) + | р(х, Г)ёх > 0.

—го

Из уравнения Пуассона следует соотношение (14), которое связывает функцию 0 с ускорением свободного падения в точке с координатой х . Отсюда следует

Фх > g0(t). (35)

По своему смыслу -фх является ускорением свободного падения в точке пространства с координатой х . Из физических соображений естественным выбором постоянной go (г) является такой, при котором ускорение свободного падения было бы равно нулю при х — +°°. Физически это связано с требованием р — 0 в этих же пределах. Однако, если полагать

Пт фх = 0, (36)

х——-го

то из соотношения (35) следует, что при х —— +го ускорение свободного падения стремится к некоторой отрицательной постоянной, равной

gгo = -4^м,

где

го

М = | р(х, Г)йх > 0

—го

есть полная масса материи в пространстве.

Поскольку р> 0, то отсюда следует, что в одномерном случае локализованное возмущение плотности создает на бесконечности постоянное поле тяготения. Выбор граничных значений напряженности гравитационного поля на бесконечности определяет в реальности значение скорости, которую приобретает сформировавшаяся структура при г — го. Если gгo = g-гo, то центр масс всей системы приобретает некоторое ускорение. Это ускорение такое, что в системе, движущейся вместе с центром масс, gГО = g-го. Здесь gГО, g -го - величины ускорения на бесконечности в системе центра масс.

В качестве примера демонстрации некоторых важных эффектов в одномерном случае рассмотрим начальное распределение плотности, которое задается простой функцией 00 (х) следующего вида:

00( х) = ММ-(гЬ(кх) + 1),

где к >0 - произвольное положительное вещественное число, величина которого обратно пропорциональна характерному пространственному

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

выборе g-гo =0, gгo = -4кОМ . При этом начальное распределение плотности и скорости будет иметь следующий вид:

р( х, 0) = —ММк— >0, и( х,0) = 0. (37)

2еЬ (кх)

В этом случае решение задачи о динамике такого возмущения дается в неявном виде как решение алгебраического трансцендентного уравнения:

Ш ^к (х + ц09ґ2)) + 1

(38)

Здесь ^0 =2кО . Решения для плотности и скорости потока можно получить из интеграла движения. Дифференцируя (38) по х и г, находим:

р = 0х =----/-------МЬ----------^, и = -2Ц0Г0. (39)

2оЬ1к(х + ^00Г ))~2МкГ ^0

На рис. 1 представлены распределения плотности в последовательные моменты времени Г = 0; 2; 4; 6; 8; 0,95 .

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

р(х,0) =----МТТ- >0, и(х,0) = 0 (40)

л(1 + к х )

представлены на рис. 2. Для этого случая имеем:

_ М ^ Мк

0 = —(к(х + ^00г )), р = 0х =-------^---------21-2—, и = -2^0Г0-

п к + кк2(х + ц0 0г2)2 -Мкг2ц0

а) б)

Рис. 1. Эволюция распределения скорости (а) и плотности (б) для начальных распределений (37)

” -2 -1 0 : X 2

а) б)

Рис. 2. Эволюция распределения скорости (а) и плотности (б) для начальных распределений (40)

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

I кМ^0 V пОМ

Ь

а для (40):

п = Ь кМц0 V ОМ'

6. Джинсовкая неустойчивость для распределений со сферической симметрией

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

0 = 7М-б(х(г, г)), (41)

4п00

где М - полная масса вещества во всем пространстве,

го

а=\&(г )йт=е(го) - е(0)

0

и

*(г, г) = г 3/2 + V 8тсО7\/о.

В этом случае функция 6(*) задает форму начального распределения функции 0 и в конечном итоге плотности, которая вычисляется в этом случае так:

р(г, г) = 3М3/2______ ))________.

4я<20г3/2 8п^л/б + М44пОд'(%(г, г))г Множитель М /(47100) в (41) выбирался из условия

го

М = 4п | рг 2ёг.

0

Выбор самой функции 6(*) должен определяться требованиями 0< 60 < го и условием монотонного неубывания этой функции с ростом г . Кроме этого, исходя из того, что по крайней мере в начальный момент времени в начале координат отсутствует сингулярность в распределении

2

плотности вещества (р = 0г / г ), зависимость 6(*) должна быть такой, что:

6(*(г ,0)) = К (г 3),

где К(^) - всюду непрерывная монотонно неубывающая функция, равная 0 при £, = 0 :

К(0) = 0.

Это обеспечивает в нуле конечное значение плотности.

Исходя из общих условий выбора представления (41) и по аналогии с одномерным случаем рассмотрим интегралы движения двух типов. Первый интеграл

0(г, г) = М (к (г3/2 +л/8пОг70)2) (42)

4п' !

-2

соответствует следующим начальным условиям для р = 0г и и при г = 0 :

00(г,0) = МП(кг3); (43)

р(г ,0)= 3Мк 3 >0, и (г ,0) = -л/вПОг-1/200 (г). (44)

4п2 (кг 3)

Здесь к >0 - как и раньше, произвольное положительное вещественное число, величина которого обратно пропорциональна характерному пространственному размеру Ь области, занятой материей: Ь = к 1, а М - полная масса системы, 60 =1.

Плотность и скорость потока (рис. 3, 4) могут быть вычислены через значения функции 0(г, г) по формулам

( г) 3Мкл/ё(1 + Ц0^7ёг-3/2) ( г) 0

р(г, Г)=-----2-----------0--------------, и (г, г) = -Ц0>/ 0,

4п7э (к*2 (г, Г)) -Мк^0Г*(г, Г)

где *(г, Г) = г3/2 + |10Гл/б , |10 = ^8пО .

Второй интеграл

М (К(г3/2 2п

0(г, Г )=—г (к (г 3/2 +78лОгл/0 )2) (45)

2п2 1 1

соответствует начальным условиям:

00(г,0) = ^(кг3); (46)

2п2

р(г,0) = М---------12->0, и(г,0) = 0. (47)

4п 1 + к2 г 6

Для этого случая 60 = п / 2 .

Анализ асимптотического поведения функций 0 и р вблизи нуля дает следующие соотношения:

0 = г 3(00 (Г) + 04 (Г )г6 + 0(г12)), р = 300 (Г) + 904 (Г )г6 + 0(г12).

Уравнение для 00 (г), следующее из (42), имеет вид

00 = (1 +л/°0г^0)2,

где ^0 =4пО . Остальные коэффициенты разложения определяются 00 (г).

Имеются два решения относительно 0+ , соответствующие двум общим решениям уравнения (42) для функции 0(г, г):

0± (г) =_________Мк_______

0 4п(1 + ^Г^Мк / 4п)

а)

б)

Рис. 3. Эволюция распределения скорости (а) и плотности (б) для начальных распределений (40) до образования сингулярности: 1 - г = 0; 2 - г = 0,5;

3 - г = 1,0; 4 - г =1,5; 5 - г = 2,0; 6 - г = 2,5; 7 - г = 3,0; 8 - г = 3,5; 9 - г = 4,0

Одно из этих решений при г >0 за конечное время

4л 1

Мк ^о

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

а) б)

Рис. 4. Распределения скорости (а) и плотности (б) для начальных распределений (40) до и после образования сингулярности: 1 - t = 0 ; 2 - t = 4,0; 3 - t = 4,8 ; 4 - t = 5,0

Заключение

Как продемонстрировано в данной работе, развитый метод построения моделей динамики сжимаемой жидкости, допускающих преобразование к интегрируемым уравнениям с помощью обобщенных подстановок Коула -Хопфа, оказывается эффективным при решении ряда задач динамики само-гравитирующей среды. В рамках этого подхода можно на основе аналитических решений рассмотреть ряд важных вопросов формирования структур в пылевых облаках, что важно для целого ряда задач астрофизики [19]. В данной работе не проводилось сравнение между формой решений, которые можно получить с помощью метода годографа [12, 19] и данным методом. Однако, как показывают результаты работы [19], с помощью метода годографа можно эффективно исследовать асимптотическое поведение решений. Данный же метод позволяет получать более тонкую информацию о точном решении для разнообразных начальных условий. Следует также отметить, что данный метод достаточно эффективно работает как в случае идеальной жидкости, так и в случае вязкой среды. Решения для вязкой среды в рамках метода годографа получить не удается. Поэтому данный метод имеет гораздо более широкую область применения.

Список литературы

1. Burgers, J. M. The nonlinear diffusion equation / J. M. Burgers. - Dordrecht, Holland: D. Reidel Publisher Company, 1974. - P. 185.

2. Hopf, E. The partial differential equation ut + uux = iwuxx / E. Hopf // Comm. Pure and Appl. Math. - 1950. - V. 3. - P. 201-230.

3. Cole, J. D. On a quasilinear parabolic equation occurring in aerodynamics / J. D. Cole // Quart. Appl. Math. - 1951. - V. 9. - № 3. - P. 225-236.

4. Уизем, Дж. Линейные и нелинейные волны / Дж. Уизем. - М. : Мир, 1978. -622 c.

5. Солитоны / под ред. Р. Буллафа, Ф. Кодри. - М. : Мир, 1983. - 408 с.

6. Захаров, В. Е. Теория солитонов: метод обратной задачи / В. Е. Захаров, С. В. Манаков, С. П. Новиков, Л. П. Питаевский. - М. : Наука, 1980. - 320 с.

7. Додд, Р. Солитоны и нелинейные волновые уравнения / Р. Додд, Дж. Эйлбек, Дж. Гиббон, Х. Моррио. - М. : Мир, 1988. - 694 с.

8. Ландау, Л. Д. Теоретическая физика. Т. VI. Гидродинамика / Л. Д. Ландау, Е. М. Лифшиц. - М. : Наука, 2003. - 732 с.

9. Кочин, Н. Е. Теоретическая гидромеханика / Н. Е. Кочин, И. А. Кибель, Н. В. Розе. - М. : Наука, 1963. - 584 с.

10. Полянин, А. Д. Системы уравнений гидродинамического типа: точные решения, преобразования, нелинейная устойчивость / А. Д. Полянин, С. Н. Аристов // Доклады академии наук. - 2009. - Т. 428, № 2. - С. 180-185.

11. Броман, Г. И. Затопленная струя Ландау: точные решения, их смысл и приложения / Г. И. Броман, О. В. Руденко // Успехи физических наук. - 2010. - Т. 180, № 1. - С. 97-104.

12. Полянин, А. Д. Методы решения нелинейных уравнений математической физики и механики / А. Д. Полянин, В. Ф. Зайцев, А. И. Журов. - М. : Физматлит, 2005. - 256 с.

13. Журавлев, В. М. О нелинейных уравнениях, связанных с уравнениями теплопроводности и Д'Аламбера с помощью подстановок типа Коула - Хопфа /

B. М. Журавлев, А. В. Никитин // Нелинейный мир. - 2007. - Т. 5, № 9. -

C. 603-611.

14. Журавлев, В. М. Нелинейные уравнения, линеаризуемые с помощью обобщенных подстановок Коула - Хопфа, и точно интегрируемые модели одномерных течений сжимаемой жидкости / В. М. Журавлев, Д. А. Зиновьев // Письма в ЖЭТФ. - 2008. - Т. 87, № 5. - С. 314-318.

15. Журавлев, В. М. Метод обобщенных подстановок Коула - Хопфа и новые примеры линеаризуемых нелинейных эволюционных уравнений / В. М. Журавлев // Теоретическая и математическая физика. - 2009. - Т. 158, № 1. - С. 58-71.

16. Журавлев, В. М. Метод обобщенных подстановок Коула - Хопфа в размерности 1+2 и интегрируемые модели двумерных течений сжимаемой жидкости /

B. М. Журавлев, Д. А. Зиновьев // Письма в ЖЭТФ. - 2008. - Т. 88, № 3. -

C. 194-197.

17. Журавлев, В. М. Инновационные технологии / В. М. Журавлев ; под ред. проф. С. В. Булярского. - Ульяновск : Изд-во УлГУ, 2010. - С. 77-93.

18. Zhuravlev, V. M. Nonlinear Waves in Self-Gravitating Compressible Fluid and Generalized Cole-Hopf Substitutions / V. M. Zhuravlev, D. A.Zinov’ev // Physics of Wave Phenomena. - 2011. - V. 19, № 4. - P. 313-317.

19. Гуревич, А. В. Крупномасштабная структура Вселенной. Аналитическая теория / А. В. Гуревич, К. П. Зыбин // Успехи физических наук. - 1995. - Т. 165, № 7. - С. 723-758.

20. Рождественский, Б. Л. Системы квазилинейных уравнений и их приложения к газовой динамике / Б. Л. Рождественский, Н. Н. Яненко. - М. : Наука, 1978. - 689 с.

21. Куликовский, А. Г. Математические методы изучения разрывных решений нелинейных гиперболических систем уравнений / А. Г. Куликовский, Е. И. Свешникова, А. П. Чугайнова. - М. : МИАН, 2011. - 120 с.

Журавлев Виктор Михайлович

доктор физико-математических наук, профессор, кафедра теоретической физики, Ульяновский государственный университет

E-mail: zhvictorm@mail.ru

Зиновьев Дмитрий Александрович ассистент, кафедра теоретической физики, Ульяновский государственный университет

E-mail: zinovjev_d@mail.ru

УДК 532.5, 52-43, 514.86 Журавлев, В. М.

Интегрируемые модели динамики сжимаемой среды в собственном поле тяготения. Метод подстановок Коула - Хопфа / В. М. Журавлев, Д. А. Зиновьев // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2012. - № 4 (24). - С. 174-190.

Zhuravlyov Viktor Mikhaylovich Doctor of physical and mathematical sciences, professor, sub-department of theoretical physics, Ulyanovsk State University

Zinovyev Dmitry Alexandrovich Assistant, sub-department of theoretical physics, Ulyanovsk State University