Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Романков, П. Г. Гидромеханические процессы химической технологии

.pdf
Скачиваний:
29
Добавлен:
22.10.2023
Размер:
10.97 Mб
Скачать

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

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

 

 

 

 

 

 

 

 

ЛИТЕРАТУРА

 

1. К и р п и ч е в

М. В. Теория

подобия. М. — Л., Изд. АН СССР, 1953,

95 с.;

 

Г у X м а н А.

А. Введение в теорию подобия. Изд. 2-е, М., «Высшая школа»,

 

1973, 295 с.;

В е н и к о в

В.

А. Теория

подобия — моделирование. М.,

«Выс­

 

шая школа»,

1966,

487 с.;

Д ь я к о н о в

Г. К. Вопросы теории подобия

в об­

2.

ласти физико-химических

процессов. М. — Л., Изд. АН СССР, 1956,

206 с.

Д р у ж и н и н

Н. И. Метод

электрогидродинамических аналогий и его при­

 

менение при исследовании фильтрации. М .— Л., Госэнергоиздат, 1959, 346 с.;

 

Т е т е л ь б а у м

И.

М. Электрическое

моделирование. М., Физматгиз,

1959,

 

319 с.; Г у т е н м а х е р

Л.

И. Электрические модели. М. — Л., Изд. АН

3.

СССР, 1949, 402 с.

 

 

 

 

 

К а ф а р о в

 

В. В. Методы кибернетики в химии и химической технологии.

 

М., «Химия»,

1968,

379 с.;

Ф р э н к с Р.

Математическое моделирование в хи­

 

мической технологии. Пер. с англ. Под ред. В. С. Торопцова. М., «Химия»,

 

1971, 279 с.;

К а ф а р о в

В.

В. Моделирование химических процессов. М.,

 

«Химия», 1968, 350

с.; Д о л е ж а л и к В. Подобие и моделирование в хими­

 

ческой технологии. Пер. с чешского. Под ред. Н. И. Гельперина. М., Гос-

4.

топтехиздат,

1960, 96 с.

Методы моделирования каталитических процессов

С л и н ь к о

М.

Г.

и др.

 

на аналоговых

и

цифровых

вычислительных машинах. Новосибирск, «На­

 

ука», 1972,

150

с.; К а ф а р о в В. В., Е р е м е н к о В. Е. В кн.: Кибернетику

5.

на службу коммунизму. М., «Энергия», 1967, с. 6—9.

 

Л е н и и В. И. Соч.,

Изд. 4-е,

т. XIV, с. 276.

Пер.

6.

В е н е д е к

П.,

Л а с л о

А.

Научные

основы химической технологии.

 

с нем. Под ред. П. Г. Романкова и М. И. Курочкиной. Л., «Химия», 1970,

 

376 с.

 

О. Proc. Lit. Phil. Soc. (Manchester), 1874, p. 14.

 

7. R e y n o l d s

 

 

8.Г у X м а II А. А. Применение теории подобия к исследованию процессов тепло-массообмена. М., «Высшая школа», 1967, 303 с.

9.L е w i s W. Mech. Eng., 1922, v. 44, p. 445—448.

10.

Л ы к о в

А.

В. Теория теплопроводности. М., «Высшая школа»,

1967,

599 с.;

 

М у ч н и к Г.

Ф.,

Р у б а ш о в И.

Б.

Методы теории теплообмена,

ч. 1,

Тепло-

И.

ароводность.

М.,

«Высшая школа», 1970, 285 с.

учения о теплообмене.

Пер.

Г р ё б е р

Г.,

Э р к С.,

Г р и г у л ь

У.

Основы

12.

с нем. Под ред. А. А. Гухмана, М. — Л.,

Издатинлит, 1958, 566 с.

Гостехтеор-

Б р и д ж м е н

П. В. Анализ размерностей. Пер.

с англ. М. — Л.,

 

издат, 1934,

120

с.; Б р а й н е с

Я.

М.

Подобие

и моделирование в химиче­

13.

ской промышленности. М. — Л., Гостоптехиздат,

1961, 220 с.

 

 

С е д о в

Л. И. Методы

подобия

и

размерностей в механике. Изд. 5-е. М„

 

