Scientific journal
Название журнала на английском
ISSN 2500-0802
ПИ №ФС77-61154

INVESTIGATION OF THE SOLUTION OF THE NAVIER – STOKES

Jakubowski E.G. 1
1 NMSU
The solution of nonlinear partial differential equations can determine the value of the dimensionless unknown functions with a high value (for example, a large Reynolds number). At the same time, they are reduced to a countable number of ordinary nonlinear differential equations. Turbulent solutions that correspond to large values ​​of the unknown function, are complex. The transition from the real solutions to complex turbulent solution is realized through the right side of the infinity of ordinary differential equations, which reduces the Navier – Stokes equations. At the same time the real solution of the Navier – Stokes equation determines the function tends to infinity. This comprehensive solution of course. In this article we defined form of the solution. The decision depends on whether multiple equilibrium position or not. And in the case of multiple coordinate position of the equilibrium observed chaotic solution, the structure of which is described in Theorem 2.3. Calculate the coefficient of drag of fluid flow in a circular pipe at various pipe wall roughness.
Navier – Stokes equations
the solution of nonlinear equations
the drag coefficient of the liquid flow

1. Учет ламинарного решения уравнения Навье – Стокса

Значение коэффициента сопротивления круглого трубопровода при произвольном числе Рейнольдса и степени шероховатости известно только из эксперимента. Предлагается с помощью комплексного решения получить решение уравнения Навье – Стокса и на основе качественных соображений определить влияние шероховатости на решение уравнения Навье – Стокса. При этом с точностью 10 % удалось построить классические графики Никурадзе зависимости коэффициента сопротивления круглого трубопровода в зависимости от произвольного числа Рейнольдса и степени шероховатости.

Проблема описания движения жидкости в турбулентном режиме не решена. Это создает трудности при расчете нефтяных, газовых трубопроводов. Кроме того, не существует теоретических методов описания движения тела в турбулентной среде. Это и описание движения летательного аппарата, движения подводной лодки или надводного корабля в турбулентном режиме. Без использования моделирования движущегося тела в аэродинамических трубах, или водных бассейнах, невозможно проектирование движущихся в вязкой среде тел. Существуют приближенные формулы аппроксимации коэффициента сопротивления трубопровода на отдельных участках числа Рейнольдса см. [1, 2]. Известны классические экспериментальные графики Никурадзе по определению коэффициента сопротивления круглого трубопровода в зависимости от числа Рейнольдса и степени шероховатости. Применяется аппроксимация конвективного члена, сводящая задачу решения уравнения Навье – Стокса к линейной с эффективной турбулентной вязкостью. Но такое преобразование искажает решение уравнения Навье – Стокса и для совпадения с экспериментом коэффициент турбулентной вязкости может иметь произвольное значение, вплоть до отрицательного. Применяется метод Галеркина, сводящий решение задачи гидродинамики к системе нелинейных обыкновенных дифференциальных уравнений. Но при турбулентном режиме эта система нелинейных уравнений имеет комплексные положения равновесия, т.е. решение комплексное. И действительно в турбулентном режиме система уравнений гидродинамики в действительной плоскости не имеет решения, решение уравнения стремится к бесконечности см. [3]. При этом комплексное решение конечно. Физический смысл комплексного решения и его колебательный характер мнимой части см. [4, 5]. Т.е. решать задачу гидродинамики в турбулентном режиме надо в комплексной плоскости. При этом турбулентное решение определяется не однозначно, а имеется счетное количество ветвей решения.

1.1. Расчет трубопровода с круглым сечением в случае несжимаемой жидкости

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

jak2.wmf,

где jak3.wmf давление в начальной и конечной части трубопровода, величина L это длина трубопровода, радиальной и угловой скоростью пренебрегаем. Внешнее воздействие равно

jak4.wmf.

Согласно формуле (1.2.2) градиент давления равен

jak5.wmf.

Получим уравнение

jak6.wmf.

Подставляем значение скорости, получим уравнение

jak7.wmf.

Умножаем это уравнение на радиус и интегрируем это уравнение по радиусу, так как система координат цилиндрическая, получим

jak8.wmf.

Для получения счетного количества решений, умножаем это уравнение на функцию jak9.wmf и интегрируем по пространству. Получаем счетное количество решений, как для гладкой, так и шероховатой поверхности. При этом, как и у уравнения Шредингера получается счетное количество решений, каждое со своей собственной энергией. При переходе из одного состояния в другое излучается дискретная энергия. Минимум собственной энергии определяет выбор решения.

Беря модуль правой части, и находя средний модуль тангенса угла отклонения, получим

jak10.wmf. (1.1.1)

В случае выбора знака минус у значения среднего модуля тангенса угла отклонения jak11.wmf, получится, что шероховатости увеличивают скорость потока, так как полная производная jak12.wmf увеличится, что не правильно, скорость потока должна уменьшаться за счет шероховатостей.

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

jak14.wmf.

Т.е. конвективный член нужно брать с учетом минуса, причем правую часть (1.1.1) брать со знаком плюс. Кроме того, надо выбирать знак плюс у среднего тангенса модуля наклона шероховатости, чтобы получилось комплексное турбулентное решение. В противном случае, решение, описывающее пульсирующий турбулентный режим не получится.

Переходя от радиуса к диаметру трубопровода и разделив на величину jak15.wmf, получим

jak16.wmf. (1.1.2)

Если микро-шероховатости jak17.wmf распределены по всей поверхности трубопровода, они находятся и на макро шероховатостях и определяют критическое число Рейнольдса и коэффициент сопротивления при числе Рейнольдса, равном 2300. Микро-шероховатости имеют молекулярную природу и определяются средним размером атома, равным среднему геометрическому между размером ядра rA, и размером орбиты Бора jak18.wmf, при расстоянии между атомами jak19.wmf, равному некоторой величине, определяемой свойствами границы трубопровода, железом, титаном и углеродом. Расстояние между атомами железа jak20.wmf, между атомами титана jak21.wmf, между атомами углерода jak22.wmf см. [7]. При этом абсолютная величина тангенса наклона высоты микро-шероховатости поверхности металла в трубопроводе определяется по формуле

jak23.wmf.

Средний тангенс наклона равен

jak24.wmf

Величина критического числа Рейнольдса относительно диаметра равна jak25.wmf. Но почему критическое число Рейнольдса для сферы равно 3⋅105. Это связано с разным определением критического числа Рейнольдса. Оно равно

jak26.wmf,

где величина jak27.wmf эффективный, гидродинамический размер тела, включая среду, a истинный геометрический размер тела, причем jak28.wmf молекулярный тангенс наклона шероховатости. Причем отношение jak29.wmf может иметь значение jak30.wmf.

