лекция 4
.pdfЛекция 4
Наименование |
|
|
|
|
|
|
|
(час) |
|
Содержание теоретической части |
|
Затраты времени |
|||||
разделов |
Номер |
|
|
|||||
|
Самостоятель- |
|||||||
|
разделов дисциплины |
|
лекции |
Лекции |
ная (внеаудитор- |
|||
дисциплины |
|
|
|
ная) работа |
||||
|
|
|
|
|
|
|
студента |
|
ОСНОВНЫЕ |
Свойства напряжений поверхностных сил. |
4 |
2 |
0,5 |
||||
УРАВНЕНИЯ |
Давление и его свойства. |
|
|
|
|
|
||
ГИДРОГАЗО- |
Нормальные и касательные напряжения. Закон |
|
|
|
||||
ДИНАМИКИ |
парности касательных напряжений. Модель |
|
|
|
||||
|
идеальной жидкости. Давление и его свойства. |
|
|
|
||||
|
Уравнение движения (количества движения или |
|
|
|
||||
|
импульса) жидкости в напряжениях. |
|
|
|
||||
|
Обобщенная гипотеза Ньютона о связи между |
|
|
|
||||
|
напряжениями |
и скоростями деформаций |
|
|
|
|||
|
(закон Стокса). |
|
|
|
|
|
|
|
|
Уравнения движения Навье-Стокса для |
|
|
|
||||
|
вязкой сжимаемой жидкости. |
|
|
|
|
|||
|
Уравнения движения идеальной жидкости в |
|
|
|
||||
|
форме Эйлера и Громеки-Лэмба. |
|
|
|
|
|||
|
Интегралы |
уравнений |
движения |
идеальной |
|
|
|
|
|
жидкости: для безвихревого неустановившегося |
|
|
|
||||
|
движения (интеграл Лагранжа или Коши- |
|
|
|
||||
|
Лагранжа), |
для |
установившегося |
движения |
|
|
|
|
|
(интеграл Бернулли). Уравнение Бернулли. |
|
|
|
||||
|
Интегральная |
форма |
закона |
сохранения |
|
|
|
|
|
количества движения (импульса) для жидкого |
|
|
|
||||
|
объема (первое уравнение Эйлера). Переход от |
|
|
|
||||
|
жидкого объема к контрольному. |
|
|
|
|
|||
|
Определение сил, действующих на твердое |
|
|
|
||||
|
тело, по состоянию потока на границах. |
|
|
|
||||
|
Полный импульс потока в сечении. Определение |
|
|
|
||||
|
усилия, действующего на стенки криволинейного |
|
|
|
||||
|
канала со стороны текущей по нему жидкости; |
|
|
|
||||
|
учет сил давления на канал со стороны |
|
|
|
||||
|
окружающей среды. |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 1 ~ |
Свойства напряжений поверхностных сил. Давление и его свойства.
Возьмем какую-либо точку M сплошной среды, проведем из этой точки координатные оси и построим на них бесконечно малый тетраэдр ABCM с ребрами dx , dy , dz . Ориентировка грани ABC произвольна и задается вектором
нормали n . Обозначим площадь грани ABC через dS , а площади трёх граней тетраэдра, лежащих в координатных плоскостях через dFx , dFy и dFz .
Рис. 4.1
Составим векторное уравнение движения жидкого тетраэдра, выражающее второй закон Ньютона. При этом учтем, что положительными принято считать напряжения, направленные в сторону внешней нормали, а орты координатных осей, являющиеся нормалями к площадкам dFx , dFy и dFz , направлены внутрь
выделенного жидкого объема.
На тетраэдр ABCM будут действовать поверхностные силы, характеризуемые напряжениями px , py , pz и pn , приложенными на
соответствующих гранях тетраэдра, а также массовая сила
dR = F dm = ρFdV = ρF 13 hdS ,
где R – напряжение массовой силы, м/с2; dm – масса в объеме тетраэдра dV , h – высота тетраэдра.
В соответствии со вторым законом Ньютона сумма сил, действующих на тетраэдр, равна произведению его массы на ускорение, то есть
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 2 ~ |
1 |
r |
r |
r |
r |
r |
= |
dυ |
dm = |
1 |
ρh |
dυ |
dS |
|
3 |
ρhFdS + pndS - pxdFx |
- pydFy |
- pzdFz |
dt |
3 |
dt |
(4.1) |
||||||
|
|
|
|
|
|
|
|
|
|
Площади граней, параллельных координатным осям, можно вычислить как площади проекций грани ABC на соответствующие координатные плоскости:
dFx |
= dS cos(n, x) = nxdS ; |
|
dFy |
= dS cos(n, y) = nydS ; |
(4.2) |
dFz |
= dS cos(n, z) = nz dS . |
|
Здесь nx = n × i = ni cos(n, x) = cos(n, x), аналогично ny , nz .
В результате, подставив (4.2) в (4.1) сокращаем все члены равенства (4.1) на dS и стягивая тетраэдр в точку (то есть полагая h → 0), получим
pn = nx px + ny py + nz pz . |
(4.3) |
С другой стороны, каждый из векторов px , py , pz может быть представлен
тремя составляющими, действующими вдоль координатных осей.
Обратимся к рис. 4.1. На грань dFx действует напряжение px , для которого
можно записать
px = iσ x + jτ xy + kτ xz ;
аналогично
py = iτ yx + jσ y + kτ yz ; |
(4.4) |
pz = iτ zx + jτ zy + kσ z . |
|
Первый индекс указывает ориентацию площадки (нормаль к ней), а второй – ось, на которую проецируется вектор. Для нормального напряжения второй индекс совпадает с первым, поэтому может опускаться.
Теперь, учитывая (4.4), спроецируем вектор |
pn на координатные оси. |
Получим: |
|
pnx = nxσ x + nyτ yx + nzτ zx |
|
pny = nxτ xy + nyσ y + nzτ zy |
(4.5) |
pnz = nxτxz + nyτ yz + nzσ z |
|
Таким образом, напряженное состояние в точке определяется совокупностью трех векторов напряжения или их девятью компонентами, определенными на трех взаимно перепендикулярных площадках. Соотношение (4.5) является определением тензора (от англ. tension - напряжение).
Компоненты образуют тензор второго ранга, которому можно поставить в соответствие матрицу:
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 3 ~ |
|
æ |
σ x |
τ yx |
τ zx |
ö |
|
|||
pni |
ç |
|
|
σ y |
|
|
÷ |
(4.6) |
|
= çτ xy |
τ zy ÷ . |
||||||||
|
ç |
τ |
xz |
τ |
yz |
σ |
z |
÷ |
|
|
è |
|
|
|
ø |
|
Диагональные компоненты (с одинаковыми индексами) называются нормальными напряжениями, а остальные компоненты – касательными напряжениями, или
напряжениями сдвига.
Тензор напряжений (4.6) зависит от координат и времени и образует тензорное поле.
Если составить уравнения моментов относительно координатных осей, можно показать, что для касательных напряжений справедливы соотношения:
τ xy =τ yx ; τ yz =τ zy ; τ zx =τ xz , |
(4.7) |
называемые законом парности касательных напряжений. Соответствующее доказательство можно найти, например, в [1, 2] из списка основной литературы. Из этого закона следует, что тензор напряжений (4.6) является симметричным. Последнее означает, что тензор напряжений содержит только шесть различных компонент. В соответствии с этим уменьшается число неизвестных в уравнениях движения.
Рассмотрим теперь модель идеальной жидкости. Идеальной жидкостью
называется изотропная сплошная среда, в которой отсутствуют касательные напряжения, то есть τij = 0 (i ¹ j ). Касательные напряжения в жидкости
возникают благодаря трению. Поэтому можно сказать, что идеальная жидкость – это жидкость, лишенная внутреннего трения.
Тогда, из (4.5) следует, что
pnx = nxσ x , |
pny |
|
= nyσ y , pnz = nzσ z . |
(4.8) |
||||||
С другой стороны, |
|
|
|
|
|
|
|
|
|
|
pnx |
= pn × i = |
|
pn |
|
× |
|
i |
|
cos(n, x) = pnnx . |
|
|
|
|
|
|
||||||
|
|
|
||||||||
Аналогично, |
= pnny , pnz = pnnz . |
|
||||||||
pny |
(4.9) |
|||||||||
Сравнивая (4.8) и (4.9), видим, что |
|
|
|
|
|
|
|
|
|
|
pn = σ x = σ y |
= σ z . |
(4.10) |
Жидкий объем в силу своего строения не может воспринимать растягивающие напряжения. Поэтому действующие на него нормальные напряжения должны быть сжимающими. Такое сжимающее нормальное
напряжение называется абсолютным давлением или просто давлением и
обозначается p . Следовательно,
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 4 ~ |
p = − pn = −σ x = −σ y = −σ z |
(4.11) |
Знак «минус» показывает, что давление всегда направлено по внутренней нормали к поверхности выделенного жидкого объема. Это первое свойство давления.
Обозначая px = −σ x , py = −σ y , pz = −σ z , можно записать p = px = py = pz .
Осюда следует второе свойство давления: в любой точке идеальной или
покоящейся вязкой жидкости, а также при движении вязкой жидкости с равномерным полем скоростей давление считается величиной положительной и не зависит от ориентации в пространстве площадки, к которой оно приложено, то есть давление – величина скалярная.
Но в общем случае p = p(x, y, z,t) , то есть давление может изменяться в
пространстве (в разных точках разное давление) и во времени (в одних и тех же точках в разные моменты времени давление может меняться).
Уравнение движения жидкости в напряжениях.
Применим к произвольному жидкому объему V , ограниченному поверхностью F , известную из теоретической механики теорему о количестве движения системы: производная по времени от количества движения системы равна сумме действующих на нее внешних сил.
Количество движения K массы жидкости в объеме V равно
K = òρWdV ,
V
где W - скорость движения центра масс жидкого объема dV . Тогда
dK |
|
d |
|
r |
d |
r |
dW |
r |
d |
(ρdV ). |
|
|
= |
|
òρWdV = ò |
|
ρWdV = òρ |
|
dV + òW |
|
|||
dt |
dt |
dt |
dt |
dt |
|||||||
|
V |
V |
V |
V |
|
Так как ρdV - масса жидкого объема dV , не изменяющаяся со временем, то есть
|
|
|
d |
(ρdV ) = 0 , |
(4.12) |
|
|
|
dt |
||
|
|
|
|
|
|
r |
d |
(ρdV ) = 0 . |
|
||
то соответственно и òW |
|
|
|||
dt |
|
||||
V |
|
|
|
|
Тогда уравнение количества движения жидкого объема V в интегральной форме запишется в виде:
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 5 ~ |
òρ |
dW |
r |
r |
(4.13) |
dt |
dV = òρRdV + ò pndF , Н |
|||
V |
V |
F |
|
Справа в (4.13) стоит сумма внешних массовых и поверхностных сил, действующих на жидкий объем V и ограничивающую его поверхность F .
Уравнение (4.13) представляет собой основное динамическое соотношение ГГД в интегральной форме. Его часто называют также уравнением импульсов, так как оно преобразуется к последнему простым умножением на dt .
Преобразуем поверхностный интеграл, входящий в уравнение (4.13), в объемный. С учетом (4.3):
é |
cos(n, x) + py |
cos(n, y) + pz |
ù |
(4.14) |
ò pndF =òë px |
cos(n, z)û dF . |
FF
Ввекторном анализе доказывается, что для произвольной функции a (векторной или скалярной) справедливы следующие соотношения:
|
r |
|
¶a |
dV ; |
||
òacos(n, x)dF = ò |
¶x |
|||||
F |
V |
|
|
|
||
|
r |
ò |
¶a |
dV |
||
òa cos(n, y)dF = |
¶y |
|||||
F |
V |
|
||||
|
r |
|
¶a |
dV |
||
òacos(n, z)dF = ò |
¶z |
|||||
F |
V |
|
|
|
Тогда равенство (4.14) с учетом (4.15) запишется в виде
r |
é |
r |
|
¶py |
|
r |
ù |
|
¶px |
|
|
¶pz |
|||||
ò pndF =òê |
|
+ |
|
+ |
|
ú dV . |
||
¶x |
¶y |
¶z |
||||||
F |
V ë |
|
|
û |
(4.15)
(4.16)
Подставляя (4.16) в (4.13), в силу произвольности объема V можем приравнять подынтегральные выражения. Разделив обе части уравнения на ρ ,
получим векторную форму дифференциального уравнения движения в напряжениях:
r |
1 |
æ |
r |
|
|
|
¶py |
|
|
|
r |
|
|
ö |
|
dW |
|
|
||||
¶p |
|
|
|
|
|
¶p |
z |
|
|
|
(4.17) |
|||||||||||
R + |
|
ç |
|
x |
+ |
|
|
|
+ |
|
|
|
÷ = |
|
|
|
||||||
ρ |
|
|
¶y |
¶z |
dt |
|
||||||||||||||||
|
è |
¶x |
|
|
|
|
|
ø |
|
|
|
|||||||||||
В проекциях на оси декартовой прямоугольной системы координат, |
||||||||||||||||||||||
учитывая (4.5), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X + |
1 |
æ |
¶σ |
x + |
|
¶τ yx |
+ |
¶τ |
ö |
= |
¶u |
; |
|
|||||||||
|
|
ç |
|
|
|
|
|
|
|
|
|
zx |
÷ |
|
|
|||||||
|
|
|
|
|
¶y |
|
|
|
|
|
¶t |
|
||||||||||
|
ρ è |
¶x |
|
|
|
|
|
¶z ø |
|
|
~ 6 ~ |
|||||||||||
С.М. Мухаметшин. Лекции по гидрогазодинамике |
Y + ρ1 X + ρ1
æ |
¶τ |
xy |
+ |
¶σ |
y |
+ |
¶τ |
zy |
ö |
= |
|||
ç |
|
|
|
|
÷ |
||||||||
|
|
|
¶y |
¶z |
|
||||||||
è |
¶x |
|
|
|
|
ø |
|
||||||
æ |
¶τ |
xz + |
|
¶τ yz |
|
+ |
¶σ |
|
ö |
= |
|||
ç |
|
|
|
|
|
|
|
z ÷ |
|||||
|
|
¶y |
|
¶z |
|||||||||
è |
¶x |
|
|
|
|
ø |
|
¶υ |
; |
(4.18) |
¶t |
|
|
¶¶wt .
Обобщенная гипотеза Ньютона о связи между напряжениями и скоростями деформаций (закон Стокса).
Обобщение гипотезы (закона) Ньютона (1.7) для жидкого объема приводит к закону трения Стокса:
напряжения пропорциональны скоростям деформаций, причем касательные – скорости угловой деформации, а нормальные, возникающие от действия сил вязкости, соответствуют скорости линейной деформации.
При этом для касательных напряжений, учитывая (4.7), можно показать, что
τ xy |
=τ yx = 2μθz |
|
æ |
¶υ |
|
+ |
¶u |
|
ö |
|
|||||||
= μ ç |
|
|
÷ |
|
|||||||||||||
|
|
|
|
è |
|
¶x |
|
|
|
|
¶y |
|
ø |
|
|||
τ yz |
=τ zy = 2μθx |
|
æ |
|
¶w |
+ |
¶υ |
ö |
(4.19) |
||||||||
= μ ç |
|
÷ |
|||||||||||||||
|
|
|
|
è |
|
¶y |
|
|
|
|
¶z |
ø |
|
||||
τ zx |
=τ xz = 2μθy |
|
æ |
¶u |
+ |
|
¶w |
ö |
|
||||||||
= μ ç |
¶z |
|
|
¶x |
÷ |
|
|||||||||||
|
|
|
|
è |
|
|
|
|
ø |
|
|||||||
Для нормальных напряжений, возникающих от действия сил вязкости, |
|||||||||||||||||
можно записать: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¢ |
= 2μεx |
= 2μ |
∂u |
; |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||||
|
σ x |
¶x |
|
|
|
|
|||||||||||
|
σ ¢y |
= 2με y |
= 2μ |
|
¶υ |
|
; |
|
|
(4.20) |
|||||||
|
|
¶y |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
¢ |
= 2μεz |
= 2μ |
|
∂w |
. |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
σ z |
|
|
¶z |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Полные нормальные напряжения в вязкой жидкости можно представить в
виде:
σ x |
= − p0 |
′ |
|
+ σ x ; |
|
||
σ y |
= − p0 + σ ′y ; |
(4.21) |
|
σ z |
= − p0 |
′ |
|
+ σ z . |
~ 7 ~ |
||
С.М. Мухаметшин. Лекции по гидрогазодинамике |
где p0 - некоторая скалярная величина, аналогичная давлению в идеальной
жидкости.
В тензорном анализе доказывается, что сумма нормальных напряжений (σ x + σ y + σ z ) инвариантна относительно поворота осей координат, то есть не
зависит от направления, как и давление в идеальной жидкости.
Для того, чтобы обеспечить совпадение свойств давления в вязкой и идеальной жидкости, принято считать, что
p = -σ x +σ y + σ z . (4.22) 3
Тогда, решая систему (4.20) – (4.22), можно найти, что
σ x |
æ ¶u |
- |
1 |
r ö |
; |
|
= - p + 2μ ç |
¶x |
3 |
divW ÷ |
|||
|
è |
|
ø |
|
||
σ y |
æ |
¶υ |
- |
|
r ö |
(4.23) |
= - p + 2μ ç |
1 divW ÷ |
|||||
|
è |
¶y |
|
3 |
ø |
|
σ z |
æ ¶w |
- |
1 |
r ö |
|
|
= - p + 2μ ç |
¶z |
3 |
divW ÷ |
|
||
|
è |
|
ø |
|
В идеальной жидкости, когда μ = 0, из (4.23) следует (4.11).
Уравнения движения Навье-Стокса для вязкой сжимаемой жидкости.
Если в уравнении движения в напряжениях (4.18) подставить зависимости (4.19) и (4.23), то получим дифференциальные уравнения неустановившегося пространственного движения вязкой сжимаемой жидкости, называемые уравнениями Навье-Стокса (1845). Не приводя преобразований, запишем их в проекциях на координатные оси (при μ = const для всего поля течения):
du |
|
|
|
1 ¶p |
|
1 |
¶ |
|
|
r |
||||
dt |
= X - |
|
|
|
¶x |
+νDu + |
3ν |
|
|
|
|
(div W |
||
ρ |
¶x |
|
|
|||||||||||
dυ |
|
|
|
1 |
|
|
¶p |
|
1 |
¶ |
|
|
r |
|
dt |
= Y - |
|
|
|
|
¶y |
+νDυ + |
3ν |
|
|
|
(div W |
||
|
ρ |
|
¶y |
|||||||||||
dw |
|
|
|
1 |
|
|
¶p |
|
1 |
¶ |
|
|
r |
|
dt |
= Z - |
|
|
|
¶z |
+νDw + |
3ν |
|
|
|
(div W |
|||
|
ρ |
|
¶z |
|
В векторной форме (4.24) примет вид:
dW |
r |
1 |
r |
1 |
r |
|
|
= R - |
|
grad p +νDW + |
|
ν grad(divW |
|
dt |
ρ |
3 |
||||
|
|
|
); ); ).
) , Н/кг.
(4.24)
(4.25)
Уравнение (4.25) представляет собой специфическую для вязкой сжимаемой
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 8 ~ |
жидкости форму второго закона Ньютона. В правой части уравнения стоят, отнесенные к единице массы, массовые силы, силы давления и силы трения, а в левой части отнесенные к единице массы силы инерции.
Уравнения движения идеальной жидкости в форме Эйлера и Громеки-Лэмба.
Из уравнений (4.25) и (4.24), полагая ν = 0, получаем дифференциальные
уравнения Эйлера (1755) движения идеальной жидкости в векторной:
и координатной формах:
Согласно (2.16)
dW dt
dWdt = Rr - ρ1 grad p .
dudt = X - ρ1 ¶¶px ; ddtυ = Y - ρ1 ¶¶py ; dwdt = Z - ρ1 ¶¶pz .
( |
r r |
) |
r ¶W |
|
|
W + |
|
||
= |
W ×Ñ |
¶t |
||
|
|
|
|
В векторном анализе доказывается, что
æ |
|
2 ö |
r |
r r |
r |
grad çW |
÷ |
= (W |
×Ñ)W + (W ´ rotW ) |
||
è |
2 |
ø |
|
|
|
Тогда можно записать, что
dW |
æ |
|
2 ö |
r |
r |
¶W |
= gradçW |
÷ |
- (W ´ rotW )+ |
||||
dt |
è |
2 |
ø |
|
|
¶t |
(4.26)
(4.27)
(4.28)
(4.29)
Подставляя (4.29) в (4.26), получим уравнение движения идеальной жидкости в форме Громеки-Лэмба (1881 г.):
r |
1 |
æ |
|
2 ö |
r |
r |
¶W |
|
|
R - |
grad p = gradçW |
÷ |
- (W ´ rotW )+ |
(4.30) |
|||||
ρ |
¶t |
||||||||
|
è |
2 |
ø |
|
|
|
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 9 ~ |
В координатной форме будем иметь
|
|
|
1 ¶p |
|
|
¶ |
|
æ |
|
2 |
|
X - |
|
= |
|
|
çW |
|
|||||
|
|
|
|
|
|
||||||
|
ρ ¶x |
|
|
|
|||||||
|
|
|
|
|
¶x è |
2 |
|||||
|
|
|
1 ¶p |
|
|
¶ |
|
æ |
|
2 |
|
Y - |
|
= |
|
|
çW |
|
|||||
|
|
|
|
|
|
||||||
ρ ¶y |
|
|
|
||||||||
|
|
|
|
¶y è |
2 |
||||||
|
|
1 ¶p |
|
|
¶ |
|
æ |
|
2 |
||
Z - |
|
= |
|
|
çW |
|
|||||
|
|
|
|
|
|
|
|
||||
|
ρ ¶z |
|
|
|
|
||||||
|
|
|
|
¶z è |
2 |
|
ö |
r |
|
÷ |
- (W ´ |
|
ø |
|
|
ö |
r |
´ |
÷ |
- (W |
|
ø |
|
|
ö |
r |
´ |
÷ |
- (W |
|
ø |
|
|
r |
+ ¶¶ut |
|
|
rotW )x |
; |
|
|
r |
+ ¶¶υt |
|
|
rotW )y |
; |
(4.31) |
r ) + ¶w rotW z ¶t .
Основное достоинство этой формы уравнений движения в том, что здесь выделен в явной форме вихрь скорости, характеризующий особенности движения жидкости как деформируемой среды.
Интегралы уравнений движения идеальной жидкости |
|
|
Пусть массовые силы имеют потенциал, |
то есть R = -gradF . |
Введем |
некоторую функцию давления P(x, y, z) такую, что |
|
|
dP = dp , то есть P = ò dp , |
(4.32) |
|
ρ |
ρ |
|
Тогда, например, ¶¶Px = ρ1 ¶¶px . При этом
|
1 |
æ |
|
|
|
ö |
|
|
|
|
|
dp = |
ç |
¶p dx + |
¶p dy + |
¶p dz ÷ |
= |
¶P dx + |
¶P dy + |
¶P dz = dP , |
|||
|
|||||||||||
ρ ρ è |
¶x |
¶y |
¶z |
ø |
|
¶x |
¶y |
¶z |
то есть получаем соотношение (4.32). |
|
|
|
|
В векторном виде имеем, что |
1 |
|
|
|
grad P = |
grad p . |
(4.33) |
||
|
||||
|
ρ |
|
Тогда уравнение Громеки-Лэмба (4.30) можно представить в виде
æ |
|
2 ö |
r |
r |
¶W |
= 0 |
gradçF + P + W |
÷ |
- (W ´ rotW )+ |
||||
è |
2 |
ø |
|
|
¶t |
(4.34) |
Рассмотрим случай неустановившегося безвихревого течения, когда
С.М. Мухаметшин. Лекции по гидрогазодинамике |
~ 10 ~ |