Проверяемый текст
[стр. 124]

124 р компоненты средней скорости воздуха (/ = 1,2, 3); / время; х)-пространственные координаты; и*— пульсационные скорости (1 = 1,2, 3); Р статическое давление воздуха; 5 энтропия; тс(Т эффективная вязкость; р динамическая вязкость; р, турбулентная вязкость; 5у дельта Кронекера; к кинетическая энергия турбулентности; Я полная энтальпия; И статическая энтальпия.
Для замыкания данной системы уравнений используется
полуэмнирическая модель турбулен тности, состоящая из двух уравнений переноса: кинетической энергии турбулентности 3(рк) 3(рМ;.1С) _ д д( дх.
дх.
ц + -^ сг ск к У дхп + 5к + +И, ди, ди, 1 ю дх.
1 рк + р, ди, ди, дх, дх.
-Р е.
(4.6) скорости диссипации кинетической энергии турбулентности д(ре) д(р и,г) д д1 + % дх,'( \ : дг— ■ 1 д*,1 + + г к 'е! И, ди ди.
-+ ' 4 ди: 2/ уох, 8хи дх, 3 ди, Рк + М,т ах.
си, / ох, -рсе2еК (4.7) ц, =рси —.
(4-8) г где си коэффициент к е модели турбулентности.
Дискретизация уравнений осущес твляется методом конечных объемов
[6 7].
Дискретизация расчетной области производится при помощи много
[стр. 1]

МОДЕЛИРОВАНИЕ АЭРОГАЗОДИНАМИЧЕСКИХ ПРОЦЕССОВ ПРИ ПРОВЕТРИВАНИИ ВЫРАБОТОК БОЛЬШОГО ПОПЕРЕЧНОГО СЕЧЕНИЯ Н.М.
Качурин, С.А.
Воробьев, А.Д.
Левин, Ф.М.
Ботов Рассмотрены аэрогазодинамические процессы при проветривании выработок большого поперечного сечения.
Обоснованы математические модели движения воздуха в выработках, имеющих большую площадь поперечного сечения.
Показано, что моделирование аэрогазодинамических процессов при проветривании выработок большого поперечного сечения основывается в общем случае на системе уравнений Рейнольдса, описывающей течение вязкого, сжимаемого теплопроводного газа в трехмерной постановке.
Ключевые слова: уравнение неразрывности, уравнение движения, математическая модель, вязкость, турбулентность, метод конечных объемов.
Моделирование аэрогазодинамических процессов при проветривании выработок большого поперечного сечения основывается в общем случае на системе уравнений Рейнольдса, описывающей течение вязкого, сжимаемого теплопроводного газа в трехмерной постановке, которая состоит из основных уравнений сохранения [1 4]: ^ + —(ри.
) = о, дt дх.]) ' д / \ д / * *\ дР „ (ри-) + (ри*иг) = -— + + (1) д дх.
дt / ди ди.
дх.
i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
V .
+ дх дх 2 ди, „ 2 „ 3^ эХ76* 3р5*к (2) д г тт\ дР —(рН)--+ дt дt дх,.
3 (ри.Я) = д дх í х^+Е дh \ i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
д дх.
и.
Г ди, ди,Л V дх, дх 2 } ди V дх.
Рг, дх.
^ + ^+ У 2 эх;зр5*к i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
_дк дх , 1 * * Н = п + — и и.
+ к, 2 ] 1 (3) (4) (5) где р плотность воздуха; и.
компоненты средней скорости воздуха (.
= 1, 2, 3); t время; х.
пространственные координаты; и*пульсационные скорости(.

= 1, 2, 3); Р статическое давление воздуха; £ энтропия; тей.
г 3 эффективная вязкость; т динамическая вязкость; т турбулентная вязкость; 5дельта Кронекера; к кинетическая энергия турбулентности; Н -полная энтальпия; h статическая энтальпия.
Для замыкания данной системы уравнений используется
полуэмпирическая модель турбулентности, состоящая из двух уравнений переноса: кинетической энергии турбулентности Э(рк) Э(ри7к) _ Э Э, Эй Эи, i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
+ Эх, Эх Эх] Эи 2 Эх 1Л А о Эк к у Эх, + ^ + Эх, 3 рк+т? Эи, i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Эх Эи,, I У Эх ре: (6) скорости диссипации кинетической энергии турбулентности Эх Эх е +— к "е! Э(ре) Э(ри7е) __Э_ Э, + __ Эи.
Эи —^+—^ ЭхЭх.
i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
т, о Эе е У Эх Эи 2 Эх, 3 рк+т, Эи, Эх Эи,, I У Эх,, + й + рСе2е i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
т _ рс к е (7) (8) где ст коэффициент к е модели турбулентности.
Дискретизация уравнений осуществляется методом конечных объемов.