Величина критического числа Рейнольдса равна jak31.wmf. Макро-шероховатости jak32.wmf более редкие и определяют коэффициент сопротивления на числах Рейнольдса в 12 и более раз больше.

Получаем условие стационарности для уравнения Навье – Стокса с учетом одного члена ряда-решения в одномерном случае

jak33.wmf.

В одномерном случае при постоянстве сечения трубопровода уравнение неразрывности выполняется тождественно. Ламинарное решение этого уравнения равно

jak34.wmf.

При внешнем давлении, равном jak35.wmf начинается комплексное решение и турбулентный режим, так как число Рейнольдса в этой точке равно критическому значению. Из эксперимента и путем проделанного вычисления имеем значение критического числа Рейнольдса для круглого трубопровода jak36.wmf. Коэффициент сопротивления трубопровода с круглым сечением определяется по формуле (в формулу подставляем перепад давления, выраженный через безразмерное давление)

jak37.wmf

Средняя скорость, входящая в число Рейнольдса, равна

jak38.wmf.

Асимптотика коэффициента сопротивления трубопровода jak39.wmf с круговым сечением для ламинарного режима вычислена верно.

jak40.wmf.

Асимптотика получена при малом числе Рейнольдса, когда конвективный член мал.

В случае большого перепада давления, получаем комплексное турбулентное решение jak41.wmf. Если считать точнее, то вклад вращательной мнимой части в поступательную скорость движения потока соответствует корню из мнимой части согласно формуле (1.1.3)

jak42.wmf (1.1.3)

причем необходимо использовать величину отношения числа Рейнольдса на корень из безразмерного давления, как величину порядка единицы в турбулентном режиме. Число Рейнольдса потока пропорционально на бесконечности давления jak43.wmf. При этом степень усреднения пропорциональна высоте шероховатости, стремящейся к нулю, хотя самая гладкая поверхность соответствует среднему модулю тангенса наклона, равному обратной величине критического числа Рейнольдса. Так как из мнимой части выражения для числа Рейнольдса потока извлекается корень четвертой степени, имеем для эффективного диаметра jak44.wmf формулу

jak45.wmf.

Величины jak46.wmf, и при бесконечно малой шероховатости сокращаются, хотя такая ситуация и не реализуется. Но эти формулы справедливы для данной схемы решения. При решении в виде ряда получится вычисляемое другое значение jak47.wmf. Самая гладкая поверхность соответствует среднему модулю тангенса наклона, равного обратному значению критического числа Рейнольдса, так как самые малые модули тангенса наклона соответствуют молекулярному уровню шероховатости. При этом эффективный диаметр меньше истинного диаметра. Высота шероховатости может стремиться к нулю, одновременно с периодом шероховатости l. Т.е. величина jak48.wmf это максимальное отношение эффективного диаметра к истинному диаметру. Для внешней задачи эффективный диаметр увеличится, и коэффициент b определится по формуле

jak50.wmf.

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

jak52.wmf.

При этом коэффициент b пропорционален

jak53.wmf,

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

Вычислим изменение диаметра трубопровода при условии равенства макро шероховатостей микро-шероховатостям из эмпирической формулы по уменьшению среднего квадрат значения диаметра, получится, что диаметр уменьшился в jak56.wmf раза. При отношении jak57.wmf, получаем уменьшение диаметра трубопровода в jak58.wmf раза.

При этом диаметр изменяется только у коэффициента пульсирующей части решения, т.е. у мнимой части, откуда и возникает множитель jak59.wmf, так как мнимый член пропорционален jak60.wmf, который усредняется. Квадратный корень jak61.wmf при этом соответствует среднему значению диаметра.

При этом учтено влияние шероховатости стенок турбулентного потока на мнимую часть числа Рейнольдса потока. Чтобы получить графики с постоянной высотой шероховатости, надо ввести эффективный средний модуль тангенса угла наклона шероховатостей. Эффективный средний модуль тангенса угла шероховатостей должно зависеть от внешнего давления jak62.wmf.

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

jak63.wmf,

где k – среднеквадратичная высота шероховатости, jak64.wmf – радиус круглого сечения трубопровода.

Формула подобрана таким образом, чтобы правильно определять зависимость числа Рейнольдса относительно внешнего давления и коэффициент сопротивления трубопровода при бесконечном числе Рейнольдса и внешнего давления

jak65.wmf

при коэффициенте сопротивления, равном величине

jak66.wmf.

Эти же коэффициенты получаются при решении с помощью ряда, но вместо величины jak67.wmf в случае ряда используется коэффициент jak68.wmf. Другое значение константы приводит к другому значению коэффициента a.

При постоянном среднем модуле тангенса угла наклона шероховатостей jak70.wmf, но разной высоте шероховатостей k, получается график отличный от графика Никурадзе.

iak.tif

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

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

jak71.wmf

Величина jak72.wmf. Влияния эффективного среднего модуля тангенса наклона шероховатостей на свойство потока зависит от числа Рейнольдса или перепада давления.

Эмпирическая формула по определению коэффициентов jak73.wmf следующая

jak74.wmf, jak75.wmf,

jak76.wmf.

При этом в начале образования мнимой части комплексного решения jak77.wmf, или начале турбулентного решения, тангенс наклона шероховатостей приближенно равен 2, и графики при разных тангенсах наклона шероховатостей совпадают.

При этом коэффициент сопротивления потока в круглом трубопроводе определяется по формуле jak78.wmf, число Рейнольдса, вычисленное на основании средней скорости течения потока, равно jak79.wmf. Коэффициент сопротивления на бесконечности давления пропорционален

jak80.wmf.

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

iak2.tif

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

Для сравнения теоретического и экспериментального графика зависимости коэффициента сопротивления от числа Рейнольдса потока приведен экспериментальный график Никурадзе справа на рис. 2. Ошибка теоретического графика относительно экспериментального примерно 10 %. Отмечу, что график, приведенный на рисунке Никурадзе в ламинарной части решения не точен, по мере приближения к критическому числу Рейнольдса появляется отличие от соотношения jak81.wmf, что отражено на теоретическом графике, и мало проявляется на экспериментальном. Экспериментальный график для логарифма давления и числа Рейнольдса изображен прямой линией, а должно быть отличие от прямой линии по мере увеличения влияния конвективного члена.

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

jak82.wmf

Причем формула построена таким образом, чтобы

jak83.wmf.

В случае ламинарного режима имеется простая формула для определения числа Рейнольдса jak84.wmf.

1.2. Алгоритм решения задачи гидродинамики при произвольной геометрии потока

