V.-Condición de Contorno Isotérmica en Sólidos Infinitos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 11

V.

- CONDICIÓN DE CONTORNO ISOTÉRMICA

EN SÓLIDOS INFINITOS
pfernandezdiez.es

V.1.- CONDUCCIÓN TRANSITORIA EN PLACA INFINITA CON CONDICIÓN DE CONTOR-


NO ISOTÉRMICA

La conducción a través de una placa plana de espesor finito L en la dirección x, y de espesor infini-
to en las otras dos, por lo que en éstas se desprecian los efectos de borde, la ecuación diferencial de la
conducción es:

∂ 2Φ 1 ∂Φ
= , con: Φ = T(x, t) - TF
∂x 2 α ∂t
La zona próxima a la periferia es de resistencia superficial despreciable por lo que Fo < 1 y desde
ahí hasta el núcleo, Fo > 1.

Sí es posible representar la distribución de temperaturas, Φ = T - TF, mediante una expresión de


la forma:

⎧ ∂Φ = ∂X θ ⇒ ∂ 2Φ = ∂ 2 X θ ⎫
⎪ ⎪
Φ = X ( x) θ ( t ) ⇒ ⎨ ∂x ∂x ∂x 2 ∂x 2 ⎬ ⇒ 1 ∂ 2 X = 1 ∂θ = - λ 2
⎪ ∂Φ = ∂θ X ⎪ X ∂x 2 α θ ∂t
⎩ ∂t ∂t ⎭

El parámetro λ2 se ha introducido por cuanto cada uno de los miembros de dicha igualdad es fun-
ción de una sola variable; esta separación de variables conduce al siguiente sistema de dos ecuaciones
diferenciales:

∂θ + α λ 2 θ = 0 ; θ = C1 e − α λ
2t

∂t