«Наука», 1965, 320 с.;

К у р о ш

А.

Г.

Теория

групп. М. — Л.,

Госте.хіпдят,

 

1953, 467 с.

 

 

 

 

 

 

 

 

 

47

14. А л а б у ж ев

П. М. и др. Теория подобия и размерностей. Моделирование.

М., «Высшая школа», 1968, 206 с.

15. Р о з е н б р о к

X., С т о р и С. Вычислительные методы для инженеров-

химиков. Пер. с англ. М., «Мир», 1968, 442 с.

16.Б о я р и н о в А. И., К а ф а р о в В. В. Методы оптимизации в химической технологии. М., «Химия», 1969, 564 с.

17.

Л е в е н ш п и л ь

О. Инженерное оформление химических процессов. Пер.

18.

с англ. Под

ред. и с дополн. М. Г. Слинько. М., «Химия», 1969, 621 с.

А р с е н ь е в

Ю. Д. Теория подобия в инженерных экономических расчетах.

19.

М., «Высшая школа», 1967, 261 с.

Н а л и м о в

В.

В., Ч е р н о в а Н. А. Статистические методы планирования

экстремальных экспериментов. М., «Наука», 1965, 340 с.

20. Р о з е н А. М. В кн.: Доклады V межвузовской конференции по физическому и математическому моделированию. Секция моделирования химико-техноло­

гических процессов.

М., МЭИ АН СССР, 1968, 173 с.;

Р о з е н А.

М., К р ы ­

л о в В. С. ТОХТ,

1967, т. I, №

3, с. 297—300; Б о р е с к о в

Г. К-,

С л и н ь ­

ко М. Г. Там же, № 1, с. 5—7.

Ю. М. Методы оптимизации химических

21. О с т р о в с к и й

Г.

М., В о л и н

реакторов. М. — Л., «Химия», 1967, 248 с.; Л и п а т о в Л. Н. Типовые

про­

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

«Химия»,

1973,

317 с.; К р о у

К.,

Т а м и л е ц

А., Х о ф ф м а н Т.

и др.

Математическое

моделирование

химических производств. М., «Мир»,

1973,

392 с.; К а ф а-

р о в В. В., ТОХТ, 1972, т. 6, № 6, с. 880—892 и 908—915.

Г л а в а 3

ПРОЦЕССЫ, СОСТАВЛЯЮЩИЕ ВНУТРЕННЮЮ ЗАДАЧУ ГИДРОДИНАМИКИ (движение жидкостей и газов по трубам и каналам)

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

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

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

УРАВНЕНИЕ НЕРАЗРЫВНОСТИ

В гидравлике для потоков газа или жидкости законы сохране­ ния представляются в виде уравнения неразрывности (сплошно­ сти) .

Уравнение неразрывности потока носит также название уравнения постоянства рас­ хода

dfiWi

df2w2

df 2w3 =

dVсек

(3- 1)

 

 

Здесь dfi,

df2, df3— элемен­

 

 

тарные

площади

сечений

I—I,

 

 

II—II,

 

III—III

потока

(рис.

 

 

3-1); Wu Wz, w3— скорости по­

 

 

тока

в

указанных

сечениях;

 

 

dVсек— объемный

расход.

 

(3-1) получим:

 

После интегрирования уравнения

 

 

 

 

 

 

Лиц = / 2а>2= / 3ш 3= Усек

(3-2)

где fi,

/2, /з — площади поперечных

сечений трубопровода,

м2; wu

wz, Щ — скорости потока,

м/с; ѴСек— расход жидкости, м3/с.

49

I

Уравнение расхода для трубопроводов круглого сечения удобно применять в следующей форме:

7 Сек = 0,785d 2w

где d — диаметр трубопровода, м; w — средняя скорость потока,

м/с.

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

пространстве, но и во времени. В частности,

скорость w элемен­

 

тарной частицы можно пред­

 

ставить в виде зависимости

 

 

 

w = ф {х, у, г, х)

(3-3)

 

где

X,

у, z — координаты; т —

 

время.

элементарный

объем

 

 

В

 

dV = dxdydz

(рис. 3-2)

за вре­

 

мя dt вдоль оси X поступит че­

 

