Равномерный и эквидистантный способы построения упорядоченной расчетной сетки для задач газовой динамики Текст научной статьи по специальности «Физика»

УДК 533.6:519.63

РАВНОМЕРНЫЙ И ЭКВИДИСТАНТНЫЙ СПОСОБЫ ПОСТРОЕНИЯ УПОРЯДОЧЕННОЙ РАСЧЕТНОЙ СЕТКИ ДЛЯ ЗАДАЧ ГАЗОВОЙ ДИНАМИКИ

©2008 А Н. Ефимов, В. Н. Матвеев

Самарский государственный аэрокосмический университет

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

Газовая динамика, исследование, погрешность, расчетная сетка, разбиение, построение

В настоящее время широкое распространение получили численные методы исследования, применение которых требует разбиения пространства на элементарные подобласти (элементы) и совместного исследования процессов в них. При этом форма элемента может быть выбрана произвольно, исходя из удобства записи уравнений для описания происходящих в нем физических процессов. На подобном подходе основаны широко распространенные в инженерной практике методы конечных элементов[1] и контрольных объемов

[2]. Но при этом возникает задача разбиения исходной геометрической модели на элементы - создание расчетной сетки. Для уменьшения погрешности расчета на форму элемента накладывается некоторый ряд требований и ограничений [6]:

- в качестве формы элемента желательно использовать треугольник или четырехугольник;

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

- для достижения наибольшей точности расчета желательно, чтобы внутренние углы при вершинах элементов лежали в пределах 45° <а <135° .

Известны два основных подхода к построению сетки.

1. Использование неупорядоченной сетки. При данном подходе сетка создается произвольным образом. В качестве формы элемента обычно используется треугольник (те^юдо в случае трехмерной задачи) ввиду того, что данная геометрическая фигура всегда является выпуклой. В то же время форма всей расчетной области может быть произ-

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

[3]. Сложность создания сетки не зависит от формы расчетной области. Основным недостатком использования неупорядоченной сетки является существенная численная диффузия [4], увеличивающая погрешность расчета. Причем эта погрешность имеет нерегулярный, вероятностный характер. Неупорядоченная сетка применяется обычно в случаях, когда сетка упорядоченного типа просто не может быть создана либо когда допускается повышенная погрешность расчета.

Использование упорядоченной сетки. В

данном случае сетка строится в соответствии с жесткими правилами. Число элементов задается точно, как и закон их распределения по расчетной области. В качестве формы элемента обычно используется четырехугольник (парадлелепипед в случае трехмерной задачи). Сеткой подобного типа легче управлять, но процесс ее создания чрезвычайно трудоемок, особенно при сложной геометрии исходной области. Этот способ построения сетки применяется в основном в случаях, когда точность результатов расчета крайне важна. Использование упорядоченной сетки желательно при расчете потоков в соплах и турбомашинах, при определении параметров газа, обтекающего профиль крыла самолета, и сверхзвуковых течений газа.

Вестник Самарского государственного аэрокосмического университета, №3, 2008

X (l ) = (l);

В процессе создания упорядоченной сетки возможно применение двух способов-алгебраического и дифференциального. Оба эти способа определяют координаты узлов элементов и отличаются типом решаемой системы уравнений. Реализация алгебраического наиболее легка, но он подходит для построения сетки достаточно простых областей. Несмотря на это, № широко используется в таких коммерческих системах, как Ашу8, 81агСВ, СБХ и др. Дифференциальный метод [5] основа на численном решении системы эллиптических уравнений, подходит для более сложных областей и используется в пакете Б2-ТигЬо[8].

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

Равномерный способ предназначен для разбиения областей со сторонами произвольной формы.

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

При использовании предложенных способов к расчетной области предъявляются следующие требования:

- количество элементов вдоль противолежащих образующих должно быть равным;

- противолежащие образующие не должны пересекаться.

Равномерный способ.

В качестве исходных данных используются следующие:

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

Y (l ) = .к, у (l),

где /-дайна сплайна от его начала до произвольной точки с координатами (х,у);

.к,х(/), ЯьуА)-сплайн-фикции; к-номер образующей.

Образующие попарно пересекаются в точках Л,Б,С,Б (рис. 1). При этом образующие СЭ и СА начинаются в точке С, образующая АВ-в точке А, образующая ВЭ-в точке Э. Тогда

в„(0) = Хс 0= Хс Г53,(0) = Хл I^0=хв

