TÓPICO 2: EQUAÇÕES DA CONTINUIDADE EM

TRANSFERÊNCIA DE MASSA

I. EQUAÇÃO DA CONTINUIDADE MÁSSICA DE UM SOLUTO A;

II. EQUAÇÕES DA CONTINUIDADE DE UM SOLUTO A EM TERMOS DA
LEI ORDINÁRIA DA DIFUSÃO;
III. EQUAÇÃO DA CONTINUIDADE MOLAR DE UM SOLUTO A;
IV. POSSÍVEIS SIMPLIFICAÇÕES DA EQUAÇÃO DA CONTINUIDADE;
V. CONDIÇÕES DE CONTORNO;
VI. APLICAÇÕES.

BIBLIOGRAFIA:
CREMASCO, M.A. Fundamentos de Transferência de Massa. Ed. Unicamp.
 I. EQUAÇÃO DA CONTINUIDADE MÁSSICA EM
TRANSFERÊNCIA DE MASSA:

 ANÁLISE PONTUAL DO FENÔMENO DE T.M. POR INTERMÉDIO DO

CONHECIMENTO DA DISTRIBUIÇÃO DE CONCENTRAÇÃO DE UM
DETERMINADO SOLUTO NO TEMPO E NO ESPAÇO;

 A PARTIR DO BALANÇO DE MASSA, ONDE A MATÉRIA FLUI ATRAVÉS DAS

FRONTEIRAS DE UM VOLUME DE CONTROLE, ELEITO NO MEIO CONTÍNUO:
 taxa de  taxa de  taxa de  taxa de 
       
entrada 
  saída 
  produção 
 acúmulo 
de massa de massa de massa  de massa
       

n A x  x REAÇÃO QUÍMICA
z

z

y
n A x  x  x
y

x
Em x: taxa de 
n A x  x    n A x  x yz
z entrada 
taxa de
   n A x  x  x yz
saída 
taxa de 
 A  xyz
,,,
  r
 produção
n A x  x  x
taxa de   A
y   xyz
acúmulo  t
x

n A x  x  x  n A x  x  n A x  x  x DERIVADA PARCIAL
x
taxa de  taxa de  
 

  n A  x  x yz  n A  x  X  n A x  x x yz 
entrada  entrada   X 

 n A x  x xyz
X
Em Y: taxa de 
z
   n A  y  y xz
entrada 
taxa de
nA y y    n A  y  y  y xz
n A  y  y  y saída 
taxa de 
y
 A  xyz
,,,
  r
 produção

y taxa de   A
  xyz
acúmulo  t
x

taxa de  taxa de 
  

y

n A  y  y xz 
entrada  entrada 
Em z: n A z  z  z taxa de 
z
   n A z  z xy
entrada 
taxa de
z    n A z  z  z xy
saída 
taxa de 
 A  xyz
,,,
  r
n A z  z
 produção

y taxa de   A
  xyz
acúmulo  t
x

taxa de  taxa de  
     n A z  z xy
entrada  entrada  z
 taxa de  taxa de  taxa de  taxa de 
       
entrada 
  saída 
  produção 
 acúmulo 
de massa de massa de massa  de massa
       

 A 
xyz  rA,,, xyz  n A x  x x yz 
t x

y
  
n A  y  y y xz  n A  z  z z xy
z
 xyz 
 A  n A, x n A, y n A, z 
    
 A r ,,,

t   x  y  z 

 A  A
    n A  rA 
,,,
   n A  rA,,,
t t
C A C A
   N A  R A 
,,,
   N A  R A,,,
t t
 EQUAÇÕES DA CONTINUIDADE EM COORDENADAS RETANGULARES:

dz

dy

 A  n A, x  n A, y  n A, z 
   
 A r ,,,

t  x y z 

C A  N A, x   N A, y  N A, z 
     R ,,,

t  x y z 
A
 EQUAÇÕES DA CONTINUIDADE EM COORDENADAS CILÍNDRICAS:

rd

 A  1 rnA,r  1 n A,  n A, z 

    ,,,

z 
r
t  r r rsen 
A

C A  1 rN A,r  1 N A,  N A, z 

    ,,,

z 
R
t  r r rsen 
A
 EQUAÇÕES DA CONTINUIDADE EM COORDENADAS ESFÉRICAS:

 A  1
 2
 
r 2

n A, r

1 senn A, 

1 n A, 
 ,,,
 A r
t  r r rsen  rsen  

C A  1
 2
 r2
N A, r

 1 senN A, 

1 N A, 
 ,, ,
 R
t  r r rsen  rsen  
A
 II. EQUAÇÃO DA CONTINUIDADE DE UM SOLUTO A EM
TERMOS DA LEI ORDINÁRIA DA DIFUSÃO:

 A
   n A  rA
,,,

t (1)

MAS: nA  j A  Av (2)

(2) EM (1):
 A
t
 
   j A   A v  rA
,,,
  A
t
 
    A v    j A  rA
,,,

FLUXO DIFUSIVO: j A   D A, B   A

 A
t
  
    A v    D AB  A  rA
,, ,

CONTRIBUIÇÃO CONTRIBUIÇÃO
ACÚMULO GERAÇÃO
CONVECTIVA DIFUSIVA
TEMOS UMA EXPRESSÃO SEMELHANTE PARA O COMPONENTE B:

 B
t
  
    B v    D AB  B  rB
,, ,

   nB
 B
    n B  rB
,, ,

t
PARA UMA MISTURA BINÁRIA:

 A  B
    n A    n B  rA  rB
,,, , ,,

t t
 PELA LEI DA
CONSERVAÇÃO
DA MASSA

 A  B
    n A    n B  rA  rB
,,, , ,,

t t

t
 
 A   B     n A  n B  0

 A  B   
n  0
nA  nB  n t
MAS: n  v

t
 
   v  0

DA ANÁLISE VETORIAL:

   
   v  v       v



t
 
 v       v  0

SE   cte
D
Dt
 
  v  0 v
III. EQUAÇÃO DA CONTINUIDADE MOLAR EM TRANSFERÊNCIA
DE MASSA:

C A
   N A  RA
,, ,

t
C B
   N B  RB,,,
t
PARA UMA MISTURA BINÁRIA:

t
 
C A  C B     N A  N B  R A,,,  RB,,,

C
    C v  R A  RB
,,, ,, ,

t
ABRINDO O DIVERGENTE:

C
t
 
 vC  C   v  R A,,,  RB,,,

DC
Dt
 
 C   v  R A  RB
,,, , ,,

EM TERMOS DA LEI ORDINÁRIA DA DIFUSÃO:

C A
t
  
   C A v    D AB C A  R A
,,,

CONTRIBUIÇÃO CONTRIBUIÇÃO
ACÚMULO GERAÇÃO
CONVECTIVA DIFUSIVA
 IV. POSSÍVEIS SIMPLIFICAÇÕES:
 QUANTIDADE DO COMPONENTE A MUITO BAIXA, DIFUNDINDO-SE NA
ESPÉCIE B, ESTAGNADA:

   
0
C A
   C A v    D AB C A  R A
,,,

t
 EM REGIME ESTACIONÁRIO:

   
0
C A
   C A v    D AB C A  R A,,,
t
 EM UM PROCESSO ONDE NÃO OCORRA REAÇÃO QUÍMICA:

C A
t
   
   C A v    D AB C A  R A,,,
0
EXEMPLO 1: A EFICÁCIA DE PRODUTOS FARMACÊUTICOS É REDUZIDA PELA
EXPOSIÇÃO PROLONGADA A ALTAS TEMPERATURAS, À LUZ E À UMIDADE. PARA
PRODUTOS CONSUMIDOS QUE SÃO SENSÍVEIS AO VAPOR D`ÁGUA E ENCONTRAM-SE
NA FORMA DE COMPRIMIDOS OU CÁPSULAS E SÃO GUARDADOS EM AMBIENTES
ÚMIDOS, EMBALAGENS BILSTER SÃO USADAS PARA LIMITAR A EXPOSIÇÃO DIRETA
DO MEDICAMENTO À CONDIÇÕES DE UMIDADE ATÉ O MOMENTO IMEDIATAMENTE
ANTERIOR A SEU USO.
CONSIDERE COMPRIMIDOS QUE ESTÃO CONTIDOS EM UMA EMBALAGEM
BLISTER COMPOSTA POR UMA FOLHA DE COBERTURA PLANA E UMA SEGUNDA FOLHA
MOLDADA, QUE POSSUI LOCAIS QUE ABRIGAM CADA COMPRIMIDO. A FOLHA
MOLDADA TEM A ESPESSURA L=50m E É FABRICADA COM UM MATERIA
POLIMÉRICO. CADA LOCAL DE COMPRIMIDO POSSUI DIÂMETRO D=5 mm E
PROFUNDIDADE h=3mm. A FOLHA DE COBERTURA É FEITA DE ALUMÍNIO. O
COEFICIENTE DE DIFUSÃO BINÁRIA É DAB=6.10-14 m2/s, ENQUANTO PODE-SE SUPOR
QUE O ALUMÍNIO SEJA IMPERMEÁVEL EM RELAÇÃO AO VAPOR D`ÁGUA. PARA
CONCENTRAÇÕES MOLARES DO VAPOR D`ÁGUA NAS SUPERFÍCIES EXTERNA E
INTERNA DE CA,S1=0,0045Kmol/m3 E CA,S2=0,0005Kmol/m3, RESPECTIVAMENTE,
DETERMINE A TAXA NA QUAL O VAPOR D`ÁGUA É TRANSFERIDO ATRAVÉS DA
PAREDE DO COMPARTIMENTO PARFA O COMPRIMIDO.
CONSIDERAÇÕES: - A FOLHA DO POLÍMERO É MUITO FINA, EM RELAÇÃO ÀS
DIMENSÕES DO COMPARTIMENTO DO COMPRIMIDO
(UNIDIRECIONAL);

- EM REGIME PERMANENTE;

- T.M. EM MEIO ESTACIONÁRIO.

 DADOS:

CV ,Z 1  0,0045 Kmol m 3 
CV ,Z 2  0,0005Kmol m 3

z
L  5  10 6 m
D  0,005m
h  0,003m
 
DVAPOR  AR  6 1014 m2 s  cte

- T.M. UNIDIMENSIONAL (Z);

- AR ESTAGNADO:

  VAPOR v  0
- EQUAÇÃO DA CONTINUIDADE EM COORDENADAS CILÍNDRICAS: REGIME
ESTACIONÁRIO

1   CV  1   CV    CV  CV

 DVA r  2  DVA    DVA   RV 
,,,

r r  r  r     z  z  t
SEM
T.M. UNIDIMENSIONAL (Z) REAÇÃO QUÍMICA
 V. CONDIÇÕES DE CONTORNO:

V.1 CONDIÇÃO INICIAL:

t  0  C A  C A0
t  0   A   A0
t  0  X A  X A0
t  0  W A  W A0

V.2 CONDIÇÕES DE CONTORNO: EM POSIÇÕES ESPECÍFICAS NO

VOLUME DE CONTROLE OU NAS FRONTEIRAS DESSE VOLUME.
BASICAMENTE AS CONDIÇÕES DE FRONTEIRA SÃO:

V.2.1 CONCENTRAÇÃO OU FRAÇÃO DO SOLUTO ESPECIFICADA NUMA

DETERMINADA FASE;

V.2.2 CONDIÇÕES DE FLUXO;

V.2.3 REAÇÃO QUÍMICA CONHECIDA.

V.2.1 CONCENTRAÇÃO OU FRAÇÃO DO SOLUTO ESPECIFICADA
NUMA DETERMINADA FASE:

VOLUME DE CONTROLE

FLUXO FLUXO

 A   AS
C A  C AS
FRONTEIRA S
FRONTEIRA S
INICIAL
FINAL
WA  WAS
X A  X AS
YA  YAS
- NA FASE GASOSA IDEAL:

YA  YAS LEI DE DALTON: PA, S  y A, S P

- NA FASE LÍQUIDA, PARA UMA SOLUÇÃO IDEAL:

X A  X AS LEI DE RAOULT: PA, S  x A,S P

VAP

EQ. DE ANTOINE:

P VAP  mmHg F
T K ln P VAP  E 
T  G 

ESPÉCIES E F G
ÁGUA 18,3036 3618,44 -46,13
BENZENO 15,9008 2788,51 -64,38
TOLUENO 16,0137 3096,52 -53,67
METANOL 18,5875 3626,55 -34,29
ETANOL 19,9119 3803,98 -41,68
- NA HIPÓTESE DO EQUILÍBRIO TERMODINÂMICO NA FRONTEIRA “S” OU
INTERFACE ENTRE AS FASES LÍQUIDA E GASOSA, CONSIDERANDO-AS
IDEAIS:

x A,S  PVAP  y A,S  P EQUAÇÃO DE RAOULT-DALTON

- SE A FASE LÍQUIDA FOR CONSTITUÍDA SOMENTE DA ESPÉCIE “A”:

P VAP
y A, S 
P

- SE A SOLUÇÃO FOR DILUÍDA:

PA,S  x A,S H
CONSTSANTE
DE
HENRY
TABELA DOS VALORES DE H PARA GASES EM ÁGUA: (Hx10-4), (PRESSÃO EM
atm):

T(C) H2 N2 O2 CO CO2
0 5,79 5,29 2,55 3,52 0,0728
10 6,36 6,68 3,27 4,42 0,1040
20 6,83 8,04 4,01 5,36 0,1420
30 7,29 9,24 4,75 6,20 0,1860

- NA CONDIÇÃO DE EQUILÍBRIO TERMODINÂMICO LÍQUIDO-VAPOR NA

FRONTEIRA OU NA INTERFACE “S” E ADMITINDO FASES IDEAIS:
- NA CONDIÇÃO DE EQUILÍBRIO TERMODINÂMICO LÍQUIDO-VAPOR NA
FRONTEIRA OU NA INTERFACE “S” E ADMITINDO FASES IDEAIS:

y A,S  mx A,S OU PA,S  m*C A,S

H
m
P
ONDE:
H
m* 
C
TAIS RELAÇÕES DE EQUILÍBRIO SÃO UTILIZADAS EM FENÔMENOS DE
DESSORÇÃO E ABSORÇÃO; NESTES FENÔMENOS O SOLUTO A ESTÁ
CONTIDO NAS FASES GASOSA E LÍQUIDA. NA VENTURA DE ESTAR
DISTRIBUÍDO E DILUÍDO NAS FASES SÓLIDO-LÍQUIDO,A RELAÇÃO DE
EQUILÍBRIO É ESCRITA ANALOGAMENTE À LEI DE HENRY SEGUNDO:

C A,1S  K P  C A, 2 S
COEFICIENTE DE DISTRIBUIÇÃO
C A,1S  K P  C A, 2 S

ESSA RELAÇÃO É ÚTIL NAS OPERAÇÕES QUE ENVOLVEM AS

FASES SÓLIDO / FLUIDO QUANDO SE DESEJA ESPECIFICAR
UMA RELAÇÃO DE EQUILÍBRIO ENTRE A CONCENTRAÇÃO DO
SOLUTO PRESENTE NO INTERIOR DO SÓLIDO E AQUELA NO
SEIO DA FASE FLUIDA, OU SEJA:

C *
 K C
A,1 P A, 2 
EXEMPLO 2: CALCULAR OS DADOS DE EQUILÍBRIO SOB A FORMA y
VERSUS x PARA O SISTEMA m-XILENO/o-XILENO À PRESSÃO TOTAL DE 0,4
atm (XA e YA SÃO AS FRAÇÕES MOLARES DE m-XILENO NAS FASES LÍQUIDA
E VAPOR, RESPECTIVAMENTE). SÃO DADAS AS PRESSÕES DE VAPOR DOS
COMPONENTES EM atm NA TABELA A SEGUIR:
EXEMPLO 3: ESCREVA A EQUAÇÃO DA CONTINUIDADE MOLAR DE A NA
FORMA JÁ SIMPLIFICADA E AS CONDIÇÕES DE CONTORNO PARA A
SEGUINTE SITUAÇÃO:
UM CERTO GÁS DIFUNDE POR UMA PELÍCULA ESTAGNADA DE AR DE
0,5 cm DE PROFUNDIDADE EM UM CAPILAR QUE CONTÉM CERTO ÁCIDO. AO
ATINGI-LO, O GÁS É ABSORVIDO INSTANTANEAMENTE. A CONCENTRAÇÃO
DO GÁS NA BOCA DO RECIPIENTE É 0,25% EM MOLS.

y A1  0,0025 z=0
N A, Z

z
N B,Z

y A2  0 z=0,5 cm

H 2 SO4 l
V.2.2 CONDIÇÕES DE FLUXO:

EM SITUAÇÕES ONDE O SOLUTO FLUI DE UMA FASE A OUTRA,

PRESSUPONDO QUE A INTERFACE NÃO OFEREÇA RESISTÊNCIA AO
TRANSPROTE DO SOLUTO  CONTINUIDADE DE FLUXO

NA FASE 1: O FLUXO DO SOLUTO SERÁ DEVIDO À DIFUSÃO;

NA FASE 2: O FLUXO DO SOLUTO SERÁ DEVIDO À CONVECÇÃO MÁSSICA.

NA FASE 1: O FLUXO DO SOLUTO SERÁ DEVIDO À DIFUSÃO;
dC A1
N A, Z Z S
  Def1
dz
d A1
n A, Z Z S
  Def1
dz

NA FASE 2: O FLUXO DO SOLUTO SERÁ DEVIDO À CONVECÇÃO MÁSSICA.


N A,Z  k m2 C A2 S  C A2  
n A,Z  k m2  A2 S   A2  
ADMITINDO QUE A INTERFACE NÃO OFEREÇA RESISTÊNCIA À MOBILIDADE
DO SOLUTO  CONTINUIDADE DO FLUXO DE MATÉRIA NA FRONTEIRA
CONSIDERADA  IGUALDADE DOS FLUXOS:

 Def1
dC A1
dz

 k m2  C A2 S  C A2  
Z S

A CONCENTRAÇÃO DE A NA FASE 2 E CONTIDA NA INTERFACE ESTÁ EM

EQUILÍBRIO TERMODINÂMICO COM A CONCENTRAÇÃO DE A NA FASE 1 E
CONTIDA NA INTERFACE “S”, POR INTERMÉDIO DA SEGUINTE RELAÇÃO:

C A,1S  K P C A, 2 S

  Def1
dC A1
dz
 k m2 
 C A1S

 KP


 C A2    Def1

dC A1
dz

k m2
KP

 C A1S  K P C A2  
Z S   Z S
MAS: K P C A, 2  C A* ,1
*S
dC A1
dz

k m2
Def1 K P

 C  C A1S
*
A1  sdC A1
dz

sk m2
Def1 K P

 C A*1  C A1S 
Z S Z S

BiM

EXEMPLO 4: UM PELLET CILINDRICO GELATINOSO DE 2 mm DE DIÂMETRO

E 20 mm DE COMPRIMENTO, CONTENDO INICIALMENTE 50 (g DE
SACAROSE)/(l DE GEL) É POSTO SUBITAMENTE EM UM TANQUE NO QUAL HÁ
UMA SOLUÇÃO AQUOSA QUE APRESENTA 8 (g DE ACAROSE)/(L DE
SOLUÇÃO). NO TANQUE HÁ ALIMENTAÇÃO E RETIRADA LENTA E CONTÍNUA
DA SOLUÇÃO. ADMITINDO QUE NÃO HOUVE TEMPO SUFICIENTE PARA O
ESTABELECIMENTO DA DIFUSÃO DA SACAROSE NO GEL EM REGIME
PERMANENTE, ASSIM COMO SUPONDO A INFLUÊCIA DA CONVECÇÃO
MÁSSICA EXTERNA NA DIFUSÃO DA SACAROSE NO INTERIOR DO GEL,
PEDE-SE A EQUAÇÃO DA CONTINUIDADE MÁSSICA DA SACAROSE JÁ
SIMPLIFICADA QUE DESCREVE A DISTRIBUIÇÃO DE CONCENTRAÇÃO NO
GEL.
V.2.3 REAÇÃO QUÍMICA CONHECIDA:

HOMOGÊNEA;
REAÇÕES QUÍMICAS:
HETEROGÊNEA.

REAÇÃO HOMOGÊNEA: OCORRE EM TODOS OS PONTOS DO VOLUME DE

CONTOLE; NESTE CASO, A DESCRIÇÃO DA REAÇÃO QUÍMICA APARECE
DIRETAMENTE COMO TERMO DA EQUAÇÃO DA CONTINUIDADE (MOLAR OU
MÁSSICA);

C A
   N A  R A,,,
t

 A
   n A  rA,,,
t
REAÇÃO HETEROGÊNEA: OCORRE NA SUPERFÍCIE DE UMA PARTÍCULA, A
QUAL É CONSIDERADA COMO UMA FORNTEIRA À REGIÃO ONDE HÁ O
TRANSPORTE DO SOLUTO; NESTE CASO, O TERMO REACIONAL APARECERÁ
COMO CONDIÇÃO DE CONTORNO E NÃO NA EQUAÇÃO DIFERENCIAL QUE
REGE O PROCESSO DE TRANSFERÊNCIA DE MASSA.

R A,,  N A,Z Z 
 k S C A  k S y A C

PARA z=:
N A, Z Z 
yA  
kS C
A TAXA DE PRODUÇÃO (OU DESAPARECIMENTO) DE UMA DETERMINADA
ESPÉCIE QUÍMICA, PRESENTE NA SOLUÇÃO, ESTÁ ASSOCIADA COM A
REAÇÃO QUE PODE OCORRER DURANTE O TRANSPORTE DO SOLUTO. POR
EXEMPLO, SE A ESPÉCIE “A” FOR GERADA POR UMA REAÇÃO DE PRIMEIRA
ORDEM E ESTIVER ORIENTADA NO SENTIDO DO FLUXO DA MATÉRIA, SEU
FLUXO DE PRODUÇÃO SERÁ:

R A,,  N A,Z Z S
 K S C A2 S

CONSIDERANDO QUE A REAÇÃO QUÍMICA OCORRA NA SUPERFÍCIE OU EM

ÁREAS RESTRITAS DE UM SÓLIDO POROSO E QUE, DEVIDO À
CONTINUIDADE DO FLUXO DE MATÉRIA, O SOLUTO DIFUNDA PELA MATRIZ,
NA FRONTEIRA S, TEREMOS A IGUALDADE ENTRE OS FLUXOS, DA SEGUNTE
MANEIRA:

dC A1 dC A1 kS
 Def  k S C A2 S  dz  CA
Def K P 1S
dz z S
z S
EXEMPLO 5: A QUEIMA DO GRAFITE (CARBONO PURO) NO AR PODE SER
DESCRITA POR MEIO DAS SEGUINTES ETAPAS:
1. O OXIGÊNIO DIFUNDE ATRAVÉS DE UMA PELÍCULA DE AR QUE ENVOLVE
A PARTÍCULA DE GRAFITE ATÉ ATINGIR A SUPERFÍCIE DO SÓLIDO;
2. HÁ O CONTATO DO O2 COM A SUPERFÍCIE DO GRAFITE,
PORPORCIINANDO A SEGUINTE EQUAÇÃO:

C s   O2 g   N 2 g   CO2 g   N 2 g 
QUE É DESCRITA PELA REAÇÃO IRREVERSÍVEL DE PRIMEIRA ORDEM:

RO,,2  k S CyO2
3. DIFUSÃO DO CO2, COMO PRODUTO DA REAÇÃO, DA SUPERFÍCIE DO
GRAFITE PARA A PELÍCULA DE AR.

ADMITINDO QUE A PARTÍCULA DE GRAFITE TENHA A FORMA ESFÉRICA,

DESEJA-SE OBTER A EQUAÇÃO DA CONTINUIDADE MOLAR QUE
DESCREVE A DISTRIBUIÇÃO DA FRAÇÃO MOLAR DO O2 NO AR, ASSIM
COMO AS CONDIÇÕES DE CONTORNO.