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

2) Межфазный перенос

.pdf
Скачиваний:
7
Добавлен:
10.10.2020
Размер:
1.12 Mб
Скачать

2.2 Межфазный перенос субстанции

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

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

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

2.2.1. Уравнения массо -, тепло- и импульсоотдачи

2.2.1.1. Локальная форма уравнений

Рассмотрим элементарный участок межфазной поверхности dF, совпадающей с плоскостью xoy. Поток субстанций направлен вдоль оси z, движение фазы - по оси x (рис.2.5).

z

wx

j дг

q тг

вг

 

iz

z

zx

x

 

 

 

dF

y

Рис 2.5. Перенос субстанций по оси z

Рассмотрим поток субстанции за счет молекулярного и турбулентного механизмов переноса:

-jizдг - диффузионный поток массы;

-вгzx - вязкий поток импульса (тензор вязких напряжений);

-qzтг - поток тепла за счет теплопроводности.

Проекция теплового потока за счет теплопроводности на ось z по закону Фурье имеет вид:

q тг

 

 

dT

 

.

(2.64)

т

 

z

 

dz

 

z 0

 

 

 

 

 

 

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

В тепловом пограничном слое δт температура среды меняется от Т г

(температура поверхности раздела фаз) доТ я (температура на внешней границы пограничного слоя, т.е. температура ядра). В ядре фазы температура не меняется. По

закону Ньютона тепловой поток qzтг может быть записан:

q тг =α( Т г

Т я ),

(2.65)

z

 

 

где α – коэффициент теплоотдачи. Коэффициент теплоотдачи зависит от многих факторов: режима движения и физических свойств среды, геометрических параметров каналов и т.д.

Аналогичным образом могут быть получены уравнения массо- и

импульсоотдачи:

 

 

 

 

 

 

 

jдг β' г μя ) β

i

(сг ся ) ,

(2.66)

iz

i

i

i

 

i i

 

 

τвг

γ(wг

wя ) ,

(2.67)

 

zx

 

x

 

 

x

 

где βi, γ – коэффициенты массо- и импульсоотдачи.

Разница значений субстанций у границы раздела фаз и в ядре фазы носит название движущей силы субстанцииотдачи.

Коэффициенты массо-, тепло-, импульсоотдачи определяются:

 

 

 

 

j дг

 

 

 

 

dM г

 

 

 

 

 

βi

 

iz

 

 

i

 

 

,

 

(сг

ся )

dF dt(сг

 

 

 

 

 

 

 

 

ся )

 

 

 

i

i

 

 

 

 

i

 

i

 

 

 

 

тг

 

 

 

 

dQг

 

 

 

 

 

α

 

qz

 

 

 

 

 

 

 

,

 

(Т г Т я )

dF dt(T г T я )

 

 

 

 

τвг

 

 

 

 

dP

г

 

 

 

 

γ

 

 

 

zx

 

 

x

 

 

 

.

(wг

wя )

dFdt(wг

 

wя )

 

 

 

x

x

 

 

 

 

x

 

x

м/с

(2.68)

Вт/(м2К)

(2.69)

кГ/ (м2с)

(2.70)

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

Коэффициенты массоотдачи рассмотрены для бинарных сред.

При ламинарном течении жидких сред вместо значения переменной в ядре потока в уравнениях (2.65) – (2.70) используют осредненное по поперечному

сечению значение. Для ламинарного режима течения модель пограничного слоя «не работает».

2.2.1.2. Интегральная форма уравнений

Интегральная форма уравнений межфазного переноса субстанций получают осреднением локальных уравнений по всей межфазной поверхности F:

M& iг

Q& г

P&xг

 

 

 

 

г

F

 

 

 

 

dM i

 

 

 

г

я

 

 

 

 

 

 

òβi (ci

ci )dF

dt

 

 

 

 

 

 

 

 

 

0

 

 

dQг

 

 

F

 

 

 

 

α(T г

T я )dF

 

 

dt

 

0

 

 

 

 

 

 

 

 

 

 

 

 

dPxг

 

 

 

F

 

 

 

 

 

γ(Wxг

Wxя )dF

 

 

 

dt

 

0

 

 

 

 

 

 

 

 

 

 

 

βi F (ciг ciя ) ,

αF (T г T я ) ,

γF (Wxг Wxя ) .

(2.71)

(2.72)

(2.73)

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

Например, независимо осредним движущую силу:

 

 

 

 

 

 

 

 

 

 

 

 

1

 

F

 

 

 

 

 

 

 

 

(T г T я )

г T я )dF ,

 

 

 

 

 

