Переход к недивергентным формам записи системы уравнений Навье-Стокса

При использовании дивергентной формы записи

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

Рассмотрим новый вектор состояния потока , взаимообратный к вектору . Перейдем в уравнении (8) к вектору по формуле: . Здесь матрица есть матрица Якоби, обратная к которой . Умножим уравнение (8) на матрицу . Получим новую недивергентную форму представления системы уравнений

(см.[7]):

. (13)

Рассмотрим более подробно производную от вектора потока. Сам вектор потока равен . Тогда

.

Например, если , то , тогда .

Теперь полученное выражение для производной вектора потока подставим в уравнение (13):

и после преобразований получаем квазилинейное уравнение второго порядка (14)

Здесь .

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

1. В случае, когда новый вектор потока равен

система уравнений

с искомым вектором преобразуется в

систему уравнений, для нахождения которой:

1) сначала находим матрицу Якоби и матрицу, обратную к ней

тогда

,

2) дальше надо находить производные от вектора потока

(см.уравнение (11)): и И после вычислений получаем:

 

3) Для нахождения матрицы В недивергентной формы записи системы имеем:

.

Аналогично находится матрица С= и вектор правых частей .

Таким образом, в одномерном случае уравнение (14) с вектором течения получается следующее Здесь , а матрица В вычислена раньше.

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

(15)

Учитывая уравнение состояния , можно заменить производную давления во втором уравнении системы (15) по формуле . Тогда система преобразуется в следующую:

(16)

Выписывая теперь матрицу коэффициентов системы (16) при первых частных производных вектора течения , получаем

.

Матрица совпадает с вычисленной раньше матрицей В.

Из системы (16) также можно выписать и другие коэффициенты:

.

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

.

Здесь коэффициенты и символы Кронекера. Другие коэффициенты недивергентной системы равны:

где диссипативная функция, вычисляемая по формуле:

2. Рассмотрим другой случай, когда новый вектор течения ищется в переменных: , . В этом случае меняются коэффициенты и . Из уравнений состояния и находим ,

. Заменяя градиент давления во втором уравнении системы (16), получаем: . В уравнении энергии заменяем производные и, используя уравнение неразрывности преобразуем его к виду:

. Тогда для нового вектора течения матрица

. Здесь = .

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

Здесь