Математический аппарат, используемый в механике жидкости

МИНИСТЕРСТВО ОБРАЗОВАНИЯ УКРАИНЫ

ОДЕССКИЙ ГОСУДАРСТВЕННЫЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

Иванов К. Ф., Сурков С. В.

МЕХАНИКА ЖИДКОСТИ И ГАЗА

КОНСПЕКТ ЛЕКЦИЙ ДЛЯ СТУДЕНТОВ МЕХАНИЧЕСКИХ И ЭНЕРГЕТИЧЕСКИХ СПЕЦИАЛЬНОСТЕЙ

ЧАСТЬ 1

ОДЕССА ОГПУ 1995

УДК 62-82(075)

Иванов К.Ф., Сурков С. В. Механика жидкости и газа. Конспект лекций для студентов механических и энергетических специальностей. Часть 1. - Одесса: ОГПУ, 1995. - 119 с.

Рецензенты:

ISBN 5-7763-1624-3

1. МАТЕМАТИЧЕСКИЙ АППАРАТ, ИСПОЛЬЗУЕМЫЙ В МЕХАНИКЕ ЖИДКОСТИ

В науке нет другого способа приобретения, как в поте лица: ни порывы, ни фантазии, ни стремление всем сердцем не заменяет труда.

А. И. Герцен

Изучение механики жидкости и газа, понимание сущности рассматриваемых физических явлений и процессов тесно связано с усвоением достаточно развитого математического аппарата, которым эта наука оперирует. Принципиально гидромеханика может излагаться как на базе векторного, так и координатного методов. Вопрос о том, какому из них отдать предпочтение, с давних пор служил источником дискуссий. Так, например, известный физик Уильям Томсон (лорд Кельвин) считал, что «векторы сберегают мел и расходуют мозг». Противником использования аппарата векторного анализа являлся и академик А.Н.Крылов, приводивший достаточно веские аргументы против его применения. Тем не менее векторное построение курса находит широчайшее применение. Одной из причин этого является общая тенденция к сокращению времени, отводимого на изучение дисциплины. В настоящем пособии не отдается решающее предпочтение ни одному из этих методов, они используются по мере необходимости с учетом конкретной ситуации и стремления наиболее простым и доступным способом донести до изучающего содержание вопроса.

Ниже приводятся некоторые, необходимые для понимания дальнейшего, сведения из векторного анализа и теория поля, в основном известные студентам из курса математики. Разумеется, что в рамках пособия они не могут претендовать на достаточную глубину и широту и носят рецептурный характер. Желающим основательно углубить свои знания в этой области можно рекомендовать книгу: Кальницкий Л.А., Добротин Д.А., Жевержеев В.Ф. Специальный курс высшей математики для втузов. М.: Высшая школа, 1976. - 389с.

Одной из важнейших особенностей механики жидкости является то, что в основу ее положена так называемая модель сплошной среды. Как известно, для описания среды, состоящей из большого числа молекул в сравнительно малом объеме (жидкости и газы) в физике широко используются два пути: феноменологический и статистический (иногда их называют корпускулярной и континуальной моделями). Феноменологический путь изучения основывается на простейших допущениях. Оставляя в стороне вопрос о строении вещества, он наделяет его такими свойствами, которые наилучшим образом устанавливают соответствие между наблюдаемыми явлениями и их описанием.

При таком подходе жидкости (газы) рассматриваются как непрерывная среда, способная делиться до бесконечности. Другими словами, жидкость (газ) представляется состоящими из достаточно малых частиц непрерывным образом заполняющих пространство. Эта среда обладает свойством инерции и наделена различными физическими свойствами. В соответствии с такой моделью все параметры жидкости (плотность, вязкость и др.) изменяются непрерывно от точки к точке, что позволяет при анализе движения среды применять математический аппарат дифференциального и интегрального исчислений хорошо разработанный для непрерывных функций.

Понятие о частицах жидкости, которым широко оперирует механика жидкости и газа, неразрывно связано с понятием о физически бесконечно малом объеме. Это объем, размеры которого пренебрежимо малы по сравнению с характерными размерами объекта, но он содержит в себе настолько много молекул, что его средние характеристики (например, плотность) становятся устойчивыми по отношению к изменению объема. Поэтому, например, фраза «объем стягивается в точку» означает, что он стремится не к нулю, а к физически бесконечно малому объему. Следует твердо усвоить, что все законы механики жидкости справедливы до тех пор, пока справедлива модель сплошной среды. Количественно это можно оценить по величине числа Кнудсена, представляющего отношение длины свободного пробега молекул l к характерному размеру течения L, т.е.

Kn = L    (1.1

Принято считать, что законы механики жидкости справедливы, если Kn < 0,01.

1.1. Векторы и операции над ними.

Полем какой-либо величины называется пространство, в каждой точке которого эта величина вполне определена. Если эта величина скаляр, т.е. характеризуется одним числом, то поле называют скалярным (поле плотности, поле температуры).

Векторным называется поле, которое характеризуется в каждой точке пространства величиной и направлением. К этому следует лишь добавить, что непременным условием, связанным с векторными величинами, является то, что они должны складываться по правилу параллелограмма. Поэтому, например, поток автомашин, движущихся по улице и характеризующийся как величиной, так и направлением не является вектором.

Единичные векторы (орты) в декартовой системе координат будем обозначать ex, ey, ez. Тогда вектор U может быть представлен

как

где а- угол между векторами.

Ясн о, что скалярное произведение обращается в нуль, если векторы u и v взаимно перпендикулярны.

1. Векторное произведение двух векторов.

В противоположность скалярному произведению, здесь первое слово указывает на то, что результат действия есть вектор.

Векторное произведение может быть записано в виде определителя третьего порядка

ez

Uz

Vz


(1.4)


r

r

r

r

ex

ey

U X v =

Ux

uy

vx

vy

Раскрывая определитель по общим правилам, получаем:

U X v = Зс (uyvz - uzvy) - ey (uxvz - uzvx ) + e- (uxvy - uyvx ) 0.5)

1.2. Операции первого порядка (дифференциальные характеристики поля).

В теории поля рассматриваются три так называемые операции первого порядка. Эти операции позволяют, выполнив определенные математические действия превратить

-    скалярную величину в векторную;

-    векторную величину в скалярную;

-    векторную - в другую векторную;

Эти операции соответственно называются: градиент, дивергенция и ротор (вихрь). Рассмотрим каждую из них.

Градиент какой-то скалярной функции р(х, у, z) есть вектор,

образующийся в результате выполнения следующих действий:

.    r др    r др    r др

grad р = е дх    + e дУ    + e д    (16)

Физически градиент есть вектор, в направлении которого функция в данной точке поля изменяется с максимальной скоростью. Дивергенцией вектора u называется выражение вида

..    r дих    dUy    duz

div    u = ^r±    + ^т^ +    ^rz    (17)

дх    ду    дz    v

Следовательно, любое векторное поле дает некоторое скалярное поле, а именно поле своей дивергенции (расходимости). Если div u = 0, то поле называют соленоидальным.

Вихрь поля (ротор) - это вектор, образующийся при выполнении операции


\ду    дz)    \ дz    дх)    V дх    ду)

Если rot u = 0, то поле называют безвихревым.

Каждая из трех операций имеет гидродинамическую интерпретацию, которая приводится в соответствующих разделах курса.

1.3. Операции второго порядка.

Операции grad р, div u, rot u, переводящие скаляр в вектор, вектор в скаляр и вектор в вектор порождают пять операций второго порядка:

-    превращение скалярной величины в векторную

grad (div u);

-    превращение векторной величины в скалярную

div (grad p);

div (rot u);

-    превращение одной векторной величины в другую

rot (grad p);

rot(rot u).

В теории поля показывается, что два из этих пяти соотношений тождественно равны нулю: div (rot u) = 0 и rot(grad p) = 0.


, г г rot u = ех


г


дuz    дuy \    Г ( д uv д u^ 1    г (дuy    дu


х


(1.8)


+ еу


+ ez


Операция div (grad p) носит название оператора Лапласа для скалярного поля и имеет вид

(    ч    д2p    д2p д2 p

div (gradр) = д2 +    + д2

дх    ду    dz


Если рассматривается векторное поле ротора (rot U), то поток этого поля представляется как

jJ(rot и) n dS (1.11)

s

1.4.2. Циркуляция вектора поля.

Пусть рассматривается векторное поле какой-то величины и. Циркуляцией вектора и вдоль контура L называют криволинейный интеграл вида

r dl


1.4. Интегральные соотношения теории поля.

1.4.1. Поток векторного поля.

Пусть dS (рис. 1.1) - элемент поверхности, а n - единичный вектор, направленный по внешней нормали.

П отоком векторного поля (например, и) называют поверхностный интеграл вида

и


(1.10)


(1.12)


Иногда этот интеграл интерпретируется как «работа» векторного поля вдоль контура L. Если циркуляция векторного поля вдоль замкнутого пути (контура) равна нулю, то поле называют потенциальным.

1.4.3. Формула Стокса.

Эта формула позволяет преобразовать криволинейный интеграл вдоль замкнутой пространственной кривой в поверхностный интеграл по поверхности, натянутой на эту кривую, т.е.


(1.13)

L    S

т.е. циркуляция вектора поля вдоль контура равна потоку вихря через поверхность, ограниченную этим контуром.



я


n dS


=fи


j и ¦ dr = jj(rot и) ¦ n dS


1.4.4. Формула Гаусса-Остроградского.

Это соотношение, часто называемое преобразованием Гаусса-Остроградского, связывает поверхностный интеграл по замкнутой поверхности с тройным интегралом по области, ограниченной этой поверхностью

Uи ¦ n dS = JJJdiv и dV

(1.15)


S    V

Формула показывает, что поток векторного поля через замкнутую поверхность равен тройному интегралу от дивергенции поля по объему, ограниченному этой поверхностью.

В механике жидкости широко используется формула, являющаяся следствием формулы Гаусса-Остроградского для скалярного поля

dS = JJJgrad (р dV    (1.16)

S    V

где р - какая-то скалярная функция.

2. ОСНОВНЫЕ ФИЗИЧЕСКИЕ СВОЙСТВА И ПАРАМЕТРЫ ЖИДКОСТИ. СИЛЫ И НАПРЯЖЕНИЯ.

Сначала это покажется сложным, но сначала все сложно.

Мусаши.

Физические свойства и параметры, характеризующие жидкость, достаточно полно изучаются в курсе физики. Поэтому в настоящем пособии рассматриваются лишь те из них, которые непосредственно связаны с явлениями и процессами, типичными для гидромеханики.

2.1. Плотность.

Под средней плотностью, либо, что то же, плотностью физически бесконечно малого объема, понимают частное от деления его массы на объем, т.е.

р=M    м>

Плотность выражается в кг/м3.

В литературе часто оперируют понятием удельного веса, т.е. частного от деления веса частицы на ее объем

G gM

Y = — = ——    (2.2)

Г V V

Как следует из (2.2), удельный вес выражается в Н/м3. Заменяя в (2.2) M/V его значением из (2.1), получаем связь между плотностью и удельным весом:

Y = pg    (2.3)

Таким образом, в международной системе (СИ) плотность воды при t =C р= 1000 кг/м3, а ее удельный вес у = 9800 Н/м3.

2.2. Вязкость.

Под вязкостью понимают свойство жидкости оказывать сопротивление перемещению ее частиц. Физической причиной вязкости является молекулярное взаимодействие. Вследствие различия в молекулярной структуре капельных жидкостей и газов различна и природа их вязкостей. В жидкостях вязкость есть проявление сил сцепления между молекулами, в газах она -результат взаимодействия, обусловленный хаотическим движением молекул. Поэтому при повышении температуры в газах вязкость увеличивается за счет более интенсивного движения молекул. Наоборот, в капельных жидкостях повышение температуры приводит к снижению вязкости, т.к. происходит увеличение среднего расстояния между молекулами.

Равновесное состояние вещества характеризуется распределением его параметров в пространстве. Если за счет какого-либо воздействия окажется, что в каком-то месте пространства возникла неравновесность, то в веществе начинает происходить механический или тепловой обмен, который стремится сгладить неравномерность. В общем случае этот обмен называют процессом переноса. В различных явлениях можно наблюдать процессы переноса энергии, массы (вещества) и количества движения. Как будет показано ниже, вязкость обусловлена процессом переноса количества движения.

Для уяснения того, как проявляются силы вязкости, рассмотрим течение жидкости в круглой трубе. Будем считать, что векторы скоростей частиц параллельны оси х. Забегая вперед, отметим, что такое течение существует в природе и носит название ламинарного.


b

a


Пользуясь чисто интуитивными

представлениями, установим вид распределения скоростей в поперечном сечении потока. Сразу же отметим, что графическое изображение распределения скоростей в поперечном сечении называют эпюрой скоростей (либо полем скоростей). Очевидно, что скорости частиц, находящихся на стенках трубы, равны нулю и возрастают по мере приближения к оси (на оси u = u max) как это показано на рис. 2.1.

Рассмотрим два слоя жидкости (a-a и b-b), расположенные на расстоянии dy. Пусть слой a-a движется со скоростью u, тогда, как следует из эпюры, слой b-b имеет скорость u+du. Таким образом, на верхней и нижней гранях прямоугольной жидкой частицы, расположенной между слоями, скорости различны, что в соответствии с законами механики должно привести к ее деформации. Заметим, что такое движение в гидромеханике называют простым сдвигом, либо течением чистого сдвига.

Взаимодействие молекул через этот элемент приводит к появлению касательной составляющей напряжения. При этом знак этой составляющей, т.е. ее направление, таково, что оно соответствует уменьшению разности скоростей по обе стороны рассматриваемого элемента. Величина силы трения, возникающая между слоями движущейся жидкости, определяется по формуле, предложенной Ньютоном и подтвержденной многочисленными и тщательно поставленными опытами нашего соотечественника профессора Н.П.Петрова. Эта формула имеет вид:

du dy

где S - площадь поверхности соприкасающихся слоев;

- - динамическая вязкость, зависящая от физической природы жидкости, ее агрегатного состояния и температуры, и практически не зависящая от давления. Динамическая вязкость выражается в Пас.

В технических приложениях часто используется не динамическая, а кинематическая вязкость, представляющая собой отношение

v =-    (2.5)

р

Кинематическая вязкость выражается в м2/с.

du

Величина — характеризует изменение скорости в dy

направлении нормали к ней, либо, если говорить об эпюре - темп изменения скорости. Иногда эту величину называют поперечным градиентом скорости.

Разделим правую и левую части (2.4) на S. Отношение Ртр / S

есть не что иное, как касательное напряжение т, т.е.

т-=Л    (2.6)

dy

Таким образом, можно сказать, что вязкость жидкости - это способность ее оказывать сопротивление касательным напряжениям.

Из (2.6) можно сделать еще один важный вывод. Если жидкость находится в состоянии покоя, то и = 0 и, следовательно, т= 0, т.е. в покоящейся жидкости силы вязкости не проявляются. Это согласуется и с обычными житейскими представлениями. Действительно, для того, чтобы ответить на вопрос о том, является ли вязкой среда, налитая в сосуд, например, стакан, стоящий на столе, необходимо либо попытаться перелить ее в другой сосуд, либо, обмакнув в нее какой-то предмет, посмотреть как она стекает с него. Смысл этих действий в том, что мы интуитивно чувствуем, что требуется наблюдать движение этой среды.

Выше было высказано предположение, что вязкость обусловлена переносом количества движения. Для того, чтобы убедиться в этом, рассмотрим формулу Ньютона с позиций физических величин, входящих в нее

му


кг


Н


кг • м


м с2м

В числителе - количество    движения, т.е.    т -    это количество движения,

переносимое через единицу поверхности в    единицу времени.

И, наконец,

A dl    установим    физический

-I Ы    U+dU    смысл поперечного

о о    о        градиента скорости, для

чего рассмотрим жидкую частицу, показанную на рис. 2.2. Вследствие разности скоростей на верхней и нижней гранях, первоначально прямоугольная частица будет деформироваться и превращаться в параллелограмм.

Отрезок dl

характеризует величину деформации за время dt, т.е. dl = du • dt,

м 2с

У


Рис. 2.2


X


du dl dl    du    tg    y    ^

тогда — =-, но — = tg y, тогда — =    . Следовательно,

dy dt • dy dy    dy    dt

поперечный градиент скорости представляет собой скорость

относительной деформации сдвига. Таким образом, касательное

напряжение в жидкости линейно зависит от скорости относительной

деформации. В этом принципиальное отличие жидкости от твердого

тела, в котором касательные напряжения зависят от величины

деформации, а не от ее скорости.

Жидкости, удовлетворяющие (2.6) называются ньютоновскими, а не подчиняющиеся этой формуле - неньютоновскими. К числу последних относятся растворы полимеров и др.

2.3. Классификация сил.

Как и в механике твердого тела, в гидромеханике силы классифицируются по разным признакам: внутренние и внешние, сосредоточенные и распределенные.

Очевидно, что в механике жидкости могут рассматриваться лишь распределенные силы, не вызывающие деформации жидкого тела. При этом они должны быть внешними по отношению к объекту. Перевод внутренних сил в категорию внешних производится известным методом (метод сечений, либо метод «замораживания»), суть которого сводится к тому, что в среде выделяется («замораживается») замкнутый объем, внешняя среда мысленно отбрасывается и ее действие заменяется действием распределенных сил. Важнейшей особенностью гидромеханики как науки является то, что в ней, помимо приведенной выше классификации, силы разделяются на массовые и поверхностные.

2.3.1. Массовые силы.

Массовыми называют силы, величина которых пропорциональна массе рассматриваемого объема. Важнейшей особенностью является то, что они действуют на все частицы жидкости. В общем случае это силы, подчиняющиеся второму закону Ньютона F = та. В проекциях на декартовы оси координат можно записать: Fx = max; Fy = may; Fz = maz. В гидромеханике вместо

ax, ay, az принято писать X, Y, Z. Поделив обе части записанных

F    F    F

выражений на массу, получим    —    = X;    —    =    Y;    — = Z.

m    m    m

Таким образом, X,    Y и    Z есть проекции    единичных массовых

сил на соответствующие координатные оси, иногда их называют напряжениями массовых сил. Если в жидкости выделить элементарный объем dV, то его масса - pdV. В общем случае массовая сила, действующая на этот объем pFdV, а главный

вектор массовых сил, действующих на весь объем, представляется как

Hjphv    (2.7)

v

2.3.2. Поверхностные силы.

В отличие от массовых, поверхностные силы действуют лишь на частицы, находящиеся на поверхности жидкого объема.

Выделим на поверхности жидкого объема элементарную площадку A S, ориентация этой площадки в пространстве задается внешней нормалью n. Обозначим через Apn поверхностную силу,

A1 r

приложенную к площадке A S. Предел отношения lim —n = pn

as^0 AS

называют напряжением поверхностной силы.

Таким образом, первое, что необходимо усвоить при рассмотрении этого вопроса - это то, что под действием внешних сил в жидкости возникают напряжения. И второе по порядку, но не менее важное по существу.


В общем случае pn не является обычным вектором. Его величина зависит от ориентации площадки в пространстве. Это означает, что если через данную точку пространства провести одинаковые по величине, но

Рис. 2.3    различно ориентированные

площадки, то действующие на них напряжения поверхностных сил

будут различны.

Физическая величина, характеризуемая в данной точке вектором pn, принимающим бесконечное множество значений в зависимости от ориентации площадки, называется тензором напряжений.

Таким образом, на площадку dS действует поверхностная сила pndS, а на всю поверхность, ограничивающую объем V

JJ PndS    (2.8)

S

Проекция pn на направление нормали называется нормальным напряжением, а проекция на площадку действия

- касательным напряжением.

2.3.3. Тензор напряжения.

Для уяснения дальнейшего необходимо подробней рассмотреть вектор


pn.

В движущейся среде мысленно выделим частицу в форме жидкого тетраэдра.

Пусть n - внешняя нормаль к четвертой (наклонной) грани тетраэдра , а площадь этой грани dS (см. рис. 2.4).

Площади других граней - соответственно dSX, dSy, dSz, т.к. их

можно рассматривать как проекции грани ABC на координатные оси. Следовательно, dSX = dS cos(n, X) = nXdS, где nX обозначает направляющий косинус. Аналогично, dSy = nydS, dSz = nzdS.

Обозначим объем тетраэдра dV, тогда действующая на него

¦

массовая сила pFdV, а массовая сила инерции p'adV, где a вектор ускорения жидкого тетраэдра. Поверхностная сила, действующая на наклонную грань - pndS. Для трех других граней можем записать:

-p xdSx = -p xnxdS

-pydSy = -pynydS

-pzdSz = -pznzdS    r r

Рис. 2.4    Знаки минус, т.к. векторы px, py

и pz направлены в стороны, противоположные координатным осям.

Запишем уравнение движения тетраэдра, которое в соответствии с общими законами механики должно иметь вид:

Масса ¦ ускорение = (результирующая массовых сил) +

+ (результирующая поверхностных сил).

И меем:

p^dV = pFdV + pndS - pXnXdS - pynydS - pznzdS

¦

Слагаемые padV и pFdV есть величины третьего порядка малости, а остальные - второго, поэтому ими можно пренебречь, что дает

pn = nxpx + nypy + nzpz    (2.9)

Из этого равенства следует, что напряжение pn при произвольной ориентации нормали n может быть определено, если известны

напряжения в той же точке для площадок, внешние нормали которых параллельны осям Ox, Оу и Oz.

kPzz


Проекции векторов px, ру и pz на

Pzy

Рху,

/

у

Рис. 2.5


Рхх


координатные оси х, у, z обозначаются:

х


Рхх Рху Pxz Р ух Р уу Р yz

Pzx Pzy Pzz Первый подстрочный индекс указывает ось,

перпендикулярную ориентации площадки, второй - ось, на которую спроектировано напряжение.

Для уяснения ориентации рассмотрим параллелепипед, выделенный в движущейся жидкости и показанный на рис. 2.5.

Из рисунка, в частности, видно, что напряжения с одинаковыми индексами являются нормальными, а с разными - касательными. В проекциях на декартовы оси координат выражение (2.9) может быть записано как

Рпх =    ПхРхх    +    ПуРух    + nzPzx

Рпу =    П хРху    +    ПуРуу    + nzPzy    (2. 10)

Pnz =    nxPxz    +    nyPyz    + nzPzz

Совокупность этих девяти составляющих компонентов напряжения образует тензор напряжения. В матричной форме он записывается в следующем виде:

Рхх

Рух

Pzx

п =

Рху

Руу

Pzy

Pxz

Pyz

Pzz

В тензорном анализе доказывается, что тензор напряжений является симметричным. Это означает, что величины, расположенные симметрично главной диагонали, равны (рух = рху;

Pxz = Pzx; Pzy = Pyz). Следовательно, для определения тензора

напряжений достаточно знать не девять, а шесть скалярных величин.

Следует учесть одно обстоятельство. Векторы напряжений px, py, pz в соотношении (2.9), носящем имя Коши, и приложенные к

координатным площадкам, не имеют объективного физического смысла, т.к. зависят от выбора системы координат. Поэтому такие величины причисляются к так называемым «квазивекторам», хотя к ним и можно применять все операции, применимые к физическим векторам.

К понятию тензора можно подойти и другим путем, который, возможно, покажется более простым. Поэтому целесообразно хотя бы кратко остановиться на нем. Для наглядности тензор можно представить как какой - то оператор , с помощью которого можно преобразовывать векторы в векторы. Упрощая и сводя математический аппарат к механическому, оператор можно представить как какую-то «машину», которая по определенным правилам перерабатывает вводимые в нее векторы. Зная принцип работы этой «машины», можем знать и вектор, который появляется на выходе. Можно записать

I    1    1

r    a    =    A    •    B

где a- входной вектор;

BI - выходной вектор;

A - оператор, который и называют тензором.

Существенное ограничение заключается в том, что оператор должен быть линейным. Определить тензор - это значит задать правила, по которым работает оператор. Для интересующихся таким подходом можно рекомендовать книгу Астарита Дж., Марручи Дж. Основы гидромеханики неньютоновских жидкостей. - М.: Мир, 1978.307с.

И в заключение еще несколько замечаний. Выше уже отмечалось, что одно из фундаментальных свойств жидкости - ее вязкость - не проявляется, если она находится в состоянии равновесия, т.е. в этом случае касательные компоненты тензора равны нулю и действуют лишь нормальные p xx, pyy, pzz,

ориентированные по внешним нормалям (см. рис. 2.5). При этом ясно, что они являются растягивающими напряжениями. Как показывает опыт, в отличие от твердого тела, которое может воспринимать как растягивающие (положительные нормальные напряжения), так и сжимающие (отрицательные нормальные напряжения) напряжения без разрыва сплошности, жидкое тело способно воспринимать лишь сжимающие усилия. Можно показать, что при отсутствии касательных напряжений pxx = pyy = pzz, из чего

следует, что нормальные напряжения в данной точке не зависят от ориентации площадки. Величины, численно равные нормальным напряжениям, но взятые с противоположным знаком, в

гидромеханике называют давлениями, либо более полно -гидростатическими давлениями. Гидростатическое давление обозначают буквой p, т.е.

Р = - Pxx = - Pyy = - Pzz Таким образом, гидростатическое давление, являясь скалярной величиной (как компонента тензора) не зависит от ориентации площадки, на которую оно действует.

Теоретическое изучение движения жидкости связано с так

w    w    n    w

называемой моделью идеальной жидкости. В этой модели жидкость рассматривается как абсолютно несжимаемая среда, неспособная сопротивляться разрывающим усилиям и обладающая абсолютной подвижностью, т.е. лишенная вязкости. Последнее исключает возникновение в ней касательных напряжений.

2.3.4. Уравнение движения в напряжениях.

Получим наиболее общее уравнение, связывающее поверхностные и массовые силы - так называемое уравнение движения в напряжениях. Для вывода уравнения проанализируем движение жидкой частицы, масса которой р dV и поверхность dS. Аналогично тому, как это было сделано для тетраэдра, можем записать уравнение движения в виде

pdLdV - pF dV + p ndS    (2.11)

Для всего движущегося объема (V), поверхность которого S, имеем

pn = nxpx + nypy + nzpz

где nx, ny, nz - направляющие косинусы.

Воспользуемся известными из векторного анализа и справедливыми для любых векторов формулами:

Я n.R dS=!0 dV

SnyR dS -Ш% dV ¦    <213>

dS - inf dV

S    V

Применяя эти формулы к тензору pn, получаем:

+


dp


y + dp z


flp -» ¦ JJJ


dV


(2.14)



du r

- PF -


'y + dp


ffl


dV = 0


z


+


P


)


dt


dx    dy dz


v


V


Но так как dV Ф 0, а объем V выбран произвольно, то

.i    * г ^l zL    ni \


dpx . dp


du i 1

= F + —


'y + dp


z


(2.15)


+


dt


dx    dy dz


Р


v


Это и есть уравнение движения в напряжениях.

В проекциях на декартовы оси координат можем записать:

г


дт


dp


xx


zx


+


)


v dx f дт


Р

1

Р

1


\


дт


xy


zy


= Y +


(2.16)


+


dx


dy

дт


v


дт


= Z +


xz


zz


+


pV dx    dy    dz )

Эта система включает в качестве неизвестных девять величин: три проекции скорости и шесть проекций напряжений. Проекции единичных массовых сил, как правило, известны из постановки задачи.


dux    1

—^ = X + — dt

dUy

dt

duz dt


dy    dz


yz + dP


dz )

Л


dp


дт


yy


yx


+


+


3. ГИДРОСТАТИКА.

Идите, идите вперед, уверенность прийдет к вам позже.

Д'Аламбер.

Гидростатика занимается изучением жидкости, находящейся в состоянии относительного покоя. Под относительным покоем понимают состояние, при котором отсутствуют перемещения частиц относительно друг друга.

В основу гидростатики положены две теоремы: равенство нулю суммы всех сил, приложенных к рассматриваемому элементу жидкости и, как следствие, равенство нулю суммы моментов этих сил относительно какой-то оси. Однако, несмотря на простоту принципов, гидростатика приводит к важным результатам и выводам.

3.1. Уравнение равновесия жидкости.

Уравнения равновесия жидкости могут быть получены из уравнений движения в напряжениях (2.16), если положить в них ux = uy = uz = 0. Кроме того, как было показано, в покоящейся

жидкости касательные напряжения не проявляются, т.е. все производные по t равны нулю. И, наконец, нормальные напряжения заменяем давлением, что дает

X-Р дР = 0;

p дх

Y--др = 0;    (3.1)

Р ду    '    '

Z-Р = 0

p dz

В векторной форме эта система может быть записана в форме

r1

F -pgrad р = 0    (3.2)

Уравнения (3.1) носят название системы дифференциальных уравнений Эйлера для гидростатики. Эта система уравнений показывает, что существует непосредственная связь между величиной гидростатического давления в точке и ее координатами. Эта связь может быть раскрыта, если проинтегрировать (3.1).

На жидкое тело могут действовать силы, имеющие различную физическую природу. Поэтому правомерна такая постановка вопроса: всегда ли под действием приложенных сил жидкость может находиться в состоянии равновесия? Для ответа на этот вопрос необходимо выполнить некоторые преобразования системы дифференциальных уравнений (3.1).

3.2. Основное уравнение гидростатики в дифференциальной форме.

Умножим каждое из уравнений, входящих в (3.1) на dx, dy и dz соответственно и просуммируем их, что даст

(Xdx + Ydy + Zdz)-1 |dp-dx +dpdy + dpdz\ = 0    (3.3)

p\dx dy dz )

Выражение, стоящее в скобках во втором члене уравнения, есть не что иное, как полный дифференциал давления - dp, поэтому можем записать

dp = p( Xdx + Ydy + Zdz)    (3.4)

Это уравнение называют основным уравнением гидростатики в дифференциальной форме. В левой части его - полный

дифференциал, поэтому и правая часть также должна быть полным дифференциалом. Следовательно, силы и плотность должны быть такими функциями х, у и z, чтобы они обращали правую часть (3.4) в полный дифференциал. Если этого не происходит, то равновесие жидкости невозможно. Другими словами, если жидкость находится в состоянии равновесия, то правая часть (3.4) является полным дифференциалом какой-то функции Ф.

Считая плотность постоянной (p= const), можем записать

Xdx + Ydy + Zdz = dФ    (3.5)

Из теоретической механики известно, что скалярное произведение силы на элементарное перемещение частицы называют элементарной работой, т.е.

