science-review.ru
1. Учет ламинарного решения уравнения Навье – Стокса
Значение коэффициента сопротивления круглого трубопровода при произвольном числе Рейнольдса и степени шероховатости известно только из эксперимента. Предлагается с помощью комплексного решения получить решение уравнения Навье – Стокса и на основе качественных соображений определить влияние шероховатости на решение уравнения Навье – Стокса. При этом с точностью 10 % удалось построить классические графики Никурадзе зависимости коэффициента сопротивления круглого трубопровода в зависимости от произвольного числа Рейнольдса и степени шероховатости.
Проблема описания движения жидкости в турбулентном режиме не решена. Это создает трудности при расчете нефтяных, газовых трубопроводов. Кроме того, не существует теоретических методов описания движения тела в турбулентной среде. Это и описание движения летательного аппарата, движения подводной лодки или надводного корабля в турбулентном режиме. Без использования моделирования движущегося тела в аэродинамических трубах, или водных бассейнах, невозможно проектирование движущихся в вязкой среде тел. Существуют приближенные формулы аппроксимации коэффициента сопротивления трубопровода на отдельных участках числа Рейнольдса см. [1, 2]. Известны классические экспериментальные графики Никурадзе по определению коэффициента сопротивления круглого трубопровода в зависимости от числа Рейнольдса и степени шероховатости. Применяется аппроксимация конвективного члена, сводящая задачу решения уравнения Навье – Стокса к линейной с эффективной турбулентной вязкостью. Но такое преобразование искажает решение уравнения Навье – Стокса и для совпадения с экспериментом коэффициент турбулентной вязкости может иметь произвольное значение, вплоть до отрицательного. Применяется метод Галеркина, сводящий решение задачи гидродинамики к системе нелинейных обыкновенных дифференциальных уравнений. Но при турбулентном режиме эта система нелинейных уравнений имеет комплексные положения равновесия, т.е. решение комплексное. И действительно в турбулентном режиме система уравнений гидродинамики в действительной плоскости не имеет решения, решение уравнения стремится к бесконечности см. [3]. При этом комплексное решение конечно. Физический смысл комплексного решения и его колебательный характер мнимой части см. [4, 5]. Т.е. решать задачу гидродинамики в турбулентном режиме надо в комплексной плоскости. При этом турбулентное решение определяется не однозначно, а имеется счетное количество ветвей решения.
1.1. Расчет трубопровода с круглым сечением в случае несжимаемой жидкости
Реализацией данного алгоритма решение является расчет трубопровода с круглым сечением. Решение задачи о трубопроводе с круглым сечением, будем искать в виде 
в цилиндрической системе координат. Так как внешнее воздействие имеется только для продольной оси
,
где 
давление в начальной и конечной части трубопровода, величина L это длина трубопровода, радиальной и угловой скоростью пренебрегаем. Внешнее воздействие равно
.
Согласно формуле (1.2.2) градиент давления равен
.
Получим уравнение
.
Подставляем значение скорости, получим уравнение
.
Умножаем это уравнение на радиус и интегрируем это уравнение по радиусу, так как система координат цилиндрическая, получим
.
Для получения счетного количества решений, умножаем это уравнение на функцию 
и интегрируем по пространству. Получаем счетное количество решений, как для гладкой, так и шероховатой поверхности. При этом, как и у уравнения Шредингера получается счетное количество решений, каждое со своей собственной энергией. При переходе из одного состояния в другое излучается дискретная энергия. Минимум собственной энергии определяет выбор решения.
Беря модуль правой части, и находя средний модуль тангенса угла отклонения, получим
. (1.1.1)
В случае выбора знака минус у значения среднего модуля тангенса угла отклонения 
, получится, что шероховатости увеличивают скорость потока, так как полная производная 
увеличится, что не правильно, скорость потока должна уменьшаться за счет шероховатостей.
При ведении турбулентной вязкости используется отрицательный знак у средней скорости, связанной с корреляционной функцией скорости процесса 
см. [1], что приводит к знаку плюс у среднего по модулю тангенса наклона шероховатости. Уравнение движения с учетом возмущения имеет вид
.
Т.е. конвективный член нужно брать с учетом минуса, причем правую часть (1.1.1) брать со знаком плюс. Кроме того, надо выбирать знак плюс у среднего тангенса модуля наклона шероховатости, чтобы получилось комплексное турбулентное решение. В противном случае, решение, описывающее пульсирующий турбулентный режим не получится.
Переходя от радиуса к диаметру трубопровода и разделив на величину 
, получим
. (1.1.2)
Если микро-шероховатости 
распределены по всей поверхности трубопровода, они находятся и на макро шероховатостях и определяют критическое число Рейнольдса и коэффициент сопротивления при числе Рейнольдса, равном 2300. Микро-шероховатости имеют молекулярную природу и определяются средним размером атома, равным среднему геометрическому между размером ядра rA, и размером орбиты Бора 
, при расстоянии между атомами 
, равному некоторой величине, определяемой свойствами границы трубопровода, железом, титаном и углеродом. Расстояние между атомами железа 
, между атомами титана 
, между атомами углерода 
см. [7]. При этом абсолютная величина тангенса наклона высоты микро-шероховатости поверхности металла в трубопроводе определяется по формуле
.
Средний тангенс наклона равен
Величина критического числа Рейнольдса относительно диаметра равна 
. Но почему критическое число Рейнольдса для сферы равно 3⋅105. Это связано с разным определением критического числа Рейнольдса. Оно равно
,
где величина 
эффективный, гидродинамический размер тела, включая среду, a истинный геометрический размер тела, причем 
молекулярный тангенс наклона шероховатости. Причем отношение 
может иметь значение 
.
Величина критического числа Рейнольдса равна 
. Макро-шероховатости 
более редкие и определяют коэффициент сопротивления на числах Рейнольдса в 12 и более раз больше.
Получаем условие стационарности для уравнения Навье – Стокса с учетом одного члена ряда-решения в одномерном случае
.
В одномерном случае при постоянстве сечения трубопровода уравнение неразрывности выполняется тождественно. Ламинарное решение этого уравнения равно
.
При внешнем давлении, равном 
начинается комплексное решение и турбулентный режим, так как число Рейнольдса в этой точке равно критическому значению. Из эксперимента и путем проделанного вычисления имеем значение критического числа Рейнольдса для круглого трубопровода 
. Коэффициент сопротивления трубопровода с круглым сечением определяется по формуле (в формулу подставляем перепад давления, выраженный через безразмерное давление)
Средняя скорость, входящая в число Рейнольдса, равна
.
Асимптотика коэффициента сопротивления трубопровода 
с круговым сечением для ламинарного режима вычислена верно.
.
Асимптотика получена при малом числе Рейнольдса, когда конвективный член мал.
В случае большого перепада давления, получаем комплексное турбулентное решение 
. Если считать точнее, то вклад вращательной мнимой части в поступательную скорость движения потока соответствует корню из мнимой части согласно формуле (1.1.3)
(1.1.3)
причем необходимо использовать величину отношения числа Рейнольдса на корень из безразмерного давления, как величину порядка единицы в турбулентном режиме. Число Рейнольдса потока пропорционально на бесконечности давления 
. При этом степень усреднения пропорциональна высоте шероховатости, стремящейся к нулю, хотя самая гладкая поверхность соответствует среднему модулю тангенса наклона, равному обратной величине критического числа Рейнольдса. Так как из мнимой части выражения для числа Рейнольдса потока извлекается корень четвертой степени, имеем для эффективного диаметра 
формулу
.
Величины 
, и при бесконечно малой шероховатости сокращаются, хотя такая ситуация и не реализуется. Но эти формулы справедливы для данной схемы решения. При решении в виде ряда получится вычисляемое другое значение 
. Самая гладкая поверхность соответствует среднему модулю тангенса наклона, равного обратному значению критического числа Рейнольдса, так как самые малые модули тангенса наклона соответствуют молекулярному уровню шероховатости. При этом эффективный диаметр меньше истинного диаметра. Высота шероховатости может стремиться к нулю, одновременно с периодом шероховатости l. Т.е. величина 
это максимальное отношение эффективного диаметра к истинному диаметру. Для внешней задачи эффективный диаметр увеличится, и коэффициент b определится по формуле
.
Для обтекания сферы при нулевой высоте шероховатости эффективный размер сферы совпадет с истинным размером сферы, т.е. коэффициенты 
. При переходе от сферического тела к цилиндрическому с круговым сечением, этот коэффициент будет равен
.
При этом коэффициент b пропорционален
,
который при нулевой макро шероховатости эффективный диаметр равен 
, т.е. при увеличении степени шероховатости эффективный диаметр уменьшается. При этом отношение тангенса наклона макро шероховатостей к микро шероховатостям больше чем величина 
.
Вычислим изменение диаметра трубопровода при условии равенства макро шероховатостей микро-шероховатостям из эмпирической формулы по уменьшению среднего квадрат значения диаметра, получится, что диаметр уменьшился в 
раза. При отношении 
, получаем уменьшение диаметра трубопровода в 
раза.
При этом диаметр изменяется только у коэффициента пульсирующей части решения, т.е. у мнимой части, откуда и возникает множитель 
, так как мнимый член пропорционален 
, который усредняется. Квадратный корень 
при этом соответствует среднему значению диаметра.
При этом учтено влияние шероховатости стенок турбулентного потока на мнимую часть числа Рейнольдса потока. Чтобы получить графики с постоянной высотой шероховатости, надо ввести эффективный средний модуль тангенса угла наклона шероховатостей. Эффективный средний модуль тангенса угла шероховатостей должно зависеть от внешнего давления 
.
Причем, на бесконечности числа Рейнольдса, или безразмерного давления имеем шероховатость, соответствующую постоянной высоте шероховатости
,
где k – среднеквадратичная высота шероховатости, 
– радиус круглого сечения трубопровода.
Формула подобрана таким образом, чтобы правильно определять зависимость числа Рейнольдса относительно внешнего давления и коэффициент сопротивления трубопровода при бесконечном числе Рейнольдса и внешнего давления
при коэффициенте сопротивления, равном величине
.
Эти же коэффициенты получаются при решении с помощью ряда, но вместо величины 
в случае ряда используется коэффициент 
. Другое значение константы приводит к другому значению коэффициента a.
При постоянном среднем модуле тангенса угла наклона шероховатостей 
, но разной высоте шероховатостей k, получается график отличный от графика Никурадзе.
Рис. 1. График зависимости коэффициента сопротивления круглого трубопровода от числа Рейнольдса при разных среднеквадратичных тангенсах наклона шероховатости
Но формула Никурадзе получена при постоянном отношении радиуса трубопровода r0 к средней высоте шероховатости k. Формула (1.1.3) содержит эффективный средний модуль тангенса угла наклона шероховатости, который выражен через отношение радиуса трубопровода к средней высоте шероховатости через безразмерное давление.
Величина 
. Влияния эффективного среднего модуля тангенса наклона шероховатостей на свойство потока зависит от числа Рейнольдса или перепада давления.
Эмпирическая формула по определению коэффициентов 
следующая
, 
,
.
При этом в начале образования мнимой части комплексного решения 
, или начале турбулентного решения, тангенс наклона шероховатостей приближенно равен 2, и графики при разных тангенсах наклона шероховатостей совпадают.
При этом коэффициент сопротивления потока в круглом трубопроводе определяется по формуле 
, число Рейнольдса, вычисленное на основании средней скорости течения потока, равно 
. Коэффициент сопротивления на бесконечности давления пропорционален
.
Приведем график решения, полученного с помощью одного члена ряда.
Рис. 2. Теоретическая и экспериментальная зависимость коэффициента сопротивления круглого трубопровода от числа Рейнольдса при разной степени шероховатости
Для сравнения теоретического и экспериментального графика зависимости коэффициента сопротивления от числа Рейнольдса потока приведен экспериментальный график Никурадзе справа на рис. 2. Ошибка теоретического графика относительно экспериментального примерно 10 %. Отмечу, что график, приведенный на рисунке Никурадзе в ламинарной части решения не точен, по мере приближения к критическому числу Рейнольдса появляется отличие от соотношения 
, что отражено на теоретическом графике, и мало проявляется на экспериментальном. Экспериментальный график для логарифма давления и числа Рейнольдса изображен прямой линией, а должно быть отличие от прямой линии по мере увеличения влияния конвективного члена.
Данный график решения получен для постоянной температуры потока для каждого сечения потока, поэтому в случае слабой зависимости кинематической вязкости от температуры формулы не изменятся. Для его реализации в случае турбулентного режима необходимо подставить в формулу приведенное давление и отношение радиуса трубопровода к высоте шероховатости
Причем формула построена таким образом, чтобы
.
В случае ламинарного режима имеется простая формула для определения числа Рейнольдса 
.
1.2. Алгоритм решения задачи гидродинамики при произвольной геометрии потока
Уравнение Навье – Стокса в декартовых координатах имеет вид
. (1.2.1)
Для этого решаем трехмерную ламинарную стационарную задачу без учета конвективного члена при заданном внешнем воздействии gl
.
Приведем эту задачу к безразмерному виду, разделив на величину 
получим безразмерное уравнение
Решением этой задачи служит функция
.
Строим решение уравнения неразрывности при внешнем воздействии, где величина 
– отклик на внешнее воздействие
. (1.2.2)
Откуда получаем уравнение по определению давления в потоке
Давление ищем в виде
.
Подставляем его в подынтегральное выражение, умножаем на величину 
, и интегрируем по пространству, получаем систему линейных уравнений
,
где имеем значения коэффициентов
где величина 
, определяется внешним воздействием. Приведем уравнение Навье – Стокса к безразмерному виду, разделив его на величину 
, получим безразмерное уравнение
.
Умножаем уравнение Навье – Стокса на сечение трубки тока, записываем уравнение вдоль ламинарного решения, вводя трубку тока, поток в которой постоянен см. [6].
.
В конвективном члене и в градиенте давления введем производную по направлению, соответствующую направлению линий потока ламинарного решения. При подстановке в уравнение решения в виде
, (1.2.3)
где величина Ss – сечение трубки тока в ламинарном режиме; величина 
– стационарное решение уравнения Навье – Стокса без конвективного члена, который для трубки тока равен нулю, как не зависящий от продольной координаты.
Причем, эти трубки тока строим при любом внешнем воздействии, проявляющемся в перепаде давления. Далее учитываем шероховатости, и при некоторых условиях получаем комплексное турбулентное решение, которое связано с влиянием квадратичного конвективного члена с малым множителем, учитывающим шероховатость, который приводит к комплексному решению при большом перепаде давления. При этом действительное решение, которое получается при другом знаке модуля среднего модуля отклонения, отбрасываем, как не определяющее колеблющееся, турбулентное решение. Причем мнимая часть решения определяет пульсации решения. В случае выбора другого знака квадратного корня при учете корреляционной функции процесса 
, где величина 
отклонение скорости от среднего значения, турбулентная вязкость становится отрицательной.
Подставим решение (1.2.3) в уравнение Навье – Стокса и проинтегрируем по трубке тока, разделив на площадь сечения трубопровода. Тогда конвективный член имеет вид
.
Введение шероховатости приводит к зависимости радиуса трубопровода 
от степени макро-шероховатости. Далее выделяем член 
связанный с шероховатостью, находим среднее значение его модуля. При этом произведем усреднение уравнения по величине s. Причем получается, что конвективный член в ламинарном режиме при гладкой поверхности, равен нулю и только при учете шероховатости он имеет значение. Получим уравнение
.
Чтобы учесть шероховатость поверхности трубопровода, и получить турбулентное решение, надо учесть средний модуль тангенса угла наклона шероховатостей. Тогда этот конвективный член будет иметь малый множитель, и конвективный член отличен от нуля и пропорционален средней величине тангенса наклона модуля шероховатостей 
при наличии члена, связанного с не постоянным сечением трубопровода 
. При этом у комплексного турбулентного решения линии тока сохранятся, соответствующие линиям тока ламинарного решения. Но будет пульсирующее решение, относительно ламинарных линий тока. Причем пульсации определяются мнимой часть скорости, причем мнимая часть решения, равная константе, означает пульсации с амплитудой, равной мнимой части скорости.
Подставим решение (1.2.3) в уравнение Навье – Стокса и проинтегрируем по трубкам тока, умножим на 
, где эта величина удовлетворяет условию 
, где 
средний модуль тангенса наклона неустранимых микро-шероховатостей, огибающая которых определяет макро-шероховатости, получим уравнение
где величины 
определены из ламинарного решения и уравнения неразрывности, функция внешнего воздействия 
задается. Получается, что микро-шероховатости, расположенные по всей длине трубопровода, определяют критическое число Рейнольдса. Эти микро-шероховатости меньше макро-шероховатостей, которые оказывают влияние на коэффициент сопротивления, при больших числах Рейнольдса. Но так как число Рейнольдса зависит от геометрии трубопровода, через его диаметр, то критическое число Рейнольдса только обратно пропорционально среднему модулю тангенсу наклона микро-шероховатостей, и зависит от геометрии трубопровода. При этом уменьшение радиуса трубопровода вызывает отрицательное значение 
и, следовательно, отсутствие комплексного, турбулентного решения на размере сужения, т.е. повышается критическое число Рейнольдса. Наоборот расширение трубопровода вызывает увеличение величины 
и, следовательно, уменьшение критического числа Рейнольдса и может вызвать более ранее комплексное решение, т.е. турбулентный режим. Причем, так как число Рейнольдса зависит от температуры через зависимость кинематической вязкости от температуры, получается, что наступление критического числа Рейнольдса находится в зависимости от температуры среды.
Координаты положения равновесия определятся из квадратного уравнения
.
При этом ламинарное решение уточняется и имеет вид 
, которое при малом перепаде давления переходит в линейное ламинарное решение 
.
При этом справедлива турбулентная формула учета шероховатостей в силу одинакового метода усреднения в турбулентном режиме
где величина 
– эффективный средний тангенс наклона шероховатостей; 
– отношение высоты шероховатости к радиусу трубопровода. Критическое число Рейнольдса 
удовлетворяет значению числа Рейнольдса, соответствующее началу комплексного решения. При этом при малом числе Рейнольдса получаем ламинарное решение. Но проблемы получения турбулентного решения на этом не кончаются. Нужно определить влияние шероховатостей поверхности, для чего пока необходимо использование экспериментальных данных. В принципе необходимо точное значение зависимости числа Рейнольдса для гладкой поверхности от макро-шероховатостей. Можно высказать предположение, что при счете с помощью ламинарного решения при подсчете турбулентного решения влияние шероховатости будет аналогично влиянию в случае круглого трубопровода с другими константами 
, которые нужно определить из значения коэффициента сопротивления на бесконечности числа Рейнольдса. Для этого необходимо, чтобы размер, входящий в число Рейнольдса определялся аналогично в круглом трубопроводе и для тела произвольной формы. Для этого необходимо определять средний радиус поверхности через площадь поверхности тела и вводить радиус сферы, у которой площадь совпадает с площадью данного тела. Но у внешней задачи имеются свои особенности, связанные с наличием кризиса сопротивления, который связан с наличием следа за телом в потоке. Этого следа нет во внутренних задачах, таких как течение в трубопроводе.
Построим решение уравнения Навье – Стокса для внешней задачи. Имеется ламинарное решение для движения сферы в жидкости при малом числе Рейнольдса. Оно определяет следующее распределение скорости см. [8]
Отметим, что это ламинарное решение справедливо с точностью 
, а не с точностью R, что связано с преобразованным видом уравнения Навье – Стокса
.
Которое имеет стационарное решение 
и нелинейным членом можно пренебречь при условии 
. В самом деле, ламинарное решение этого уравнения 
и нелинейный член равен 
, что эквивалентно 
. Последнее равенство выполняется в ламинарном режиме. Условие наступления комплексного, турбулентного решения является равенство нулю дискриминанта этого уравнения при числе Рейнольдса равном критическому.
При этом определится зависимость давления от параметров потока
.
Определим скорость жидкости в декартовой системе координат
где 
.
Направление линий тока определяется по формуле
Значение единичного вектора вдоль линии тока равно
.
Вдоль линии тока в декартовой системе координат скорость потока равна константе. При этом равняется нулю первая производная вдоль линии тока и не равна нулю вторая производная вдоль этого же направления.
Найдется значение скорости вдоль линии тока из условия, что градиент скорости ортогонален линии тока. В самом деле, вдоль линии тока изменение скорости равно нулю, и максимально изменяется в плоскости, перпендикулярной линии тока. При этом в случае гладкой поверхности нелинейный член равен нулю. Поэтому выполняется
,
где Us удовлетворяет (1.2.3). При этом Us первый интеграл дифференциального уравнения
не является первым интегралом дифференциального уравнения 
, при этом выполняются условия 
. При вычислении оператора Лапласа вдоль линии тока 
используется направление и дифференциальное уравнение
.
Задание параметров 
на поверхности сферы определит кривую линии тока.
А при вычислении первого интеграла используется дифференциальное уравнение
.
При этом скорость равна проекции ламинарных решений на линию тока
, (1.2.3)
причем, направление 
– это направление линий тока.
Решение уравнения Навье – Стокса ищем в полученном из предыдущих рассуждений виде
Эта формула определяет постоянное значение скорости вдоль траектории 
, заданной 
, равенство нулю нелинейного члена и в общем случае не равенство нулю второй производной от скорости вдоль траектории. Причем, дифференцирование осуществляется по величине s, и полагаем, вдоль траектории движения выполняется равенство 
. Подставляем найденное значение скорости Us в уравнение Навье – Стокса, записанное вдоль линии тока получим
При этом вдоль направления линии тока 
выполняется
,
при этом величина 
, так как вторая производная вдоль направления линии тока не равна нулю.
Интегрируя по величине s, получим нелинейное уравнение
.
Это уравнение имеет стационарное решение
.
Также существует комплексное решение при большом перепаде давления
,
где величина g определяется по формуле
.
Получается такая степень s, так как
.
Если для внутренней задачи шероховатость уменьшает эффективное сечение потока, то для внешней задачи должно увеличивать эффективное сечение сопротивления потока.
Величина 
. Величина k – средний модуль высоты шероховатостей. Влияния эффективного среднего модуля тангенса наклона шероховатостей на свойство потока зависит от числа Рейнольдса или перепада давления.
Эмпирическая формула по определению коэффициентов 
следующая
,
,
.
Решение с помощью ламинарного члена определяет основное состояние турбулентного потока и является однозначным. При этом комплексное решение является бесконечное число раз дифференцируемой функцией по координате и по внешнему безразмерному давлению, кроме одного значения давления, соответствующего критическому режиму. Это значение давления определяет критическое число Рейнольдса, перехода в турбулентный комплексный режим. Но решение в виде ряда определяет при некоторых условиях счетное количество решений задачи Коши, имея счетное количество точек ветвления.
1.3. Построение решения в виде ряда для трубопровода с круговым сечением
Можно привести уравнение Навье – Стокса к безразмерному виду для ряда, являющегося решением
Запишем уравнение Навье – Стокса в безразмерном виде, умножая на величину 
, и, умножая уравнение неразрывности на величину 
, получим уравнения
Имеем значение переменных
, 
,
,
где величина a – радиус трубопровода, 
– плотность и кинематическая вязкость жидкости, 
– скорость и давление в потоке. Величина давления равна
,
где 
– начальное и конечное давление в потоке. Имеется связь
,
которая является определением давления 
. Величина скорости определяется по формуле
.
Величина a(z) определяет высоту шероховатости, при постоянном радиусе трубопровода. Уравнение неразрывности удовлетворяется тождественно с точностью
,
так как радиальная и угловая скорость равна нулю в силу отсутствия внешнего радиального и углового воздействия и нулевых начальных условиях. Причем физический смысл
.
Критическое число Рейнольдса имеет иной вид по отношению к решению раздела 1.1, так как используется другая опорная функция. С помощью метода Галеркина приводим задачу к системе обыкновенных дифференциальных уравнений.
Удалось решить уравнение Навье – Стокса с точность до 0.6 % при постоянном среднеквадратичном значении тангенса наклона шероховатостей. Строилось решение той же задачи на основе алгоритма (1.3.1). Получился следующий график (рис. 3).
Рис. 3
При этом формулы, по которым считался данный алгоритм с турбулентным решением в виде ряда
(1.3.1)
где величина 
определяется по формуле (1.3.2). При этом, так как учитывается множество членов, отражающих разную степень завихренности, корень из мнимой части комплексного числа Рейнольдса потока не извлекается. Степень шероховатости учитывалась из условия одинаковости коэффициента сопротивления при разных формулах для решения, изменялся только коэффициент a, так как на бесконечности числа Рейнольдса величина X является константой, определяемой формулой численного метода и необходимо, чтобы
откуда определится коэффициент 
при условии средней гладкости 
. Здесьвеличина 
является корнем уравнения
Решение этой системы нелинейных уравнений ищется в виде 
, причем для коэффициента 
получаем квадратное уравнение
(1.3.2)
При этом величина 
определяется по формуле 
.
Решение этой алгебраической системы уравнений имеет N+1 корень. Докажем это. Допустим, нашлось одно решение этой системы уравнений 
. Алгоритм его построения описан выше по тексту. Составим решение 
. Подставим это решение в уравнение по определению координат положения равновесия, получим уравнение
.
Получим уравнение
Выберем начальное приближение так, чтобы определитель этой линейной системы уравнений равнялся нулю. Для этого определим 
, и найдем из равенства нулю определителя значение a. Тогда следующее приближение определит значения 
с точностью до множителя. Вычислим этот множитель из условия 
. Это можно сделать N способами, определяя корень полинома N степени. В результате получим n+1 решение алгебраического уравнения.
Получение счетного количества (
) координат положения равновесия эквивалентно получению счетного количества решений для уравнения Шредингера. Эквивалентность уравнения Шредингера и Навье – Стокса докажем в данной статье далее по тексту.
При этом основное решение 
определится однозначным образом. Это доказывается алгоритмом его построения и невозможности выделить из решений 
другую константу, удовлетворяющую уравнению по определению координат положения равновесия. Отметим, так как коэффициенты 
меньше, чем линейный член при малом перепаде давления, значение 
соответствует линейной части решения. Продолжая это решение при большом перепаде давления, получим основной член решения задачи гидродинамики.
График на рис.3 отличается от графика Никурадзе на рис.4. Асимптотика при больших и малых значениях числа Рейнольдса совпадает с графиком Никурадзе. Но график построен при постоянном тангенсе наклона шероховатости, а график Никурадзе при постоянном отношении радиуса трубопровода к средней высоте шероховатости. Построим значение эффективного тангенса наклона шероховатостей, по средней высоте шероховатости.
причем, эквивалентная степень шероховатости зависит от безразмерного давления. Эквивалентную степень шероховатости можно получить при расчете динамики процесса течения жидкости на микро шероховатостях постоянной высоты. Величина критического давления определяется из формулы 
. Эмпирическая формула по определению коэффициентов 
следующая
, 
, 
,
где 
– коэффициенты a, образующие множитель при мнимой части решения для решения в виде одного члена и в виде ряда соответственно.
При этом в начале образования мнимой части комплексного решения, или начале турбулентного решения, тангенс наклона шероховатостей приближенно равен 2, и графики при разных тангенсах наклона шероховатостей совпадают.
Рис. 4
График построен для отношения радиуса трубопровода r0 к высоте шероховатостей 
равной 
. Точность аппроксимации относительно графиков Никурадзе примерно 10 %.
2. Построение решения обыкновенных дифференциальных уравнений
Рассмотрим систему нелинейных автономных дифференциальных уравнений
. (2.1)
К системе нелинейных дифференциальных уравнений сводится система уравнений Навье – Стокса и уравнение неразрывности.
, (2.2)
где величина трехмерной скорости определяется по формуле
.
Причем функция 
определяется в виде синуса. Тогда коэффициенты 
для непрерывной функции убывают с ростом индекса как величины не большие 
и возможна редукция, т.е. вместо бесконечного числа членов используется конечное число членов.
Исследуются комплексные и действительные решения задачи Коши этого дифференциального уравнения в случае действительных и комплексных начальных условиях, при действительном аргументе t. Начальные условия имеют вид 
, где величина t0 соответствует начальному моменту интегрирования, а величина 
в общем случае комплексная. Причем в случае действительных значениях 
, правая часть (2.1) действительна.
Систему дифференциальных уравнений (2.2) можно представить при не кратных положениях равновесия путем подстановки 
. При этом положения равновесия системы (2.2) 
перейдут в положения равновесия 
. При этом определятся собственные числа и собственные векторы линеаризованной системы (2.2).
Система уравнений (2.2) запишется в виде
. (2.3а)
Систему уравнений (2.3а) можно записать в виде
, (2.3б)
где введен не обращающийся в ноль множитель 
, который равен
.
При подстановке этого множителя в (2.3б), получим (2.3а). Покажем, что этот множитель в ноль не обращается. Величины 
удовлетворяют условию 
, где величина S конечна.
При условии 
имеем конечный предел
где произвели сокращение множителя 
, числитель дроби в ноль не обращается, так как рассматриваются не совпадающие корни, являющиеся координатами положения равновесия. Показали, что этот множитель в ноль не обращается при не совпадающих координатах положения равновесия. Величины 
удовлетворяют условию
,
где величина S конечна.
При этом дифференциальное уравнение можно записать в виде
, (2.4)
где 
– стремящаяся к бесконечности функция при условии стремления координат к положению равновесия. В случае решения в действительной плоскости это монотонная функция.
Лемма 2. Необходимым и достаточным условием стремления неизвестной функции к устойчивым координатам положения равновесия является условие 
, причем 
.
Имеем соотношения
(2.5)
при 
и значит 
, как интеграл от константы. Справедлива и обратная теорема, при условии 
, реализуется одно из устойчивых положений равновесия. Это следует из вида решения, при условии 
имеется отрицательная действительная часть у числа 
в формуле (2.6) согласно лемме 4, и решение стремится к координате положения равновесия 
в формуле (2.4). При этом величина времени стремится к бесконечности.
Лемма 3. Решением дифференциального уравнения (2.1) является функция 
, удовлетворяющая формуле (2.6).
Для получения (2.6) разделим уравнение (2.4) на произведение множителей 
и умножим (2.4) на величину 
. Раскладываем полученную дробь на сумму простых дробей и их интегрируем. Получаем уравнения
.
Здесь в случае излучения звуковой энергии в интервале 
получим разные ветви логарифма.
Потенцируя полученное выражение, получим (2.6)
(2.6)
где все значения координат положения равновесия не кратные и не зависят от процесса излучения, происходящего в интервале 
. Причем в случае ламинарного действительного решения излучения не будет, а в случае турбулентного, сопровождающегося излучением решения будет переход энергии. В самом деле, наличие излучения определяет комплексное решение, которое описывает турбулентный пульсирующий режим. При этом турбулентный режим при перестройке решения сопровождается звуковым шумом. На координаты положения равновесия, определяющие стационарное решение, экспоненциальный множитель не оказывает влияние.
Наличие S совокупностей не кратных координат положения равновесия налагает следующие ограничения на систему нелинейных уравнений. Правая часть системы дифференциальных уравнений может быть представлена в виде
,
где если величины 
не равны нулю, то они имеют множитель 
, иначе уравнение не имеет координат положения равновесия. Причем функция 
, имеющая S совокупностей координат положения равновесия должна обязательно иметь множители с одним или несколькими n, входящих в произведение 
.
Тогда нет определенной ветви вычисления координаты положения равновесия и под координатой 
можно понимать любую из координат 
. Т.е. совокупности координат перепутаны и могут состоять из любой комбинации и тогда обязательно имеются устойчивые положения равновесия.
Когда система нелинейных уравнений по определению координат положения равновесия, сводится к одному уравнению, относительно одной неизвестной, а остальные неизвестные координаты положения равновесия, являются функцией этой определяемой неизвестной, координаты положения равновесия не перепутаны, а являются отдельными совокупностями, тогда возможны случаи, когда все координаты положения равновесия неустойчивы.
При этом необходимо уточнить понятие кратного корня. Если в одном уравнении корни одинарны, то они не являются кратными, хотя и перепутаны как координаты положения равновесия. При этом нельзя добиться равенства нулю определителя линеаризованной эквивалентной системы дифференциальных уравнений см. текст перед теоремой 2, и значит решение не хаотическое. Если в одном уравнении имеются одинаковые корни, то они являются кратными и можно добиться равенства нулю определителя линеаризованной эквивалентной системы дифференциальных уравнений, см. текст перед теоремой 2, и решения хаотические.
Лемма 4. Сумма коэффициентов 
по индексу s равна нулю, т.е. 
.
Для доказательства этого тождества рассмотрим полином S–1 степени относительно y
,
В точках положения равновесия 
полином удовлетворяет 
. В силу единственности полинома степени S–1, проходящего через S точек, получаем 
, так как это значение удовлетворяет точкам аппроксимации. Распишем формулу для полинома, равного единице, разделив его на произведение 
, получим
полагая, 
получим тождество 
, в случае, если имеется S+1 положение равновесия.
В случае если разлагается дробь
,
где 
полином степени S–1. Значение коэффициентов изменится, но свойство 
останется, 
. Докажем это. Для чего рассмотрим сумму
.
Эта сумма равна 
. Распишем формулу для полинома, равного 
, разделив его на произведение 
, получим
Полагая 
, получим тождество 
, в случае, если имеется S+1 положение равновесия.
Но чтобы реализовать решение, надо знать положения равновесия этой системы нелинейных уравнений. Кроме того, положения равновесия могут быть кратные, что изменяет ход решения, оно становится случайным, или хаотическим. Но удается доказать важную теорему.
Теорема 1. Рассматривается задача Коши при произвольных действительных начальных условиях для системы нормальных нелинейных обыкновенных дифференциальных уравнений (2.1). Случай вырожденного решения задачи Коши – положения равновесия, не рассматривается. В случае, если у системы (2.1) имеются комплексно-сопряженные положения равновесия с действительной частью, то при конечном аргументе t действительное решение задачи Коши системы (2.1) при действительных начальных условиях стремится к бесконечности. Потом это решение переходит в комплексное решение, стремясь к положению равновесия, в случае, если комплексные координаты положения равновесия имеют действительную часть. При этом правую часть (2.1) считаем регулярной функцией, действительной при действительных аргументах. Она имеет конечное число не кратных положений равновесия.
Доказательство.
Если решать систему (2.2) при не кратных положениях равновесия, то получим согласно с леммой 3
(2.7)
где 
выделенное комплексное положение равновесия, 
остальные положения равновесия. Коэффициенты 
удовлетворяют 
в соответствии с леммой 4. При этом в сумме 
величина действительной части 
в случае комплексного значения 
участвует дважды и в силу того, что все числа 
удовлетворяют условию 
, имеем формулу 
.
Обоснуем формулу (2.3). Для этого два комплексно сопряженных члена решения преобразуем (для упрощения записи индекс l опускаем)
, (2.8)
где 
. После интегрирования (2.8) по аргументу x, получим формулу (2.7)
.
Решение равняется
,
где
При этом величина 
действительная в силу существования комплексно-сопряженных положений равновесия. Т.е. имеем равенство при условии 
при конечном t
. (2.9)
Это уравнение имеет решение, стремящееся к бесконечности.
При этом решение дифференциального уравнения при росте 
согласно лемме 2, может иметь комплексные корни
.
При этом, так как справедливо 
согласно лемме 4, и положения равновесия имеют действительную часть, имеются числа с отрицательной действительной частью 
, значит, имеется сходимость к одному из положений равновесия. Действительное решение будет стремиться к бесконечности, причем нарушатся условия существования и единственности задачи Коши. При этом при бесконечности 
согласно лемме 2 неизвестная функция будет стремиться к одному из положений равновесия. Это положение равновесия не может быть действительным, так как действительное решение бесконечно. Значит, решение будет иметь точку ветвления и стремиться к комплексному положению равновесия. Значит, при комплексных положениях равновесия получается конечное комплексное решение при изменении 
. Т.е. в некоторой точке начнется комплексное решение.
Конец доказательства.
Приведем пример, описывающий это свойство дифференциального уравнения, переход к комплексному решению. Так, для дифференциального уравнения может возникнуть комплексное решение, вместо бесконечного действительного решения
.
Причем положения равновесия чисто мнимые 
, и значит, решение может не стремиться к положению равновесия. Причем действительное решение этого дифференциального уравнения быстро стремиться к бесконечности 
.
Используя неявную схему решения, получим следующее уравнение
.
Разрешая относительно неизвестной функции x, получим неявную схему
.
Эта неявная схема с постоянным шагом правильно описывает стремление решения к бесконечности. При счете с уменьшенным шагом она определяет большее значение переменной t и значит, определяет большее значение неизвестной функции. Т.е. правильно описывает решение дифференциального уравнения до бесконечности решения. Когда бесконечность достигнута, при условии 
определится конечное комплексное решение. Численный счет этого уравнения подтвердил правильность проведенного анализа решения.
Причем комплексное решение обладает новыми свойствами, оно сложным образом вращается вокруг положения равновесия. При этом действительное решение стремится к бесконечности, т.е. правая часть дифференциального уравнения стремится к бесконечности, и нарушаются условия существования и единственности решения задачи Коши, и возникает дополнительное комплексное решение.
Решение с комплексными начальными данными определится формулой
при любом t. Т.е. приближенно имеем
При этом знаменатель этой дроби в ноль не обращается.
Т.е. конечного решения задачи в действительной плоскости не существует. А в комплексной плоскости имеется конечное непрерывное решение в случае не кратных положениях равновесия.
Но возникает вопрос, каков физический смысл мнимой части комплексного решения.
3. Физический смысл комплексного решения
Опишем физический смысл комплексного турбулентного решения. Итак, рассмотрим действительное решение системы обыкновенных дифференциальных уравнений 
. Пусть начальные данные имеют среднее 
и дисперсию 
(дисперсия начальных данных в случае уравнения Навье – Стокса определяется шероховатостью поверхности или не точно заданными начальными данными). Тогда для дисперсии решения имеем
.
Значит имеем
. (3.1)
Приведу формулировку обратной теоремы Пифагора. Для всякой тройки положительных чисел a,b и c, такой, что 
, существует прямоугольный треугольник с катетами a и b и гипотенузой c. Значит, математическое ожидание и среднеквадратическое отклонение образуют катеты, а гипотенузой является корень из среднего квадрата величины. Т.е величина среднего 
ортогональна среднеквадратическому отклонению 
, которое образует мнимую часть координаты тела. Таким образом, полученное в результате усреднения во времени декартово пространство с колебательной скоростью высокой частоты (период колебания меньше времени измерения) становится комплексным пространством. Т.е. в случае большой дисперсии величины действительного пространства, его нужно рассматривать как комплексное трехмерное пространство, где мнимая часть соответствует среднеквадратическому отклонению. При этом имеется следующая связь между переменными
,
причем комплексное число a выбирается из условия, чтобы мнимая часть имела положительное или отрицательное значение. Этому удовлетворяет среднеквадратичное отклонение. Но иногда среднеквадратичное отклонение положительно, например, в случае диэлектрической проницаемости, где вмешиваются положительные и отрицательные заряды. Тогда имеем формулу
,
где действительная часть пропорциональна положительному среднеквадратичному отклонению диполя, а проводимость пропорциональна среднему значению. Но зато проводимость делится на частоту, которая имеет положительный и отрицательный знак.
Следовательно, алгоритм нахождения среднего решения, или среднего в фазовом пространстве решения, и его дисперсию сводится к нахождению комплексного решения. Среднее решение соответствует действительной части решения, а квадрат комплексной части соответствует дисперсии решения. Таков физический смысл комплексного решения, действительная часть – это среднее решение, а мнимая часть – это его среднеквадратическое отклонение. Причем действительная и мнимая часть ортогональны, и образуют комплексное пространство. В самом деле, согласно обратной теореме Пифагора в силу формулы (3.1) математическое ожидание и среднеквадратическое отклонение образуют катеты, а средний квадрат является гипотенузой.
Отметим, что при расчете движения потока при учете одного члена ряда, надо брать квадратный корень из мнимой части, так как вычисляется поступательная скорость. Мнимая часть соответствует квадратному корню из колебательной части безразмерной скорости. Ситуация аналогична вычислению отклонения при случайном выборе шага вперед или назад с вероятностью ½ и положение точки после N шагов определяется значением 
. Действительная часть решения и мнимая часть решения расположены на разных осях комплексного пространства. Но если усреднять мнимую безразмерную часть, получим
причем, решение равно модулю последней величины, и при разных шероховатостях мнимую часть решения надо умножать на усредняющий множитель. При этом, если вычисляются все коэффициенты ряда из системы нелинейных уравнений, извлекать квадратный корень из мнимой части не надо. Необходимо складывать комплексные величины и брать модуль суммы.
Покажем, что мнимая часть комплексной производной от координаты в фазовом пространстве дифференциального уравнения, которую назовем скоростью, образует пульсирующее перемещение координаты в фазовом пространстве, т.е. в пространстве переменных 
.
Используется в качестве переменных среднее значение, ведь на молекулярном уровне среда не гладкая.
Лемма 5. Комплексное решение определяет колеблющуюся пульсирующую функцию координат перемещения потока.
Мнимая часть скорости соответствует скорости вращения в фазовом пространстве. Так как известен радиус вращения, то можно определить и частоту вращения. В плоскости вращения комплексную скорость с постоянным радиусом вращения и постоянной частотой можно представить в виде 
. В случае переменной по пространству стационарной скорости эту формулу можно представить локально в одной плоскости в виде
причем частота зависит от времени, так как смещение фазы обеспечивается гармоническими колебаниями в соседних точках. Сумма гармонических колебаний с разными частотами, зависящими от времени, определяет пульсирующий режим в фазовом пространстве, при стационарной комплексной скорости. Т.е. получается, что комплексная скорость описывает пульсирующие во времени координаты точек фазового пространства. Ситуация аналогична наличию нескольких стационарных вихрей, описывающих пульсирующее вращение потока.
Лемма 6. Трехмерную скорость потока можно представить в виде
Причем, скорости определяются в виде интеграла от касательного ускорения, по формуле
Интеграл от нормального ускорения определяет нормальную компоненту скорости, по формуле
.
При этом величина локальной скорости 
. Но проинтегрированная относительно центростремительного ускорения скорость отлична от нуля 
, обращаясь при постоянной скорости частицы и постоянном радиусе кривизны, за период
,
где величина R – радиус кривизны, в ноль при той же помеченной начальной точке. При переменной скорости частицы за время, когда один из интегралов 
, которое, при конечном радиусе кривизны одного знака траектории, конечно и равно
,
так как касательное направление 
, при вращении меняет знак.
При этом вклад в поступательную часть комплексной скорости за один оборот вращения равен
.
Причем, если эту связь записать в безразмерном виде, получим
.
Причем имеем значение нормальной компоненты скорости 
, где величина 
соответствует замкнутой траектории скорости вращения, а величина 
вклада в поступательную скорость. Откуда имеем величину дополнительного вклада в поступательную скорость за счет мнимой части числа Рейнольдса 
. Откуда имеем дополнительный вклад в поступательную скорость 
.
При этом в случае точного решения в виде ряда складываются члены ряда, и результирующая скорость усредняется и она равна корню из модуля скорости. Если же используется один член ряда, то он не усредняется и из мнимой части надо извлечь корень, и использовать модуль скорости с измененной мнимой частью.
При этом тело сместится относительно помеченной начальной точки. Чтобы смещение было существенным радиус кривизны должен менять свой знак. Причем, когда этот период мал, по сравнению с временем процесса, это вращение воспринимается как мнимое среднеквадратичное отклонение скорости. Отметим, что тангенциальное ускорение и нормальное ускорение образуют скорость, которая направлена по касательной к траектории частицы. Величины 
это тангенциальные и нормальные орты. Тангенциальное ускорение определяется по формуле
.
Направление скоростей 
ортогонально и их сумма приводит к приращению модуля скорости движения
,
так как
.
Дифференцируемые по времени компоненты этих проекций определяют тангенциальное и нормальное ускорение. При этом вводится понятие тангенциальной и нормальной скорости, которые в декартовом пространстве не ортогональны 
, но в шестимерном комплексном пространстве ортогональны и их модуль комплексного вектора 
равен
.
Это доказывается представлением 
и вычислением модуля как произведения комплексно сопряженных векторов с учетом ортогональности шести действительных ортов.
Итак, получено решение уравнений Навье – Стокса в случае не кратных положений равновесия. Оно имеет вид
.
,
где величины 
координаты положения равновесия.
Ламинарное решение соответствует координатам положения равновесия, структура турбулентного решения имеет вид
,
где 
известная определенная непрерывная функция, величина tn определится из начальных условий, причем 
. Причем решение содержит множество полюсов, которые для действительного решения при действительных начальных данных приводят к бесконечности. Комплексное решение при действительном времени и комплексных начальных условиях, определяющих комплексное значение tn, конечно. При этом комплексное решение стремится к комплексным координатам положения равновесия.
При этом формула
. (3.2)
Может содержать точки ветвления, в которых решение непрерывным образом переходит в другую ветвь решения. Это не противоречит теореме единственности решения задачи Коши, так как левая часть дифференциального уравнения стремится в точке ветвления к бесконечности. Производная от правой части обыкновенного дифференциального уравнения также стремится в точке ветвления к бесконечности. Получается точка не продолжаемого решения. Но это решение можно продолжить с помощью формулы (3.2).
4. Построение решения в случае кратных положений равновесия
Но как же строить решение в случае кратных координат положения равновесия. При этом одно из собственных чисел равно нулю и в формуле (2.5) собственное число равно нулю и значит, hl может не стремиться к бесконечности, что нарушает теорему о сходимости к комплексному положению равновесия. Для построения решения докажем теоремы 2, 3.
Решение задачи Коши для обыкновенных дифференциальных уравнений описывает детерминированное решение. Показано, что в случае кратных значений положения равновесия может иметься хаотическое решение. При этом малым изменениям начальных данных соответствует большое изменение решения.
Ставится задача нахождение условий, когда система нелинейных уравнений имеет хаотическое решение. Т.е. перескок от одного положения равновесия к другому. Эта система уравнений описывает скачкообразное, турбулентное решение уравнений гидродинамики, когда локальная скорость потока не предсказуема во времени. Оказывается, что это условие наступает, при кратных положениях равновесия системы обыкновенных нелинейных уравнений, к которым сводится задача гидродинамики.
Рассмотрим систему нелинейных автономных дифференциальных уравнений, имеющих двукратный корень
. (4.1)
Осуществим вырожденное преобразование системы дифференциальных уравнений с кратным корнем (кратный корень является кратной координатой положения равновесия) с помощью неизвестных функций 
. Подставим в дифференциальное уравнение, получим
. (4.2)
Будем рассматривать для l уравнения координату 
независимой, а остальные координаты как функции времени, являющиеся решением дифференциального уравнения. Система уравнений разбивается на N независимых неавтономных уравнений, которые в случае, когда решение равно координатам положения равновесия, правая часть дифференциального уравнения не зависит от времени. При этом связь между решениями относительно преобразования 
одинакова со связью между положениями равновесия (связь вырожденна в кратном положении равновесия, но при этом для остальных произвольных значениях xl получается единственное значение yl). При этом справедливо следующее равенство
.
Т.е. допустим, матрица уравнения (4.1) 
не вырождена, но система (4.1) имеет двукратное положение равновесия. Тогда матрица уравнения (4.2), полученного из уравнения (4.1), с помощью преобразования 
равна (допустим кратный корень имеет последнее уравнение системы (4.1), тогда при дифференцировании по 
в точке двукратного положения равновесия останется только последнее уравнение, остальные уравнения имеют неопределенность, которая равна константе и при дифференцировании дают ноль)
.
При этом первый член второй строки этой формулы является константой в точке положения равновесия, так как
и выбором преобразования 
можно добиться, что матрица 
может быть вырождена в точке кратных координат положения равновесия системы (3.1).
Но преобразование переменных 
не произвольно. Преобразование таково, что действительным или комплексным корням 
должны соответствовать действительные или комплексные корни 
. Строить систему (4.2) необходимо, так как уравнение (4.2) имеет не кратные положения равновесия, соответствующие положениям равновесия системы (4.1). По координатам системы (4.2) однозначно определяются соответствующие положения равновесия системы (4.1). При этом добиваемся, чтобы определитель линеаризованной системы уравнений (4.2) равнялся нулю.
При условии, что содержится кратный корень 
, являющийся одинаковым решением s и p ветви решения алгебраического уравнения, эквивалентное N уравнение для переменной xN имеет вид
.
При этом формула учитывает несколько ветвей координат положения равновесия. Значит, частная производная в положении равновесия имеет вид
.
При этом эквивалентный определитель системы дифференциальных уравнений равен нулю. Но формула для линеаризованной системы уравнений учитывает одну ветвь координат положения равновесия, упуская нулевое значение для определителя 
в случае кратных положений равновесия при приближении xN к значениям 
или 
. При линеаризации системы нелинейных уравнений учитывается либо 
, либо 
, а надо учитывать обе одновременно.
При двукратном положении равновесия определитель линеаризованной системы не обязательно равен нулю, например, аттрактор Лоренца, имеет двукратные положения равновесия, но его определитель линеаризованной системы не равен нулю см. [1]. Приведем пример существования хаотического решения при кратных корнях. Это аттрактор Лоренца
Эта система имеет положения равновесия
При этом корень 
двукратный. Причем, определитель линеаризованной системы этих уравнений не вырожденный. Эта задача в хаотическом режиме вращается вокруг одного из положений равновесия. Потом решение системы нелинейных уравнений резко перескакивает к другому положению равновесия и вращается вокруг второго положения равновесия. Это связано с тем, что хотя положения равновесия кратные, определитель линеаризованной системы уравнений у аттрактора Лоренца не равен нулю. Но у системы линеаризованных уравнений, построенных с помощью двукратного корня по мере приближения к точке, соответствующей 
, в случае 
согласно теореме 2 образует скачок от этого положения равновесия. При этом происходит произвольный перескок на решение с положением равновесия
.
Причем возможно колебание между кратными положениями равновесия, так как не кратное положение равновесия неустойчиво.
В книге [6], описан сценарий рождения аттрактора Лоренца через неполный двойной гомоклинический каскад бифуркаций. В этой книге считается, что условием хаотического решения является наличие в системе Лоренца седло – узла и двух седло – фокусов, откуда вытекает возможность существования в ней различных гомоклинических и гетероклинических контуров особых точек и связанных с ними каскадов бифуркаций. Терминологию и объяснение обозначений см. [6].
В предлагаемой статье определен простой критерий существования хаотического решения. Наличие хаотических решений связано с наличием кратного положения равновесия системы уравнений.
Теорема 2. Пусть 
координаты одного из положений равновесия системы (2.1). В случае двукратного значения положения равновесия 
, линейное приближение решения дифференциального уравнения с устойчивым положением равновесия 
, определяет сходимость к положению равновесия при условии 
, где матрица 
, это собственные векторы матрицы линеаризованной системы. В случае противоположного условия 
положение равновесия не достижимо и получается многозначное решение. Предполагается, что выполняется 
.
Доказательство.
В силу существования двукратного корня, равного 
уравнение (2.1) можно привести к виду (4.3) с вырожденной матрицей линеаризованной системы
(4.3)
Т.е. определитель матрицы 
равен нулю в этом представлении.
Решаем линеаризованное уравнение (4.4)
(4.4)
Это уравнение имеет решение в окрестности положения равновесия
,
так как 
, где собственные векторы 
и собственные числа 
определятся из системы уравнений
Причем наблюдается приближение 
к значению 
. Т.е. положение равновесия не обязательно достижимо в случае кратных корней, даже для системы с 
, так как выполняется условие 
. Найдем условие сходимости к кратному корню. Имеем соотношение
.
При этом для определения константы cN при дискретном вычислении решения имеем следующее рекуррентное соотношение 
при достижении устойчивыми координатами положения равновесия. Т.е условием сходимости решения в направлении 
является неравенство 
. При условии 
наблюдается отсутствие сходимости в одном направлении 
, остальные cl не растут. Скачок решения осуществляется мгновенно на произвольное значение величины cN, по формуле
,
в момент времени, когда по устойчивым направлениям достигнуто положение равновесия. При этом решение увеличится в разы
. (4.5)
Но в силу приближенности решения, получится не бесконечность, а переход к другому положению равновесия. Точное решение при численном счете получается мгновенно при условии 
и имеет значение, зависящее от малости величины h. Т.е. происходит произвольный скачок. Остальные возможно комплексные направления 
собственного вектора, соответствуют устойчивому собственному числу положения равновесия.
Причем в случае численной схемы получится в случае рационального значения k в формуле (4.5) конечное число состояний, а в общем случае при иррациональном значении k получится счетное число состояний. Численный счет в случае k целого определит единственное решение, так как период мнимой части фазы умножается на целое число.
Если же наблюдается одно положительное собственное число при остальных отрицательных, то имеем
и рано или поздно растущий член будет иметь большее значение, и будет удаление от положения равновесия.
Получается, что в случае одного нулевого собственного числа и остальных отрицательных собственных числах, наблюдается приближение к положению равновесия вдоль комплексных собственных векторов с отрицательными действительными частями собственных чисел. При определенных условиях имеется удаление решения от положений равновесия по комплексному направлению собственных векторов с нулевым собственным числом.
Конец доказательства.
Опишем структуру возможно хаотического решения, которое является турбулентным при кратных положениях равновесия. В случае решения уравнения Навье – Стокса в виде ряда кратные корни, это свойство комплексного решения, или турбулентного решения.
Теорема 3. В случае системы (4.1) с двукратным корнем положения равновесия решение задачи Коши в комплексной плоскости для системы дифференциальных уравнений (4.1) с действительными и комплексными положениями равновесия будет определяться рядом Лорана с существенно особой точкой и, следовательно, при приближении к особенности может носить хаотический характер. Решение может иметь полюсы со знаменателем
,
где K+2 количество положений равновесия. Положений равновесия должно быть больше двух. В общем случае решение задачи не единственно, а имеется счетное количество решений. Причем реализуется состояние с наименьшей энергией. Аргумент решения t является действительным.
Доказательство.
В случае двукратного корня aN систему дифференциальных уравнений (2.1) можно представить в виде
, (4.6)
где введен не обращающийся в ноль множитель 
, который равен
.
При подстановке этого множителя в (4.6), получим (4.2). Покажем, что этот множитель в ноль не обращается. Величины 
удовлетворяют условию
,
где величина K конечна.
При условии 
имеем конечный предел
При условии 
тоже имеем конечный предел
.
В случае совпадающих корней сокращается множитель 
. При этом получается не обращающийся в ноль множитель 
.
При этом дифференциальное уравнение можно записать в виде
, (4.7)
где 
– стремящаяся к бесконечности функция при условии приближения к положениям равновесия. В случае решения в действительной плоскости это монотонная функция.
Запишем дифференциальное уравнение с кратным значением положения равновесия
Величина 
в ноль не обращается. Решение этого дифференциального уравнения имеет вид
(4.8)
– координаты положения равновесия, величина 
– значение двух кратного положения равновесия. Вычислим коэффициенты этого разложения
Имеем три возможных варианта разложения функции в ряд. Либо в точке 
имеется существенная особая точка, и решение стремится к произвольным точкам при приближении к особенности, либо справедливо
(4.9)
либо решение регулярно. При этом величина 
конечна. Допустим справедливо (4.9). Определим решение этого дифференциального уравнения, для чего подставим часть решения в виде полюса 
в уравнение
. (4.10)
При этом оказывается, что приравниваются члены
,
что приводит к равенству 
и для множителя справедливо 
. При этом решение при условии 
имеет сложную структуру и не имеет однозначного решения в комплексной плоскости, т.е. в точке 
имеется ветвление решения и значит, возможен переход к комплексному решению. Имеем 
. Т.е. решение в этом случае является турбулентным, комплексным и имеющим счетное количество ветвей.
Допустим, решение регулярно. Разложим решение в окрестности точки 
по положительным степеням 
, т.е. 
. Тогда решение представимо в виде
.
Т.е. справедлива формула
.
Т.е. предположение о регулярности решения не подтверждается, решение имеет бесконечное число членов с отрицательной степенью 
. Значит, решение имеет существенную особую точку при условии 
.
При этом решение с существенной особой точкой определяется из уравнения (4.6). Т.е. по мере приближения к положению равновесия 
координаты 
происходит скачок решения, который является произвольным, так как решение в виде ряда Лорана содержит существенную особую точку и при приближении к особой точке, значение функции 
произвольно. При достижении положения равновесия остальными координатами 
, в случае, если положения равновесия кратные, наблюдается скачок решения на сходимость к другому положению равновесия. В результате решения в окрестности кратного положения равновесия, получаются ряды Лорана с существенной особой точкой.
Причем коэффициенты ряда Лорана могут оказаться многозначной функцией, так как согласно теореме 2 решение многозначно при кратных координатах положения равновесия, а значит в турбулентном режиме. Отметим, что ряд Лорана принимает любое комплексное значение в окрестности существенно особой точки, кроме быть может, одного значения, при комплексном аргументе.
При этом многозначное решение имеет вид ряда Лорана и может иметь существенно особую точку, а может быть многозначное решение имеет вид
, (4.11)
где величина t1 определится из начальных условий, а величины tk имеют постоянное значение. В случае ряда Лорана реализуется хаотическое решение, которое соответствует случаю 
теоремы 2 и происходит отскок от положения равновесия. В случае решения вида (4.9) реализуется случай 
, и положение равновесия достигается в комплексной плоскости. В первом случае решение является многозначной функцией. Во втором случае наблюдается стремление к положению равновесия с помощью многозначной функции при бесконечном значении V.
Конец доказательства.
В случае двукратного положения равновесия по переменной xN решение уравнения Навье – Стокса имеет вид
где величины 
однократные координаты положения равновесия, а величины 
двукратные положения равновесия.
Ситуация аналогична уравнению Шредингера, когда имеется в общем случае счетное количество решений. Это не удивительно, ведь уравнение Шредингера сводится к уравнению Навье – Стокса. Докажем это. Для чего запишем уравнение Шредингера и преобразуем
.
Разделив на массу 
, получим уравнение
.
Получим уравнение в частных производных, взяв градиент от обеих частей уравнения, введем действительную скорость по формуле 
.
.
Подставляя значение скорости в преобразованное уравнение Шредингера, получим
.
Получим трехмерное уравнение Навье – Стокса с давлением, соответствующим потенциалу. Но задача гидродинамики отличается от уравнения Навье – Стокса, полученного из уравнения Шредингера, уравнением неразрывности.
При этом можно провести аналогию между ламинарным однозначным режимом и классическим однозначным описанием тел. Между турбулентным режимом, имеющим счетное количество решений и квантовым описанием частиц, имеющих счетное количество решений. В случае турбулентного комплексного и ламинарного действительного режима между ними имеется граница, критическое число Рейнольдса. Аналогичная граница имеется между классическим и квантовым описанием частиц, соответствующая переходу энергии состояния из отрицательного в положительное состояние. В свою очередь уравнение Навье – Стокса должно иметь дискретные уровни энергии состояния турбулентного потока, должны реализовываться переходы между этими состояниями с выделением энергии или поглощением энергии.
Границу между классическим описанием частиц и квантовым можно провести, это переход к комплексному квантовому числу, или бесконечности главного квантового числа атома водорода. При этом бесконечное квантовое число атома водорода через нулевое значение выражения 
, где n главное квантовое число, становится мнимым, непрерывным. Ламинарному решению задачи гидродинамики, при котором имеется однозначное решение, соответствует волновая функция свободного движения, которая однозначна, при непрерывной энергии. Причем при большом квантовом числе система является квазиклассической, т.е. близко к граничному квантовому числу, равному бесконечности, система почти классическая.
Причем, имеется граница, между свободным решением, и решением, описывающим связанные состояния. Это нулевое значение энергии, также как и в нелинейных уравнениях в частных производных, существует граница между турбулентным комплексным решением и ламинарным действительным решением.
Библиографическая ссылка
Якубовский E.Г. ИССЛЕДОВАНИЕ РЕШЕНИЯ УРАВНЕНИЯ НАВЬЕ – СТОКСА // Научное обозрение. Реферативный журнал. – 2016. – № 1. – С. 46-80;URL: https://abstract.science-review.ru/ru/article/view?id=632 (дата обращения: 21.12.2024).