close

Вход

Забыли?

вход по аккаунту

?

Динамика адсорбции с учетом химической реакции при очистке стоков плавательных бассейнов.

код для вставкиСкачать
УДК 544.35.03:665.7
О. Р. Каратаев, С. В. Анаников
ДИНАМИКА АДСОРБЦИИ С УЧЕТОМ ХИМИЧЕСКОЙ РЕАКЦИИ
ПРИ ОЧИСТКЕ СТОКОВ ПЛАВАТЕЛЬНЫХ БАССЕЙНОВ
Ключевые слова: модель, адсорбция, массоперенос, преобразование Лапласа, оригинал, изображение.
В статье рассматривается модель адсорбции хлорорганических соединений на цеолитах с учетом химической
реакции. Задача решается для случая, когда основное сопротивление массопереносу сосредоточено во внешней
фазе. Решения получены методом преобразования Лапласа. Они могут быть распространены на аналогичные
задачи при кусочно-линейной аппроксимации нелинейной изотермы адсорбции, осложнённой химической
реакцией первого порядка.
Keywords: model, adsorption, masstransfer, transformation of Laplace, original, representation.
In the article is descrebed model adsorption of chlororganic compaunds on ceolits with chemical reaction of first
order. It is taked into consideration masstransfer without calculation longitudinal diffusion (masstransfer). Solutions of
task was obtained by method trasformation of Laplace and can be disseminate on analogy tasks at part-linear
approximation unlinear isotherm of adsorption at chemical reaction of first order.
потоке на расстоянии x в момент времени  ; W скорость потока; K - коэффициент массообмена;
(1   )
K c - константа скорости реакции;  
- ко-
В настоящей работе аналитически решается
задача динамики адсорбции с использованием линейной изотермы наиболее общего вида с учетом
изменения концентраций целевого компонента за
счет химической реакции.
Как уже отмечалось [1, 2] задача динамики
адсорбции основывается на уравнении баланса массы целевого компонента (адсорбтива) для бесконечно малого элемента слоя, уравнениях кинетики адсорбции и изотермы адсорбции. Решение данной
задачи позволит при кусочно-линейной аппроксимации нелинейной изотермы адсорбции выяснить
физическую сущность изучаемого процесса и без
проблем описать любую часть изотермы сложной
формы. Кстати, следует заметить, что синтетические
цеолиты, которые часто применяются в процессах
очистки, имеют слегка выпуклую изотерму адсорбции близкую к уравнению прямой [3].
Общая система дифференциальных уравнений, описывающих динамику адсорбции, в одномерном потоке без учета продольной диффузии (с
учетом только внешнего массопереноса), с отрицательным источником массы, имеет вид
C a ( x , )
C( x , )
C( x , )


W



x
  K c C( x , ), x  0,   0 ,
(1)


C a ( x , )
 K C( x , )  C* ( C a ) , x  0,  0.

при краевых условиях
C( x ,0 )  0 , x  0 ,

эффициент;  - доля свободного сечения адсорбента (постоянная по объему порозность неподвижного
слоя); C* - концентрация целевого компонента в
потоке равновесная со средним содержанием адсорбитива C a в слое.
Первые два слагаемых в левой части уравнения (1) представляют собой скорость изменения
массы целевого компонента в зазорах между частицами и внутри частиц соответственно. Третье слагаемое соответствует приращению массы целевого
компонента за счет конвективного переноса с потоком. Слагаемое в правой части уравнения (1) учитывает изменение массы целевого компонента за счет
протекания химической реакции.
Краевые условия (3) - (5) выражают следующее.
В начальном сечении неподвижного слоя в
произвольный момент времени  концентрация
целевого компонента постоянны и равна C0 : условие (5); в начальный момент времени   0 неподвижный слой свободен от адсорбируемого вещества:
условия (3), (4).
В целях обобщения задачи принимается,
что коэффициенты A, B , K могут принимать как
положительные, так и отрицательные значения.
При K  0 имеет место сток вещества, при
K  0 источник массы за счет химической реакции.
Задача (1) - (5) решается с использованием
одностороннего преобразования Лапласа по переменной  .
Преобразования (2) с учетом выражения (6),
а также (1) с учетом (2) и (6) позволяют записать
C( x , )
C( x , )
W


x
C ( x , ) B 

  K C( x , )  a
   K C( x , ),
A
A