Уравнение Навье – Стокса в декартовых координатах имеет вид

jak85.wmf. (1.2.1)

Для этого решаем трехмерную ламинарную стационарную задачу без учета конвективного члена при заданном внешнем воздействии gl

jak86.wmf.

Приведем эту задачу к безразмерному виду, разделив на величину jak87.wmf получим безразмерное уравнение

jak88.wmf

Решением этой задачи служит функция

jak89.wmf.

Строим решение уравнения неразрывности при внешнем воздействии, где величина jak90.wmf – отклик на внешнее воздействие

jak91.wmf. (1.2.2)

Откуда получаем уравнение по определению давления в потоке

jak92.wmf

Давление ищем в виде

jak93.wmf.

Подставляем его в подынтегральное выражение, умножаем на величину jak94.wmf, и интегрируем по пространству, получаем систему линейных уравнений

jak95.wmf,

где имеем значения коэффициентов

jak96.wmf

где величина jak97.wmf, определяется внешним воздействием. Приведем уравнение Навье – Стокса к безразмерному виду, разделив его на величину jak98.wmf, получим безразмерное уравнение

jak99.wmf.

Умножаем уравнение Навье – Стокса на сечение трубки тока, записываем уравнение вдоль ламинарного решения, вводя трубку тока, поток в которой постоянен см. [6].

jak100.wmf.

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

jak101.wmf, (1.2.3)

где величина Ss – сечение трубки тока в ламинарном режиме; величина jak102.wmf – стационарное решение уравнения Навье – Стокса без конвективного члена, который для трубки тока равен нулю, как не зависящий от продольной координаты.

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

Подставим решение (1.2.3) в уравнение Навье – Стокса и проинтегрируем по трубке тока, разделив на площадь сечения трубопровода. Тогда конвективный член имеет вид

jak105.wmf.

Введение шероховатости приводит к зависимости радиуса трубопровода jak106.wmf от степени макро-шероховатости. Далее выделяем член jak107.wmf связанный с шероховатостью, находим среднее значение его модуля. При этом произведем усреднение уравнения по величине s. Причем получается, что конвективный член в ламинарном режиме при гладкой поверхности, равен нулю и только при учете шероховатости он имеет значение. Получим уравнение

jak108.wmf.

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

Подставим решение (1.2.3) в уравнение Навье – Стокса и проинтегрируем по трубкам тока, умножим на jak111.wmf, где эта величина удовлетворяет условию jak112.wmf, где jak113.wmf средний модуль тангенса наклона неустранимых микро-шероховатостей, огибающая которых определяет макро-шероховатости, получим уравнение

jak114.wmf

jak115.wmf

где величины jak116.wmf определены из ламинарного решения и уравнения неразрывности, функция внешнего воздействия jak117.wmf задается. Получается, что микро-шероховатости, расположенные по всей длине трубопровода, определяют критическое число Рейнольдса. Эти микро-шероховатости меньше макро-шероховатостей, которые оказывают влияние на коэффициент сопротивления, при больших числах Рейнольдса. Но так как число Рейнольдса зависит от геометрии трубопровода, через его диаметр, то критическое число Рейнольдса только обратно пропорционально среднему модулю тангенсу наклона микро-шероховатостей, и зависит от геометрии трубопровода. При этом уменьшение радиуса трубопровода вызывает отрицательное значение jak118.wmf и, следовательно, отсутствие комплексного, турбулентного решения на размере сужения, т.е. повышается критическое число Рейнольдса. Наоборот расширение трубопровода вызывает увеличение величины jak119.wmf и, следовательно, уменьшение критического числа Рейнольдса и может вызвать более ранее комплексное решение, т.е. турбулентный режим. Причем, так как число Рейнольдса зависит от температуры через зависимость кинематической вязкости от температуры, получается, что наступление критического числа Рейнольдса находится в зависимости от температуры среды.

Координаты положения равновесия определятся из квадратного уравнения

jak120.wmf.

При этом ламинарное решение уточняется и имеет вид jak121.wmf, которое при малом перепаде давления переходит в линейное ламинарное решение jak122.wmf.

При этом справедлива турбулентная формула учета шероховатостей в силу одинакового метода усреднения в турбулентном режиме

jak123.wmf

где величина jak124.wmf – эффективный средний тангенс наклона шероховатостей; jak125.wmf – отношение высоты шероховатости к радиусу трубопровода. Критическое число Рейнольдса jak126.wmf удовлетворяет значению числа Рейнольдса, соответствующее началу комплексного решения. При этом при малом числе Рейнольдса получаем ламинарное решение. Но проблемы получения турбулентного решения на этом не кончаются. Нужно определить влияние шероховатостей поверхности, для чего пока необходимо использование экспериментальных данных. В принципе необходимо точное значение зависимости числа Рейнольдса для гладкой поверхности от макро-шероховатостей. Можно высказать предположение, что при счете с помощью ламинарного решения при подсчете турбулентного решения влияние шероховатости будет аналогично влиянию в случае круглого трубопровода с другими константами jak127.wmf, которые нужно определить из значения коэффициента сопротивления на бесконечности числа Рейнольдса. Для этого необходимо, чтобы размер, входящий в число Рейнольдса определялся аналогично в круглом трубопроводе и для тела произвольной формы. Для этого необходимо определять средний радиус поверхности через площадь поверхности тела и вводить радиус сферы, у которой площадь совпадает с площадью данного тела. Но у внешней задачи имеются свои особенности, связанные с наличием кризиса сопротивления, который связан с наличием следа за телом в потоке. Этого следа нет во внутренних задачах, таких как течение в трубопроводе.

Построим решение уравнения Навье – Стокса для внешней задачи. Имеется ламинарное решение для движения сферы в жидкости при малом числе Рейнольдса. Оно определяет следующее распределение скорости см. [8]

jak128.wmf

Отметим, что это ламинарное решение справедливо с точностью jak129.wmf, а не с точностью R, что связано с преобразованным видом уравнения Навье – Стокса

jak131.wmf.

Которое имеет стационарное решение jak132.wmf и нелинейным членом можно пренебречь при условии jak133.wmf. В самом деле, ламинарное решение этого уравнения jak134.wmf и нелинейный член равен jak135.wmf, что эквивалентно jak136.wmf. Последнее равенство выполняется в ламинарном режиме. Условие наступления комплексного, турбулентного решения является равенство нулю дискриминанта этого уравнения при числе Рейнольдса равном критическому.

При этом определится зависимость давления от параметров потока

jak137.wmf.

Определим скорость жидкости в декартовой системе координат

jak138.wmf

где jak139.wmf.