fxdx + fydy + fzdz    (3.6)

Силы, работа которых не зависит от пути движения, а только от начального и конечного положений, называют потенциальными. При этом для того, чтобы работа силы не зависела от пути движения, необходимо и достаточно, чтобы выражение для элементарной работы, т.е. (3.6), было полным дифференциалом некоторой скалярной функции P, называемой силовой. Взятая с противоположным знаком, она называется потенциалом. Таким образом, рассмотренную выше функцию можно назвать силовой функцией, а (3.4) представить как

dp = pdФ    (3.7)

Из чего следует, что несжимаемая жидкость может находиться в равновесии только под действием сил, имеющих потенциал.

3.3. Эквипотенциальные поверхности и поверхности равного давления.

Поверхности, в каждой точке которых Ф = const, называют эквипотенциальными. Частным случаем эквипотенциальной поверхности является поверхность равного давления, т.е. поверхность, в каждой точке которой р = const. В этом случае dp = 0 и (3.4) принимает вид

p(Xdx + Ydy + Zdz) = 0 Но плотность p^ 0, и, следовательно,

Xdx + Ydy + Zdz = 0    (3.8)

Уравнение (3.8) называют уравнением поверхности равного давления. Если из массовых сил на жидкость действует только сила тяжести, то X = Y = 0; Z = -g (знак минус, т.к. сила тяжести

ориентирована в сторону, противоположную оси z); - gdz = 0 и z = const, т.е. в покоящейся жидкости любая горизонтальная плоскость есть поверхность равного давления.

3.4. Равновесие однородной несжимаемой жидкости в поле сил тяжести. Закон Паскаля. Гидростатический закон распределения давления.

Проинтегрируем основное уравнение гидростатики (3.4) в предположении, что р— const (жидкость несжимаема) и считая, что из    массовых сил    действует только сила тяжести. Как показано    выше,

в этом случае    X    —    Y — О, Z —    -д, т.е. dp — -pgdz, и после

интегрирования

p — -pgz + C    (3.9)

Рис. 3.1


где C - произвольная постоянная. Для ее нахождения используем следующее граничное условие (см. рис. 3.1): при z — z0 p — po. Из (3.9)

следует, что

C — p О + pgzo

И после подстановки

p — po + pg(z> - z)    (3.10)

Как видно из рис. 3.1, разность    (z0 - z)

- глубина погружения рассматриваемой частицы, которую будем обозначать буквой h, т.е.

p — po + pgh    (3.11)

Полученное уравнение выражает известный из курса физики закон Паскаля: давление, приложенное к свободной поверхности, передается во все точки без изменения.

Поскольку любое правильное физическое уравнение должно быть размерностно однородным, то ясно, что член pgh должен выражаться в единицах давления, т.е. в паскалях (Па - Н/м2). Эту величину называют избыточным давлением. Она может быть как положительной, так и отрицательной. Такая трактовка приводит нас к понятию абсолютного давления, которое в соответствии с (3.11) может быть представлено как сумма барометрического (атмосферного) давления и избыточного, т.е.

pабс. — pат. ± pизб.    (3.12)

Отрицательное избыточное давление называют вакуумом.

Вернемся вновь к уравнению (3.10). После деления обеих его частей на pg получаем

pg

pg


В таком виде все его члены выражаются в единицах длины и носят название напоров. Величина z характеризует положение жидкой частицы над произвольно выбираемой горизонтальной плоскостью

P

пьезометрический


отсчета, т.е. z - это геометрический напор;

рд

напор. Сумму этих величин z + Р- называют гидростатическим

р—

напором. Чтобы уяснить физический смысл этих величин, рассмотрим простую схему, показанную на рис. 3.2.

Представим

Рис. 3.2


герметично закрытый сосуд, заполненный жидкостью, находящейся под давлением. Выберем в этом сосуде две произвольно

расположенные точки A и B и, опять-таки произвольно,

горизонтальную плоскость O-O, которую назовем плоскостью отсчета. Координаты частиц, расположенных в точках A и B будут zA и

zB. В соответствии со сказанным выше, величины zA и zB выражают геометрический напор. Введем теперь через крышку сосуда в точки A и B сообщенные с атмосферой стеклянные трубки. Эти трубки называют пьезометрами. Поскольку по условию жидкость находится под давлением, то она начнет подниматься по пьезометрам. Не представляет труда и ответ на вопрос о том, когда прекратится подъем. Очевидно, что это произойдет в тот момент, когда высота столба жидкости уравновесит давление в рассматриваемой точке. Это и есть пьезометрическая высота, либо пьезометрический напор.

Соотношение (3.13) справедливо для любых произвольно выбранных частиц покоящейся жидкости, поэтому в общем виде его

можно записать как z + = const, т.е. для любых точек жидкости

р-

гидростатический напор одинаков. Следовательно, уровни в пьезометрах установятся на одной и той же высоте (плоскость C-C

на рис. 3.2). Уравнение (3.13) выражает так называемый гидростатический закон распределения давления.

3.5. Определение силы давления жидкости на поверхности тел.

Задача сводится к нахождению силы давления жидкости на поверхности стенок, ограничивающих ее.

Рассмотрим криволинейную поверхность АВ произвольной формы, площадь которой S (рис. 3.3). Выделим на ней элементарную площадку dS, пусть h - орт внешней нормали. Сила, действующая на эту площадку

dF = phdS

Рис. 3.3


где р - гидростатическое давление в центре площадки. Обычно в технических приложениях интерес представляет лишь сила, возникающая от избыточного давления. Имея в виду, что р = pgh, получаем

dF = pghhdS    (3.14)

На всю площадь действует сила

F = JJpghndS    (3.15)

Запишем это выражение в проекциях на оси координат, что дает

= pgjjh cos(h,x)dS    (3.16)

Fz = pg^hcos(h,z)dS

(3.17)


Для удобства изобразим отдельно элементарную площадку (см. рис. 3.4). рисунка следует, что


Из


dS • cos(/?,x) = dSe

dScos(h,z) = dS3 где dSe - вертикальная, и dSa -горизонтальная проекции dS. Таким образом

f — pgЯл ¦ d^ (s.18)

Fz — pgjjh • dSs    (S.19)

S

Рис. S.4    Рассмотрим горизонтальную

составляющую.

Из механики известно, что интеграл (3.18) есть статический момент площади, равный произведению Лц т Se, где hц т. -

координата центра тяжести вертикальной проекции.

Следовательно,

Fx — pghц т.Se    (3.20)

т.е. горизонтальная составляющая равна произведению вертикальной проекции стенки на гидростатическое давление в центре тяжести этой проекции.

Определим теперь вертикальную составляющую силы, для чего воспользуемся следствием из формулы Гаусса-Остроградского (см. ф-лу 1.16)

jjpn dS — jjjgradpdV

S    V

I

Из уравнения равновесия (3.2) имеем pF — grad p, т.е.

j]Jgrad pdV — j jjpFdV

V    V

I    _

Вертикальная проекция единичной массовой силы F — Z — g (знак

плюс, т.к. в данном случае ось z ориентирована вниз). Следовательно,

Fz W\pgdV — pg JjJdV — pgV    (S.21)

VV

V носит название объема тела давления. Таким образом, вертикальная составляющая равна весу жидкости, заключенному в

объеме тела давления. Для




нахождения этого объема следует использовать формальное правило: тело давления - это объем, образованный криволинейной стенкой, ее проекцией на свободную поверхность(либо на продолжение свободной поверхности) и вертикальными проектирующими плоскостями. На рис. 3.5 показаны примеры определения тел давлений для двух случаев.

Как следует из рисунка, тело давления может быть как Рис. 3.5    положительным, так и

отрицательным (фиктивным).

3.5.1. Плоская поверхность.

Этот случай можно рассматривать как частный предыдущего, но можно получить и более удобное соотношение. Действительно, общее выражение для силы давления имеет вид (3.15), но так как поверхность плоская, то ориентация нормали для всех ее точек остается одинаковой, и, сл едовательно,

F = р-n hцт S    (3.22)

I

Из формулы (3.22) следует, что F направлена по нормали к стенке, поэтому можно записать

F = Р-h ц т S    (3.23)

Следовательно, сила давления на плоскую поверхность равна произведению ее площади на гидростатическое давление в центре тяжести этой поверхности. Следует отметить, что задачи, связанные с определением сил давления на поверхности, играют исключительно важную роль в гидротехнической практике. Применительно к энергетике и машиностроению круг этих задач заметно сужается и ограничивается, главным образом, расчетом болтовых соединений люков различных резервуаров, находящихся под давлением.


4. КИНЕМАТИКА.

Приобретение любого познания всегда полезно для ума, ибо он сможет впоследствии отвергнуть бесполезное и сохранить хорошее. Ведь ни одну вещь нельзя ни любить, ни ненавидеть, если сначала ее не познать.

Леонардо да Винчи.

Кинематика занимается изучением движения жидкости, не интересуясь причинами, которые его вызвали. По образному выражению Н.Е.Жуковского, кинематика изучает «геометрию движения». Принципиально можно пойти двумя путями. По первому из них изучается движение каждой отдельной жидкой частицы. Чтобы выделить ее, в начальный момент времени to отмечаются ее координаты xo, yo и zo. Движение считается определенным, если в каждый момент времени для каждой частицы известны уравнения, описывающие ее путь во времени, т.е. известны параметрические уравнения траекторий всех частиц. Этот путь предложен Лагранжем. По методу Эйлера изучается изменение скорости и других параметров в точках пространства x, y, z.

В настоящем пособии используется главным образом метод Эйлера. Для желающих глубже разобраться в этом вопросе можно рекомендовать книгу: Федяевский К.К., Войткунский, Фаддеев Ю.И. Гидромеханика. - М.: Судостроение, 1968. - 567 с.

4.1. Установившееся и неустановившееся движения

жидкости.

Установившемся (стационарным) называют движение, при котором основные параметры потока (скорость, давление, плотность) в данной точке пространства не изменяются с течением времени, т.е.

и = f (x, y, z); p = f (x, y, z); p= f (x, y, z)    (4.1)

Если это условие не соблюдается и параметры в точке меняются с

течением времени


и = f (x, y, z, t);

p = f (x, y, z, t);    (4.2)

p= f (x, y, z, t)

движение называют неустановившимся (нестационарным).

В этих формулировках следует обратить внимание на то, что

Рис. 4.1

речь идет о параметрах в точке. Чтобы уяснить это, рассмотрим канал, показанный на рис. 4.1. В гидромеханике такие каналы, в которых площадь сечения уменьшается по ходу потока, называют конфузорами. Исходя из чисто интуитивных представлений ясно, что скорость течения по ходу канала будет возрастать. Возникает вопрос, может ли быть установившемся движение в таком канале. Очевидно, может, если параметры в точках A и B не будут изменяться с течением времени. Определение вида движения не требует, чтобы параметры в точках А, В и С были одинаковы.

4.2. Уравнение неразрывности (сплошности).

Уравнение неразрывности либо сплошности выражает один из фундаментальных законов природы - закон сохранения массы

применительно к жидкой среде.


Рассмотрим объем V, ограниченный поверхностью S (рис.

4.2). Выделим элемент поверхности dS. Пусть n - орт внешней нормали, а и - вектор скорости. Через выделенный элемент dS в единицу времени внутрь объема проникает масса жидкости

-ри • ndS.

(знак минус, т.к. направления и и n противоположны). Секундная масса, проникающая в объем через всю поверхность,

- §Jpu • П dS.

S

С другой стороны, приток жидкости в объем приводит к изменению ее массы. При этом, поскольку выделенный объем является постоянным, то изменение массы может происходить только за счет изменения ее плотности. Скорость изменения массы можно представить как

dt ЯК

V

либо с учетом того, что V = const, можно записать

ШтГ-

V

Очевидно, что изменение массы внутри объема должно быть равно массе, поступившей в него извне, т.е.

ЯШ -Я'dS

V    S

Применяя преобразование Гаусса-Остроградского, получим:

dp

Шм = ~ IfJdiv (pu )dV, лиео

dt

V    V


ш


— + div (p u) dt }

V L

Равенство нулю интеграла возможно лишь при условии

д + div (pru) — 0.    (4.S)

Это и есть уравнение неразрывности. Поскольку при выводе его не делалось никаких ограничений, то оно справедливо как для установившегося, так и для неустановившегося движений сжимаемой и несжимаемой жидкости. Уравнение (4.3) относится к числу фундаментальных уравнений механики жидкости.

Рассмотрим некоторые частные случаи. При установившемся движении все производные по времени равны нулю, что следует из самого определения этого понятия, поэтому

Установим физический смысл этого соотношения. Частные

dux dUy duz производные , dy-,    характеризуют скорость

относительного удлинения (укорочения) жидкой частицы. Если этот процесс происходит одновременно вдоль всех координатных осей, то он приводит к объемному расширению либо сжатию частицы. Ясно, что если частица удлиняется вдоль осей x и у, то она должна укорачиваться относительно оси z. Другими словами, хотя бы одна из производных, входящих в (4.6), должна быть отрицательна, т.к. в противном случае соотношение не может быть равным нулю.

Как уже отмечалось в 1.1, поле, в котором div U — 0, носит название соленоидального.

dV — 0.

4.3. Линии тока и траектории.

Линией тока называется кривая, обладающая тем свойством, что в данный момент времени векторы скоростей в любой ее точке совпадают по направлению с касательными.

В векторной форме это условие может быть записано как и х dS = 0, т.е. векторное произведение должно быть равно нулю. Это, как известно (см. формулу 1.4), может быть записано в виде определителя

ex

ey

ez

Ux

Uy

Uz

= 0

(4.7)

dx

dy

dz

Раскрывая определитель, получаем дифференциальное уравнение линии тока в виде

dx dy dz

(4.8)


и,

и.


Под траекторией понимается след, оставленный движущейся частицей в пространстве. Дифференциальное уравнение траектории

dx dy dz ~    (4.9)

= dt

и

и


и


x " y " z

Из сопоставления (4.8) и (4.9) следует, что в общем случае, т.е. при неустановившемся движении, линии тока и траектории не совпадают.

4.4. Трубка тока (поверхность тока)

В движущейся жидкости наметим бесконечно малый замкнутый контур, и через все точки его периметра проведем линии тока (рис.

4.3).

Рис. 4.3


Образованная таким образом поверхность носит название трубки либо поверхности тока. Ясно также, что поскольку контур намечался в пространстве, занятом движущейся жидкостью, то какая-то часть ее должна находиться и внутри

поверхности тока.

4.5. Струйная модель потока.

Струйная модель потока введена в рассмотрение Л.Эйлером. Основу этой модели составляет понятие о струйке (либо элементарной струйке), под которой понимают жидкость, протекающую внутри трубки тока. Если вспомнить, что границами боковой поверхности трубки тока являются линии тока, т.е. линии, к которым касателен вектор скорости частиц, которые в данный момент времени находятся в ней, то ясно, что ни одна частица не может проникнуть извне в струйку, либо, наоборот, выйти из нее через боковую поверхность. Действительно, вектор скорости частицы, пытающейся, например, проникнуть в струйку извне, должен быть ориентирован к ее границе под каким-то углом, а на самой границе - линии тока - касателен (рис. 4.4).

Из сказанного следует, что струйка ведет себя как трубка с непроницаемыми стенками.


Поперечное сечение струйки мало, поэтому можно допустить, что в пределах сечения все частицы движутся с Рис. 4.4    одинаковыми    скоростями

либо, что то же, эпюра скоростей в сечении представляет собой цилиндр для трехмерной струйки либо прямоугольник - для плоской (двумерной).

На рис. 4.4 показаны эпюры для двух произвольно выбранных сечений плоской струйки. Заметим лишь, что равномерность распределения скоростей в сечении, т.е. движение всех частиц, находящихся в нем, с одной и той же скоростью, вовсе не означает, что в другом сечении эти скорости должны быть такими же, т.е., не обязательно, чтобы tv1 = и2 (см. рис. 4.4).

Совокупность струек, заполняющих поперечное сечение канала конечных размеров, образует поток. Если представить, что соломинка для коктейля - это струйка, то пучок таких соломинок -поток.

4.6. Уравнение неразрывности для струйки.

Первое свойство струйки, говорящее о том, что боковая поверхность непроницаема для частиц, по существу выражает закон сохранения секундной массы. Действительно, если через сечение 1-1 в единицу времени вошла масса dmv то за то же время через сечение 2-2 должна выйти масса dm2, равная dmv Массу жидкости, протекающую через поперечное сечение струйки в единицу времени называют элементарным массовым расходом и обозначают dQm.

Легко убедиться в том, что dQm = ри dA, где dA - площадь поперечного сечения струйки. Действительно, выражая параметры, входящие в это соотношение через единицы физических величин, получаем

кг м 2 кг

— м2 ^ —. м3 с    с

Из сказанного выше следует, что

p1u 1dA 1 — p2 u2dA2    (4.10)

Это и есть уравнение неразрывности для струйки. Если жидкость

несжимаема, т.е. p— const, то p1 — p2 и

u 1dA 1u 2dA 2    (4.11)

При этом произведение udA выражает элементарный объемный расход - dQ.

4.7. Ускорение жидкой частицы.

Запишем выражение для проекции ускорения жидкой частицы на какую-либо координатную ось, например, x. Имеем

я — dux. ax    dt

Для нахождения этой величины следует учесть, что проекция

скорости ux (как и две другие проекции) является функцией

координат x, у, z, которые, в свою очередь, в общем случае зависят

от времени t. Представим величину dux в виде полного

дифференциала

. dux    dux . dux . dux .

dux —^x dt + dx dy + dz

x dt    dx    dy    dz

dx dy

Разделим обе части на dt. Имея в виду, что    — ux,    — uy и


u —— dt    ux,dt

dz

dt = uz, получим

dux    dux    dux    dux    1Л

ax    + ux dx + uy d + uz d    (412)

Аналогичные соотношения можно записать и для двух других компонент.

Выражение (4.12) носит название полной либо субстанциональной производной. Установим смысл величин,

dux

входящих в нее. Производная ^ ~ проекция локального ускорения,

которое характеризует изменение скорости во времени в данной точке пространства. Локальное ускорение обусловлено нестационарностью процесса. Из чего следует, что если движение стационарное (установившееся), то локальное ускорение

отсутствует, т.е. ^ = 0. Три остальных члена (4.12) - проекции

конвективного ускорения, которое возникает при переходе частицы от одной точки пространства к другой, оно обусловлено неравномерностью скоростного поля, т.е. неравномерным распределением скоростей.

4.8. Анализ движения жидкой частицы.

Движение жидкой частицы является более сложным, чем движение твердого тела, которое, как известно из механики, может быть поступательным и вращательным. Особенностью жидкости и ее частиц, как уже неоднократно отмечалось, является легкая деформируемость. Поэтому помимо поступательного и вращательного, жидкая частица может участвовать и в деформационном движении. Это положение и составляет суть так называемой первой теоремы Гельмгольца, к рассмотрению которой мы и приступаем. Оценивая значение работы Г.Гельмгольца, основоположник отечественной аэродинамики Н.Е.Жуковский писал, что «современная гидродинамика своим развитием обязана главным образом Гельмгольцу». Важнейшим достоинством приводимых ниже выкладок и рассуждений является то, что они раскрывают физический смысл и вносят ясность в ряд казалось бы совершенно абстрактных понятий. Выкладки эти достаточно просты, но требуют внимания. Поэтому нужно запастись определенной долей терпения и помнить, что достигаемое понимание сути явлений безусловно оправдает эти затраты труда.

Рассмотрим жидкую частицу в форме прямоугольного параллелепипеда (рис. 4.5). Длина его ребер dx, dy, dz. Деформация такой жидкой частицы может быть как линейной (ребра удлиняются и укорачиваются), так и угловой (грани скашиваются). Удобней рассмотреть каждый из этих видов раздельно. Начнем с

Рис. 4.5


угловых деформаций.

4.8.1. Угловые деформации.

Из рис. 4.5 следует, что угловая деформация (скашивание) может возникнуть из-за разности скоростей, перпендикулярных ребрам (частично этот вопрос уже обсуждался в разделе 2.2). Для

A B

D

A

dx


упрощения целесообразно ограничиться лишь одной гранью, показанной на рис. 4.6.

Пусть компоненты скорости в точке A равны uX, uy, uz. Найдем скорости в точке

B, считая, что движение

установившееся и,

следовательно, все

производные по t равны нулю.

Приращение компоненты

скорости при переходе из одной точки пространства в другую можно

представить как u+du. Так для проекции uX можем записать

Ux + duX, где

.    duX .    duX .

duX =    dx + dy +


с


X


dX

duX


dy


dz


(4.13)


+


dz


Аналогичные выражения можно записать и для других проекций.

Рис. 4.6    Рассмотрим приращение

uX при переходе от точки A к точке B. В этом случае dX = dz = 0, т.е.

duX .

+ ^^ dy


uX = uX + duX = uX


dy


'(A)


(B)


(A)


Предположим, что за время dt за счет разности скоростей в точках A и B ребро займет положение AB'.

Аналогично рассуждая относительно скорости uy в точках A и

D получим:

Точка A: uy (по условию)

duv ,

ТОчка D:    uy( D) = uy( A) +dtdX

За счет разности этих скоростей точка D займет позицию D'. Таким образом

duX

dy


uX - uX    =

dy

dUy

dX

(A)

(B)


dX


U y(D)    U y(A)


Путь, проходимый точкой B за время dt в положение B', определяет величину скашивания, которую можно найти как

BB' = dy- dy dt

Угловая деформация характеризуется тангенсом угла da. При

этом

, ч BB дих    .

tg(da) =-=    —- dt * da

AB    д y

(имея в виду, что AB = dy).

Вследствие малости угла da можно считать, что tg(da) * da.

Аналогично,

Полное скашивание первоначально прямого угла A определяется как сумма


(4.14)

Здесь следует обратить внимание на одно весьма существенное обстоятельство: рассматриваемое перемещение ребер вызвано не только деформацией, но и вращением частицы. Действительно, если бы грань только деформировалась без вращения, то ребра повернулись бы на одинаковый угол навстречу друг другу. Наоборот, если бы происходило только вращение, то ребра поворачивались бы на одинаковый угол в направлении вращения. Следовательно, в общем случае движение элемента можно рассматривать как сумму деформационного и вращательного движений, и таким образом определить da и dp. Рассмотрим

деформацию прямого угла A, считая, что вращение происходит против часовой стрелки. Чисто деформационное движение будем характеризовать углами dy, а чисто вращательное - ds.

ykdy


Из рис. 4.7 следует, что

dp= dy+ ds либо da+ dp= 2 dy, откуда

Вычитая, получим Рис. 4.7

Таким образом, деформация характеризуется полусуммой углов, а вращение - полуразностью. Имея в виду (4.14), можем записать:

г


1


dy =


dt


(4.17)


J


duv du.


dy    dx


+


Скорость угловой деформации, происходящей вокруг оси z

dY

dt


1

2


(4.18)

(4.19)

(4.20)


Yz


И по аналогии


Yy =

Yx =


V


J


dz    dx

du. + duy


du„ dux +


duv du. +


V dy    dx


dy    dz


J


о есть угловая скорость вращения жидкой

ds

dt


Выражение


частицы. Проекции угловых скоростей

о = —

оо = -

о = -

1

( duz

duy Л

2

V dy

dz J

1

(dux

duz

2

V dz

dx J

1

I

u

V

dux ]

2

v dx

dy J

Соотношения (4.21-4.23) играют исключительно важную роль в механике жидкости. Они устанавливают связь между угловой и поступательной скоростями жидкой частицы. Вопрос о знаках чисто

(4.21)

(4.22)

(4.23)


w n    w

условный. В гидромеханике поворот против часовой стрелки считается положительным, по часовой - отрицательным.

В векторной форме выражение для угловой скорости может быть записано как

о = exox + eyoy + ezoz    (4.24)

Заменяя ox, oy и oz их выражениями по (4.21-4.23) получаем:

г 1

о = — 2


дих    d Uz

dz дх


дх


ду


j


v


Г диг    диу

v dy dz


(4.25)

Сопоставляя выражение в квадратных скобках с формулой (1.8) видим их полную идентичность, поэтому можем записать:

r 1    ,    r

о = 2 rot и


2    (4.26)

либо

rot и = 2 о    (4.27)

Формула (4.27) раскрывает гидромеханический смысл вихря (ротора) векторного поля. Если и характеризует поле мгновенных скоростей, то векторное поле rot и представляет собой поле удвоенных угловых скоростей частиц жидкости этого поля.

4.8.2. Линейные деформации.

Очевидно, что линейные

с

о


Г


dy


dx


z-


Ю"


о


A


D


х


дх-ах

дх


Рис. 4.8


(A)


деформации частицы (рис. 4.8) могут возникнуть в результате различия в скоростях, совпадающих с направлением ребер. Как и ранее, компоненты скорости в точке A - их, иу, и

Вдоль оси х :

Точка A: и


'(A)


Точка D: их = их +


М D)


Разность скоростей,


вызывающая удлинение ребра AD: DD" за время dt

дих

дх


dx. Удлинение частицы


дих дх

dt = ds

дх    х


^х dt


(4.28)

(4.29)


Относительное удлинение

DD" AD


DD" =


Скорость относительного удлинения

dsx = дих dt    дх

(4.30)


Аналогично для других осей

duz

du


y


y    dyz    dz

Если процесс происходит одновременно вдоль всех осей, то это приводит к объемному расширению либо сжатию частицы. Таким образом, объемная деформация сводится к изменению первоначального объема параллелепипеда dV = dx dy dz на величину SV = SVx + SVy + SVz за счет растяжения либо сжатия ребер. При этом SVx = DD" dydz, и с учетом (4.28)

duy


du


duz

dz


x dVdt. Аналогично SVy = dVdt и SVz =

y    dy    z


SVx = Л x    dx

Таким образом


dVdt.


duv duv du..


S V =


dVdt


+


+


dx    dy    dz

Скоростью относительной объемной деформации назовем отношение изменения объема к его первоначальному объему и скорости деформации, т.е.


V


dux    du.


SV


duz    r

z = div u


+


+


drV dt dx dy dz Если div U = 0, то это означает, что SV = 0, т.е. деформация жидкой частицы происходит без изменения ее объема. В этом и заключается гидромеханический смысл равенства нулю дивергенции.

Полученную выше связь между поступательной и вращательной скоростями жидкой частицы можно получить и более коротким путем, представляющим определенный интерес. Разные подходы к одному и тому же вопросу способствуют углубленному пониманию. Поэтому рассмотрим этот путь.

