
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfвектор v = gradcp. Функция ср называется |
потенциалом |
скоростей. |
|
Для потенциального поля уравнение (5.12) |
приобретает |
вид |
|
-^-grad<p + |
grad(-y-) = F — ^grad/> . |
(8.5) |
|
Будем опять иметь |
в виду баротропную жидкость. |
Положим |
р (р)
тогда как нетрудно проверить,
1р grad/?=gradP
и уравнение (8.5) переписывается так:
F = g r a d ( - | j - + - £ + p ) . |
(8.6) |
Отсюда видно, что F есть потенциальный вектор, иными сло вами, безвихревое движение возможно только тогда, когда объем ные силы имеют потенциал, т. е. F = —grad V. Теперь (8.6) приво дится к виду
grad ( - £ - + - £ + ^ + / > ) = 0 , |
(8.7) |
показывающему, что выражение в скобках не зависит от коорди нат и является только функцией времени t
|
dt |
V+P=W(f). |
(8.8) |
||
|
|
|
|
|
|
Здесь |
—произвольная |
функция. Выражение |
(8.8) назы |
||
вается интегралом Коши. |
|
|
|
|
|
Если жидкость несжимаема, то уравнение неразрывности при |
|||||
обретает вид |
|
|
|
|
|
|
дх°- 1 ду2 |
1 |
дг* |
-О, |
(8.9) |
|
|
|
и решение задачи о безвихревом движении сводится к отысканию одной функции ф, удовлетворяющей уравнению Лапласа (8.9), краевым и начальным условиям. Гидродинамическое давление опре делится из (8.8), а произвольная функция х¥ (t) может быть най дена, если задана зависимость р от t в одной точке поля. Если, кроме того, движение установившееся, а объемные силы сводятся к одной силе тяжести, то интеграл Коши принимает вид
z + ^ + - | - = c o n s t , |
(8.10) |
внешне совпадающий с интегралом Бернулли (8.4). Однако между выражениями (8.4) и (8.10) есть существенная разница. Первое из них не предполагает, что движение безвихревое, и в общем слу-
чае оно правильно только для отдельных линий тока. Второе же относится ко всему потоку, который обязательно должен быть без
вихревым. |
|
|
|
Приведем |
теперь без |
доказательства |
некоторые теоремы |
•о вихрях. |
|
|
|
Т е о р е м а |
Л а г р а н ж а. |
Если объемные |
силы имеют потен |
циал и если в какой-то момент в некоторой части идеальной баро-
тропной жидкости нет |
вихрей (rotv = 0), то их не было |
раньше и |
не будет позже в этой |
части жидкости (подчеркнем, что |
речь идет |
об определенной массе жидкости, а не об определенной части пространства).
Введем понятие о вихревой линии и вихревой трубке. Вихре вой линией называется линия тока вихря скорости. Вихревая трубка есть трубка, образованная вихревыми линиями, проведен ными через некоторую замкнутую кривую. Интенсивностью вихре вой трубки называется поток вектора вихря через сечение вихревой трубки. Для вихревых линий и вихревых трубок имеет место сле
дующая |
т е о р е м а |
Г е л ь м г о л ь ц а . |
Если |
объемные |
силы имеют потенциал, а жидкость идеаль |
на и баротропна, то частицы жидкости, образующие в некоторый момент вихревую линию, образуют вихревую линию во все после
дующие |
моменты, |
а |
интенсивность любой вихревой |
трубки |
во все время движения |
остается постоянной. |
|
||
Эту теорему называют также теоремой о сохранении |
вихревых |
|||
линий и |
вихревых |
трубок. Более общая теорема, не связанная |
с идеальностью и баротропностью жидкости и потенциальностью объемных сил, принадлежит А. А. Фридману: необходимое и доста точное условие сохранения векторных линий и интенсивности век
торных трубок вектора Е |
(в нашем |
случае £ = rotu) состоит в том, |
|
что г е л ь м г о л ь ц и а н |
этого |
вектора |
|
h e l m E = - ^ |
- ( E , |
V)v + E d i v v |
должен быть равен нулю.
Из теорем Лагранжа и Гельмгольца вытекает, что при движе нии идеальной баротропной жидкости в поле потенциальных сил вихри не могут ни возникать (если их не было), ни исчезать (если они существовали).
Возникновение вихрей связано с разными причинами, основные из них следующие:
1) небаротропность жидкости. Если в начальный момент изоба рические поверхности (поверхности равного давления) не совпа дают с изостерическими (поверхностями равной плотности), то в жидкости возникают вихри, тем более интенсивные, чем выше гра диенты плотности и давления и чем больше угол между векторами градиентов. Это главная причина образования вихрей в атмосфере и вызываемых ими циркуляционных течений;
4* |
51 |
2) иепотеї-іциальпость действующих на жидкость сил, в первую очередь сил внутреннего трения (вязкости). Мы увидим в следую щих главах, что движение жидкости с внутренним трением вообще ие может быть потенциальным. Непотенциалыюсть вектора объем ных сил F, имеющая место в том случае, когда, кроме силы тяже сти, приходится учитывать еще силу Кориолиса, связанную с вра щением Земли, порождает вихри, даже если внутренним трением
вжидкости можно пренебречь;
3)индуцирование вихрей движущимися в жидкости телами, даже когда жидкость идеальная. Третья причина связана с внеш ней задачей гидромеханики, не рассматриваемой в этой книге.
Подчеркнем следующие обстоятельства, характеризующие иде альную жидкость как модель реальной жидкости.
1. Идеальных жидкостей в природе не существует, решения гид ромеханических задач для идеальной жидкости, вообще говоря, не есть пределы, к которым стремятся решения соответствующих задач для вязких жидкостей при стремлении вязкости к нулю.
2. Поскольку не существует жидкостей без внутреннего трения, не существует и безвихревых движений, которые присущи только идеальным жидкостям. Реально можно говорить только о движе ниях, обладающих слабой завихренностью, т. е. мало отличаю щихся от безвихревых.
Несмотря на эти обстоятельства, гидродинамика идеальной жидкости дает очень хорошее приближение к действительности для многих (хотя и далеко не для всех) действительных волновых II вихревых движений. Там, где скорость течения идеальной жидкости как функция координат делает скачок (у стенок, на поверхности раздела двух жидкостей), в реальной жидкости наблюдается не прерывное, но с очень большим градиентом в достаточно тонком слое изменение соответствующих компонентов скорости. Учитывая сказанное, в первую очередь мы рассмотрим именно движения идеальной жидкости — простейшей из всех жидких сред. Здесь не ставится цель дать полную теорию внутренних задач гидродина мики идеальной жидкости. Мы рассмотрим только простейшие вол новые движения и движения с поперечной циркуляцией, чтобы оце нить возможности гидродинамической и гидравлической идеали зации.
Глава III
Н Е К О Т О Р Ы Е Д В И Ж Е Н И Я И Д Е А Л Ь Н О Й Ж И Д К О С Т И
9. Простейшие волновые движения жидкости постоянной плотности
Чтобы выяснить одну из важнейших причин возникновения волн, рассмотрим действие мгновенных сил на жидкость. Для этого на пишем уравнение Эйлера (5.12) в таком виде
|
- 3 T = F — \ - & * * Р - |
|
(9-І) |
||||
|
Положим |
|
|
|
|
|
|
|
F = F 0 + |
F», |
р=Ро+Р*, |
|
|
|
|
где |
и р%— составляющие |
вектора |
F и скаляра |
р, отличные |
от |
||
нуля только в течение некоторого интервала времени от |
^ = 0 |
до |
|||||
t=x, |
но принимающие в течение |
этого интервала |
очень |
большие |
|||
значения, так что пределы |
|
|
|
|
|
|
|
|
lim \F%dt=l, |
|
lim |
\p^dt—v: |
|
|
|
конечны. Эти пределы называются соответственно импульсом мас совых сил и импульсом сил давления. Интегрируя уравнение (9.1) по t от ^ = 0 до t = x, переходя к пределу при т—>-0 и замечая, что для конечных Fo и ро соответствующие пределы равны нулю, полу чаем:
v - v 0 = I - g r a d ( - ^ ) . |
(9.2> |
Здесь Vo—скорость в момент £ = 0 до начала действия |
мгновен |
ных сил. |
|
Наибольший интерес представляет случай 1=0, когда мгновен |
|
ные силы есть только силы давления. Тогда |
|
v - v 0 = - g r a d ( - ^ - ) . |
(9.3> |
Из этого выражения можно сделать следующие выводы. |
|
1. Если до приложения импульса давлений движение |
жидкости |
было безвихревым, т. е. было vo = grad<p, то после приложения им
пульса оно останется |
безвихревым с потенциалом q> — л/р + Су |
где С — произвольная |
постоянная. |
2. Безвихревое движение может возникнуть из состояния покоя (УО = 0) в результате приложения импульса давления. Но никаким импульсом давления нельзя вызвать вихревого движения покоя щейся жидкости.
3. Безвихревое движение с потенциалом скорости ф может быть полностью остановлено во всей жидкости приложением импульса давления я==рф, по никаким импульсом нельзя уничтожить вихре вое движение.
Таким образом, приложение импульса давлений, например к свободной поверхности покоящейся жидкости, приводит эту жид кость в состояние безвихревого движения с некоторым потенциа лом скорости <р, удовлетворяющим уравнению (8.9)
дх* ^ ду2 Т" дгч — и -
Если р = const, то из уравнения (8.8) получаем
Будем предполагать, что движение совершается в поле силы тя жести без участия других массовых сил. Направим ось z верти кально вверх, тогда V=gz. Если движение достаточно медленное, то членом v2[2 в выражении (9.5) можно пренебречь и записать его так:
Р _ |
d?i |
•gz, |
(9.6) |
|
р |
dt |
|||
|
|
|||
где |
|
|
|
|
Так как уравнение (9.4) |
сохраняет силу, если заменить |
в нем <р |
на фі, то индекс единица при ф в формуле (9.9) |
можно отбросить |
и написать ее в виде |
|
f — Я - - * * - |
<9-7> |
Нормальная составляющая скорости на твердых поверхностях, ограничивающих жидкость, равна нулю, следовательно, на этих по верхностях имеет место следующее краевое условие:
0-8)
где п — нормаль к поверхности. На свободной поверхности жид кости давление после окончания действия импульса постоянно, и его, не нарушая общности, можно считать равным нулю, т. е. на свободной поверхности имеем следующее краевое условие:
•&+gz=0. |
(9.9) |
Пусть уравнение свободной поверхности есть z — \(x, у , t). Крае вое условие (9.9) можно записать так:
г д<? (х, у, z, |
t) |
|
dt |
- Ц ( , . , . о + # = 0 - |
( 9 Л 0 > |
Плоскость Оху совместим |
со свободной поверхностью |
жидкости |
в положении равновесия. Мы ограничимся волнами, высота кото
рых очень мала по сравнению с длиной, |
тогда в |
(9.10) |
можна |
||||||||
брать значение dcp/dt |
при 2 = 0, а не при z = £, т. е. |
|
|
||||||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
ді |
|
|
1 |
(х, |
у, 0, |
t) |
|
|
(9.11> |
|
|
dt |
|
|
g |
|
д& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
но так как |
|
|
dx |
|
|
rfy |
|
dz |
|
|
|
|
|
|
|
|
W — |
|
|
|
|||
|
|
и=—гг, |
' |
- V = —тг , |
—rr, |
|
|
||||
|
|
|
dt |
|
dt ' |
|
dt |
' |
|
|
|
то на свободной |
поверхности |
|
|
|
|
|
|
||||
dl |
dS , |
dS |
dx |
, .55 |
dy |
dt |
. |
d5 |
, „, й |
/ Q 1 Г Л |
®=1Г=-эт+-1Г-^гт-^г--^г=^г+и -йГ+г , ^г- (9Л2>
Так как высота волн мала по сравнению с их длиной, то плос кость, касательная к свободной поверхности, мало отличается or горизонтальной плоскости и, следовательно, дЦдх и дЦду в (9.12) пренебрежимо малы. Значит на свободной поверхности
dz |
|
dt |
|
|
и краевое условие (9.11) можно переписать так: |
|
|||
дер |
1 |
д 2 ? |
(9.13) |
|
~~дТ~ |
g~ ' ~d~W |
|||
|
||||
при z = 0. |
|
|
|
|
Имея в виду волновые движения, |
мы не будем формулировать |
начальных условий, а будем искать периодические решения урав
нения (9.4) |
при краевых условиях |
(9.8) |
и (9.13). Для |
этого по |
|
ложим |
cp=cos (аг!-(-є) Ф (л:, |
у, |
z). |
(9.14) |
|
|
|||||
Подставляя значение <р в (9.4), |
(9.8) |
и |
(9.13), видим, что функ |
||
ция Ф должна удовлетворять уравнению |
|
|
|
||
|
^ " + ^ + |
^ = |
0 |
|
0М5> |
и краевым |
условиям |
|
|
|
|
|
- Т и - * |
|
|
(9.16) |
на твердых границах и
(9-17)
ПрИ 2 = 0.
Мы ограничимся плоским случаем, когда движение каждой час тицы происходит параллельно плоскости Oxz, а скорость v и давле ние р не зависят от координаты у . В плоском случае потенциал
скорости |
(9.14) есть |
<p=COS(a* + E)i]>(jc, z). |
|
(9.18) |
|||
|
|
|
|||||
Уравнение (9.15) |
записывается так: |
|
|
|
|||
|
|
<?2ф |
, |
<Э2ф |
О, |
|
(9.19) |
|
|
дх* |
~ |
дг* |
|
||
|
|
|
|
|
|||
а краевые условия |
(9.16) и (9.17) остаются |
без |
изменения. |
||||
Будем рассматривать водоем постоянной глубины h. Непосред |
|||||||
ственной |
подстановкой можно |
проверить, |
что |
уравнение (9.19) |
|||
имеет частное решение вида |
|
|
|
|
|
||
|
|
Ф(л-, z ) = P |
(2) sin £ (* -- /)), |
(9.20) |
|||
•если Р(г) |
удовлетворяет дифференциальному уравнению |
||||||
|
|
Р" |
|
(z)-k2P(z)=0. |
|
|
|
Из этого уравнения |
|
|
|
|
|
||
|
|
P(z)=Clek2~\-C2e-'iZ. |
|
|
|
||
Согласно краевому условию |
(9.16), должно быть дФ/дг = 0 при |
z — —h, т. е.
C,Ae* A - C 2 te - * f t =0 . Это условие удовлетворяется, если принять
т. е. |
|
P{z)=Cohk{z+h). |
(9.21) |
Из условия же (9.17) получаем |
|
e2=gk th kit |
|
Следовательно, в данном случае |
|
ср=С ch к (z+A)sin к {х - т)) cos (ot+e). |
(9.22) |
Интегрируя уравнение (9.11), найдем уравнение свободной по верхности
Здесь ар (л:)—произвольная функция. Но, согласно сделанному выбору системы координат, гр (х) = 0 и, следовательно,
е = - - у |
t ] L o = а s i n / г - ^ s i n ^ + є > - |
<9 -2 3 > |
где |
|
|
|
a = - ^ - c h A A |
(9.24) |
— амплитуда волны. Різ (9.23) видно, что в точках с абсциссами, удовлетворяющими условию кх— ц=пп (п — любое целое число или нуль), т. е.
хЪ '
вертикальные перемещения свободной поверхности отсутствуют.
Такие точки называются узлами. В |
точках, в которых kx — т) = |
= я я + я/2, амплитуды вертикальных |
перемещений равны а, т. е. |
достигают максимума. Эти точки называются пучностями. Волны, характеризующиеся наличием узлов и пучностей, носят название
стоячих волн. Длина |
такой волны K = 2n/k, |
а период Т = 2л/а. |
Проекции скорости частиц на координатные оси равны: |
||
a==~^-==-^-=Ck |
ch k (z 4-/г) cos k (x — -q) cos^-J-s), |
|
w—-^~—~^~ |
= Ck sh k (z-\-h)sink(x |
— y\) cos(at-\-z). |
Так как рассматриваются волны малой амплитуды, то в правых частях этих уравнений вследствие малых смещений частиц можно принять х и г равными тем значениям х = хо и г=го, которые час тица занимала в равновесном положении. При этом условии после интегрирования получаем следующие уравнения движения час тицы:
|
|
|
Ck |
|
|
|
|
|
|
х — х0=—— ch k(zQ-\-h) cos k(x0 |
— т}) sin (at-{-є), |
|
|
||||
|
|
|
Ck |
|
|
|
|
|
Отсюда z — z0— |
|
sh |
k (z0 - j - h) sink (x0 |
— 7j) sin (at-f- є). |
|
|
||
|
|
z |
z°- = |
|
thk(z0+h)tgk(x0-rl), |
|
|
|
|
|
|
••*o |
|
|
|
|
|
т. е. траектории |
частиц |
есть прямые линии. При этом |
в |
узлах |
||||
sin£(xo — т)) = 0 и |
частицы |
совершают |
колебательные |
движения |
||||
точно в |
горизонтальном |
направлении, |
в пучностях |
же, |
где |
|||
cos£(xo — т)) =0 , |
частицы |
колеблются вертикально. Амплитуда ко |
||||||
лебаний |
есть |
|
|
|
|
|
|
|
Ck , |
, ,ч |
sh k (z0 + h) |
В случае бесконечной глубины (/г-н»со) эта амплитуда будет aehz\ Отсюда видно, что волны с глубиной очень быстро затухают:
например, |
на глубине, равной длине |
волны, т. |
е. при Zo =—Х = |
=—2n/k, |
амплитуда будет а е _ 2 я , т. е. |
примерно |
в 500 раз меньше |
амплитуды на свободной поверхности. Иными словами, когда глу бина жидкости велика по сравнению с длиной волны, волнение но сит ярко выраженный поверхностный характер.
Мы считали до сих пор, что водоем простирается в направлении оси х неограниченно. В этом случае длины волн могут быть лю
быми. Если же он при |
х = 0 |
и |
x = L ограничен |
вертикальными |
|
стенками, нормальными |
к оси х, |
то, согласно условию |
(9.16), дол |
||
жно быть дФ/дх — О при х = 0 и при x = L , или |
|
|
|||
COS£T) = 0, |
COS k (L — v))=0. |
|
|
||
Отсюда следует |
|
|
|
|
|
= |
\-iriK, |
£(/,_ 7)) = -|-4-/даГ |
|
|
|
и значит |
|
|
|
|
|
где т — любое натуральное число или нуль. |
|
|
|||
Вернемся к случаю |
водоема, |
неограниченного |
в |
направлении |
оси л-. Нетрудно видеть, что всем уравнениям и краевым условиям
удовлетворяет не только потенциал |
(9.22), но и потенциал |
|
c?=Cch k(z-\-h)cosk(x |
— -q) sin (at-[-в). |
(9.25) |
В силу линейности уравнений и краевых условий разность потен
циалов |
(9.22) и |
(9.25), т. е. функция |
|
|
|
||||
|
|
|
|
<P=Cchfc(z + A)sin(fce —о/) |
|
|
(9.26) |
||
(здесь |
для простоты |
принято |
г| = е = 0 ) , также |
есть частное |
реше |
||||
ние |
рассматриваемой |
задачи. |
Уравнение свободной |
поверхности |
|||||
в этом случае будет, по аналогии с (9.23), |
|
|
|
||||||
|
|
|
|
|
£ = a c o s ( / b c — at). |
|
|
(9.27) |
|
|
Постоянные |
а |
и |
С по-прежнему связаны |
формулой |
(9.24). |
|||
Из |
(9.27) видно, |
что профиль |
волны, имеющий |
вид |
косинусоиды, |
||||
перемещается со скоростью |
|
|
|
|
Такие волны носят название прогрессивных волн. Интересно от метить два крайних случая:
1) если глубина жидкости очень велика по сравнению с длиной волны, то по соотношению
и m t h — г - = 1
имеем |
|
|
|
|
|
' - / - g |
r ; |
(9-29) |
|
2) если глубина жидкости очень мала по сравнению с длиной |
||||
волны, то |
|
|
|
|
,, |
|
2кк |
2тсЛ. |
|
t h |
— |
~ |
— |
> |
и тогда |
|
|
|
|
|
c--=Vgh, |
|
(9.30) |
т. е. скорость перемещения волн не зависит от их длины. Скорости и траектории жидких частиц находятся, как и в пре
дыдущем случае, по формулам: |
|
|
|
|
||
u=^=-^=Ck |
ch k{z+h) |
cos (kx - |
at), |
(9.31) |
||
•w=-^-=-^ |
= |
Ckshk{z^-h)sin{kx-at), |
|
(9.32) |
||
|
— |
ch k(z0+h) |
cos (kx0-at), |
|
(9.33) |
|
Ck |
|
|
|
|
|
|
z — zQ =—^— sh k (z0-f-/z)sin (kx0 — at). |
|
(9.34) |
||||
Исключая отсюда |
время |
и выражая |
С через |
амплитуду волн |
||
а, получаем уравнение |
эллипса |
|
|
|
|
|
( ^ ) 2 + ( ^ Н ' |
|
(9-35> |
||||
где |
|
|
|
|
|
|
|
А = а - ch k (г 0 + |
К) |
|
|
|
|
|
|
sh kh |
|
|
|
|
|
|
sh k (z0 + |
h) |
|
|
|
|
|
sh kh |
|
|
|
|
Таким образом, частицы двигаются |
по эллипсам с полуосями А |
|||||
и В. Для волн, распространяющихся в положительную |
сторону |
|||||
оси Ох, частицы вращаются |
по эллипсам |
по часовой стрелке. Как |
и в предыдущем случае, смещения частиц очень быстро умень шаются с погружением под свободную поверхность.
Обратимся к более общему случаю прогрессивных |
волн — вол |
нам на поверхности раздела двух жидкостей (рис. 8). |
Плоскость |
Оху совместим с поверхностью раздела жидкостей в состоянии равновесия. Пусть движение совершается между двумя параллель ными плоскостями z=—hi и z = hi, причем жидкости текут парал лельно оси Ох со скоростями Ui и ІІ2- Очевидно, что плотность Р2 верхней жидкости должна быть не больше плотности pi нижней