∂2 X
+ λ2 X = 0 ; X = C2 sen ( λ x) + C3 cos ( λ x )
∂x 2
Fig V.1.- Placa plana infinita
La solución general es:
2 2
Φ = C1 e - α λ t {C2 sen ( λ x ) + C3 cos ( λ x )} = e - α λ t { B1 sen ( λ x ) + B2 cos ( λ x )}
pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-121
Condición inicial, para:{ 0t =≤ x0 ≤ L ⇒ Φ = Φ ( x , 0 ) = f(x) ó Φ = T - T 0 0 F

Condiciones de contorno, para: {


t> 0 ⎧ en: x = 0 ; Φ = 0 = T ( 0 , t )
-∞ ≤ x ≤ +∞ ⇒ ⎨⎩ en: x = L ; Φ = 0 = T ( L , t )

La condición de contorno se presenta para el caso límite de considerar un valor muy grande del
coeficiente de transferencia térmica por convección, {metales líquidos}, por lo que la resistencia térmi-
ca de la capa de convección es despreciable y la temperatura de la superficie del cuerpo en el tiempo t
es idéntica a la temperatura del fluido, situación a la que se debe llegar en un tiempo muy pequeño
(condición de contorno isotérmica).

Aplicando las condiciones de contorno a la ecuación diferencial:

⎧ x = 0 , Φ = 0 ⇒ B2 = 0
Para: ⎨ 2
⎩ x = L , Φ = 0 ⇒ 0 = e - α λ t B1 sen ( λ L ) ⇒ sen ( λ L ) = 0

y considerando que para cualquier valor de t, se satisface por un número infinito de valores del pará-
metro (λ L) se puede poner:

sen ( λ n L) = 0 ⇒ λn = π n ; n = 1, 2 , 3,...
L
Para cada valor de n se obtiene uno de λn siendo éstos los valores característicos del problema, por
lo que la solución para la distribución de temperaturas es un desarrollo en serie de la forma:


∑ e- α λ t
2
Φ= n Bn sen ( λ n x)
n=1


Aplicando la condición inicial, t = 0, Φ = f(x), resulta: Φ = f(x) = ∑ Bn sen ( λ n x)
n=1

En una serie infinita de funciones de la forma:

sen ( λ1 x ), sen ( λ 2 x ), sen ( λ3 x ), ..., sen ( λn x )

éstas son ortogonales, cuando se cumple que:



L
∫0 sen ( λi x ) sen ( λ j x ) dx = 0 , con: i ≠ j

y tiene un valor determinado en un instante considerado.

€ Si la distribución de temperaturas Φ = f(x) es una función arbitraria se puede poner en función de


una combinación lineal de funciones ortogonales, en la forma:


f ( x ) = B1 sen ( λ 1 x ) + B2 sen ( λ 2 x ) + ... + Bn sen ( λ n x ) + ... = ∑ Bn sen ( λ n x )
n=1

en la que los valores de Bi son constantes a determinar.

Si la serie anterior es convergente e integrable, y la multiplicamos por, sen ( λ n x ) , se obtiene:


L L L
∫0 f ( x ) sen ( λn x ) dx = B1 ∫0 sen ( λ1 x ) sen ( λn x ) dx + ... + Bn ∫0 sen 2 ( λn x ) dx
pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-122


Por definición de ortogonalidad se hacen cero todas las integrales del segundo miembro, menos la
correspondiente al coeficiente Bn, por lo que:

L
∫0 f ( x ) sen ( λn x ) dx
Bn = L
=
∫0 sen 2 ( λn x ) dx

L L sen ( λn L) cos ( λn L) ⎧ π n ⎫ L 2 L


= ∫ 0
sen 2 ( λn x) dx =
2
-
2 λn
= ⎨ λn =
⎩
⎬ =
L ⎭ 2
=
L ∫ 0
f(x) sen ( λn x) dx

obteniéndose la siguiente distribución de temperaturas:

€ πn
2 ∞
-α λn2 t
L 2 ∞ - ( )2 α t πnx L πnx
Φ =
L

n=1
e sen (λn x) ∫ 0
f(x) sen (λn x) dx =
L

n=1
e L sen
L ∫ 0
f(x) sen
L
dx

que depende del tiempo t, de la posición x y de la distribución de temperaturas inicial f(x), que es co-
nocida.

Si se considera f(x) = Φ0, el valor de Bn es:

2 L
Bn =
L ∫ 0
f (x) sen ( λn x) dx =

2 Φ0 L 2 Φ0 L πn 2 Φ0 4 Φ0
=
L ∫ 0
sen ( λn x) dx =
L ∫ 0
sen (
L
x) dx =
πn
{1 - (-1)n } =
πn

por lo que:

Φ 2 €∞ -(π n)2 α t πnx L πnx 2 ∞


{1 - (-1)n } sen ( λn x) 4 ∞ sen ( λn x )
Φ0
= ∑
L n=1
e L sen
L ∫ 0
sen
L
dx =
π

n=1
e-λ2nα t
n
= ∑
π n=1,3 ,..
e-λ2nα t
n

Si se pone en función del nº de Fo y de un parámetro adimensional de la posición ξ = x :


L


Φ 2 2 L πnx 4 ∞ 2 sen ( π n ξ )
= ∑ e-(π n) Fo
sen ( π n ξ ) ∫0 sen dx = ∑ e-(π n) Fo
Φ0 L n=1
L π n=1,3.. n

Si la distribución de temperaturas fuese, por ejemplo, de la forma, f(x) = a x, con, a = Cte, se obtie-
€ ne:

2 L 2a L 2aL
Bn = ∫0 f ( x ) sen ( λn x ) dx = ∫0 x sen ( λn x ) dx = ... = (-1)n+1
L L πn

(-1 )n+1
Φ = 2 aL
2

π ∑ e-( π n) Fo sen ( π n ξ ) n
€ n=1

cuya representación gráfica exponemos en la Fig V.2.

El flujo térmico en el instante t es:

∂Φ 2 Φ0 k ∞
2 1 - (-1)n π n 4 Φ0 k ∞ 2 4 Φ0 k ∞ 2
q=-k( )x=0 = ∑ e-(π n) Fo
= ∑ e-(π n) Fo = ∑ e- λ n α t
∂x π n=1
n L L n=1,3... L n=1,3...

pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-123


Estas series convergen muy rápidamente a menos que el nº de Fo sea muy pequeño.

Para Fo > 0,2 sólo es necesario conservar el primer término de la serie, cometiéndose un error me-
nor del 2%. Para valores muy pequeños del nº de Fo < 0,2, (es decir, poco después de que la placa se
haya sumergido en el líquido), la serie converge lentamente y es necesario conservar los suficientes
términos del desarrollo para obtener un resultado exacto o hacer uso de la gráfica V.3.

Fig V.2.- Desarrollo temporal de la distribución de temperaturas en una placa plana infinita, en régimen transitorio,
con temperatura inicial de la forma, f(x)= at

Fig V.3.- Abaco para la determinación gráfica de la temperatura y de Q/Q0


en una placa plana infinita en régimen transitorio, con condición de contorno isotérmica, para Fo < 0,2

k ( T0 - TF )
Para Fo→ 0, lím qFo → 0 = , que es la expresión que se encontrará para el sólido semiin-
π αt
finito, y que indica que la temperatura sólo cambia en una posición muy delgada cerca de la superfi-
cie, comportándose en esta zona como un sólido semiinfinito, mientras que la temperatura en el inte-
rior de la placa permanece constante; el proceso de la transmisión de calor se limita a esta delgada re-
gión y el espesor de la placa no afecta en absoluto, por cuanto L no aparece en la ecuación.
pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-124
El calor Q que abandona la placa por las dos caras, de superficie de contacto con el fluido (2 A), en
el intervalo de tiempo (0 ÷ t), es:


t 4 (T0 - TF ) k A t 2
Q = 2A ∫0 q dt =
π
∫ 0 ∑ e-(π n) Fo
{ 1 - (-1)n } dt =
n=1


4 (T0 - TF ) k A L2 8 (T0 - TF ) k A L ∞ 1 - e-( π n) 2 Fo
= ∑ ( e-(π n)2 Fo - 1) { 1 - (-1)n } = ∑ =
π n=1
( π n )2 α α n=1,3...
( π n) 2


1 - e-( π n ) 2 Fo 8 (T0 - TF ) ρ c p A ∞ 1 - e- λ2n α t
= 8 (T0 - TF ) ρ c p A L ∑ = ∑
n=1,3...
( π n )2 L n=1,3... λ2n

⎧
El calor Q0 almacenado en la placa es, para ⎨ x ≥ 0 , de la forma:
⎩T = TF

Q0 = L A ρ c p Φ 0 = L A ρ c p (T0 - TF )

2 €
Q ∞
1 - e -λ n α t
y la relación: = 8 ∑
€ Q0 n=1,3... ( π n )2

El calor almacenado en la placa en el intervalo t = 0 a t = t, es:

€ 
Q = L A ρ c p (TFinal − T0 )


 1 L TFinal - TF 8 ∞
e-λ2n α t

Φ Final =
L
∫ 0 Φ(x) dx ;
T0 - TF
=
π2

n2
n=1,3...

OTROS CASOS:

a) Pared plana infinita inicialmente a T0 que experimenta cambios instantáneos de temperatura en
las superficies a T1 y T2 . El problema se resuelve mediante la superposición de un problema estacio-
nario y otro transitorio, de forma que:

T ( x , t ) = T1 ( x ) - T2 ( x , t ) ó Φ ( x , t ) = Φ 1 ( x ) - Φ 2 ( x , t )
∂ 2 T1 ∂T2 ∂ 2 T2
=0 ; =α
∂x 2 ∂t dx 2
T = T0 ; 0 < x < L ; t = 0
T = T1 ; x = 0 ; t > 0
T = T2 ; x = L ; t > 0

∞ 2
2 T1 - (-1) n T2 2 2 T0 e-λn α t ( 2 n + 1) π
T ( x,t ) = ∑ { 1 - e-λn α t + } sen ( λ n x ) ; λn =
π n λn L
n=1
.....................................................................................................................................................................................................

€ b) Pared plana infinita inicialmente a T0; una de sus superficies experimenta un aislamiento térmi-
co y la otra intercambia calor con un fluido exterior.

pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-125


T = T0 ; 0 < x < L ; t = 0
∂T
=0 ; x = 0 ; t > 0
∂x
∂T hcF
= - a1 T = - T ; x=L ; t>0
∂x k
∂Φ hcF
= - a1 Φ = - Φ ; x=L ; t> 0
∂x k

∞ 2
Φ(x, t) cos ( λ n x ) e-λn α t ( 2 n + 1) π
= 2 a1 T0 ∑ ; λn =
Φ0
n=1 { L ( λ2n + a12 ) + a1 } cos ( λ n L ) L
.....................................................................................................................................................................................................

€ V.2.- CONDUCCIÓN TRANSITORIA EN PARED CILÍNDRICA INFINITA CON CONDICIÓN


DE CONTORNO ISOTÉRMICA.

Se considerará que en el cilindro infinito no existen efectos de borde, debido a que se le supone
longitud infinita y, por lo tanto, que la conducción del calor se verifica en la dirección radial.

La ecuación diferencial que rige el proceso, en coordenadas cilíndricas, es:

∂2T ∂ 2Φ
α 1 ∂ ( r ∂T ) = ∂T ; α( +1
∂T
)=
∂T
; α( + 1
∂Φ
)=
∂Φ
r dr ∂r ∂t dr 2 r ∂r ∂t dr 2 r ∂r ∂t