Пусть жидкая частица вращается вокруг оси z с угловой скоростью (Oz. Запишем выражение для ротора в проекциях на оси координат (см. формулу 1.8). Имеем:

rotxU =



duz

y

u

d

dy

dz

dux

duz

dz

dx

y

u

d

X

u

d

dx

dy


rotyu =


rotzu =


rotzU = 0

z    дх

Аналогично для двух дру гих компонент

rotXu = 2ох; rotyU = 2оу

Либо в векторной форме

r 1 r о = 2 rot и

что полностью совпадает с (4.26).

Движение, при котором rot U Ф 0 называют вихревым, при rot U = 0 - безвихревым либо потенциальным. Из чего следует, что если течение вихревое, то движение жидких частиц происходит с вращением.


Рассмотрим точку M на жидкой частице (рис. 4.10).

Линейная скорость этой частицы и = 0z х r. Запишем выражения для проекций скоростей на оси координат:

их = °zy; и, = OzX;

uz = 0

ди,

Откуда находим -д-

Таким образом r    ди


5. ВИХРЕВОЕ ДВИЖЕНИЕ ЖИДКОСТИ

Вихревое движение широко распространено как в природе, так и в разного рода технических устройствах. Поэтому изучение его закономерностей представляет несомненный практический интерес. Вращательное движение жидких частиц характеризуется вихрем скорости


r 1 r о = 2 rot и


Это означает, что в каждой точке пространства вращение жидких частиц может быть охарактеризовано этим вектором. Его модуль

о = ^ о + о2у + о\    (5.2)

Движение, при котором величина вихря скорости не равна нулю, т.е. rot U ф 0, называют вихревым. При условии rot U = 0 движение безвихревое либо потенциальное.



дих

ду


о


- Оо


z


дих    2 О

— = 2 oz ду z


(5.1)


5.1. Кинематика вихревого движения.

Кинематические понятия для вихревого движения можно получить по аналогии с общими понятиями кинематики. В основу кинематики вихревого движения положено представление о вихревой линии, которое аналогично понятию линии тока. Вихревой называется линия, в каждой точке которой в данный момент времени вектор вихря скорости совпадает с касательной (рис. 5.1). Другими словами, вихревая линия - это мгновенная ось вращения частиц жидкости, которые в данный момент времени расположены на ней.

По аналогии с дифференциальным уравнением линии тока можно записать

dx _ dy _ dz

(5.3)

СОх    (Dy COz

Вихревая трубка - аналог трубки (поверхности) тока. Это поверхность, образованная вихревыми линиями, проведенными через все точки бесконечно малого замкнутого контура. Вихревая нить - аналог струйки - это жидкость, заключенная в вихревой трубке. Если вихревая трубка имеет конечные размеры, то частицы, заполняющие ее и находящиеся во вращательном движении, образуют вихревой шнур.

5.2. Интенсивность вихря.


Понятие интенсивности вихря достаточно абстрактно и вводится чисто математически. Напомним, что потоком векторного поля называют интеграл вида


и hdA    (5.4)

А

Поскольку вихрь скорости (ротор) есть вектор, то вместо и можно подставить rot й, что и приводит нас к понятию интенсивности вихря, т.е. интенсивность вихря - это поток вектора вихря

/'= ffrot и hdA    (5.5)

Рис. 5.1    JJ

А

Можно использовать и другую форму записи: rot и • h = rot п й]

i = jJrotn?f dA    (5.6)

A

Имея в виду, что rot й = 2 'со, можем записать

A    A

Воспользуемся формулой Гаусса-Остроградского и перейдем от интеграла по поверхности к интегралу по объему. Имеем:


®dV = 2

+


\


y    d®


л


i = 2


®ndA = 2


dV.


+


J


v    dx    dy    dz

A    V    V    7

Раскроем выражение, стоящее    под    знаком    интеграла,    имея в виду,

что проекции вектора вихря имеют    вид:


® = —

= -

1

Г duz

u

d

1

2

V dy

d

J

1

(dUx

duz']

2

d

N

1

d

X

1

i

u

d

V

dUx']

2

d

X

d

->2

du


Имеем

1 2


f п2 d U..


->2

d u,


->2

du


->2

du


->2

du


= 0.


+


+


dy dx    dx dz    dz dy dx dy

Следовательно, можно записать

JJ®ndA = 0    (5.8)

A

Заметим, что это выражение по структуре напоминает уравнение неразрывности.


dx dz    dy dz


j


¦    Применим    (5.8)    к вихревому шнуру

(рис. 5.2). На боковой поверхности соп = О, так как 'со направлен по касательной к поверхности. Поэтому можем записать

-    +    jjo)n2dA2 = 0;

А1 А2

А-\ А2 Если допустить, что в пределах сечения соп = const, то

= ^n2^2    (5-9)

Либо в общем случае

соА = const    (5.10)

т.е. это своеобразное «уравнение неразрывности». Полученный результат носит название теоремы Гельмгольца о Рис. 5.2    вихрях,    которую    можно    сформулировать

следующим образом: интенсивность вихревого шнура на всей его протяженности остается постоянной. Из выражения (5.10) следует и другой весьма важный вывод, сделанный Г.Гельмгольцем в 1855 г. в работе «Об интегралах уравнений, соответствующих вихревым движениям».Так как произведение соА остается неизменным, то уменьшение площади сечения шнура должно приводить к увеличению угловой скорости вращения частиц. При А = 0 со - оо, что физически невозможно. Следовательно, вихрь не может зарождаться либо оканчиваться в толще жидкости. Окончательно развившись, он должен замкнуться либо на твердую поверхность, либо сам на себя, т.е. образовать вихревое кольцо. Подробное описание этого явления можно найти в книге: Фабрикант Н.Я. Аэродинамика. - М.: Наука, 1964. - 814 с.

Понятие об интенсивности является весьма важным, но, к сожалению, непосредственное определение этой величины экспериментальным путем связано с непреодолимыми трудностями. Кроме того, если пытаться распространить это понятие на вихри конечных размеров, то по аналогии со средней скоростью пришлось бы вводить понятие о средней угловой скорости, что связано с определенными трудностями чисто математического характера. Поэтому гидромеханика избрала другой путь, заменив это понятие другим, более удобным для целей практики. К рассмотрению этого понятия, называемого циркуляцией скорости, мы и приступим.

5.3. Циркуляция скорости.

Для введения понятия о циркуляции скорости в настоящем пособии используется методика Н.Я.Фабриканта, приведенная в упомянутой выше книге. Несомненным преимуществом ее является то, что в отличие от других она позволяет ввести понятие циркуляции не чисто математически, а исходя из достаточно простых и ясных физических предпосылок.

Рассмотрим крыловой профиль, находящийся в потоке газа (воздуха). Как известно, на профиль в этом случае будет действовать подъемная сила (см. рис. 5.3). Физически наличие этой силы можно объяснить лишь тем, _ что давление под профилем (р^


Рис. 5.3    больше, а давление над

профилем (р2) меньше, чем давление на каком-то удалении от него, которое мы обозначим рх. Это позволяет утверждать, что под крыловым профилем скорость и| < иж, а над ним и2 > иж. В данном случае иж - скорость невозмущенного потока.

Вычтем теперь из скоростей и\и и2 скорость иж, т.е. и\ - иж и и2- и о,,. Это действие приводит нас к понятию потока возмущения, т.е. движения, которое возникает в среде из-за того, что в нее внесено инородное тело, т.е., по существу, это реакция потока, обусловленная в рассматриваемом случае тем, что в ней появился крыловой профиль. Установим теперь направление потоков возмущения. Под профилем и\ < иж, и он направлен против скорости и0о, над профилем - наоборот. В результате появляется циркуляционный поток, направленный по часовой стрелке, как это показано на рис. 5.3. Теперь необходимо охарактеризовать этот поток количественно. Именно с этой целью вводится понятие циркуляции скорости по замкнутому контуру.

Рассмотрим замкнутый контур С, показанный на рис. 5.4. Пусть в произвольной точке М скорость равна и. Составим скалярное произведение и • dl, где dl - направленный элемент дуги.

Циркуляцией скорости называют контурный интеграл вида

Г = d u • dl



(5.11)


= ju

Обратим внимание на структуру этого соотношения. Оно построено аналогично выражению для работы, поэтому иногда говорят, что циркуляция - это своеобразная «работа» вектора скорости.

Имея в виду, что u(ux uyuz) и

dl (dx, dy, dz), по правилу скалярного произведения получим

Г = uxdx + uydy + uzdz) (5.12)


Для плоского течения:

Г = <f( uxdx + uydy)

(5.13)


В конце предыдущего раздела утверждалось, что понятие циркуляции является более удобным, чем интенсивность вихря. Действительно, из (5.13) следует, что для определения циркуляции достаточно знать проекции скорости, нахождение которых не связано с существенными трудностями. Однако остается пока открытым вопрос о том, существует ли связь между циркуляцией и интенсивностью вихря. Ответ на него дает теорема Стокса.

5.4. Теорема Стокса.

В движущейся жидкости рассматриваем вихревое поле и выделяем в нем малый замкнутый контур со сторонами dx и dy (рис.

5.5). Пусть в начале координат скорости будут ux и uy. Запишем

выражение для элементарной циркуляции по этому контуру, имея в виду, что поток двумерный: d Г = uxdx + uydy.


Рассмотрим контур OABC. Если вдоль OA скорость ux, то вдоль CB ее приращение составит

dux

~ ay, и аналогично dy

Риг. Я Я

duv .

вдоль AB - ~dx dx . Это следует из выражения для полного

дифференциала скорости, например, dux = dux dx + dux dy.

dx dy

Запишем теперь выражение для элементарной циркуляции вдоль контура OABCO. Имеем:

г


du


y


dr = u„dx +


dx


dy -


uy +Л

y dx


ux +


V


du^ dy] dx - uydy dy yJ    y '


Раскрывая скобки и выполнив сокращения, получаем

dx dy = 2®zdA

dr =


dx    dy


Из чего следует, что циркуляция по бесконечно малому замкнутому контуру равна интенсивности вихря, пронизывающего этот контур.

Этот вывод легко обобщить и на случай произвольной кривой конечных размеров (см., например, Аржаников Н.С. и Мальцев В.Н. Аэродинамика. - М.: Оборонгиз, 1956 - 483 с.; упомянутую выше книгу Н.Я.Фабриканта).

Таким образом, можем записать:

Г = 2 jJ^dA = i

(5.14)


A

Это и есть формула Стокса, показывающая, что циркуляция по произвольному контуру равна сумме интенсивностей (напряжений) вихрей, пронизывающих поверхность, натянутую на контур.

6. ПОТЕНЦИАЛЬНОЕ ДВИЖЕНИЕ ЖИДКОСТИ

Как уже отмечалось, условием потенциаль ности движения является равенство нулю вихря скорости, т.е. rot и = 0. Физически это означает, что движение жидкости происходит без вращения частиц. Как будет показано, потенциальное движение играет исключительно важную роль в механике жидкости.

6.1. Потенциал скорости.

Сущность теоремы Стокса, по существу, сводится к утверждению о равенстве числовых значений интенсивности вихря и циркуляции, т.е. i = Г, либо

я-

rot U • n dA = Г


i =


A

С другой стороны, для потенциального потока по его определению rot U = 0, т.е. в потенциальном поле циркуляция по замкнутому контуру равна нулю.

Запишем выражения для проекций угловых скоростей.

Из сказанного выше следует, что для безвихревого (потенци


0. Следовательно, в этом случае


ального) движения юх =

duz    du


ю,


duz dUy дих


dy dzdz dxdx dy Эти соотношения позволяют существенным образом упростить вычисления компонент скорости их, иу и uz.

Рассмотрим выражение

uxdx + и ydy + uzdz    (а)

Оно построено аналогично известному из механики твердого тела выражению для элементарной работы. Зададимся вопросом, в каком случае (а) является полным дифференциалом. Напомним, что если выражение для работы является полным дифференциалом, то силы называются консервативными или имеющими потенциал. Ответ на


^ = —

ю

ю, = —

1

r duz

и

d

1

2

V dy

dz j

1