Направление линий тока определяется по формуле

jak140.wmf

Значение единичного вектора вдоль линии тока равно

jak141.wmf.

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

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

jak142.wmf,

где Us удовлетворяет (1.2.3). При этом Us первый интеграл дифференциального уравнения

jak143.wmf

не является первым интегралом дифференциального уравнения jak144.wmf, при этом выполняются условия jak145.wmf. При вычислении оператора Лапласа вдоль линии тока jak146.wmf используется направление и дифференциальное уравнение

jak147.wmf.

Задание параметров jak148.wmf на поверхности сферы определит кривую линии тока.

А при вычислении первого интеграла используется дифференциальное уравнение

jak149.wmf.

При этом скорость равна проекции ламинарных решений на линию тока

jak150.wmf, (1.2.3)

причем, направление jak151.wmf – это направление линий тока.

Решение уравнения Навье – Стокса ищем в полученном из предыдущих рассуждений виде

jak152.wmf

Эта формула определяет постоянное значение скорости вдоль траектории jak153.wmf, заданной jak154.wmf, равенство нулю нелинейного члена и в общем случае не равенство нулю второй производной от скорости вдоль траектории. Причем, дифференцирование осуществляется по величине s, и полагаем, вдоль траектории движения выполняется равенство jak156.wmf. Подставляем найденное значение скорости Us в уравнение Навье – Стокса, записанное вдоль линии тока получим

jak157.wmf

При этом вдоль направления линии тока jak158.wmf выполняется

jak159.wmf,

при этом величина jak160.wmf, так как вторая производная вдоль направления линии тока не равна нулю.

Интегрируя по величине s, получим нелинейное уравнение

jak162.wmf.

Это уравнение имеет стационарное решение

jak163.wmf.

Также существует комплексное решение при большом перепаде давления

jak164.wmf,

где величина g определяется по формуле

jak165.wmf.

Получается такая степень s, так как

jak166.wmf.

Если для внутренней задачи шероховатость уменьшает эффективное сечение потока, то для внешней задачи должно увеличивать эффективное сечение сопротивления потока.

jak167.wmf

Величина jak168.wmf. Величина k – средний модуль высоты шероховатостей. Влияния эффективного среднего модуля тангенса наклона шероховатостей на свойство потока зависит от числа Рейнольдса или перепада давления.

Эмпирическая формула по определению коэффициентов jak169.wmf следующая

jak170.wmf,

jak171.wmf,

jak172.wmf.

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

1.3. Построение решения в виде ряда для трубопровода с круговым сечением

Можно привести уравнение Навье – Стокса к безразмерному виду для ряда, являющегося решением

jak173.wmf

Запишем уравнение Навье – Стокса в безразмерном виде, умножая на величину jak174.wmf, и, умножая уравнение неразрывности на величину jak175.wmf, получим уравнения

jak176.wmf

Имеем значение переменных

jak177.wmf, jak178.wmf,

jak179.wmf

jak180.wmf,

где величина a – радиус трубопровода, jak181.wmf – плотность и кинематическая вязкость жидкости, jak182.wmf – скорость и давление в потоке. Величина давления равна

jak183.wmf,

где jak184.wmf – начальное и конечное давление в потоке. Имеется связь

jak185.wmf,

которая является определением давления jak186.wmf. Величина скорости определяется по формуле

jak187.wmf.

Величина a(z) определяет высоту шероховатости, при постоянном радиусе трубопровода. Уравнение неразрывности удовлетворяется тождественно с точностью

jak188.wmf,

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

jak189.wmf.

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

Удалось решить уравнение Навье – Стокса с точность до 0.6 % при постоянном среднеквадратичном значении тангенса наклона шероховатостей. Строилось решение той же задачи на основе алгоритма (1.3.1). Получился следующий график (рис. 3).

iak3.tif

Рис. 3

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

jak190.wmf (1.3.1)

где величина jak191.wmf определяется по формуле (1.3.2). При этом, так как учитывается множество членов, отражающих разную степень завихренности, корень из мнимой части комплексного числа Рейнольдса потока не извлекается. Степень шероховатости учитывалась из условия одинаковости коэффициента сопротивления при разных формулах для решения, изменялся только коэффициент a, так как на бесконечности числа Рейнольдса величина X является константой, определяемой формулой численного метода и необходимо, чтобы

jak192.wmf

откуда определится коэффициент jak193.wmf при условии средней гладкости jak194.wmf. Здесьвеличина jak195.wmf является корнем уравнения

jak196.wmf

jak197.wmf

Решение этой системы нелинейных уравнений ищется в виде jak198.wmf, причем для коэффициента jak199.wmf получаем квадратное уравнение

jak200.wmf (1.3.2)

При этом величина jak201.wmf определяется по формуле jak202.wmf.

Решение этой алгебраической системы уравнений имеет N+1 корень. Докажем это. Допустим, нашлось одно решение этой системы уравнений jak203.wmf. Алгоритм его построения описан выше по тексту. Составим решение jak204.wmf. Подставим это решение в уравнение по определению координат положения равновесия, получим уравнение

jak205.wmf.

Получим уравнение

jak206.wmf

Выберем начальное приближение так, чтобы определитель этой линейной системы уравнений равнялся нулю. Для этого определим jak207.wmf, и найдем из равенства нулю определителя значение a. Тогда следующее приближение определит значения jak209.wmf с точностью до множителя. Вычислим этот множитель из условия jak210.wmf. Это можно сделать N способами, определяя корень полинома N степени. В результате получим n+1 решение алгебраического уравнения.

Получение счетного количества (jak211.wmf) координат положения равновесия эквивалентно получению счетного количества решений для уравнения Шредингера. Эквивалентность уравнения Шредингера и Навье – Стокса докажем в данной статье далее по тексту.

При этом основное решение jak212.wmf определится однозначным образом. Это доказывается алгоритмом его построения и невозможности выделить из решений jak213.wmf другую константу, удовлетворяющую уравнению по определению координат положения равновесия. Отметим, так как коэффициенты jak214.wmf меньше, чем линейный член при малом перепаде давления, значение jak215.wmf соответствует линейной части решения. Продолжая это решение при большом перепаде давления, получим основной член решения задачи гидродинамики.

График на рис.3 отличается от графика Никурадзе на рис.4. Асимптотика при больших и малых значениях числа Рейнольдса совпадает с графиком Никурадзе. Но график построен при постоянном тангенсе наклона шероховатости, а график Никурадзе при постоянном отношении радиуса трубопровода к средней высоте шероховатости. Построим значение эффективного тангенса наклона шероховатостей, по средней высоте шероховатости.

jak216.wmf

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