con: Φ = T - TF

Si representamos Φ en la forma, Φ = R(r) θ(t), resulta:

∂Φ 2 ∂ 2R
=
∂R
θ ; ∂ Φ2
= θ ; ∂Φ
=
∂θ
R
∂r ∂r dr dr 2 ∂t ∂t

2 1 ( ∂ 2 R + 1 ∂R ) = 1 ∂θ = - λ 2
α(∂ R θ + θ ∂R ) = ∂θ R ;
Fig V.4.- Cilindro infinito dr 2 r ∂r ∂t R dr 2 r ∂r α θ ∂t

Las ecuaciones diferenciales ordinarias resultantes y sus soluciones son:

d2R dR
+ 1 = - λ2 R ⇒ R = B1 J0 ( λ r ) + B2 Y0 ( λ r )
dr 2 r dr
d θ = - λ 2α dt ⇒ θ = B3 e - λ α t
2

Como el cilindro no puede admitir en su eje (r = 0), una solución infinita, por cuanto Y0(0) = ∞, re-
sulta que B2 tiene que ser cero, por lo que se obtiene una ecuación de la forma:

R = B1 J 0 ( λ r )

La solución general que proporciona la distribución de temperaturas de la forma:


2 2
Φ = B3 e - λ α t B1 J 0 ( λ r ) = B e - λ α t J 0 ( λ r )

en la que B y λ son constantes, que habrá que determinar mediante las condiciones de contorno;
J0(λ r) es la función de Bessel de primera especie y orden cero.

pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-126


t=0
La condición inicial es: ⎧⎨ ⇒ Φ = f (r ) ó Φ 0
⎩ 0 ≤ r ≤ R

La condición de contorno para un cambio brusco de la temperatura en la superficie lateral del ci-
lindro infinito a: Φ = T - TF, es:

2
t > 0 ; Φ r=R = 0 = B e - λ α t J 0 ( λ R ) ; J 0 ( λ R ) = 0 ; J0 ( λ n R) = 0

que se tiene que cumplir para cualquier valor de t.

Raices de las funciones de Bessel: J0(x) = 0, y J1(x) = 0


J0 ( x ) = 0 2,4048 5,5201 8,6537 11,7915 14,9309 18,0711
J1 ( x) = 0 3,8317 7,0156 10,1735 13,3237 16,4706 19,6159

Los valores de λn se calculan como raíces de la ecuación, J0(λnR) = 0, con, n = 1, 2, 3,... obteniéndo-
se una serie de valores que conforman un desarrollo en serie para la distribución de temperaturas, de
la forma:


Bn e -λnα t J 0 ( λ n r )
2
Φ= ∑
n=1

⎧ t = 0
Aplicando la condición inicial ⎨ , resulta:
⎩ F = f(r)

f ( r ) = B1 J 0 ( λ 1 r ) + B2 J0 ( λ 2 r ) + ... + Bn J0 ( λ n r ) = ∑ Bn J 0 ( λn r )
n=1