(T

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tг - T я F .

 

 

 

 

 

 

dQт

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для случая модели гидродинамической структуры потока

вытеснения (МИВ) получено:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tвх Tвых

,

 

 

(T г T я ) T

 

 

 

 

 

 

 

 

 

 

 

ср

 

 

 

 

 

Tвх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вых

где (T г T я

) T

,

(T г T я

) T .

 

 

 

 

вх

вх

 

 

вых

 

 

 

 

вых

 

 

 

 

Использование диффузионной (MD) и ячеечной моделей (МЯ) более сложным зависимостям для Тя(F) и (T г T я ) .

Т

Тявых МИВ

(2.74)

(2.75)

идеального

(2.76)

приводит к

я МИС

Тявых MD,

Рис 2.6. Изменение температуры в ядре потока по длине аппарата для различных моделей

Как видно из рис. 2.6., максимальную среднюю движущую силу и, соответственно, эффективность теплообменного аппарата обеспечит структура потока, соответствующая МИВ, а минимальную – соответствующая модели идеального смещения (МИС). Диффузионная и ячеечная модель дают промежуточные результаты.

2.2.2 Уравнения массо-, тепло- и импульсопередачи

2.2.2.1 Локальная форма уравнений

Рассмотрим перенос субстанций из фазы 1 через межфазную поверхность в фазу 2 за счет молекулярного и турбулентного механизмов. Примем, что сопротивлением переносу субстанции со стороны межфазной поверхности можно пренебречь. Это равносильно предположению об установлении равновесия на границе раздела фаз, т.е.:

г

г

, T г T г , wг wг .

(2.77)

i1

i2

1

2

1

2

 

y

2 фаза

межфазная

поверхность

1 фаза

x

Рис 2.7. Схема межфазного переноса субстанций.

Предположим μi1i2, тогда;

j дг i1 я

г

iy

 

i1

i1

.

j дг

i2

г

я

 

iy

i2

i2

 

Разделим уравнения на βi1 и βi2 соответственно и их сложим:

 

 

1

 

1

 

 

 

 

j дг

 

 

 

 

я

я

,

 

 

iy

 

i1

 

i2

 

i1

i2

 

 

 

 

 

 

 

 

j дг

 

1

 

1

 

1

я

я

K

 

я

я

.

 

 

 

 

 

 

iд

 

 

 

iy

 

i1

 

i2

 

 

i1

i2

 

i1

i2

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.78)

(2.79)

Здесь K- коэффициент массопередачи, (μi1я μi2я ) - движущая сила

массопередачи. Уравнение (2.79) носит название уравнения массопередачи. Химические потенциалы неидеальных (реальных) систем определить

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

Аналогичным образом могут быть получены уравнения тепло- и импульсопередачи:

 

qтг

K

 

(T я T я ), K

 

 

 

 

 

 

 

 

1

 

 

,

(2.80)

 

т

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

1

2

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α1

 

 

α 2

 

 

 

если Т12.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τвг

K

 

(wя

wя

) , K

 

 

 

 

1

 

 

 

 

,

 

 

 

(2.81)

г

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yx

 

x1

x2

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ1

 

γ2

 

 

 

 

 

 

если wx1>wx2.

Здесь Кт и Кг – коэффициенты тепло- и импульсопередачи. Коэффициенты в соотношениях (2.79) – (2.81) могут быть представлены иначе:

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

Ki д

βi1

 

βi2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

1

 

,

,

(2.82)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kт

 

 

 

 

α1

 

 

α2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

Kг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ1

 

 

γ2

 

 

где

 

1

,

1

,

1

- сопротивления

 

 

 

массо-, тепло-,

импульсопередачи

 

 

 

 

 

 

 

 

K

Kт

Kг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(межфазные

сопротивления), а

1

,

1

и

1

 

-

 

сопротивления

массо-, тепло- и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

βi

α

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

импульсоотдачи (фазовые сопротивления).

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

1

 

1

 

1

r

,

(2.83)

 

 

 

Kт

 

α1

 

ст

 

 

α2

 

 

где rст – термическое сопротивление стенки.

Профили wx, Т, μi в процессе переноса субстанции через границу раздела фаз, не обладающую сопротивлением, приведены на рис. 2.8.

y

wяx2

Тя2

μяi2

2 фаза

δг2

δт2

δд2

i

дт

q

тг

 

вг

 

 

 

 

y

yx

 

 

 

ij

 

 

 

wг

Т г

μiг

 

 

 

 

межфазная

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поверхность

 

δг1

δТ1

δд1

 

 

 

 

 

1 фаза

wяx1

Тя1

μяi1

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

Здесь δ – толщина пограничных слоев.

