книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов
.pdfТогда |
|
|
2 • t mi сек |
|
|
([ |
км) = |
1000 м |
-0,004 м!сек. |
||
|
|
|
500-103 |
м |
|
Б течение суток это |
приведет |
к опусканию |
слоя |
воздуха на |
350 м и адиабатическому повышению температуры в нем па 3,5°С, что может вызвать полное рассеяние или значительное уменьшение облачности.
П р и м е р |
2. Линии тока параллельны друг другу, |
причем ве |
||||||
личина скорости в направлении движения убывает на |
1 м/сек на |
|||||||
каждые 100 км. |
В этом случае |
|
|
|
||||
div v |
|
dv |
\ v |
— 1 |
м/сек |
10-5 секг\ |
||
|
иг |
|
10 0 - |
10п м |
||||
|
|
|
|
|
||||
Таким образом, |
в этом случае на высоте 1 км имеют место вос |
|||||||
ходящие движения |
воздуха со скоростью порядка 1 см/сек. |
|||||||
П р име р |
3. |
Скорость ветра равна |
10 м/сек п вдоль линии тока |
|||||
не меняется. Линии тока, расположенные на расстоянии |
1 км друг |
|||||||
от друга, расходятся |
на 6° (0,1 радиана). В этом случае |
|||||||
div t'i = |
c |
оВ |
|
А? |
0.1 |
- 10_3 сек~\ |
||
ап |
|
v —т— = 10 м/сек ■—^— • |
||||||
|
|
|
А/; |
1 0 я м |
|
|
||
|
|
|
|
vz (1 |
км) — - -1 |
м/сек. |
|
|
т. е. на высоте 1 км имеют место нисходящие |
движения воздуха |
|||||||
со скоростью 1 м/сек. |
|
|
|
|
|
|||
|
|
|
|
|
ГЛАВА IV |
|
|
|
ЛОКАЛЬНАЯ, |
КОНВЕКТИВНАЯ И ИНДИВИДУАЛЬНАЯ |
|||||||
|
|
|
|
ПРОИЗВОДНЫЕ |
|
|
||
Здесь мы |
ставим |
своей |
задачей установить |
физический смысл |
и найти выражения производных по времени от плотности, темпе ратуры, проекций скорости и других скалярных характеристик жидкости, для которых введем унифицированное обозначение
/= / (х. У. х, /).
§1. ГРАДИЕНТ СКАЛЯРНЫХ ВЕЛИЧИН
Предварительно напомним, что эквискалярной поверхностью называется поверхность равных значений скалярной величины, т. е. / (х, у, г) =const. Если построены эквмекалярлые поверхности, соответствующие значениям функции /о, /о+Д/, /о+2 Д/ и т. д., т. с. такие, у которых переход от одной! поверхности к соседней сопро вождается изменением f на одну и ту же величину А/, то густота
'50 |
.а |
расположения поверхностей характеризует собой быстроту изме нения /: чем ближе друг к другу лежат поверхности, тем быстрее
меняется /.
Поместим теперь начало системы координат на какои-лиоо эквпскалярной поверхности (рис. 4.1), направив одну ось п к ней по нормали. Тогда две другие Л и к расположатся в соприкасаю
сь/ |
df |
щепря плоскости, за счет чего |
dT~ —в (в направлении |
этих осей значение / не меняется). |
|
П |
|
Запишем производную от / по любому направлению / в данной системе координат:
df |
|
df |
да |
| |
Of |
О/ , |
df |
d i\ |
||
dl. |
|
dn |
dl |
^ |
d f |
d! |
1 |
d t2 |
dl |
|
В указанной |
системе |
|
координат это эквивалентно равенству |
|||||||
|
|
df |
|
|
df |
|
л |
|
|
|
|
|
S |
r |
= |
- d T cos{',-‘>■ |
|
||||
d/i |
(/?, |
л |
|
|
|
|
|
|
|
(4.1) |
где - щ - — cos |
/), |
|
|
|
|
|
|
|
||
.Максимальное |
изменение |
/ |
будет |
иметь место в том случае, |
||||||
если / совпадает с п (cos |
(п, |
/) = 1 ). |
|
|
|
|
||||
Величина — |
|
называется |
градиентом |
скалярной величины /. |
Это вектор, характеризующий изменение / при смещении на еди ницу длины в направлении нормали. Как мы только что установи ли, модуль градиента равен производной функции по направлению ее быстрейшего изменения п.
4* |
51 |
Найдем теперь проекции градиента на осп любым образом за данной прямоугольной системы координат xyz. Для этого укажем, что формула (4.1) представляет собой запись проекции градиента на любое направление /. Если, в частности, / совпадает с осью х (I hs.v), то получим х-ю проекцию градиента:
Of Of
дх дп cos п, X)
Аналогично по осям у и 2 , соответственно:
/ |
QJ |
\ |
of . |
/ |
Of |
\ |
of |
\ |
дп |
/у |
ду ' |
\ |
дп |
/, |
dz |
Таким образом, |
|
|
|
|
|
|
|
|
Of |
|
Of |
> |
Of |
|
> Of |
|
дп |
1 |
дх |
' 1 |
д\) |
' |
(4.2) |
|
dz |
Модуль градиента:
I |
df |
(4.3) |
|
||
|
ду |
|
|
|
Отметим, что градиент часто еще обозначается как grad/. Через посредство оператора Гамильтона он может быть записан в виде V/. Таким образом
|
|
о]_ |
grad / = V/. |
|
|
Оп |
|
|
|
|
|
В качестве |
примера |
запишем градиент для случая, когда |
|
/ = vx или / ~ |
р. Ясно, |
что: |
|
> Op |
■> |
op |
_ |
> op |
grad p■ ' ~0x |
1 |
dy |
1 |
N ! |
У0 Vy |
->Ov у |
|
-> Ov.- |
|
grad Ux~ l Ox |
7 |
~dv |
+ k -oY- |
Аналогично:
Op |
I op |
Op |
op |
On |
~1~0x ’ |
~dy ’ |
dz |
O V у |
\Ovx |
O'Vy |
о v , | |
On |
\ dx ’ |
Ov ’ |
dz 1 |
dvy |
f О vv |
Ovy |
0 v y ) |
Ov, |
| |
Оч). |
o v z |
О v, j |
дп |
( дх ‘ |
dy ' |
dz j ’ |
On |
| |
Ox ’ |
dy ’ |
dz j ' |
52
■ Таким образом, можно констатировать, что градиент вектора
скорости grad и выражается через посредство совокупности девяти величин, т. е. является тензором второго ранга:
d v r |
d v x |
d v x |
дх |
ду |
dz |
d v y |
d v y |
d v |
дх |
ду |
dz |
dvz |
д v, |
dvz |
дх |
ду |
dz |
Отметим, что скалярное поле при переходе к градиенту поро ждает векторное, а векторное - тензорное с рангом два.
§2. ФИЗИЧЕСКИЙ СМЫСЛ ЛОКАЛЬНОЙ, КОНВЕКТИВНОЙ
ИИНДИВИДУАЛЬНОЙ ПРОИЗВОДНЫХ
Как известно, |
у, |
полный |
дифференциал какой-либо |
скалярной |
|||||||||
функции f = f (х, |
х, /) записывается в виде |
|
|
|
|||||||||
|
<//= |
df |
<и |
Of . |
|
Of |
dv |
!- |
ov |
|
|||
|
|
|
— |
dx |
- 1— —- |
— d z , |
|||||||
|
|
|
|
|
|
ах |
|
' |
dv |
- |
■ |
dz |
|
а полная производная по времени |
|
|
|
|
|
|
|||||||
|
d f |
df |
^ |
df |
dx |
| |
|
df |
dy |
, |
df |
dz |
|
|
dt |
dt |
|
1 |
dx |
dt |
|
|
dy |
dt |
' |
dz |
dt |
|
i |
|
|
dx |
|
dv |
|
|
dz |
|
\ |
|
|
или |
^поскольку |
- |
|
— гл |
dt |
= |
'V. |
~dt |
- |
) |
|
||
|
d f |
|
|
df |
|
df |
-Г ?'V |
rtf |
|
df |
|
||
|
dt |
|
~ |
dt |
^ |
l '- dx |
dy |
‘' r-- ~dz~ |
|
||||
Рассмотрим |
смысл |
каждого |
пз |
слагаемых |
в выражении (4.5). |
||||||||
1 . |
Локальная |
производная. |
Частная |
производная |
по времени |
——представляет собой изменение величины / за единицу вре
мени при фиксированных координатах, т. е. дает скорость изме нения свойства в неподвижной точке пространства. В силу этого она называется локальной. Знак производной характеризует собой воз-
( df |
\ |
I df |
\ |
растание |
|
или убывание l - ^ - C O ] / в0 времени. |
Показания приборов на станциях гидрометсети соответствуют измерениям локальной производной. Так, например, если имеются
53
обработанные наблюдения самописцев давления или температуры
dp |
ОТ |
(рис. 4.2), то для приближенного определения |
или v ^ до |
статочно взять на барограмме или термограмме приращение соот ветствующей величины за малый промежуток времени ЛI и разде
|
лить |
это |
приращение |
на |
||||
|
\/. |
Очевидно, |
|
что |
баро |
|||
|
метрическую |
|
тенденцию, |
|||||
|
т. е. приращение |
давления |
||||||
|
на станции |
за |
|
3 часа |
мы |
|||
|
можем |
рассматривать |
как |
|||||
|
приближенное |
значение |
ло |
|||||
|
кальной производной давле |
|||||||
|
ния, когда за единицу вре |
|||||||
Рис. -1.2. |
мени |
принят |
промежуток |
|||||
|
времени, |
равный |
3 |
часам. |
В случае установившегося движения, по определению которого, характеристики жидкости в любой зафиксированной точке про странства не претерпевают изменения.
2 . Конвективная производная. Это наименование предопреде ляется физическим смыслом рассматриваемой производной. Дей-
O '/ |
df |
df |
ствптельно, под ней подразумевается сумма vx~ ^ -f |
+ |
v -~dz- |
Of
При этом ^ представляет собой приращение f на единицу длины
df
в направлении оси Ox; vx — равно приращению / на расстоянии
их в том же направлении, т. е. приращению, обусловленному сме щением частицы за единицу времени на расстояние vx в направле- df
шш оси Ох. Подобно этому v,j |
выражает собой приращение /, |
обусловленное смещением частицы в направлении Оу, и аналогич- df
ное истолкование получает vz |
Сумма всех трех слагаемых равна |
приращению /, обусловленному одновременным смещением в на правлении каждой координатной оси за единицу времени.
Легко видеть, что конвективную производную можно рассма тривать как скалярное произведение двух векторов:
с = 1’х i-f-vu j -j-vz k,
51
grail/: |
df |
df ->• |
|
dx |
d\ |
||
|
n записывать ее в виде: (и-grad/).
Нетрудно и этому выражению дать физическое нстолковайне,- ибо
> |
-•> |
■> |
. |
—г |
(y-grad/) = |
[t>| j grad/ ( |
cos {v, |
grad /) = |
|z>[ grad/f, |
де grad//есть проекция grad/ на направление линии тока. В силу |
|||
Г а >- s |
|
|
|
этого |
|
|
|
|
|
df |
, |
|
|
(Ifо--ggrad/)raci i) = |»||-- ^ |
|
причем |
df |
определяет собой приращение / |
при смещении на еди- |
д/ |
df
вицу длины по линии тока. Поэтому —у представляет собой
приращение / при смещении по линии тока на отрезок | у|, т. е. приращение, обусловленное перемещением частицы по линии тока за единицу времени. Вывод тождественный предыдущему.
Конвективная производная обращается в пуль в следующих случаях:
1) когда о—О, что соответствует отсутствию перемещения ча стицы;
2 ) когда g r a d /= 0, что соответствует условиям, когда во всех точках, близких к данной точке, / имеет одну и ту же величину (так как производная в данной точке по любому направлению равна нулю). При этом изменение положения частицы в простран
стве само по себе не вызывает изменения /:
- v
3) когда o_Lgracl/. Это имеет место в том случае, когда частица смещается по эквпскалярнон поверхности функции /, и в этом слу чае изменение положения частицы не сопровождается изменением величины / в этой частице.
Наконец, |
отметим, |
что, если |
угол |
между |
v и grad/ |
острый |
||||
-> -А |
|
то / вдоль линии |
тока |
возрастает, |
в |
против |
||||
(cos(y, g r a d /) > 0), |
||||||||||
ном случае, наоборот, падает. |
|
|
Для |
выявления |
ее |
физиче |
||||
3. Индивидуальная |
производная *. |
|||||||||
ского смысла перепишем (4.5) |
в виде |
|
|
|
|
|
||||
df |
d f |
( |
df |
|
df |
|
df |
|
(4.6) |
|
df |
~ |
df |
Г'-1' |
dx |
+ Vy ~d7 |
+ |
~dz |
|
||
|
|
* Она еще носит название субстанциональной.
оо
Ясно, что изменение свойства / в точке может вызываться, во-первых, тем, что в нее непрерывно поступает жидкость из дру
гих точек. |
(Это отражает конвективная производная). Так, напри |
||||
мер, если |
в нее пришла чаетпна, |
находившаяся в |
момент i — 1 |
||
в другой точке, то она принесла |
в |
рассматриваемую значение |
|||
/= /, которое имелось у нее в исходном положении. |
Во-вторых, за |
||||
время движения свойства в самой |
частице |
также |
меняются за |
||
счет взаимодействия с окружающей средой. |
(Так, если происходит |
||||
поглощение тепла, то температура |
ее |
повышается). |
Этот послед |
ний фактор, как показывает (4.6), и отражает индивидуальная производная. Иначе говоря, она характеризует изменение рассма триваемой величины в фиксированной частице, которое происхо дит в процессе ес движения за единицу времени. Указанное изме нение оценивается по показаниям прибора, связанного с движу щейся массой.
Итак, согласно |
(4.5), d f |
равна сумме локальной и конвектив |
ной производных. |
В частном |
случае стационарного движения |
имеет место равенство индивидуальной и конвективной производ
ных, |
а в тех случаях, |
когда конвективная |
обращается в нуль, |
||||
dT |
— |
оТ |
п |
|
|
|
|
|
|
При этом изменение положения частицы в пространстве |
|||||
не вызывает приращения / в точке. |
|
|
|
||||
|
В заключение заметим, что с учетом модификаций записи кон- |
||||||
вективнои |
производной |
величина |
«7 |
может |
Л |
||
|
быть представлена |
либо в виде (4.5), либо как:
df df dt dt
или
dt dt
' (и-grad/), |
(4.7) |
->л |
(4.8) |
grad/| cos (u, grad/). |
§ 3. ПРИМЕРЫ
Пр и м е р 1. Определение барометрической тенденции по дан ным судовых наблюдений.
Выражение для индивидуальной производной оказывается справедливым не только применительно к движущейся частице данной среды, но и в том случае, когда нужно записать изменение во времени характеристики среды в каком-либо объекте, переме щающемся относительно среды произвольным образом, например
в самолете, на корабле и т. п. В этом случае в формуле (4.7) v представляет собой вектор скорости движущегося объекта, а в
5(3
формуле (4,5 ) фигурируют проекции этого вектора на оси коор динат. Такое обобщение позволяет применить эти формулы к ре шению практической задачи, с которой встречаются при анализе синоптических материалов.
Барометрической тенденцией называется изменение атмосфер ного давления в данном пункте за 3 час, т. е. локальная производ ная давления, при вычислении которой за единицу времени при нято 3 час. Задача ставится следующим образом. Известны ско рость и направление движения корабля, а также распределение давления на уровне моря (карта изобар). Зная изменение давле ния за 3 час, показываемое барометром па корабле, определить величину барометрической тенденции в пункте, где находится ко рабль.
Для определения локальной производной воспользуемся равен ством
др |
dp |
|
dt |
------- (t>-grad/>). |
|
|
|
|
Индивидуальная производная |
равна приращению давления |
за единицу времени, показываемому барометром на корабле. Конвективная производная определяется следующим выраже-
н нем:
|
|
|
|
(у-grad р) = |
tv • gradjр) |
OP |
COS а. |
||||
|
|
|
|
дпх |
|||||||
где gradiр — горизонтальный |
градиент давления; а — угол между |
||||||||||
курсом |
корабля и направлением градиента давления (рис. 4.3). |
||||||||||
|
Скорость |
корабля |
|
обычно |
|
|
|
||||
указывается в |
узлах. |
Поскольку |
|
|
|
||||||
1 |
узел— 1 |
миля/час=1,85 |
км/час, |
|
|
|
|||||
то |
v |
(км/час) — 1,85 |
v |
(узлов)-, |
|
|
|
||||
др |
|
Др |
|
|
с помощью |
|
|
|
|||
-— |
-г— определяется |
|
|
|
|||||||
дпх |
Дпх |
у |
карты, |
где |
вели- |
|
, |
|
|||
синоптической |
|
________ |
|||||||||
чине |
Ар |
соответствует |
|
разность |
________ дгаа |
||||||
давления на двух соседних изо |
|
Р||С- |
1-П| |
||||||||
барах, |
равная |
5 мб, а |
A/i| |
пред- |
|
||||||
ставляет |
собой |
расстояние |
между |
расстоянию |
на карте в санти |
||||||
изобарами |
в километрах, |
равное |
|||||||||
метрах, умноженному на 100. |
|
|
|
|
57
Таким образом,
|
* |
|
|
|
|
о |
|
|
|
|
(у-grad р) ==t.85 v (узлов) |
\/ii (см) ■100 cos а — |
|
||||||
|
|
— 0,090 |
с (узлов)cos а |
|
|
|
|||
|
|
|
\п\(см) |
(мб/час), |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
о (узлов) |
|
|
|
|
|
|
|
|
0,093 Д//1 (см) |
С О > |
rj_ |
||
|
|
|
|
|
|
|
|
||
и окончательно |
|
|
|
|
а (узлов) |
|
|||
dp |
I |
мб |
\ |
dp |
мо |
0,28 |
cos а |
||
dt |
\ 3 |
час |
' |
dt |
3 час |
:\п1 (см) |
|||
П р и м е р |
2. |
Изобары |
на синоптической |
карте представляют |
собой параллельные линии, расположенные в широтном направле нии па расстоянии 1 см друг от друга; давление растет с юга на север (рис. 4.4). Корабль в точке А движется в северо-западном направлении со скоростью 30 узлов, причем барометр на нем пока зывает рост давления на 2 мб за 3 часа. Определим величину ба рометрической тенденции в точке А:
|
др |
2 |
- 0,28 |
30 |
|
|
|
мб |
\ |
|||
|
'dt |
— cos 45°= --4 |
|
3 |
час |
) |
||||||
|
|
|
|
|
|
|
|
|||||
Таким |
образом, |
действительная |
величииа |
барометри чсекои |
||||||||
тенденции |
равна |
не |
+ 2 - |
мб |
как показывает оарометр на ко |
|||||||
час |
||||||||||||
|
|
мо |
|
3 |
|
|
|
|
|
|
||
рабле, а |
—4 |
.Легко |
понять, что отличие объясняется тем,' |
|||||||||
|
^ |
час |
|
|
что |
за 3 час |
корабль смещается |
|||||
|
q |
|
|
|
в сторону высокого давления на |
|||||||
|
|
|
|
|
|
90 |
миль или |
на |
63 мили (118 /сиг) |
|||
|
|
|
|
|
|
по нормали к изобарам. Смеще |
||||||
grad(jp |
|
|
|
|
ние корабля само но себе приво |
|||||||
|
|
|
|
|
|
дит к росту давления на 6 мб за |
||||||
|
~ v d / v |
|
|
3 час и этот |
рост |
перекрывает |
||||||
|
|
|
|
|
|
отрицательную |
барометрическую |
|||||
|
|
|
|
|
ДTlf тенденцию, т. е. общее падение |
|||||||
|
|
|
|
|
|
давления |
в |
|
рассматриваемом |
|||
|
|
|
|
|
|
районе. В силу этого барометр |
||||||
|
|
|
|
|
|
на |
корабле |
показывает величину |
Ю__ разности 6—4 = + 2 мб за 3 час.
Различие между показаниями барометра на корабле и истин ной величиной барометрической
тенденции, очевидно, оказывается тем большим, чем больше скорорость корабля, чем сильнее его курс отклоняется от направле ния изобар, чем гуще располагаются изобары.
58
§ 4. УСКОРЕНИЕ
1. Запись ускорения в декартовых координатах и векторно форме. Если последовательно положить f ^ v x, f ' = v y, f = vz, то, согласно (4.5), получим, соответственно, производные от проекций скорости на осп координат, т. с. проекции ускорения. Они имеют вид:
d vx |
— |
d v v |
-Ц |
|
<)v* |
I |
r |
d b \ |
, |
r |
C)V-V |
||
d t |
--- |
|
|
dy |
|
||||||||
|
dt |
|
r |
v d x + |
-v |
1 -• |
dz |
||||||
d vy |
|
di\, |
|
|
d vy |
|
|
O vy |
|
|
d vy |
||
d t |
~ |
dt |
|
4 - |
vx |
H t |
+ |
|
dy |
|
|
(4.9) |
|
d v. |
|
d v , |
+ |
r x |
d v, , |
vv |
d v . |
|
|
|
|||
d t |
~ |
dt |
|
■4 |
f- |
dv |
|
|
|
||||
|
|
|
d x |
|
|
|
|
|
|||||
В тензорных обозначениях три последние зависимости можно |
|||||||||||||
записать более компактно, в виде одной формулы: |
|
||||||||||||
|
|
|
|
d V: |
д V, |
|
|
д V; |
|
|
(4-9') |
||
|
|
|
|
d t |
|
d t |
' |
Vj |
д Xj |
|
|||
|
|
|
|
|
|
|
|||||||
Из приведенных формул мы видим, |
что |
|
проекции ускорения, |
||||||||||
а следовательно, |
и |
само |
ускорение, |
слагается |
из локального и |
конвективного ускорений. Первое может быть оценено по показа ниям прибора, способного непрерывно регистрировать скорость и неподвижно укрепленного в данной точке пространства. Приме рами таких приборов могут служить: в метеорологии — анеморумбограф, в гидрологии — гидрометрическая вертушка. Индиви дуальное изменение скорости, т. е. полное ускорение частицы, легко получить через посредство показаний регистратора скорости, не подвижно скрепленного с движущейся частицей. По разности этих величин можно определить и конвективную составляющую.
Умножая |
1-е, 2-е и 3-е уравнения |
(4.9) соответственно па /, /, к |
||||
и складывая |
результаты, |
получим |
выражение для ускорения |
|||
в векторной форме: |
|
|
|
|
|
|
|
|
d v |
dv |
4-(i>-y) v. |
(4.10) |
|
|
|
d l |
dt |
|
|
|
Здесь (ц-v) представляет |
собой |
скалярное |
произведение вектора |
|||
скорости па оператор Гамильтона, т. е. |
|
|||||
(y-V)= |
О , |
d |
|
d |
d |
|
Ж "г |
"у Н у |
|
d z ~~ v> d Xj ’ |
|||
|
|
|
59