Para que ésto sea así, es necesario que las funciones , J0 ( λ 1 r ), J0 ( λ 2 r ), ... , J 0 ( λ n r ), formen un
agrupamiento ortogonal en el intervalo, 0 ≤ r ≤ R, respecto a un factor ponderal r, de forma que:

R
∫0 r J 0 ( λir ) J 0 ( λ j r ) dr = 0, con: i ≠ j

Si la serie es convergente e integrable, se tiene que:

€ R
∫ 0 r f (r ) J 0 ( λnr ) dr =

R R R
= B1 ∫ r J 0 ( λ1r ) J 0 ( λnr ) dr + B2 ∫ r J 0 ( λ 2r ) J 0 ( λnr ) dr + ... + Bn ∫ r J 02 ( λnr ) dr + ...
0 0 0

Por definición de ortogonalidad, todas las integrales del segundo miembro a excepción de la últi-
ma, son cero, es decir:

R
R R ∫0 r f (r ) J 0 ( λnr ) dr
∫0 r f (r ) J 0 ( λnr ) dr = Bn ∫ r J 02 ( λnr ) dr ⇒ Bn = R
=
0
∫0 r J 02 ( λnr ) dr
R
R R2 R 2 J 12 ( λn R) ∫0 r f (r ) J 0 ( λnr ) dr
= ∫0 r J 02 ( λnr ) dr = { J 02 ( λn R) + J 12 ( λn R)} = =
2 2 R 2 J 12 ( λn R)

2
obteniéndose la siguiente distribución de temperaturas:


pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-127
R
∞ ∫0 r f (r ) J 0 ( λnr ) dr 2 2 ∞ 2 J 0 ( λn R) R
Φ= ∑ 2
e- λn α t J 0 ( λnr ) = ∑ e- λn a t ∫0 r f (r ) J 0 ( λnr ) dr
n=1 R J 12 ( λn R) R 2
n=1 J 12 ( λn R)
2
en la que J1(λnR) es la función de Bessel de primera especie y orden uno.

€ Si la distribución de la temperatura inicial es uniforme: Φ0 = T0 - TF

f ( r ) = Φ ( r , 0 ) = Φ 0 = Cte, para: 0 ≤ r ≤ R

Teniendo en cuenta que J0 (λnR) = 0, resulta:

R R
∫0 r J 0 ( λnr ) dr = J ( λ R)
λn 1 n

Φ = T - TF = 2 e - λn α t
2 J 0 ( λn r )
Φ0 T0 - TF R ∑ λ n J 1 ( λ n R)
n=1

siendo la temperatura Tc en el eje del cilindro (r = 0) para: Φc = Tc - TF:

Φc T - TF e -λ n α t
2

= c = 2 ∑
Φ0 T0 - TF R n=1
λ n J1 ( λ n R)

El flujo de calor es:

∂ J ( λ r) = - λ J ( λ r)
0 n n 1 n
∂r 2 Φ 0 ∞ - λ2n α t
q = - k ∂Φ 〉 r= R = = - k ∑ e
- λ2n α t - λ n J 1 ( λ n R ) = - 2 Φ 0
∞ ∞
∂r ∂Φ 〉 2
e- λ n α t R n=1
r= R = 2 Φ 0 ∑ e ∑
∂r n=1
λ n R J 1 ( λn R ) R n=1

y el calor Q que llega a (r = R), en el intervalo (0 ÷ t) es:


2
t t

-λn2α t

1 - e -λ n α t
Q= A ∫0 q dt = - 4 π k L Φ 0 ∫0 ∑ e dt = ... = - 4 π k L Φ 0 ∑
n=1 n=1 λ2n α

El calor Q0 almacenado inicialmente en el cilindro a Φ0 es:

€ Q0 π R2 k Φ
= π R 2 ρ cp Φ0 = 0
L α

Fig V.5.- Desarrollo temporal de la distribución de temperaturas en conducción transitoria,


en un cilindro infinito, con temperatura inicial constante
pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-128
y la fracción de energía perdida:

t ∂Φ
Q -k 2πR ∫ ( ) dt
0 ∂r r= R = - 2 α t ∂Φ
= ∫0 ( ) dt =
∂r r= R
J 0 ( λn R) = 0 =
Q0 π R2 ρ c p Φ0 R Φ0