(2)
(3)
C a ( x ,0 )  0, x  0 ,
(4)
C( 0, ),   0.
(5)
Отличие от [1] в данной задаче состоит в
учете отрицательного источника массы за счет химической реакции первого порядка.
Уравнение изотермы адсорбции
C a ( x , )  A C * ( C a )  B .
(6)
Здесь C a ( x , ) - концентрация адсорбированного вещества в сорбенте в сечении x в момент
времени  ; C( x , ) - концентрация адсорбтива в
59
x  0 ,  0
Замена в решениях (15), (16) констант
A1 , B1 их значениями приводит к выражениям
(1а)
C a ( x , )
C ( x , ) B 

 K C( x , )  a
 ,

A
A

x  0,  0.
(2а)
К уравнениям (1а), (2а) и граничным условиям (3) - (5) применяется прямое преобразование
Лапласа [4]
F( x,s ) 
 
 
 
  s K s
K
C0

 c  x 
exp   
K W  
s
W
  W
s
 
A
 
 
K B
1


A s 2  a1s  a0

C( x , )

F ( x , s )  C( x , )e  s d ,
(7)
0

C a ( x , )

Ф( x , s )  C a ( x , )e  s d ,
(8)
 
 
 
  s K s
K
c
 x ;
 exp   

K W  
W
  W
s
 
 
A
 
0
C( x , )
s F ( x , s )  C( x ,0 ),
(9)

C a ( x , )
s Ф( x , s )  C a ( x ,0 ).
(10)

обозначает переход от ориЗдесь знак
гинала к изображению и наоборот [4].
Если умножить обе части уравнений (1а),
Ф( x , s ) 
(17)
K 2 B
1

2
A s  a1s  a0

 
 



K c   
1 
s K s


1




x 
exp


K W  
K 
W
W

s
s
 
 
A
A 
  
1


K
s( s  )
A


 
 


  s K s

K
B K 
 K C0 exp  
 c  x 
, (18)
A 
 W W s  K W  

 
 


A
 
K  A K  A K c
K K
где a1 
, a0  c .
A
A
В соотношениях (17), (18) необходимо преобразовать экспоненциальную функцию для получения табличных изображений.
Эти изображения будут использованы при
свертке оригиналов для получения окончательного
решения. Необходимо также выполнить преобразование отдельных членов, входящих в выражения
(17), (18).
 
 
 
  s K s
K
 c  x 
exp   
K W  
W
  W
s
 
 
A
 


 K 2 x 
 x 
 exp  s  exp 

 W 
 AW ( s  K ) 

A 
x

 exp  (  K  K c )  .
W

Решение в изображениях (17), (18) позволяют независимо друг от друга переходит к оригиналам. Поэтому можно осуществить последовательный переход к оригиналам сначала для функции
F ( x , s ), а затем - Ф( x , s ) или наоборот.
Преобразование (17) дает
(2а) и условие (5) на e  s и проинтегрировать каждый член уравнений по  в пределах от 0 до  , то
после незначительных преобразований с учетом (7) (10), можно получить
F ( x , s )

s F( x,s )  W
x
K
KB1
 (  K  K c ) F( x,s ) 
,
(11)
Ф( x , s ) 
A
A s
AK
KB 1
Ф( x , s ) 
F( x,s ) 
,
(12)
AS  K
As  K s
C0
.
(*)
s
Если подставить Ф( x , s ) из (12) в (11), то
совместно с граничным условием (*) получается
задача Коши для обыкновенного дифференциального уравнения первого порядка
dF ( x , s )
 A1 F ( x , s )  B1  0 ,
(13)
dx
C
(14)
F ( 0, s )  0 ,
s
K
Ks
K B
s
.
где A1 

 c , B1 
K
K
W
W
( s  )W
AW ( s  )
A
A
Здесь буква s выполняет роль параметра.
Решением (13) с условием (14) будет функция

1 
 A1 C 0

 B1  e  A1 x  .
F( x,s ) 
(15)
 B1  
A1 
 s


F ( 0, s ) 
Подстановка (15) в соотношение (14) дает
K B1
1
Ф( x , s ) 
1  e  A1 x 
K
A1
(s )
A
BK
1

K C 0 e  A1 x 

.
(16)

K
A 
s( s  ) 
A

K B
1

A s 2  a1s  a0

60
F( x,S ) 
K B
1

