ytS:;—_ Изв. Сарат. ун-та. Нов. сер. 2013. Т. 13. Сер. Математика. Механика. Информатика, вып. 1, ч. 2
УДК 517.51
СОХРАНЕНИЕ ОРИЕНТАЦИИ СИМПЛЕКСА ПРИ КВАЗИИЗОМЕТРИЧНОМ ОТОБРАЖЕНИИ
А. В. Болучевская
Волгоградский государственный университет E-mail: a.v.boluch@gmail.com
Статья посвящена проблеме сохранения ориентации симплекса при квазиизометричном отображении в . Данная проблема возникает в задачах построения расчетных сеток с помощью отображений различных классов. Найдены условия на отображение, обеспечивающие сохранение ориентации.
Ключевые слова: симплекс, ориентация, квазиизометричное отображение, триангуляция.
On the Quasiisometric Mapping Preserving Simplex Orientation A. V. Boluchevskaya
The paper concerns simplex orientation preserving under the quasiisometric mapping. This problem arises from the problem of mesh generation using different kinds of mappings. We find the conditions for the quasiisometric mapping to be simplex orientation preserving.
Key words: simplex, orientation, quasiisometric mapping, triangulation.
ПОСТАНОВКА ЗАДАЧИ
При численном решении различных задач существует множество способов построения расчетных сеток. Один из таких способов — отображение некоторой стандартной сетки на заданную область. Однако при этом важно контролировать искажение исходных ячеек, чтобы не допустить их вырождения при отображении. Этот вопрос исследуется, например, в работах [1-7].
Для частного случая нерегулярных сеток — триангуляций — в [5-7] показано, что одним из наиболее важных условий сохранения триангуляции является сохранение ориентации каждого симплекса. В [6, 7] найдены условия на отображения некоторых классов, обеспечивающие сохранение ориентации треугольника на плоскости. Кроме того, в [7] поставлена задача нахождения аналогичных условий для сохранения ориентации симплексов в пространстве Мп. Решению этой задачи для квазиизометричных отображений посвящена данная статья.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ
Пусть задана область Б с Мп и ориентированный симплекс 5 (<= Б, образованный точками Ро,Р2,...Рп е Б.
Для набора векторов , £2,... £п введем обозначение
det(£i ,6 ,...£n )
£11 £21
^12 £22
£1n £2n
£n1 £n2
£nn
где £ = (£¿1 ,£г2,... ,£т), г = 1,... ,п.
Будем говорить, что симплекс 5 имеет положительную ориентацию, если точки Р0,Р1, Р2
.,Pn
заданы так, что
det(P1 - Po, P2 - Po,..., Pn - Po) > 0,
и отрицательную ориентацию, когда
ёв1;(Р1 - Ро,Р2 - Ро,... ,Рп - Ро) < 0.
Пусть отображение /: Б ^ Б*, Б* с Мп дифференцируемо п. в. и квазиизометрично [7], т.е. существуют такие положительные Ь, 1, I < Ь, что для любых х, у е Б выполнено
- у|<|/(х) - /(у)|< Ь|х - у|.
Пусть / дифференцируемо в точке Р0. Через 3$ (Р0) обозначим якобиан / в этой точке, через и(£) — модуль непрерывности [8] дифференциала /, йт — длина максимальной стороны симплекса 1 т
5 и д(т) = — Г и(£) М, т > 0. Тогда выполнена следующая теорема. т о
© Болучевская А. В., 2013
А В. Болучевская. Сохранение ориентации симплекса при тазтзометртном отображения
Теорема. Если якобиан отображения / в точке р0 симплекса Б удовлетворяет неравенству
п
•™ яА¥) - ^ -
где V — объем симплекса Б, то симплекс Б' с вершинами Р0 = /(Р0), Р{ = /(Р1), • • •, РП = /(Рп) имеет ту же ориентацию, что и симплекс Б.
Доказательство. Предположим, что симплекс Б имеет положительную ориентацию. Величина ёе1;(Р1 — Р0, Р2 — Р0, • • •, Рп — Р0) равна увеличенному в п! раз ориентированному объему симплекса, натянутого на векторы Р1 — Р0, Р2 — Р0, • • • ,Рп — Р0. Следовательно, если V — объем симплекса Б, то в силу его положительной ориентации получим:
п! ■ V = det(Pl — Р0, Р2 — Р0, • • •, Рп — Р0)•
Пусть V' обозначает ориентированный объем симплекса Б'. Необходимо показать, что V' > 0. Для вершин этого сиплекса имеем:
Рк — Р0 = /(Рк) — /(Р0) = ¿Ро/(Рк — Р0) + Н(Рк), к = 1, 2, • • •, п,
где Н(Рк) = /(Рк) — /(Р0) — ¿Ро/(Рк — Р0), ¿Ро/ — дифференциал / в точке Р0. Тогда
V' = п^еЧР' — Р0 ,Р2 — Р0, • • • ,РП — Р0) =
= П^Ц^ / (Р1 — Р0),^Ро / (Р2 — Р0), • • • ,^Ро / (Рп — Р0)) + Пу det (Н (Р1), • • • ,Н (Рп)) + 1 п
+ - £ det(¿Ро/(Р1 — Р0), • • •, ¿Ро/(Рг1-1 — Р0), Н(Рп), П! ¿1=1
¿Ро /(Рг1 +1 — Р0),•••, ¿Ро /(Рп — Р0)) +
1 п —1 п
+щЕ Е deЧdPо / (Р1 — Р0 ),•••,¿Ро / (Р<1—1 — Р0),Н (Р<1),
¿1 =1 ¿2=»1 +1
¿Ро /(Рг1 +1 — Р0),•••, ¿Ро /(Рг2 — 1 — Р0), Н(Р<2 ), ¿Ро /^ + 1 — Р0),•••,
.. п—2 п —1 п
¿Ро/(Рп — Р0)) + ^Е Е Е deЧ¿Ро/(Р1 — Р0),•••, ¿Ро/(Р*1—1 — Р0),Н(Р<1),
¿1 =1 ¿2=»1 +1 ¿3=^2 +1
¿Ро / (P¿1 +1 — Р0 ),•••, ¿Ро / ^¿2 — 1 — Р0), Н ^¿2 ), ¿Ро /№2 + 1 — Р0),•••, ¿Ро /(P¿3 — 1 — Р0), Н(P¿3 ), ¿Ро /(P¿3 + 1 — Р0),•••,
1 2 п—1 п
¿Ро/(Рп — Р0)) + ■ ■ ■ + п Е ■ ■ ■ Е Е det(¿Ро/(Р1 — Р0),•••, ¿Ро/(P¿1—1 — Р0),
¿1=1 ¿„^=¿^-3 + 1 ¿п-1 =¿^-2 +1
Н(P¿1 ), • • • , ¿Ро /^ —1 — Р0), Н^ ), • • • , ¿Ро /(Рп — Р0))
Отсюда в силу неравенства Адамара [9] для определителей получим:
1 1 п V' > V.//(Р0) — п|Н(Р1)| ■ • • • ■ |Н(Рп)| — - £ |^Ро/(Р1 — Р0)| ■ • • • ■ |^Ро/(P¿l—1 — Р0)| ■ |Н(P¿l)|х
¿1 =1
.. п —1 п
Х|^о/(P¿l +1 — Р0)| ■ • • • ■ |^Ро/(Рп — Р0)| — ^ Е 1^Ро/(Р1 — Р0)1 ■ • • • ■ 1^Ро/^¿1 —1 — Р0)|Х
¿1 =1 ¿2 =¿1 +1
Х|Н^¿1 )| ■ |^Ро/^¿1 +1 — Р0)| ■ • • • ■ Мро/^¿2 — 1 — Р0)| ■ |Н(^^¿2)| ■ МРо/(^^¿2 + 1 — Р0)| Х • • •
.. 2 п —1 п
х|^Ро/(Рп — Р0 )| — ••• — ^^ ••• Е Е 1^Ро / (Р1 — Р0 ^ ••• ■ |Н (Рб1 )| Х •••
¿1 =1 ¿п-2 =¿^-3 +1 ¿^-^¿п-2+1
Х|Н ^ ••• ■ |^Ро / (Рп — Р0 )|-Математика 21
Изв. Сарат. ун-та. Нов. сер. 2013. Т. 13. Сер. Математика. Механика. Информатика, вып. 1, ч. 2
Из квазиизометричности / имеем:
Ир0/(Р - Ро)| =
Нш
е^0
/(Ро + £(Р - Ро)) - /(Ро)
< Ь|Р - Ро г = 1, 2,..., п.
Кроме того, для всех г = 1, 2,..., п выполнено [6]
|Н(Р)| _ |/(РО - /(Ро) - ¿ро/(Р - Ро)| < 1
|Р<—Ро |
|Р - Ро|
|Рг - Ро|
|Рг - Ро |
Ввиду монотонности функции <^(£):
|Н(Рг) |
|Рг - Ро |
Используя полученные оценки, имеем:
..и ..и и
V' > (Ро) - -дп(^т) П |Рг - Ро| - ) П |Рг - Ро | -
— • 1 — • • 1 -1
г=1 гх =1 г=1
и-1 и и
— •
Рп-У(^Щ |Рг - Ро|
гх =1 ¿2 =¿1 +1 2 и — 1
г=1
- ■■■ £ £ ьди—)П |Рг - Ро|
о| =
г1 =1 гта-2=гИ-з + 1 ¿п-1 =гИ-2 + 1 и
г=1
г=1
— , ти — 2,.2,1 М- - 1)
V,// (Ро ) - - П |Рг - Ро| ( —1 ^¿т) Ц + Ри —2д2 (¿т )
2!
+
+ьи—3д3 (¿т)—(— 1)(— 2) + ■ ■ ■ + Ьди—)—(— 1) •••2 + ди (¿т) ) =
3!
(— - 1)!
(Ро) - -7 П |Рг - Ро | £ .!( —• ■)! Ри —г ^(¿т).
—• г!(— - г)!
г=1
г=1
Отсюда,
V' > (Ро) - £((1 + ^
ии
1 П |Рг - Ро |.
г=1
Из условия теоремы получаем, что V' > 0. □
V
Замечание. Если величина —- — отграничена от нуля, 7/ (Ро) > 0, то при достаточно
малых выполнено
П |Рг - Ро|
г=1
(Ро) - Ш + ^
ии
1 П |Рг - Ро| > 0,
г=1
поскольку 1 +
ь
- 1 ^ 0 при т ^ 0.
Следствие. Если для отображения / выполнено неравенство
V
П |Рг - Ро|
г=1
> М ь
> —• VI
1+д^') -1),
то симплекс Б' имеет ту же ориентацию, что и симплекс Б. Доказательство. Заметим, что
(Ро) (Ро )| = (7/ (Ро ))2 = Л1 (ро) ■ ... ■ Ли (ро) > Аи (ро)
£
о
и
и
и
и
и
и
и
Научный отдел
а В. Болучевская. Сохранение ориентации симплекса при тазтзометртном отображении
где Ai(po), • • •, An(po) — собственные числа оператора JM(Po)JM(Po), JM(Po) — матрица Якоби отображения f в точке po, A(po) = min (Aj(po)}.
Оценим A(po). Выберем x E S и построим ориентированную спрямляемую кривую y(s) с параметром s, 0 < s < M, где M — длина кривой и po = y(0),x = y(M). Причем касательная к кривой в каждой точке коллинеарна собственному вектору с минимальным собственным значением оператора
JT (Y (s))Jm (Y (s)).
В силу квазиизометричности f имеем:
M
M
1/2
1|x - po| < |f (x) - f (po |d7(s)/(7 '(s))| ds < /ds
1/2
|dY(s)f(Y'(s))| ds <
M
1/2
1/2
< M1/2 / J (y(s))JM(y(s))Y'(s), Y'(s))ds = M1/2 / Л(7(s))|y'(s)|2ds
где Л(Y(s)) — минимальное из собственных чисел оператора JM(y(s))JM(y(s))-Отсюда,
/ M \ 1/2
l MM(s)) ds
| x _ Po l
Если x ^ p0, то -—m— ^ 1, и l < Л1/2(р0). Тогда (Jf (P0))2 > 12n и из теоремы
V' > VP _ — 1 + r n! V V L
g(dm)
1 П |Pi _ Pol
i=1
Отсюда получаем требуемое.
Работа выполнена при финансовой поддержке РФФИ (проект 11-01-97021-р_поволжье_а). Библиографический список
□
1. Азаренок Б. Н. О построении подвижных адаптивных пространственных сеток. М. : Вычислительный центр им. А. А. Дородницына РАН, 2007. 50 с. [Azarenok B. N. On generation of dynamically adaptive spatial grids. Moscow, 2007. 50 p.]
2. Миклюков В. М. Геометрический анализ. Дифференциальные формы, почти-решения, почти-квазиконформные отображения. Волгоград : Изд-во ВолГУ, 2007. 532 с. [Miklukov V. M. Geometric analysis. Differential forms, almost-solutions, almost-quasiconformal mappings. Volgograd, 2007. 532 p.]
3. Ушакова О. В. Условия невырожденности трехмерных ячеек. Формула для объема ячеек // Журн. вы-числ. мат. и мат. физ. 2001. Т. 41, № 6. С. 881-894. [Ushakova O. V. Nondegeneracy conditions for three-dimensional cells and a formula for the cell's volume // Comp. Math. and Math. Phys. 2001. Vol. 41, № 6. P. 832-845.]
4. Бобылев Н. А., Иваненко С. А., Казунин А. В. О кусочно-гладких гомеоморфных отображениях ограниченных областей и их приложения к теории сеток // Журн. вычисл. мат. и мат. физ. 2003. Т. 43, № 6. С. 808-817. [Bobylev N. A., Kazunin A. V., Ivanenko S. A. Piecewise smooth homeomorphisms of bounded domains and their applications to the theory of
grids // Comp. Math. and Math. Phys. 2003. Vol. 43, № 6. P. 772-781.]
5. Прохорова М. В. Проблемы гомеоморфизма, возникающие в теории сеток // Тр. ИММ УрО РАН. 2008. Т. 14, № 1. С. 112-129. [Prohorova M. V. Homeomorphism problems arising in the grid theory // Trudy IMM UrO RAN. 2008. Vol. 14, № 1. P. 112-129]
6. Клячин В. А. О гомеоморфизмах, сохраняющих триангуляцию // Записки семинара «Сверхмедленные процессы». 2009. № 4. С. 169-182. [Klyachin V. A. On homeomorphisms preserving triangulation // Zapiski seminara «Sverhmedlennye processy». 2009. № 4. P. 169-182.]
7. Миклюков В. М. Введение в негладкий анализ. Волгоград : Изд-во ВолГУ, Лаб. «Сверхмедленные процессы», 2008. [Miklukov V. M. Introduction to non-smooth analysis. Volgograd, 2008. 424 p.]
8. Кудрявцев Л. Д. Курс математического анализа : в 3 т. Т. 2. М. : Дрофа, 2004. 720 с. [Kudryavcev L. D. Mathematical analysis. Vol. 2. Moscow : Drofa, 2004. 720 p.]
9. Гантмахер Ф. Р. Теория матриц. М. : Наука, Гл. ред. физ.-мат. лит., 1966. 576 c. [Gantmacher F. R. The Theory of Matrices. AMS Chelsea Publishing : Reprinted by Amer. Math. Soc., 2000. 660 p.]
o
o
o
o
n
n
Математика