Нахождение приближенных аналитических решений линеаризованного по скорости уравнения Навье-Стокса в сфероидальной системе координат

Размещено на http://www.allbest.ru/

Нахождение приближенных аналитических решений линеаризованного по скорости уравнения Навье-Стокса в сфероидальной системе координат

дифференциальный уравнение задача ньютоновский

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

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

Частицы, входящие в состав реальных систем (газ, жидкость), могут иметь форму поверхности, отличной от сферической, например, сфероидальную (вытянутый или сплюснутый эллипсоид вращения). В данной работе получено решение линеаризованного по скорости уравнения Навье-Стокса в сфероидальной системе координат. При этом предполагается степенной вид зависимости коэффициентов вязкости , теплопроводности , и плотности газообразной среды от температуры (, , , , , ). Здесь и далее индексы «g» и «« будем относить соответственно к газообразной среде и частице; индексом «« обозначены параметры газообразной среды на бесконечности в невозмущенном потоке, а индексом «« - значения физических величин, взятых при средней температуре поверхности частицы, равной .

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

Предполагается также, что коэффициент теплопроводности частицы по величине много больше коэффициента теплопроводности газа (). Это допущение приводит к тому, что в коэффициенте динамической вязкости можно пренебречь зависимостью от угла в системе частица - газообразная среда (слабая угловая асимметрия распределения температуры). В этом случае можно считать, что . Здесь , ; и определяются из решения тепловой задачи. Это допущение позволяет рассматривать гидродинамическую часть отдельно от тепловой части, а связь между ними осуществляется через граничные условия.

В рамках сформулированных выше допущений система газодинамических уравнений, описывающая распределения полей скорости , давления и температур вне и внутри частицы, имеет вид [2, 3]

, , (1)

, , (2)

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

Система (1), (2) решается со следующими граничными условиями в сфероидальной системе координат :

, , , ,

, ,

, , .

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

Исследования показали, что при поиске выражений для компонент массовой скорости в виде

, , (3)

где и - некоторые произвольные функции, зависящие от радиальной координаты (связь между функциями и определяется из уравнения непрерывности (2) - см. формулу (6)), задачу можно свести к обыкновенному неоднородному дифференциальному уравнению третьего порядка для функции , решение которого можно получить в виде обобщенных степенных рядов [4-8].

В выражении (3) введены следующие обозначения: - в случае вытянутого сфероида (, формула (4)) и - в случае сплюснутого сфероида (, формула (5)); и - полуоси сфероида; , , , - коэффициенты Ламэ. Декартовы координаты связаны со сфероидальными координатами соотношениями (4) - (5)

, , , (4)

, , , (5)

, где , . (6)

Функция , входящая в (6), находится из решения тепловой задачи и имеет вид:

, (7)

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

С учетом выражения (7) формула для коэффициента динамической вязкости принимает вид

. (8)

Подставляя (8), (6) в (1), получаем следующее уравнение для функции в системе координат сплюснутого сфероида

(9)

с краевыми условиями

, , , , (10)

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

В уравнении (9) перейдем к новой переменной , в результате имеем

. (11)

Здесь , , , , , , ; , .

Точка =0 для однородного уравнения, соответствующего неоднородному уравнению (11), является регулярной особой точкой [4, 5]. Поэтому решение однородного дифференциального уравнения будем искать с помощью обобщенных степенных рядов. Для этого разложим в степенные ряды в окрестности нуля функции, входящие в уравнение (11). При этом основную сложность составляет поиск разложения в ряд функции L. Построим степенной ряд таким образом, чтобы он начинался с 1 и аппроксимировал нашу функцию L с достаточно высокой точностью. Ряд вида удовлетворяет всем выше перечисленным требованиям, причем погрешность приближения не превышает 5%. Ниже приведены графики, показывающие с какой точностью полученный ряд аппроксимирует исходную функцию L при различных значениях температуры (количество и значения коэффициентов для каждой температуры различны). Сплошной линией изображена исходная функция L, а точками - аппроксимирующая ее функция.

После разложения функций уравнения (11) в ряды, имеем следующее уравнение:

, . (12)

Рекуррентные соотношения для коэффициентов , получаемых при перемножении рядов для функции и L в степени m (m=1,2,3), имеют громоздкий вид, поэтому в данной статье они приведены не будут.

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

, . (13)

Вычисляя производные, получим:

, ,

.

Подставляя (13) в однородное уравнение, соответствующее неоднородному уравнению (12), и приравнивая коэффициенты при , получаем определяющее уравнение

,

корни которого равны соответственно , , .