2 2
e- λn α t 1 - e- λ nα t
∞ ∞ ∞
4α t - λn2α t 4α
= ∫0 ∑ e dt = - ∑ 〉 t0 = 4 ∑
€ R2 n=1 R2 n=1 λ2nα n=1 ( λn R) 2

V.3.- CONDUCCIÓN TRANSITORIA EN UNA ESFERA CON CONDICIÓN DE CONTORNO


ISOTÉRMICA. €

Para estudiar la distribución de temperaturas en una esfera que se calienta, o enfría instantánea-
∂ 2Φ
mente en su superficie, se parte de: α ( 2 + 2
∂Φ ∂Φ
)= , en la que r es el radio correspondiente a
∂r r ∂r ∂t
un punto cualquiera de la esfera.

La solución general es de la forma:


2
αt sen( λ r )
Φ = B e-λ
λr
⎧ t = 0 ; r ≤ R ; Φ ( r ,0 ) = Φ = f ( r ) ó Φ ( r,0 ) = Φ 0
con las siguientes condiciones de contorno: ⎨
⎩ t > 0 ; Φ r=R = 0
Aplicando la segunda condición se obtiene:

2 sen ( λ R)
Φ = 0 = B e- λ α t ⇒ sen ( λ R) = 0
λR

que se debe cumplir para cualquier valor de t, por lo que habrá ∞ soluciones, deduciéndose de ella la
relación que define las λ en la forma:

sen ( λ n R) = 0 ⇒ λn = π n
R

sen ( λ n r )
La solución general de la distribución de temperaturas es: Φ = ∑ Bn e - λ nα t
2

n=1
λ nr

⎧t = 0
Para determinar el valor de la constante Bn aplicamos la primera condición de contorno: ⎨
⎩Φ = f (r )

sen ( λ n r ) sen ( λ 1 r ) sen ( λ 2 r ) sen ( λ n r )
f ( r ) = ∑ Bn = B1 + B2 + ... + Bn
n=1
λn r λ1 r λ 2r λn r

y haciendo uso de las propiedades de las funciones ortogonales se tiene:

R
∫0 f (r ) ( λnr ) sen ( λnr ) dr =

R λn r R λnr R
= B1 ∫ sen ( λ1r ) sen ( λnr ) dr + ... + Bn ∫ sen 2 ( λnr ) dr = Bn ∫ sen 2 ( λnr ) dr
0 λ1r 0 λnr 0

despejando Bn y teniendo en cuenta que, sen (λn R) = 0, resulta:


pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-129
R
∫0 f (r ) ( λnr ) sen ( λnr ) dr 2 R
Bn = R
=
R
∫0 f (r ) ( λnr ) sen ( λnr ) dr
∫0 sen 2 ( λnr ) dr

Los valores de Bn son distintos para cada superficie equipotencial, definida por r.

€ La ecuación general que proporciona la distribución de temperaturas de las superficies isotermas


de radio r a lo largo del tiempo, es de la forma:


2 2 sen ( λnr ) R
Φ= ∑ e- λn α t ∫0 f (r ) ( λnr ) sen ( λnr ) dr
R n=1
λnr

Para: f ( r ) = Φ 0 se tiene:
€ ∞
Φ = 2 sen ( λ n r ) - ( λ n r ) cos ( λ n r ) + sen ( λ n r ) R
e -λ nα t
2
∑ { }0 =
Φ0 R n=1
λn r λn

∞ sen ( π n r ) ∞ sen ( π n r )
-λ 2nα t R = cos ( π n ) = (-1 ) n = 2 R - 2 αt R
= 2 ∑ - cos ( π n ) e
πnr ∑ -(-1 ) n e n
π n=1
λ
nr
n=1
R
Si la temperatura en el centro de la esfera (r = 0) es TC, sale una indeterminación de la forma
(0/0), por lo que aplicando la regla de L’Hôpital se obtiene:


ΦC
= - 2 ∑ (-1 ) n e- λnα t
2

Φ0 n=1

La disipación de calor adimensional es:

