CAPÍTULO II

MÉTODO DE SEPARACIÓN DE VARIABLES O MÉTODO DE FOURIER

El método de separación de variables, MSV, es una herramienta poderosa para encontrar la

solución de problemas de la física-matemática, en su proceso se obtienen varias ecuaciones
diferenciales ordinarias, EDO, a partir de la ecuación diferencial parcial, EDP, con lo cual, se
simplifica la solución del problema. Más adelante, se presenta un ejemplo que describe de forma
explícita este método.
Etapas generales del método de separación de variables
Requisitos para la aplicación del MSV

1. Que el problema descrito por la EDP sea homogeneo

2. Que se cuente con al menos una segunda derivada dentro del problema y para ella exista un
par de condiciones homogeneas evaluadas en distinto punto
3. Contar con al menos UNA condicion No-homogenea en el problema

2.1 Transferencia de Calor Unidimensional en Estado No Estacionario (Problema que genera

una serie de Fourier de senos).

Suponer una barra idealmente aislada para llevar a cabo transferencia de calor solo en
dirección x, dentro de la cual existe una distribución de temperatura que se puede expresar como
una función de la posición. En el instante t = 0, los extremos x = 0 y x = l se colocan a temperatura
cero y son mantenidos a esa temperatura, con lo cual se tendrá una disipación del calor hacia los
extremos. Se desea evaluar la distribución de temperatura como una función de la posición y del
tiempo, T(x,t).

Este proceso lo gobierna la ecuación de calor

 T  x, t   2T  x, t 
 a2 , 0  x  L, t  0 . (2.1.1)H
t  x2
Con condición inicial
T  x, 0   f  x  , 0  x  L , (2.1.2)NH

y condiciones frontera
T 0, t   0, t  0 , (2.1.3)H

T  L, t   0, t  0 . (2.1.4)H

El problema anterior consta de una EDP de segundo orden homogénea (2.1.1), una condición inicial,
CI, no homogénea (2.1.2) y dos condiciones frontera, CF, homogéneas, (2.1.3 y 4). Estas
características definen un problema homogéneo. Hay que notar que la EDP contiene igual número
de derivadas que de condiciones, por lo cual se tiene un problema totalmente determinado. El
método de separación de variables (MSV) es específico para este tipo de problemas y no se puede
aplicar a problemas que No cumplan con los requisitos anteriores.
La suposición fundamental del MSV es que la variable dependiente se considera igual a un producto
de funciones, en donde cada una de estas, depende solo de una de las variables independiente. En
este caso
T  x ,t   X  x  t  , (2.1.5)
esta solución propuesta debe satisfacer la EDP (2.1.1) y las CF (2.1.3) y (2.1.4).
Por lo tanto, al sustituir la ecuación (2.1.5) en (2.1.1) se obtiene
( ) ( ) ( ) ( ) ( ) ( )
=𝑎  𝑋(𝑥) = 𝑎 𝜃(𝑡) 

la cual al ser dividida  a 2 t X  x  es

( ) ( ) ( ) ( ) ( ) ( )
=  =
( ) ( ) ( ) ( ) ( ) ( )

1 1 d t  1 d 2 X x 
  2
a  t  dt
2
X x  dx 2

y en donde la igualdad a  2 (constante de separación por determinar) se justifica por la

independencia de cada uno de los términos en la ecuación, ya que el lado izquierdo de la ecuación
es independiente de x y el derecho es independiente de t.
La ecuación anterior proporciona dos EDO’s,
 d t 
 2 a 2 t  (2.1.6)
dt

∫ = −𝜆 𝑎 ∫ 𝑑𝑡  𝐿𝑛(𝜃) = −𝜆 𝑎 𝑡 + 𝐿𝑛𝐶  solución (2.1.8)

d 2 X x 
2
 2 X  x  . (2.1.7)
dx
Su ec’n auxiliar es 𝑚 + 𝜆 = 0  𝑚 = ±𝜆𝑖  solución (2.1.9)
Las soluciones de las ecuaciones (2.1.6) y (2.1.7) son:

t   C1e   a t
2 2
(2.1.8)

X  x   C 2 cos( x )  C 3 senx  (2.1.9)

Para encontrar el valor de las constantes de integración y de separación, es necesario obtener
condiciones para cada una de las funciones independientes. Esto se hace mediante la sustitución de
(2.1.5) en (2.1.3) y (2.1.4),
T 0, t    t X 0  0  X 0  0 (2.1.10)

T  L, t     t  X  L   0  X  L   0 (2.1.11)

La obtención de (2.1.10) y (2.1.11) se justifica a causa de que  t  debe ser diferente de cero, ya
que forma parte, en producto, de la solución propuesta.

Aplicando la condición (2.1.10) en (2.1.9) se obtiene

1 0

X 0   0  C 2 cos0   C 3 sen0 

de donde se determina que C 2  0 y por lo tanto de la solución (2.1.9) queda solo

X  x   C 3 senx  . (2.1.12)
Ahora, aplicando la condición (2.1.11) en la ecuación (2.1.12) se tiene
X  L   0  C3 sen   L  (2.1.13)

Como C 3 no puede ser igual cero debido a que X(x) sería igual cero y en consecuencia la solución
buscada también seria igual a cero, entonces se tiene que
sen   L   0 (2.1.14)

La ecuación anterior es llamada ecuación característica y se cumple cuando el producto λL es

múltiplo de  , o sea
 L  n , n  1, 2,3,......
por lo tanto, se pueden determinar los valores de λ,
n
 , n  1, 2,3,...... (2.1.15)
L
Los valores de  son llamados valores característicos y siempre se obtienen como las raíces
positivas de la ecuación característica. Con lo anterior las soluciones (2.1.8) y (2.1.12) para θ(t) y
X(x) son:
2
 n  2
  a t
  t   C1e  L 
n  1, 2, 3,..... , (2.1.16)

 n 
X  x   C3 sen   x, n  1, 2,3,..... , (2.1.17)
 L 
y de acuerdo a la solución propuesta (2.1.5) se tiene
2
 n 
 n    a2 t
T  x, t   C1C3 sen  xe  L  , n  1, 2,3,.... , (2.1.18)
 L 
El principio de superposición establece que la solución de una ecuación diferencial es la
combinación lineal de soluciones linealmente independientes, por lo tanto, la solución es
2
 n 

 n    a2 t
T  x, t    An sen  xe  L  (2.1.19)
n 1  L 
donde las An , n  1,2,3,..... son llamadas constantes de combinación lineal. La determinación de
estas constantes se efectúa aplicando la condición inicial (la condición no homogénea), por lo tanto,
aplicando (2.1.2) en (2.1.19) se tiene

 n 
f  x    An sen x (2.1.20)
n 1  L 
La ecuación anterior indica que f x  puede ser expresada como una serie de Fourier de senos y

que los valores de las An , n = 1,2,3,…, deben satisfacer dicha igualdad. Para determinar dichos
valores se aplica la propiedad de ortogonalidad de funciones (similar a la ortogonalidad de vectores).
La ortogonalización se efectúa multiplicando toda la ecuación, por la misma función característica,
con diferente valor característico y se integra en el dominio (0 ≤ x ≤ L).

L  m   L  n   m 
 f  x  sen  x dx   An  sen  x sen  x  dx, m  1, 2, 3,... (2.1.21)
0
 L  n 1
0
 L   L 
Hay que notar, que en la anterior ecuación se define un vector de ecuaciones, una ecuación para
cada valor de m. Se puede mostrar por integración directa, que la integral del lado derecho de la
ecuación (2.1.21) es igual a cero si m  n e igual a L / 2 si m  n , por lo tanto, todos los términos
de la sumatoria exceptuando uno son cero, en cada ecuación.
Integración de (2.1.21) si m  n
𝑛𝜋 𝑚𝜋
𝑠𝑒𝑛 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥
𝐿 𝐿

1 (𝑛 − 𝑚)𝜋 1 (𝑛 + 𝑚)𝜋
= 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 𝑥 − 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 𝑥
2 − 𝐿 2 + 𝐿
𝐿 𝐿 𝐿 𝐿
1 (𝑛 − 𝑚)𝜋 1 (𝑛 + 𝑚)𝜋
= 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 𝐿 − 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 𝐿
2 − 𝐿 2 + 𝐿
𝐿 𝐿 𝐿 𝐿
1 (𝑛 − 𝑚)𝜋 1 (𝑛 + 𝑚)𝜋
− 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 0 + 𝑛𝜋 𝑚𝜋 𝑠𝑒𝑛 0 =0
2 − 𝐿 2 + 𝐿
𝐿 𝐿 𝐿 𝐿
Integración de (2.1.21) si m  n

𝑛𝜋 𝑥 1 2𝑛𝜋 𝐿 1 0 1 𝐿
𝑠𝑒𝑛 𝑥 𝑑𝑥 = − 𝑛𝜋 𝑠𝑒𝑛 𝑥 = − 𝑛𝜋 𝑠𝑒𝑛(2𝑛𝜋) − + 𝑛𝜋 𝑠𝑒𝑛(0) =
𝐿 2 4 𝐿 2 4 2 4 2
𝐿 𝐿 𝐿

De esta manera se puede despejar para An y obtener

2 L  n 
An   f  x  sen  x dx, n  1, 2,3,.... (2.1.22)
L 0
 L 
Al sustituir los valores de las constante An de la serie de Fourier de senos en la ecuación (2.1.19)
se tiene la solución completamente determinada para el problema planteado
2
 n 
2  L  n   n    a2 t
T  x, t     f  x  sen  x dx  sen  xe  L  (2.1.23)
L n 1 0  L   L 
Esta serie es convergente por lo que se puede truncar la sumatoria hasta un valor determinado de n,
dependiendo de la precisión deseada.
2
 n 
2  L  n   n    a2 t
T  x, t     f  z  sen  z dz  sen  xe  L 
L n 1 0  L   L 

Ejemplo
Como caso particular se considerará una barra de Fe de 10 cm de longitud y difusividad térmica a2
= 0.15 cgs sujeta a una temperatura inicial dada por la ecuación
 T(x,0) = 100, 0 < x <L.
Sustituyendo el perfil térmico anterior en la ecuación (2.1.22) se obtiene que

An  200
1   1 
n

n
así que

f  x   100  200

1   1  sen n x 
n

n 1 n  l 
La expresión anterior indica calcular la sumatoria desde n = 1 hasta n = , pero esto además de
imposible, es innecesario, debido a que la serie es convergente, es decir, que conforme aumenta el
valor de n, los valores de An son más pequeños, hasta que se pueden despreciar. En la Tabla 1 se
muestran algunos valores de An.

Tabla 2.1. Valores de An.

n An N An n An N An
1 127.3240 11 14.1471 21 7.48964 201 0.633453
3 42.4413 13 11.5749 25 6.70126 501 0.25414
5 25.4648 15 9.79415 51 2.49655 801 0.158956
9 18.1891 17 8.48826 101 1.26063

Al sustituir los valores de An, en la solución se obtiene que la temperatura está dada por

1   1n  2  n

2

 sen  n

 a
T  x, t   200 
  t
 xe  L 

n 1 n  L 

1   1n   n 
2

   n  0.15 10  t
T  x, t   200 sen  xe
n 1 n  10 

Hay que notar, que la serie escrita en esta ecuación converge más rápidamente que la constante de
combinación lineal, a causa de que el factor exponencial dependiente del tiempo es negativo y
contiene n2 el cual se incrementa en valor rápidamente, por ejemplo, para t = 10 seg, se tiene que

T x,10   109.80 sen 0.314 x   11.20 sen 0.3942 x   0.63 sen 1.571x   0.01 sen 2.199 x  
 8.8  10 5 sen 2.827 x   1.9  10  7 sen 3.456 x   1.3  10 10 sen 4.084 x   
en este caso, el 5º término de la serie ya es despreciable a pesar que A5 es aproximadamente seis
veces menor que A1 . En particular para x = 5 cm, el cálculo de la temperatura conduce a la serie:

T 5,1  109.80  11.1979  0.628886  0.0128638 8.766  10 5  1.922  10 7  1.334  10 10    99.2215

En la figura 2, se muestra el comportamiento de la temperatura en la barra a lo largo de tiempo,

como se observa, al inicio se tiene el perfil de temperatura parabólico, f(x).

Figura 2.2. Transferencia de calor en una barra unidireccional en estado

no estacionario
Pasos para la solución de un problema por medio del método de Fourier (MSV)
1. Verificar que se pueda aplicar el MSV (a. Problema debe ser homogéneo; b.
Contar con un par de condiciones homogéneas prescritas en la misma
variable y en diferente punto de la misma; c. Contar con al menos una
condición No homogénea).
2. Si el problema No es homogéneo, hay la posibilidad de hacerlo homogéneo
mediante un cambio de variable. 2 bis. Si el problema Es Homogéneo,
proseguir con el paso #3
3. Proponer el cambio de variable que hace la separación (Ej. 𝑇(𝑥, 𝑦, 𝑡) =
𝑀(𝑥)𝑁(𝑦)𝜃(𝑡)).
4. Sustituir en la EDP para obtener el mismo número de EDO’s que de variables
independientes, separando variables y sustituir en las condiciones
homogéneas.
5. Resolver las EDO’s.
6. Aplicar las condiciones homogéneas  se obtiene el problema Sturm-
Liouville (con ecuación característica, valores característicos y función
característica).
7. Regresar a la variable original.
8. Aplicamos principio de superposición de soluciones.
9. Aplicamos condición No-Homogénea.
10. Ortogonalizar la serie para encontrar las constantes de combinación lineal.
11. Sustituir las constantes en la respuesta antes obtenida (en 8).

2.2 Problema de Sturm-Liouville.

El método de separación de variables es el indicado para resolver analíticamente los

problemas de la física - matemática homogéneos. Como ya se vio en el ejemplo anterior, el objetivo
de este método es la generación de varias EDO a partir de la EDP, y la obtención de condiciones
para cada EDO, a partir de las CF homogéneas. Para que el MSV pueda aplicarse, al menos una de
las EDO’s obtenidas y sus condiciones deben formar un problema de Sturm-Liouville. Lo cual
asegura la existencia de funciones ortogonales y la posibilidad de determinar las constantes de
combinación lineal.
El problema de Sturm-Liouville está definido por la ecuación diferencial
dy  x  
d 

dx 
r x 
dx    
 q  x   2 p  x  y  x   0 (2.2.1)

y dos condiciones homogéneas

dy a 
A1 y a   A2  0, (2.2.2)
dx
dy b 
B1 y b   B2  0. (2.2.3)
dx
Dependiendo del valor particular de las constantes A1, A2, B1 y B2, se obtienen condiciones de
primero, segundo y tercer tipo, homogéneas.

Sean yn(x) y ym(x) dos funciones características que se corresponden a dos valores característicos n
y m y que satisfacen al problema de Sturm Liouville. Sustituyendo estas dos funciones y valores
característicos en la ecuación (2.2.1) se tiene lo siguiente.
dy  x  
d 

dx 
 
r  x  n   q  x   2n p x  y n  x   0
dx 
dy  x  
d 

dx 

r  x  m   q  x   2m p  x  y m  x   0
dx 

Multiplicando la primera ecuación por ym(x), la segunda por yn(x) y restando las ecuaciones
resultantes se obtiene

d  dyn  x   d  dym  x  
ym  x  r  x    yn  x   r  x      m   n  p  x  yn  x  y m  x   0
2 2

dx  dx  dx  dx 

Integrando la anterior en el dominio, a  x  b

 2
m   2n  pxy x y x dx   y x  dxd r x dydxx  dx   y x dxd r x dydxx dx
b
n m
b
m
n
b
n
m
a a
  a
 

dy  x  dy  x  dy m  x 
b

 r x  y m x  n   r x  n
b
dx
dx a a dx dx
dy m  x  dy  x  dy n  x 
b

 r x  y n x    r x  m
b
dx
dx a a dx dx

y m b  y n b  y m a  y n a 
 r b  dy m b  dy n b   r a  dy m a  dy n a  (2.2.4)
dx dx dx dx
Pero, las funciones características también deben satisfacer las condiciones (2.2.2) y (2.2.3), por lo
cual se tiene
dy n a  dy n a 
A1 y n a   A2 0 y n a 
dx  dx  0
dy a  dy m a 
A1 y m a   A2 m 0 y m a 
dx dx

dy n b  dy n b 
B1 y n b   B2 0 y n b 
dx  dx  0
dy b  dy m b 
B1 y m b   b2 m 0 y m b 
dx dx

La, igualdad a cero de los determinantes anteriores se explica a causa de que un sistema de
ecuaciones homogéneo tiene solución si y solo si su determinante es igual a cero. Por lo tanto, la
ecuación (2.2.4) se reduce a

 2
n   2m  px y x y x dx  0
b

a
n m (2.2.5)

Para n  m se obtiene que

 px  y x  y x dx  0
b
n m n   m (2.2.6)
a

La ecuación anterior se le conoce como integral de ortogonalidad con función de peso, p(x).
Una generalización del concepto anterior se obtiene al considerar la expansión de una función
seccionalmente continua, f(x), en una serie generalizada de Fourier

f  x    Ann  x  , (2.2.7)
n 1

donde n  x  es una función característica obtenida de un problema de Sturm Liouville, con

propiedad de ortogonalidad. Las constantes de combinación lineal, An, se pueden encontrar al
multiplicar la ecuación (2.2.7) por la función de peso, p(x), y la función característica con otro valor
característico, λm(x), para luego, integrar en el dominio, a  x  b ,
 
p  x  f  x  m  x dx   An  p x  n  x  m  x dx
b b
 , m  1,2,3,  (2.2.8)
a a
n 1

de acuerdo a la ecuación (2.2.6), todas las integrales del lado derecho son cero, a excepción de una,
cuando m = n, por lo tanto, la constante de combinación lineal es

 px  f x  x dx
b
n
An  a
, n  1,2,3,  (2.2.9)
n x 
2

en donde n  x  , es la magnitud o norma de la función característica, n  x  , definida por

 n x    px
b

a
n
2
x dx 
1
2
, n  1,2,3,  (2.1.10)

hay que notar la semejanza entre la norma definida en la ecuación (2.2.10) y la norma para un
conjunto de vectores n-dimensional con componentes vk, dada por
n
v  v
k 1
2
k . (2.2.11)

3. Problema que genera una Serie de Fourier de Cosenos

El siguiente ejemplo genera la serie de Fourier de Cosenos.

Una sustancia disuelta con concentración inicial constante, C0, encerrada entre los planos x = 0 y x
= h, se empieza a difundir a t = 0, en un solvente localizado entre los planos x = h y x = L. Se desea
determinar el proceso de compensación de concentraciones, suponiendo que los planos x = 0 y x =
L son impermeables a la sustancia.

C=C0 C= 0

x x=h x=L

La ecuación de difusión que gobierna este proceso es

 C  x, t   2 C  x, t 
D , 0  x  L, t  0 , (2.3.1)H
t  x2
Con la CI (seccionalmente continua)
 C 0 xh
C  x, 0    0 , (2.3.2)NH
0 h xL
y CF (fronteras aisladas)
 C 0 ,t 
 0, t  0 , (2.3.3)H
x
 C  L, t 
 0, t  0 . (2.3.4)H
x
SOLUCIÓN: Matemáticamente la ecuación (2.3.1) es una EDP de segundo orden, homogénea, y las
CF son de segundo tipo homogéneas, además de una condición No homogénea, con lo cual se tiene
un problema homogéneo que puede ser resuelto mediante el MSV.
Se propone solución
C  x ,t   X  x  t  , (2.3.5)
con la cual, y de forma similar al procedimiento seguido en el problema anterior, tenemos: (2.3.5)
en (2.3.1)
( ) ( ) ( ) ( ) ( ) ( )
=𝐷  𝑋(𝑥) = 𝐷𝜃(𝑡)  ÷ 𝐷𝑋(𝑥)𝜃(𝑡) 
( ) ( )
= como se separaron ya las variables, igualo la anterior al parámetro de
( ) ( )

separación 𝜆 con signo contrario al termino donde se tenga el par de condiciones homogéneas
( ) ( ) ( ) ( )
= = −𝜆  = −𝜆 & = −𝜆 reacomodando
( ) ( ) ( ) ( )

( ) ( )
+ 𝐷𝜃(𝑡)𝜆 = 0 & + 𝑋(𝑥)𝜆 = 0

y se llega a las soluciones

t   C1e   D t ,
2
(2.3.6)
X  x   C 2 cosx   C 3 senx  . (2.3.7)
Mediante la sustitución de (2.3.5) en (2.3.3) y (2.3.4) se obtiene
dX 0 
 0, (2.3.8)
dx
dX  L 
y 0. (2.3.9)
dx
Aplicando la condición (2.3.8) en la ecuación (2.3.7) se tiene