рез

грань

dy dz

количество

 

массы жидкости (газа),равное

Рис. 3-2. К выводу уравнения нераз­

pWxdVdr [величину рwx следует

рывности потока.

рассматривать

как

плотность

 

х-компоненты

количества дви­

жения в кг-м/(с-м3)]. За это же время dx из элементарного объе­ ма в направлении оси х на расстоянии х + dx выйдет количество массы, равное

, д wx ) ,

9wx + ~ ^ ~ d x

Накопление массы жидкости в элементарном объеме в направ­ лении оси X составит:

рwx dx dy dz dx pwx -f d dx 1 dy dz dx = d dx dy dz dx (3-4)

Аналогично в направлении оси у накопление массы будет равно

d{pwy)

<Э(рwz)

— ^ — d y d x d z d г,

а в направлении оси z -------— dzdxdydx .

По закону сохранения массы за время dx суммарное накопле­ ние массы жидкости (газа), движущейся по всем трем направле­ ниям элементарного объема dV должно быть равно убыли массы дрдх dx dy dz dx в этом объеме.

Таким образом, материальный баланс неустановившегося по­ тока сжимаемой жидкости (газа) можно выразить зависимостью,

50

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

др

д ( р wx) |

d ( p W y ) '

д (pwz)

(3-5)

d r =

дх

ду

dz

 

Величины, входящие в правую часть уравнения (3-5), представ­ ляют собой дивергенцию (расхождение) плотности потока.

Уравнение (3-5) можно привести к виду

I dwx

dwу

dwz \

dp

dp

dp

dp

(3-6)

P V âx

ây

dz )

W* dx

Wy dy

Wz dz

âx

 

Поскольку плотность жидкости р =

/' (х,

у, z, т), т. е. р зависит

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

 

dp =

-^2. ух

dy

уу _), ^£. yz

dx

yr

(3-7)

v

âx

J dz

 

 

 

Полная производная по времени будет

 

 

 

 

dp_

dp_'dx_

, _dp

 

 

 

 

 

dx

dx dx

dy

dx

dz

dx

dx

’ '

и при dx/dx = dy/dx = dz/dx = 0 скорость изменения плотности в данной точке

dp _ âp_

(3-9)

dx dx

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

В векторно-аналитической форме уравнение (3-5) будет иметь вид:

div (ри>) + — = 0

(3-10)

Для установившегося (стационарного) потока dp/dx — 0, т. е. плотность жидкости (газа) не зависит от времени, и уравнение неразрывности можно записать в виде

â{pwx)

t

â(pwy)

_

d( pwz)

л

fn ііч

dx

+

dy

+

dz

 