(dux

duz л

2

V dz

dx J

1

i

и

d

V

dux '

2

v dx

dy J

dUx


(6.1)


у.


поставленный вопрос был дан Алесисом Клодом Клеро (с жизнью и деятельностью этого удивительного ученого можно познакомиться по превосходной книге: Идельсон Н.И. Этюды по истории небесной механики. - М.: Наука, 1975. - 494 с.)

Клеро показал, что выражение типа (а) является полным дифференциалом, если обеспечивается равенство накрест взятых производных. Соотношения (6.1) как раз и удовлетворяют этому требованию, т.е. взятые накрест производные в (а) дают соотношения (6.1). Таким образом, при потенциальном движении выражение (а) является полным дифференциалом какой-то функции p, и

По предложению Гельмгольца функцию p называют потенциалом скорости.

Таким образом, всякому движению жидкости, происходящему без вращения частиц, соответствует свой потенциал скорости. Справедливо и обратное утверждение: если существует потенциал скорости, то движение происходит без вращения частиц.

Соотношения (6.4) можно получить и другим путем. Поскольку разные подходы к одному и тому же вопросу способствуют углубленному его пониманию, то получим эти же соотношения, используя другую методику.

Как уже отмечалось, условием потенциальности является rot u = 0. С другой стороны, как показано при рассмотрении операций второго порядка, операция ротора над градиентом какой-то скалярной функции тождественно равна нулю, т.е.

rot grad p = 0 Сопоставляя эти соотношения, можем записать

u = grad p    (6.5)

Это означает, что вектор скорости можно рассматривать как градиент какой-то скалярной функции p. Раскроем значения u и grad p. Имеем    rrrr

u = exux + eyuy + ejuz;

т.е. вновь приходим к соотношениям (6.4).

Пока что остается открытым вопрос о необходимости и целесообразности введения понятия о потенциале скорости. Чтобы разобраться в этом, следует иметь в виду, что к числу центральных задач гидромеханики относится определение сил, действующих на тела, обтекаемые потоками жидкости либо газа. Решение этих задач непосредственно связано с необходимостью расчета поля скоростей, т.е. определением проекций скоростей (их, иу, и7) в каждой его точке.

Из выражений (6.4) непосредственно следует, что все три компоненты скорости могут быть определены, если известна лишь одна величина - потенциал скорости. Таким образом, знание потенциала скорости существенно упрощает расчет поля. Однако немедленно возникает следующая проблема - как же найти потенциал скорости течения. Чтобы решить ее, необходимо прежде всего уяснить некоторые свойства, присущие потенциалу.

6.2. Уравнение Лапласа.

Операция дивергенции над градиентом скалярной функции приводит к оператору Лапласа. Если в качестве скалярной функции использовать потенциал скорости, то можно записать

..    .    ^2    д2ф    д2ф    д2ф

Выражения (6.7) и (6.8) носят название уравнения Лапласа. Таким образом, для нахождения потенциала скорости необходимо проинтегрировать уравнение Лапласа. Любая функция, удовлетворяющая этому уравнению, носит название гармонической. Следовательно, потенциал скорости является гармонической функцией. Как любое дифференциальное уравнение, уравнение Лапласа имеет бесчисленное множество решений, поэтому для того, чтобы однозначно определить потенциал скорости, необходимо задать граничные условия. Для задач, связанных с обтеканием тел, так называемых внешних задач гидромеханики, такими условиями являются un = 0 и

u = ux.

Первое условие характеризует безотрывность течения (равенство нулю нормальной компоненты скорости). Второе - показывает, что вдали от тела распределение скоростей известно.

Поверхности (либо линии для двумерных потоков), в каждой точке которых p = const, называются эквипотенциальными.

6.3. Циркуляция скорости в потенциальном поле.

Рассмотрим плоский (двумерный) поток. Выделим в нем произвольную кривую (рис. 6.1) и запишем выражение для циркуляции вдоль этой кривой

B    B    B

г = Juxdx + uydy =fdpdx + y-dy = Jdp = pa -pb (69)

A    A

A


т.е. циркуляция вдоль кривой не зависит от ее формы, а определяется лишь разностью потенциалов в ее конечных точках. Если кривая замкнута, то очевидно, что pB = pA и Г = 0, т.е. циркуляция по замкнутому контуру в потенциальном поле равна нулю.

B


6.4. Функция тока плоского течения.

В практических задачах гидромеханики двумерных потоков широчайшее применение находит понятие о функции тока. Рассмотрим двумерный поток и ограничимся несжимаемой жидкостью.

Как было показано, дифференциальное уравнение линии тока имеет вид

dx dy

u,

либо

uxdy - uydx = 0    (6.10)

Запишем уравнение неразрывности для этого случая, которое будет иметь вид

du d uy

+-L = 0

dx    dy

Аналогично тому, как это делалось при рассмотрении потенциала скорости, поставим вопрос об условиях необходимых и достаточных для того, чтобы выражение (6.10) являлось полным дифференциалом какой-то скалярной функции. Применим к (6.10) условия Клеро (равенство взятых накрест производных). Имеем:

du


y


= 0.


и


d uv du.


+


dx    dy    dx    dy

Но это есть не что иное, как уравнение неразрывности (6.11) для плоского потока, которое удовлетворяется всегда, если только движение существует. Следовательно, можно записать:

dy = uxdy - uydx    (6.12)

где у носит название функции тока. С другой стороны, поскольку, как показано выше, dy является полным дифференциалом, то можно записать:

(6.13)


ux =


uy =


Из чего следует, что если функция тока течения известна, то можно определить компоненты скорости в любой точке пространства. Сопоставляя (6.10) и (6.12) приходим к выводу, что если частица движется вдоль линии тока, то функция тока остается постоянной (при у= const, dy = 0 и (6.12) превращается в (6.10)). Проверим теперь, является ли функция тока гармонической функцией, т.е. удовлетворяет ли она уравнению Лапласа.

Для плоского потенциального течения юг = 0, но

г


duy du.


dy

dy


= 0, откуда


Из (6.14) ux


ю = -


V dx dy dx


dy


dx


dy


и uy =


следовательно


rdyу


d

dx


dy dx y


d


откуда

дх ду

Таким образом, функция тока, как и потенциал скорости, является гармонической функцией. И еще одно важное обстоятельство. Если потенциал скорости существует только в потенциальном потоке, то функция тока этим условием не ограничена. Это объясняется тем, что уравнение неразрывности, которое используется для получения этого понятия, справедливо как для вихревого, так и для безвихревого движений.

6.5. Гидромеханический смысл функции тока.

Рис. 6.2


Установим гидромеханический смысл функции тока, для чего проведем две достаточно близко расположенные линии тока (рис. 6.2). Вычислим объемный расход жидкости, протекающий между ними, для чего разложим вектор скорости частицы и на две составляющие их и иу, что позволит

представить расход как сумму dQ = dQx + dQy, при этом

dQx = uxdy и dQv = -uvdx

(рис. 6.2).

dQ = uxdy

Uydx

в


в


Q = [(uxdy - uydx) = (dty = Vb-Va    (6.15)

A    A

т.е. разность значений функций тока на двух смежных линиях тока равна объемному расходу между ними.

6.6. Связь потенциала скорости и функции тока.

Связь между этими параметрами может быть легко установлена, если записать полученные выше выражения для проекций скоростей

д(р _ fj    д(р

иу =

их =


ду

дх


откуда


д(р _ дц/ _ д(р _ дц/

дх дуду дх

(6.16)

Эти соотношения играют чрезвычайно важную роль в механике жидкости и носят название соотношений Коши-Римана. Более подробно они будут

рассмотрены ниже. Пока же ограничимся тем, что перемножим их. Это дает

дф_д<р_ = _дщ_ду_

дх ду ду дх

Из математики известно, что выражения типа (6.17) свидетельствуют о взаимной ортогональности кривых. Следовательно, линии тока и эквипотенциальные линии образуют сетку взаимно ортогональных кривых, которая носит название гидродинамической сетки движения. Примерный ее вид показан на рис. 6.3.

Рис. 6.3

6.7. Методы расчета потенциальных потоков.

Как уже отмечалось, для нахождения потенциала скорости необходимо проинтегрировать уравнение Лапласа при заданных граничных условиях. Задача эта достаточно сложна. Поэтому в теории потенциальных течений особый интерес представляют случаи, которые дают точные значения функций тока и потенциала скорости без интегрирования уравнения Лапласа. Общая идея такого подхода сводится к следующему: задаются какой-то функцией, которая заведомо удовлетворяет уравнению Лапласа и выясняют, что представляет собой гидродинамическая сетка движения. Эту методику рассмотрим на ряде простейших примеров.

Пример 6.1. Пусть выражение для потенциала скорости имеет вид = ах + by, где а и Ь - действительные числа.

Найдем компоненты скорости. Имеем


= а и и у =

Вторые производные равны нулю, т.е. уравнение Лапласа удовлетворяется. Так как ux = а и uy = b, то из этого следует, что

поток движется с постоянной скоростью

=4ат+ь

2 2

их + иу


и =


Выясним, что представляют собой линии тока. Дифференциальное уравнение линий тока

dy = uxdy - uydx = ady - bdx.

И после интегрирования

у = ay - bx    (6.18)


Приравнивая (6.18) какой-то постоянной, получаем семейство линий тока - параллельных прямых, наклоненных к оси под углом b

tg а = — (см. рис. 6.4). Действи-a

тельно, для линии тока можем записать:

dx = dy; dy =    =    Ь

их~ uydx ~ их~ a'

Пример 6.2. Потенциал скорости задан выражнием

(р = a (x2 - y2)

где a - действительное число. Необходимо найти линии тока этого течения.

Прежде всего проверим, удовлетворяет ли руравнению Лап

d2p


d2p dy2


dp

dy


ласа. Имеем    = 2ax


= -2ay;


= 2a:


= -2a


2


dx


dx


d2p


d2p


+


dx    dy

т.е. уравнение Лапласа удовлетворяется. Выясним, какое же движение описывается этой функцией, для чего установим вид функции тока.

dy= uxdy - uydx = 2ax dy + 2aydx =

= 2a(xdy + ydx) = 2a d (xy)


= 2a - 2a = 0,


Следовательно, ц/ = \2ad(xy) = 2 ах у (произвольная постоянная

в данном случае нас не интересует).

Для нахождения линии тока приравняем ц/ какой-то постоянной величине 2ах у = const либо ху = const. Следовательно, линии тока - гиперболы, для которых оси х и у - асимптоты. На рис. 6.5 показаны линии тока для верхней половины. Если считать, что оси ко-1нат являются твердыми стен-, то получим картину обтекания потоком прямого угла.

Существует ряд простейших течений, для которых потенциалы скорости могут быть получены аналитическим путем. Эти течения играют заметную роль в гидромеханике, и поэтому их рассмотрение представляет несомненный интерес.

Пример 6.3. Источник (сток) на плоскости. Ограничимся плоской задачей. Интересующиеся объемной (трехмерной) задачей могут найти ее в книге: Талиев В.Н. Аэродинамика вентиляции. - М.:

Изд. по строительству и архитектуре, 1954. - 287 с.

Под источником (стоком) на плоскости понимают точку, из которой происходит истечение (либо втекание) жидкости. Пусть точка О на рис. 6.6 представляет плоский источник, из которого, как из цен-

тра, проведем несколько концентрических окружностей. Запишем уравнение неразрывности для цилинд-

Рис. 6.6

В декартовой системе координат


(6.19)

В рассматриваемом случае удобней использовать цилиндрическую систему координат. Увязка систем может быть получена, исходя из рис. 6.7. Для цилиндрической системы

1 dp dp u° = 1 д; ur = д (621)

Рис. 6.7


Вывод этих соотношений можно найти в книге: Аржаников Н.С., Мальцев В.Н. Аэродинамика. -М.: Оборонгиз, 1956. - 483 с. Из (6.21) следует, что ur не зависит от полярного угла, поэтому можно за-dp

Приравнивая это


писать ur = —j—. r dr

выражение (6.19), получим


Q = dp

Q dr


dF, откуда dp= 2n r И после интегрирования

2nr


Q_

2n


ln r


p=


(6.22)


Из (6.22) следует, что эквипотенциальные линии источника представляют собой окружности. Формулу (6.22) можно записать и в следующей форме

p = 2nln -\/x2 + y2    (6.23)

Для нахождения функции тока удобней использовать декарто ву систему координат. При этом (6.19) принимает вид:

Q 1

(6.24)

ur =


2nJ x 2 + y2

С другой стороны, из рис. 6.7 следует:

ux = ur cosO = ur x/r

Таким образом

ux =


Аналогично


Дифференциальное уравнение функции тока

dy= uxdy - uydx

Подстановка значений ux и uy в (6.25) дает

Q xdy - ydx

dy =


x2 + y2

2п


f l'

v x


АуЛ


V x


x


т.е. xdy - ydx = x d Из знаменателя (6.26) выносим за скобки x2, при этом


Выполним некоторые преобразования. Дифференциал от ча-

xdy - ydx    .    .    2


стного имеет вид d


2 2 2 x + у = x

1 +

f у)

2

V x


Таким образом, (6.26) принимает вид

d

Q

2п


dy =

f L_2

1 +


V x J

f L'

v x


d


y)

V x J


Q

2п


Q_

2n



arctg


и y


f L'2

V x


1 +


Но с другой стороны — = tg 9, т.е. arctg(tg9) = 9, и

x

Q_

2n


9


(6.27)


y


В полярной системе координат (6.27) представляет собой семейство прямых, проходящих через начало координат. Для стока потенциал скорости и функция тока имеют те же выражения, но с противоположными знаками, т.е.


Q


Q


ln r и y =


9


(6.28)


р =


2п    г    ~    2п

Иногда Q называют мощностью (обильностью) источника.

6.8. Наложение потенциальных потоков.

Предположим, что имеются два потока с известными потенциалами скорости 1 и 2, удовлетворяющими уравнению Лапласа. Из теории линейных дифференциальных уравнений, к которым принадлежит и уравнение Лапласа, известно, что сумма частных решений этих уравнений также является их решением. Другими словами, это означает, что потенциал р, образованный как 1 + р2, также бу-


дет удовлетворять уравнению Лапласа, т.е. будет описывать какой-то новый поток, имеющий потенциал (р. Из этого следует, что можно получить новый поток путем сложения (наложения) уже известных. Следует обратить внимание на то, что собственно наложение потоков здесь не производится, а речь идет о сложении потенциалов скорости уже известных течений.

Скорость в каждой точке нового потока является суммой скоростей первоначальных потоков. Задача нахождения нового течения может быть решена как графически, так и аналитически.

Рассмотрим сначала графический метод. Общий подход сводится к следующему. Необходимо построить линии тока течений в одинаковом масштабе, что при достаточной густоте линий тока при пересечении дает фигуру, близкую к параллелограмму (рис. 6.8).

Рис. 6.8


Отрезки АВ и AD в каком-то масштабе представляют скорости течения, их результирующая определяется как диагональ параллелограмма (АС). Для построения такой сетки необходимо соблюсти следующее условие: расход между соседними линиями тока обоих течений должен быть одинаков.

В качестве примера рассмотрим картину течения, образующуюся при наложении плоско-го параллельного потока на сток (рис. 6.9). Как следует из рис. 6.9, частицы жидкости в новом течении будут двигаться по кривым, направленным к стоку.


Задача, как отмечалось выше, может быть решена и аналитически. В этом случае должны быть известны и у/ обоих течений.

Пример 6.4. Выполним сложение источника и стока с одинаковыми расходами, симметрично расположенными относительно начала координат на расстоянии а (см. рис. 6.10).

Потенциалы скорости:

Q

источника и -


Рис. 6.9

(Рс =


In ru ; стока


271


Q .

¦f- In rc.

Z.71


Выбираем произвольную точку М с координатами х и у. Потенциал скорости в этой точке <р = <ри + с, т.е.

Q    Q


(In ги - In гс) = — In

2п    с)    2    71    г„


Выполним некоторые преобразования этого соотношения. Из треугольников МИх и МСх получаем:


г„ =    +    (х2    +    а2)


2 + (х2 - а2)


=


Следовательно, потенциал скорости нового течения

Q


2 71

\

2

У

+

(*2

+ а2)

2

У

+

(*2

-I


In


(6.29)



Существенно больший интерес представляет функция тока.

Как было показано, щи = -Q- 6и и ц/с- -$-6с.

?.7t    ?.7t

Аналогично предыдущему

Q /    \

?=?u + ?c = jz\8u ~ 0о)

С другой стороны, из рис. 6.10 следует, что вс = 6+ ви, откуда 9и~ 9С = т.е. у/= -S-в. При этом условию у/= const (т.е. ли-

^ JL

нии тока) соответствует в - const. Таким образом, линии тока нового течения представляют собой окружности, проходящие через источник и сток.

Рассмотрим теперь картину, образующуюся при сближении источника и стока.

Пример 6.5. Забегая несколько вперед отметим, что получаемое при сближении источника и стока течение называется диполем.

В чем особенность рассматриваемой задачи? Если просто предположить, что расстояние a = 0, то ru = rc, и p и / тождественно равны нулю. Поэтому рассмотрим другой предельный случай. Пусть при 2а ^ 0 расход Q ^ да, но так, что произведение 2а • Q =

= const = M, где M носит название момента диполя. Таким образом,

Q = M    (6.30)

При этом потенциал скорости диполя

M Injy + (x2 + а2) - InJy2 + (x2 - а2)

MIn д|

pd = — lim v


Разберемся теперь в том, что представляет собой выражение, стоящее под знаком предела. Знаменатель можно рассматривать как приращение независимого переменного, а числитель - как соответствующее приращение функции. Действительно, рассмотрим функцию

In д/y2 + x2 . Придадим x значение x + a и x - a. Если теперь из значения функции, соответствующей x + a, вычесть ее значение при x-a, то получим числитель. Разность значений независимого переменного (x + a) - (x - a) = 2a есть знаменатель. Таким образом,

мы должны вычислить предел отношения приращения функции к приращению независимого переменного при стремлении последнего к нулю. Как известно, в математике такой предел называют производной функции, т.е.

M д .    Г~2    2

p    ln^x + y

Дифференцирование легко выполняется методом подстановок. Пусть u = д/x2 + y2 ; u = x2 + у2. Тогда z = In u; zf =u';

u

u ' =г (u *) 2V uv '

x

Уд —


M


x


M


2nx2 + у2

Из чего следует, что линии тока и эквипотенциальные линии - окруж ности, касающиеся осей Ox и Оу в начале координат (рис. 6.11). Действительно, придавая функции тока постоянные значения, получаем:

x2 + у2 = Су

M


2nx2 + у2 (6.31)

Действуя аналогичным образом, можно показать, что

У


2Vx2 + y2Jx2 + y2 x2 + у

т.е. -d- lnJx2 + y2 = 2 x 2. dx ^    * x2 + y2


2jx2+f

2x    1    x


Таким образом:

Рд


Имеем: u' —


где С — -


(6.32)



x2 + у2 - Су +

С2 = С2 4    4    ;

сГ\ = С

4 J 4

2


2

x +


Су +

f С_

2

f С_

+

У-тт

V

V 2 J

а это и есть уравнения окружностей с разными центрами.

Рис. 6.11

6.9. Бесциркуляционное обтекание круглого цилиндра.

Продолжим рассмотрение метода наложения потоков. Полученное в примере 6.5 течение, называемое диполем, на первый взгляд носит достаточно абстрактный характер. Однако, как будет показано ниже, такая точка зрения не совсем справедлива. Используя понятие диполя, можно получить весьма интересные и полезные для практических приложений результаты. Для подтверждения этого проанализируем течение, возникающее при наложении прямолинейного поступательного потока на диполь с центром, расположенным в начале координат. Прямолинейный поток движется вдоль оси Ох со скоростью, равной единице, т.е. их = uQ = const; иу = 0. Потенциал

скорости

d(p= uxdx + Uydy

и п = UqX с точностью до произвольной постоянной.

Функция тока dy/= uxdy - uydx = u0dy и y/= u0y. Если, как

принято в условии, и0 = 1, то п = х и ц/ - у. Примем для упроще-

..    х

ния выкладок момент диполя М = 2ж, тогда д -    2    и

X + у

у

Уд =--2-2 ¦ Складывая потенциалы и функции тока, получаем

X + у

Найдем линии тока, для чего приравняем функцию тока постоянной: у/ = у--2~^——2 = С, откуда

X + у

у (х2 + у2) - l] = С(х2 + у2)    (6.33)

Из чего следует, что линии тока течения представляют семейство кривых третьего порядка. Найдем нулевую линию тока, т.е. линию, для которой С = 0. Это дает два уравнения:

у = 0 их + у =1,


т.е. линия тока представляет собой ось х-ов и окружность единичного радиуса с центром в начале координат (см. рис.

6.12). Это позволяет рассматривать окружность как твердую границу и течение вне ее, что приводит к задаче обтекания бесконечно длинного цилиндра. Покажем, что на доста-Рис 6 12    точно    большом    удалении    от    ци-

линдра скорость направлена вдоль оси x и равна u= 1. Найдем проекции скоростей ux и uy.


\


2 2 0 2 x + y - 2x

(x2 + у2 )2


dp d


= 1 +


Имеем: uv =


x +


2 2 x + y J


dx    dx


У


x2 - y2;

x2 + y2


Откуда ux = 1


\


2xy

¦


dp d


аналогично uy =


x +


2 2 x + y


dy    dy


Для дальнейшего удобно перейти к полярным координатам, имея в виду, что x = r cose и y = r sin e. Подстановка этих значений в выражения для ux и uy дает:

, cos в- sin2в ux = 1--12--(6.34)


2 sin в cos#


(6.35)

= 0,


uy =


2


Перейдем к пределу. При r ^ да получаем ux^ = 1 и uy(x)


т.е. то, что и требовалось доказать.

Точки B и A, показанные на рис. 6.12, являются так называемыми особыми либо критическими точками, т.к. скорость в них обращается в нуль. Покажем, что это действительно так, для чего запишем выражение для потенциала скорости в полярных координатах:


r cose


= r cose+ 1Л


p x +--5-tz

x2 + y2


г


cose


(6.36)


r +

У    rJ


p


Найдем проекции скорости в произвольной точке на произвольной линии тока (рис. 6.13). Имеем:

dp d f 1^ „Л    1Л


cose = 1

у


cose


r +

у rJ


ur =


r2J


dr

1 dp r de


dr

f

1 +

у


r


\


1


cose = - 1 + 2de    У    r    j


sin e.


ue =


r2 J


¦    На поверхности цилинд-

ра г = 1 и иг = 0, т.е. обтекание безотрывно. Компонента ив = -2 sin в. В общем случае, когда и00*1,

ив = -2 их sin в (6.37) Знак «минус» указывает на то, что направление скорости на верхней половине цилиндра противоположно положительному направлению от-Рис 113    счета угла в. В точках в и Л

(sin в = 0) скорости равны нулю, т.е. действительно эти точки являются критическими.

6.10. Применение теории функций комплексного переменного к изучению плоских потоков идеальной жидкости.

Рассматриваемый ниже метод относится к числу наиболее эффективных способов анализа плоских потоков. Вернемся к полученным выше (см. 6.15) соотношениям Коши-Римана. Они показывают, что комплексная комбинация этих двух функций (т.е. и у/) от действительных переменных х и у, т.е. (р(х,у) + /у/(х,у), является аналитической функцией комплексного переменного z = х + iy Другими словами, эти условия показывают, что существует функция комплексной переменной W(z) либо просто И/, вещественная и мнимая

части которой ср и у/ соответственно, т.е.

l/l/(z) = + / у/ (либо W = (р+ /у/).

Функция W{z) называется аналитической в данной точке, если она дифференцируема как в самой точке , так и в некоторой ее окрестности. В гидромеханике функция W(z) называется комплексным потенциалом. Следует отметить, что теория аналитических функций является одной из наиболее разработанных ветвей классической математики. Обстоятельное изучение этого материала далеко выходит за рамки курса. Ограниченный объем данного пособия позволяет привести лишь весьма краткие сведения, необходимые для уяснения самой общей идеи метода. При необходимости подробное и обстоятельное изложение его можно найти в книге: Фильчаков П.Ф. Приближенные методы конформных отображений. - К.: Наукова думка, 1964. - 530 с.

Интересующиеся приложениями теории функций комплексного переменного для решения технических задач, в частности, задач гидромеханики, могут обратиться к книге: Лаврентьев М. А., Шабат Б. В. Методы теории функций комплексного переменного. - М.: Наука,

1987. - 688 с.

Как показывается в теории функций комплексного переменного, производная от комплексного потенциала W = ф + iy по комплексному же переменному z = x + iy имеет вид

дф . дш

~^Т + in = Ux - iuv дх дх x y


dW    "    -•    (6.38)

dz

Это выражение называется комплексной скоростью. Модуль этой величины дает саму скорость, т.е.

dW =^j U2 + Uj; = U    (6.39)

dz

Рассмотрим некоторые примеры.

Пример 6.6. Пусть течение задано комплексным потенциалом W = az2, где a - действительное число. Имея в виду, что W = ф + iy и z = х + iy, можно записать:

Ф + iy = a(х + iy)2 = ax2 + 2aixy- ay2 = a(x2 - y2) + i2axy

Разделяя действительную и мнимую части, получаем:

Ф = a(x2 - y2) и ш = 2axy.

Этот поток рассмотрен выше в примере 6.2. Обратим лишь внимание на то, что с помощью комплексного потенциала результат достигается более коротким путем.

Найдем комплексную скорость. Имеем:

dW

= Ux - iUy;

dz


uv = ^ = 2ax ;


dW


dz dW


x2 + y2 = 2ar,


U=


dz


дx    yдy

= 2ax + i 2ay = 2a (x + iy) = 2az; Uy = 2a^


дФ

u„ = —— = -2ay;


=    ,    lUx +


т.е. частицы движутся по гиперболическим линиям тока со скоростью u = 2ar.

6.11. Конформные отображения.

Геометрические преобразования, при которых величины углов между любыми двумя линиями, содержащимися в преобразуемой фигуре, не изменяются, называются конформными преобразованиями или отображениями. Широкое применение конформные отображения находят в гидромеханике. Обсудим лишь общую идею метода.

Рассмотрим две координатные сетки на плоскостях комплексных переменных zx + iy и Z — %+ in (рис. 6.14).

В плоскости z имеется какая-то фигура (A), которую необходимо отобразить на плоскость Z Эта операция может быть выполнена при одном непременном условии: должно быть известно соотношение, устанавливающее связь Z и z, т.е. Z — f (z). Эта зависимость

носит название отображающей функции. Предположим, что она нам известна. Тогда, задавшись какой-то произвольной точкой на контуре A, например 1, можно вычислить z1, и подставив это значение в отображающую функцию, найти значение Z и соответствующую точку на плоскости Z (1'). Повторив эти операции для точек 2, 3 и т.д., найдем 2', 3', ... . В результате этих действий получим контур B на плоскости Z, т.е. контур A отобразился в контур B. Такое преобразование получило название конформного. В теории функций комплексного переменного доказывается, что модуль производной р' (z) характеризует

изменение линейных размеров области при преобразовании, а аргумент ее определяет угол поворота радиуса-вектора. При этом преобразование, осуществляемое аналитической функцией, сохраняет эти углы во всех точках, где производная отображающей функции отличается от нуля. Теперь вопрос может быть сформулирован таким образом: какие же практические преимущества можно получить, используя метод конформных отображений?

Рис. 6.14

Остановимся лишь на одном, но крайне важном случае. Как известно, одной из главных задач расчета крыля является определение его подъемной силы. Для ее нахождения необходимо знать скорости частиц в каждой точке потока, обтекающего крыло. Крыловой профиль - достаточно сложная фигура, и рассчитать скорости теоретическим путем не представляется возможным. Но, как было показано выше, расчет легко выполняется для цилиндра. Поэтому задача была бы решена, если бы удалось заменить обтекание крылового профиля обтеканием цилиндра. Это можно сделать с помощью конформного отображения.

Рассмотрим рис. 6.15. Конформно отобразив фигуру, заштрихованную на рис. 6.15а (внешность профиля) на заштрихованную фигуру рис. 6.156 (внешность окружности) мы сводим задачу обтекания профиля к задаче обтекания цилиндра. Рассчитав скорость в любой точке цилиндра, обратным переходом можно найти скорость в соответствующей ей точке профиля.

Рис. 6.15

Нахождение вида отображающей функции, позволяющей осуществить требуемое конкретными условиями рассматриваемой задачи конформное отображение, является отдельным специальным вопросом. Решение рассмотренной выше задачи было найдено Н.Е.Жуковским. Отображающая функция в этом случае имеет вид

С 2\

г

(6.40)


Z +

V

и носит название функции Жуковского.

7. ГИДРОДИНАМИКА ИДЕАЛЬНОЙ ЖИДКОСТИ

В механике жидкости понятию «гидродинамика» придается весьма широкий смысл. В настоящем пособии этот термин будет использоваться в его классическом значении, как раздел курса, который, в отличие от кинематики, рассматриваюшей движение жидкости без учета причин, обусловивших его, изучает как само движение, так и причины, приводящие к его возникновению. Движение жидкости вызывается действием сил, а если иметь в виду, что давление есть частное от деления силы на площадь, то можно считать, что причиной возникновения движения частиц с какими-то скоростями является разность (перепад) давлений. Таким образом, для расчета течений необходимо иметь уравнение, связывающее давление в точке со скоростью движения частицы.

7.1 Уравнения движения идеальной жидкости.

Уравнения движения идеальной жидкости можно получить из уравнений движения в напряжениях, положив в них все производные от т равными нулю и заменив нормальные напряжения давлениями, имея в виду, что pxx = pyy = pzz = - p. Таким образом, уравнения гидродинамики принимают вид

Система (7.1) называется системой дифференциальных уравнений Эйлера для гидродинамики, она связывает давления и скорости в движущейся жидкости. Следует помнить, что выражения в правой части уравнений системы являются полными либо субстанциональными производными. Наличие конвективных членов ускорения приводит к тому, что система является нелинейной, содержащей четыре неизвестных: три проекции скорости и давление. Проекции единичных массовых сил обычно известны из постановки задачи.

Три уравнения (7.1) плюс уравнение неразрывности образуют замкнутую систему.

7.2 Преобразование Громеки-Лэмба.

Рассмотрение теоремы Гельмгольца о движении жидкой частицы показывает, что жидкость, как любое материальное тело, может участвовать в поступательном и вращательном движениях.

Следует обратить внимание на то, что для совершения работы в современных технических устройствах может использоваться только энергия поступательного движения. Энергия же вращательного (вихревого) движения полностью теряется, рассеивается в окружающей среде, превращаясь в теплоту.

Система уравнений Эйлера (7.4) не учитывает факт существования этих двух движений, что в определенной степени обедняет ее. Поэтому целесообразно использовать преобразование, позволяющее учесть эту особенность движения жидких частиц, называемое преобразованием Громеки-Лэмба. Формально оно сводится к тому, что в выражение для ускорения вводятся члены, характеризующие вращение жидких частиц.

Рассмотрим лишь одну компоненту:

dUx


dUx


dUx


(7.3)


+ Ux


+ Uy


+ Uz


dt ' ~x dx    ~y dy    ~z dz

Прибавим и вычтем в конвективной части ускорения выражение


ах =


dUx


duy

dUz

dx


+ Uz

Uy


y dx

Скомпонуем члены с учетом знаков

r


\


dUy

dx

dUs


1 d


1    dU_

2    dx


dU


dUx

dx


/ 2 2 2\

(Ux + Uy + Uz) =


U


+


U


+ U-.


2 dx

dU


dx


J


d Uz

dx


dUx

' dz


U


+ U.


U


+ Uz


y


dy


dx


V


d Uv    dUz


= - U,


+ Uz


dx


dy


dz    dx


V


V


Выражения в скобках есть не что иное, как удвоенные компоненты вихря (oz и coy, т.е. можем записать

-2Uyyz + 2Uzyy = 2(Uzyy - Uywz)

Подставляя полученные значения в (7.3) имеем

(7.4)


ax =


+


U yO z


zy y


d Uv 1 dU2

+


dt    2 dx


и по аналогии

7.3 Уравнение движения в форме Громеки-Лэмба.

Если в (7.2) в правую часть подставить ускорение в виде (7.7)

либо (7.8), то это приводит к уравнению движения в форме Громеки-

Лэмба. Для установившегося движения имеем

r 1    u2    r r

F --grad p — grad + 2wx u    (7.9)

Выполним некоторые преобразования (7.9).

В разделе гидростатики было введено понятие о скалярной функции Ф, называемой силовой. Было показано, что

Xdx + Ydy + Zdz — dФ    (7.10)

Поскольку эта функция является полным дифференциалом, то можно записать

дФ . дФ . дФ . dФ —    dx dy + dz    (7.11)

дх dy J dz    v 1

Сопоставляя (7.10) и (7.11), получаем

дФ — X; дФ — y д — z

дх    дy    дz    v 1

¦ _

С другой стороны вектор F, проекциями которого являются X, Y, и Z

F — exX + eyY + ezZ    (7.13)

Из (7.12) и (7.13) следует, что

r r дФ г дФ г дФ дх + ey дy + Г дz С учетом (7.14) выражение (7.9) принимает вид

F — ex^r + ey    + ez^ — grad Ф    (7.14)


2 p

J


Следует иметь в виду, что эта форма записи справедлива лишь для несжимаемой жидкости, т.е. при условии p= const. И, наконец, уравнению движения (7.15) можно придать более удобную для анализа форму, умножив скалярно его левую и правую части на произвольный направленный отрезок

dl = ex dx + eydy + e dz    (7.16)

Опуская эту операцию, которую обучающийся при желании может выполнить самостоятельно, приведем лишь конечный результат

2

u

2

Л

dx

dy

dz

+ p - ф

= -2

ux

uy

uz

P

(Dx

Dy

Dz


d


(7.17)


v


7.4 Интегрирование уравнения движения для установившегося течения

Интегрирование уравнения движения (7.17) возможно лишь в случае, когда его правая часть равна нулю. Из теории определителей известно, что признаками равенства нулю являются: равенство нулю какой-либо строки или пропорциональность элементов одной строки элементам другой.

Исходя из физического смысла имеем четыре возможных случая:

(7.18)

(7.19)

(7.20)

(7.21)


x =

Dy =

Dz =

ux

= uy

= uz

Dx

Dy

Dz

dx

= dy

= dz

Dx

Dy

Dz

dx

= dy

= dz

ux

uz

\


uy


Для любого из них можем записать

2

u p


d


= 0


+ — - Ф P у


2


v


И после интегрирования


+ ^-Ф = с

u2


(7.22)


2 P

Если из массовых сил действует только сила тяжести, то, как показано в разделе гидростатики,

ф = - gz

и (7.22) принимает вид

p    U2

gz + P +    =    C    (7.23)

Еще раз обратим внимание на то, что вид уравнения (7.23) одинаков вне зависимости от того, какой из четырех случаев равенства нулю определителя рассматривается. Однако смысл интеграла и область его применения различны. Именно поэтому следует разобраться в этом вопросе подробней.

Первый случай, как известно, является признаком потенциальности движения. Интеграл (7.23) в этом случае называют интегралом Коши-Лагранжа. Он справедлив для любых точек жидкости, движущейся без вращения частиц, т.е. потенциально.

Второй случай является признаком коллинеарности вектора вихря и вектора скорости. Это весьма редкий случай так называемого винтового движения.

Третий случай характеризует движение жидкой частицы вдоль вихревой линии, а четвертый - движение вдоль линии тока. Интеграл

(7.23) при этом носит название интеграла Бернулли. Он справедлив как для потенциального, так и для вихревого движений. Именно этот случай и будет интересовать нас в дальнейшем.

7.5 Упрощенный вывод уравнения Бернулли.

В ряде пособий и учебников рассматривается упрощенный вывод уравнения Бернулли. Поэтому с целью расширения и углубления представления об этом основополагающем уравнении механики жидкости представляется целесообразным рассмотреть и этот подход. В основу его положено принимаемое без каких-либо доказательств положение о том, что рассматривается жидкая частица, движущаяся вдоль линии тока. После чего производится преобразование системы дифференциальных уравнений Эйлера (7.1) путем умножения каждой из его проекций соответственно на dx, dy и dz и почленного их сложения аналогично тому, как это делалось в гидростатике. Это преобразование уже рассматривалось в случае, когда из массовых сил действуют лишь силы тяжести (см. раздел «Гидростатика»). Оно

приводит к соотношению: -gdz - dpp. Поэтому рассмотрим лишь

правую часть. Имеем

dUx.dx + dUzdv + dU±dz dt + dt ay + dt d

dx    dy    dz

Считая, что df    = Ux;    dtf    = Uy; df    = Uz,    можем записать:

111

Uxdx + Uydy + Uzdz = 2 dU2x + 2 dU 2 + 2 dU\ =

Это выражение называют уравнением Бернулли в дифференциальной форме. При условии p= const (для несжимаемой жидкости) интегрирование его дает


Таким образом


либо


(7.24)


(7.25)

т.е. соотношение (7.23).

Очевидно, для обеспечения математической строгости следовало бы доказать, что вдоль линии тока проекции вектора скорости могут быть представлены не как частные, а как полные производные от соответствующих координат частицы. Но при этом вывод уравнения Бернулли утратил бы свою простоту.

7.6 Энергетический смысл уравнения Бернулли

Прежде чем приступить к анализу физического содержания полученного соотношения, следует вспомнить одно важное обстоятельство. При введении понятия о струйке было показано (см. раздел «Кинематика»), что одним из ее свойств является равномерное распределение скоростей в пределах любого ее поперечного сечения. Это означает, что соотношение (7.25) остается справедливым для любой линии тока, проходящей внутри струйки. Поэтому уравнение


(7.25) можно назвать уравнением Бернулли для струйки идеальной жидкости. Для двух произвольных поперечных сечений струйки можно записать

(7.26)

Выясним физический смысл величин, входящих в уравнение Бернулли. Любое правильное физическое соотношение размерност-но однородно, т.е. все его члены имеют одинаковую размерность, поэтому достаточно рассмотреть один из его членов. Наиболее удобно

обратиться к третьему - u/2. Эта величина выражается в м22. Умножим и разделим числитель и знаменатель на кг, что дает:

г / /2


с2 кг &кг    кг    кг

Из чего следует, что каждый член уравнения выражает энергию, отнесенную к единице массы, т.е. удельную энергию. Это позволяет придать уравнению Бернулли энергетический смысл. Первые два члена выражают удельную потенциальную энергию (положения - gz и Я


), а третий - удельную кинетическую энергию. Следо


кг - м ¦ м    Н м Дж

-» -о-->-->


вательно, полная удельная энергия в любом сечении струйки остается неизменной. Другими словами, уравнение Бернулли выражает закон сохранения энергии в ее простейшей форме - форме сохранения механической энергии.

7.7 Уравнение Бернулли в форме напоров.

В практических приложениях широко используется другая форма уравнения Бернулли - форма напоров. Разделив обе части уравнения (7.26) на ускорение свободного падения д, получаем


полного напора (трубке Пито). Принцип действия этого устройства легко

уясняется из рис. 7.1.

Сумма двух первых членов носит название гидростатического напора, а трех - полного либо гидродинамического напора. Таким образом, уравнению Бернулли придается геометрическое толкование, которое сводится к следующему. Сумма трех высот: геометрической

(z), пьезометрической (у^д)и скоростной (и/2д) есть величина по-


рд 2д    ' рд ' 2д    (7    27)

Каждый член (7.27) имеет линейную размерность и выражает напор, под которым в общем случае понимают высоту столба жидкости, уравновешивающую давление в данной точке. Таким образом, z - геометрический напор, характеризующий положение жидкой частицы над какой-то произвольной плоскостью, называемой плоскостью

Я


Zi + Pi + ^ = Z2 + P2 + oI


давления


отсчета;

пьезометрическии на


ке;

Рис. 7.1


пор - высота столба жидкости, уравновешивающая давление в данной точ


скоростной напор, пред



ставляющий собой высоту столба жидкости в так называемой трубке


стоянная вдоль струйки. Либо, что то же самое, полный либо гидродинамический напор при движении вдоль струйки остается неизменным. Сказанное иллюстрируется рис. 7.2, который иногда называют диаграммой уравнения Бернулли.

На рис. 7.2 N-N - напорная линия; O-O - плоскость (линия) отсчета; P-P - пьезометрическая линия, лежащая ниже напорной на величину скоростного напора в данном сечении.

8. ГИДРОДИНАМИКА ВЯЗКОЙ ЖИДКОСТИ

8.1 Модель вязкой жидкости

Приступая к рассмотрению движения вязкой жидкости, необходимо прежде всего уяснить терминологию , т. е. смысл, вкладываемый в понятие «вязкая жидкость». С математических позиций необходимо установить вид функциональной зависимости для напряжений, либо, другими словами, сформировать модель вязкой жидкости. В дальнейшем под вязкой мы будем понимать жидкость, удовлетворяющую трем гипотезам: линейности, однородности и изотропности.

8.1.1. Гипотеза линейности .

Применим закон Ньютона к жидкости, движущейся параллельно плоскости xOy (рис. 8.1), что дает

диГ Tzx =>‘-дГ

Воспользуемся результатом, полученным при рассмотрении теоремы Гельмгольца о движении жидкой частицы. Согласно теореме, скорость угловой деформации относительно оси y

d Uv dU. + '

Yy

Рис. 8.1


1


dz    dx j

Так как движение происходит в плоскости xOy, то Uz = 0 и

и, следовательно, касательное напряжение

Tzx = 2m Yy    (8.1)

Полученный результат иллюстрирует так называемый закон трения Стокса. Согласно этому закону, напряжения, возникающие в жидкости, в отличие от твердого тела, пропорциональны не величинам, а скоростям деформаций, и связаны с ними линейной зависимостью. При этом коэффициент пропорциональности остается неизменным и равным

Кроме того, согласно закону Стокса касательные напряжения, как показано выше, пропорциональны скоростям угловой деформа-

dUx dUy

y


w    w    I    X    ¦

ции, а нормальные - скорости линейной деформации, т.е.

dxdy

dUz dz

(8.2)


+


xy


yx


V


J


Таким образом, можем записать

f dU,


dx    dy


= T


и т.д.

Рассмотрим теперь нормальные напряжения, возникающие от сил вязкости. Согласно закону Стокса, их можно записать в виде так называемых девиаторов напряжения, имеющих вид:

dUx

Cxx = 2^

dx

dUy

dy


2 duz

°zz = 2V~dZT

Полные нормальные напряжения отличаются тем, что помимо записанных выше в любой, как в вязкой, так и в невязкой жидкости, действуют и статические давления. Другими словами

о dux

pxx = p + 2j dx

Pyy =- P + 2^-dyy    (84)

о duz pzz = p + 2j dz

Выполним следующую операцию: из утроенной величины pxx вычтем сумму (Pxx + Pyy + Pzz). Это дает:

d uy

3Pxx - (Pxx + Pyy + Pzz) = -3P + 6

dx

f о..    d,,

dux    duy    duz

x +-- + z


- dux    n    .. г

= 6j—- -2jdiv u dx


-3p + 2j

V dx dy dz J


откуда найдем

P    2 dux 2 div u I Pxx + Pyy + Pzz

Pxx = 2^дТГ - 3 jdiv u +-y-

В качестве давления в вязкой жидкости принимают среднее

Pxx + Pvv + Pzz

арифметическое, т.е. P =--3-. И, следовательно,

Pzz = -P + 2^-^zz - 3 jdiv u

Для несжимаемой жидкости div u = 0, и выражения упрощаются.

8.1.2.    Гипотеза однородности

Предполагается, что вид линейной зависимости между напряжениями и скоростями деформаций одинаков для всех точек пространства.

8.1.3.    Гипотеза изотропности

Вязкая жидкость предполагается изотропной, т.е. ее свойства в любом направлении одинаковы.

8.2 Уравнение движения вязкой жидкости. (уравнение Навье-Стокса)

18 марта 1822 года в докладе, представленном Французской академии наук, Клод Луи Навье писал о полученных им уравнениях: «Хотя уравнения основаны на гипотезе Ньютона о том, что касательные напряжения пропорциональны скорости деформации, никак нельзя сказать, что они не выражают ничего нового».

Уравнения движения вязкой жидкости можно получить из уравнений движения в напряжениях (2.16), выполнив некоторые преобразования. Рассмотрим лишь одну проекцию этих уравнений:

Г


dT


dUx

dt


= X +1

Р


dp


dT


yx


xx


zx


+


+


dx


dy


dz


J


Как было показано при рассмотрении модели вязкой жидкости, нормальные напряжения

о dUx 2    ..    r

Pxx = -p + 2Md- - 3Mdiv u

Для упрощения задачи будем считать жидкость несжимаемой (div U = 0), тогда

Г    о    ..    \


dp + 2m ^Ux


dp


d


о dUx

p + 2m— dx


xx


(8.6)


dx    dx


dx


dx


V


r


Касательное напряжение Tyx = м


+


V dx \


dy


J


->2

d Uy

dy dx


d2Ux

dy2


dT


d Uy dU


d

dy


yx


(8.7)


+


+ M'


= M


= M


dy

аналогично

dT


dx


dy


J


2

d Ux 2

dz


2

d Uz

dz dx


d

dy


dUv dU.

+ '


zx


(8.8)


+ м


м


м


V


dz


dz    dx


J


Суммируя (8.6), (8.7) и (8.8) и группируя члены, получаем:

( п2 d U.


2

d U,


2

d U


2

x + d Ux


dp

dx


d U.


+ м


+


+ м


+


+


dy"


dx    dy dx    dz dx


dx'


dz


V


J


V


J


Третий член можно записать в виде:


d Uv    d Uv    dU.


d

dx


d r div U


+


+


м


= M


dx


dy    dz


dx


V


но жидкость несжимаема, и div U = 0. Таким образом получаем:

222

d ux +d ux +d ux

dux    1 dp и

x = x---— +


v dx2    dy2    dz2 j


dt    p    dx p

Выражение в скобках есть ни что иное, как оператор Лапласа

222

d ux    d ux d ux

A +    A    +    A


(8.10)


v dx2    dy2    dz2 J


V2ru, а P = V. Окончательно получаем:

dux dux    dux dux ХУг 1 dp

x + u^^r^ + u»—^ + uz^r2L = X---+


dt    x dx    y dy    z dz    p dx

+ V


Аналогично можно расписать и две другие проекции. Полученная система уравнений движения вязкой жидкости и носит название системы уравнений Навье-Стокса.

В векторной форме можно записать

r r 1    r

a = F -pgrad p + vV2u    (8.11)

Как следует из (8.11), это уравнение отличается от уравнения движения идеальной жидкости дополнительным членом (vV2u), учитывающим действие сил вязкого трения.

Целью гидродинамического расчета является нахождение полей скоростей и давлений, т.е. в результате расчета должны быть найдены четыре величины: ux, uy, uz и p. Принципиально это оказывается возможным, так как три уравнения Навье-Стокса (в проекциях) плюс уравнение неразрывности образуют замкнутую систему. Плотность и вязкость, входящие в них, считаются известными, а проекции массовых сил (X, Y, Z) задаются условиями конкретной задачи.

С чисто математических позиций уравнения Навье-Стокса относится к классу нелинейных дифференциальных уравнений в частных производных второго порядка. Одно из наиболее неприятных из их свойств - нелинейность, обусловленная наличием конвективных членов ускорения. Следует отметить, что до настоящего времени вследствие практически непреодолимых математических трудностей не получено ни одного общего решения уравнений Навье-Стокса в их полном виде, т.е. при сохранении всех конвективных членов и всех членов, учитывающих вязкость. Известны лишь отдельные частные решения.

Одним из основных граничных условий при интегрировании является условие «прилипания», т.е. равенство нулю скорости жидкости на стенке.

9. ОДНОМЕРНЫЕ ТЕЧЕНИЯ НЕСЖИМАЕМОЙ ЖИДКОСТИ (основы гидравлики).

Одномерными называются течения, в которых основные параметры потока зависят лишь от одной координаты, направление которой совпадает с направлением вектора скорости. Использование одномерных течений позволяет достаточно просто решать многие важные прикладные задачи. Раздел механики жидкости, изучающий одномерные течения, называют гидравликой.

9.1 Расход потока и средняя скорость.

Для решения широкого круга прикладных инженерных задач плодотворной явилась введенная Эйлером так называемая струйная модель потока. Согласно этой модели поток представляется состоящим из бесконечного множества струек жидкости. При рассмотрении потока поперечные сечения в нем выбираются так, чтобы пересекающие их линии тока были нормальны к ним. В этом случае сечение потока называется «живым». Очевидно, что если линии тока параллельны, то живое сечение будет плоским.

Ранее, в разделе «Кинематика», было показано, что элементарный объемный расход несжимаемой жидкости может быть определен как

dQ = udA

(9.1)


где u -скорость в сечении струйки, dA - площадь ее поперечного сечения.

n    w    w

В соответствии со струйной моделью расход потока

A

Рассмотрим движение жидкости в трубе круглого поперечного сечения. В силу тормозящего действия сил вязкого трения распределение скоростей в поперечном сечении трубопровода (эпюра скорости) будет иметь вид, показанный на рис.

Рис. 9.1


v


9.1. Для удобства перейдем к цилиндрическим координатам (r, в),

где в - полярный угол.

В этой системе

dA = rdr de Подставляя (9.3) в (9.2) получаем

0 0

Имея в виду, что

2п,имеем

0

R

Q = 2п J u (r) r dr

(9.5)


0

Запись u (r) обозначает, что местные скорости в сечении трубы изменяются по радиусу. Другими словами, u (r) описывает закон изме

нения скорости, т.е. является математическим описанием эпюры. Следовательно, для того, чтобы вычислить расход по (9.5), необходимо знать уравнение эпюры скорости, которое, как правило, неизвестно. Поскольку расход является важнейшим параметром, знание которого требуется при проведении любых гидравлических расчетов, необходимо найти путь, позволяющий преодолеть возникшее затруднение.

Рассмотрим, как решается эта задача в механике жидкости. С чисто математических позиций интеграл в правой части выражает объем эпюры скорости. Представим теперь, что при неизменном расходе Q в силу каких-то причин жидкость потеряла вязкость. Это, очевидно, приведет к тому, что эпюра начнет перестраиваться и, так как исчезнут силы вязкого трения, то все частицы жидкости будут двигаться с какой-то одинаковой скоростью v (см. рис. 9.1), а так как по условию расход остается тем же, то объем новой эпюры равен объему старой. При этом условии u (r) = v = const, и из (9.5) получаем

R

R

0 0


R2

r dr = 2nv — = vA 2


(9.6)

Скорость v, введенная таким образом носит название средней либо среднерасходной скорости. Следовательно, формально средняя скорость может быть определена как фиктивная скорость, с которой должны были бы двигаться все частицы жидкости для того, чтобы расход был равен его истинному значению.

С физической точки зрения использование понятия средней скорости, одинаковой для всех частиц жидкости в сечении, позволяет свести задачу о движении жидкости в трубах и каналах к одномерной.

9.2 Слабодеформированные потоки и их свойства.

Рассмотрим движение в трубе с несколько иных позиций. Если считать его установившимся, то все производные по времени, входящие в уравнение движения, равны нулю. Если исходить из одномерной модели, то равны нулю и компоненты скорости Uy и Uz. При

dUx ~

этом из уравнения неразрывности следует, что ~d^ = 0. Применительно к этому случаю система дифференциальных уравнений На-вье-Стокса принимает вид:

С п2    ~2

d Ux +d Ux


V 1 dp

X---— + v


= 0


2


р dx


dy    dz


V


Y -- dp = 0    (9.7)

Р dy

Z - 1    =    0

р dz

Последние два уравнения (9.7) совпадают с уравнениями гидростатики, а это означает, что в плоскости поперечного сечения движущейся жидкости давления распределены по гидростатическому закону

z + ppg = const (либо gz + р = const)    (9.8)

Этот вывод приближенно справедлив для слабодеформированных потоков. Под слабодеформированными понимают потоки, у которых угол расхождения линий тока мал, а радиус кривизны - велик, т.е. понятие это носит скорее качественный, чем количественный характер.

9.3 Уравнение Бернулли для потока вязкой жидкости.

Как уже отмечалось, уравнение Навье-Стокса в подавляющем большинстве случаев не поддаются интегрированию. Вместе с тем, практическая деятельность, связанная с необходимостью использования законов движения жидких сред, настоятельно требовала разработки инженерных методов расчета.

Одним из путей решения этой задачи, оказавшимся наиболее плодотворным, явился путь обобщения уравнения Бернулли, т.е. распространения его на поток вязкой жидкости. В основу этого метода, как уже отмечалось, положена струйная модель - представление

о потоке как о бесконечно большой сумме струек, протекающих через сечение.

Исходим из того, что движение установившееся и в рассматриваемом сечении поток слабо деформирован. Определим энергию, проносимую секундной массой струйки через сечение (т.е. мощность струйки ). Эта величина может быть найдена как произведение пол-

2

ной удельной энергии струйки (gz + p + —) на ее массовый расход

p 2

(pu dA). В справедливости этого легко убедиться непосредственно. Действительно, удельная энергия -Дж/кг, массовый расход - кг/с, их произведение

Дж кг = Дж кг с ~ с

Таким образом

2

p


dN =


)


Секундная энергия (мощность) потока в соответствии со струйной моделью

2


= Вт


pudA


(9.9)


v


pudA

N =


(9.10)


P

A


либо

N = pjjfgz + — udA + Pp JJu3dA

(9.11)


A    P    A

Так как поток слабодеформированный, то gz + р = const и первый интеграл принимает вид

р( gz + -! J udA = p( gz + P Q

(9.12)


v Py A    v    p)

Г p]

gz + -

py


q+pp jj u

A


N =


dA


(9.13)


p


По физическому смыслу второй член в (9.13) представляет собой кинетическую энергию секундной массы.

Поскольку мы ограничимся одномерным представлением, то в (9.13) необходимо ввести среднюю скорость. Поступим следующим

образом: разделим обе части уравнения на массовый расход pQ, т.е. отнесем это соотношение, как и уравнение Бернулли для струйки,

Дж • с Дж

--> —— и,


N_

pQ


кг • с


кг


к единице массы (N ^Дж/с; pQ ^кг/с;

i= N

следовательно, E = —~ -удельная энергия.


pQ

Таким образом, имеем


E = gz + Р + —u

dA


(9.14)


У p 2Q JJ

A

Разделив и умножив третий член на квадрат средней скорости v2, с учетом того, что Q = vA , получим

E = gz + p + —-3— I u3dA    (9.15)

p 2 v A A

A

'a I u


dA = a; тогда


A


Обозначим выражение 3

v A


p    av2

E = gz + — + p

(9.16)


Величина a носит название коэффициента кинетической энергии, корректива скорости либо коэффициента Кориолиса. Физический смысл этой величины будет раскрыт позже.

Разделив обе части (9.16) на ускорение свободного падения g, выразим это соотношение в единицах длины, т.е. в форме напоров

E ы p    av

— = H = z + — + g

(9.17)


pg    2g

Рассмотрим движение потока вязкой жидкости в канале (рис. 9.2) от сечения 1-1 к сечению 2-2. Обозначим удельную энергию потока в сечении 1-1 через E', а в 2-2 - E2.


Так как жидкость вязкая, то процесс ее перемещения сопро-

вождается диссипацией энергии, т.е. какая-то ее часть расходуется на преодоление сил внутреннего трения и превращается в тепло, следовательно, E2 < E'. Поэтому баланс энергии для выбранных сечений должен быть записан в виде

E' = E2 +Ае    (9.18)

где Аe - потери энергии.

Раскрывая значения E1 и E2, получаем:

22 Р1    1 V1    Р2    СС2У2    А    /Г\ АГ\\

gz1 + — +    = gz2 + — +    +    Ае    (9.19)

р 2 р 2 Это и есть энергетическая форма уравнения Бернулли для потока вязкой жидкости.

В практических приложениях чаще используют уравнение Бернулли, выраженное в напорах

22 Р1    ^1^1    Р2    a2w2    .1    .    ...

z1 + —    + -L-L = z2    + — +    —2-2 + Ah    (9.20)

1    pg    2g    2    pg    2g    K '

Ае

где — = Ah - потери напора. g

Для газовых потоков (без учета сжимаемости), а также при расчетах систем гидравлического привода обычно используют уравнение Бернулли в форме давлений

paV    pa2vj    .

pgz1 + Р1    = pgz2 + Р2 +    + Ap    (9.21)

где Ap - потери давления.

Обычно в упомянутых системах член pgz оказывается пренебрежимо малым по сравнению с остальными. В этих случаях (9.21) принимает вид:

Р1 + О- = Р2 + РТ2 + Ap    (9.22)

9.4 Физический смысл коэффициента Кориолиса.

Как уже упоминалось, коэффициент a носит название коэффициента кинетической энергии, корректива скорости, коэффициента Кориолиса. Выясним физический смысл этой величины.

Как уже отмечалось выше, второй член в уравнении (9.13) представляет собой кинетическую энергию секундной массы потока, определяемую истинным распределением скоростей в сечении, т.е.

ЕГ =р JJuZdA    (9.23)

A

Если бы скорости в сечении были бы распределены равномерно, то u = v = const (v - средняя скорость потока), и кинетическая энергия потока была бы

^=2    v'dA=2v3 ЯdA=Рг    (924)

A    A

Разделив (9.23) на (9.24), получим:

Следовательно, коэффициент Кориолиса представляет собой отношение кинетической энергии потока, вычисленной по истинному распределению скоростей, к кинетической энергии, определенной по средней скорости.

Для уяснения вопроса рассмотрим гипотетический «поток», состоящий из двух струек, скорости которых u1 = 2м/с и u2 = 4м/с и вычислим коэффициент Кориолиса.

Истинная кинетическая энергия (сумма кинетических энергий струек)

Еист = И2 + ul = 4 + 16 = 10 м2

и а =    *    = -9, т. е. а >1 (истинная кинетическая энергия больше

Ек

средней).

Легко убедится, что чем больше неравномерность распределения скоростей, тем больше коэффициент Кориолиса. Так, если

u1 = 2 м/с, а u2 = 6 м/с, то а = ^4. Очевидно, что минимальное значение а = 1 будет при равномерном распределении скоростей. Дей-

м2

ствительно, пусть u1 = u2 = 4 м/с, тогда Еист = 16-2- и

с2

м2

Е^р = 16—2. Следовательно, можно утверждать, что а корректирует

к    с2

ошибку, возникающую при вычислении кинетической энергии при замене истинного распределения скоростей условным равномерным.

Забегая несколько вперед, отметим, что в природе существует два принципиально отличающихся режима течения жидкости: ламинарный и турбулентный. При ламинарном течении в трубах ал = 2, при турбулентном ат = 1, 02. ..1,04. Это позволяет утверждать, что в турбулентном потоке скорости в поперечном сечении распределены существенно равномерней, чем в ламинарном (эпюра турбулентного потока более «наполненная», ближе к прямоугольной по сравнению с эпюрой ламинарного потока).

Подведем некоторые итоги. Использование струйной модели потока и сведение его к одномерному путем введения представления

о средней скорости позволяют получить одно из основных уравнений гидродинамики - уравнение Бернулли для потока вязкой жидкости. Принципиально, с помощью этого уравнения можно рассчитать движение жидкости в каналах при установившемся течении и условии, что в выбранных сечениях поток слабодеформированный либо параллельно-струйный. Однако, для полного решения задачи необходимо уметь определять потери напора (Ah), возникающие при движении жидкости в каналах. Эта далеко не простая задача и будет являться предметом дальнейшего рассмотрения.

10. КЛАССИФИКАЦИЯ ТЕЧЕНИЙ ЖИДКОСТИ. УСТОЙЧИВОСТЬ ДВИЖЕНИЯ.

В 80-х годах прошлого столетия работы, связанные с изучением сопротивления движению жидкости при течении в трубах, зашли в тупик. Опыты одних исследователей (немецкий инженер-строитель Г.Хаген, французский врач Ж.Пуазейль) показали, что сопротивление линейно зависит от скорости. В то же время не менее тщательные и точные опыты французского инженера А.Дарси свидетельствовали, что сопротивление пропорционально квадрату скорости. Возникшее противоречие тормозило развитие инженерной практики и требовало разрешения.

Наблюдения, выполненные Г.Хагеном еще в 1855 г. показали, что характер движения в трубе изменяется при достижении каких-то определенных условий. На это же со всей определенностью было указано в 1870 году нашим соотечественником проф. Н.Н.Петровым при разработке им теории гидродинамической смазки. Эта гипотеза нашла блестящее подтверждение в опытах английского физика Осборна Рейнольдса, результаты которых были опубликованы в 18831884 годах и имели далеко идущие последствия для всей механики жидкости.

Идея опытов отличалась ясностью и предельной простотой. В стеклянную трубу, скорость движения воды в которой могла регулироваться, Рейнольдс вводил струйки красителя. При малых скоростях струйки двигались параллельно оси трубы и вся картина представлялась неподвижной. При увеличении скорости воды за счет открытия крана картина изменялась, струйка красителя сначала приобретала синусоидальную форму, а дальнейшее увеличение скорости приводило к ее размыву, что свидетельствовало о беспорядочном движении.

n    w    w    w    w ^

Первый режим - спокойный, слоистый без перемешивания частиц был назван ламинарным. Второй - бурный, хаотичный, приводящий к перемешиванию частиц, позднее по предложению У. Томсона (Лорда Кельвина) получил название турбулентного. Как истинный ученый, Рейнольдс не остановился на констатации факта. Он предположил, что увеличении скорости потока приводит к возникновению каких-то возмущений, дестабилизирующих его структуру. Если понимать под устойчивостью способность потока подавлять возникающие в нем малые возмущения, то переход к турбулентному режиму может рассматриваться как потеря устойчивости. При этом из двух категорий сил, действующих на жидкие частицы, вязкого трения и инерции, первые играют стабилизирующую роль, а вторые - дестабилизирующую. Таким образом, отношение этих сил может служить критерием (мерой) устойчивости потока, т.е.

«    силы    инерции

Мера устойчивости =-—-

^ 7    силы    вязкого трения

Такой подход позволяет получить и количественную меру. Действительно, сила инерции F = та. Массу можно представить как произведение плотности на объем, но объем пропорционален кубу линейных размеров, т. е. т — pl . Ускорение есть изменение скорости в

единицу времени а = Ц. Таким образом

Fu„ — ^    (10.1)

По смыслу It есть скорость, следовательно,

Рин —pl2 и2    (10.2)

Сила вязкого трения (по Ньютону)

du dy

5тр =MdUS    (10.3)


Действуя аналогично предыдущему, получаем

Fтр М~и1 — М U 1 и безразмерный комплекс, характеризующий устойчивость, приобретает вид

^ = рМ-    (10.4)

Fтр    м

В дальнейшем это соотношение получило название числа Рейнольдса, т.е.

Re = рМ    (10.5)

М

где u - характерная скорость течения; l - характерный линейный размер.

Оригинальное толкование этого комплекса дано самим Рейнольдсом. Он писал: «Жидкость можно уподобить отряду воинов, ламинарное течение - монолитному походному строю, турбулентное -беспорядочному движению. Скорость жидкости - скорость отряда, диаметр трубы - величина отряда. Вязкость - дисциплина, а плотность - вооружение. Чем больше отряд, чем быстрее его движение и тяжелей вооружение, тем раньше распадается строй».

Для круглых труб характерный размер - диаметр, характерной скоростью является средняя скорость. С учетом этого, имея в виду,

что ^/р = v, выражение (10.5) принимает вид

Re = —    (10.6)

v

При течении в каналах некруглого сечения в качестве характерного размера принимают так называемый гидравлический радиус

A

R = —    (10.7)

П

где A - площадь поперечного сечения канала; П - смоченный периметр (часть периметра, находящаяся в контакте с жидкостью).

Для круглых труб при напорном движении A = П^, П = nd и

d

R = 4, т.е. гидравлический радиус в два раза меньше геометрического.

Одним из наиболее существенных результатов, обнаруженных в опытах Рейнольдса являлось то, что переход от ламинарного течения к турбулентному происходил при одном и том же численном значении введенного им критерия устойчивости, названного впоследствии критическим значением числа Рейнольдса (ReKp). По данным

многочисленных опытов в круглых трубах ReKp « 2300. Это так называемое нижнее критическое число Рейнольдса, которое получают, если не принимать специальных мер по стабилизации потока. При принятии мер, переход к турбулентному течению можно существенно затянуть. При выполнении технических расчетов принято считать, что если число Рейнольдса, вычисленное по фактическим значениям параметров, меньше критического, то режим ламинарный, и наоборот.

11. ЗАКОНОМЕРНОСТИ ЛАМИНАРНОГО РЕЖИМА ТЕЧЕНИЯ В КРУГЛЫХ ТРУБАХ.

При рассмотрении уравнений движения вязкой жидкости (уравнений Навье-Стокса) отмечалось, что интегрирование их в большинстве случаев связано с непреодолимыми математическими трудностями. Однако известны и исключения. К числу их относится ламинарное течение между параллельными пластинами, одна из которых движется с какой-то скоростью u. Это так называемое течение Куэтта.

Рассмотрение закономерностей этого течения можно найти в оригинально построенном современном американском курсе прикладной гидродинамики: Дейли Дж., Харлеман Д. Механика жидкости. -М.: Энергия, 1971. - 480 с.

Другим примером, интересующим нас в данном случае, является установившееся течение в круглой трубе, происходящее под действием постоянного перепада давлений - течение Пуазейля.

Профессор медицины Жан Пуазейль (1799-1869 гг.) во введении к своему трактату «Движение жидкостей в трубах малого диаметра» писал: «Я начал свои исследования потому, что прогресс в физиологии требовал определения законов движения жидкости в трубах малого диаметра (порядка 0,1 мм). Конечно, Дю Буа, Жирар, Навье и другие уже исследовали эти проблемы, однако они нуждаются в дальнейшем аналитическом и экспериментальном изучении, что было необходимо для надежного согласования теории с экспериментом». Опыты, выполненные Пуазейлем с трубкой диаметром 0,14 мм согласовывались с полученным им соотношением до тех пор, пока длина трубки составляла 51 мм; при уменьшении длины эта зависимость не соблюдалась. Этот факт и объясняется переходом от ламинарного к турбулентному режиму течения.

Как отмечалось выше, закономерности ламинарного течения в трубах можно получить путем прямого интегрирования уравнений Навье-Стокса. Решение задачи таким методом можно найти в книге: Аржаников Н.С., Мальцев В.Н. Аэродинамика. -М.: Изд-во оборонной промышл., 1956. -483 с.

В данном пособии мы воспользуемся другим способом, позволяющим получить более ясные физические представления.

Рассматриваем установившееся ламинарное течение в горизонтальной трубе, происходящее под действием постоянного перепада давления. Радиус трубопровода - R.

Двумя сечениями, отстоящими на расстоянии I друг от друга, выделим отсек трубопровода, и в нем цилиндр радиуса r. Составим уравнение движения. Так как течение установившееся, то сумма проекций на ось всех сил, действующих на цилиндр, должна быть равна нулю. Другими словами, активные силы, приводящие частицы жидкости в движение, должны быть равны силам сопротивления.


Активные силы: p1A - p2A = ApA = nr2Ap.

Силы сопротивления: 2nrI т.

Таким образом, nr Ap = 2nrI т и

Apr

(11.1)

т =


21


Из (11.1), в частности, следует, что касательные напряжения

Т =


Знак «минус» потому, что направления отсчета y и r противоположны.

Приравнивая (11.1) и (11.2), получаем

Ap r du

(11.3)

21 M'dr Либо после разделения переменных

Ap

du = -

r dr


(11.4)


2pI

и после интегрирования

Ap r2

+ C


(11.5)

u = -


2^1 2

Произвольную постоянную интегрирования находим из граничных условий: при r = R u = 0 (условие прилипания), и

Следовательно,


и = ^Р (R2 - r2)

4^Л    '


либо


ApR2

4^/


2


1 -


(11.7)


и =


R2


v


J


Максимальная скорость движения частиц будет на оси трубы, т.е. при r = 0, а ее величина


ApR"


(11.8)


и


max


4^/


Подставляя (11.8) в (11.7) получим


,2 Л


1


(11.9)


и = и


max


R2


V


J


Из чего следует, что в поперечном сечении трубы скорости распределены по параболическому закону, т.е. эпюра скорости представляет собой параболоид вращения.

Выражение (11.9) можно представить в виде


и


= 1 -


(11.10)


R2


и


max


Из чего следует, что отношение скорости в любой точке к скорости на оси не зависит от расхода, рода жидкости и материала стенок трубы: при всех значениях Re < Rew оно одинаково.

Определим расход, протекающий через трубопровод. При введении понятия о средней скорости было показано, что

R


Q = 2п j и (r)r dr


(11.11)


где и (r) - уравнение эпюры скорости. Воспользуемся (11.6), что дает Q = 2nApR


J(R2 - r2)


r dr


(11.12)


4^/


Выполнив интегрирование и имея в виду (11.8), можно полу


чить


Q = nu m2aX R = J u max A = VA Из чего следует, что

u max = 2V    (11.13)

Раскрывая значение umax по (11.8), получаем выражение для определения потерь давления при ламинарном режиме течения в круглой трубе

АР = ^    (11.14)

Либо, заменяя радиус диаметром,

АР = 32^2/V    (11.15)

Полученное соотношение носит название формулы Хагена-Пуазейля. Для потерь напора с учетом того, что Ар = pgAh, формула принимает вид

А h =    (11.16)

pgd

Важнейший вывод, следующий из этого соотношения, можно сформулировать так: потери давления (напора) при ламинарном течении в круглых трубах линейно зависят от средней скорости.

Выполним некоторые формальные преобразования формулы Хагена-Пуазейля, которые окажутся полезными в дальнейшем. Умножим числитель и знаменатель (11.16) на    2v, что дает

.    32 • 2v I V2    64v I V2    64 I V2

Ah =---=---=------(11.17)

vd d 2g vd d 2g    Red 2g

Таким образом, можем записать, что    в формуле    Ah =    f (vm)

при ламинарном течении m = 1.

12. ОСНОВНЫЕ ЗАКОНОМЕРНОСТИ ТУРБУЛЕНТНОГО ДВИЖЕНИЯ.

12.1. Общие сведения.

^    W    ^    W    W    W

Теория турбулентных течений представляет собой важнейший для практики, но и наиболее сложный раздел гидродинамики.

Как уже отмечалось, первые серьезные исследования перехода к турбулентности были выполнены О. Рейнольдсом в 1883 году.

Им же со ссылкой на Стокса был предложен ответ: «Общей причиной изменения стационарного течения на завихряющееся является то обстоятельство, что при некоторых условиях стационарное движение становится неустойчивым, так что бесконечно малые возмущения могут привести его к переходу в волнистое движение». «Волнистое движение», так первоначально было названо турбулентное движение Рейнольдсом. К сожалению, исследование бесконечно малых возмущений не дало критических значений, близких к наблюдавшимся в опытах.

Основной, определяющей чертой турбулентного движения является его хаотичность. Это означает, что скорость (и другие параметры) в любой точке потока зависят от времени. Более того, эти флуктуации скорости в данной точке также являются хаотическими.

Подробный исторический обзор развития теории турбулентности можно найти в капитальном двухтомном труде известных советских специалистов А.С. Монина и А.М. Яглома «Статистическая гидромеханика» (ч.1. -М.: Наука, 1968. -639 с.)

В настоящем пособии мы ограничимся лишь самыми общими сведениями, в какой-то мере поясняющими сложные и еще не до конца понятые вопросы, связанные с турбулентным движением.

Впервые гипотеза о физическом механизме турбулентного перемешивания была высказана английским ученым Л. Ричардсоном в 1922 г. Условно турбулентное движение принято рассматривать как совокупное движение отдельных структур, называемых молями либо вихрями, совершающими как поступательное, так и вращательное движение. По Ричардсону развитая турбулентность представляет собой иерархию «вихрей». При зарождении вихри имеют большие размеры, соизмеримые с размерами канала. Затем за счет потери устойчивости они распадаются на более мелкие, передавая при этом им свою энергию. Возникает каскадный процесс, в котором энергия осредненного потока последовательно передается вихрям все более мелких масштабов. В конечном итоге образуются вихри минимального масштаба, которые далее не разрушаются. При этом нижний размер вихря (турбулентного образования) определяется вязкостью среды. В самых малых вихрях кинетическая энергия турбулентности за счет сил вязкого трения превращается в тепло, т.е. происходит диссипация энергии. Это указывает на необратимый характер процесса.

Из сказанного ясно, что турбулентное движение по своей фи-

W    Л    W

зической природе является движением неустановившимся. С другой стороны, непосредственные измерения свидетельствуют, что при турбулентном характере потока в нем можно выделить основную, так называемую регулярную часть, на которую накладывается случайная часть движения.

На рис. 12.1 показан типичный вид экспериментально снятой зависимости проекции скорости в какой-то точке потока от времени при сохраняющихся неизменными граничных условиях.

Как следует из графика, особенностью этого процесса является его непериодичность, при этом


U'x = Ux - (Ux ), где (ux^ - осредненная скорость, представляющая регулярную часть; и'х - пульсаци-онная скорость, разность между мгновенным и регулярным значением скорости.

Аналогичные соотношения можно записать и для других компонент.

Таким образом, осредненная скорость - это какое-то устойчивое значение, вокруг которого происходит изменение рассматриваемой проекции скорости (в данном случае ). Все сказанное в равной мере относится и к другим параметрам, в частности, к давлению.

Наиболее важной характеристикой течения при его расчете является поле скоростей. Но, как показано выше, в любой точке потока при турбулентном течении скорость выступает как случайная величина, что исключает возможность записи начальных условий для системы дифференциальных уравнений Навье-Стокса, т.е. оказывается невозможной математическая постановка задачи. Именно это и приводит к необходимости перехода к какому-то осредненному описанию, использующему не истинные, а осредненные величины скоростей и давлений. Осреднение скоростей и давлений производится путем интегрирования функций их (х, y, z, t), uy(х, y, z, t), uz(x, y, z, t),

p (x, y, z, t) по промежутку времени T (см. рис. 12.1), величина которого намного больше так называемого характерного времени турбулентных пульсаций. Это время определяется как частное от деления масштаба / на скорость турбулентных пульсаций. Под масштабом турбулентных пульсаций понимают расстояние, на котором пульсации претерпевают заметное изменение. Так, например, при турбулентном движении в трубах наибольший масштаб пульсаций равен диаметру трубы. Таким образом, осредненная компонента скорости, например, (ux)

t0

Аналогичное соотношение можно записать и для давления. При этом, поскольку флуктуации (пульсации) имеют как положительный так и отрицательный знак, то

1 ^0+T

(ux) = - Ju'x(t)dt = 0    (12.2)

t0

Ясно также, что ^u'x2^ Ф 0. Если в данной точке потока

(uх2^ = (u'y2^ = (u'z2), то турбулентность называют изотропной, а

если это условие соблюдается во всех точках, то она называется еще и однородной.

12.2. Уравнения Рейнольдса.

Как уже отмечалось, сложность турбулентного движения делает невозможным строгое рассмотрение течений при заданных граничных условиях. Одной из возможных альтернатив является переход от истинной картины, детали которой нам неизвестны, к рассмотрению осредненного турбулентного течения, т.е., по существу, замена принципиально неустановившегося движения на квазиустановив-шееся. Этот переход был предложен О.Рейнольдсом. Суть его сводится к тому, что в уравнениях движения вязкой жидкости (уравнениях Навье-Стокса) и уравнении неразрывности истинные значения параметров по определенным правилам заменяются их осредненными значениями. Получаемая таким образом новая система уравнений носит название уравнений Рейнольдса. Вывод этих уравнений выходит за рамки настоящего курса. Интересующиеся могут найти его в ряде учебных пособий, в частности, Федяевский К.К., Войткунский Я.И., Фаддеев Ю.И. Гидромеханика. - Л.: Судостроение, 1968. - 567 с.

Наиболее существенным результатом этой операции является то, что вследствие нелинейности уравнений Навье-Стокса в уравнениях Рейнольдса появляются дополнительные члены, которые получили название напряжений Рейнольдса. Для наиболее простого плоскопараллельного течения эти напряжения имеют вид:

тт = р( u х uy)    (12.3)

(угловые скобки - символ осреднения).

Таким образом, в осредненном турбулентном потоке к обычным вязкостным напряжениям добавляются напряжения, зависящие от пульсации скорости. Физически это объясняется тем, что между

разными участками турбулентного потока происходит обмен количеством движения, обусловленный перемешиванием частиц. Перенос количества движения вызывает дополнительное торможение либо ускорение отдельных масс жидкости, т.е. приводит к возникновению турбулентных напряжений.

Поскольку исходная система уравнений являлась замкнутой (четыре уравнения и четыре неизвестных - ux, uy, uz, p), то появление дополнительных членов в уравнениях Рейнольдса приводит к тому, что она превращается в незамкнутую. Возникает новая проблема «замыкания системы уравнений Рейнольдса».

12.3. Полуэмпирические теории турбулентности.

Современная теория турбулентности не располагает возможностями теоретическим путем получить уравнения для определения напряжений Рейнольдса. Поэтому единственным способом, позволяющим замкнуть систему, является привлечение полуэмпирических соотношений, связывающих эти напряжения с осредненными по времени компонентами скорости (^ux), (uy^ и (uz).

Один из первых исследователей турбулентности, Ж.Буссинеск, предложил выражать турбулентные напряжения аналогично закону трения Ньютона, т.е.


(12.4)

где п -турбулентная вязкость.

В отличие от физической, турбулентная вязкость характеризует не физические свойства жидкости, а статистические свойства пульсационного движения. Поэтому она не является постоянной величиной, а может изменяться как в пространстве, так и во времени. Важно также отметить, что даже на небольших удалениях от твердых границ турбулентная вязкость существенно превосходит физическую

(П >> м).

В целом для турбулентного потока можно записать

Однако представление Буссинеска не приводит к решению задачи, т.к., к сожалению, отсутствуют прямые методы определения турбулентной вязкости.

Первого заметного успеха в этом направлении добился Л.Прандтль в 1925 году, предложив так называемую теорию пути перемешивания (смешения).

n    w    w

В основе ее лежит аналогия с кинетической теорией газов и предположение о том, что путь смешения зависит от условий течения. В соответствии с гипотезой Прандтля, каждый турбулентный моль (вихрь) жидкости переносит некоторое количество движения, которое сохраняется постоянным на пути перемешивания. Другими словами, длина пути перемешивания в известной мере аналогична длине свободного пробега молекул в кинетической теории газов, и определяет путь, который проходит моль жидкости, прежде чем он перемешается с другими жидкими образованиями и передаст свой импульс.

Допустив далее, что вертикальная и горизонтальная компоненты пульсационной скорости (u'x и иy) являются величинами одного

порядка, Прандтль получил формулу для определения турбулентного напряжения в виде

Тт = Р/" f dy]    ('26)

V dy J

где /п - длина пути перемешивания.

Угловые скобки вокруг и, символизирующие операцию осреднения, для упрощения записи опущены.

Интересующиеся выводом формулы Прандтля, могут найти его в книгах: Аржаников Н.С., Мальцев В.Н. Аэродинамика. - М.: Изд-во оборонной промышл., 1956. - 483 с., либо Шлихтинг Г. Теория пограничного слоя. - М.: Наука, 1974. - 711 с.

На первый взгляд может показаться, что формула Прандтля не имеет каких-либо существенных преимуществ по сравнению с формулой Буссинеска, и единственным результатом является замена одной не поддающейся вычислению величины п другой - /п. Однако это не так, поскольку величину /п оценить значительно проще, чем п. В частности, /п не может быть больше размера канала и должна стремиться к нулю вблизи твердой стенки (поперечное движение у стенки невозможно).

12.4. Турбулентное течение в трубах.

Расчет турбулентного течения в трубах относится к широко распространенным инженерным задачам. Одним из важных элементов расчета является нахождение закона распределения осреднен-ных скоростей в поперечном сечении трубы.

По Прандтлю, поток в трубах при турбулентном течении условно разбивается на две области (двухслойная модель Прандтля): турбулентное ядро, в котором определяющими являются напряжения

м    W    WWW/    W    W    п

Рейнольдса, и тонкий вязкий подслой (ламинарный подслой по Пран-

дтлю либо пристенный слой) вблизи стенки, в котором влияние турбулентности пренебрежимо мало, а касательные напряжения обусловлены физической вязкостью в соответствии с законом трения Ньютона.

На рис. 12.2 приведен примерный вид поля осредненных скоростей (эпюра скорости) при турбулентном течении в трубопроводе. Следует обратить внимание на ее большую наполненность (большую равномерность) по сравнению с ламинарным течением. Это объясняется тем, что вследствие перемешивания частиц за счет турбулентных пульсаций происходит обмен количеством движения и, как следствие, более равномерное распределение скоростей в поперечном сечении.

В непосредственной близости от стенки в пределах пристенного слоя решающее влияние на течение оказывают жесткость стенки, ее непроницаемость и эффект прилипания частиц. На самой стенке справедливы условия:


ux = uy = uz = 0;

u'x = uy = uy = 0; / u'x uy) = 0.

Таким образом, для области в пределах вязкого подслоя можно записать:

t= T0 = p

du

dy


(12.7)


где t0 - касательное напряжение на стенке. Интегрирование (12.7) дает

pu = yT0 + C при y = 0, u = 0 и C = 0. Таким образом,

(12.8)


Имея в виду, что p


(12.9)


To

u =y

p

vp, после подстановки получаем

T y


P v


u =


Из чего следует, что в пределах подслоя скорость изменяется по линейному закону. Величина Тт/ имеет размерность квадрата скорости, поэтому корень квадратный из нее, т.е.

1ТТ = uT    (12.10)

Р

называют динамической скоростью либо скоростью трения. Из выражения для напряжений Рейнольдса (см. 12.3) следует, что

= pu'xuy) и

тт

p ^u'x u'y) = ut    (12.11)

Таким образом, динамическая скорость является мерой интенсивности турбулентного пульсационного движения, т.е. мерой интенсивности переноса количества движения.

Подставляя (12.10) в (12.9), получаем

u = uT —    (12.12)

Оценим толщину вязкого подслоя. На его границе y = 8, и

(12.12) можно придать вид

u = ut8

(12.13)

uT    V

В правой части стоит выражение, аналогичное числу Рейнольдса. Согласно тщательным опытам ближайшего сотрудника Л.Прандтля, Никурадзе, эта величина приближенно равна 11,6; тогда

8= 11,6    (12.14)

u Т

Очевидно, что этим соотношением можно воспользоваться лишь в случае, если известна динамическая скорость. Для ее нахождения необходимо увязать ее с параметрами осредненного потока, что является решаемой задачей.

Чтобы завершить вопрос о турбулентном течении в трубах, установим закон распределения осредненных скоростей в ядре потока. В этой области определяющую роль играют турбулентные касательные напряжения, и, следовательно, можно воспользоваться формулой Прандтля (см. 12.6). Однако для того, чтобы продвинуться дальше, необходимо принять дополнительные допущения. Они оказываются достаточно грубыми, и единственным их оправданием является то, что результаты, к которым они приводят, достаточно хорошо согласуются с экспериментальными данными.

Первое допущение связано с длиной пути перемешивания. Согласно наиболее простой гипотезе, принадлежащей Л.Прандтлю,

ln = кУ    (12.15)

где к - какая-то величина, называемая постоянной Кармана. Выполненные измерения показывают, что к « 0,4. Более поздние исследования показали, что зависимость (12.15) справедлива лишь в пристенной части турбулентного ядра потока.

Вторым является допущение о касательных напряжениях. Следует полагать, что принципиально они являются величинами переменными. Однако, если рассматривать область, расположенную достаточно близко к стенке, то здесь величина касательного напряжения изменяется незначительно, и можно принять ее равной касательному напряжению на стенке, т.е. тт = т0.

При этих допущениях формула Прандтля принимает вид

либо

2 2f du]2


y -Т

V dy J

Извлекая квадратный корень и разделяя переменные, получаем

и после интегрирования

(12.16)

т.е. скорости в ядре потока распределены по логарифмическому закону.

Произвольную постоянную интегрирования можно найти из граничных условий на оси трубы: при y = R и = иmax, и

C = иmax - — In R. После подстановки и простых преобразований

- и    1    R

max-- = - In—    (12.17)

U т    К    y

Строго говоря, соотношение (12.17) выводится для плоских труб, но опыт показывает, что оно оказывается справедливым и для круглых, и подтверждает экспериментально установленный факт о независимости распределения скорости от причин, обусловливающих возникновение касательных напряжений (вязкости, шероховатости).

Выражение (12.17) иногда называют законом дефекта скорости.

Использование двухслойной модели, т.е. разделение потока на ядро и пристенный слой, приводит к специфической классификации стенок труб. Если толщина пристенного слоя больше выступов шероховатости, трубы называют гидравлически гладкими, в противном случае - шероховатыми.

Завершая раздел, обратим внимание на следующие обстоятельства. Как отмечалось, для получения закона распределения скоростей в поперечном сечении трубопровода использовались простейшие гипотезы: постоянство касательных напряжений в ядре потока т = т0) и линейная зависимость для длины пути перемешивания (1п = Ky). Легко показать, что первая из них не согласуется с реальностью при рассмотрении течения в трубах. Действительно, выделим в трубе цилиндрический элемент жидкости длиной I и радиусом г, на который действует постоянный перепад давления Ар. Сила

давления на этот элемент Арп r , а сила трения 2пг1т. Приравняв эти силы, получаем

2I

Ар = — т    (12.18)

А для всей трубы длиной I и радиусом R

Ар = R T0    (12.19)

где т0 - напряжение на стенке.

Поскольку Ар = const по условию, то приравняв (12.18) и (12.19), с учетом того, что r = R - y,

г    y\

1 - y

(12.20)


т = т0

'V    RJ

т.е. касательные напряжения по сечению не постоянны, а изменяются по линейному закону, и лишь на достаточно малом расстоянии от

стенки (RR << 1) можно считать, что т = т0.

Вторая гипотеза также не согласуется с данными опытов. На рис. 12.3 приведены графики, характеризующие распределение длины пути перемешивания в поперечном сечении круглой трубы по данным опытов Никурадзе (кружки) и по формулам, предложенным различными авторами. В соответствии с результатами экспериментов, значение 1п достигает максимума на оси трубы. Из графика следует, что гипотеза Прандтля (прямая 1) неприемлема.

Существенно отличаются от опытной и кривые, полученные другими авторами: Карманом (кривая 2), Конаковым (кривая 4), Саткевичем (кривая 5). Достаточно близка к эксперименту кривая Альтшуля (кривая 3), описывающая длину пути перемешивания с помощью формулы

2у

Рис. 12.3


у 1

2

R )

(12.21)


I n =к1

R

В последнее время Д.Н.Васильевым получена аппроксимирующая зависимость, практически точно совпадающая с данными опыта и имеющая вид

3

к


ln =


R


(12.22)


1 - 1 -


R )


Использование этого соотношения с учетом линейного распределения касательных напряжений по сечению трубы приводит к закону распределения скоростей, соответствующему гиперболическому тангенсу. Вывод этого соотношения можно найти в книге: Павленко В. Г. Основы механики жидкости. - Л.: Судостроение, 1988. - 240 с.

Существуют и другие подходы к этому непростому вопросу. Так, например, А.Д.Альтшуль считает, что разделение потока на две области является грубой схематизацией, носящей искусственный характер. Не оправдана с теоретических позиций гипотеза о ламинарном подслое, как об области, в которой отсутствуют пульсации. Пульсации проникают и в этот слой, но следуют там особым закономерностям. Слабо обосновано и то, что в ядре потока физическая вязкость не играет никакой роли. На базе этих представлений автором разработана полуэмпирическая теория, рассматривающая турбулентный поток в трубе как единое целое, без разделения его на ядро и ламинарный подслой. Достаточно ясное и подробное изложение этой теории можно найти в книге: Альтшуль А. Д. Гидравлические сопротивления. - М.: Недра, 1970. - 215 с.

Полуэмпирические теории неоднократно подвергались серьезной критике. Главные возражения обычно касались выводов, связанных с особенностями структуры турбулентности. Тем не менее, они широко распространены из-за их простоты и удобства, хотя получаемые результаты достаточно грубы и приближенны. Поэтому нельзя не согласиться с мыслью, высказанной крупнейшим специалистом в области механики жидкости Х. Раузом (Механика жидкости. - М.:

Стройиздат, 1967. - 390 с.): «При сравнении простоты соотношений со сложностью явления вызывает удивление степень их полезности, а не их недостатки».

12.5. Степенные законы распределения скоростей.

Логарифмический закон распределения скоростей хорошо подтверждается результатами эксперимента, но вызывает определенные трудности при численных расчетах. Поэтому в последние годы получают распространение степенные зависимости вида

'n


и


U max


V R J

Главным достоинством этих формул является их простота, а


(12.23)


ry\


недостатком - зависимость показателя степени от числа Рейнольдса. Поэтому степенной закон нельзя рассматривать как универсальный.

3    6

В диапазоне изменения чисел Re = 4 • 10 ...3 • 10 показатель степени 1/n меняется в пределах от 1/6 до 1/10.

Следует отметить, что ни логарифмический, ни степенной законы не удовлетворяют условию равенства нулю производной от скорости на оси симметрии потока.

12.6. Потери давления (напора) при турбулентном

течении в трубах.

Напомним, что рассмотрение закономерностей как ламинарного, так и турбулентного течений в трубах помимо чисто познавательных целей преследовало и цели сугубо практические: получить соотношения, позволяющие определять потери давления (напора) в трубопроводных сетях при выполнении инженерных расчетов. Для ламинарного течения эта задача решается с помощью формулы Хаге-на-Пуазейля. Из рассмотрения закономерностей турбулентного течения становится ясным, что вследствие его чрезвычайной сложности получение аналогичного соотношения чисто теоретическим путем практически невозможно. Поэтому, основываясь на уже известных положениях, установим хотя бы общую структуру необходимой формулы.

Как было показано, выражение для турбулентных касательных напряжений (напряжений Рейнольдса) имеет вид

Т=р их и ’y)

Это с большой долей уверенности позволяет утверждать, что существует связь между средней скоростью и касательным напряжением на стенке трубы вида

т0 = kpv2    (12.24)

где k - коэффициент пропорциональности.

С другой стороны, из условия равновесия движущегося под действием постоянного перепада давления жидкого цилиндра длиной

I (см. 12.19)

2I

Ар = R Т0

После замены радиуса диаметром и подстановки т0

I 2

Ар = 4k ^ pv    (12.25)

либо

I pv2

Ар = 8k    (12.26)

pv 2

В такой форме записи выражение имеет четкий физический смысл. Это так называемое динамическое давление потока, обусловленное средней скоростью, либо кинетическая энергия потока, заключенная в единице объема.

Обозначим величину 8k = Я и назовем ее гидравлическим коэффициентом трения, тогда

, 2 I pv

Ар = Я-^-—    (12.27)

d 2

либо

I v2

А^ = Я----(12.28)

d 2g

Полученное соотношение носит название формулы Дарси. Более строго это соотношение будет получено методом анализа размерностей.

Отметим попутно, что если в преобразованной формуле Хаге-

64

на-Пуазейля (см. 11.17) обозначить величину буквой Я, то она

превращается в формулу Дарси. В этом смысле формула Дарси может быть названа универсальной, т.е. пригодной как для ламинарного, так и для турбулентного течений. В последнем случае открытым остается вопрос о нахождении гидравлического коэффициента трения, который, как следует из всего сказанного выше, может быть решен только экспериментальным путем.

В заключение отметим, что хотя поставленная главная проблема и оказалась теоретически неразрешимой, полученные результаты позволяют найти решения ряда частных задач, имеющих важное практическое значение.

13. ОСНОВЫ ТЕОРИИ ПОДОБИЯ И МОДЕЛИРОВАНИЯ

При рассмотрении различных разделов, связанных с движением жидких сред, неоднократно приходилось сталкиваться с процессами и явлениями, которые в силу своей сложности не позволяют получить аналитические решения, необходимые для инженерной практики. Вместе с тем переход от качественных суждений к количественным соотношениям играет ведущую роль в творческой деятельности человека.

Рассматриваемые в настоящем пособии вопросы непосредственно связаны с методологией научного познания. Однако, этот аспект, безусловно важный с познавательных позиций, далеко выходит за рамки курса, поэтому в настоящем пособии мы ограничимся лишь технической стороной.

Принципиально , процесс познания человеком природы можно условно разделить на две стадии: анализ и синтез. На первой стадии, т.е. на стадии анализа, изучаемый объект мысленно расчленяется на более простые составные части, выделяются свойства и связи.

На этапе синтеза происходит их соединение с целью воссоздания единого целого. Этап завершается построением математической модели, которая с какой-то степенью приближения описывает поведение изучаемого объекта. Обычно математическая модель представляет систему либо системы дифференциальных уравнений. Что же касается степени приближения модели, то она обусловлена теми упрощающими предпосылками, которые положены в основу. Здесь важную роль играет так называемый фактор неопределенности. Суть его сводится к тому, что с усложнением математической модели за счет более полного учета влияющих факторов уменьшается возможность получения точного, имеющего практическое значение представления. Другими словами, неопределенность решения возрастает по мере углубленного анализа реальной задачи.

Так, например, система дифференциальных уравнений Эйлера для гидродинамики является математической моделью, описывающей движение идеальной жидкости. Усложнение модели за счет учета сил вязкого трения приводит к системе дифференциальных уравнений Навье-Стокса.

Если модель разрешима, т.е. уравнения могут быть проинтегрированы любым путем, то можно считать, что решена и поставленная конкретная задача. Полученные результаты сопоставляются с теми, что наблюдаются в природе. Если они близки, то это означает, что модель правильно отражает поведение и свойства реального объекта, если нет, нужно ввести какие-то дополнительные факторы, не учтенные ранее, т.е. улучшить ее. Все это, конечно, не означает, что этот процесс идет легко и просто. Он может быть связан с преодолением огромных трудностей как математического, так и вычислительного характера. Новые проблемы возникают в двух случаях: несмотря на все усилия уравнения, составляющие математическую модель, проинтегрировать не удается; изучаемое явление оказывается столь сложным, что не поддается математическому описанию.

В качестве примера первого случая можно привести уравнения Навье-Стокса, которые не могут быть проинтегрированы для большинства важных для практики случаев. Очевидно, что единственным в этих условиях способом решения задачи является эксперимент на физической модели, под которой понимается уменьшенный (либо увеличенный) реальный объект исследования. При этом сразу возникают три вопроса: как спроектировать и построить модель, какие величины необходимо измерять при проведении опытов, и как перенести результаты опытов, полученных на модели на натурный объект. На эти вопросы и отвечает теория подобия, являющаяся основой современного физического эксперимента. Прежде чем приступить к в ее рассмотрению, необходимо уяснить, что же понимается под подобием? Одно из наиболее удачных определений этого понятия принадлежит академику Л.И.Седову: «Подобными называются такие явления (процессы), когда по характеристикам одного из них можно получить характеристики другого простым пересчетом, аналогичным переходу от одной системы единиц к другой».

В общем случае различают три вида подобия: геометрическое, кинематическое и динамическое. Наиболее простым является подобие геометрическое, требующее, чтобы линейные размеры натуры и модели находились в постоянном соотношении, другими словами, модель повторяет натуру в каком-то масштабе.

Это требование можно записать в виде

L м _ к L L

*“ н

где kL - масштабный множитель.

Для площадей (S) и объемов (V)

_ к2; _ к3 S L' V L

нн

Можно отметить, что правила геометрического подобия были известны еще Джонотану Свифту, который отмечал, что в стране, в которую попал Гулливер в одном из своих путешествий, он

обнаружил существа, превосходившие его по росту в 12 раз, по площади - в 144 раза и по объему - в 1728 раз.

Применительно к физическим явлениям элементарные представления геометрического подобия расширяются и распространяются на все величины, характеризующие данный процесс. Если учесть, что они могут изменяться как во времени, так и в пространстве, образуя поля, то возникает понятие о временном подобии и подобии полей, называемое кинематическим подобием.

В механике жидкости оно сводится к подобию полей скоростей в потоках, движущихся в геометрически подобных каналах.

И наконец, имея в виду, что механическое движение происходит под действием сил, вводится понятие динамического подобия, которое требует, чтобы в соответствующих точках натуры и модели силы находились в постоянном соотношении.

Рассмотрим простейший пример. Известно, что движение любой механической системы подчиняется закону Ньютона

F = md-    (13.1)

Для двух подобных систем можно записать

du1    du2

F - m dT и F2m2 t

Разделив первое на второе получим:

F1 m1 du1 dt2    _ F1 m1 u1 t2

— - —112 либо — - —1L

F2 m2 dt2 dt1    F2 m2 t2 t1

Имея в виду, что m - pV = pL имеем

F1 p1L1U1t2

F2    p2L2U 2t1

По смыслу L/t есть скорость, поэтому

F1 p1L1U1 F2    p2L2 U 2

(13.2)


либо

F1    -    F2

(13.3)


2 2 2 2 p1L 1U1 p2L 2 U 2

Очевидно, что полученные комплексы безразмерны. Таким образом, для двух подобных систем сохраняется

F_

числовое равенство безразмерных комплексов f 2 2 . Кратко это


pL2U ‘

условие можно записать так:    ,    2    2    =    idem    . В честь Ньютона этот

pL2u 2

комплекс обозначается двумя первыми буквами его фамилии, т.е.

Ne =-^~2    (13.4)

pL u

и называют числом подобия Ньютона, а выражение Ne = idem -основным законом динамического подобия механических систем (законом Ньютона).

Величины L и и, входящие в (13.4), называются определяющим линейным размером и определяющей скоростью. При проведении опытов они выбираются экспериментатором произвольно, исходя из удобства их измерения.

Полученные результаты заслуживают того, чтобы остановиться и сделать кое-какие полезные выводы. Во первых, они позволяют ответить на один из поставленных выше вопросов: как спроектировать и построить модель. Ответ очевиден: так, чтобы она была геометрически подобна натуре.

Во-вторых, из сказанного следует, что для обеспечения динамического подобия не требуется, чтобы все величины, определяющие характер процесса в натурном объекте, были численно равны аналогичным величинам в модели. Достаточным является равенство безразмерных комплексов, составленных из этих величин для натуры и модели, называемых числами подобия.

Какие преимущества дает такой подход в практическом плане?

Из математической статистики известно, что число опытов, которое необходимо поставить для того, чтобы получить закономерность, достоверно описывающую какое-то физическое явление, определяется из соотношения:

N = crk    (13.5)

где с - число экспериментальных точек, которое необходимо снять

для обеспечения представительности опыта (cmin = 5); k - число величин, подлежащих варьированию в опытах.

Таким образом, минимальное число опытов

N = 5k    (13.6)

Если в опытах варьируется число Ньютона (например, за счет изменения скорости), то k = 1 и N = 5, но если изучать влияние каждой из величин (p, и, L), то k = 3 и число опытов N = 125. Следовательно, использование числа подобия в качестве своеобразной «обобщенной переменной» позволяет уменьшить

число необходимых опытов в 25 раз, а если для надежности принять а _ 10, то в 100 раз.

И наконец, в-третьих, можно ответить на вопрос о том, какие величины следует измерять в опытах и как переносить результаты на натурный объект. Так как при проведении опытов необходимо обеспечить равенство чисел подобия натуры и модели, то ясно, что измерению подлежат лишь те величины, которые входят в эти числа.

По результатам измерений можно вычислить числа подобия модели и, исходя из равенства их числам подобия натуры, произвести пересчет.

Остается открытым вопрос, который, по существу, является центральным. Как же найти числа подобия, характеризующие изучаемый процесс либо явление? Очевидно, что только ответ на него открывает путь для практической реализации теории подобия.

13.1. Инспекционный анализ дифференциальных

уравнений.

Исходим из того, что математическая модель процесса нам известна, но она не может быть проинтегрирована. В этом случае числа подобия могут быть найдены методом, который по предложению известного американского математика и гидродинамика Г.Биркгофа назван инспекционным анализом. Как следует из названия, метод заключается в организованном по определенным правилам «инспектировании» дифференциальных уравнений, которое должно выявить числа подобия, позволяющие моделировать процесс. Отметим лишь, что этот метод не является единственным. Интересующиеся другими подходами могут найти их в книге Я.М. Брайниса «Подобие и моделирование в химической и нефтехимической технологии».-М.: Гостоптехиздат, 1961. - 219 с.

Базой инспекционного анализа является положение, рассматриваемое как постулат и сводящееся к следующему.

Если две системы описываются одинаковыми дифференциальными уравнениями и имеют одинаковые граничные условия, и если значения всех параметров в этих уравнениях и граничных уравнениях равны, то эти две системы подобны, при условии существования единственности решения.

Желающие познакомиться с подробным анализом всех элементов, входящих в это утверждение могут обратиться к превосходной книге А.А.Гухмана «Введение в теорию подобия». - М.: Высшая школа, 1963. - 253 с.

В данном пособии мы примем его как постулат, опуская все обоснования. Вторая его половина относится к так называемым условиям однозначности. Это крайне важное понятие требует более внимательного рассмотрения.

Любое дифференциальное уравнение описывает целый класс явлений, т.е. решение их многозначно. Так, например, то же уравнение Навье-Стокса, к которому мы уже неоднократно возвращались, может описывать движение жидкости в каналах, реках и океанах, движение атмосферных масс воздуха и т.п. Инженера интересует конкретное явление данного класса. Поэтому из множества возможных решений требуется лишь одно, соответствующее изучаемому явлению. Этого можно добится, если при постановке задачи ввести дополнительные так называемые условия однозначности, которые включают:

-    данные о физических свойствах среды (плотность, вязкость);

-    сведения о начальном состоянии системы (начальные условия);

-    данные о поведении системы на её границах (граничные условия).

Инспекционный анализ представляет собой определенный алгоритм, включающий два этапа: на первом из них отношение дифференциальных величин заменяются отношениями самих переменных, на втором - уравнение приводится к безразмерному виду путем деления всех его членов на один из них, выбранный произвольно.

Метод наиболее просто усвоить, обратившись к рассмотрению конкретного примера. Имея в виду, что в механике жидкости основными соотношениями, описывающими движение вязких сред, являются уравнения Навье-Стокса, целесообразно воспользоваться именно ими. Рассмотрим одну из проекций в декартовой системе координат. В данном случае безразлично какую, так как структура уравнений одинакова, что обеспечит и одинаковость получаемых результатов.

В проекции на ось х-ов имеем

дих д их дих    дих    1 др

+ v


v дх ду    д^Т j

Будем считать, что из массовых сил действует только сила тяжести, т.е. X = g cos в (cose учитывает знак). С учетом этого и после умножения всех членов уравнения на плотность получим

dp

dx

д

тр


+


В такой форме записи каждый из членов выражает силу, отнесенную к единице объема. При этом

Fu1 и Fu 2 - силы инерции;

Fm - сила тяжести;

Fd - сила давления;

F р - сила вязкого трения.

Действуя по алгоритму, заменим дифференциальные соотношения отношениями величин. Имеем:

2

F pU; F pU; F пег

Fu1    4.    ;    Fu 2    /    ;    — PQ;


t

с Ар    

Fd — —; f р — vP


dx


Приводим эти соотношения к безразмерному виду, приняв в качестве делителя один из комплексов. Как отмечалось выше, он может быть выбран произвольно. Пусть им будет Fu 2, т.е. силы инерции. Получаем:

Fu 1 = L = Sh - это так называемый критерий гомохронности


tu


u2


либо число подобия Струхаля.


L

_d_

dx


u

V P —2


V


У


2

QL ^    и    г

т

u2


2, обратная величина — = Fr - число Фруда

и    QL

отношение сил инерции к силам тяжести.

Ар


д


= Еи - число Эйлера, отношение сил давления к


Fu2    P u

силам инерции.


v ^    uL    _

обратная величина — = Re - уже известное нам


т р


uL


v


u 2


число Рейнольдса - отношение сил инерции к силам вязкого трения.

Следует отметить, что вопрос о правильности интерпретации чисел подобия как отношения сил ставился рядом исследователей.

По некоторым сведениям, еще Прандтль высказывал мнение о том, что число Рейнольдса не всегда равно отношению силы инерции к силе внутреннего трения. Более точным и правильным является утверждение, что если две системы геометрически подобны и течение в них происходят при одинаковых числах Рейнольдса, то отношение сил инерции к силам трения для обоих потоков одинаково. Подробности, связанные с таким подходом можно найти в книге С. Клайна «Подобие и приближенные методы». - М.: Мир, 1968. - 302 с.

Таким образом, при моделировании гидромеханических явлений необходимо использовать числа подобия Струхаля, Фруда, Рейнольдса и Эйлера.

Анализируя величины, входящие в числа подобия, легко заметить, что они составлены из параметров, входящих в условия однозначности. Эти числа подобия называют определяющими. Экспериментатор, разумеется, в определенных пределах, может изменять их величину ( менять скорость, геометрические размеры, вязкость). В число Эйлера входит величина Ар - перепад давления (потеря давления), которая, как правило, является искомой. Другими словами, величина числа Эйлера является следствием (результатом) процесса. Числа подобия такого рода называются неопределяющими. С чисто математических позиций сказанное можно представить в виде

Eu - f (Sh, Fr, Re)    (13.7)

Если изучается установившееся движение, при котором параметры в точке не изменяются с течением времени, то из рассмотрения выпадает число Струхаля и

Eu - f (Re, Fr)    (13.8)

Следовательно, при моделировании гидромеханических явлений в данном случае должны соблюдаться следующие условия, обеспечивающие динамическое подобие:

Fr„ - FrH; ReM - ReH

либо

UM    _    UH . UML M    UHL H

(13.9)

дмL м    днL H    VM    VH

Если при проведении опытов удается соблюсти эти требования, то подобие называется полным. Однако в реальных условиях добиться этого достаточно трудно, а иногда и просто невозможно. Поэтому обычно ограничиваются частичным подобием. Анализируя сущность явления экспериментатор устанавливает какие из сил (тяжести,

трения) играют определяющую роль в исследуемом процессе и моделирует только их. В этом случае при установившемся движении зависимость (13.9) распадается на две

Еи = f(Re) и Еи = f(Fr)    (13.10)

из которых и выбирается определяющая.

Дополнительно отметим, что для сжимаемых сред в число определяющих чисел подобия помимо полученных выше входит и число Маха.

Остается открытым лишь вопрос о кинематическом подобии. Опыт многочисленных исследований показывает, что для его решения не требуется каких-либо специальных мер. Если системы динамически подобны и течение происходит в геометрически подобных каналах, то кинематическое подобие обеспечивается автоматически.

Кратко остановимся еще на двух вопросах, носящих принципиальный характер. Первый из них связан с понятием геометрического подобия «в большом» и «малом». В начале раздела было показано, что геометрическое подобие натуры и модели может быть легко реализовано. Подобие геометрических границ объектов относится к подобию «в большом». Вместе с тем стенки каналов как натуры, так и модели имеют какую-то шероховатость. Очевидно, что моделирование шероховатости практически невозможно, и геометрическое подобие «в малом» недостижимо.

Второй вопрос связан с так называемым «масштабным эффектом». Суть его в том, что моделирование, основанное на классических принципах теории подобия, не обеспечивает масштабный переход. Это означает, что эффективность различного рода промышленных технологических аппаратов оказывается ниже той, которая должна была бы быть по результатам, полученным пересчетом с модельных испытаний. Более того, она ухудшается по мере увеличения размеров аппаратов. Это вынуждает исследователей отказываться от испытаний на моделях и переходу к испытаниям на объектах, построенных в натуральную величину, что резко повышает стоимость эксперимента, а при создании особо крупных аппаратов такой подход вообще невозможно реализовать. Исследования, выполненные в последние годы, показали, что в основе масштабного эффекта лежат чисто гидродинамические явления: неравномерность распределения потоков по сечению аппарата, увеличение масштаба турбулентности и т.п., что позволяет найти способы устранения этого эффекта. Достаточно полное изложение теории можно найти в книге под ред. А.М. Розена «Масштабный переход в химической технологии: разработка

промышленных аппаратов методом гидродинамического моделирования»: Химия, 1980. - 320 с.

13.2. Понятие об автомодельности.

Автомодельность - кардинальное понятие теории подобия, принципиальное содержание которого сводится к так называемому вырождению чисел подобия. Формальным признаком её служит выпадение чисел подобия как аргументов, входящих в функциональную зависимость.

Обстоятельное рассмотрение этого вопроса можно найти в книге А.А.Гухмана «Применение теории подобия к исследованию процессов тепломассообмена. Процессы переноса в движущейся среде». - М.: Высшая шкала,1967. - 302 с.

Мы же ограничимся лишь кратким рассмотрением содержания этого понятия без уяснения которого невозможна грамотная постановка эксперимента.

Для простоты будем считать, что в интересующем исследователя процессе определяющими является силы вязкого трения т.е. зависимость (13.8) имеет вид Eu - f (Re). График этой

зависимости устанавливается экспериментально, и часто имеет вид, показанный на рис. 13.1.

Как следует из рисунка, при увеличении числа Рейнольдса в опытах зависимость Eu - f (Re)

Рис. 13.1    Исчезновение (вырож

дение) числа Рейнольдса означает отсутствие предпосылок для подобия . Очевидно, механизм процесса таков, что не надо никаких условий для подобия и все процессы такого типа автоматически подобны между собой. Этот случай и называется автомодельностью. На рис.13.1 автомодельная область обозначена римской цифрой II.


ослабевает и при некотором конкретном для каждого случая значении числа Re, называемого граничным (Regp) происходит

Eu


«вырождение», т.е. число Эйлера перестает зависеть от Re.

Re—-


В общем случае под автомодельной понимают область, в которой неопределяющее число подобия перестает зависеть от определяющего (либо определяющих).

Проведение опытов в этой области существенно упрощается. Действительно, если в области I экспериментатор должен заботиться о том, чтобы ReM = ReH, что далеко не всегда возможно, то в автомодельной области достаточно, чтобы ReM было больше Re^. Нужно лишь помнить, что какого-то универсального значения Re^ не

существует, оно всегда зависит от природы изучаемого объекта, в частности, от его формы. Поэтому, как правило, задачей первого этапа экспериментального исследования является нахождение граничного значения определяющего числа подобия.

Таким образом, приведенные сведения показывают, что если в результате анализа изучаемого явления удается составить его математическую модель, то принципиально задача постановки эксперимента может считаться разрешенной. К сожалению, возможность аналитического описания является скорее исключением, чем правилом. Поэтому целью следующего раздела является ознакомление со стратегией исследователя при возникновении такой ситуации.

13.3. Анализ размерностей.

Следует подчеркнуть, что конечная цель в рассматриваемом случае остается прежней: нахождение чисел подобия, по которым следует вести моделирование, но решается она при существенно меньшем объеме информации о характере процесса.

Для уяснения дальнейшего кратко рассмотрим некоторые основополагающие понятия. Обстоятельное изложение можно найти в книге А.Н.Лебедева «Моделирование в научно-технических исследованиях». - М.: Радио и связь. 1989. -224 с.

Любой материальный объект обладает рядом свойств, которые допускают количественное выражение. При этом каждое из свойств характеризуется размером определенной физической величины. Единицы некоторых физических величин можно выбирать произвольно, и с их помощью представлять единицы всех остальных. Физические единицы, выбираемые произвольно, называют основными. В международной системе (применительно к механике) это - килограмм, метр и секунда. Остальные величины, выраженные через эти три, называют производными.

Основная единица может обозначаться либо символом соответствующей величины, либо специальным символом.

Например, единицы длины - L, единицы массы - M, единица времени

- T. Либо, единица длины - метр (м), единица массы - килограмм (кг), единица времени - секунда (с).

Под размерностью понимают символическое выражение (иногда его называют формулой) в виде степенного одночлена, связывающее производную величину с основными. Общий вид этой закономерности имеет вид

dim X = LxMyTz    (13.11)

где x, y, z- показатели размерности.

Например,размерность скорости

dim v = L • T_1 Для безразмерной величины все показатели x = y = z = 0, и, следовательно, dim X = 1.

Два следующих утверждения достаточно ясны и не нуждаются в каких-либо специальных доказательствах.

Отношение размеров двух объектов является величиной постоянной вне зависимости от того, в каких единицах они выражаются. Так, например, если отношение площади, занимаемой окнами, к площади стен составляет 0,2, то этот результат останется неизменным, если сами площади выражать в мм2, м2 или км2.

Второе положение можно сформулировать следующим образом. Любое правильное физическое соотношение должно быть размерностно однородным. Это означает, что все члены, входящие как в правую, так и в левую его части должны иметь одинаковую размерность. Это простое правило четко реализуется в житейском обиходе. Все осознают, что метры можно складывать только с метрами и никак не с килограммами или с секундами. Нужно четко представлять, что правило остается справедливым и при рассмотрении даже самых сложных уравнений.

Метод анализа размерностей базируется на так называемой п-теореме (читается: пи-теорема). п-теорема устанавливает связь между функцией, выраженной через размерные параметры, и функцией в безразмерной форме. Более полно теорема может сформулирована так:

Любая функциональная зависимость между размерными величинами может быть представлена в виде зависимости между N безразмерными комплексами (числами п), составленными из этих величин. Число этих комплексов N = m - n, где n - число основных единиц. Как уже отмечалось выше, в гидромеханике n = 3 (кг, м, с).

Пусть, например, величина А является функцией пяти размерных величин (m = 5), т.е.

Из п-теоремы следует, что эта зависимость может быть преобразована в зависимость, содержащую два числа п

(N - m - n - - 5 - 3 - 2)

A - f1(n1,n2)    (13.13)

где п1 и п2 - безразмерные комплексы, составленные из размерных величин.

Эту теорему иногда приписывают Бэкингему и называют п-теоремой Бэкингема. В действительности в её разработку внесли вклад многие крупные ученые, в том числе Фурье, Рябушинский, Рэлей.

Доказательство теоремы выходит за рамки курса. При необходимости оно может быть найдено в книге Л.И.Седова «Методы подобия и размерностей в механике» - М.: Наука, 1972. - 440 с. Подробное обоснование метода приводится и в книге В. А. Веникова и Г.В.Веникова «Теория подобия и моделирования» - М.: Высшая школа, 1984. -439 с. Особенностью этой книги является то, что помимо вопросов, связанных с подобием, в нее включены сведения о методике постановки эксперимента и обработки его результатов.

Использование анализа размерностей для решения конкретных практических задач связано с необходимостью составления функциональной зависимости вида (13.12), которая на следующем этапе обрабатывается специальными приемами, приводящими в конечном итоге к получению чисел п (чисел подобия).

Основным, носящим творческий характер, является первый этап, так как получаемые результаты зависят от того, насколько правильно и полно представление исследователя о физической природе процесса. Другими словами, насколько функциональная зависимость (13.12) правильно и полно учитывает все параметры, влияющие на изучаемый процесс. Любая ошибка здесь неизбежно приводит к ошибочным выводам. В истории науки известна так называемая «ошибка Рэлея». Суть ее в том, что изучая задачу о теплообмене при турбулентном течении, Рэлей не учел влияние вязкости потока, т.е. не включил её в зависимость (13.12). В результате в конечные соотношения, полученные им, не вошло число подобия Рейнольдса, играющее исключительно важную роль в теплообмене.

Для уяснения сущности метода рассмотрим пример, иллюстрирующий как общий подход к задаче, так и способ получения чисел подобия.

Необходимо установить вид зависимости, позволяющий определить потери давления либо напора при турбулентном течении в круглых трубах.

Напомним, что эта задача уже рассматривалась в разделе

12.6. Поэтому представляет несомненный интерес установить, как она может быть разрешена с помощью анализа размерностей и дает ли это решение какую-то новую информацию.

Ясно, что падение давления вдоль трубы, обусловленное затратами энергии на преодоление сил вязкого трения обратно пропорционально её длине, поэтому с целью сокращения числа переменных целесообразно рассматривать не Ар, а Ар / I, т.е. потери давления на единицу длины трубы. Напомним, что отношение Ah / I, где Ah - потери напора, носит название гидравлического уклона.

Из представлений о физической сущности процесса можно предположить что возникающие потери должны зависеть: от средней скорости течения рабочей среды (v); от размера трубопровода, определяемого его диаметром (d); от физических свойств транспортируемой среды, характеризуемых её плотностью (p) и вязкостью (р); и, наконец, разумно считать, что потери должны быть как-то связаны с состоянием внутренней поверхностью трубы, т.е. с шероховатостью (к) ее стенок. Таким образом, зависимость (13.12) в рассматриваемом случае имеет вид

либо

V    I    J


(13.14)

На этом и заканчивается первый и, нужно подчеркнуть, наиболее ответственный этап анализа размерностей.

В соответствии с п-теоремой, число влияющих параметров, входящих в зависимость, m = 6. Следовательно, число безразмерных комплексов N = m - n = 6 - 3 = 3, т.е. после соответствующей обработки (13.14) должна принять вид

) = 0


f3 (


^, п3


(13.15)


Существует несколько способов нахождения чисел п. Мы воспользуемся методом, предложенным Рэлеем.

Основным достоинством его является то, что он представляет собой своеобразный алгоритм, приводящий к решению задачи.

Из параметров, входящих в (13.15) необходимо выбрать три любых, но так, чтобы в них входили основные единицы, т.е. метр, килограмм и секунда. Пусть ими будут v, d, p. Легко убедиться, что они удовлетворяют поставленному требованию.

Образуются числа п в виде степенных одночленов из выбранных параметров, умноженных на один из оставшихся в (13.14)

п = vx1 dy'pz' Ар;    (13.16)

п2 = vx 2 dy 2 Pz2 ju;    (13.17)

п3 = vx 3 dy 3 PZ3 k;    (13.18)

Теперь задача сводится к нахождению всех показателей степеней. При этом они должны быть подобраны так, чтобы числа п были безразмерны.

Для решения этой задачи определим прежде всего размерности всех параметров:

dim v = L • T_1; dim d = L ; dim p= M • L~z

кг • м • с кг    ..    -i-,--i

Вязкость [u] — Па • с — —22——-, т.е. dim u = ML T

с • м мс

Ар . l

кг• м    кг    ..    Ар    _-2, -2

—— —2—2-—— —22, и di m — — MT L .

с • м • м    м • с    l


Параметр


И, наконец, dim k = L.

Таким образом, размерности чисел п будут

dim п = (LT)x 1 Lyi (ML -3)zi MT-2L -2

либо

dim п = Lx1 +yi-3zi-2T-x 1 -2MZl +1 Аналогично два других

dim п2 = Lx2 +y2-3Z2-1T-x2-1MZ2+1

dim п3 = Lx3 +Уз-3z3 +1!-x3MZ3 В начале раздела 13.3 уже отмечалось, что для любой безразмерной величины показатели размерности x = y = z = 0. Поэтому, например, для числа п1 можем записать Lxi + yi-3zi-2t—xi-2M zi 1 = L ^T ^M0 Приравнивая показатели степеней, получаем три уравнения с тремя неизвестными

x 1 + y1 - 3z-, - 2 = 0 - x 1 - 2 = 0 z1 + 1 = 0 Откуда находим x 1 = -2; z1 = -1; y1 = 1.

Подставляя эти значения в ( 13.6), получаем

Ap d

27 pv l

Действуя аналогично, легко показать, что

vd p vd

k


П2 -    H ~

- Re и п3 - -j.

v    3    d


- 0


(13.20)


3


Таким образом, зависимость (13.15) принимает вид

( Ap d    кл

—^ —, Re, -\pv2 Г    ' dj


есть неопределяющее число подобия (число Эйлера),

Ap

pv 2


Так как


то (13.20) можно записать как функциональную зависимость

^ % - f {Re к]- о


pv l


либо


k] l    2

Ap - f Re, - -pv [    dJ d


dj


V


(13.21)


Следует иметь в виду, что анализ размерностей не дает и

принципиально не может дать каких-то числовых значений в

получаемых с его помощью соотношениях. Поэтому он должен

завершаться анализом результатов и при необходимости их

корректировкой, исходя из общих физических представлений.

Рассмотрим с этих позиций выражение (13.21). В правую его часть

входит квадрат скорости, но эта запись не выражает ничего, кроме

того, что скорость возводится в квадрат. Однако, если поделить эту

v2

величину на два, т.е. т, то как известно из гидромеханики, она

^    W    I    W    W    W

приобретает важный физический смысл: удельной кинетической

v 2

энергии, а — динамическое давление, обусловленное средней скоростью. С учетом этого (13.21) целесообразно записать в виде

( к] l v

Ap - f Re, k - p—

[ dJ dH 2

(    k]

Если теперь, как в (12.26), обозначить f Re, — буквой Л, то

[    dJ

приходим к формуле Дарси

либо

I v2

Ah = Я----(13.24)

d 2g

где Я - гидравлический коэффициент трения, который, как следует из (13.22), является функцией числа Рейнольдса и относительной шероховатости (k/d). Вид этой зависимости может быть найден только экспериментальным путем.

ОГЛАВЛЕНИЕ

1.    МАТЕМАТИЧЕСКИЙ АППАРАТ, ИСПОЛЬЗУЕМЫЙ В

Учебное издание

Ким Федорович Иванов Сергей Владимирович Сурков

Механика жидкости и газа (конспект лекций) Часть 1

Редактор    Т.И.Лучнева

Корректор    Л.А.Гречанова

Подписано к печати    .    Формат    60x84/16.

Бумага газетная. Печать офсетная    усл. печ. л.

уч.-изд. л. Тираж    экз. Заказ №

Одесский государственный политехнический университет 270044, Одесса, пр. Шевченко, 1


1.    Абдурахманов И.М., Алишаее М.Г. Плоская стационарная фильтрация в пласте, разделенном прямолинейной трещиной // Изв. АН СССР. Механика жидкости и газа-1973 - + 4.- С. 173-177.

2.    Азиз X., Сеттари Э. Математическое моделирование пластовых систем. - М.: Недра, 1982. - 408 с.

3.    Антипов Д.М., Ибрагимов А.И., Панфилов М.Б. Модель сопряженного течения жидкости в пласте и внутри горизонтальной скважины // Изв. РАН. Механика жидкости и газа. - 1996. - + 5. - С.

112-117.

4.    Афанасьева А.В., Горбунов А.Т., Шустеф И.Н. Заводнение нефтяных месторождений при высоких давлениях нагнетания.- М.: Недра, 1975. - 216 с.

5.    Бан А., Богомолова А.Ф., Максимов В.А., Николаевский В.Н., Оганджанянц В.Г., Рыжик В.М. Влияние свойств горных пород на движение в них жидкости.- М.: Гостоптехиздат, 1962.- 276 с.

6.    Баренблатт Г.И. О некоторых задачах теории упругости, возникающих при исследовании механизма гидравлического разрыва пласта // Прикл. матем. и механика.- 1956.- Т. 20.- + 4.- С. 475-486.

7.    Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в пористых пластах.- М.: Недра, 1984.- 208 с.

8.    Басниев К.С., Дмитриев Н.М. Обобщенный закон Дарси для анизотропных пористых сред // Изв. вузов. Нефть и газ.- 1986.- + 5.- С. 54-59.

9.    Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидродинамика.- М.: Недра, 1993.- 416 с.

10.    Гайдуков В.Н. Анализ результатов гидроразрыва пласта на месторождениях ОАО “Варьеганнефтегаз” // Сб. докладов научно-технич. совещания по обмену опытом и выработке основных направлений деятельности в области повышения эффективности разработки месторождений ОАО “Варьеганнефтегаз”. - М.: НК “СИДАНКО”.- 1998. -С. 146-160.

11.    Глова В.Н., Латышев В.Н. Результаты гидроразрыва пласта на месторождениях ОАО “Пурнефтегаз”// Нефтяное хозяйство. -1996.- + 1.- С. 15-18.

12.    Голубева О.В. Обобщение теоремы об окружности на фильтрационные течения (К вопросу о течениях в кусочно-неоднородных грунтах) II Изв. АН СССР. Механика жидкости и газа.-1966.- № 1.- С.

113-116.

13.    Голубева О.В., Шпилевой А.Я. О плоской фильтрации в средах с прерывно изменяющейся проницаемостью вдоль кривых второго порядка II Изв. АН СССР. Механика жидкости и газа. - 1967.- № 2.- С. 174-179.

14. Голф-Рахт Т.Д. Основы нефтепромысловой геологии и разработки трещиноватых коллекторов.- М.: Недра, 1986.- 608 с.

15.    Гордеев Ю.Н. Автомодельные задачи о развитии протяженной трещины гидроразрыва в проницаемой среде II Изв. РАН. Механика жидкости и газа.- 1992.- № 2.- С. 91-101.

16. Гордеев Ю.Н., Ентов В.М. О распределении давления в окрестности растущей трещины II Прикл. матем. и механика.- 1997. - Т. 61. - Вып.

6.- С. 1060-1064.

17.    Гордеев Ю.Н., Зазовский А.Ф. Автомодельное решение задачи о глубокопроникающем гидравлическом разрыве пласта II Изв. АН СССР. Механика твердого тела.- 1991.- № 5.- С. 119-131.

18.    Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений.- М.: Физматгиз, 1962. - 1100 с.

19.    Григорян С.С. Об осреднении физических полей II Доклады АН СССР. - 1980. - Т. 254. - № 4. - С. 846-850.

20.    Гуревич А.В., Крылов А.Л., Топор Д.Н. Решение плоских задач гидродинамики пористых сред вблизи разрывных нарушений методом комплексного потенциала II Докл. АН СССР.- 1988. - Т. 298.- № 4.- С. 846-850.

21.    Гусев С.В., Бриллиант Л.С., Янин А.Н. Результаты широкомасштабного применения ГРП на месторождениях Западной Сибири II Материалы совещания “Разработка нефтяных и нефтегазовых месторождений” (г. Альметьевск, 1995 г.).- М.: ВНИИОЭНГ.- 1996.- С. 291-303.

22.    Гусев С.В., Коваль Я.Г., Кольчугин И.С. Анализ эффективности гидроразрыва пластов на месторождениях ПО “Юганскнефтегаз” II Нефтяное хозяйство.- 1991.- № 7.- С. 15-18.

23.    Данилов В.Л., Ка— Р.М. Гидродинамические расчеты взаимного вытеснения жидкостей в пористой среде. - М.: Недра, 1980.- 264 с.

24.    Данилов В.Л. Методы установления в прикладных обратных задачах потенциала. - М.: Наука, 1996. - 248 с.

25.    Джеффрис Г., Свирлс Б. Методы математической физики.- М.: Мир, 1969. - Вып. 1. - 424 с.

26.    Джеффрис Г., Свирлс Б. Методы математической физики.- М.: Мир, 1970. - Вып. 2. - 352 с.

27.    Дмитриев Н.М., Максимов В.М. Определяющие уравнения двухфазной фильтрации в анизотропных пористых средах II Изв. РАН.

Механика жидкости и газа. - 1998.- № 2. -    С. 87-94.

28.    Ентов В.М., Мурзенко В.В. Стационарная фильтрация однородной жидкости в элементе разработки нефтяного пласта с трещиной гидроразрыва // Изв. РАН. Механика жидкости и газа. - 1994.- № 1.- С. 104-112.

29.    Жданов С.А., Константинов С.В. Проектирование и применение гидроразрыва пласта в системе скважин // Нефтяное хозяйство.-

1995.- № 9.- С. 24-25.

30.    Желтое Ю.П. Деформации горных пород. - М.: Недра, 1966. - 198 с.

31.    Желтое Ю.П. Механика нефтегазоносного пласта.- М.: Недра, 1975.

-    207 с.

32.    Желтое Ю.П., Христианоеич С.А. О гидравлическом разрыве

нефтеносного пласта // Изв. АН СССР. ОТН. - 1955. -    №    5.    -    С.    3

41.

33.    Зазоеский А.Ф., Лемеха А.В., Федоренко Р.П. Нелинейные эффекты при циркуляции жидкости в трещине гидроразрыва // Ин-т прикл. матем. Препринт.- М., 1987. - № 210. - 30 с.

34.    Зазоеский А.Ф., Лемеха А.В., Федоренко Р.П. Численный анализ упругогидродинамической задачи о циркуляции жидкости в трещине гидроразрыва // Ин-т прикл. матем. Препринт.- М., 1987. - № 1.

-    29 с.

35.    Зазоеский А.Ф., Тодуа Г.Т. О стационарном притоке жидкости к скважине с вертикальной трещиной гидроразрыва большой протяженности // Изв. АН СССР. Механика жидкости и газа. - 1990. - № 4. - С. 107-116.

36.    Ильин А.Ф. Астраханское газоконденсатное месторождение: состояние и проблемы освоения // Газовая промышленность. - 1991. - № 11. -С. 10-11.

37.    Кадет В.В., Селякое В.И. Фильтрация флюида в среде, содержащей эллиптическую трещину гидроразрыва // Изв. вузов. Нефть и газ.

-    1988.- № 5. - С. 54-60.

38.    Канееская Р.Д. Асимптотический анализ влияния капиллярных и гравитационных сил на двумерный фильтрационный перенос двухфазных систем // Изв. АН СССР. Механика жидкости и газа. - 1988. -№ 4. - С. 88-95.

39.    Канееская Р.Д. Асимптотический анализ моделей двухфазной фильтрации и определение областей их применимости // Сб. науч. тр. ВНИИ. - Вып. 102. - М., 1988. - С. 108-122.

40.    Канееская Р.Д. Влияние неполноты вытеснения нефти водой в отдельных пропластках на вид модифицированных фазовых проницаемостей слоистого пласта // Сб. науч. тр. ВНИИ. - Вып. 103. - М., 1988.

-    С. 110-121.

41.    Каневская Р.Д. Зарубежный и отечественный опыт применения гидроразрыва пласта. - М.: ВНИИОЭНГ, 1998. - 40 с.

42.    Каневская Р.Д. К вопросу о расчете процесса вытеснения нефти водой из системы изолированных пропластков // Сб. науч. тр. ВНИИ.

- Вып. 95.- М., 1986. - С. 31-42.

43.    Каневская Р.Д., Кац Р.М. Аналитические решения задач о притоке жидкости к скважине с вертикальной трещиной гидроразрыва и их использование в численных моделях фильтрации // Изв. РАН. Механика жидкости и газа.- 1996. - № 6. - С. 69-80.

44.    Каневская Р.Д., Кац Р.М. Оценка эффективности гидроразрыва добывающих и нагнетательных скважин при различных системах заводнения пласта // Нефтяное хозяйство. - 1998. - № 6. - С. 34-37.

45.    Каневская Р.Д. Методы введения модифицированных фазовых проницаемостей // Сб. науч. тр. ВНИИ. - Вып. 94. - М., 1986. - С. 45-52.

46.    Каневская Р.Д. О влиянии направления трещин гидроразрыва на динамику обводнения продукции скважин // Нефтяное хозяйство. -1999. - № 2. - С. 26-29.

47.    Каневская Р.Д. О притоке жидкости к скважине с вертикальной трещиной гидроразрыва в кусочно-однородном анизотропном пласте // Изв. РАН. Механика жидкости и газа.- 1999. - № 2. - С. 64-71.

48.    Каневская Р.Д. О притоке жидкости к скважине с трещиной гидроразрыва в слоистом пласте // Нефтепромысловое дело. - 1999. - № 1. -С. 6-8.

49.    Каневская Р.Д. Оценка влияния гидроразрыва на дебит скважин в неоднородных коллекторах // Геология, геофизика и разработка нефтяных месторождений.- 1999. - № 2. - С. 26-28.

50.    Кац Р.М., Андриасов А.Р. Математическая модель трехфазной фильтрации в трещиновато-пористой среде // Сб. науч. тр. ВНИИ. -М., 1986. - Вып. 95. - С. 61-66.

51.    Кац Р.М., Каневская Р.Д. Проектирование разработки нефтяных месторождений с применением гидроразрыва пластов (ГРП) на основе современных компьютерных технологий // Материалы научнотехнической конференции “Повышение уровня добычи нефти на месторождениях ОАО ННГ в 1997-2005 гг.” (г. Ноябрьск, 1997 г.). - М.: ВНИИОЭНГ, 1999. - С. 248-259.

52.    Кокорин А.А., Заболотнов А.Р. Особенности разработки юрских залежей нефти Нижневартовского района с применением гидроразрыва пласта // Нефтяное хозяйство. - 1997. - № 10. - С. 54-57.

53.    Константинов С.В., Гусев В.И. Техника и технология проведения гидравлического разрыва пластов за рубежом. - М.: ВНИИОЭНГ. Обзорная информация. - Сер. Нефтепромысловое дело. - 1985. - 61 с.

54.    Константинов С.В., Лесик Н.П., Гусев В.И., Борисов Ю.П.

Глубокопроникающий гидравлический разрыв пласта - метод интенсификации разработки низкопроницаемых коллекторов // Нефтяное хозяйство. - 1987. - № 5. - С. 22-25.

55.    Константиное С.В., Матееее Ю.М. Результаты опытных работ по гидроразрыву канадской фирмы “Фракмастер” на месторождениях ПО “Юганскнефтегаз” // Нефтяное хозяйство. - 1989. - № 6. - С. 20-26.

56.    Копаее А.В. Метод функциональных уравнений в задачах фильтрации в слоистой среде // Изв. РАН. Механика жидкости и газа. - 1997. -№ 5. - С. 81-89.

57.    Криеоносое И.В., Чарный И.А. Расчет дебитов скважин с трещиноватой призабойной зоной пласта // Нефтяное хозяйство. - 1955. - №

7. - С. 40-47.

58.    Кристиан М., Сокол С., Константинеску А. Увеличение продуктивности и приемистости скважин. - М.: Недра, 1985.- 185 с.

59.    Кричлоу Г.Б. Современная разработка нефтяных месторождений -проблемы моделирования. - М.: Недра, 1979. - 303 с.

60.    Курамшин Р.М., Иеаное С.В., Кузьмичее Н.Д. Эффективность проведения гидроразрыва пласта на месторождениях Ноябрьского района // Нефтяное хозяйство. - 1997. - № 12. - С. 58-60.

61.    Кураное И.Ф., Шехтман Ю.М. Определение дебита скважины при наличии горизонтальной трещины с заполнителем // Нефтяное хозяйство. -1961.- № 9. - С. 37-39.

62.    Курант Р. Уравнения с частными производными. - М.: Мир, 1964. -830 с.

63.    Курбаное А.К. О некоторых обобщениях уравнений фильтрации двухфазной жидкости // Науч.-техн. сб. ВНИИ. - Вып. 15. - М., 1961.

-    С. 32-38.

64.    Курбаное А.К., Атаное Г.А. К вопросу о вытеснении нефти водой из неоднородного пласта // Нефть и газ Тюмени. - 1974. - Вып. 13. -С. 36-38.

65.    Лисоеский Н.Н., Филиппое В.П. Состояние разработки нефтяных месторождений России и задачи по дальнейшему ее усовершенствованию // Материалы совещания “Разработка нефтяных и нефтегазовых месторождений” (г. Альметьевск, 1995 г.). - М.: ВНИИОЭНГ. - 1996.

-    С. 3-18.

66.    Логиное Б.Г., Блажееич В.А. Гидравлический разрыв пластов. - М.: Недра, 1966. - 148 с.

67.    Максимое М.М., Рыбицкая Л.П. Математическое моделирование процессов разработки нефтяных месторождений. - М.: Недра, 1976. -264 с.

68.    Максимоеич Г.К. Гидравлический разрыв нефтяных пластов. - М.: Гостоптехиздат, 1957. - 98 с.

69.    Малышев А.Г., Малышев Г.А., Журба В.Н., Сальникова Н.Н.

Анализ технологии проведения ГРП на месторождениях ОАО “Сургутнефтегаз” II Нефтяное хозяйство. -1997. - № 9. - С. 46-51.

70.    Марчук Г.И. Методы вычислительной математики. - М.: Наука, 1980.

-    536 с.

71.    Маскет М. Течение однородных жидкостей в пористой среде. - М.: Гостоптехиздат, 1949. - 628 с.

72.    Методическое руководство по проектированию разработки нефтяных месторождений с применением гидроразрыва пластов на основе современных компьютерных технологий. РД 153-39.2-032-98 I Кац Р.М., Каневская Р.Д. и др.- М.: Минтопэнерго РФ, 1998. - 70 с.

73.    Муравьев И.М., Го Шан-пин. Об эффективности проведения массового гидравлического разрыва пласта II Нефтяное хозяйство. -1958.

-    № 4. - С. 39-44.

74.    Мурзенко В.В. Аналитические решения задач стационарного течения жидкости в пластах с трещинами гидроразрыва II Изв. АН СССР. Механика жидкости и газа. - 1994.- № 2.- С. 74-82.

75.    Николаевский В.Н., Басниев К.С., Горбунов А.Т., Зотов Г.А.

Механика насыщенных пористых сред. - М.: Недра, 1970.- 336 с.

76.    Николаевский В.Н. Геомеханика и флюидодинамика. - М.: Недра, 1996. - 448 с.

77.    Николаевский В.Н. Применение гидравлического разрыва пласта на месторождении Умбаки II Нефтяное хозяйство.-1958. - № 4. - С. 50

53.

78. Пивень В.Ф. Двумерная фильтрация в слоях переменной проводимости, моделируемой гармонической функцией координат II Изв. РАН. Механика жидкости и газа.- 1995.- № 3.- С. 102-112.

79.    Пивень В.Ф. О двумерной фильтрации в    слоях с    прерывно изменяющейся проводимостью вдоль кривых    второго    порядка II    Изв.

РАН. Механика жидкости и газа. - 1993. -    № 1. - С.    120-128.

80.    Пискунов Н.С. Разрыв пласта и влияние    разрыва    на процесс    экс

плуатации месторождений II Труды ВНИИ. - Вып.ХУЪ- М.: Гостоптехиздат, 1958. - С. 3-24.

81. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. - М.: Наука, 1977. - 664 с.

82.    Регламент составления проектных технологических документов на разработку нефтяных и газонефтяных месторождений. РД 153-39-007

96. - М.: Минтопэнерго РФ, 1996. - 204 с.

83.    Реутов В.А. Гидравлический разрыв пласта II Итоги науки и техники. Механика деформируемого твердого тела.- М.: ВИНИТИ, 1989. -Т. 20. - С. 84-188.

84.    Реутов В.А. Гидравлический разрыв пласта: условия образования трещин, их практическое определение и использование // Итоги науки и техники. Разработка нефтяных и газовых месторождений.- М.: ВИНИТИ, 1991.- Т. 23. - С. 73-153.

85.    Седов Л.И. Механика сплошной среды. - М.: Наука, 1976.- 536 с.

86.    Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Добыча нефти / Под ред. Ш.К. Гима-тудинова. - М.: Недра, 1983. - 456 с.

87.    Токунов В.И., Рылов Е.А., Поляков Г.А., Кунавин В.В., Булда ков С.В., Бусыгин Г.И. Гидроразрыв пласта на Астраханском ГКМ // Газовая промышленность.- 1998.- № 3.- С. 47-48.

88.    Усачев П.М. Гидравлический разрыв пласта.- М.: Недра.- 1986.

89.    Халимов Э.М., Леви Б.И., Дзюба В.И., Пономарев С.А. Технология повышения нефтеотдачи пластов.- М.: Недра, 1984.- 271 с.

90.    Холодовский С.Е. О решении задач нестационарной фильтрации в средах с протяженной трещиной (завесой) // Изв. РАН. Механика жидкости и газа.- 1995.- № 6.- С. 95-98.

91.    Холодовский С.Е. О фильтрации жидкости в кусочно-однородных средах с трещинами и завесами, расположенными на одной прямой // Изв. РАН. Механика жидкости и газа.- 1996. - № 6. - С. 85-91.

92.    Чарный И.А. Подземная гидрогазодинамика.- М.: Гостоптехиздат, 1963. - 346 с.

93.    Чарный И.А. Подземная гидромеханика.- М.-Л.: Гостехиздат, 1948.

-    196 с.

94.    Чикеш М., Чубрич С. Гидравлический разрыв пласта как средство оптимизации добычи нефти из юрских залежей нефти Западной Сибири // Сб. докладов научно-технич. совещания по обмену опытом и выработке основных направлений деятельности в области повышения эффективности разработки месторождений ОАО “Варьеганнефтегаз”.

-    М.: НК “СИДАНКО”. - 1998. - С. 266-295.

95.    Швидлер М.И. Приток жидкости к скважине с трещиной в призабойной зоне // Изв. АН СССР. ОТН. - 1955.- № 11.- С. 95-100.

96.    Швидлер М.И. Статистическая гидродинамика пористых сред. - М.: Недра, 1985.- 288 с.

97.    Шелепов В.В., Таранин В.В. Анализ применения ГРП на Пов-ховском месторождении // Материалы совещания “Разработка нефтяных и нефтегазовых месторождений” (г. Альметьевск, 1995 г.). - М.: ВНИИОЭНГ, 1996. - С. 303-308.

98.    Шехтман Ю.М. Приток жидкости к горизонтальной осесимметричной трещине с заполнителем // Изв. АН СССР. ОТН. Механика и машиностроение. -1959.- № 5.- С. 53-57.

99.    Шехтман Ю.М. Приток жидкости к горизонтальной трещине при конечной мощности пласта // Изв. АН СССР. ОТН. Механика и маши-

100.    Шехтман Ю.М. Приток жидкости к одиночной вертикальной трещине с заполнителем // Изв. АН СССР. ОТН.- 1957.- № 7. - С. 146

149.

101.    Шпурое И.В., Разуменко В.Е., Горее В.Г., Шарифуллин Ф.А.

Анализ эффективности разработки залежей нефти Самотлорского месторождения с применением гидроразрыва пласта // Нефтяное хозяйство. - 1997.- № 10.- С. 50-53.

102.    Щелкачее В.Н. Разработка нефтеводоносных пластов на упругом режиме. - М.: Гостоптехиздат, 1959. - 468 с.

103.    Щурое В.И., Трубина А.Ф. Решение при помощи метода электролитического моделирования задачи о притоке жидкости к скважине при наличии трещины в пласте // Труды ВНИИ.- Вып. XVI. - М.: Гос-топтехиздат, 1958.- С. 86-105.

104.    Abass H.H., Hedayati S., Meadows D.L. Nonplanar fracture propagation from a horizontal wellbore: experimental study // SPE Prod. & Fac. -1996. - V.11. - N 3. - P. 133-137.

105.    Abou-Sayed I.S., Schueler S.K., Ehrl E., Hendriks W.P. Multiple hydraulic fracture stimulation in a deep, horizontal tight-gas well // J. Petrol. Technol. - 1996.- N 2.- P. 163-168.

106.    Agarwal R.G., Carter R.D., Pollock C.B. Evaluation and performance prediction of low permeability gas wells stimulated by massive hydraulic fracturing // J. Petrol. Technol.- 1979. -V. 31.- N 3. - P. 362-372.

107.    Agarwal R.G., Carter R.D., Pollock C.B. Type curves for evaluation and performance prediction of low permeability gas wells stimulated by massive hydraulic fracturing // J. Petrol. Technol. - 1979. - V. 31. - N 5. -P. 651-654.

108.    Ahmed U., Khelkar S., Schatz J.F. Mini-frac: an aid to formation in-situ stress and permeability measurements // Proc., 24th Symposium on rock mechanics.- 1983.- P. 319-326.

109.    Al-Hashim H., Kissami M, Al-Yousef H.Y. Effect of multiple hydraulic fractures on gas-well performance // J. Petrol. Technol. - 1993. - V. 45.-N 6.- P. 558-563.

110.    Anderson S.A. Exploring reservoirs with horizontal wells: the Maersk experience offshore // Offshore. -1991.- V. 51.- N 2. - P. 23.

111.    Appah D. Application of the theory of diffuse set to optimize hydraulic fracturing // J. of Petrol. Science and Eng. - 1994.- V.11.- P. 335-340.

112.    Baumgarthner W.E. Fracture stimulation of a horizontal well in a deep tight gas reservoir: a case history from offshore the Netherlands // Paper SPE 26765.- 1993.

113.    Beier R.A. Pressure-transient model for a vertically fractured well in a fractal reservoir // SPE Form. Eval.- 1994.- V.9.- N 2.- P.122-128.

114.    Ben Naceur K, Roeglers J.-C. Design of fracturing treatments in multilayered formations // SPE Prod. Eng. - 1990. - N 1.- P. 21-26.

115. Bennett C.O., Raghavan R, Reynolds A.C. Analysis of finite-conductivity fractures intercepting multilayer reservoirs // SPE Form. Eval.

- 1986. - N 3. - P. 259-274.

116. Bennett C.O., Reynolds A.C, Raghavan R, Elbel J.E. Performance of finite-conductivity, vertically fractured wells in single-layer reservoirs // SPE Form. Eval.- 1986.- N 4.- P. 399-412.

117.    Bennett C.O, Rosato N.D., Reynolds A.C., Raghavan R. Influence of fracture heterogeneity and wing lenght on the response of vertically fractured wells // Soc. Petrol. Eng. Journal.- 1983.- V. 23.- N 2.- P. 219230.

118.    Berumen S., Tiab D. Interpretation of stress damage on fracture conductivity // J. Petrol. Science and Eng. - 1997. - V. 17.- P. 71-85.

119.    Blanco E.R. Hydraulic fracturing requires extensive disciplinary interaction // Oil and Gas J. -1990. - N 12. - P. 112-118.

120.    Bostic J.N., Agarwal R.G., Carter R.D. Combined analysis of post-fracturing performance and pressure buildup data for evaluating an MHF gas well // J. Petrol. Technol. - 1980. - V. 32. - N 10.- P. 1711-1719.

121.    Camacho R.G., Raghavan R, Reynolds A.C. Response of wells producing layered reservoirs: unequal fracture lenght // SPE Form. Eval. -1987.- V. 2.- N 1. - P. 9-28.

122.    Carter R.D. Type curves for finite radial and linear gas-flow systems: constant terminal pressure case // Paper SPE 12917.- 1984.

123.    Chen C.-C, Raghavan R. Modeling a fractured well in a composite reservoir // SPE Form. Eval.- 1995.- V. 10.- N 4.- P. 241-246.

124.    Chu W, Shank G.D. A new model for a fractured well in a radial composite reservoir // SPE Form. Eval.- 1993.- V. 8.- N 3.- P. 225- 233.

125.    Cinco-Ley H. Evaluation of hydraulic fracturing by transient pressure analysis methods // Paper SPE 10043.- 1982.

126.    Cinco-Ley H, Samaniego V.F. Transient pressure analysis for fractured wells // J. Petrol. Technol. - 1981.- V. 33.- N 9. - P. 1749-1766.

127.    Cinco-Ley H, Samaniego V.F. Transient pressure analysis: finite conductivity fracture case versus damaged fracture case // Paper SPE 10179.1981.

128. Cinco-Ley H, Samaniego V.F., Dominguez N. Transient pressure behavior for a well with a finite-conductivity vertical fracture // Soc. Petrol. Eng. Journal.- 1978.- V. 18.- N 4.- P. 253-264.

129.    Cipolla C.L., Meehan D.N., Stevens P.L. Hydraulic fracturing performance in the Moxa Arch Frontier formation // SPE Prod. and Fac.-1996.- V. 11.- N 4.- P. 216-222.

130.    Clark J.B. Hydraulic process for increasing productivity of wells // Trans.

AIME.- 1949.- V. 186.- P. 1-8.

131.    Clark K.K. Transient pressure testing of fractured water injection wells // J. Petrol. Technol.- 1968.- V. 20.- N 6. - P. 639-643.

132.    Cleary M.P. Comprehensive design formulae for hydraulic fracturing // Paper SPE 9259.- 1980.

133.    Cooke C.E. Conductivity of fracture proppants in multiple layers // J. Petrol. Technol. - 1973.- V. 25.- N 9. - P. 1101- 1107.

134.    Daneshy A.A. On the design of vertical hydraulic fractures // J. Petrol. Technol. - 1973.- N 1.- P. 83-97.

135. Dusterhoft R.G., Chapman B.J. Fracturing high-permeability reservoirs increases productivity // Oil and Gas J. - 1994. - N 20. - P. 40-44.

136. Dyes A.B., Kemp C.E., Caudle B.H. Effect of fractures on sweep-out pattern // Trans. AIME.- 1958.- V. 213.- P. 245-249.

137. Economides M.J., Nolte K.G. Reservoir Stimulation.-Prentice Hall, Eglewood Cliffs, New Jersey 07632. - 1989.- 430 pp.

138.    Ehlig-Economides C., Ramey H.J.Jr. Transient rate decline analysis for wells produced at constant pressure // Soc. Petrol. Eng. Journal. -1981.- V. 21.- N 1.- P. 5-20.

139.    Ehlig-Economides C.A., Economides M.G. Pressure transient analysis in an elongated linear flow system // Soc. Petrol. Eng. Journal. - 1985.-V. 25.- N 6.- P. 839-847.

140.    Elbel J.L. Considerations for optimum fracture geometry design // SPE Prod. Eng.- 1988.- V. 3.- N 3.- P. 323-327.

141.    Elbel J.L. Designing fractures for efficient reserve recovery // Paper SPE 15231.- 1986.

142.    Elkins L.F., Skov A.M. Determination of fracture orientation from pressure interference // Trans. AIME. - 1960. - V. 219. - P. 301-304.

143.    Fairchild N.R.Jr, Paugh L Advanced stimulation technologies cut costs and increase production in Medina stimulation // Paper SPE 37362.1996.

144.    Fast C.R., Holmar G.B., Covlin R.J. A study of the application of MHF to the tight muddy “J” formation, Wattenberg field, Adams and Weld Counties, Colorado // Paper SPE 5624.- 1975.

145. Fetkovich M.J. Decline curve analysis using type curves //    J. Petrol. Technol. - 1980.- V. 32.- N 6.- P. 1065-1077.

146.    Geertsma J., de Klerk F. A rapid method of predicting width and extent of hydraulically induced fractures // J. Petrol. Technol. - 1969. - N 12.- P. 1571-1581.

147.    Gidley J.L, Holditch S.A, Nierode D.E, Veatch R.W. Recent advances in hydraulic fracturing. - Monograph Series. SPE of AIME. Richardson.- TX, 1989.- V. 12.

148. Gordeev Yu.N., Zazovsky A.F. Self-similar solution for deep-penetrating hydraulic fracture propagation // Transport in porous media. -1992. - V. 7. - N 3. - P. 283-304.

149. Gringarten A.C., Ramey H.J.Jr. The use of source and Green’s functions in solving unsteady flow problems in reservoirs // Soc. Petrol. Eng. Journal.- 1973.- V. 13.- N 5. - P. 285-296.

150.    Gringarten A.C., Ramey H.J.Jr. Unsteady-state pressure distributions created by a well with a single horizontal fracture, partial penetration or restricted entry // Soc. Petrol. Eng. Journal.- 1974.- V. 14.- N 4.- P. 413426.

151. Gringarten A.C., Ramey H.J.Jr., Raghavan R. Applied pressure analysis for fractured wells // J. Petrol. Technol. - 1975.- V. 27.- N 7.-P. 887-892.

152. Gringarten A.C., Ramey H.J.Jr., Raghavan R. Unsteady-state pressure distributions created by a well with a single infinite-conductivity vertical fracture // Soc. Petrol. Eng. Journal. - 1974. -V. 14. - N 4. - P. 347-360.

153. Gruber N.G., Anderson H.A. Carbonated hydrocarbons for improved gas well fracturing results // J. Canad. Petrol. Technol.-1996.-V. 35.- N

8.- P. 15-24.

154. Guppy K.H., Cinco-Ley H., Ramey H.J., Samaniego V.F. Non-Darcy flow in wells with finite-conductivity vertical fractures // Soc. Petrol. Eng. Journal.- 1982.- V. 22.- N 5.- P. 681-698.

155. Guppy K.H., Cinco-Ley H, Ramey H.J. Pressure building analysis of fractured wells producing at high flow rates // J. Petrol. Technol. - 1982.-V. 34. - N 11. - P. 2656-2666.

156. Hannah R.R., Park E.I., Porter D.A., Black J.W. Combination fracturing/gravel-packing completion technique on the Amberjack, Mississippi Canion 109 field // SPE Prod. and Fac.- 1994.- V. 9.- N 4.- P. 262-266.

157.    Hartsock J.H., Warren J.E. The effect of horizontal hydraulic fracturing on well performance // J. Petrol. Technol. - 1961.- V. 13. - N 10.- P. 10511056.

158.    Hickey J.W, Brown W.E., Crittenden S.J. The comparative effectiveness of propping agents in the Red Fork Formation of the Anadarko Basin // Paper SPE 10132.- 1981.

159.    Hill A.D., Zhu D., Wang Y. The effect of wormholing on the fluid-loss coefficient in acid fracturing //SPE Prod. and Fac. - 1995. - V. 10.- N 4.- P. 257-263.

160. Holditch S.A., Jennings J.W, Neuse S.H., Wyman R.E. The optimization of well spacing and fracture length in low-permeability gas reservoirs // Paper SPE 7496.- 1978.

161. Holditch S.A, Laufer W.L. The analysis of fractured gas wells using reservoir simulation // Paper SPE 7473.- 1978.

162.    Holditch S.A., Morse R.A. The effects of non-Darcy flow on the behavior hydraulically fractured gas wells // J. Petrol. Technol.- 1976.- V. 28.-N 10.- P. 1169-1178.

163. Holditch S.A, Robinson B.M, Whitehead W.S. Pre-fracture and post-fracture formation evaluations necessary to characterize the 3D Shape of hydraulic fracture // Paper SPE 14086.- 1986.

164. Howard G.C., Fast C.R. Hydraulic fracturing. - SPE Monograph Series.- 1970.- V. 2.- 203 pp.

165.    Howwell J.C., Thomas B.D. Evaluation of injection well stimulation as applied to a large micellar-polymer project // Paper SPE 7180.- 1978.

166.    Hubbert M.K., Willis D.G. Mechanics of hydraulic fracturing // Trans. AIME. - 1957.- V. 210.- P. 153-168.

167. Kanevskaya R.D., Kats R.M. New analytical solutions of fluid flow problems in hydraulic fracturing // Proceedings of the 1998 International Gas Research Conference. San Diego, CA, USA, Nov. 8-11.- 1998.- V. 1.- P. 75-86.

168.    Kozik H.G., Holditch S.A. A case history for massive hydraulic fracturing the Cotton Valley Lime matrix, Fallon and Personville Fields // J. Petrol. Technol.- 1981.- N 2.- P. 229-244.

169.    Lacy L.L. Comparison of hydraulic-fracture orientation techniques // SPE Form. Eval.- 1987.- V. 2.- N 1.- P. 66-76.

170.    Lee W.J., Holditch S.A. Fracture evaluation with pressure transient testing in low-permeability gas reservoir // J. Petrol. Technol. - 1981.- V. 33.-N 9.- P. 1776- 1792.

171.    Mader D. Hydraulic proppant fracturing and gravel packing. Developments in petroleum science. - Elsevier Science Publishers, 1989.- V. 26.1240 pp.

172. Malekzadeh D., Khan F.U, Day J.J. Analysis of pressure behaviour of hydraulically fractured vertical wells by the effective hydraulic fracture lenght concept // J. Canad. Petrol. Technol.- 1996.- V. 35.- N 3.- P. 36

43.

173. Martins J.P, Leung K.H, Jackson M.R, Stewart D.R, Carr A.H.

Tip screen out fracturing applied to the Ravensprun South gas field development // SPE Prod. Eng.- 1992.- V. 7.- N 3.- P. 252-258.

174.    Martins J.R., Murray L.R., Clifford P.J., Maclelland W.G., Hanna M.F, Sharp Jr. J.W. Produced-water reinjection and fracturing in Prad-hoe Bay // SPE Res. Eng.- 1995.-V. 10.- N 3.- P. 176-182.

175. Mayerhofer M.J, Ehlig-Economides C.A, Economides M.G. Pressure-transient analysis of fracture-calibration tests // J. Petrol. Technol.

- 1995.- V. 47.- N 3.- P. 229-234.

176.    McDaniel R.R., Willingham J.R. The effect of various proppants and

proppant mixtures on fracture permeability // Paper SPE 7573.- 1978.

177. McGuire W.J, Sikora V.J. The effect of vertical fractures on well productivity // Trans. AIME.- 1960.- V. 219.- P. 401-403.

178.    Meehan D.N. Optimization of fracture lenght and well spacing in heterogeneous reservoirs // SPE Prod. and Fac. - 1995. - V. 10. - N 2.- P. 82-88.

179.    Meese C.A., Mullen M.E, Barree R.D. Offshore hydraulic fracturing technique // J. Petrol. Technol.- 1994.- V. 46.- N 3. - P. 226-229.

180.    Meng H-Z, Brown K.E. Coupling of production forecasting, fracture geometry requirements and treatment scheduling in the optimum hydraulic fracture design // Paper SPE 16435.- 1987.

181. Meng H-Z., Proano E.A, Buhidma I.M, Mach J.M. Production systems analysis of vertically fractured wells // Paper SPE/DOE 10842. -1 982.

182.    Millheim K.K, Cichowicz L. Testing and analysing low-permeability fractured gas wells // J. Petrol. Technol.- 1968.-V. 20.- N 2.- P. 193-198.

183.    Morse R.A., Von Gonten W.D. Productivity of vertically fractured wells prior to the stabilized flow // Paper SPE 3631.- 1971.

184. Mousli N.A, Raghavan R, Cinco-Ley H, Samaniego V.F. The influence of vertical fractures intercepting active and observation wells on interference tests // Soc. Petrol. Eng. Journal.- 1982.- V. 22.- N 6.- P. 933-944.

185.    Mukherjee R, Paoli B.F., McDonald T., Cartaya H, Anderson J.A. Successful control of fracture height growth by placement of artificial barrier // SPE Prod. and Fac.- 1995.-V. 10.- N 2.- P. 89-95.

186. Mullen M.E, Norman W.D, Wine J.D, Stewart B.R.Investigation of

height grawth in frac-pack completions // Paper SPE 36458.- 1996.

187.    Nelson R.A. Fractured reservoirs: turning knowledge into prac-tice // J. Petrol. Technol. - 1987.- V. 39.- N 4. - P. 407-414.

188.    Nghiem L.X. Modeling infinite-conductivity vertical fractures with source and sink terms // Soc. Petrol. Eng. Journal.- 1983.- V. 23.- N 4. - P. 633644.

189.    Nghiem L.X, Forsyth P.A.Jr, Be hie A. A fully implicit hydraulic fracture model // J. Petrol. Technol.- 1984.- V. 36.- N 6.- P. 1191-1198.

190. Nolte K.G., Smith M.B. Interpretation of fracturing pressures // J. Petrol. Technol. -1981.- N 9.- P. 1767-1775.

191.    Nordgren R.P. Propagation of vertical hydraulic fracture // Soc. Petrol. Eng. Journal.- 1972.- V. 12.- N 4.- P. 306-314.

192.    Osman M.E. Pressure analysis of a fractured well in multilayered reservoirs // J. Petrol. Science and Eng.- 1993.- V. 9.- P. 49-66.

193.    Osman M.E, Abou-Kassem J.H. Effect of boundary conditions on pressure behavior of finite-conductivity fractures in bounded stratified reservoirs // J. Petrol. Science and Eng. - 1996.- V. 15.- P. 291-307.

194. Osman M.E., Abou-Kassem J.H. A new method for pressure test analysis of a vertically fractured well producing commingled zones in bounded square reservoirs // J. Petrol. Science and Eng.- 1997.- V. 18.- P. 131-145.

195. Overbey Jr. W.K., Yost II A.B, Wilkins D.A. Inducing miltiple hydraulic fractures from a horizontal wellbore // Paper 18249.- 1988.

196.    Ozkan E, Raghavan R. New solutions for well-test-analysis problems: Part 1 - Analytical considerations // SPE Form. Eval.- 1991.- V. 6.- N 3.-P. 359-368.

197.    Ozkan E, Raghavan R. New solutions for well-test-analysis problems: Part 2 -Computational considerations and applications // SPE Form. Eval.-

1991.- V. 6.- N 3.- P. 369-377.

198.    Peaceman D.W. Interpretation of well-block pressures in numerical reservoir simulation // Soc. Petrol. Eng. Journal. - 1978.- V. 18. - N 3. -P. 183-194.

199.    Peaceman D.W. Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability // Soc. Petrol. Eng. Journal.- 1983.- V. 23.- N 3.- P. 531-543.

200.    Peaceman D.W. Interpretation of well-block pressures in numerical reservoir Simulation: Part 3 - Off-center and multiple wells within a well-block // SPE Res. Eng. - 1990.- V. 5. - N 2. - P. 227-232.

201. Pearson C, Clonts M., Vaughn N. Use of longitudinally fractured horizontal wells in a multizone sandstone formation // Paper SPE 36454.

1996.

202. Pearson C.M, Bond A.J., Eck M.E., Schmidt J.H. Results of stress-oriented and aligned perforating in fracturing deviated wells // J. Petrol. Technol. - 1992.- V. 44.- N 1.- P. 10-18.

203. Perkins T.K, Gonzalez J.A. The effect of thermoelastic stresses on injection well fracturing // Soc. Petrol. Eng. J. - 1985.- V. 25. - N 2.- P. 78

87.

204. Perkins T.K, Kern L.R. Widths of hydraulic fracturing // J. Petrol. Technol. - 1961.- N 9.- P. 937-949.

205. Poollen H.K. van, Tinsley J.M, Saunders C.D. Hydraulic fracturing: fracture flow capacity as well productivity // Trans. AIME.- 1958.- V.

213.- P. 91-95.

206.    Prats M. Effect of vertical fractures on reservoir behavior - incompressible fluid case // Soc. Petrol. Eng. Journal.- 1961.- V. 1.- N 2.- P. 105

118.

207.    Prats M, Hazebroek P., Strickler W.R.Effect of vertical fractures on reservoir behavior - compressible fluid case // Soc. Petrol. Eng. Journal-1962.- V. 2.- N 2.- P. 87-94.

208.    Raghavan R.Some practical considerations in the analysis of pressure data // J. Petrol. Technol. - 1976. - V. 28. - N 10.- P. 1256-1268.

209.    Raghavan R, Joshi S.D. Productivity of multiple drainholes or fractured horizontal wells // SPE Form. Eval. - 1993.- V. 8.- N 1.- P. 11-16.

210.    Raghavan R, Reynolds A.C. Analysis of pressure buildup data following a short flow period // J. Petrol. Technol.- 1982. - V. 34. - N 4.- P. 904-916.

211.    Raghavan R, Uraiet A., Thomas G.W. Vertical fracture height: effect on transient flow behavior // Soc. Petrol. Eng. Journal.- 1978.- V. 18.- N

4.- P. 265-277.

212.    Rahim Z, Holditch S.A. Using a three-dimensional concept in a twodimensional model to predict accurate hydraulic fracture dimensions // Journal of Petroleum Science & Engineering.- 1995.- N 13.- P. 15 - 27.

213. Robinson B.M, Holditch S.A, Whitehead W.S, Peterson R.E. Hydraulic fracturing research in East Texas: third GRI staged field experiment // J. Petrol. Technol. - 1992.- V. 44.- N 1.- P. 78-87.

214.    Rodriguez F, Horne R.N., Cinco-Ley H. Partially penetrating vertical fractures: pressure transient behavior of a finite-conductivity fracture // Paper SPE 13057.- 1984.

215.    Roodhard L.P, Fokker P.A., Davies D.R, Shlyapobersky J, Wong G.K. Frac-and-pack stimulation: application, design, and field experience // J. Petrol. Technol. - 1994.- V. 46.- N 3.- P. 230-238.

216.    Russell D.G., Truitt N.E. Transient pressure behavior in vertically fractured reservoirs // J. Petrol. Technol. - 1964. - V. 16. - N 10. - P. 11591170.

217.    Schechter R.S. Oil well stimulation.- Prentice Hall. Englewood Cliffs. NJ, 1992. - 278 p.

218.    Settari A. A new general model of fluid loss in hydraulic fracturing // Soc. Petrol. Eng. Journal. - 1985. - V. 25. - N 4. - P. 491-501.

219.    Settari A. Simulation of hydraulic fracturing processes // Soc. Petrol. Eng. Journal. - 1980. - V. 20. - N 4.- P. 487-500.

220. Settari A, Bachman R.C., Hovern K.A. Paulsen S.G. Productivity of fractured gas-condensate wells: a case study of the Smorbukk field // SPE Res. Eng. - 1996. - V. 11. - N 4. - P. 236-244.

221.    Settari A, Cleary M.P. Three-dimensional simulation of hydraulic fracturing // J. Petrol. Technol. - 1984. - V. 36. - N 7. - P. 1177-1190.

222.    Settari A., Ito Y, Jha K.N. Coupling of a fracture mechanics model and a thermal reservoir simulator for tar sands // J. Canad. Petrol. Technol. -

1992. - V. 31. - N 9. - P. 20-27.

223. Settari A., Puchir P.J, Bachman R.C. Partially decoupled modeling of hydraulic fracturing processes // SPE Prod. Eng. - 1990.- V. 5.- N 1.-P. 37-44.

224. Simonson E.R, Abou Sayed A.S, Clifton R.J. Containment of massive hydraulic fractures // Soc. Petrol. Eng. Journal. - 1978.- N 1.- P. 27

32.

225.    Sleefe G.E., Warpinski N.R., Engler B.P. The use of broadband microseisms for hydraulic-fracture mapping // SPE Form. Eval.- 1995.- V.

10.- N 4.- P. 233-239.

226.    Smith M.B., Hannah R.R. High-permeability fracturing: the evolution of a technology // J. Petrol. Technol. - 1996. - V. 48.- N 6. - P. 628-633.

227.    Soliman M. Design and analysis of a fracture with changing conductivity // J. Canad. Petrol. Technol. - 1986. - V. 25. - N 5. - P. 62-67.

228.    Soliman M. Fracture conductivity distributions studied // Oil and Gas J. -1986.- V. 84.- N 6.- P. 89-93.

229. Soliman M, Hunt J.L.II, Azari M. Fracturing gorizontal wells in gas reservoirs // Paper SPE 35260.- 1996.

230.    Soliman M, Hunt J.L.II, El-Rabaa W. Fracturing aspects of gorizontal wells // J. Petrol. Technol. - 1990.- V. 42.- N 8. - P. 966-973.

231.    Strubhar M.K. Multiple, vertical fractures from an inclined well-bore - A field experience // J. Petrol. Technol. - 1975. - N 5. - P. 641-647.

232. Sullivan R.B., Lee W.J, Holditch S.A. Pressure transient response in multilayer gas reservoirs containing hydraulic fractures // Paper SPE/DOE 16399. - 1987.

233.    Suprunowicz R, Butler R.M. The effect of vertical fractures upon the performance of horizontal wells when coning can occur // J. Canad. Petrol. Technol. - 1996. - V. 35. - N 5.- P. 45 - 54.

234.    Tiab D. Analysis of pressure and pressure derivative without type-curve matching: Vertically fractured wells in closed systems // J. Petrol. Science and Eng. - 1994. - V. 11. - P. 323-333.

235.    Tirant P. le, Gay L. Manuel de fracturation hydraulique. - Paris: Tech-nip, 1972.- 334 pp.

236.    Tongchun Yi, Peden J.M. A comprehensive model of fluid loss in hydraulic fracturing // SPE Prod. and Fac. -1994. - V. 9. - N 4.- P. 267272.

237. Torok J.S, Advani S.H. Non-Newtonian fluid flow in a reservoir - an application to hydraulic fracturing // Trans. ASME. Journal Energy Resour. Technol.- 1987.- V. 109.- N 1.- P. 6-10.

238.    Tudor R, Poleschuk A. Low-viscosity, low-temperature fracture fluids // J. Canad. Petrol. Technol.- 1996.- V. 35. - N 7.- P. 31-36.

239.    Underwood P.J, Kerley L. Evaluation of selective vs. рoint-source perforating for hydraulic fracturing // Paper SPE 36480. - 1996.

240.    Uraiet A., Raghavan R, Thomas G.W. Determination of the orientation of a vertical fracture by interference tests // J. Petrol. Technol. -1977.- V. 29.- N 1.- P. 73-80.

241.    Vogel J.V. Inflow performance relationships for solution-gas drive wells // J. Petrol. Technol. - 1968. - V. 20. - N 1.- P. 83-92.

242.    Voneiff G.W, Holditch S.A. Economic assessment of applying advances in fracturing technology // J. Petrol. Technol. - 1994. - V. 46.- N 1.- P. 5157.

243. Voneiff G.W., Robinson B.M., Holditch S.A. The effects of unbroken fracture fluid on gas well performance // SPE Prod. and Fac. - 1996. - V. 11.-N 4.- P. 223-229.

244. Wagner P.B, Economides M.J, Ehlig-Economides C.A. Custom type-curve generation for pressure transient analysis of elongated linear flow systems // Revue de l’Inst. Fran. du Petrole. - 1984. - V. 39.- N 6.-P. 717-733.

245. Warpinski N.R., Moschovidis Z.A., Parker C.D., Abou-Sayed ,S.

Comparison study of hydraulic fracturing models - Test case: GRI staged field experiment No. 3 // SPE Prod. & Fac.- 1994.- N 1.- P. 7 - 18.

246. Wattenbarger R.A., Ramey H.J.Jr. Well test interpretation of vertically fractured gas wells // J. Petrol. Technol. - 1969.- V. 21.- N 5. - P. 625-632.

247.    Williams D.B., Sherrard D.W, Lin C.Y. Impact of inducing fractures in Pradhoe Bay // J. Petrol. Technol. - 1989. -V. 41. - N 10.- P. 10961102.

248. Wong D.W, Harrington A.B, Cinco-Ley H. Application of the pressure derivative function in the pressure transient testing of fractured wells // SPE Form. Eval. - 1986. - V. 1.- N 3. - P. 470-480.

ёеаёеа aaiOeAiieo

Борисенко З.Г. Методика геометризации резервуаров и залежей нефти и газа. — М.: Недра, 1980.

Викторин В.Д., Лыков Н.А. Разработка нефтяных и газовых месторождений, приуроченных к карбонатным коллекторам. — М.: Недра, 1982.

Гиматудинов Ш.К., Ширковский А.И. Физика нефтяного и газового пласта. — М.: Недра, 1981.

Гутман И.С. Методы подсчета запасов нефти и газа. — М.: Недра, 1985.

Дахнов В.Н. Интерпретация результатов геофизических исследований разрезов скважин. — М.: Недра, 1980.

Дементьев Л.Ф. Математические методы и ЭВМ в нефтегазовой геологии. — М.: Недра, 1983.

Жданов М.А. Нефтегазопромысловая геология и подсчет запасов нефти и газа. — М.: Недра, 1981.

Коротаев Ю.П., Закиров С.Н. Теория и проектирование разработки газовых и газоконденсатных месторождений. — М.: Недра, 1981.

Нефтегазопромысловая геология. Терминологический справочник: Под ред. М.М. Ивановой / 2-е изд. — М.: АО "ТВАНТ", 1994.

Справочник по нефтегазопромысловой геологии: Под ред. Н.Е. Быкова, М.И. Максимова, А.Я. Фурсова. — М.: Недра, 1981.

Спутник нефтегазопромыслового геолога: Под ред. И.П. Чоловского. — М.: Недра, 1978.

Сургучев М.Л. Вторичные и третичные методы увеличения нефтеотдачи пластов. — М.: Недра, 1985.

Чоловский И.П., Тимофеев В.А., Брагин Ю.И. Методы геологопромыслового контроля разработки нефтяных и газовых месторождений. 2-е изд. — Элиста: АПП "Джангар", 1996.

Иванова М.М., Дементьев Л.Ф., Чоловский И.П. Нефтегазопромысловая геология и геологические основы разработки месторождений нефти и газа. 2-е изд.— М.: Недра, 1992.

Гавура В.Е. Геология и разработка нефтяных и газонефтяных месторождений. — М.: ВНИИОЭНГ, 1995.

Геология и разработка крупнейших и уникальных нефтяных и нефтегазовых месторождений России / Под ред. В.Е. Гавуры (2 тома). — М.: ВНИИОЭНГ, 1 996.

Разработка нефтяных и нефтегазовых месторождений / Материалы совещания в г. Альметьевске в сентябре 1995 г. — М.: ВНИИОЭНГ, 1996.

Щелкачев В.Н. Избранные труды. 2 тома. — М.: Недра, 1990.

Дьяконов Д.И., Леонтьев Е.И., Кузнецов Г.С. Общий курс геофизических исследований скважин. — М.: Недра, 1984.

Кудинов В.И., Сучков Б.М. Новые технологии повышения добычи нефти. — Самарское книжное изд-во, 1998.

Регламент составления проектных технологических документов на разработку нефтяных и нефтегазовых месторождений. — М.: Минтопэнерго России, 1996.

Каналин В.Г., Вачин С.Б., Токарев М.А. Нефтегазопромысловая геология и гидрогеология. — М.: Недра, 1997.

Токарев М.А. Комплексный геолого-промысловый контроль за текущей нефтеотдачей при вытеснении нефти водой. — М.: Недра, 1990.

Фурсов А.Я. Оптимизация изученности нефтяных месторождений. — М.: Недра, 1985.

Панов Г.Е., Петряшин Л.Ф., Лысянский Г.Н. Охрана окружающей среды на предприятиях нефтяной и газовой промышленности. — М.: Недра, 1986.

Булыгин В.Я., Булыгин Д.В. Имитация разработки залежей нефти. — М.: Недра, 1 990.

Лысенко В.Д. Разработка нефтяных месторождений. Теория и практика. — М.: Недра, 1996.

Лебединец Н.П. Изучение и разработка нефтяных месторождений с трещиноватыми коллекторами. — М.: Наука, 1997.

Производительность скважин  »
Библиотека »