Université des Sciences et de la Technologie

Houari Boumediene

Faculté de Génie Mécanique et de Génie des Procédés

Filière Génie des Procédés (GP)
1ère année Master :
Pétrochimie

TD

Phénomènes de Transfert II

2020 / 2021
DMS 4: Diffusion avec réaction chimique homogène
Un soluté A diffuse à travers la surface d’un liquide, au niveau de laquelle
sa concentration molaire est cA0, vers le liquide B (inerte) avec lequel il réagit
dans un réacteur.

La diffusion est supposée unidirectionnelle et unidimensionnelle et la

réaction est du premier ordre par rapport au soluté. Ce dernier est supposé être très
dilué dans le liquide.

Dans le réacteur fonctionnant en régime stationnaire, la densité de flux

molaire de A à une profondeur L chute de moitié par rapport à celle relative à
l’interface.

1. Donner l’expression de la concentration du soluté à une profondeur z

de l’interface en fonction de la diffusivité DA,B et de la constante de vitesse k.

2. Exprimer la densité du flux molaire à l’interface en fonction des

données du problème.
Hypothèses :

 Ect. unidirectionnel et unidimensionnel.

 Régime stationnaire.
 A est très dilué dans B (yA → 0).
 B  inerte
 Réaction chimique irréversible d’ordre 1

CA(z) = ?
 
N A   c DAB grad y A  y A  N j
c A
N A   DAB
j
1
N A   c DAB
y A
 y A N A  N B 
z
z

Or d’après l’équation du bilan

c A N Ax N Ay N Az
    RA
t x y z
 N A
 kC A 2
z
Disparition de A

 c A 
De 1 et 2   D AB   kCA
z  z 
 2c A
 D AB  kCA  0
z 2
 2c A k
 CA  0 Equation différentielle
z 2
DAB du 2nd ordre

2 2 2 k
L’équation caractéristique: r a 0 avec a 
D AB

La solution aura la forme: C A   e az   e az

 C A ( z )  A Ch( az )  B Sh( az )
Calculons les constantes A et B:
CL1: z  0
 C A ( z  0 )  C A0
1
CL2: zL
 N A zL
 NA z 0
2

CL1: C A ( 0 )  A Ch( 0 )  B Sh( 0 )  C A0 A  C A0

1 0
dC A 1 dC A
CL2: 
dz zL 2 dz z 0

a A Sh( aL )  a B Ch( aL )  a A Sh( 0 )  a B Ch( 0 )

1
2
0 1
2C A0 Sh( aL )
B
1  2Ch( aL )
 2C A0 Sh( aL )
C A ( z )  C A0 Ch( az )  Sh( az )
1  2Ch( aL )

C A( z ) 
C A0
Ch( az )  2Ch( aL ) Ch( az )  Sh( aL ) Sh( az )
1  2Ch( aL )

C A( z ) 
C A0
Ch( az )  2Ch( aL ) Ch( az )  Sh( aL ) Sh( az )
1  2Ch( aL )
Ch( aL  az )
C A( z ) 
C A0
Ch( az )  2Ch( aL  az )
1  2Ch( aL )

Posons: kL2
b1  a L 
D AB
C A0    z  
C A( z )  Ch( az )  2Ch b1 1   
1  2Ch( b1 )    L  
La densité de flux molaire à l’interface (z = 0)
c
N A ( z )   DAB A
z
  
N A ( z )   DAB 
C A0
Ch( az )  2Ch( aL  az )
z  1  2Ch( aL ) 

N A ( z )   DAB
C A0
aSh( az )  2aSh aL  az 
1  2Ch( b1 )

N A ( z )  DAB
aC A0
Sh( az )  2Sh aL  az 
1  2Ch( aL )
 Sh( az )  2Sh aL  az 
N A ( z )  C A0 kDAB  
 1  2Ch( aL ) 

 2Sh aL  
z=0 N A ( 0 )  C A0 kDAB  
1  2Ch( aL ) 