Если сопротивление одной из фаз, например первой, гораздо больше второй, то последним можно пренебречь:

1

 

1

,

1

 

1

,

1

 

1

.

(2.84)

 

 

 

 

 

 

Kig

 

βi1

Kт

 

α1

Kг

 

γ1

 

Из (2.84) следует, что Kβi1, Kт α1, Kг γ при βi1 << βi2, α1 << α2, γ1 << γ2.

Интенсификация процессов переноса требует увеличения коэффициентов субстанциипередачи. Для этого необходимо увеличить наименьший коэффициент субстанцииотдачи.

2.2.2.2 Интегральная форма уравнений

Осреднив локальные уравнения межфазного переноса субстанций по участку поверхности F можно получить интегральную форму уравнений:

 

 

 

 

 

 

 

 

г

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

г

 

 

 

dM i

 

 

дг

 

 

 

 

 

 

 

 

 

 

я

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M i

 

dt

 

 

jiy

dF KF i1

μ i2 ) ,

(2.85)

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

& г

 

 

dQг

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тг

 

 

 

 

 

 

 

 

 

я

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

dt

qy

dF K т F (T1

T2 ) ,

(2.86)

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dPг

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&г

 

 

 

 

вг

 

 

 

 

 

 

 

 

 

я

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Px

 

 

 

dt

τyxdF K г F (Wx1

Wx2 ) .

(2.87)

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вобщем случае при одновременном изменении кинетического коэффициента

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

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

-прямоток (движение фаз в одном направлении), -противоток (движение фаз в противоположных направлениях), -перекрестный ток, -смешанный ток.

1

1

 

1

2

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

прямоток

противоток

перекрестный

 

 

смешанный

 

 

 

ток

 

 

 

 

 

ток

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

2.3. Моделирование технологических процессов

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

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

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

Объект-оригинал – объект, свойство которого подлежат изучению методом моделирования.

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

Наибольшее распространение в инженерной практике получила математическое и физическое проектирование.

2.3.1. Математическое моделирование

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

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

Любая модель неполно отображает оригинал. Поэтому следующим этапом моделирования является проверка адекватности модели - соответствия ее моделируемому объекту. Это достигается путем сопоставления результатов моделирования с численным либо физическим экспериментом.

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

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

Итак, этапы математического моделирования:

-составление математической модели;

-идентификация модели;

-проверка адекватности модели, при необходимости коррекция;

-использование модели для описания объекта-оригинала.

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

2.3.2. Физическое моделирование

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

-какую модель использовать (форма, размер, среда),

-какие характеристики измерять,

-как перенести результаты исследования с модели на объект.

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

2.3.2.1. Теория подобия

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

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

Геометрическое подобие – постоянство отношения всех сходственных линейных размеров модели и оригинала.

 

 

 

 

 

l m

 

l m

 

 

 

 

 

 

 

 

 

 

1

 

2

K

l

const,

(2.88)

 

 

 

 

 

 

 

 

 

 

 

l 0

 

l 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

где

l m

и

l 0

- сходственные линейные размеры модели

и объекта; K

l

-

 

1

 

1

 

 

 

 

 

 

 

 

 

константа геометрического подобия.

Временное подобие (гомохранность) – постоянство отношения сходственных интервалов времени модели и оригинала:

t m

 

t m

 

 

 

1

 

2

K

t

const .

(2.89)

t 0

t 0

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

Если Kt =1, то имеем синхронность.

Подобие физических величин – постоянство отношения физических величин для модели и оригинала в сходственных точках в сходственные моменты времени:

ρm

Kρ ,

μm

Kμ ,

λm

Kλ .

(2.90)

ρ0

μ0

λ0

 

 

 

 

Подобие модели и объекта предполагают подобие полей физических величин:

wm K w - гидродинамическое подобие (подобие полей скоростей); w0

Т m Kт - тепловое подобие (подобие полей температуры);

Т 0

c m K c - концентрационное подобие (подобие полей концентраций). c 0

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

Подобие граничных условий – постоянство отношения соответствующих величин на границах модели и оригинала.

Константа подобия – отношения одноименных величин модели и оригинала. Они постоянны для различных сходственных точек подобных систем.

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

Симплексы подобия – инварианты подобия, представляющие собой отношения однородных величин:

l m

 

l m

 

 

 

1

 

2

idem Г

l

.

(2.91)

l 0

l 0

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

Критерии подобия – инварианты подобия, представляющие собой отношения разнородных, сложных величин:

Например,

wl м

wl 0

idem .

(2.92)

м

0

 

 

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

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

Теоремы подобия:

1. Подобные объекты характеризуются численно равными критериями подобия.