Файл: Зак М.А. Неклассические проблемы механики сплошных сред.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 26.07.2024
Просмотров: 94
Скачиваний: 0
ный разрыв скорости vB. Далее |
6[Э„] = 2 (p'vn' + p' [v/]) • б [vs'j-f- |
||
+ 2 (p"vn"+p" [vB"]) -6 [vB"] = 0. |
Учитывая, что [Оаз'^—^ [«э'] и |
||
полагая |
[иВ1/]= [г'в2/]= Ь В1//]= [ив2//]= 0, приходим |
к соотношению |
|
{х2в- ( 2 |
+*v) |
Ы |
=+О, *Ц[ай] Ф 0. |
Отсюда получаем скорость распространения поперечной волны рассматриваемой поверхности в направлении скорости разрыва ив и условие «размывания» этой поверхности:
1 __ |
Р ив |
| |
^п1 _i_ |
~l f |
^“111 |
P,p,,U"B |
,, |
2нв У р'р" |
'в — |
Р'+Р" |
“Г |
2 - |
у |
4 |
(Р'+Р")3’ |
п1<" |
Р'+Р" • |
Здесь vni — проекция vn на направление скорости ив. Последнее условие упрощается, если разделяемые области жидкости имеют одинаковые плотности: uni<wB. Меняя ролями вихревую и потен циальную составляющие скорости, получаем аналогичные ре зультаты для поверхности тангенциального разрыва потенциаль ной составляющей скорости:
} ___ |
Р "« П |
, «В1 |
_L |
l / |
и 3В1 |
р ' р " « ап |
/ |
2 « п W |
p |
" |
'п |
р' + р" |
2 |
— |
V |
4 |
р' + р" ’ |
v>il ^ |
р' + |
р" |
’ |
а при равенстве р' = р" последнее условие принимает вид oBi< u n.
Здесь ип— тангенциальный разрыв |
скорости vn, uBi — проекция |
|
vBна направление скорости ип. |
|
|
10. |
О единственности решения уравнений движения идеальн |
|
жидкости. |
Использование принципа |
наименьшего принуждения |
(II—1—4) |
может оказаться эффективным и в том случае, когда |
классические уравнения движения сплошной среды при физиче
ски оправданных граничных условиях хотя и имеют |
решение |
в классе непрерывных функций, но это решение не |
является |
единственным. Примером тому могут служить задачи об обтека нии твердого тела идеальной несжимаемой жидкостью.
Рассмотрим вначале задачу о плоском стационарном безот рывном обтекании твердого тела безграничным потоком идеаль ной (невесомой) несжимаемой жидкости.
Известно (см. [16]), что такая задача при заданной по вели чине и направлению скорости на бесконечности v<*>имеет бесчис ленное множество решений, зависящих от выбора величины цир куляции скорости Г. И только в том случае, когда на задней кромке обтекаемого тела расположена угловая точка, можно, ис ходя из постулата Жуковского — Чаплыгина, определить величи ну циркуляции Г, а следовательно, выделить единственное реше ние задачи. При отсутствии на профиле угловой точки постулат уже не применим и для выделения единственного решения при ходится привлекать дополнительные предположения.
Обратимся к принципу наименьшего принуждения в форме (II—1—4), согласно которому истинное движение жидкости ми нимизирует функционал
100
Ф(т^, r) = j аЧа, |
(III—6—25) |
а |
|
где а — часть плоскости, занятая жидкостью. Тогда искомая, цир куляция Г* определится из условия
Ф ^ , Г*) = inf Ф (■ово, Г).
Отсюда, как частный случай, вытекает постулат Жуковско го— Чаплыгина; действительно, при наличии на обтекаемом те ле угловой точки циркуляция Г* должна принять такое значе ние, чтобы интеграл в (25) не расходился, т. е. ускорение не об
ратилось бы в бесконечность в угловой точке. |
для |
обтека |
|||
Найдем, используя (25), единственное решение |
|||||
ния круга с радиусом R. Согласно [16] для комплексно-сопряжен |
|||||
ной скорости в этом случае имеем выражение |
|
|
|||
ц = К |
I (l - |
4 9 |
+ 2Wz’ z = rel?’ >'2 = ХЛ' + У2 > R2- (III—6—26) |
||
Тогда |
|а |2= |и |2 |
\\dv/dz\2 и, следовательно-, |
|
|
|
|
О2“ |
|
1 |
2Ь |
л |
Ф = J J |
|
||||
| a fd y rd r — 32*з/?4 |
4 |
vОО • |
R о
Из равенства дФ/дГ = 0 следует, что принуждение Ф достига ет экстремума лишь при Г =0, причем из неравенства 52Ф/дГ2>0 при Г = 0 следует, что этот экстремум является минимумом. Сле довательно, задача об обтекании круга имеет единственное реше ние в классе непрерывных функций, причем это решение соответ ствует безциркуляционному обтеканию.
Докажем существование единственного решения задачи об обтекании произвольного плоского профиля без угловых точек в классе непрерывных функций. Для этого воспользуемся выра жением для комплексно-сопряженной скорости, получаемой из (26) соответствующим конформным преобразованием. Учитывая, что циркуляция Г по-прежнему будет линейно входить в это вы
ражение, получаем Ф = А4Г4+ А3Г3-f А2Г2+ А4Г + А0 > О,
причем коэффициенты А* зависят от v«, и от вида обтекаемого профиля.
Пусть А4=т^ 0. Тогда уравнение <ЗФ/дГ=0 является кубическим относительно Г и, следовательно, имеет по крайней мере один вещественный корень. Этот корень должен соответствовать ми нимуму Ф, так как Ф->-оо при Г-^-оо. Если же уравнение <ЗФ/5Г=0 имеет три вещественных корня, то два из них соответ ствуют минимумам Ф, и из этих корней следует выбрать тот, при котором принуждение Ф наименьшее.
Пусть А4= 0. Тогда и А3 = 0; в противном случае Ф<0 при Г-9— оо, что невозможно. Если А2ф 0, то из равенства дФ/дТ = 0
101
получим Г= —А\/2А2. Этот корень соответствует минимуму Ф, так как Ф->-оо при |Г|-*-оо. Если А2 = 0, то и i4j = 0; в противном случае Ф <0 при Г-»— оо, что невозможно; но тогда принужде ние Ф окажется вообще не зависящим от циркуляции Г. В этом случае вопрос о единственности решения остается открытым.
Таким образом, существование единственного решения дока зано.
Заметим, что решения с безотрывным обтеканием, получае мые при условии divv = 0, сохраняют смысл лишь при р^О , так как в действительности div v^O (см. гл. II, § 3, п. 3). В качестве иллюстрации рассмотрим полуокружность, описываемую урав нением x2+ y2 = R2 (х'<0) и обтекаемую потоком, имеющим на бесконечности скорость Vco, параллельную оси х. Нетрудно прове рить, что комплексно-сопряженная скорость
г , = 1'у ~ | ( 1 ------§ |
" ) ’ |
z — r e 1'р, |
0 |
< г < |
оо, |
|||
удовлетворяет и уравнению |
Эйлера |
(II—3—5), |
и |
уравнению |
||||
div и= 0, и граничным условиям, |
причем обтекание является без |
|||||||
отрывным. Действительно, при \ z \ ^ R |
приходим к |
|
(26), т. е. |
|||||
к полю скоростей, получающемуся при |
обтекании |
круга. При |
||||||
\z \ ^ R образуется область, |
в которой жидкость циркулирует по |
|||||||
замкнутым траекториям (диполь). Однако |
это решение все же |
|||||||
непригодно, так как и->-оо при |
|2 |-ЯЗ и, в соответствии с урав |
|||||||
нением Бернулли, р < 0 при |
|г|->-0для любого конечного, давле |
|||||||
ния жидкости на бесконечности. |
|
|
|
|
|
|
||
Обратимся теперь к случаю отрывных обтеканий. Положим, |
||||||||
что уравнение контура |
обтекаемого |
тела |
имеет |
вид г (s) = |
||||
= r (s + nl), л = 1, 2 ,..., |
где 5 — дуговая |
координата, |
I — длина |
контура. Введем вспомогательный контур г* (s), охватывающий тело так, что |г * |^ |г |. Очевидно, что любое решение задачи об обтекании контура г* (s) является одновременно и решением исходной задачи, но уже в более широком классе функций, допу скающем разрывы поля скоростей вдоль некоторых линий (в дан ном случае вдоль контура г* (s)). Поэтому имеет место оценка ф * ^ ф , где Ф* определяет принуждение в случае отрывного об текания. Варьируя г* (s) и величину тангенциального разрыва скорости вдоль этого контура так, чтобы достигнуть inf ф*, при дем к искомому отрывному обтеканию. В частном случае может, конечно, оказаться, что г* (s) = г (s), и тогда имеет место безот рывное обтекание/
Особый интерес при исследовании отрывного обтекания пред ставляет случай, когда минимизирующим для Ф* оказывается тангенциальный разрыв скорости вдоль г* (s), равный величине самой скорости на этом контуре; при этом между контурами г* (s) и г (s) образуется застойная зона (ц=0), а поэтому в силу уравнения Бернулли величина скорости вдоль контура г* (s) должна быть постоянной. Последнее обстоятельство позволяет
102
находить такого рода отрывные обтекания посредством уравне ния Эйлера (II—3—5) с учетом равенства divv = 0 и граничных условий v„ = 0, |w |= const при г = г* ( 5), где v„ — составляющая скорости жидкости, нормальная к контуру г* (s). Однако и в этом случае единственное решение можно выделить лишь то гда, когда известны заранее точки отрыва (например, при обте кании тела с угловыми точками). В противном случае следует обратиться к принципу наименьшего принуждения и. искать точ ки отрыва г+ из условия
ф (г+. = inf Ф (г+*, г>те),
где г+ * — возможные точки отрыва на контуре г (s).
Непосредственное практическое использование решений с от рывным обтеканием осложнено тем, что, как это следует из пре дыдущего пункта, линии тангенциальных разрывов скорости не устойчивы; следовательно, описанные выше решения будут не устойчивы в классе кусочно-непрерывных функций, т. е. сколь угодно малые случайные возмущения скоростей, перпендикуляр ные к линии разрыва, приведут к возникновению «всюду разрыв ного» поля скоростей. Действительно, в соответствии с (II—3—5) недетерминированные связями движения жидкости со скоростями vr (пульсационные движения) могут появиться или в силу специфических граничных условий, как в случае пе ресечения двух потоков (гл. II, § 3, п. 3), или в силу «всюду раз рывных» начальных условий, которые в рассматриваемом случае имеют место из-за неустойчивости линии тангенциального
разрыва скорости. Расширим класс функций, |
в |
котором |
отыскивается движение жидкости в застойной |
зоне, допу |
|
стив появление в ней недетерминированных связями |
скоростей |
vr (пульсационных скоростей) в соответствии с (I—3—23). В ка честве конкретного примера возьмем пластинку, обтекаемую по током, перпендикулярным к ее плоскости. Точками отрыва в этом случае являются концы пластинки, а застойная зона, в со-' ответствии с [16], простирается от тыльной стороны пластинки до бесконечности, неограниченно расширяясь. Для описания явле ний, возникающих после потери устойчивости границы застойной зюны, можно воспользоваться уравнениями, вытекающими из ми
нимизации функционала (II—1—7) |
(см. гл. II, |
§ 3, п. 3): |
1 |
|
|
р -^ - = — grad р, d iv [v cl; = 0, |
v = v(r, |
£) при ^= О, |
о |
|
|
считая, что р = р (г, t), v= v (г,t,Q , O sS ^^l, т. е. полагая, что состояние каждой точки пространства характеризуется множест вом скоростей, образующих континуум v (£). Начальные значе ния v(£) обязаны своим происхождением случайным факторам, благодаря которым оказалась нарушенной устойчивость границы
103
застойной зоны, и, следовательно, при выборе этих значений дол жны быть привлечены статистические соображения.
11. Движение гибкой нити в идеальной жидкости. Проект руя уравнения движения гибкой нити в виде (II—2—9) с учетом сил инерции на оси естественного трехгранника, получаем
dv°, . |
0 |
- |
ди |
ди |
дТ |
с |
(III—6—27) |
- з г + |
^ з - |
“Н dt |
и I f |
|
|
||
р(~дГ~ + шз^01— |
+-2w3«j + (ри2— ?) S3 = |
^2. |
(Ш—6—28) |
||||
|
p ( ^ f + |
“ |
ш^ ° ‘ ~ |
2шз») = Ft. |
(Ill—6—29) |
||
Здесь p — плотность, отнесенная к единице длины; |
Т — натяже |
||||||
ние; v° — переносная скорость нити вместе с ее |
фиксированной |
формой; и — скорость течения нити вдоль фиксированной формы; со — угловая скорость вращения элемента нити; F — внешние силы, отнесенные к единице длины; Qb Q3 — кручение и кривиз на формы нити.
К этим уравнениям следует добавить кинематические соотно шения совместности перемещений и поворотов в проекциях на те же оси:
dvi° |
T)°3S 3, |
dv°2 __ : -y03Q i _ |
^ |
|
|
= |
c)0 |
|
|
|
|
дф |
|
|
|
|
|
|
dv°3 |
■= — V°2 Q i — |
cd2, |
(III—6—30) |
|
|
dii |
||||
|
|
|
|
|
|
дь>1 |
dQi _ |
UJ2W"I3> |
|
^ |
|
ddi |
dt |
city |
* |
||
|
ди>з |
dQ3 |
Ш2'Х| 1 > |
(III—6—31) |
|
|
d6 |
dt |
|||
|
|
|
|
а также связи между течением и, удлинением нити f, натяжени ем Г и начальной плотностью р0:
ди f = |
Ё1. |
df_ |
Т = ? ( / ) , Ро= Р/, |
/= -§-• |
(Ш—6—32) |
|
J |
dt |
и дЬ ’ |
||||
Нетрудно проверить, |
что система (27) — (32) |
замкнута, если |
||||
задана внешняя сила F. Эту силу разобьем |
на |
составляющие: |
||||
F=F(1)+ F ^+ F (3\ где F<9— «мертвые» силы, например силы веса |
||||||
F(9=—p0gk (k — единичный орт вертикальной |
оси |
декартовой |
системы координат); F® — следящие силы негидродинамическо го происхождения, например пондеромоторные силы; будем по лагать, что F(2) = F<2>(oi°, v2°, t>3°); F<3>— силы гидродинамическо го происхождения.
Если считать жидкость идеальной и искать решение задачи о взаимодействии ее с нитью в классе дифференцируемых функ ций, то единственной силой гидродинамического взаимодействия
104