Дискретизация расчетной области производится при помощи многоблочных,
неортогональных, адаптивных, структурированных сеток.
Каждая подобласть представляется в виде трехмерной матрицы сеточных узлов & j, ф, где 1 £ j, k £ ID.
В каждом сеточном узле определены все зависимые переменные.
Для описания распределения узлов внутри сеточной подобласти вводится понятие потокового элемента, который по своей сути является конечным элементов и на котором определены функции формы конечного элемента.
На рис.1 представлен гексаэдрический потоковый (конечный) элемент, определяемый восемью узлами.
Для потокового элемента с номерами узлов (ij,k)-(i+1j+1,k+1) область элемента может быть определена в терминах локальной неортогональной системы координат s, ^ u в виде записи через функции формы конечного элемента: х^,,,и) = N1х— + N2х.+1-,к + NзXi-+l,k + N4х.+1-+1,к + + N5Xi-,k+1 + N6Xi+l-,k+1 + N7Xi-+l,k+1 + N8Xi+l-+l,k+l; у^,г,и) = N^1--^ + N2Уi+l-,k + ^у^^ + N4Уi+l-+l,k + + N5Уi-,k+l + N6Уi+l-,k+l + Nту.+1М1 + NgУi+l-+l,k+l'} 2(5,г,и) = N12,^ + N22+1,^ + N32^+1^ + N42+1+1,к + + N52^+1 + N2+1^+1 + N72^+1^+1 + N^+1^+1^+1', где функции формы потокового элемента N записываются как N1 = 1/8(1 s)(1 г)(1 и); N2 = 1/8(1 + s)(1 г)(1 и); N3 = 1/8(1 5)(1 + г)(1 u); N,1 = 1/8(1 + s)(1 + г)(1 и); N5 = 1/8(1 5)(1 г)(1 + и); N6 = 1/8(1 + 5)(1 г)(1 + и); i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
N7 = 1/8(1 я)(1 + 1 + и); N8 = 1/8(1 + в)(1 + ¡)(1 + и).
Рис.
1.
Гексаэдрический потоковый (конечный) элемент Диапазон изменения локальных неортогональных координат 5, г, и: -1 <5, г, и<1.
Численная схема представляет собой полностью консервативный метод контрольных объемов.
Для определения контрольных объемов внутри расчетной области, дискретизированной регулярными массивами потоковых элементов, вводится понятие октантов.
Октант (квадрант в двумерной постановке) определен так, что каждый потоковый элемент состоит из восьми октантов, каждый из которых связан с одним из восьми узлов потокового элемента.
Для каждого потокового элемента октанты определены через диапазон изменения локальных координат s, ^ и следующим образом: 1) -1 < 5 < 0, -1 < г < 0, -1 < и < 0, 2) 0 < 5 < 1, -1 < г < 0, -1 < и < 0, 3) -1 < 5 < 0, 0 < г < 1, -1 < и < 0, 4) 0 < 5 < 1, 0 < г < 1, -1 < и < 0, 5) -1 < 5 < 0, -1 < г < 0, 0 < и < 1, 6) 0 < 5 < 1, -1 < г < 0, 0 < и < 1, 7) -1 < 5 < 0, 0 < г < 1, 0 < и < 1, 8) 0 < 5 < 1, 0 < г < 1, 0 < и < 1.
На рис.
2 представлен потоковый элемент, подразделенный на октанты.
Шесть поверхностей каждого из восьми октантов потокового элемента подразделяются на две группы так, что три поверхности октанта совпадают с поверхностями потокового элемента, а три оставшиеся находятся внутри потокового элемента.
Последняя группа формирует поверхности контрольного объема и связанные с ними точки интегрирования, на которых определяется интеграл по поверхности.
Рис.
2.
Подразделение потокового элемента на октанты (один октант вынесен) i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
После применения процедуры интегрирования с использованием теоремы Гаусса, которая предусмотрена в методе контрольных объемов, уравнения (1), (2), (3), (6), (7) примут вид: т р а и I +1 Р = (9) Э' .
и я Э-I ри^и I + ри.и^.
=-1Pdnl + SU1dи + Эг + I т Эи Эи V дх] Эх.
2 ди, 3' eff Эх, ~1] V ~1У 2 i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
dn I"\\РШи1 Ы^ I + IРи1Шп Эг =I Эг ' ЭТ т Э^ л 1— + —-Эх Рг.
Эх V з х з У + < и.
т Эи Эи, Эх Эх V з г 2 Эи i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
2 ГЭх, 5'3 ^к +1 SEd и Эk ч Эх dnJ.; —1 Iрк dи +1ри.
к dnJ = I Эг +1 и {\ / тг О Эк Эх i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
к у .
+1S ^ и+ тг Эи.
Эи.
—+ —3 V Эх.
Эх.
Эи 2 У Эхз 3 рк+тг Эи, Эх Эи,, I У Эх -Э I ре dи +1ри.е dnJ = I О i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Л Эе е У Эх.
dnJ +1 £е d и + (10) (11) dи-1 ре dи; (12) +К к 'е1 тг Эи Эи, + V Эх.
Эх.
Эи 2 i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Эх.
31 рк + тг Эи, Эх Эи,, I У Эх р^е 2е )d и, (13) где се1, се2 коэффициент к е модели турбулентности; dni произведение компоненты вектора внешней нормали на площадь октанта в декартовой системе координат.
Интеграл по поверхности отвечает за интегрирование потоков консервативных величин, интеграл по объему учитывает действие источника.
На рис.
3 представлен конечный объем и потоковый элемент.
Каждый конечный объем определяется двадцатью четырьмя билинейными поверхностями в трехмерной постановке или восемью линейными сегментами в двумерной постановке.
Интегральные уравнения (9) (13) записаны для каждого контрольного объема, полученного путем соединения середин противоположных сторон в каждом элементе.
Непрерывный интеграл по поверхности переводится в дискретную форму и оценивается через так называемые точки интегрирования.
Место положения точек интегрирования для одного потокового элемента в двумерном представлении показано на рис.
4.
и и 5 5 i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Г s и 5 5 5 и 5 и и и и 5 5 г i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Рис.
4.
Расположение точек интегрирования на двумерном билинейном элементе В трехмерной постановке потоковый элемент состоит из восьми октантов и двадцати четырех поверхностей, содержащих точки интегрирования.
В дискретной форме интегральные уравнения (9) (13) могут быть записаны в следующем виде: Уо1 .0 \ дг + ф (14) pVol Г (Л иг иг + = -^(ГАп, \р + s*Vol + 'р 'р +1 i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
и eff Э и Э и —L + —1 dXj Эх 2 ди, „ 2 „ --LI ,.,.--О--DO К 3^eff Эх, v 3 у А/7.
(15) ip pVol н-н At 0 Л = 1 i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
^ЭГ 4 Эh Эх Prt Эх Vol У V Л An r Р _ р° ^ ip j V j J J At ip ip U: i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
l^eff Гд и Эм;Л —L + —JЭх, Эх, V j 1 2 du, „ 2 „ Э£ .
^ J (16) pVol о Л к-К At +Z(P «¿4)»*=I 11+ ik о Эк i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Эх, к у гр Д/7 + к Vol + + Эг/ Эг/ —L + —¿ Эх Эхг Э и 2 Эл0 3V Р к+Ц, Эи, Эх, Эг/,.
i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
/ У Эх, Vol peVol (17) pVol — 8° ^ V У у с M-, i + — о Л Эе * у Эх/ Ли.
+ SeVol + i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
'Р + К С, el Эм Э м —L + —¿ Эх, Эх, V j 1 Э г/.
2 Эх7 3 v рк + ц, Эи, Эх Эг/, i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
/ у Эх,, РСе2£ (Vol, (18) где Vol величина контрольного объема; нижний символ «ip» точка интегрирования, суммирование производится по всем точкам интегрирования; Arij произведение компоненты вектора внешней нормали на площадь грани; At шаг по времени; верхний символ 0 определяет параметр, взятый на старом временном уровне; верхняя черта над членом, учитывающем влияние источника, обозначает осредненное значение на контрольном объеме.
Следует отметить: фундаментальное преимущество метода конечных объемов означает то, что потоки в точке интегрирования на соприкасающихся поверхностях соседних контрольных объемов равны, т.е.
поток, истекающий из одного контрольного объема и втекающий в прилегающий объем, идентичен.
Для аппроксимации по времени используется схема Эйлера первого порядка точности.
Такая аппроксимация не накладывает жестких ограничений на размер шага по времени Ко/ (19) М Стандартный подход метода конечных элементов (через функции формы конечного элемента) используется для оценки производных для всех диффузионных членов.
Например, производная по направлению х в точке интегрирования ¡р Эф Эх гр п ЭХ Ф„ • (20) 'р Суммирование производится по всем функциям формы элемента.
Производная от функции формы в декартовой системе координат может быть выражена в виде локальных производных функции формы через якобиан матрицы преобразования координат.
i Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Для повышения устойчивости схемы оценка градиента производится в точках интегрирования, расположенных на пересечении поверхности интегрирования с ребром потокового элемента (линейно-линейная интерполяция).
Применение такого подхода снижает степень аппроксимации ЭА7Э{.у, /, и) до первого порядка точности как для ортогонального, так и для искривленного элемента.
В противном случае при оценке градиента в точках интегрирования, расположенных согласно рис.
4 (три линейная интерполяция), можно обеспечить второй порядок точности для ортогонального элемента и первый для искривленного.
В случае качественной сеточной дискретизации можно использовать метод со вторым порядком интерполяции диффузионных членов, что снижает схемную вязкость и позволяет более точно и физически обоснованно описывать происходящие процессы.

[Back]