Ы 9=Гс ’{.‘У0) =Гс ’[5Ц9=ул к (9=Ъ

Количество элементов вдоль образующих Щ: И! = N3 и N2= N4.

о Закон распределения узлов

вдоль каждой образующей, определяющий расстояние вдоль к-го сплайна от его начала до п-ого узла, описывается зависимостью /(п ) = Тк (п). Здесь п-номер узла вдоль образующей.

Функция /(п) обладает следующими свойствами:

/( 0 ) = 0;

/(N ) = Ьк, где Ьк- длина к-ой образующей.

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

/( п ) = —п

Процесс построения расчетной сетки заключается в создании матрицы Я узлов размером П х V и определении ее эле ментов:

о . Я 0 . п, 0 • \ 0^

Я = Я . т . Я п, т • Я, п, т

Л V • • ЯпV • Яп,V у

где к=Щг + 1, V=N2+1, {X

Я .=■

I, ]

-элемент матрицы, кото-

рый соответствует узлу, заданному вектором координат. Определение элементов матрицы Я проходит в три этапа. На первом этапе на-

ходятся элементы Я0>0,Яп>0,Я0^,Яп^.- координаты узлов, расположенных в точках С,Б,Л и Б.

Точке С соответствует элемент Хх ( 0)1

Я0,0 -

Яп,0 -

Яс.=

( 0)\ К (Ь)]

К (ь А

^3,х (0)1 . ( 0)\

точке Б- элемент

точке А- элемент

а точке Б- элемент

Кх (Ь3 )}

^ [.зу ( ьъ )\

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

Ххв(п))

Я 0 =

п ,0

п)).

>,1 < п<к- 1,

а б

Рис. 1. Создание расчетной сетки равномерным способом: а-исследуемая область, б- сетка

,т

К,, -

ЯЛ ,т

► ,1 < т<, — 1.

•,1< п<Ь — 1,

•,1< т<, — 1.

на образующей СА- узлы с координатами

'.2х(Т2 (т))1

.-ЛТЛ т))\

на образующей ЛБ- узлы с координатами

'Хз,х(Тз( п))' .3 ,г(Тз( п))\

на образующей ЭВ- узлы с координатами МТ4( т))

Х4,(Т.( т)\

На третьем этапе создаются внутренние узлы Яп, т. Для этого доится восьмиугольник АБВОЭНСЕ с прямолинейными сторонами (рис.2), для чего выбираются точки, координаты которых совпадают с координатами имеющихся узлов:

Н = Яп 0 ,Е = Я0 т,^ = Яп ,0 = Я, т. Значение

п, 0 7 0, т7 п, V 7 п, т

п изменяется в диапазоне 1 < п < й — 1 и определяет положение точек Н и V, в то время как т изменяется в диапазоне 1 < т < V — 1 и

определяет положение точек Е и О. Суммарное количество восьмиугольников равно (Л - 2 - 2).

Строятся четыре дополнительные точки N1,K1,N2,K2- точки пересечения прямых, проходящих через стороны многоугольника АБВОЭНСЕ:

СНгЛЕ,К=ЕБгНО

(ис. 2).

Координаты точки Е, которая является точкой пересечения двух прямых АВ и СЭ фис. 3), заданных координатами двух точек, могут быть найдены по формулам

(Ул-е-Хл)~{Ус-/-Хс)

Уе=-

при ^ /. Если g=f , то следует принять

Ге= 106, Хе=Уе-^Ус-^Хс ,

где

g=■

Х„-Хл

У -У

g = 106,

если У„Ф У л

если УБ=УЛ

/ =ХГ—хС, если У0 ФУС

У -У

■'D С

если Уа= УС

[/ = 106,

В дальнейшем строятся четыре дополнительные ТОЧКИ І1,І2,І3,І4 фис.2), являющиеся точками пересечения соответственно пересечением следующих прямых: Т1Н и

КіЕ, КіЕ И Ы2Е, Ы2Е и КО, К2ви ЫН. Координаты точек І1,І2,І3,І4 определяются по тому же правилу что и координаты точки Е на рис. 3.

Рис.2. К определению дополнительных точек І1^І6

и Я т

точки пересечения двух прямых

Методом линейной интерполяции определяются координаты дополнительных точек /5 и 16. Точка 15 расположена на отрезке, проходящем между точками 11 и 12, а точка /6- между точками /3,/4. Их координаты могут быть найдены по формулам:

хг = хг -і і--\+х \

п

н:

У = УІ,

X = х

і6 1 з

Уі = У І

її,

1-п

н

+У •п;

