Demostracion Matematica de La Ecuacion de Difusividad en Condiciones de Fronteras
Demostracion Matematica de La Ecuacion de Difusividad en Condiciones de Fronteras
Demostracion Matematica de La Ecuacion de Difusividad en Condiciones de Fronteras
U(x,s)
L u(x,t) = (4)
[ x ] x
Derivando parcialmente con respecto a t, la transformada de Laplace sera:
Z
u(x,t)
L u(x,t) = est dt (5)
[ t ]
0 t
Se reconoce:
u
= ut
t
Reemplazando: Z
L u(x,t) = est ut dt
[ ]
t 0
du
Luego se convierte a derivada ordinaria: ut = dt
Reemplazando: Z
L u(x,t) = est du
[ ]
t 0
Integracion por partes:
Z
L u(x,t) = est u(x,t) |
0 +s est dt
[ ]
t 0
Z
st
L u(x,t) = lim (e u(x,t) ) u(x,0) + s est dt
[ ] t
t 0
1
De la misma manera que en la ecuacion (4), se obtiene la transformada de Laplace para la segunda derivada
parcial de u(x,t) con respecto a x
2 U(x,s)
L 2 u(x,t) = (7)
[ x2 ] x2
Luego de realizar las demostraciones pertinentes para obtener las propiedades a utilizar, se pasa a resolver
la Ecuacion de difusividad:
2 p 1 p Ct p
2
+ = (8)
r r r k t
Para facilitar las operaciones, se utiliza la ecuacion adimensional:
2 pD 1 pD pD
2 + = (9)
rD rD rD tD
Siendo pD una funcion de variables rD y tD , se lleva la ecuacion (9) al espacio de Laplace:
1
L 2 pD + L pD = L[ pD ] (10)
[
r 2
D
] rD [ rD ] tD
1
L[p(1,tD ) ] = U(1,s) = (16)
s
La solucion para la ecuacion (15) en el espacio de Laplace es la siguiente combinacion lineal1 :
U(rD ,s) = A0 I0 (rD s) + B0 K0 (rD s) (17)
1 A y B :Coeficientes; I : Funcion modificada de Bessel de primera especie y orden cero; K : Funcion de Bessel de segunda
0 0 0 0
especie y de orden cero
2 I : Funcion modificada de Bessel de primera especie y de primer orden; K : Funcion de Bessel de segunda especie y de primer
1 1
orden
2
Los coeficientes luego de despejar las ecuaciones (18) y (19) son:
K1 (reD s)
A0 = (20)
s[I0 ( s)K1 (reD s) + K0 ( s)I1 (reD s)]
I1 (reD s)
B0 = (21)
s[I0 ( s)K1 (reD s) + K0 ( s)I1 (reD s)]
Se define qD :
pD
qD = rD (22)
rD rD =1
Definiendo la transformada de Laplace de qD :
U(rD ,s)
L[qD (rD ,t) ] = M(rD ,s) = (23)
rD rD =1
L[qD (rD ,t) ] = M(rD ,s) = [A0 I1 ( s) + B0 K1 ( s)] (24)
2 )
min(i, N N
N
+i
X k 2 (2k)!
V (i) = (1) 2
N
(26)
2 k !k! (k 1)! (i k)! (2k i)!
k=floor( i+1
2 )
3 Las ecuaciones estan escritas de acuerdo al codigo de programacion usado en wxM axima
3
4
Figure 1: qD(1) y qD(2) en escala cartesiana