Заметим, что разность корней равна целому числу. Следовательно, согласно общей теории решения дифференциальных уравнений в виде обобщенных степенных рядов (методом Фробениуса) во всех остальных решениях, кроме первого, соответствующего , появляется дополнительное слагаемое, содержащее множитель , умноженный на первое решение [4, 5].

Таким образом, учитывая значения корней определяющего уравнения, система линейно независимых решений однородного уравнения, соответствующего неоднородному уравнению (9), имеет вид:

, ,

. (14)

Из вида правой части неоднородного уравнения (9) следует, что его частным решением будет функция

, . (15)

Подставляя (14), (15) в уравнение (12) и используя метод неопределенных коэффициентов, получим следующую рекуррентную формулу для нахождения коэффициентов :

,

,

.

Через обозначена целая часть числа .

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

, , , ,

, , , ,

,

,

.

Заметим, что выбор постоянных , и осуществляется таким образом, чтобы функции , и стремились к соответствующим функциям для эллипсоида вращения при малых относительных перепадах температуры [2]. Четвертое решение уравнения (12), вычисляемое по формуле (14), не удовлетворяет краевому условию (10) , поэтому его явный вид мы приводить не будем. Таким образом, мы получили общее решение уравнения (12)

.

Функция формально удовлетворяет уравнению (12) по построению. Можно показать, что ряды, определяющие функции , i=1,2,3,4 равномерно сходятся при , причем радиус сходимости рядов равен единице. Постоянные интегрирования , , и однозначно определяются из граничных условий (10). Очевидно, , . Для определения постоянных , имеем следующую систему уравнений:

Эта система имеет единственное решение, поскольку ее главный определитель отличен от нуля в силу линейной независимости решений , и . Таким образом, имеем:

, . (16)

В результате проведенного исследования доказана следующая теорема.

Теорема. Функция с коэффициентами, определяемыми формулой (16), является единственным решением уравнения (9), удовлетворяющим краевым условиям (10).

Заметим, что полученное решение справедливо при . Если интересующая нас область не принадлежит указанному интервалу, то можно найти аналитическое продолжение полученного решения.

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

Заключение

Известный в научной литературе метод решения обыкновенных дифференциальных уравнений с помощью обобщенных степенных рядов (метод Фробениуса) был применен для нахождения аналитического решения линеаризованного по скорости уравнения Навье-Стокса с учетом зависимости коэффициента динамической вязкости от температуры. Специальным выбором выражений для компонент массовой скорости при условии, что коэффициент теплопроводности частицы по величине много больше коэффициента теплопроводности газообразной среды (слабая угловая асимметрия распределения температуры в системе частица-газ), система уравнений Навье-Стокса сводится к неоднородному дифференциальному уравнению третьего порядка с изолированной особой точкой. Решение этого уравнения находилось в виде обобщенных степенных рядов. Доказана теорема существования и единственности полученного решения.

Литература

1. Ладыженская О.А. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Наука, 1970. 288 с. 1. Ladyzhenskaya O.A. Mathematical questions of dynamics of a viscous incompressible liquid. М.: Nauka, 1970. 288 pp.

2. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. М.: Мир, 1976. 630 с. 2. Happel J., Brenner H. Low Reynolds number hydrodynamics. М.: Mir, 1976. 630 pp.

3. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. М.: Наука, 1986. 736 с. 3. Landau L.D., Lifshits E.M. Hydrodynamics. М.: Nauka, 1986. 736 pp.

4. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука. 1961. 703 с. 4. Kamke E. Reference book on the ordinary differential equations. М.: Nauka. 1961. 703 pp.

5. Коддингтон Э.А., Левинсон Н. Теория обыкновенных дифференциальных уравнений. М.: Иностр. лит., 1958. 474 с. 5. Koddington E.A., Levinson N. The theory of the ordinary differential equations. М.: Inostr. Lit., 1958. 474 pp.

6. Привалов И.И. Введение в теорию функций комплексного переменного. М.: Наука, 1967. 444 с. 6. Privalov I.I. Introduction in the theory of functions of the complex variable. М.: Nauka, 1967. 444 pp.

7. Двайт Г.Б. Таблицы интегралов и другие математические формулы. М.: Наука, 1978. 222 с. 7. Dvayt G.B. Tables of integrals and other mathematical formulas. М.: Nauka, 1978. 222 pp.

8. Бретшнайдер Ст. Свойства газов и жидкостей. Инженерные методы расчета. М.: Химия, 1966. 535 с. 8. Bretshnaider St. Properties of gases and liquids. Engineering methods of calculation. М.: Chemistry, 1966. 535 pp.

Размещено на Allbest.ru