І2 н

/

1-

п

+ Х1

п

Т

1-п\+У,-п.

к) 4 к Искомый узел Яп тн^одатся на отрезке, соединяющем точки /5 и /6. Его координаты могут быть найдены методом линейной интерполяции по формуле

Я =

п ,т

ХТ

1-

т

+ X1

т

V

1 т

После заполнения матрицы Я можно приступить к формированию элементов сетки. Процесс формирования элементов сетки описан ниже.

Эквидистантный способ

Исходные данные идентичны используемым в равномерном способе. Вместе с тем, к расчетной области предъявляется ряд дополнительных требований (рис. 4):

- образующие ЛБ и СБ должны быть эквидистантны;

- при построении сетки между эквидистантными кривыми обычно трудно подобрать такие функции распределения узлов по эквидистантным кривым Т1 (п) и Т3 (п), чтобы ребра элементов были нормальны к образующим. Поэтому предполагается, что функция Т3(п) распределена узлов по образующей АВ неизвестна;

- образующие области взаимно перпендикулярны в точках пересечения ЛС±ЛБ, ЛС1СП, БП1ЛБ, БП1СВ длины образующих АС и ВЭ равны. Для данных образующих заданы одинаковые функции распределения узлов Т2(т) = Т4(т), где 0 < т < V. Это обуславливает одинаковое распределение соответствующих узлов вдоль данных образующих и равенство длин сторон соответствующих элементов, расположенных как на образующих, так и внутри области.

X

Создание матрицы Я с помощью эквидистантного метода осуществляется в два этапа.

На первом этапе определяются элементы матрицы

Rn,0 = '

где 0 < n<h.

КИ")).

На втором этапе создаются остальные элементы матрицы Япт в диапазоне

0 < п<к, 1 < т < V послойно от ЛС к ББ. Для этого восстанавливаются нормали к образующей СБ в каждом ее узле с ко ординатами Яп,0. Затем определяются координаты

узлов, лежащих на этих нормалях от Яп 1 до

Яп ^, по формулам

XRn 0 + Tl(m) COS(a) YR + Tl{m) sin{a)

где

a = arctg

dshy(iЯ

1 y dl

dShx(l)

dS, (l)

, при —--------------Ф 0

dl

и

dl _ dShx (l)

n

a = — при

2 dl

= 0.

Создание элементов сетки на базе полученных узлов

По полученным узлам (как в случае использования равномерного, так и эквидистантного способов) создается матрица элементов Е размером (к-1)х(V-1). Каждый элемент матрицы Е соответствует элементу

сетки. Элемент сетки представляет из себя четырехугольник, вершинами которого являются узлы. Нумерация узлов зависит от конкретного пакета программ, для которого строится сетка фис. 5). В пакетах Ansys, StarCD принята нумерация узлов против часовой стрелки, в пакетах Fluent и CFX -послойная нумерация узлов. В зависимости от варианта нумерации узлов элемент матрицы записывается в виде

Е =

m n

Rn

п.

R

R

)

R

"n+1, m

"n+1, m+1

'n, m+1

для варианта 1,

и

E =

m n

Rn

R

R

}

R

n+1, m

'n, m+1

n+1, m+1

для варианта 2.

Полученную матрицу узлов Я можно в случае необходимости использовать также и для создания треугольной упорядоченной сетки. Для этого каждый плоский элемент разбивается на два треугольных.

Разбиение выполняется между противолежащими узлами. Существует два варианта подобного разбиения (рис. 6) - отрезками АС и ВЭ. Искомым вариантом разбиения является вариант с наименьшей длиной секущего отрезка [3], в данном случае вариант а.

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

Рис. 5. Варианты нумерации узлов плоского элемента: a -Ansys, StarCD; б - Fluent, CFX

C

а б

Рис. 6. Варианты разбиения четырехугольного элемента на треугольные

Для этого расчетная область делится вспомогательными образующими на четырехугольные зоны, в которых создается расчетная сетка с использованием равномерного или эквидистантоного способов. При этом эквидистантный способ используется для построения сетки вблизи твердых стенок в области повышенных нормальных градиентов скоростей, так как полученная сетка лучше удовлетворяет критериям качества [6]. Между тем, требование эквидистантности двух противолежащих образующих в общем случае не позволяет использовать данный способ для построения сетки всей расчетной области. Оставшаяся расчетная область разбивается сеткой с использованием равномерного способа, не веющего подобных ограничений.

