Разделы сайта

Читаемое

Обновления Mar-2024

Промышленность Ижоры -->  Динамика жидкости: уравнения 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 [ 91 ] 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

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

15.4.L Квазихарактеристическое поведение

Из уравнений (15.84)-(15.86) следует, что вдоль направлений X и Z существует лишь конвективный перенос; вдоль направления у - конвекция и диффузия. Формально (т. е. согласно процедуре, описанной в гл. 2) уравнения (15.84) - (15.86)

Область зависимости

Область влияния


Поверхность гела

Рис. 15.21. Области влияния и зависимости.

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

д д ,

(15.87)

Квазихарактеристическое направление определяется проекцией линии тока на плоскость, касательную к поверхности тела, т. е. dz/dx = w/u. Но w/u изменяется поперек пограничного слоя. Поскольку при заданных (х, г) возмущения в направлении у распространяются по всей области, необходимо определить области зависимости и влияния (рис. 15.21). Возмущения в некоторой точке Р внутри пограничного слоя будут воздействовать на всю область, расположенную вниз по потоку и ограниченную плоскостями, проходящими через нормаль к телу, проведенную через точку Р; причем одна плоскость проходит через предельную линию тока на твердой поверхности, а другая



через проекцию линии тока, расположенную на внешней границе пограничного слоя.

Гиперболичность задачи в плоскости (дс, z) существенно ограничивает область, в которой по заданным начальным данным может быть найдено решение. Если начальные данные Uo(y) и Woiy) определены на поверхности Si, перпендикулярной стенке (рис. 15.22), решение вниз по потоку может быть


Рис. 15.22. Начальные и граничные данные.

получено лишь в области, ограниченной квазихарактеристиками S2 и S3 (и твердой стенкой и внешней границей пограничного слоя). Квазихарактеристические направления S2 и S3 определяются либо направлением предельных линий тока на стенке, либо направлением линий тока на внешней границе слоя. Попытки определить решение слева от S2 или справа от S3 (если смотреть против потока) нарушают принцип влияния Рэтца [Krause, 1973], поскольку локальное решение в этих областях определяется данными, лежащими вне плоскости Si. Как отмечено в п. 2.2.3, если определить дополнительные граничные условия, например на поверхностях S4 или S5, решение можно получить в любой области вниз по потоку, ограниченной начальными и граничными условиями.

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

Чтобы численное решение было устойчиво, должно выполняться условие Куранта -Фридрихса -Леви (КФЛ) (п. 9.1.2),



согласно которому расчетная область зависимости в плоскости (x,z) должна включать в себя физическую область зависимости, определяемую уравнениями в частных производных. Это существенно ограничивает относительные размеры шагов Ах и Az [Krause, 1973].

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



Рис. 15.23. Точки сетки, задействованные в схеме; (а) схема Кранка - Николсона, (Ь) схема загзаг Краузе.

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

Прямое обобщение схемы Кранка - Николсона на случай трех пространственных переменных состоит в центральной дискретизации в точке {k-1/2, п+1/2) в плоскости z (рис. 15.23) и в построении двух трехдиагональных систем разностных уравнений для и и w из уравнений х- и г-компонент импульса вдоль линий сетки (k, n+l), т. е. в направлении у. Из условия КФЛ следует, что скорость w должна быть положительна и Ах ограничено условием \wAx/uAz\l. Узлы сетки для трехмерной схемы Кранка - Николсона показаны на рис. 15.23.

Проход по Z осуществляется в направлении увеличения k. Если W везде меньше нуля, проход по z осуществляется в направлении уменьшающихся значений k. Таким образом, при-



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 [ 91 ] 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка