Потенциальное обтекание эллиптического параболоида

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УЧЕНЫЕ ЗАПИСКИ Ц, А Г И
Том XVIII
19 87
№ 3
УДК 533.6. 011. 32
ПОТЕНЦИАЛЬНОЕ ОБТЕКАНИЕ ЭЛЛИПТИЧЕСКОГО ПАРАБОЛОИДА
И. В. Петухов
Рассмотрена замкнутая параболоидальная система координат, в которую включены предельные координатные поверхности: параболические пластины и пластины с параболическим вырезом. Контуры этих поверхностей являются геометрическим местом точек закругления координатных параболоидов. Описано течение, индуцированное параболическим источником.
Параболоидальная система использована при построении точного решения задачи потенциального обтекания несжимаемой жидкостью эллиптического параболоида и параболической пластины, произвольно ориентированных к набегающему потоку. Рассмотрена общая картина течения. На каждом из координатных параболоидов линии уровня, ортогональные линиям тока, образуются сечениями параболоида плоскостями, одна из которых в общем случае касается параболоида в критической точке поля касательной составляющей скорости, а при поперечном обтекании проходит через его ось.
Точные аналитические решения пространственных задач гидродинамики используются при апробировании численных методов, а также при разработке программ численного расчета пространственного пограничного слоя. Большую роль в таком применении играет известное решение задачи потенциального обтекания трехосного эллипсоида произвольной ориентации. Соответствующее решение для эллиптического параболоида представляет интерес как предельный случай решения для сильно вытянутого эллипсоида.
Известно решение осесимметричной задачи обтекания параболоида вращения, полученное методом конформных отображений [1]. Для построения решения задачи потенциального обтекания параболоида произвольной формы и ориентации естественно использовать парабо-лоидальную систему координат. Ниже решение задачи получено на основе уравнения Лапласа в параболоидальных координатах. При определении потенциала продольного обтекания использован параболический источник. Решение выражается в элементарных функциях.
1. Параболоидальная система координат. Параболоидальные координаты Я, |я, V определяются как корни уравнения
(1. 1)
— оо& lt-Х & lt- Ех & lt- (а & lt-?у & lt- V & lt--{-оо- Ех& lt-Еу с& gt-0, (1. 2)
где х, у, г — прямоугольные координаты- Ех, Еу, с — константы.
Рассматриваемая система является замкнутой, поскольку в (1. 2) включены знаки равенства. Знакам равенства отвечают предельные поверхности параболоидальной системы. Система имеет два фокуса 2= А" /у на оси г:
= Л = -|г& gt- Л-Л =1 й = УЁ^Жх (1. 3)
и характеризуется параметром d/c. При заданных d п с величина Ех выбирается произвольно, Ey=Ex+d2. При х = Я, v уравнение (1. 1) описывает семейства эллиптических параболоидов противоположной ориентации, а при х = ц, — семейство гиперболических параболоидов. Сечения семейств к = Х, ц, v плоскостями х=0 и у = 0 (главные сечения) образуют семейства софокусных парабол с фокусами z=fy и
z=fx-
Основные соотношения для параболоидальной системы имеют вид
X -{- р. -j- v — Ех-- Еу-- 2 cz ,
Xjj, _j. JJ. V + vX = ЕхЕу — 2с (Ех + Ey) z — с2(хъ + у2), Xpv = - 2сЕх Е z — c3 (Е х2 + Ех у2) —
(1. 4)
{Ех-)(Ех-?.)(Ех-ч) 9 (Еу-)(Е,-р)(Еу-1)
Х '- сЦЕх-Еу) ' У ~ c*(Ev-Ex) 5
-1 = (?г~х) (Еу-Х) = (Ех — IX) (Еу — ri
С2 (ц _ X) (V _ X) ' С2 (v — ?х) (X — (J.) '
h (Ех — у) (Еу — ¦-)
С2 (X — v) (?л — v) '
где кх = (2Щ~1, Нх — параметры Ламэ (Ах& gt-0),
Здесь и ниже — орт координатной линии g, направленный в сторону возрастания g. Параметры А* удовлетворяют соотношениям
Н1, Я, ^ п Ах., ^ А2
= 0, — *• ^ ±-^-=0,
Ey — l Ev
hl+ h -f h = c2
Для векторов кривизны k* = (д*°1д%)1Нх координатных линий % имеем
Ла сг° hz х° k% =--------1- ----- - x, о, x = {X, ja, v}.
x — a x — т
Здесь И ниже {X, (X, v] - циклические перестановки ИЗ X, [А, V. Координатные линии х являются линиями кривизны для поверхностей (1. 1) (теорема Дюпена для триортогональных поверхностей [2]), т. е. орты образуют главные направления, а составляющие векторов kA совпадают с главными кривизнами для этих поверхностей. Следовательно, имеем
h h kx, a-z--------------- «G, X -- {X, [A, v) ,
О — X T — X
где kx,» И x — кривизны нормальных сечений поверхности х = = const, касательных к а° и т°, в проекции на нормаль х°. В точках закругления'-(омбилических точках), — и = получим
X = ?х = Ех, [а = v = Еу. (1. 8)
Кривые (1. 8) описывают геометрическое место точек закругления эллиптических параболоидов x = v и х = Я. Из (1. 3) — (1. 5) следует, что к этим кривым принадлежат и фокусы параболоидальной системы
Х = р = Ех, ч = Еу — у. = v = Er X = Ех. (1. 9)
2. Предельные поверхности и особые кривые параболоидальной
системы. Предельные поверхности описываются уравнениями
~к=Ех р = Ех, Еу- у = Еу.
Обозначим
= -¦/* - s = x, у, (2. 1)
Es — х
где в соответствии с (1. 2)
dxs = -j-l при Es — /& gt--0- & amp-Xf = - 1 при Es — х& lt-0. (2. 2)
Тогда уравнение (1. 1) при % = ES можно записать в виде
+^?±?k = _& amp-xiS*2- ifS = {x, у}, (2. 3)
Et — Es с*
где s* - параметр уравнения предельной поверхности.
Из (1. 5) следует
s"2=a_bue. (Es-«)(Bs--). X)0jX={X)[i)V}5 (2. 4)
сЗ (Es — Et)
т. е. s* - известная функция координат а, х поверхности % - Es. Каждая
из предельных поверхностей расположена в одной из главных плоскос-
тей х, у = 0 и является двусторонней, т. е. состоит из верхней и нижней стороны поверхности. Обе стороны имеют общую границу (контур предельной поверхности). Выражение
в-О [5*] - 5 = Х, у
ниже означает, что
й= + 0 при в*& gt-0 (верхняя сторона поверхности),
5 = --0 при 5*& lt-0 (нижняя сторона поверхности),
5 = 0 при 5* = 0 (контур предельной поверхности).
Перейдем к более детальному описанию предельных поверхностей. Семейство эллиптических параболоидов % = % меняется от поверхностей, приближающихся к плоскостям 2сг = -X при Я,-«-°о, до параболической пластины в плоскости х = 0
X = Ех, х = О [х*],
2 сг + Ех «с2
= - а: *2, 0л* 2& lt-оо. (2. 5)
Семейство гиперболических параболоидов х = ^ меняется от пластины с параболическим вырезом в плоскости х = 0
Р — х — О[**], Ех + 2, 0& lt-**2 & lt- оо (2. 6)
до пластины с параболическим вырезом в плоскости у = 0
Xі 2 сг + Еу
Р = ?"у = 0[у*],.
Ех Еу
+
с2
2, 0& lt-У*
(2. 7)
Семейство эллиптических гиперболоидов к=у меняется от параболической пластины в плоскости у = 0
V = Еу, у = 0[у*,
%С2 |
Л----'-?г-^ = У*2& gt- 0& lt-у**<-оо (2. 8)
до поверхностей, приближающихся к плоскостям 2сг =-V при ^-& gt- + оо. Общими границами областей (2. 5), (2. 6) и (2. 7), (2. 8) являются особые кривые (1. 8) параболоидальной системы
X = (а = Ех, х — 0,
— у = Еу, у = 0,
Еу Ех
2 сг + Ех -1------------------------------- = 0
+ ¦
С2
2 сг + Еу & quot-с2
= 0.
(2. 9)
Эти кривые представляют собой геометрическое место точек закругления на эллиптических параболоидах к = у, X. Вершина любой из парабол (2. 9) совпадает с фокусом другой.
Выражения (1. 4) — (1. 7) на предельных поверхностях (2. 1) -(2. 3) за вычетом их контуров имеют вид
х = Е5, а х = - 2сг, ах = - 2с Еі г — с2 і2,
8 = ±о, ?V* = -(?,-о) (Е,-т),
Л* = 0,
_Е, -а Л?___________- Т
С2 х — а ' с3 а — х
(2. 10)
ЗС°=-|-4
¦Ь Г
т
К
г г
Е, — т
(2. 11)
где /, 5= {л, у}- о, т, х = {X, ц, V}, причем х = X, ?а при э — Х] х = р., V при 5 = г/. Верхний или нижний знак «+ «, «-» отвечает верхней или нижней стороне поверхности К=Еа.
На особых кривых (2. 9)
Е$у 5
О,
Ег
+ **. ±Ё? = 0- 5, Ь={х, у) (2. 12)
получается однопараметрическая зависимость для Аа, Ат и а°, т0. Выражения (1. 4) — (1. 7) на этих кривых можно представить в виде
Е3 — 2сг,
5 = 0, с2Р- - {Е (- Е$(Е (- х),
Е5 X
Е (-
х
X
5/
и° г°
— * с И
Е (~Е$
Ас^рсоэх, Ат = р эш у, 0& lt-х<-^/2,
=+ «° ЭШ X — р° С08 X, = + 5° соэ X ~ Р° ЭШ X,
(2. 13)
(2,14)
ю»
+ 5° БШ 2х + р° сов 2х,
(2. 15)
где t, ж = {*, з/}- о, т, х = {X, ?л, V), причем / - V при 5 = л: — х = Х при 5 = & gt-'-- р° - орт нормали к контуру (2. 12) в плоскости контура- щ°-орт, ортогональный к контуру.
Обозначим
т — о, — о — г2 соэ2 х& gt- х'-
Е = г2 вт2/.
Соотношения (2. 14) получаются путем предельного перехода в (1. 6), (1. 7) при Х=СОП81, е-*-0. Этот переход соответствует движению к кон-^УРУ по лучу с ортом о& gt-°. Верхний или нижний знак «+ «, «-» в (2. 14) отвечает расположению в° в полуплоскости в& gt-0 или 5& lt-0 в соответствии со знаком в (2. 15). При х=0& gt- я/2 выражения (2. 13), (2. 14) сводятся к соответствующим выражениям из (2. 10), (2. 11). Однопараметрическая зависимость (2. 14) отражает тот факт, что в точках закругления поверхности любые два взаимно ортогональные касательные направления могут считаться главными.
В фокусах (1. 9) предельных поверхностей (2. 5), (2. 8) имеем из (2. 13)
¦а-«=т *=е» * =
5 = 0, Ь = + 0, 2сг — - Е$,
А* - 0, = +
где і, «= {х, у} х = Х при і = х- к = v при 1 = у- верхний или нижний знак «+ «, «^-» отвечает верхней или нижней стороне предельной поверхности к = Еі. Выражения (2. 14), (2. 15), где можно положить р = с, р° = г°, сохраняют прежний вид.
3. Параболический источник. Уравнение Лапласа в параболоидаль-ной системе координат имеет вид
* 1 (¦w f) ¦+ ¦* ¦i [v~f^ fi (V/lw f) ¦= °'- pk (v — Ji.) УЩ, /v = - (X — V) У- f (?a), p., = (?J. — X) yjQC), /(x) = (Ex — x) (Ey — x) — x = X, ?л, v.
Его решением Ф = Фj (X) является
У/(Ц^ = К? = const,
Фг = - 2k In {УЁ~^1 + У Еу — X). Потенциальное течение со скоростью ui = uiX°, где
1 rf4& gt-i _ 2kc _____
(3. 1)
(3. 2)
и,=
Hi dl у (?А- X) (V — X) '
индуцируется параболическим источником (стоком),распределенным на параболической пластине (2. 5). Учитывая, что Х0 = 4-л:0 при % = ЕХ, х=±0^[см. (2. 11)], а также (2. 4), получим на поверхности источника: и1 = + игХ0, где
2к у2 2сг 4- Ех
111= -=г, Х*2=- * х
У (Еу-Ех)я**'- Еу-Ех с*
Потенциал Фе
Фе = Ф1 при л& gt- + 0, Фе = -Фг при — О
описывает обтекание пластины (2. 6) с параболическим вырезом, в частности, течение жидкости в слое, ограниченном двумя любыми гиперболическими параболоидами семейства [г = сопз1:. Скорость протекания через параболический вырез % = ЕХ равна ие = - И/Лг°.
Аналогичное решение для эллиптического источника рассмотрено
в [1].
4. Потенциальное обтекание эллиптического параболоида и параболической пластины. Ниже, принимая? ж& gt-0, обозначим
а0 = /Я7, Ь0 = УЁ~
Пусть э — решение уравнения Лапласа (3. 1) -такое, что
= т (4. 1)
и пусть некоторая функция Е=Е (Я) удовлетворяет уравнению
= т ¦§ + г (Х)А = 0. (4. 2)
Тогда, как можно убедиться, функция Ф»
Ф, = «Я (Х)
также удовлетворяет уравнению Лапласа. При этом условию непроте-кания дф8/дх = 0 на параболоиде (1. 1) при Я = 0, а0& gt-0 удовлетворяет любое решение уравнения Лапласа следующего вида:
+ И0 =--------------Е0. (4. 3)
& quot-о ёоУ/о
Здесь и ниже с1о=й{0) для любой функции ??(X).
Элементарными решениями уравнения Лапласа являются х, у, г. Из (1. 5) следует, что условию (4. 1) удовлетворяют э = х, у. Пусть а0& gt-0. Используя (4. 1) — (4. 3) при эфх, у и полагая Е=А, В, получим для потенциалов Фж, Фу поперечных обтеканий параболоида
ф* «* 4--4-•*. ф^у + -~у& gt-
& amp-0 Ао
Л (X) — 2а°ь°
В (к)
(4. 4)
В эллипсоидальной системе координат решения вида (4. 3) при з = х, у, г определяют суммарный потенциал обтекания эллипсоида произвольной ориентации [1]. В соответствии с этим выражения (4. 4) могут быть получены предельным переходом в решении для эллипсоида. Однако такую процедуру не удается провести для потенциала Ф2 продольного обтекания.
При определении Ф2 используем решение (3. 2) для параболического источника
а0Ь0 У а2 — X + VЬ%-
Ф, = г& lt- = -1п-----------0-^то--------- • (4. 5)
Нетрудно проверить, что дФг/дХ = 0 при Х = 0. Можно убедиться, что частные производные по х, у, г от вторых слагаемых в выражениях для Фж, Фу, Фг стремятся к нулю при X-*--оо вдоль любой координатной
ЛИНИИ Х (|Л, v = const).
Рассматривая внешнюю задачу обтекания эллиптического параболоида
х2 у2 2г
-?г + тг Н-----==0, ?& gt-0 & gt- а0 & gt- 0, с & gt- 0 (4. 6)
ао Ьо с
или параболической пластины (2. 5)
, х = 0[л: *], ~х*2, 0& lt-х*2 & lt-оо, 6ц& gt-а0 = 0, с& gt-0 (4. 7)
Ьо с
со скоростью набегающего потока
— их1их== чхх° +. Т>-3,0 + Т22Г°. Иоо=|"оо|,. (4. 8)
выделим семейство параболоидов (1. 1) при х = А& lt-0×2 у2 2 сг + X
+ -?5-= 0, — оо& lt-Х<-0,
и = йо — X, Ь = V Ьо — X
Решение (4. 4), (4. 5) задачи (4. 6), (4. 8) представим в виде
— «/"00 = grad Ф,
. I х У а0Ьа, а + Ь
ф = /1(1,т + т, т)+1г (2−2(.),
а0 + Ь0
к = а + Ь0 -^ + *9 =Ь + а0
а + Ь
ао + Ьй, а --Ъ
(4. 9)
(4. 10)
где и — вектор скорости.
На параболоиде (4. 6) имеем: гг0=0, к0 = а0 -}- Ь0. Перепишем (4. 9) в других обозначениях:
х*2 + у*2 +
Ус2,
X
а
У
У = ~Г
г* =¦
(4−11)
и в соответствующей форме представим решение (4. 10)
Ф = к{'-(хх* + у у*) + с^(-г* - 2*), 2* = 2-/с. (4. 12)
Из (1. 5) следует, что х*, у* не меняются вдоль координатной линии к:.
. Решение в форме (4. 12) включает решение задачи (4. 7), (4. 8) для параболической пластины. В этом случае имеем: 2 г = 0, !г = Ь. На поверхности пластины при заданных у, г (или у*, г*) и знаке стороны поверхности параметр х* определяется из (4. 7)
х* = + у -у2 — 2г* при л: = + 0.
Если Ф рассматривается как функция от х, у, г, то К определяется как наименьший корень уравнения (4. 9). Очевидно, что все выражения (4. 8) — (4. 12) справедливы и для параболоида вращения (4. 6) при Ьо=а0.
5. Общая картина течения. Ниже потенциал Ф будем рассматривать в функциональной зависимости Ф (х, у, г, Я), где Я- параметр, т. е. Ф будем рассматривать на параболоидах семейства (4. 9). Выделим параболоид (4. 9) при фиксированном Я& lt-0 и рассмотрим векторное поле касательной составляющей скорости ю^/"оо = -gradFФ, где? га^ - градиент на поверхности У7 параболоида. Отбрасывая постоянное сла-
гаемое будем иметь из (4. 10), (4. 12)
— иР1ик = ёгас!/? Ф/т, ФР = Ф + Тг гь фр = (ц*, г) = (д, г*), йФР = (д*, сіг) = (д, сіг*),
4* = - Т^° + 4- ТуУ + ї,*0, 9,"Ат, д?0 + Ат),^0 + С’г, г0, —
(5. 1)
где г* = х*лг° +_У*3)° = 2*2Г°, дифференциалы взяты вдоль Л
Пусть ухт^О. Тогда имеется критическая точка Р поля а^, в которой с1Фр = 0 в любом направлении на / Из (5. 1) видно, что эта точка определяется вектором q*, ортогональным к Б в Р. Сравнивая д* с вектором нормали от* к параболоиду
хх у V г°
т* --?г +• -ТГ + -, аJ о* '- с
запишем Фр' в виде
Фр = С тг (от*, г) = с т2 (-4 X* + * + г*), (5. 2)
где индексом р обозначены значения величин в точке Р,
* Ь 1х * А
X = - -, у=-Жш (5. 3)
Значение определяется получим
из (4. 11). Из условия
2,2,2 л 7* + ь + Тг = 1
2
7*
1 +Ур
*2
,. *2
= 1
(5. 4)
Выражения (5. 3) определяют положение критической точки в зависимости от направления набегающего потока (4. 8). При заданном положении Р направление в» определится из (5. 3), (5. 4). Точки Р полей ир составляют геометрическое место точек (кривую Ьр), в которых линии тока поля и ортогональны параболоидам (4. 9). В общем случае кривая Ьр не совпадает с присоединенной линией тока поля и. Выраже-. ния (5. 1) -(5. 4) [в (5. 1, (5. 2)-вторые выражения для Ф^] справедливы и на поверхности параболической пластины.
На любом параболоиде семейства (4. 9) картина касательного течения определяется положением критической точки Р (направлением вектора q* при уг = 0). При поперечном обтекании у2 = 0 критические точки Р полей ир отсутствуют. Из (5. 1), (5. 2) непосредственно следует, что на указанном параболоиде:
а) с точностью до постоянных множителя и слагаемого потенциал Ф равен проекции радиуса-вектора на направление нормали в точке Р (на направление вектора д* при у2 = 0) —
б) линии уровня Ф = const, ортогональные линиям тока, образуются сечениями параболоида плоскостями, параллельными касательной плоскости в точке Р (плоскостями, ортогональными вектору q* при
Yz = 0) —
в) картина касательного течения на параболоиде Я& lt-0 совпадает с картиной течения на непроницаемом параболоиде той же формы и с тем же положением точки Р (с тем же направлением вектора q* при
Yz = 0).
В а), б) исключается случай пластины, а = % = 0.
Полученные результаты могут найти применение при апробировании численных методов решения пространственных задач гидродинамики и пограничного слоя.
ЛИТЕРАТУРА
1. Милн-Томсон JI. М. Теоретическая гидродинамика. — М.: Мир,
1964.
2. Рашевский П. К. Курс дифференциальной геометрии. — М.: Гостехиздат, 1956.
Рукопись поступила 28/XII 1985 г.

ПоказатьСвернуть
Заполнить форму текущей работой