jak219.wmf, jak220.wmf, jak221.wmf jak222.wmf,

где jak223.wmf – коэффициенты a, образующие множитель при мнимой части решения для решения в виде одного члена и в виде ряда соответственно.

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

iak4.tif

Рис. 4

График построен для отношения радиуса трубопровода r0 к высоте шероховатостей jak225.wmf равной jak226.wmf. Точность аппроксимации относительно графиков Никурадзе примерно 10 %.

2. Построение решения обыкновенных дифференциальных уравнений

Рассмотрим систему нелинейных автономных дифференциальных уравнений

jak227.wmf. (2.1)

К системе нелинейных дифференциальных уравнений сводится система уравнений Навье – Стокса и уравнение неразрывности.

jak228.wmf, (2.2)

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

jak229.wmf.

Причем функция jak230.wmf определяется в виде синуса. Тогда коэффициенты jak231.wmf для непрерывной функции убывают с ростом индекса как величины не большие jak232.wmf и возможна редукция, т.е. вместо бесконечного числа членов используется конечное число членов.

Исследуются комплексные и действительные решения задачи Коши этого дифференциального уравнения в случае действительных и комплексных начальных условиях, при действительном аргументе t. Начальные условия имеют вид jak233.wmf, где величина t0 соответствует начальному моменту интегрирования, а величина jak234.wmf в общем случае комплексная. Причем в случае действительных значениях jak235.wmf, правая часть (2.1) действительна.

Систему дифференциальных уравнений (2.2) можно представить при не кратных положениях равновесия путем подстановки jak236.wmf. При этом положения равновесия системы (2.2) jak237.wmf перейдут в положения равновесия jak238.wmf. При этом определятся собственные числа и собственные векторы линеаризованной системы (2.2).

jak239.wmf

Система уравнений (2.2) запишется в виде

jak240.wmf. (2.3а)

Систему уравнений (2.3а) можно записать в виде

jak241.wmf, (2.3б)

где введен не обращающийся в ноль множитель jak242.wmf, который равен

jak243.wmf.

При подстановке этого множителя в (2.3б), получим (2.3а). Покажем, что этот множитель в ноль не обращается. Величины jak244.wmf удовлетворяют условию jak245.wmf, где величина S конечна.

При условии jak246.wmf имеем конечный предел

jak247.wmf

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

jak250.wmf,

где величина S конечна.

При этом дифференциальное уравнение можно записать в виде

jak251.wmf, (2.4)

где jak252.wmf – стремящаяся к бесконечности функция при условии стремления координат к положению равновесия. В случае решения в действительной плоскости это монотонная функция.

Лемма 2. Необходимым и достаточным условием стремления неизвестной функции к устойчивым координатам положения равновесия является условие jak253.wmf, причем jak254.wmf.

Имеем соотношения

jak255.wmf (2.5)

при jak256.wmf и значит jak257.wmf, как интеграл от константы. Справедлива и обратная теорема, при условии jak258.wmf, реализуется одно из устойчивых положений равновесия. Это следует из вида решения, при условии jak259.wmf имеется отрицательная действительная часть у числа jak260.wmf в формуле (2.6) согласно лемме 4, и решение стремится к координате положения равновесия jak261.wmf в формуле (2.4). При этом величина времени стремится к бесконечности.

Лемма 3. Решением дифференциального уравнения (2.1) является функция jak262.wmf, удовлетворяющая формуле (2.6).

Для получения (2.6) разделим уравнение (2.4) на произведение множителей jak263.wmf и умножим (2.4) на величину jak264.wmf. Раскладываем полученную дробь на сумму простых дробей и их интегрируем. Получаем уравнения

jak265.wmf.

Здесь в случае излучения звуковой энергии в интервале jak266.wmf получим разные ветви логарифма.

Потенцируя полученное выражение, получим (2.6)

jak267.wmf (2.6)

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

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

jak269.wmf,

где если величины jak270.wmf не равны нулю, то они имеют множитель jak271.wmf, иначе уравнение не имеет координат положения равновесия. Причем функция jak272.wmf, имеющая S совокупностей координат положения равновесия должна обязательно иметь множители с одним или несколькими n, входящих в произведение jak275.wmf.

Тогда нет определенной ветви вычисления координаты положения равновесия и под координатой jak276.wmf можно понимать любую из координат jak277.wmf. Т.е. совокупности координат перепутаны и могут состоять из любой комбинации и тогда обязательно имеются устойчивые положения равновесия.

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

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

Лемма 4. Сумма коэффициентов jak278.wmf по индексу s равна нулю, т.е. jak280.wmf.

Для доказательства этого тождества рассмотрим полином S–1 степени относительно y

jak281.wmf,

В точках положения равновесия jak282.wmf полином удовлетворяет jak283.wmf. В силу единственности полинома степени S–1, проходящего через S точек, получаем jak284.wmf, так как это значение удовлетворяет точкам аппроксимации. Распишем формулу для полинома, равного единице, разделив его на произведение jak285.wmf, получим

jak286.wmf

полагая, jak287.wmf получим тождество jak288.wmf, в случае, если имеется S+1 положение равновесия.

В случае если разлагается дробь

jak289.wmf,

где jak290.wmf полином степени S–1. Значение коэффициентов изменится, но свойство jak291.wmf останется, jak292.wmf. Докажем это. Для чего рассмотрим сумму

jak293.wmf.

Эта сумма равна jak294.wmf. Распишем формулу для полинома, равного jak295.wmf, разделив его на произведение jak296.wmf, получим

jak297.wmf

Полагая jak298.wmf, получим тождество jak299.wmf, в случае, если имеется S+1 положение равновесия.

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

Теорема 1. Рассматривается задача Коши при произвольных действительных начальных условиях для системы нормальных нелинейных обыкновенных дифференциальных уравнений (2.1). Случай вырожденного решения задачи Коши – положения равновесия, не рассматривается. В случае, если у системы (2.1) имеются комплексно-сопряженные положения равновесия с действительной частью, то при конечном аргументе t действительное решение задачи Коши системы (2.1) при действительных начальных условиях стремится к бесконечности. Потом это решение переходит в комплексное решение, стремясь к положению равновесия, в случае, если комплексные координаты положения равновесия имеют действительную часть. При этом правую часть (2.1) считаем регулярной функцией, действительной при действительных аргументах. Она имеет конечное число не кратных положений равновесия.

Доказательство.

Если решать систему (2.2) при не кратных положениях равновесия, то получим согласно с леммой 3

jak301.wmf (2.7)