2
A s  a1s  a0

 
 C  x s   K 2 x 1
0
e W exp

  AW s  K
 s
 

A
Для перехода к оригиналам от изображений
в данной задаче базовыми являются следующие соответствия [4, 5]
1
 
 
  1 
 
 
 
s 2  a1s  a0
a

 1
 1 e 2 sin( a 2  ) при a 2  0 ,
 a 2

a
 1
 1
e 2 sh(  a 2  ) при a 2  0;

  a 2
s
x
C  s K B
1
 0e W 

2
s
A s  a1s  a0
e




x

s
W exp


K 2 x 1
K
 AW
s

A

K B
1
e
A s 2  a1s  a0

 
 
  1 
 
 
 
s 2  a1s  a0
a

 1
 1 e 2 sin( a 2  
 a2

 a0

 arctg a 2 ) при a  0,
2

a
 1


2

a1
 1

e 2 sh(  a 2  

a2

 a
0


 a2
 Arth
) при a 2  0;
a1



2
x
x 
s ( K  K c )
W e
W
(17а)
.


Преобразование в (18) с раскрытием скобок
и последующим умножением числителя и знаменателя у второго и третьего членов в правой части выK
ражения на ( s  ) и выделение множителя
A
1
можно получить
K
( s  )2
A
Ф( x , s ) 
K 2 B
1
1

2
K
A s  a1s  a0
(s )
A

x
2
1
s
 K B   W s


e
2
s  a1s  a0 ( s  K ) 2
 A

A


 K 2 x 1  K 3 B   x s
1

e W
 exp

2
2
K
 AW
A
s
a
s
a


1
0

s


A



x
 K 2 x 1 
 s
1
  K C0 e W 

exp
K
K
 AW
s 
( s  )2

A
A



 K 2 x 1 
1


exp
K 2
K
 AW
(s )
s 

A
A

 
  K 2 x 1
exp
  AW s  K
 
A

,


(21)
 K 2 x 
где I1  2
  - бесселева функция первого
AW




рода мнимого аргумента первого порядка;


 K 2 x 1 
1


exp
K
K 2
 AW
s 
(s )

A
A

K
  
AW
K 2 x 
e A I1  2

 ,
(22)
AW


K 2 x


x
x

 s
0 при 0    W ,
e W
(23)

s
1 при   x ;

W
K 2C0  W s 1
e
A
s
1

K 2
(s )
A


 K 2 x 1  (  K  K c ) x
W 
 e
 exp
K 
 AW
s 

A 

BK
1
.

K
A
s( s  )
A
K
 
A
(20)
 
 
  1
 
 
 
 K 2 x
I1  2

AW

AW


K 2 x
e
x

(19)
(18а)
1
s
61
K
A
e

K

A
;
(24)
x

0 при   W ,

 x 

a1 t  
W


 

x 

 1
2
e
sin  a 2  t   
 a
 W 

 2 x
W


K (  t )

K 2 x

A

 AW (   t ) e




2
 I 2 K  x (   t )  dt
1



AW




x
при   , a 2  0 ,
W

 x 

a1 t  

W

 

x 

 1
2
e
sh   a 2  t   
 a
 W 

2 x


W

K (  t )


K 2 x
A
e




A
W
(
t
)



 K 2 x



 dt
2

(
t
)
I


 1

AW




x
при   , a  0 ;
2

W
1
1
2
s  a1s  a0 s  K
A
a
t
K

1


 (  t )
 1
e 2 sin a 2 t e A
dt
 a2
0

при a 2  0 ,

at
K
 (  t )
 1   1
e 2 sh  a 2 t e A
dt

  a2 0

при a 2  0;
K
A A  A
 e