t ∂Φ
Q - k 4 π R2 ∫ ( ) dt 3k
0 ∂r r= R t ∂Φ
Q0
= 3
=-
R ρ c Φ
∫ 0 ( ∂r )r= Rdt =
4πR p 0
ρ c p Φ0
3
πnr πn πnr
6k t
∞ cos ( ) (n r ) - sen ( )n
n - λn2 α t R R R
= - ∫ ∑ (-1 ) e 〉 r =R dt =

π ρ c p 0 n=1 n 2 R2
2
6k t

n - λ n2 α t (-1)n π n 2 6k t

- λn2 α t -6k ∞
e- λn α t
= ∫0 ∑ (-1) e 2 2
dt = 2 ∫0 ∑ e dt = 2
∑ 〉 t0 =
€ π ρ cp n=1 n R ρ cp R n=1 ρ cp R n=1 λ2n α

2 2 2 2
-6k ∞
e- λ n α t - 1 -6k ∞
e- λ nα t - 1 6k ∞
1 - e- λ n α t 6 ∞
1 - e- λ n α t
= ∑ = ∑ = ∑ = ∑
€ ρ c p R2 n=1 λ2n α ρ c p R2 n=1 (
πn 2
) α
2
π α ρ cp n=1 n2 π2 n=1 n2
R

V.5.- TRANSMISIÓN DE CALOR POR CONDUCCIÓN EN RÉGIMEN TRANSITORIO CON


€ GENERACIÓN DE CALOR E

a) Pared plana infinita con generación de calor E; condición de contorno isotérmica.

⎧⎪ t = 0 ; Φ = Φ 0 ; 0 ≤ x ≤ L
Condiciones de contorno : ⎨ t > 0 ⎧ x = 0 ; Φ 0 = 0
⎪ ⎨
⎩ ⎩ E = Cte

pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-130


2 ∞ sen ( λ n x )
Φ ( x , t ) = 4 E3 L ( 1 - e - λ n α t ) , con: n = 1, 3 , 5 , ... ; λ n = π n
2

π k
∑ n3 L
n=1
...........................................................................................................................................................
b) Cilindro infinito con generación de calor E; condición de contorno isotérmica.

⎧⎪ t = 0 ; Φ = Φ 0 ; 0 ≤ r ≤ R
Condiciones de contorno: ⎨ t > 0 ⎧ r = R ; Φ 0 = 0
⎪ ⎨
⎩ ⎩ E = Cte
∞ J0 ( λ n r )
Φ ( x, t ) = 4 E (1 - e - λ nα t ) , con λ n raices de: J 0 ( λ n R ) = 0
2

Rk ∑ λ 3n J1 ( λ n R )
n=1

...........................................................................................................................................................
c) Esfera con generación de calor E; condición de contorno isotérmica.
⎧⎪ t = 0 ; Φ = Φ 0 ; 0 ≤ r ≤ R
Condicioners de contorno: ⎨ t > 0 ⎧ r = R ; Φ 0 = 0
⎪ ⎨
⎩ ⎩ E = Cte

2E sen ( λ n r ) 2 πn
Φ ( r,t ) = ∑ (1 - e- λn α t ) (-1) n ; λn =
k
n=1 λ3n r R
...........................................................................................................................................................

€ Valores de las funciones de Bessel J0(x) y J1(x)