где jak302.wmf выделенное комплексное положение равновесия, jak303.wmf остальные положения равновесия. Коэффициенты jak304.wmf удовлетворяют jak305.wmf в соответствии с леммой 4. При этом в сумме jak306.wmf величина действительной части jak307.wmf в случае комплексного значения jak308.wmf участвует дважды и в силу того, что все числа jak309.wmf удовлетворяют условию jak310.wmf, имеем формулу jak311.wmf.

Обоснуем формулу (2.3). Для этого два комплексно сопряженных члена решения преобразуем (для упрощения записи индекс l опускаем)

jak313.wmf, (2.8)

где jak314.wmf. После интегрирования (2.8) по аргументу x, получим формулу (2.7)

jak316.wmf.

Решение равняется

jak317.wmf,

где

jak318.wmf

При этом величина jak319.wmf действительная в силу существования комплексно-сопряженных положений равновесия. Т.е. имеем равенство при условии jak320.wmf при конечном t

jak322.wmf. (2.9)

Это уравнение имеет решение, стремящееся к бесконечности.

При этом решение дифференциального уравнения при росте jak323.wmf согласно лемме 2, может иметь комплексные корни

jak324.wmf.

При этом, так как справедливо jak325.wmf согласно лемме 4, и положения равновесия имеют действительную часть, имеются числа с отрицательной действительной частью jak326.wmf, значит, имеется сходимость к одному из положений равновесия. Действительное решение будет стремиться к бесконечности, причем нарушатся условия существования и единственности задачи Коши. При этом при бесконечности jak327.wmf согласно лемме 2 неизвестная функция будет стремиться к одному из положений равновесия. Это положение равновесия не может быть действительным, так как действительное решение бесконечно. Значит, решение будет иметь точку ветвления и стремиться к комплексному положению равновесия. Значит, при комплексных положениях равновесия получается конечное комплексное решение при изменении jak328.wmf. Т.е. в некоторой точке начнется комплексное решение.

Конец доказательства.

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

jak329.wmf.

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

Используя неявную схему решения, получим следующее уравнение

jak332.wmf.

Разрешая относительно неизвестной функции x, получим неявную схему

jak334.wmf.

Эта неявная схема с постоянным шагом правильно описывает стремление решения к бесконечности. При счете с уменьшенным шагом она определяет большее значение переменной t и значит, определяет большее значение неизвестной функции. Т.е. правильно описывает решение дифференциального уравнения до бесконечности решения. Когда бесконечность достигнута, при условии jak336.wmf определится конечное комплексное решение. Численный счет этого уравнения подтвердил правильность проведенного анализа решения.

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

Решение с комплексными начальными данными определится формулой

jak337.wmf

при любом t. Т.е. приближенно имеем

jak339.wmf

При этом знаменатель этой дроби в ноль не обращается.

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

Но возникает вопрос, каков физический смысл мнимой части комплексного решения.

3. Физический смысл комплексного решения

Опишем физический смысл комплексного турбулентного решения. Итак, рассмотрим действительное решение системы обыкновенных дифференциальных уравнений jak340.wmf. Пусть начальные данные имеют среднее jak341.wmf и дисперсию jak342.wmf (дисперсия начальных данных в случае уравнения Навье – Стокса определяется шероховатостью поверхности или не точно заданными начальными данными). Тогда для дисперсии решения имеем

jak343.wmf.

Значит имеем

jak344.wmf. (3.1)

Приведу формулировку обратной теоремы Пифагора. Для всякой тройки положительных чисел a,b и c, такой, что jak345.wmf, существует прямоугольный треугольник с катетами a и b и гипотенузой c. Значит, математическое ожидание и среднеквадратическое отклонение образуют катеты, а гипотенузой является корень из среднего квадрата величины. Т.е величина среднего jak346.wmf ортогональна среднеквадратическому отклонению jak347.wmf, которое образует мнимую часть координаты тела. Таким образом, полученное в результате усреднения во времени декартово пространство с колебательной скоростью высокой частоты (период колебания меньше времени измерения) становится комплексным пространством. Т.е. в случае большой дисперсии величины действительного пространства, его нужно рассматривать как комплексное трехмерное пространство, где мнимая часть соответствует среднеквадратическому отклонению. При этом имеется следующая связь между переменными

jak348.wmf,

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

jak350.wmf,

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

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

Отметим, что при расчете движения потока при учете одного члена ряда, надо брать квадратный корень из мнимой части, так как вычисляется поступательная скорость. Мнимая часть соответствует квадратному корню из колебательной части безразмерной скорости. Ситуация аналогична вычислению отклонения при случайном выборе шага вперед или назад с вероятностью ½ и положение точки после N шагов определяется значением jak351.wmf. Действительная часть решения и мнимая часть решения расположены на разных осях комплексного пространства. Но если усреднять мнимую безразмерную часть, получим

jak352.wmf

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

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

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

Лемма 5. Комплексное решение определяет колеблющуюся пульсирующую функцию координат перемещения потока.

Мнимая часть скорости соответствует скорости вращения в фазовом пространстве. Так как известен радиус вращения, то можно определить и частоту вращения. В плоскости вращения комплексную скорость с постоянным радиусом вращения и постоянной частотой можно представить в виде jak354.wmf. В случае переменной по пространству стационарной скорости эту формулу можно представить локально в одной плоскости в виде

jak355.wmf

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

Лемма 6. Трехмерную скорость потока можно представить в виде

jak356.wmf

Причем, скорости определяются в виде интеграла от касательного ускорения, по формуле

jak357.wmf

Интеграл от нормального ускорения определяет нормальную компоненту скорости, по формуле

jak358.wmf.

При этом величина локальной скорости jak359.wmf. Но проинтегрированная относительно центростремительного ускорения скорость отлична от нуля jak360.wmf, обращаясь при постоянной скорости частицы и постоянном радиусе кривизны, за период

jak361.wmf,

где величина R – радиус кривизны, в ноль при той же помеченной начальной точке. При переменной скорости частицы за время, когда один из интегралов jak362.wmf, которое, при конечном радиусе кривизны одного знака траектории, конечно и равно

jak363.wmf,

так как касательное направление jak364.wmf, при вращении меняет знак.

При этом вклад в поступательную часть комплексной скорости за один оборот вращения равен

jak365.wmf.

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

jak366.wmf.

Причем имеем значение нормальной компоненты скорости jak367.wmf, где величина jak368.wmf соответствует замкнутой траектории скорости вращения, а величина jak369.wmf вклада в поступательную скорость. Откуда имеем величину дополнительного вклада в поступательную скорость за счет мнимой части числа Рейнольдса jak370.wmf. Откуда имеем дополнительный вклад в поступательную скорость jak371.wmf.

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

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