.
K K
1
(25)
K
)
A
С целью получения окончательных формул
(оригиналов) от изображений (17а), (18а) выполнялись операции свертывания оригиналов от произведения изображений с учетом таблицы оригиналов
для отдельных членов, представленных соответствиями (19) - (25).
Операция свертывания оригиналов выполнялась по известному соотношению [4]
f1 ( s ) f 2 ( s )
F1 (  )* F2 (  ) 
s( s 





 F1 ( t )F2 (   t )dt  F1(   t )F2 ( t )dt .
0
(26)
0
Кроме этого использовалась первая теорема
смещения [4]
0 при   a ,
e a s f ( s )
(27)

 F (   a ) при   a .
В результате были получены следующие
выражения
e

x
s
W
1
2
s  a1s  a0
x

0 при   W ,

x 


a1  
W




x 

 1 e
2
sin  a 2   
 a
W 


 2
x

при   , a 2  0 ,
W

x 


a1  

W



x 

 1
2
sh   a 2   
 a e
W 


2


x
при   , a 2  0;
W

x
 

 s 
e W   K 2 x 1  
 1
exp
K
s   AW
s  

A  
 
x

0 при   W ,


K

K 2 x  A t  K 2 x  ;
t dt
e
I1  2

AW


 x AWt


W
 
 
x
  K 2 x 1  
 s
1
W
  1
e
exp
s 2  a1s  a0   AW s  K  
 
A  


(28)




x
s


1

s  a1s  a0 ( s  K )2
A


 K 2 x 1 

 exp
K
 AW
s



A

e W
(29)
62
1

2
(30)
(31)
x

0 при   W ,

 x 

a1 t  
W


 
 1
2
e

 a
2 x

 a0 W





a2 
x


 sin  a 2  t  W   arctg a  



 1


2 



K
 (  t )
AW


(  t ) e A

2
 K x

 K 2 x




 I1 2 AW (   t )  dt




x

 при   , a 2  0 ,
W

x 


a1 t  

W



 1
2

e

a
2
x
 

a0 W





 a2 

x

 sh   a 2  t    Arth a  
 W


 1 


2 

K
 (  t )

AW
A



(
t
)
e


2
K
x



 K 2 x


(   t )  dt
 I1 2


AW




 при   x , a  0;
2

W
x

0 при   W ,

 x 

a1 t  
W


 
 1
2

e
 a
2 x

W




 sin  a 2  t  x  

 W 


K

 (  t )
AW

(  t ) e A

 K 2 x

 K 2 x


 I1 2
(   t )  dt


AW




x

 при   , a 2  0 ,
W

 x 

a1 t  

W

 
 1
2

e
 a
2
x


W

 sh   a  t  x  
2



 W 


K

 (  t )
 AW (   t ) e A

 K 2 x




2
 I1 2 K  x (   t )  dt



AW




x
 при   , a 2  0;
W





e

x
s
W
1
1
2
2
s  a1s  a0 
K
s  
A

 
  K 2 x 1
 exp
  AW s  K
 
A
(32)
e



 K 2 x 
1
exp 

K
 AW ( s  K ) 
( s  )2

A
A 
x

0 при   W ,

x 


K   
W


x 
 AW 
A

   e

2
W


K

x



  K 2 x
x
(   t )  при   ;
 I1  2


AW
W
 


x
s
W

 
 
  1
 
 
 
e

x
s
W
s
63
(33)

 K 2 x
1
exp 
K
 AW ( s  K
( s  )2

A
A



)

(34)



x
0 при   ,
W


K
 t
AW

A 
t
e

2
x K x
W

  K 2 x 
x
 I1  2
t  dt при   ;
AW 
W
 

 


a

x
 e
K 




 sin a 2 t dt  B1  e A .




x
При   , a 2  0
W

 
   2
K x 

e
 C0 1 
AW t
 
x
  W



K B
A a2
e
K 2B
A  a2
A a2
K
t
A I
 K  a1x  



e  A 2W 
x
a
 1
2
e



sh  a 2  ,
K 
A

a K 
 1  t
e  2 A
e



.

a
 1
e 2
x
W
 2
A

x

sin  a 2  t 
 W

  AW
(  t ) 

  K 2 x
K
x 
  
A W 
I
 K 2 x

x   C0 K 2



2





1

AW 
W 
A



K 
 
2

AW
K
x


t e A I1  2
t dt  
2
AW  

K x



K



 B 1  e A


x
При   , a 2  0
W
x
(  K  K c )
W
e


sin a 2  
C( x , ) 
F( x,s )

 K 2 x  

t dt  
1 2
AW  


 

a K 
 1  t
e  2 A



0
 
   2
K x 

e
 C0 1 
AW t
 
x
  W


K 2 x

AW (   t )

W

x

 sin  a 2  t 
 W



 K 2 x
AW 
x
 I1 2
(   t )  dt  K C 0
   
2


AW
W

K  x


A  a2
C( x , ) 
F( x,s )

W
K B
K B

a K
 1  t
e  2 A
0

  a1  K t
0
K




2 
 sh  a t dt  B1  e A




x
При   , a 2  0
W
e
K 
A
 K  a1x 

 K 2 x
3




K
B
(   t )  dt 
e  a 2W  
 I1 2
2


AW
A a2



C a ( x , ) 
A a2




a 2  AW
x

 sin  a 2  t    arctg
(  t ) 

2
a
 W

 1 K x
2 


Ф( x , s )
K 2B

 K  a1 x  
 a1 K 
 2



 K B   a 2W   2  A t
e
e

 sin a 2 t dt  
 A a2
x

a0
W

1
K B  2
e
sin a 2 t ,
A a2
K   a1  K t
K 2B  A
C a ( x , ) 
e
e  2 A 
A a2
Ф( x , s )
C a ( x , ) 

C( x , ) 
C( x , ) 
A a2
Ф( x , s )
F( x,s )
F( x,s )
K B


x

x   (  K  Kc )W

sin  a 2    e
;
W  




(35)
Наконец, можно, в соответствии с выполненными преобразованиями, записать окончательное решение исходной задачи
x
При   , a 2  0
W

x 

a1  
W


2
e


2
  K  x
(   t )  dt 
  I1 2

AW
 


K B
A  a2
K B
A  a2
K
t
A I
 K  a1x  



e  A 2W 
x
W
64


.

a
 1
e 2


sh  a 2  

 K 2 x  

t dt  
1 2
AW  


 

a K 
 1  t
e  2 A
K 2 x

AW (   t )


x    K 2 x

 sh   a 2  t   I1 2
(   t )  dt 

AW
 W  



x



a1  
x

 W


K B
x   (  K  Kc )W

2
;
e
sh   a 2    e

W  
A  a2




Ф( x , s )
C a ( x , ) 
K 2B
A  a2
e

K 
A

a K 
 1  t
e  2 A


0

 K  a1 x  
a K 
 2



 1  t
 K B
a 2W 

 sh  a 2 t dt  
e
e  2 A 
a
A  2
x

a
W
0






x

 sh   a 2  t 
W



 a2

  Arth
a

 1
2

 AW
(  t ) 

2
 K x

Литература
 K  a1x 
 K 2 x

K 3 B   a  2W 


 I1 2
(   t ) dt 
e



AW
A2  a2


  a1  K t
x
 e
 2
A
K
x 
  
A W 
I
 K 2 x
x   C 0 K 2


2



 

1

AW 
W 
A



K 



2

AW
K x  
t e A I1  2
t dt  

2
AW  

x K  x



W
x
K 

(  K  K c )



W
e
 B 1  e A .




При переходе в пространство оригиналов
для функции F ( x , s ) - использовались соответствия
(19), (23), (28) - (30), а для функции Ф( x , s ) - соответствия (31) - (35). При этом были преобразованы
произведения экспоненциальных функций, имеющиеся в представленных оригиналах.
Следует отметить, что в решениях для
x

определенные интегралы могут быть вычисW
лены в конечном виде.
e

1. С.В. Анаников, Вестн. Казан. технол. ун-та, 15, 8, 247253 (2012).
2. С.В. Анаников, Вестн. Казан. технол. ун-та, 15, 10, 247250 (2012).
3. Л.М. Никитина Термодинамические параметры и коэффициенты массопереноса во влажных материалах. Энергия, Москва, 1968. 500 с.
4. Г. Дёч, Руководство к практическому применению преобразования Лапласа. ГИФМЛ, Москва, 1958. 208 с.
5. В.А. Диткин, Справочник по операционному исчислению. ГИТТЛ, Москва - Ленинград, 1951. 256 с.

x   AW

sh   a 2  t  
(  t ) 
 W   K 2 x

W

 K 2 x
AW 
x
 I1 2
(   t )  dt  K C 0
   
2


AW
W


K x




© О. Р. Каратаев - к.т.н., доцент каф. машиноведения КНИТУ, oskar_karataev@mail.ru; С. В. Анаников – д-р техн. наук,
проф. каф. химической кибернетики КНИТУ; ananikovsv@rambler.ru.
65
Документ
Категория
Без категории
Просмотров
4
Размер файла
224 Кб
Теги
динамика, очистки, бассейном, адсорбция, реакций, стоков, плавательных, химические, учетом
1/--страниц
Пожаловаться на содержимое документа