293
{ 
° 
»
В пределах граничных элементов Г.Ф с узлами 1, } ‚ пронумерованными
против часовой стрелки, скорость записывается аналогично
о9) = { ® (0,0Е )и = М, } 
(6.74)
где о'(0)= 2 о„Ф (О) = о,Ф 
(О)+ оФ® (0), 0 ЕГ.
Координаты \)Ш неизвестны и являются значениями искомой функции ско-
рости в узлах сетки:
Эт = о(ХШ,ХШ)= от Е,
где 0 = {Пп‚п = ], ЫШ}= {П1‚ 0,,‚..., Эц } - множество узлов сетки.
Нормальную производную скорости на каждом граничном элементе
примем постоянной:
дА 
дА
Р (О) = (
3 ) (О),О 
е 10 ,г = 1, М,
П
ое[2А 
] @)- ^°( 
- ^“,ое1
Кинематическую вязкость в пределах каждого конечного элемента тоже
будем считать постоянной:
У(0) = /° (0),0е50,г= 1, М, }
где у®) (О) = у) = сопя(, 
О е Г.
Перейдем к преобразованию уравнения (6.64) к виду, удобному для чис-
ленной реализации. Полагая функцию фФ равной базисным функциям в узлах
сетки:
Ф (О 
) = {:(0)1 = 1, М }.
получаем систему уравнений для области $ ‚ число уравнений в которых опре-
деляется числом узлов в области (М,)
[|(У(О)вга4 (О)вта%(О) - Ф:(О)в1о )4 
=
(6.75)
= [9:(О)(О) 
2° (Ох о‚О е 5 =1, М,
Т 
П
Преобразуем каждое слагаемое уравнения (6.75) в отдельности. Начнем с пер-
вого слагаемого. Так как область $ является объединением конечных элемен-
().
тов 5°: