www.volsu.ru
МАТЕМАТИКА
DOI: http://dx.doi.oгg/10.15688/jvolsu1.2015.4.1
УДК 517.951, 519.632 ББК 22.161, 22.19
ТРИАНГУЛЯЦИЯ
_____ ___________О
ПРОСТРАНСТВЕННЫХ ЭЛЕМЕНТАРНЫХ ОБЛАСТЕЙ 1
Алексей Александрович Клячин
Доктор физико-математических наук, заведующий кафедрой математического анализа и теории функций, Волгоградский государственный университет klyachin-aa@yandex.ru, matf@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Анжелика Юрьевна Беленикина
Студентка Института математики и информационных технологий,
Волгоградский государственный университет
matf@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
га к к к ■х к
ч
Аннотация. В работе приводится построение триангуляции пространственных областей специального вида и дается нижняя оценка минимального
§ двугранного угла ее тетраэдров.
Ключевые слова: триангуляция, тетраэдр, двугранный угол, элементарная область, разбиение области, условие Липшица.
Введение
й Пусть задан конечный набор точек {Д}™! в пространстве И3. Триангуляцией
^ данного набора точек называется совокупность невырожденных тетраэдров {Т.,-}м
<
х удовлетворяющих условиям:
? 1) любая точка Рг является вершиной хотя бы одного тетраэдра ^;
3 } = !'
^ 2) каждый тетраэдр Т^ содержит только четыре точки из данного набора, явля-
@ ющиеся вершинами этого тетраэдра.
Пусть П ограниченная область в R3. Триангуляцией области П называется триангуляция произвольного конечного набора точек, лежащего в замыкании области П.
В настоящее время для плоских областей имеется довольно большое число различных алгоритмов триангуляции, описание которых можно найти, например, в работах [1; 5; 8-10]. Многие из приведенных там алгоритмов подходят, после соответствующих изменений, и для трехмерных областей. Описание некоторых прямых и итерационных алгоритмов в пространственном случае имеется в работах [2;3]. Как правило, в перечисленных работах не исследуется математически строго качество построенной триангуляции и не выводится никаких количественных оценок ее числовых характеристик.
В настоящей работе мы описываем достаточно простой в реализации алгоритм триангуляции пространственной области специального вида, который относится к так называемым прямым методам. Для оценки качества построенной триангуляции мы используем минимальный двугранный угол, вычисленный среди всех тетраэдров триангуляции. Именно мы приводим нижнюю оценку этого угла через геометрические характеристики области и параметры ее разбиения на ячейки. Стоит отметить, что условие пустой сферы, то есть условие Делоне, не обеспечивает качественной триангуляции (в отличие от триангуляции на плоскости). Например, в работе [6] строится пример триангуляции пространственной области, удовлетворяющей условию Делоне, но не аппроксимирующей с достаточной точностью градиент гладкой функции. Это ограничивает область применения таких триангуляций в различных вычислительных задачах, например, при численном решении уравнений с частными производными. В то же время более сильные условия на триангуляцию, которые обеспечивают необходимую точность вычисления площади, дают равномерную сходимость приближенных решений уравнения минимальной поверхности [4]. Для пространственных триангуляций (см., например, [7]) отделимость от нуля двугранных углов в тетраэдрах позволяет показать, что погрешность вычисления градиента гладкой функции стремится к нулю при стремлении к нулю максимального диаметра тетраэдров.
1. Основные результаты
Рассмотрим элементарную область П С R3, которая имеет вид
П = {(x,y,z) : а < х < b, с < у < d, <р(х, у) < z < ф(х, у)} ,
где ip(x,y) и ф(х,у) — заданные на прямоугольнике [a,b] х [c,d] функции, удовлетворяющие условию Липшица.
Далее полагаем
Mo = inf (ф(х,у) - <р(х,у)) > 0,
(х,у)£[а,Ь] X [c,d\
Ml = sup (ф(х,у) - (р(х,у)) < ж.
(x,y)£[a,b] X [c,d]
Будем считать, что для функций <р(х,у) и ф(х,у) найдется постоянная L < ж такая, что выполняются неравенства
) - ф",у")| < L^(x> - х'')2 + (у' - у1 ')2, 1Ф(х',у') - ф(х",у")1 < L^j(х' - х")2 + (у' - у'')2 ISSN 2222-8896. Вестн. Волгогр. гос. ун-та. Сер. 1, Мат. Физ. 2015. № 4 (29) 7
для любых х',х" Е [а,Ь],у',у" Е [с,в\. Далее опишем построение триангуляции области П. Пусть а = хо < х\ < х2 < ... < хп = Ь — разбиение отрезка [а, Ь], с = у0 < у1 < х2 < < ... < ут = в, — разбиение отрезка [с,в\. Положим
¡Т(х,у) = тф(х,у) + (1 - т)p(x,y), т Е [0, 1].
Рис. 1. Пример элементарной области О
Разобьем отрезок [0,1] точками 0 = т0 < Т1 < т2 < ... < тк = 1 ив области П рассмотрим сетку, определяемую системой точек
М]1(хг,Уз) = (хг,Уэ,/тг(х^у^)), г = 0,...,п, ] = 0,...,т, I = 0,...,к.
Ясно, что все точки (х^,у^Е П. Зафиксируем индексы г,],1, где г = 0,...,п — 1, ] = 0,...,га — 1, I = 0,...,к — 1 и рассмотрим ячейку, задаваемую точками А]1, А^^, А^+и, А^+ц, А^+ъ А+щ+ъ А]+п+1, Аг+1^+ц+1. Данная ячейка определяется следующим образом
ПЩ = {(х,У,г) : Хг <х< Хг+и Уз <у< yj+l, /т (х,у) < г < /тг+1 (х,у)}.
В каждой такой подобласти П^г построим триангуляцию, состоящую из шести тетраэдров. Способ выбора этих тетраэдров для параллелепипеда указан в работе [3]. Мы используем ту же схему.
\
1] 1+1
1+1 Рк1 "К
/
ч+11 41
Рис. 2. Схема построения тетраэдров в одной части ячейки Приведем эти тетраэдры, выписав четыре вершины каждого из них. Итак, полагаем
= (А^1, Ащ+ц,Aij+ll+l,Ai+lj+ll), т2з1 = (Aijl+l, Aij+ll+l, Аг+1^+ц+1 ,Ai+lj+ll),
T'iji — {Aiji,Aiß+i,Ai+ij+u,Aij+u+i), T?3l — (Aiji, Ai+ij+u, Ai+iji, Ai+iji+i), Tiji — {Aiji+i, Ai+iji+i, Ai+ij+u+i, Ai+ij+u), Tifl — (Aiji, Ai+ij+u,Aiji+i,Ai+iji+i). Пробегая все допустимые значения индексов i,j,l, получим триангуляцию всего множества точек [Aiji}, г — 0, ...,п, j — 0, ...,т, I — 0,...,к.
Для формулировки основного результата введем следующие обозначения. Полагаем
hi — min (xi+i — Xi), hl — max (xi+i — xi),
x 0<i<n-i fx 0<i<ra-i
hl — ^min !(^j + i — yj), hl — ^maX !(^j + i — yj), 0<]<m-1 " 0<]<m-i
t1 — min (т+i — Tz), t2 — max (т+i — t).
0<l<k-i 0<1<к-Г
Оценка углов будет проводиться одинаково в каждой ячейки П^. Более того, ее достаточно получить только для тетраэдров Tijl,T2jl,T3jl, так как для остальных трех рассуждения будут такими же.
Для начала отметим следующее. Пусть заданы два вектора
vi — (ai, ßi, 1),V2 — (a-2, ß2,1). Тогда синус угла 0 между ними удовлетворяет неравенству
.2 (ai — a2)2 + (ßi — ß2)2 Sin 0 " (1 + a2 + ßi)(1 + a2 + ß2).
Итак, рассмотрим первый тетраэдр Т^. Его грани лежат на плоскостях, нормали которых равны векторам
щ — (1, 0, 0), V2 — (0,1, 0), щ —(Zi+ij+il — Zij+il, Zij+il — Zijl, —Л ,
v Xi+i — Xi Vj+i — yj )
^ _ / zí+1j+1l — Zjj+U+l Zjj+U+l — Zjji — Л
4 V xí+i - xí ' Vj+I - Vi ' )
Обозначим углы между соответствующими векторами 0l2, 0l3, 0l4, 023, 024, 034- Тогда, используя оценку (1), получаем для угла 034 неравенство
sin2 034 >
' Zii + 1l + 1-Zij+1l \ 2 I ( Zjj+1l + 1-Zjjl + 1 N 2
í Zij + 1l + 1-Zjj+1l\ I í Zjj+1l + 1-Zjjl + 1 A
> _y Xi+1-Xi ) У Vj+1-Vj )_ >
( 1 I í Zi+1j+1l-Zij+1l \ 2 I í Zij+1l-Zijl\ j / 1 I í Zi+1j+1l-Zij+1l+1\2 I í Zij+1l+1 -Zijl\ 2 I
^ У Xi+1-Xi ) У Vj+1-Vj ) J y y Xi+1-Xi ) У Vj+1-Vj ) J
(£ )2 M2 + L2
>-T-—-2-Л • (2)
(1 + 2L2^1 + (i + §MX) + (i + §MX) J
Аналогично приходим к нижним оценкам остальных углов
sin2 01 3 > ^ , sin2 014 > ---1--2 , sin2 012 = 1,
1 + L + L)
Sin2 023 > -+-J-2 , Sin2 024 > -у-1-^.
1 + L ^{ч^ + L)
Из углов других тетраэдров отличительными будут углы для граней, находящихся на плоскости, ортогональной вектору v5 = {j^+^^j!., —1, , то есть на диагональной плоскости. Оценка угла 045 между векторами v4 и v5 будет такой
1
2
sin 045 >
й М1 + L) +(Ч М1 + L)
2 •
Остальные углы в оставшихся тетраэдрах оцениваются через те же величины. В итоге приходим к утверждению, которое показывает, что при определенных условиях на разбиения отрезков [а, Ь] и [с, ¿] при измельчении триангуляции двугранные углы построенных тетраэдров будут отделены от нуля некоторым положительным числом.
Теорема 1. Для любого двугранного угла 0 всякого тетраэдра построенной триангуляции выполнено неравенство
mm{ (|У М2 + ^ (|у М2 + ^ 1 + 2L4
sin2 0 > k }
(1 + 2L2)[1 + (IMl + L) 2 + (|M1 + L)
В частности, если n = m = k и Xi = a + i(b — a)/n, yj = с + j(d — c)/m, т = l/k, то
sin2 0
mm{ (&)2 + L2, (Ж)2 + L2,1 + 2b2}
(1 + 2L2)(1 + (^ + L)2 + (^ + L)2) Наконец, для куба имеем sin2 0 > 1/3.
ПРИМЕЧАНИЕ
1 Работа выполнена при финансовой поддержке РФФИ, проект № 15-41-02517-р_по-волжье_а.
СПИСОК ЛИТЕРАТУРЫ
1. Алейников, С. М. Алгоритм генерации сетки в методе граничных элементов для плоских областей / С. М. Алейников, А. А. Седаев // Математическое моделирование. — 1995. — № 7 (7). — С. 81-93.
2. Галанин, М. П. Разработка и реализация алгоритмов трехмерной триангуляции сложных пространственных областей: итерационные методы / М. П. Галанин, И. А. Щеглов // Препринт ИПМ им. М.В. Келдыша РАН, 009. — 2006. — С. 1-32.
3. Галанин, М. П. Разработка и реализация алгоритмов трехмерной триангуляции сложных пространственных областей: прямые методы / М. П. Галанин, И. А. Щеглов // Препринт ИПМ им. М.В. Келдыша РАН, 010. — 2006. — С. 1-32.
4. Клячин, А. А. О равномерной сходимости кусочно-линейных решений уравнения минимальной поверхности / А. А. Клячин, М. А. Гацунаев // Уфимский математический журнал. — 2014. — № 6 (3). — С. 3-16.
5. Клячин, А. А. Равномерная триангуляция плоских областей / А. А. Клячин // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2011. — № 15 (2). — С. 43-49.
6. Клячин, В. А. О многомерном аналоге примера Шварца / В. А. Клячин // Изв. РАН. Сер. математическая. — 2012. — № 76:4. — С. 41-48.
7. Клячин, В. А. Триангуляция Делоне многомерных поверхностей и ее аппроксима-ционные свойства / В. А. Клячин, А. А. Широкий // Изв. вузов. Математика. — 2012. — № 1. — С. 31-39.
8. Немировский, Ю. В. Автоматизированная триангуляция многосвязных областей со сгущением и разрежением узлов / Ю. В. Немировский, С. Ф. Пятаев // Вычислительные технологии. — 2000. — № 5 (2). — С. 82-91.
9. Скворцов, А. В. Алгоритмы построения триангуляции с ограничениями / А. В. Скворцов // Вычислительные методы и программирование. — 2002. — № 3. — С. 82-92.
10. Скворцов, А. В. Обзор алгоритмов построения триангуляции Делоне / А. В. Скворцов // Вычислительные методы и программирование. — 2002. — № 3. — С. 14-39.
REFERENCES
1. Aleynikov S.M., Sedaev A.A. Algoritm generatsii setki v metode granichnykh elementov dlya ploskikh oblastey [Algorithm of mesh formation in the boundary element method for plane domains]. Matеmatichеskoе modеlirovaniе [Mathematical Models and Computer Simulations], 1995, no. 7 (7), pp. 81-93.
2. Galanin M.P., Shcheglov I.A. Razrabotka i realizatsiya algoritmov trekhmernoy triangulyatsii slozhnykh prostranstvennykh oblastey: iteratsionnye metody [Development and implementation of algorithms for three-dimensional triangulation of complex spatial domains: iterative methods]. P^print IPM im. M.V. ^ldysha RAN, 009, 2006, pp. 1-32.
3. Galanin M.P., Shcheglov I.A. Razrabotka i realizatsiya algoritmov trekhmernoy triangulyatsii slozhnykh prostranstvennykh oblastey: pryamye metody [Development and implementation of algorithms for three-dimensional triangulation of complex spatial domains: direct methods]. P^print IPM im. M.V. ^ldysha RAN, 010, 2006, pp. 1-32.
4. Klyachin A.A., Gatsunaev M.A. O ravnomernoy skhodimosti kusochno-lineynykh resheniy uravneniya minimalnoy poverkhnosti [On uniform convergence of piecewise linear solutions of the minimal surface equation]. Ufimskiy matеmatichеskiy zhurnal [Ufa Mathematical Journal], 2014, no. 6 (3), pp. 3-16.
5. Klyachin A.A. Ravnomernaya triangulyatsiya ploskikh oblastey [Uniform triangulation of planar domains]. Vеstnik Volgogradskogo gosudars^nnogo univеrsitеta. Sеriya 1, Matеmatika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics],
2011, no. 15 (2), pp. 43-49.
6. Klyachin V.A. O mnogomernom analoge primera Shvartsa [A multidimensional analogue of Schwartz example]. Izv. RAN. Sеr matеmatichеskaya [Izvestiya: Mathematics],
2012, no. 76:4, pp. 41-48.
7. Klyachin V.A., Shirokiy A.A. Triangulyatsiya Delone mnogomernykh poverkhnostey i ee approksimatsionnye svoystva [Delaunay triangulation of multidimensional surfaces and its approximation properties]. Izv. vuzov. Matеmatika [Russian Mathematics], 2012, no. 1, pp. 31-39.
8. Nemirovskiy Yu.V., Pyataev S.F. Avtomatizirovannaya triangulyatsiya mnogosvyaznykh oblastey so sgushcheniem i razrezheniem uzlov [Automated triangulation of multiply connected domains with concentration and rarefaction of nodes]. Vychislitеlnyе tеkhnologii [Computational Technologies], 2000, no. 5 (2), pp. 82-91.
9. Skvortsov A.V. Algoritmy postroeniya triangulyatsii s ogranicheniyami [Algorithms for constructing a triangulation with restrictions]. Vychislitеlnyе mеtody i programmirovaniе, 2002, no. 3, pp. 82-92.
10. Skvortsov A.V. Obzor algoritmov postroeniya triangulyatsii Delone [Review of algorithms for constructing Delaunay triangulation]. Vychislitеlnyе mеtody i programmirovaniе, 2002, no. 3, pp. 14-39.
TRIANGULATION OF SPATIAL ELEMENTARY DOMAINS
Alеksеy Aleksandrovich Klyachin
Doctor of Physical and Mathematical Sciences, Head of Department of Mathematical Analysis and Function Theory, Volgograd State University klyachin-aa@yandex.ru, matf@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Anzhеlika Yuryevna Bеlеnikina
Student, Institute of Mathematics and Information Technologies,
Volgograd State University
matf@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Abstract. We consider a domain Q c R3 that has the form
Q = {(x,y,z) : a < x < b, c < y < d, <p(x, y) < z < ^>(x, y)} ,
where <p(x,y) and 4>(x,y) are given functions in rectangle [a,b] x [c,d] which satisfy Lipschitz condition. Let a = x0 < x1 < x2 < ... < xn = b be a partition of the segment [a,b] and c = y0 < yi < x2 < ... < ym = d be a partition of the segment [c,d]. We put
fT{x,y) = x^{x,y) + (1 - x)p{x,y), T G [0, 1].
We divide the segment [0,1] by points 0 = t0 < t1 < t2 < ... < Tk = 1 and consider the grid in the domain Q defined points
Aijl(xi,yj,Ziji) = (xi,yj,/t,(Xi,yj)), i = 0,...,n, j = 0,...,m, I = 0,...,k.
In this paper we built a triangulation of the region Q of nodes A^i such that a decrease in the fineness of the partition, and under certain conditions, the dihedral angles are separated from zero to some positive constant.
Key words: triangulation, tetrahedron, dihedral angle, elementary domain, partition of domain, Lipschitz condition.