jak373.wmf.

Направление скоростей jak374.wmf ортогонально и их сумма приводит к приращению модуля скорости движения

jak375.wmf,

так как

jak376.wmf.

Дифференцируемые по времени компоненты этих проекций определяют тангенциальное и нормальное ускорение. При этом вводится понятие тангенциальной и нормальной скорости, которые в декартовом пространстве не ортогональны jak377.wmf, но в шестимерном комплексном пространстве ортогональны и их модуль комплексного вектора jak378.wmf равен

jak379.wmf.

Это доказывается представлением jak380.wmf и вычислением модуля как произведения комплексно сопряженных векторов с учетом ортогональности шести действительных ортов.

Итак, получено решение уравнений Навье – Стокса в случае не кратных положений равновесия. Оно имеет вид

jak381.wmf.

jak382.wmf,

где величины jak383.wmf координаты положения равновесия.

Ламинарное решение соответствует координатам положения равновесия, структура турбулентного решения имеет вид

jak384.wmf,

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

При этом формула

jak387.wmf. (3.2)

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

4. Построение решения в случае кратных положений равновесия

Но как же строить решение в случае кратных координат положения равновесия. При этом одно из собственных чисел равно нулю и в формуле (2.5) собственное число равно нулю и значит, hl может не стремиться к бесконечности, что нарушает теорему о сходимости к комплексному положению равновесия. Для построения решения докажем теоремы 2, 3.

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

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

Рассмотрим систему нелинейных автономных дифференциальных уравнений, имеющих двукратный корень

jak388.wmf. (4.1)

Осуществим вырожденное преобразование системы дифференциальных уравнений с кратным корнем (кратный корень является кратной координатой положения равновесия) с помощью неизвестных функций jak389.wmf. Подставим в дифференциальное уравнение, получим

jak390.wmf. (4.2)

Будем рассматривать для l уравнения координату jak392.wmf независимой, а остальные координаты как функции времени, являющиеся решением дифференциального уравнения. Система уравнений разбивается на N независимых неавтономных уравнений, которые в случае, когда решение равно координатам положения равновесия, правая часть дифференциального уравнения не зависит от времени. При этом связь между решениями относительно преобразования jak394.wmf одинакова со связью между положениями равновесия (связь вырожденна в кратном положении равновесия, но при этом для остальных произвольных значениях xl получается единственное значение yl). При этом справедливо следующее равенство

jak395.wmf.

Т.е. допустим, матрица уравнения (4.1) jak396.wmf не вырождена, но система (4.1) имеет двукратное положение равновесия. Тогда матрица уравнения (4.2), полученного из уравнения (4.1), с помощью преобразования jak397.wmf равна (допустим кратный корень имеет последнее уравнение системы (4.1), тогда при дифференцировании по jak398.wmf в точке двукратного положения равновесия останется только последнее уравнение, остальные уравнения имеют неопределенность, которая равна константе и при дифференцировании дают ноль)

jak399.wmf.

При этом первый член второй строки этой формулы является константой в точке положения равновесия, так как

jak400.wmf

и выбором преобразования jak401.wmf можно добиться, что матрица jak402.wmf может быть вырождена в точке кратных координат положения равновесия системы (3.1).

Но преобразование переменных jak403.wmf не произвольно. Преобразование таково, что действительным или комплексным корням jak404.wmf должны соответствовать действительные или комплексные корни jak405.wmf. Строить систему (4.2) необходимо, так как уравнение (4.2) имеет не кратные положения равновесия, соответствующие положениям равновесия системы (4.1). По координатам системы (4.2) однозначно определяются соответствующие положения равновесия системы (4.1). При этом добиваемся, чтобы определитель линеаризованной системы уравнений (4.2) равнялся нулю.

При условии, что содержится кратный корень jak406.wmf, являющийся одинаковым решением s и p ветви решения алгебраического уравнения, эквивалентное N уравнение для переменной xN имеет вид

jak408.wmf.

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

jak409.wmf.

При этом эквивалентный определитель системы дифференциальных уравнений равен нулю. Но формула для линеаризованной системы уравнений учитывает одну ветвь координат положения равновесия, упуская нулевое значение для определителя jak410.wmf в случае кратных положений равновесия при приближении xN к значениям jak411.wmf или jak412.wmf. При линеаризации системы нелинейных уравнений учитывается либо jak413.wmf, либо jak414.wmf, а надо учитывать обе одновременно.

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

jak415.wmf

Эта система имеет положения равновесия

jak416.wmf

При этом корень jak417.wmf двукратный. Причем, определитель линеаризованной системы этих уравнений не вырожденный. Эта задача в хаотическом режиме вращается вокруг одного из положений равновесия. Потом решение системы нелинейных уравнений резко перескакивает к другому положению равновесия и вращается вокруг второго положения равновесия. Это связано с тем, что хотя положения равновесия кратные, определитель линеаризованной системы уравнений у аттрактора Лоренца не равен нулю. Но у системы линеаризованных уравнений, построенных с помощью двукратного корня по мере приближения к точке, соответствующей jak418.wmf, в случае jak419.wmf согласно теореме 2 образует скачок от этого положения равновесия. При этом происходит произвольный перескок на решение с положением равновесия

jak420.wmf.

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

В книге [6], описан сценарий рождения аттрактора Лоренца через неполный двойной гомоклинический каскад бифуркаций. В этой книге считается, что условием хаотического решения является наличие в системе Лоренца седло – узла и двух седло – фокусов, откуда вытекает возможность существования в ней различных гомоклинических и гетероклинических контуров особых точек и связанных с ними каскадов бифуркаций. Терминологию и объяснение обозначений см. [6].

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

Теорема 2. Пусть jak421.wmf координаты одного из положений равновесия системы (2.1). В случае двукратного значения положения равновесия jak422.wmf, линейное приближение решения дифференциального уравнения с устойчивым положением равновесия jak423.wmf, определяет сходимость к положению равновесия при условии jak424.wmf, где матрица jak425.wmf, это собственные векторы матрицы линеаризованной системы. В случае противоположного условия jak426.wmf положение равновесия не достижимо и получается многозначное решение. Предполагается, что выполняется jak427.wmf.

Доказательство.

В силу существования двукратного корня, равного jak428.wmf уравнение (2.1) можно привести к виду (4.3) с вырожденной матрицей линеаризованной системы

jak429.wmf (4.3)

Т.е. определитель матрицы jak430.wmf равен нулю в этом представлении.

Решаем линеаризованное уравнение (4.4)

jak431.wmf (4.4)

Это уравнение имеет решение в окрестности положения равновесия

jak432.wmf,