(3' '

или в векторной форме:

div(p® )=«0

(3-12)

В

потоке

несжимаемой

жидкости плотность

постоянна

(р =

const) и

уравнение неразрывности упростится:

 

 

 

dwx

dw,,

dwz

(3-13)

 

 

dx ^

dy

' dz

или

 

 

 

div (pw) =

p div w

(3-14)

 

 

61

Для одномерного потока сжимаемой жидкости (газа), направ­ ленного вдоль оси X и проходящего через сечение /, уравнение неразрывности можно выразить в виде

d (P WxU.+ f £р = о

(3-15)

д х ^ ' д х

 

к ’

или для установившегося потока

 

 

d(pwx)f

0

(3-16)

дх

 

 

 

Откуда

 

(3-17)

pwxf —const

Таким образом, для двух произвольных сечении потока

 

= Р2®г/2

(3-18)

Для несжимаемой жидкости выражения (3-17) и (3-18)

упро-

стятся:

 

(3-19)

wf = const

 

Отсюда Wifi = w2f2 = . . . , т. е. мы

получим уравнение постоян­

ства расхода (3-2), которое представляет наиболее употребитель­ ную в технике форму уравнения неразрывности.

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

УРАВНЕНИЯ ДВИЖЕНИЯ

При изучении законов равновесия и движения жидкостей и га­ зов (а также взаимодействия движущихся жидкостей и газов с об­ текаемыми имй" твердыми телами) в гидромеханике жидкость или газ ** рассматривают как непрерывно распределенную в простран­ стве сплошную среду. Движение жидкости или газа осуществляет­ ся под действием массовых (объемных) и поверхностных сил.

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

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

*Идеальной жидкостью называется жидкость, в которой отсутствует внут­ реннее трение.

**Это допущение неприменимо к разреженным газам.

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

виях

равновесия касательные напряжения * в

жидкости равны

нулю

и поверхностные силы выражаются лишь

силами давления

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

Таким образом, при выводе уравнений движения жидкости в основу берется второй закон движения Ньютона, согласно кото­ рому изменение скорости движения во времени пропорционально

действующей силе и

имеет одинаковое

 

с нею направление. Другими словами,

 

векториальная сумма всех сил, действую­

 

щих на

выделенный

элемент жидкости

 

(рис. 3-3), равна произведению его мас­

 

сы на ускорение (принцип Д’Аламбера).

 

Баланс действующих в потоке сил вы­

 

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

 

жидкости уравнениями Эйлера, а в слу­

 

чае движения реальной жидкости — урав­

 

нениями Навье — Стокса.

 

 

Рассмотрим самый общий случай и

 

составим

уравнения

движения вязкой

 

сжимаемой жидкости.

силы тяжести

сов­

Рис. 3-3. К выводу уравне­

Если

направление

ния движения жидкости.

падает с направлением оси х (рис.

3-3),

 

то проекция на ось х массовых и поверхностных сил, действующих на элементарный объем dV = dxdydz жидкости, может быть полу­

чена с учетом силы внутреннего трения по уравнению

Ньютона

(для одномерного потока)

 

 

То

dwx

(3-20)

dz

 

 

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

Проекция силы трения на ось х будет иметь вид:

То + ^т0 + - | j - dxj dy dz

(3-21)

Тогда движение элементарного объема под действием сил тяже­ сти, давления и трения вдоль оси х может быть описано уравне­ нием

p^

= p* * “ 4 £ + ^V2* * + T ^ (dlvw)

(3'22)

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

S3

движение параллельно оси у

ÜWy

dp

 

I

д

 

(3-23)

Р Ч Г = Му -

т

+ PV2^

+ 3 llW

(div w)

 

и оси z

dp

,

. 1

<?

,

 

Dwz

(3-24)

p -dT -= p ^

- i

+ ^ +

i

^

(dlvw)

 

В левой части уравнений

(3-22) — (3-24)

производная Dw/dx

представляет собой полную производную скорости по времени:

Dwx

dwx

+ wx

dwx

dWy

 

dw2

dx

dx

dx

+ Wy ~дГ +- Wz

dz

Dw,j

dWy

+

dwx

dWy

+ wz

dwz

dx

dx

dx

+ wy~ d f

dz

Dwz

dwz

 

dwx

dWy

 

dwz

dx

 

 

dx

+ wy dy

+ wz ~dz

Величины gx, gy и gz — компоненты ускорения массы в направ­ лении отдельных осей координат (например, gy = g cos ß, где ß — угол между направлением действия силы тяжести и осью у).

Влияние вязкости р и сжимаемости жидкости учитывается по­ следними двумя членами правой части уравнений (при условии, что вязкость во всей области течения постоянна, т. е. р = const).

Оператор Лапласа обозначен символом V2 (набла) и представ­

ляет собой сумму вторых производных *:

 

ѵ2 =

д2

д2

 

dx2

âz2

Сжимаемость реальной жидкости учитывается последним чле­ ном правой части уравнений. Дивергенция вектора скорости обо­

значает:

 

âWx

dWy

 

 

 

dlу w

 

 

 

dx

dy

dz

 

 

 

 

В векторной форме уравнения

(3-22) —(3-24) сводятся к одному

уравнению:

 

 

 

 

 

dw

= Р * . grad р + рѴ2и> + -д- Р grad (diva>)

(3-25)

dx

 

 

 

 

 

где R — вектор напряженности массовых сил (для силы тяжести R

равно ускорению свободного падения g);

р — давление; р — вяз­

кость; р — плотность; т — время;

w — вектор скорости.

 

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

Оператор Лапласа часто обозначают также через Д,

64

Для неустановившегося движения несжимаемой жидкости (р=

=

const)

последние

члены

в

уравнениях

(3-22)— (3-25)

отсут­

ствуют, так как div "ш=

0,

и уравнения движения упрощаются:

 

 

р

dx

=

Pgx -

дх + цѴ2шX

 

 

 

 

Dwy

 

 

др

 

(3-26)

 

 

р

dx

=

PBy -

—--- (- рѴ2ш

 

 

ду

У

 

 

 

р

Dwz

 

 

+ рѴ2аі>Z

 

 

 

dx

=

Р£г -

 

или

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

grad р + рѴ2а>

(3-27)

 

Уравнение (3-27) может быть выражено в форме, не содержа­

щей давления *:

 

 

 

 

 

 

 

 

 

rot w =

rot

+

rot [w rot w] +

vV2 rot w

(3-28)

 

 

OX

 

 

 

 

 

 

 

где

V =

p/p — кинематический

коэффициент вязкости; rot яу— ро­

тор (или вихрь) вектора скорости.

 

ѵ = 0.

 

В поле тяготения

rot R =

0,

для идеальной жидкости

Тогда уравнение (3-28) упростится:

 

 

 

 

 

 

rot w — rot [ги rot а>]

 

(3-29)

 

Уравнения движения

жидкости Навье — Стокса (3-22) — (3-24)

или (3-25) совместно с уравнением неразрывности (3-5) или (3-10) дают возможность решить основную задачу гидродинамики — оп­ ределить поля скоростей, давления и плотности в жидкости, дви­ жущейся под действием заданных внешних сил. Решение возможно для идеальной несжимаемой жидкости или для изотермического движения вязкой жидкости, когда плотность и вязкость зависят только от давления и вид этих зависимостей известен: р — f(p) и

рі = f'(p). В этом случае

возможно определение следующих зави­

симостей:

 

 

 

wx = fi(x, у, z, т);

Wy =

f2{x, y, z, т); wz = /

3(х, у, г, т)

P = fi (х, у, г,

т); р = U (х, у, г, х)

 

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

* Так как в уравнения движения жидкости давление входит под знак диф­ ференциала, то следует отметить, что для гидромеханических процессов харак­ терными являются не абсолютные давления, а разность давлений в двух какихлибо точках потока. Таким образом, уравнения движения жидкости позволяют исследовать вопрос о распределении скоростей в потоке вне связи с распре­ делением давлений, причем обратное исследование поля давлений без предвари­ тельного установления поля скоростей невозможно [1].

55

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

Решение уравнений движения Навье — Стокса получено только для некоторых простейших случаев одномерного или двумерного потока, например, для течения вязкой жидкости по прямой трубе (задача Пуазейля), для течения между двумя плоскими парал­ лельными стенками, одна из которых неподвижна (задача Куэтта), а также при обтекании неподвижной тонкой пластинки (в этом

случае уравнения Навье — Стокса оказывается возможным

заме­

нить более простыми уравнениями пограничного слоя).

не от­

Следует отметить также, что уравнения Навье — Стокса

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

РЕЖИМЫ ДВИЖЕНИЯ ЖИДКОСТИ (ГАЗА)

Экспериментально установлено [2], что существуют два режима движения потока — ламинарный и турбулентный. Ламинарное, или слоистое, движение наблюдается при малых скоростях или в тру­ бах малого диаметра. При ламинарном движении слои жидкости скользят один относительно другого не перемешиваясь. В условиях установившегося движения скорость w при ламинарном режиме постоянна в каждой точке потока, т. е. w = f(x, у, z).

Рис. 3-4. Пульсации скорости потока (о>ср = w ± Ада).

При высоких скоростях отдельные слои (струйки) потока бес­ порядочно перемешиваются между собой и в каждой точке потока (даже в условиях установившегося движения) имеют место бы­ стрые изменения скорости пульсации около некоторого ее среднего значения wcv. Такой режим движения называется турбулентным.

Интенсивность пульсаций служит мерой турбулентности потока. Пульсации или пульсационные скорости, представляющие собой отклонения действительной мгновенной скорости w от среднего (по времени) значения скорости потока шср (рис. 3-4), можно разло­ жить на отдельные составляющие в направлении осей координат:

Соседние файлы в папке книги из ГПНТБ