x J0 ( x ) J1 ( x ) x J0 ( x) J1 ( x ) x J0 ( x ) J1 ( x ) x J0 ( x ) J1 ( x )
0,00 1,0000 0,0000 1,65 0,4267 0,5743 3,30 -0,3443 0,2206 5,10 -0,1443 -0,3371
0,05 0,9993 0,0249 1,70 0,3979 0,5777 3,35 -0,3548 0,2000 5,20 -0,1103 -0,3432
0,10 0,9975 0,0499 1,75 0,3690 0,5801 3,40 -0,3643 0,1792 5,30 -0,0758 -0,3460
0,15 0,9943 0,0747 1,80 0,3399 0,5815 3,45 -0,3727 0,1583 5,40 -0,0412 -0,3453
0,20 0,9900 0,0995 1,85 0,3109 0,5818 3,55 -0,3864 0,1164 5,50 -0,0068 -0,3414
0,25 0,9844 0,1240 1,90 0,2818 0,5811 3,65 -0,3960 0,0745 5,60 0,0270 -0,3343
0,30 0,9776 0,1483 1,95 0,2528 0,5794 3,70 -0,3923 0,0538 5,70 0,0599 -0,3241
0,35 0,9696 0,1723 2,00 0,2238 0,5767 3,75 -0,4014 0,0332 5,80 0,0917 -0,3110
0,40 0,9604 0,1960 2,05 0,1951 0,5730 3,80 -0,4025 0,0128 5,90 0,1220 -0,2951
0,45 0,9500 0,2193 2,10 0,1666 0,5682 3,85 -0,4026 -0,0073 6,00 0,1506 -0,2767
0,50 0,9384 0,2422 2,15 0,1383 0,5626 3,90 -0,4018 -0,0272 6,10 0,1773 -0,2559
0,55 0,9257 0,2647 2,20 0,1103 0,5559 3,95 -0,3999 -0,0468 6,20 0,2017 -0,2329
0,60 0,9120 0,2867 2,25 0,0827 0,5483 4,00 -0,3971 -0,0660 6,30 0,2238 -0,2081
0,65 0,8971 0,3081 2,30 0,0555 0,5398 4,05 -0,3933 -0,0848 6,40 0,2433 -0,1816
0,70 0,8812 0,3290 2,35 0,0287 0,5304 4,10 -0,3886 -0,1032 6,50 0,2601 -0,1538
0,75 0,8642 0,3492 2,40 0,0025 0,5201 4,15 -0,3830 -0,1212 6,6 0,2740 -0,1250
0,80 0,8462 0,3688 2,45 -0,0232 0,5090 4,20 -0,3765 -0,1386 6,8 0,2931 -0,0625
0,85 0,8273 0,3877 2,50 -0,0483 0,4970 4,25 -0,3692 -0,1555 7,0 0,3001 -0,0047
0,90 0,8075 0,4059 2,55 -0,0729 0,4843 4,30 -0,3610 -0,1719 7,2 0,2951 0,0543
0,95 0,7867 0,4233 2,60 -0,0968 0,4708 4,35 -0,3520 -0,1876 7,4 0,2786 0,1096
1,00 0,7652 0,4400 2,65 -0,1199 0,4565 4,40 -0,3422 -0,2027 7,6 0,2516 0,1592
1,05 0,7428 0,4559 2,70 -0,1424 0,4416 4,45 -0,3370 -0,2172 7,8 0,2554 0,2014
1,10 0,7196 0,4709 2,75 -0,1641 0,4259 4,50 -0,3205 -0,2310 8,0 0,1716 0,2346
1,15 0,6957 0,4950 2,80 -0,1850 0,4097 4,55 -0,3086 -0,2441 8,2 0,1222 0,2580
1,20 0,6711 0,4982 2,85 -0,2051 0,3928 4,60 -0,2961 -0,2565 8,4 0,0692 0,2708
1,25 0,6459 0,5106 2,90 -0,2243 0,3754 4,65 -0,2830 -0,2681 8,6 0,0146 0,2727
1,30 0,6200 0,5220 2,95 -0,2426 0,3574 4,70 -0,2693 -0,2790 8,8 -0,0392 0,2641
1,35 0,5937 0,5324 3,00 -0,2600 0,3390 4,75 -0,2551 -0,2891 9,0 -0,0903 0,2453
1,40 0,5668 0,5419 3,05 -0,2765 0,3201 4,80 -0,2404 -0,2985 9,2 -0,1368 0,2147
1,45 0,5395 0,5504 3,10 -0,2920 0,3009 4,85 -0,2252 -0,3070 9,4 -0,1768 0,1816
1,50 0,5118 0,5579 3,15 -0,3066 0,2812 4,90 -0,2097 -0,3146 9,6 -0,2090 0,1395
1,55 0,4837 0,5644 3,20 -0,3201 0,2613 4,95 0,1938 0,3215 9,8 -0,2323 0,0928
1,60 0,4554 0,5699 3,25 -0,3327 0,2411 5,00 -0,1776 -0,3275 10,0 -0,2459 0,0435

pfernandezdiez.es Condición isotérmica sólidos infinitos.V.-131

También podría gustarte