так как jak433.wmf, где собственные векторы jak434.wmf и собственные числа jak435.wmf определятся из системы уравнений

jak436.wmf

Причем наблюдается приближение jak437.wmf к значению jak438.wmf. Т.е. положение равновесия не обязательно достижимо в случае кратных корней, даже для системы с jak439.wmf, так как выполняется условие jak440.wmf. Найдем условие сходимости к кратному корню. Имеем соотношение

jak441.wmf.

При этом для определения константы cN при дискретном вычислении решения имеем следующее рекуррентное соотношение jak442.wmf при достижении устойчивыми координатами положения равновесия. Т.е условием сходимости решения в направлении jak443.wmf является неравенство jak444.wmf. При условии jak445.wmf наблюдается отсутствие сходимости в одном направлении jak446.wmf, остальные cl не растут. Скачок решения осуществляется мгновенно на произвольное значение величины cN, по формуле

jak447.wmf,

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

jak448.wmf. (4.5)

Но в силу приближенности решения, получится не бесконечность, а переход к другому положению равновесия. Точное решение при численном счете получается мгновенно при условии jak449.wmf и имеет значение, зависящее от малости величины h. Т.е. происходит произвольный скачок. Остальные возможно комплексные направления jak450.wmfсобственного вектора, соответствуют устойчивому собственному числу положения равновесия.

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

Если же наблюдается одно положительное собственное число при остальных отрицательных, то имеем

jak452.wmf

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

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

Конец доказательства.

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

Теорема 3. В случае системы (4.1) с двукратным корнем положения равновесия решение задачи Коши в комплексной плоскости для системы дифференциальных уравнений (4.1) с действительными и комплексными положениями равновесия будет определяться рядом Лорана с существенно особой точкой и, следовательно, при приближении к особенности может носить хаотический характер. Решение может иметь полюсы со знаменателем

jak453.wmf,

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

Доказательство.

В случае двукратного корня aN систему дифференциальных уравнений (2.1) можно представить в виде

jak454.wmf, (4.6)

где введен не обращающийся в ноль множитель jak455.wmf, который равен

jak456.wmf.

При подстановке этого множителя в (4.6), получим (4.2). Покажем, что этот множитель в ноль не обращается. Величины jak457.wmf удовлетворяют условию

jak458.wmf,

где величина K конечна.

При условии jak460.wmf имеем конечный предел

jak461.wmf

При условии jak462.wmf тоже имеем конечный предел

jak463.wmf.

В случае совпадающих корней сокращается множитель jak464.wmf. При этом получается не обращающийся в ноль множитель jak465.wmf.

При этом дифференциальное уравнение можно записать в виде

jak466.wmf, (4.7)

где jak467.wmf – стремящаяся к бесконечности функция при условии приближения к положениям равновесия. В случае решения в действительной плоскости это монотонная функция.

Запишем дифференциальное уравнение с кратным значением положения равновесия

jak468.wmf

Величина jak469.wmf в ноль не обращается. Решение этого дифференциального уравнения имеет вид

jak470.wmf (4.8)

jak471.wmf – координаты положения равновесия, величина jak472.wmf – значение двух кратного положения равновесия. Вычислим коэффициенты этого разложения

jak473.wmf

jak474.wmf

Имеем три возможных варианта разложения функции в ряд. Либо в точке jak475.wmf имеется существенная особая точка, и решение стремится к произвольным точкам при приближении к особенности, либо справедливо

jak476.wmf (4.9)

либо решение регулярно. При этом величина jak477.wmf конечна. Допустим справедливо (4.9). Определим решение этого дифференциального уравнения, для чего подставим часть решения в виде полюса jak478.wmf в уравнение

jak479.wmf. (4.10)

При этом оказывается, что приравниваются члены

jak480.wmf,

что приводит к равенству jak481.wmf и для множителя справедливо jak482.wmf. При этом решение при условии jak483.wmf имеет сложную структуру и не имеет однозначного решения в комплексной плоскости, т.е. в точке jak484.wmf имеется ветвление решения и значит, возможен переход к комплексному решению. Имеем jak485.wmf. Т.е. решение в этом случае является турбулентным, комплексным и имеющим счетное количество ветвей.

Допустим, решение регулярно. Разложим решение в окрестности точки jak486.wmf по положительным степеням jak487.wmf, т.е. jak488.wmf. Тогда решение представимо в виде

jak489.wmf.

Т.е. справедлива формула

jak490.wmf.

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

При этом решение с существенной особой точкой определяется из уравнения (4.6). Т.е. по мере приближения к положению равновесия jak493.wmf координаты jak494.wmf происходит скачок решения, который является произвольным, так как решение в виде ряда Лорана содержит существенную особую точку и при приближении к особой точке, значение функции jak495.wmf произвольно. При достижении положения равновесия остальными координатами jak496.wmf, в случае, если положения равновесия кратные, наблюдается скачок решения на сходимость к другому положению равновесия. В результате решения в окрестности кратного положения равновесия, получаются ряды Лорана с существенной особой точкой.

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

При этом многозначное решение имеет вид ряда Лорана и может иметь существенно особую точку, а может быть многозначное решение имеет вид

jak497.wmf, (4.11)

где величина t1 определится из начальных условий, а величины tk имеют постоянное значение. В случае ряда Лорана реализуется хаотическое решение, которое соответствует случаю jak498.wmf теоремы 2 и происходит отскок от положения равновесия. В случае решения вида (4.9) реализуется случай jak499.wmf, и положение равновесия достигается в комплексной плоскости. В первом случае решение является многозначной функцией. Во втором случае наблюдается стремление к положению равновесия с помощью многозначной функции при бесконечном значении V.

Конец доказательства.

В случае двукратного положения равновесия по переменной xN решение уравнения Навье – Стокса имеет вид

jak500.wmf

jak501.wmf

где величины jak502.wmf однократные координаты положения равновесия, а величины jak503.wmf двукратные положения равновесия.

Ситуация аналогична уравнению Шредингера, когда имеется в общем случае счетное количество решений. Это не удивительно, ведь уравнение Шредингера сводится к уравнению Навье – Стокса. Докажем это. Для чего запишем уравнение Шредингера и преобразуем

jak504.wmf.

Разделив на массу jak505.wmf, получим уравнение

jak506.wmf.

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

jak508.wmf.

Подставляя значение скорости в преобразованное уравнение Шредингера, получим

jak509.wmf.

Получим трехмерное уравнение Навье – Стокса с давлением, соответствующим потенциалу. Но задача гидродинамики отличается от уравнения Навье – Стокса, полученного из уравнения Шредингера, уравнением неразрывности.

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

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

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