0 1
dX 0
 0  C 2 sen0  C 3  cos0  ,
dx
de donde, C 3  0 y/o   0 .

Primero, si C 3  0 , la ecuación (2.3.7) es

X x   C2 cosx  , (2.3.10)
Con la condición (2.3.9) en la anterior se obtiene
dX  l 
 0  C2  sen   L  ,
dx
de donde se encuentra la ecuación característica
sen   L   0 ,

la cual se cumple para los valores característicos

n
 , n  1, 2,3...... (2.3.11)
L
Por lo tanto, la solución de (2.3.10) es
 n 
X  x   C2 cos  x  , n  1, 2,3,..... (2.3.12)
 L 
y de (2.3.6) es
2
 n 
  Dt
  t   C1e  L 
, n  1, 2, 3,.... (2.3.13)

De acuerdo a la ecuación (2.3.5) se tiene

2
 n 
 n    Dt
C  x, t   C1C2 cos  xe  L  , n  1, 2, 3,.... , (2.3.14)
 L 
Segundo, si   0 , (2.3.12) y (2.3.13) son
X  x   C2 1 ,
t   C1 1 ,
por lo cual, se tiene otra solución
C  x ,t   C1C2 1 . (2.3.15)
La combinación lineal de las soluciones (2.3.14) y (2.3.15) es
2
 n 

 n    Dt
C  x, t   A0   An cos  xe  L  . (2.3.16)
n 1  L 
Aplicando la condición inicial (2.3.2) se obtiene

 n 
C  x, 0   A0   An cos  x. (2.3.17)
n 0  L 
La ecuación (2.3.17) representa la expansión en serie de Fourier de cosenos de C(x,0).
Ortogonalizando la anterior se tiene

 m   m  
 m   n 
x  dx   An  cos 
L L L
 C  x, 0  cos  x dx   A0 cos  x  cos  x  dx, n  0,1, 2,3,....
0
 L  0
 L  n 0
0
 L   L 
(2.3.18)

al sustituir el valor de la concentración inicial en la integral de la izquierda de la ecuación (2.3.18),

ésta puede ser separada en dos integrales,

L  m  h  m 
 C  x, 0  cos  L x  dx  
0 0
C0 cos 
 L
x  dx

L  m  CL  m 
   0  cos  x  dx  0 sen  h  , m  0,1, 2,3,...
h
 L  m  L 

La integral de la derecha de la ecuación (2.3.18) es igual a cero si m  n, a L/2 si m = n  0, y si m

= n = 0 la integral es igual a L.
L
L  m  L   m 
La integral intermedia 0
A0 cos 
 L
x  dx  A0
 m  sen  L
 
x   0
 0
Por lo tanto, de la ecuación (2.3.18) se obtiene que
2C0  n 
An  sen  h  , n  1, 2,3,....
n  L 
Y para obtener el valor de A0, es necesario hacer otra ortogonalización usando la función
característica =1

 n 
C  x, 0 dx   A0 dx   An  cos 
L L L
0 0
n0
0
 L
x dx

L

L   n 
0dx  A0  x 0   An
h L
 C0 dx  
L

n  sen  L x 
0 h
n 0   0

C0 h  A0 L

h
A0  C0
L
Al sustituir las anteriores en la ecuación (2.3.16) se tiene
 n
sen h  n 
L cos  n

h 2C    Dt
C  x, t   C0  0
L 
 n0 n

 L
xe  L 

(2.3.19)

Hay que notar que la solución (2.3.19), es la suma de una solución en estado estacionario y una
solución transitoria, esta última tiende a cero cuando t tiende a infinito.
En la Figura 2.3 se presenta el comportamiento para el caso particular en que l
h  3 . Hay que notar
que la solución en estado estacionario es igual a 1/3 de la concentración inicial, debido a la selección
del factor l
h  3.

t0

1 

 t1

0.8  t2
Concentración, C(x,t)/Co

 t3
0.6 t4

t5

0.4  tiempo t
 t6 


t  t6
tiempo   t4 
 t3 t5
0.2
 t2
 t1
0 
t0

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Posición, x/l

Figura 2.3. Transferencia de masa unidireccional en estado no

estacionario dentro de un recipiente inicialmente con
dos secciones, hl  3 .

C<C0 C>0

x x=h x=L