
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfможет быть представлен как произведение скаляра X на единичный тензор
|
0 |
0 1 |
' 1 |
0 |
0 |
0 |
X |
0 • = \ - |
0 |
1 |
0 |
0 |
0 |
1 |
.0 |
0 |
1. |
т. е. он эквивалентен |
скаляру. |
|
|
|
Заметим, во избежание недоразумений, что тензор не определя ется однозначно тремя своими инвариантами. Чтобы вполне опре делить, например, симметричный тензор, нужно, кроме трех инва риантов, задать еще три главных направления, а они также зада ются тремя числами: двумя направляющими косинусами одного из направлений, например е( 1 ) , и одним направляющим косинусом од ного из двух других главных направлений, т. е. мы приходим к тому, что симметричный тензор определяется шестеркой чисел, как и в случае задания его в форме (3.7), а меняется только смысл этих чисел.
4. Криволинейные координаты
Вернемся к прежним обозначениям прямолинейных прямоуголь ных координат х, у, г. Если заданы три различные дифференцируе мые однозначные функции этих координат
|
Яі = |
ЯЛх, |
У, z), |
qi=q2{x, |
у, z), |
q3=q3(x, |
у, z) |
(4.1) |
|
таких, что обратные |
функции |
|
|
|
|
||||
|
x=x(qu |
q2, |
<7з), |
y=y(qu |
q2, q3), |
z=z(qu |
q2, q3) |
(4.2) |
|
также |
однозначны |
|
и дифференцируемы, то между |
тройками |
чисел |
||||
х, у, z |
и qi, qz, qz |
устанавливается |
взаимно |
однозначное соответст |
вие. Поэтому положение точки в пространстве можно определить заданием трех чисел qi, q%, qz, которые называются ее криволиней ными координатами. Очевидно, выбор функций qi, q2, qz не может быть вполне произвольным: требование взаимной однозначности со ответствия между криволинейными и прямолинейными прямоуголь ными координатами означает, что якобиан
dqi |
<>Я\ |
dqt |
|
|
дх |
ду |
dz |
|
|
dq2 |
|
dqo |
(4.3) |
|
дх |
ду |
dz |
||
|
||||
ддз |
<Э<7з |
дЯз |
|
|
дх |
ду |
dz |
|
|
должен отличаться от нуля. |
|
|
|
|
Три поверхности |
|
|
|
|
qx (х, у, £) = const, q2(x, |
у, z)=const, q3(x, у, z)=const |
|
называются координатными поверхностями. Линия пересечения двух координатных поверхностей есть координатная линия. Коор динатная линия <7i образуется пересечением координатных поверх ностей 92 = const и <7з = const; вдоль этой линии меняется только координата qu а координаты q2 и q3 сохраняют постоянные зна чения.
Введем единичные векторы Єї, е 2 , ез, направленные по касатель ным к координатным линиям в точке М (рис. 4). Вектор производ ной dr/dqi в точке М имеет, очевидно, направление Єї, ибо qi и q3 при дифференцировании сохраняют постоянные значения. По этому
где Li — длина вектора dr/dqi, и вообще
-*L.=Lfil(i=h |
2, 3). |
(4.4) |
|
||
А так |
как |
|
|
|
|
дг |
дх |
ду |
dz |
|
|
|
|
|
dqt |
|
|
TO |
дх |
ду |
|
|
|
|
|
|
|
||
|
|
dqt |
|
|
Система криволинейных ко |
|
|
|
|
(4.5) |
ординат. |
|
|
|
|
|
|
Величины |
L l t Li, L 3 |
называются |
коэффициентами Ламэ. |
||
В данном |
параграфе |
мы |
будем |
рассматривать только ортого |
нальные криволинейные координаты, т. е. такие, что в каждой точке
пространства направления |
e i , е 2 , е 3 |
взаимно |
перпендикулярны или |
иначе |
|
|
|
(е„ |
еА ,)=0, |
іфк. |
(4.6) |
Необходимое и достаточное условие ортогональности криволи нейных координат есть, очевидно, условие ортогональности векто ров dr/dqi и dr/dqk при 1фк, т. е.
дх |
дх |
, ду |
ду , dz |
dz |
Q |
(4.7) |
|
дЯі |
dqu ~ |
dqt |
dqk |
dqb |
~~ |
||
|
Однако в решении конкретных задач обычно не приходится проверять ортогональность так или иначе выбранной системы коор динат, а наоборот, выбирать систему координат, адекватную задаче из условия ортогональности координатных линий. В таких случаях условия (4.7) выполняются автоматически и проверять их не при ходится.
Вектор F в системе ортогональных криволинейных координат можно представить разложением
F^Ffr+F^+Fte |
(4.8) |
Числа Fx, F2, F3 называются криволинейными составляющими вектора или проекциями вектора F на оси криволинейных коорди нат. Беря в качестве вектора F вектор а'г, получаем из (4.8) и (4.4)
dr=-~ |
dqx + ^ d q 2 + -^-dq3=Lx |
dqxQx - f Z.2 dq2e2 + L 3 |
dq3e3. |
|
Таким образом, составляющие вектора dr есть |
(4.9) |
|||
|
||||
|
ds—L.dq, |
{І=Л, |
2, 3). |
(4.10) |
В конкретных задачах определение dsi по dqx обычно легко вы полняется по свойствам выбранной системы координат, поэтому коэффициенты Ламэ большей частью определяются из формулы
(4.10), а не из формулы (4.5). |
|
|
|
|
|
|
|
|
|||
Пусть |
N — точка, бесконечно |
близкая |
к М, |
так |
что |
dr=MN |
|||||
(рис. 5). Проведем |
через точки М и N по три координатных поверх |
||||||||||
|
|
|
|
ности. Эти |
поверхности образуют |
||||||
|
|
|
УУ бесконечно |
малый |
криволинейный |
||||||
|
|
|
1 |
параллелепипед |
с |
ребрами |
ds^ |
||||
|
|
|
|
Площади doi, do*, daa |
его |
граней, |
|||||
|
|
|
|
нормальных |
к |
направлениям |
е ь |
||||
|
|
|
|
ег, е3 , соответственно |
будут |
|
|
||||
|
|
|
|
do1=L2L3dq2dq3, |
|
|
|
do2—L3Lxdq3dqx, |
|||
|
|
|
|
|
da3=LxL2dqxdq2. |
|
(4.11) |
||||
Рис 5. |
К |
определению элемента |
а |
объем dV определится по формуле |
|||||||
объема |
н |
элемента |
поверхности |
|
dV—LxL2L3dqxdq2dq3. |
|
(4.12) |
||||
в криволинейных координатах. |
|
|
Выразим теперь основные дифференциальные операции в кри волинейных координатах.
По определению grad Ф, проекция этого вектора на ось qi есть дФ/dsi или, согласно (4.10),
|
|
дФ |
1 |
дФ |
|
|
|
|
|
|
dsi |
Li |
dqi |
|
|
|
|
Следовательно, |
|
|
|
|
|
|
|
|
ь |
Ц |
dqx |
' L 2 |
dq2 1 L 3 |
dq3 |
4 |
' |
|
Отсюда видно, |
что для |
оператора |
Гамильтона |
мы |
имеем |
вы |
||
ражение |
|
|
|
|
|
|
|
|
V |
Ц |
dqx |
' L 2 |
dq2 |
Т I 3 dq3 |
' |
|
|
|
|
|
|
|
|
|
|
Для выражения других операций поля в криволинейных коор динатах существуют три пути:
1) можно пользоваться формулами (2.16), (2.24) и т. д., вы ражающими эти операции через оператор Гамильтона; при этом
необходимо учитывать, что единичные векторы криволинейных ко ординат е ь Є2, ез сами есть функции координат; последнее делает этот путь очень громоздким;
2)можно исходить из определений операций поля, т. е. из вы ражений (2.16), (2.24) и т. д., рассматривая их в криволинейных координатах так же, как это было сделано в § 2 и 3 для прямоли нейных координат;
3)можно апеллировать к общей теории тензоров.
Любой из этих путей позволяет получить следующие выражения для основных дифференциальных операций поля в криволинейных координатах.
1. Криволинейные составляющие wt ускорения w точки, движу щейся со скоростью v — обобщение формул (2.13) для случая F = v:
|
З |
|
|
З |
о |
|
|
З |
|
n , , - * L i . V A |
|
dvi |
V |
vk |
dLk |
і V |
v(Ok dLt |
||
2. Расхождение |
вектора |
F — обобщение |
формулы |
(2.16) |
|||||
d i v F |
1 |
Г d{L2L3Fx) |
|
. d{L3LxF2) |
, d(LxL2F3) |
) |
|||
|
LXL2L3 |
\ |
dqx |
|
|
dq2 |
"r" |
dq3 |
f' |
3. Криволинейные составляющие вихря вектора F — обобщение формулы (2.22)
|
1 |
f |
d(F3L3) |
(rotF)1 = = |
|
{ |
dq2 |
|
1 |
f |
d(FxLx) |
(rot F ) 2 = |
L3LX |
[ |
dq3 |
(rot F);3" |
1 |
f |
d(F2L2) |
Lit* |
|
dqi |
д {F2L2) дЧз
д (F3L3) dq\
d{FxLx) dq3
4. Оператор |
Лапласа |
от |
|
скаляра |
Ф — обобщение |
формулы |
|||
(2.26) |
|
|
|
|
|
|
|
|
|
у , 2 ф г ^ |
1 |
Г Д / L2L3 |
|
_ дФ \ • |
д ( L3LX |
_ дФ \ . |
|||
|
LXL2L3 |
\ dqi { |
L x |
|
dqi j |
" 1 |
" dq2 [ L 2 |
dq2 |
J ' |
|
|
і |
д |
I LX 2 |
дФ |
|
|
||
|
|
' |
dq3 |
\ |
L 3 |
dq3 |
|
|
5. Криволинейные составляющие вектора div T — расхождения тензора Т с криволинейными компонентами Fih — обобщение фор мулы (3.21)
|
|
3 |
|
|
|
|
|
|
(Л\\т Т\ |
1 |
V I |
1 |
д |
I L\LiLbLi |
77 \ |
Fkk |
dLk \ |
3 Заказ № 428
Глава |
II |
О С Н О В Н Ы Е |
УР А В Н Е Н И Я Г И Д Р О М Е Х А Н И К И
5.Основные уравнения механики сплошной среды
Будем рассматривать некоторую (не обязательно жидкую) дви жущуюся сплошную среду, пронизывающую в своем движении про странство, в котором задана какая-либо неподвижная система ко
ординат. Эту систему для простоты можно считать |
прямолинейной |
||
прямоугольной, не нарушая этим общности |
построений. Выделим |
||
в пространстве тетраэдр (рис.6) |
|||
МАВС, |
ребра |
MA, MB, |
МС |
которого |
параллельны |
осям |
координат х, у, z соответст венно. Если 5 — площадь грани
ABC, а |
п — внешняя |
нормаль |
к этой |
плоскости, то |
площа |
ди граней МВС, MAC и МАВ
будут равны соответственно
5 cos
X •
Рис. 6. К выводу динамического урав нения механики сплошной среды.
Силы, действующие на грани ABC, МВС, MAC и МАВ тетра эдра со стороны остальной среды могут быть выражены так:
phSr p*S cos |
py S cos |
p^Scos I |
где pn, p*, py, pz — средние значения сил, приходящихся на единицу площади соответствующих граней. Эти средние значения есть не которые векторы, направления которых не совпадают в общем слу чае с направлениями нормалей к граням. Из условия равновесия сил, действующих на тетраэдр, имеем
pnS— p.v5cos1 |
py Scos |
pzS |
cos |
) + V q = 0 , |
|
где q — вектор |
объемных сил (например; |
силы |
тяжести) |
и сил |
|
инерции на единицу объема тетраэдра; |
V — его полный объем. Раз |
||||
делив это равенство на S, устремим все ребра тетраэдра |
к нулю. |
||||
Учитывая, что при этом отношение V/S |
стремится |
к нулю |
(числи- |
тель — малая |
третьего порядка, |
а знаменатель — второго), полу |
||
чаем |
> |
|
|
|
PH = |
PA-COS ( i O ) + Py cos (nTy)-f-p;Cos |
(tCz). |
||
Векторы p.v, ру, p z можно представить в виде |
|
|||
|
Px = |
Wxx + |
iPxy+^Pxz, |
|
|
Py = |
iPyx + iPyy + kpv2, |
|
|
|
Pz = |
iPzx + |
iPzy+kpz2- |
|
Соста вляющие рхх, рху, Pxz вектора р.-с єсть |
силы, |
действующие |
|
на единицу |
поверхности площадки, перпендикулярной |
оси х в на |
|
правлениях |
осей х, у, z, т. е. рхх — нормальное, |
а рху |
и pXz — каса |
тельные напряжения по этой площадке. |
|
|
Таким образом, в силу определения, данного в § 3, напряженное состояние в произвольной точке М пространства занятого сплош
ной средой есть тензор |
|
|
|
Рхх |
Рху Pxz' |
|
|
П = / V |
Руу |
PyZ • , |
(5.1) |
Pzx |
Pzy |
Pzz. |
|
называемый тензором напряжений.
Перейдем к выводу основного уравнения механики сплошной среды. Выделим для этого некоторый объем V без пустот, ограни ченный поверхностью S, неподвижный относительно системы коор динат.
Сумма сил, приложенных |
к субстанции, |
заполняющей |
этот |
объем в любой момент, включая силы инерции, должна быть |
равна |
||
нулю |
|
|
|
J p(F-w)dV |
+ <$pndS=0. |
|
(5.2) |
V |
s |
|
|
Здесь р — плотность среды; |
w — ускорение |
ее частиц; F — уско |
рение объемных сил (например, силы тяжести). Момент указан ных сил относительно любой точки пространства, например, отно
сительно начала |
координат |
0, |
также |
должен |
быть равен |
нулю. |
|||
А так как момент силы |
Р есть |
векторное произведение этой |
силы |
||||||
на плечо (радиус-вектор) |
г точки |
ее приложения, |
то это второе ус |
||||||
ловие записывается так: |
|
|
|
|
|
|
|
|
|
J |
р [г, F - |
w] d V+ |
j> [г, |
р„] dS=0. |
|
(5.3) |
|||
V |
|
|
|
S |
|
|
|
|
|
По формуле |
(3.23) |
|
|
|
|
|
|
|
|
|
Ф Ря ^ |
= |
J dlv ПrfV. |
|
|
(5.4) |
|||
|
5 |
|
|
|
v |
|
|
|
|
Поэтому и в силу произвольности |
объема |
V, |
уравнение |
(5.2) |
|||||
будет иметь вид |
|
|
|
|
|
|
|
|
|
|
p ( F - w ) + divn=0 . |
|
|
(5.5) |
3* |
-35 |
Это есть первое основное уравнение механики сплошной среды. Условие (5.3) равенства нулю момента приложенных к объему сил показывает, что П в уравнении (5.5) есть симметричный тен зор. Действительно, согласно формулам (3.23) и (3.21), а также вследствие того, что векторное произведение дифференцируется по тем же правилам, что и обычное, имеем
|
|
|
У \ |
|
дх |
* |
ду |
' |
дг |
|
|
|
|
i l l ' |
дх |
1 |
ду |
|
|
|
+ . дх ' P v |
+ " дудт' |
+ |
||
+ dzдг |
p]}dV=U[r, |
|
|
divn] + |
[i, p,] + [j, |
p,l + |
[k, |
pz]}d V. |
||||
Теперь уравнение |
(5.3) |
можно записать так: |
|
|
||||||||
J [г, P ( F - w ) + d i v n ] d \ / + j ( i ) |
p,] + |
[j, p,] + [k, |
|
p,]}dV=0. |
||||||||
V |
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.6) |
Из |
(5.5) следует, что первый интеграл |
в (5.6) равен нулю. Атак |
||||||||||
как объем V произволен, то уравнение (5.6) сводится |
к |
уравнению |
||||||||||
|
|
|
[i. |
P.v] + [J, |
Py] + [k, |
Р , ] = 0 . |
|
(5.7) |
||||
Но |
уравнение |
(5.7) |
и |
есть |
условие |
симметричности тензора. |
||||||
Чтобы в этом убедиться, умножим |
(5.7) |
векторно на |
произвольный |
|||||||||
вектор |
а. Учитывая, |
что |
согласно |
формулам |
(1.3) и (1.4) |
|||||||
|
|
|
[і, |
р * ] Х а = р Л * . а ) - і ( р Л |
а) |
|
|
и что аналогичные соотношения имеют место для других членов
(5.7), получаем |
|
|
РЛі» a)+py (j, a)+p,(k, |
а ) - і ( р х , |
a ) - j ( p , , a ) - k ( p z , a ) = 0 , |
откуда, согласно формулам |
(3.14) и |
(3.15), имеем |
(Пс , а ) - ( П , |
а ) = 0 . |
А так как а — произвольный вектор, то это значит, что П |
С = П, |
т. е. тензор П совпадает со своим сопряженным, и значит он |
сим |
метричен
Рху==РуХ> Pyz==Pzy> Pzx===Pxz-
Если рассматривается жидкая (капельная или газообразная) среда, то, положив в формуле (2.12) F = v, получим следующее вы ражение для ускорения:
Основное уравнение (5.5) для этого случая будет иметь вид
- g - + ( v , V v ) = F + - I - d l v I I . |
(5.9) |
|
Следует заметить, что |
оператор (v, V ) часто |
обозначают как |
(v, grad), тогда уравнение |
(5.9) будет иметь вид |
|
- ^ - + (v, grad)v = F + -£- divH.
Важный частный случай жидкой среды есть идеальная жид
кость — среда, неспособная |
к |
восприятию касательных напряже |
|
ний. В идеальной |
жидкости |
р х у |
= pyx = pVz = Pzy = pzx = Pxz:=0, а век |
торы р п , рл;, ру, |
pz нормальны |
к граням тетраэдра соответственно |
ABC, МВС, МАС, МАВ |
(рис. 6), при этом |
|
|
Vx = |
Wxx, Py = j/>yy. |
Pz=kPzz, |
|
P„=iA™cos ( n T x ) - f j^yycos (тСу) + кргг cos (пГ*)- |
(5.10) |
||
Проектируя равенство (5.10) последовательно на направления |
|||
х, у, г, получаем |
|
|
|
|
Pn=Pxx=Pyy=Pzz, |
|
(5-11) |
т. е. нормальное напряжение не зависит от направления площадки, по которой оно действует, иными словами, в данном случае тензор П есть шаровой тензор
\-р |
0 |
01 |
||
п = |
|
О -р |
о |
, |
( |
о |
0 |
~р\ |
|
где — р — инвариантное по |
отношению |
к |
ориентации площадки |
значение нормального напряжения. Величина р есть скаляр, на зываемый гидродинамическим давлением или просто давлением. Знак минус при р появился потому, что сила, действующая на любую грань тетраэдра, направлена всегда навстречу внешней нор мали к этой грани.
Если |
П — шаровой тензор, |
то, согласно |
формулам |
для |
дивер |
||
генции |
тензора, |
d i v l l = gradp |
и уравнение |
(5.9) |
для |
идеальной |
|
жидкости будет выглядеть так: |
|
|
|
|
|
||
|
|
- 5 T + ( v , V ) v = F — J - g r a d j * . |
|
|
(5.12) |
||
Это уравнение |
называется уравнением Эйлера |
(точнее, |
уравне |
ниями Эйлера называются три проекции этого векторного уравне ния на три координатные оси).
Если жидкость не идеальна, и в ней могут существовать ка сательные напряжения, то соотношение (5.11) уже не имеет места.
Ho pxx+Pyy + Pzz есть инвариант тензора |
П, т. е. скаляр. Положим |
в этом случае |
|
Р = --Ъ-(Рхх+Руу+Ргг) |
(5-13) |
и введем следующие обозначения для компонентов тензора П:
Рхх=—Р+*хх, |
Руу=—Р+*уу. |
Ргг=—Р |
+ 1г*> |
Рху ==/^'y.v= = Zxy= = ^у і-1 Pyz' |
Pzy ' ^уг= = |
^гуі |
|
|
/'г.ї —Агг= = xzx |
==тл"г • |
|
Теперь уравнение (5.9) можно записать в виде
d v ' (г>, g r a d ) x » = F — f g r a d p + - f div7\ |
(5.14) |
где
^xy
7=={ XyX Lyy
г-
''zx ^y ••zz
— тензор, порожденный отклонениями свойств жидкости от иде альной. Скаляр р в этом уравнении по-прежнему называется дав лением.
Второе основное уравнение сплошной среды дает закон сохране ния массы. Рассмотрим действие этого закона в условиях проте кания жидкости через тот же объем V, который имелся в виду при выводе первого основного уравнения. За время dt через поверх ность 5, ограничивающую этот объем, вытекает наружу масса жид кости
<J>(n, |
pv)dSdt. |
s |
|
С другой стороны, за то же время масса жидкости в объеме V изменяется на
1-%-dVdt. v
Общее приращение массы за время dt есть сумма этих двух вы ражений. По закону сохранения вещества это приращение равно нулю, следовательно,
j-|*-rfVcW+<j>(n, pv)dSdt=0.
Сокращая это уравнение на dt и преобразуя затем поверхност ный интеграл в объемный по формуле (2.27), получаем
J ^ + d i v P v = 0 . |
(5-15) |
Это уравнение называется уравнением неразрывности. Оно мо жет быть записано также в форме
д? |
•V • grad p-j-p div v = 0 , |
|
dt |
||
|
или, согласно формуле (2.10),
1 |
RFP ' - d i v v = 0 . |
(5.16) |
p |
dt |
|
Обращаясь к физическому смыслу формул (2.9) и (2.10), за ключаем, что dp/dt характеризует изменение плотности жидкой ча стицы при ее перемещении в потоке. Если жидкость несжимаема, то плотность каждой ее частицы остается постоянной, т. е.
- ^ = = 4 - + v - g r a d p = 0 . |
(5.17) |
и уравнение неразрывности приобретает вид |
|
div v = 0 . |
|
При этом, конечно, ни dp/dt, ни grad р, вообще говоря, |
в нуль |
не обращаются, ибо неизменность плотности движущихся |
жидких |
частиц не означает неизменность плотности жидкости в данной точке пространства, если плотность этих частиц почему-либо раз лична.
По характеру изменения плотности мы будем различать сле дующие виды жидкостей.
А. Несжимаемые жидкости:
1)жидкость постоянной плотности;
2)общий случай несжимаемой жидкости. Б. Сжимаемые жидкости:
1)баротропная жидкость. Плотность ее зависит только от дав ления р = р(р). Простейший пример — упругая капельная жидкость, подчиняющаяся при сжатии закону Гука;
2)бароклинные жидкости (в частности, все газы). Плотность
еезависит не только от давления.
В механике деление жидкостей по этому принципу характери зует не столько физические свойства жидкостей, сколько свойства математических моделей, применяемых для описания их движения: одна и та же жидкость в различных задачах может рассматри ваться, например, как несжимаемая и как бароклинная.
6. Уравнения Навье—Стокса
Векторное уравнение (5.14), эквивалентное трем скалярным, и скалярное уравнение (5.15) образуют систему четырех скалярных уравнений, в общем случае с одиннадцатью неизвестными: и, v, w (так мы обозначим здесь компоненты вектора скорости v по осям X, у, Z), р, р, Гхх, tyy, Tzz, txy, Тх г , Xyz. Общие ЗЭКОНЫ МеХЭНИКИ, ПрИВОдящие к этим уравнениям, никаких других независимых уравнений