Capitulo N°5 Parte B

CAPITULO 5
MÉTODOS NUMÉRICOS APLICADOS

A LA CONDUCCIÓN DE CALOR
PARTE B:
5-6.1.- APLICACIÓN EN PROBLEMAS DE
TRANSFERENCIA DE CALOR POR
CONDUCCIÓN
La ecuación de LAPLACE, es un modelo que se aplica para determinar el perfil de
temperaturas en una o dos direcciones, en coordenadas rectangulares, cilíndricas o
esféricas, en modelos isotrópicos, sin generación de energía calorífica interna, con
propiedades físico químicas constantes, en estado estacionario.
En coordenadas rectangulares la ecuación del Laplaciano bidimensional será la
siguiente:
2 2
2 ∂ T ∂T
∇ T= 2
+ 2 =0
∂x ∂ y
Aplicando diferencias finitas en el punto m, n en la cruceta, y considerando
diferencias centrales, se tiene:
] ]
2 2 2 2
∂T ∂ T ∂ T ∂ T
2
+ 2= 2 + 2
=0
∂x ∂ y ∂ x m, n ∂y m, n
T m +1 ,n−2 T m ,n +T m−1 ,n T m , n+1−2 T m ,n +T m ,n−1

¿ 2
+ 2
=0
∆x ∆y
Donde los intervalos ∆x y ∆y son diferentes.
Para resolver este sistema de ecuaciones algebraicas, deben conocerse las
condiciones de frontera (cuatro condiciones limite). Una aproximación a la solución
de esta ecuación diferencial parcial en dos dimensiones, es suponer que ∆x = ∆y.
Esto significa que la ecuación para los puntos interiores tendrá la forma.
0=T m+1 ,n +T m−1 ,n−4 T m ,n +T m−1 ,n +T m , n−1
5-7.- ECUACIONES DIFERENCIALES CON

CONDUCCION UNIDIMENSIONAL
Resolver una ecuación diferencial parcial u ordinaria por un método numérico,
significa discretizar la ecuación diferencial en nodos, aplicando una red o malla.
pág. 1
Una ecuación diferencial ordinaria de una dimensión (⟶ x), emplea solo nodos en
m desde 0 que es la condición inicial hasta un valor próximo a M -1 = L, aplicado a
placas planas de longitud L.
En ecuaciones diferenciales parciales de segundo orden, se aplican diferencias
numéricas centrales.
Generalmente la dificultad más importante es la interpretación de las condiciones
de frontera. Para nodos desde m = 1 hasta m = M – 2, se emplea la diferencia
central para los puntos internos.
Los problemas o modelos de conducción de calor, consideran ecuaciones
diferenciales de segundo orden, en el caso unidimensional puede aplicarse la
ecuación diferencial parcial.
∂
k ( )
∂T
∂x ∂x
+ ė gen=ρ c
∂T
∂t
(5-23)
En modelos con conductividad térmica constante, en estado estacionario o cuasi

estacionario, las ecuaciones a emplear serán:
ECUACION DE BIOT – FOURIER
∂ T ė gen 1 ∂T
2
+ = (5-24)
2
∂x k α ∂t
ECUACION DE POISSON
∂ T ė gen
2
+ =0 (5-25)
2
∂x k
ECUACION DE LA DIFUSION O SEGUNDA LEY DE FOURIER
2
∂ T 1 ∂T (5-26)
2
=
∂x α ∂t
ECUACION DE LAPLACE
2
d T (5-27)
2
=0
dx
5-8.- CONDICIONES DE FRONTERA

La aplicación de los nodos, en las condiciones de frontera, generan una ecuación
algebraica con valores fijos. Estas ecuaciones son diferentes que un nodo interior.
Para determinar la expresión en estos nodos debe emplearse un balance de
energía o la condición de frontera descrita en el capítulo 2. Además, en estos
puntos, se considera un volumen de control en ∆x/2.
Las condiciones de frontera principales que se aplican son las siguientes:
5-8.1.- CONDICIÓN DE FRONTERA INICIAL – x = 0.

