Выпускная квалификационная работа Двухфазное течение в узком канале с источниками на границе (стр. 6 )

Как видно из Таблицы 1 при и условие устойчивости не выполняется при числах Махаи . Расчет был проведен и при данных значениях решение “разваливалось”. Результаты остальных расчетов приведены на рисунках 2-38. На каждом графике представлены три кривые которые соответствуют начальному распределению, распределению после отскока от стенки и промежуточному значению.

Рис.3:Распределение давления при числе Маха = 1.1,

Рис.4:Распределение давления при числе Маха = 2,

Рис.5:Распределение давления при числе Маха = 5,

Рис.6:Распределение плотности при числе Маха = 1.1,

Рис.7:Распределение плотности при числе Маха = 2,

Рис.8:Распределение плотности при числе Маха = 5,

Рис.9:Распределение скорости при числе Маха = 1.1,

Рис.10:Распределение скорости при числе Маха = 2,

Рис.11:Распределение скорости при числе Маха = 5,

Рис.12:Распределение энергии при числе Маха = 1.1,

Рис.13:Распределение энергии при числе Маха = 2,

Рис.14:Распределение энергии при числе Маха = 5,

Рис.15:Распределение давления при числе Маха = 1.1,

Рис.16:Распределение давления при числе Маха = 2,

Рис.17:Распределение давления при числе Маха = 5,

Рис.18:Распределение плотности при числе Маха = 1.1,

Рис.19:Распределение плотности при числе Маха = 2,

Рис.20:Распределение плотности при числе Маха = 5,

Рис.21:Распределение скорости при числе Маха = 1.1,

Рис.22:Распределение скорости при числе Маха = 2,

Рис.23:Распределение скорости при числе Маха = 5,

Рис.24:Распределение энергии при числе Маха = 1.1,

Рис.25:Распределение энергии при числе Маха = 2,

Рис.26:Распределение энергии при числе Маха = 5,

Рис.27:Распределение давления при числе Маха = 1.1,

Рис.28:Распределение давления при числе Маха = 2,

Рис.29:Распределение давления при числе Маха = 5,

Рис.30:Распределение плотности при числе Маха = 1.1,

Рис.31:Распределение плотности при числе Маха = 2,

Рис.32:Распределение плотности при числе Маха = 5,

Рис.33:Распределение скорости при числе Маха = 1.1,

Рис.34:Распределение скорости при числе Маха = 2,

Рис.35:Распределение скорости при числе Маха = 5,

Рис.36:Распределение энергии при числе Маха = 1.1,

Рис.37:Распределение энергии при числе Маха = 2,

Рис.38:Распределение энергии при числе Маха = 5,

2.3                ИСКУССТВЕННАЯ ВЯЗКОСТЬ

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

- тензор напряжений;

- гидростатическое давление;