Приведенные в статье способы построения сетки реализованы в препроцессоре ТигЬоМе8к и используются для построения расчетной сетки межлопаточных каналов турбин. При этом эквидистантный способ используется для построения сетки в области пограничного слоя, а с помощью равномерного способа создается сетка в области ядра потока.

Библиографический список

1. Сабоннадьер Ж.-К ,Кулон Ж.-Л. Метод конечных элементов и САПР / Пер. с франц. - М.: Мир, 1989. -190 с.

2. С.Патанкар. Численные методы решения задач теплообмена и динамики жидкости - М.: Энергоатомиздат, 1984. - 152с

3. М. Щербак. Алгоритм автоматической генерации двумерной конечно-элементной сетки.

rhttp://home.onego.ru/~scherbak/1.31.05.2005.

4. Эффективное решение конвективного уравнения переноса современными численными методами: методические указания по изучению спецкурса «Решение задач на ЭВМ» / Перм. ун-т; Сост. А.В. Штраубе. -Пермь, 2003. - 28 с.

5. Основы проектирования турбин авиадвигателей /А.В.Деревянко, В.АЖ^авлев, В.В.Зикеев и др. Шодред С.З.Копелева - М.: Машиностроение, 1988. -328 с.

6. ANSYS Modeling and Meshing Guide. [http: / / strelka. ftf2.tsu. ru/~sid/Linux/ Ansys_do c/ /Ansys html/guide 55/g-mod/GMODToc.htm1. 20.06.2005.

7. Боглаев ЮЛ. Вычислительная математика и программирование: учеб. пособие для студентов втузов. - М.: Высш.шк., 1990 -544с.

8. EZ-Turbo.[http://www.cd-adapco.com/ news/16/es-turbo.htm1. 20.06.2005.

References

1. Sabonadier J-C., Kulon J-L. Finish Elements Method and CAD. Transl. from French. Moscow: “Mir”, 1989.

2. S.Patankar. Numerical solution methods of heat exchange problems and fluid dynamics. Moscow: “Energoatomizdat, 1984.

3. Scherbak M. Algorithm of 2D FE grid automatic generation. [http://home. one-go.ru./~scherbak/]. 31.05.2005.

4. Schtraube A.V. Effective solution of convective equation by modern computing: learners guide for “Problem solving using computer”. Perm: Perm’s University, 2003.

5. Derevyanko A.V., Zhuravlev V.A., Zik-heev V.V. Aircraft Turbines Design Basis. Moscow: ”Mashinostroenie”, 1988.

6. ANSYS Modelling and Meshing Guide. [http: / / strelka. ftf2.tsu. ru/~sid/Linux/

Ansys_do c/______/Ansys_html/ guide_55/g-

mod/GMODToc.htm]. 20.06.2005.

7. Bogkaev Y.P Calculus Mathematics and Programming. Tutorial. Moscow: ”Vysshaya Shkola”, 1990.

8. EZ-Turbo. [http://www.cd-adapco.com/ news/16/es-turbo.htm] 20.06.2005,

UNIFORM AND EQUIDISTANT WAYS OF STRUCTURED COMPUTATION GRIGD CREATION FOR CFD PROBLEMS

©2008 A. N. Yefimow, V. N. Matveev

Samara State Aerospace University

Two ways of structured plain algebraic computation grid for tetragonal regions with curvilinear sides are introduced. Uniform way can be used in tetragonal regions with arbitrary side shape. Equidistant way is used in regions where two opposite sides are equidistant.

Gas dynamics, research, inaccuracy, computational grid, meshing, modeling

Информация об авторах

Ефимов Алексей Николаевич, аспирант кафедры теории двигателей летательных аппаратов Самарского государственного аэрокосмического университета. E-mail: tdla@ssau.ru. Область на^ных интересов - турбомашины, газовая динамика, системы охлаждения.

Матвеев Валерий Николаевич, доктор технических наук, профессор, заведующий кафедрой теории двигателей летательных аппаратов Самарского государственного аэрокосмического университета. E-mail: tdla@ssau.ru. Область научных интересов: рабочие процессы микротурбомашин.

Efimov Alexei Nikolaevich, postgraduate of Samara State Aerospace University “Theory of Aircraft Engines” department. E-mail: tdla@ssau.ru. Area of research: turbomachines, gas dynamics, cooling systems.

Matveev Valeriy Nikolaevich, Doctor of Engineering Science, professor, head of SSAU “Theory of Aircraft Engines” department. E-mail: tdla@ssau.ru. Area of research: operation processes in microturbomachines.