En este punto en m = 0 Tm = T0. Se tiene los siguientes casos
pág. 2
A.- Temperatura especifica. – Significa un valor definido de temperatura uniforme
y constante
x=0 m=0 T (0) = T0
B.- Con flujo de calor q0 (w/m2), conducción de calor, y generación interna de
calor.
( T 1−T 0 ) ė o A ∆ x
q̇ o A +kA + =0
∆x 2
C.- Superficie aislada, conducción del calor y generación interna.
( T 1−T 0 ) ėo A ∆ x
kA + =0
∆x 2
D.- Con convección, conducción y generación interna.
( T 1−T 0 ) ė o A ∆ x
hA ( T ∞ −T 0 ) +kA + =0
∆x 2
E.- Con radiación, conducción, y generación interna
( T 1−T 0 ) ė o A ∆ x
εσA ( T alr −T 0 ) + kA
4 4
+ =0
∆x 2
F.- Con convección, radiación, conducción y generación interna.
( T 1−T 0 ) ė o A ∆ x
h A ( T ∞−T 0 ) + εσA ( T alr −T 0 ) + kA
4 4
+ =0
∆x 2
G.- Con convección, radiación, flujo de calor, conducción, y generación
interna
( T 1−T 0 ) ė o A ∆ x
q̇ 0 A+ h A ( T ∞−T 0 ) +εσA ( T alr −T 0 ) + kA
4 4
+ =0
∆x 2
5-8.2.- CONDICIÓN DE FRONTERA FINAL - x = L.

En estas condiciones, el número de nodos a emplear, determina el nodo final, pero
este se considera contando desde el cero. En una ecuación diferencial con 30
nodos, el nodo inicial está en m = 0, la condición inicial, y los nodos interiores son m
= 1,2,3,4 …… 28, son los nodos internos y el ultimo nodo 29 es la condición final.
Por lo tanto, en n=29 está el valor de T o valor superficial T L (29) y TL-1 (28) son los
dos valores que se deben tomar en cuenta, para aplicar en la condición final.
A.- Temperatura especifica. - Se tiene un valor definido
x=L m=L T(x) = T L
pág. 3
B.- Con flujo de calor q0, conducción, y generación interna
( T L−T L−1 ) ėo A ∆ x

q̇ o A−kA + =0
∆x 2
C.- Con frontera aislada, conducción, y generación interna
( T L−T L−1 ) ė o A ∆ x
−kA + =0
∆x 2
D.- Con convección, conducción y generación interna
( T 1−T 0 ) ėo A ∆ x
−kA + =h A ( T L −T ∞ )
∆x 2
E.- Con radiación, conducción y generación interna
( T 1−T 0 ) ėo A ∆ x
=εσA ( T L −T alr )
4 4
−kA +
∆x 2
E.- Con convección, radiación, conducción y generación interna
( T 1−T 0 ) ėo A ∆ x
=h A ( T L −T ∞ ) +εσA ( T L −T alr )
4 4
−kA +
∆x 2
F.- Con convección, radiación, flujo de calor, conducción, y generación
interna
( T 1−T 0 ) ėo A ∆ x
=h A ( T L −T ∞ ) +εσA ( T L −T alr )
4 4
q̇ 0 A−kA +
∆x 2
En todas estas ecuaciones debe eliminarse el término del área transversal A, todas
las ecuaciones con flujo de calor, convección, radiación deben ser acomodadas a
una expresión donde solo se considere la temperatura y los demás términos en
relación de las propiedades físicas de k, h o e0.
pág. 4

Capitulo N°5 Parte B

Cargado por

Copyright:

Formatos disponibles

Capitulo N°5 Parte B

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Capitulo N°5 Parte B

Cargado por

Copyright:

Formatos disponibles

CAPITULO 5

MÉTODOS NUMÉRICOS APLICADOS

T m +1 ,n−2 T m ,n +T m−1 ,n T m , n+1−2 T m ,n +T m ,n−1

0=T m+1 ,n +T m−1 ,n−4 T m ,n +T m−1 ,n +T m , n−1

5-7.- ECUACIONES DIFERENCIALES CON

En modelos con conductividad térmica constante, en estado estacionario o cuasi

5-8.- CONDICIONES DE FRONTERA

5-8.1.- CONDICIÓN DE FRONTERA INICIAL – x = 0.

5-8.2.- CONDICIÓN DE FRONTERA FINAL - x = L.

( T L−T L−1 ) ėo A ∆ x

También podría gustarte