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

Читаемое

Обновления Oct-2018

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

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

где pdpi - отношение давлений (PRAT в программе SHOCK). Из уравнений (14.169) также следуют граничные условия при х = О yl 1.0. Решите задачу при NX = 101, NT = 170 и DT = 0.100 по схеме Лакса - Вендроффа с искусственной вязкостью. Общий характер решения можно сравнить с численным решением Сода [Sod, 1978]. Точное решение этой задачи можно найти у Липмана и Рошко [Liepmann, Roshko, 1957]. Наибольшую сложность представляют расчет контактного разрыва и ударной волны.

14.9. Получите решение задачи о распространении умеренной {pjp2=2.b, рис. 14.18) и сильной (pjp2 = 10, рис. 14.25) ударных волн по схеме FCT при следующих значениях диффузионных и антидиффузионных параметров в (14.93):

(1) т1о = 0.125, Л1=0. Л2 = 0,

(2) т1о = 0.500, Л1=0, Л2 = 0,

(3) Ло=1/3, Л1 = 1/3, Л2==-1/6,

(4) Ло=1/6, Л1 = 1А Л2 = -1/3. Сравните профили ударных волн, полученные в различных вариантах.

14.10. Примените схему FCT к решению задачи 14.8 об ударной трубе. Сравните полученные решения с решением, найденным по схеме Лакса - Вендроффа с искусственной вязкостью (задача 14.8). Особое внимание уделите профилям ударной волны и контактного разрыва.

Трансзвуковые невязкие течения (§ 14.3)

14.11. Течение у тонкого крыла, задаваемого уравнением у = т{х - \ ), -1 Jc 1, в области, изображенной на рис. 14.34, описывается трансзвуковым уравнением малых возмущений (14.130). Граничные условия приведены на рис. 14.34. Получите численное решение при т = 0.1, == 0.8 и 0.9


Рис. 14.34. Модельная задача для уравнения потенциала.

на сетке 41 (NX)X 21 (NY). Для конечно-разностной дискретизации используйте формулу (14.135). Для расчета используйте итерационный алгоритм SOR.

14.12. Решите задачу, представленную на рис. 14.34, используя дискретное представление уравнения полного потенциала (14.129), (14.138), (14.139),



Для введения диссипативного механизма в сверхзвуковой области используйте (14.143). Граничное условие Фу = fx на АВ можно трактовать как вдув газа с заданным профилем скорости.

Для решения дискретных уравнений используйте метод приближенной факторизации, подобный (14.156), (14.157), но с использованием метода установления (14.137). Получите решение на сетке 41 (NX)X 21 (NY) при т=0.1 и Моо = 0.8 и 0.9. Сравните полученные решения с решениями задачи 14.11.

14.13. При помощи центральных разностей проведите дискретизацию уравнения Лапласа дФ1дх + дН/ду = О в прямоугольной области О 1, О 1/1 с граничными условиями 0(0, г/) = со8(0.5яг/), 0(1, /) = = ехр(я/2)со5(0.5яг/), ф(х, 0)= ехр(0.5ял:), ф(х, 1) = 0.

Получите решение дискретных уравнений

(1) по схеме приближенной факторизации (14.156), (14.157),

(2) по схеме приближенной факторизации в сочетании с многосеточным подходом (14.156), (14.157), (14.162)

и сравните эффективность методов, принимая во внимание число итераций необходимое для получения сходящегося решения и примерное число операций. Точное решение этой задачи Фех = ехр(0.5ял:)со8(0.5яг/).

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



Глава 15 Течения в пограничном слое

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

u=u=0


Рис. 15.1. Течение в пограничном слое.

{рис. 15.1). Уравнение для нормальной составляющей импульса может быть заменено условием постоянства давления в этом направлении. Если распределение давления известно, то уравнения, описывающие такое течение, перестают быть эллиптическими, что позволяет разработать весьма эффективные одно-проходовые маршевые алгоритмы (в направлении х на рис. 15.1) для их решения.

В пограничных слоях имеют место большие градиенты скорости в направлении, нормальном к поверхности. Поэтому желательно использовать преобразования координат, для которых эти градиенты стали бы меньше в новой системе координат. Наиболее эффективные преобразования описаны в § 15.2 и 15.3. Кроме того, для получения хорошего разрешения вблизи стенки используются сетки, шаг которых геометрически возрастает в направлении нормали к поверхности (п. 15.1.2).

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



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 - 2018 Prom Izhora
При копировании текстов приветствуется обратная ссылка