Tarea Transferencia de Calor
Tarea Transferencia de Calor
Tarea Transferencia de Calor
Facultad de Ingeniería
E.A.P. Ingeniería En Energía
AUTOR:
Zavaleta Jara Gerson
2
DEDICATORIA
3
INTRODUCCIÓN
En la siguiente parte nos ocupamos de los problemas en los que la
temperatura y el flujo de calor puede ser tratado como funciones de una
sola variable. Muchos problemas prácticos entran en esta categoría,
pero cuando los límites de un sistema son irregulares o cuando la
temperatura a lo largo de un límite es no uniforme, un unidimensional
tratamiento ya no puede ser satisfactorio. En tales casos, la temperatura
es una función de dos o posiblemente incluso tres coordenadas. El flujo
de calor a través de una sección de la esquina donde dos o tres paredes
se unen, la conducción de calor a través de las paredes de un cilindro
corto y hueco, y la pérdida de calor desde una enterrados tubería son
ejemplos típicos de esta clase de problema.
Ahora vamos a considerar algunos métodos para el análisis de la
conducción en dos y sistemas tridimensionales. El énfasis será puesto
en problemas bidimensionales porque son menos engorroso para
resolver, sin embargo, que ilustran los métodos básicos de análisis para
sistemas tridimensionales. La conducción de calor en dos y tres
dimensiones sistemas pueden ser tratados por analítico, gráfico,
analógico, numérico y métodos computacionales. Para algunos casos,
"factores de forma" también están disponibles.
Nosotros tendremos en cuenta analítica, gráfica, y los métodos de factor
4
OBJETIVOS
Generales:
5
1.1. Solución Analítica
El objetivo de cualquier análisis de la transferencia de calor es predecir la
tasa de flujo de calor, la temperatura distribución, o ambos. De acuerdo
con la Ecuación,
𝛛𝟐 𝐓 𝛛𝟐 𝐓 𝛛𝟐 𝐓 …(1)
𝟐
+ 𝟐 + 𝟐 − 𝛁𝟐𝐓 = 𝟎
𝛛𝐱 𝛛𝐲 𝛛𝐳
En un sistema de dos dimensiones sin fuentes de calor la ecuación de
conducción general que rige la distribución de temperaturas en el estado
estacionario es:
𝛛𝟐 𝐓 𝛛𝟐 𝐓 …(2)
+ =𝟎
𝛛𝐱 𝟐 𝛛𝐲 𝟐
Si la conductividad térmica es uniforme. La solución de la ecuación (2)
dará T(x, y), la temperatura como una función de las dos coordenadas
espaciales x y y. Los componentes del flujo de calor por unidad de área o
flujo de calor 𝑞 ′′ en las direcciones x y y (𝑞𝑥′′ , y 𝑞𝑦′′ , respectivamente) se
pueden obtener a partir de la ley de Fourier:
𝑞 𝜕𝑇
𝑞𝑥′′ = ( ) = −𝑘
𝐴 𝑥 𝜕𝑥
𝑞 𝜕𝑇
𝑞𝑦′′ = ( ) = −𝑘
𝐴 𝑦 𝜕𝑦
Hay que señalar que mientras que la temperatura es un escalar, el flujo de
calor depende del gradiente de temperatura y por lo tanto es un vector. El
flujo de calor 𝑞 ′′ en un punto dado x, y es la resultante de las componentes
𝑞𝑥′′ y 𝑞𝑦′′ en ese punto y está dirigido directamente a la isoterma, como se
muestra en la Fig. 1. Por tanto, si se conoce la distribución de temperatura
en un sistema, la tasa de flujo de calor se puede calcular fácilmente. Por
tanto, los análisis de transferencia de calor suelen concentrarse en la
determinar el campo de temperatura.
6
de separación de variables. Este enfoque se ilustrará aplicándolo a un
problema relativamente simple.
𝑻 = 𝑿𝒀 …(3)
𝟏 𝒅𝟐 𝑿 𝟏 𝒅𝟐 𝒀
− = …(4)
𝑿 𝒅𝒙𝟐 𝒀 𝒅𝒚𝟐
Ahora las variables están separadas. El lado izquierdo es una función de
solo x, mientras que el lado derecho es una función solo de y. Dado que
ningún lado puede cambiar ya que x y y varían, ambos deben ser igual a
7
una constante, digamos 𝜆2 . Por lo tanto se llega a las dos ecuaciones
diferenciales
𝒅𝟐 𝑿
+ 𝝀𝟐 𝑿 = 𝟎 …(5)
𝒅𝒙𝟐
𝒅𝟐 𝒀
𝟐
− 𝝀𝟐 𝒀 = 𝟎 …(6)
𝒅𝒚
𝑋 = 𝐴 cos 𝜆𝑥 + 𝐵 sin 𝜆𝑥
La solución general de la ecuación (6) es
𝑌 = 𝐶𝑒 −𝜆𝑦 + 𝐷𝑒 𝜆𝑦
Por lo tanto, la ecuación (3),
8
𝑇 = 0 𝑒𝑛 𝑦=0
𝑇=0 𝑒𝑛 𝑥=0
𝑇=0 𝑒𝑛 𝑥=𝐿
𝑇 = 𝑇𝑚 sin(𝜋𝑥/𝐿) 𝑒𝑛 𝑦=𝑏
Sustituyendo estas condiciones en la ecuación (7) para T, de la primera
condición se obtiene
𝐴(𝐶𝑒 −𝜆𝑦 + 𝐷𝑒 𝜆𝑦 ) = 0
y de la tercera condición
9
𝒔𝒊𝒏𝒉(𝝅𝒚/𝑳) 𝝅𝒙 …(9)
𝑻(𝒙, 𝒚) = 𝑻𝒎 𝐬𝐢𝐧
𝒔𝒊𝒏𝒉(𝝅𝒃/𝑳) 𝑳
El campo de temperatura correspondiente se muestra en la Fig. 1.2. Las
líneas continuas son isotermas, y las líneas discontinuas son líneas de
flujo de calor. Cabe señalar que las líneas que indica la dirección del flujo
de calor son perpendiculares a las isotermas.
∂2 T ∂2 T ∂2 T
+ + =0
∂x 2 ∂y 2 ∂z 2
Separando las variables e integrando las ecuaciones diferenciales totales
resultantes sujetas a las condiciones de frontera dadas.
10
1.2. Método gráfico
El método gráfico presentado en esta sección puede producir rápidamente
una estimación razonablemente buena de la distribución de la temperatura
y el flujo de calor en sistemas bidimensionales geométricamente
complejos, pero su aplicación está limita a los problemas con fronteras
isotérmicas y aisladas. El objeto de una solución gráfica es construir una
red que consta de las isotermas (líneas de temperatura constante) y líneas
de flujo constante (líneas de flujo de calor constante). Las líneas de flujo
son análogas a las líneas de corriente en un flujo potencial, es decir, que
son tangentes a la dirección de flujo de calor en cualquier punto. En
consecuencia, el calor no puede fluir a través de las líneas de flujo
constante. Las isotermas son análogas a las líneas de potencial constante,
y los flujos de calor perpendiculares a ellas. Así pues, las líneas de
temperatura constante y las líneas de flujo de calor constante se
intersecan a ángulos rectos. Para obtener la distribución de temperatura
primero se elabora un modelo a escala y luego se trazan a mano alzada
las isotermas y las líneas de flujo, con el método de prueba y error, hasta
que formen una red de cuadrados curvilíneos.
11
FIGURA 1.3. Trazo de una red de cuadrados curvilíneos para
una sección en esquina: (a) modelo a escala; (b) trazo de flujo;
(c) cuadrado curvilíneo común.
∆𝑇
∆𝑞 = −𝑘(∆𝑙 × 1) = −𝑘∆𝑇
∆𝑙
𝑇2 − 𝑇1
∆𝑇 =
𝑁
12
de los cuadrados en una red de cuadrados curvilíneos. Por lo tanto, la tasa
total de transferencia de calor se puede escribir como
𝒏=𝑴
𝑴 𝑴
𝒒 = ∑ ∆𝒒𝒏 = 𝒌(𝑻𝟐 − 𝑻𝟏 ) = 𝒌∆𝑻𝒈𝒍𝒐𝒃𝒂𝒍 …(11)
𝑵 𝑵
𝒏=𝟏
𝒒 = 𝒌𝑺∆𝑻𝒈𝒍𝒐𝒃𝒂𝒍 …(12)
𝑐𝑜𝑠ℎ−1 𝑥 = ln(𝑥 ± √𝑥 2 − 1)
13
TABLA 1
Factor de forma conductivos
Descripción del
Bosquejo simbólico Factor de forma S
sistema
2𝜋𝐿
Cilindro isoterma de 𝐿≫𝑟
𝑐𝑜𝑠ℎ−1 (𝐷/𝑟)
radio r inmerso en un
medio semi-infinito
cuya superficie es 2𝜋𝐿
𝐿 ≫ 𝑟 ; 𝐷 > 3𝑟
isoterma 𝐷
𝐿𝑛 ( )
𝑟
Esfera isoterma de
radio r inmersa en un 4𝜋𝑟
medio infinito
Esfera isoterma de
radio r
inmersa en un medio
semi-infinito cuya 4𝜋𝑟
superficie isotérmica 1 − 𝑟/2𝐷
∆𝑻
= 𝑻𝒔𝒖𝒑 − 𝑻𝒄𝒂𝒎𝒑𝒐 𝒍𝒆𝒋𝒂𝒏𝒐
Cubo inmerso en un
medio infinito, lado L 8.24 𝐿
14
Cilindro isotermo de 2𝜋𝐿
radio r situado en un 𝐿𝑛(2𝐿/𝑟)
medio semi-infinito
como se muestra 𝐿 ≫ 2𝑟
Paralelepípedo
rectangular isotermo
𝑏 −0.59 𝑏 −0.078
inmerso en un medio 1.685𝐿 [log (1 + )] ( )
𝑎 𝑐
semi-infinito cuya
superficie es isoterma
𝐴
Pared plana 𝐿
Flujo de calor unidimensional
2𝜋𝐿
Cilindro hueco, longitud 𝐿𝑛(𝑟0 /𝑟𝑖 )
L
𝐿 ≫𝑟
4𝜋𝑟0 𝑟𝑖
Esfera hueca
𝑟0 − 𝑟𝑖
15
Disco delgado 4𝑟
horizontal inmerso en
8𝑟
un medio semi-infinito
cuya superficie es 4𝜋𝑟
isoterma 𝜋/2 − 𝑡𝑎𝑛−1 (𝑟/2𝐷)
Semiesfera inmersa en
un medio semi-infinito
2𝜋𝑟
∆𝑻 = 𝑻𝒆𝒔𝒇𝒆𝒓 − 𝑻𝒄𝒂𝒎𝒑𝒐 𝒍𝒆𝒋
Esfera isoterma
inmersa en un medio 4𝜋𝑟
semi-infinito cuya 1 + 𝑟/2𝐷
superficie está aislada
4𝜋𝑟2
Dos esferas isotermas 𝑟2 (𝑟1 /𝐷)4 2𝑟
[1 − 2 ] − 2
inmersas en un medio 𝑟1 1 − (𝑟2 /𝐷) 𝐷
infinito
D > 5𝑟𝑚𝑎𝑥
𝜋𝑊 𝑫=𝟎
𝐿𝑛(4𝑊/𝐿)
𝑾>𝑳
Placa rectangular
delgada de longitud L, 2𝜋𝑊 𝑫≫𝑾
inmersa en un medio 𝐿𝑛(4𝑊/𝐿) 𝑾>𝑳
semi-infinito cuya
superficie es isoterma 2𝜋𝑊 𝑾≫𝑳
𝐿𝑛(2𝜋𝐷/𝐿) 𝑫 > 𝟐𝑾
16
4𝜋𝑟
𝜋
Discos paralelos [ − 𝑡𝑎𝑛−1 (𝑟/𝐷)]
2
inmersos en un medio
infinito
𝐷 > 5𝑟
−1 (𝑟/𝐷)𝑒𝑛
𝑡𝑎𝑛 𝑟𝑎𝑑𝑖𝑎𝑛𝑒𝑠
2𝜋𝐿
𝑟 2 + 𝑟22 − 𝐷2
Cilindros excéntricos 𝑐𝑜𝑠ℎ−1 ( 1 )
de longitud L 2𝑟1 𝑟2
𝐿 ≫ 𝑟2
2𝜋𝐿
Cilindro centrado en un
prisma cuadrado de
𝐿𝑛(0.54𝑊/𝑟)
longitud L
𝐿≫𝑊
17
Para una pared tridimensional, como la de un horno, se utilizan por
separados factores de forma para calcular el flujo de calor a través de las
secciones de las aristas y las esquinas. Si todas las dimensiones
interiores son mayores que un quinto de espesor de la pared
𝐴
𝑆𝑝𝑎𝑟𝑒𝑑 = 𝑆𝑎𝑟𝑖𝑠𝑡𝑎 = 0.54𝐷 𝑆𝑒𝑠𝑞𝑢𝑖𝑛𝑎 = 0.15𝐿
𝐿
Donde: 𝐴 = á𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑝𝑎𝑟𝑒𝑑
𝐿 = 𝑒𝑠𝑝𝑒𝑠𝑜𝑟 𝑑𝑒 𝑙𝑎 𝑝𝑎𝑟𝑒𝑑
𝐷 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑎𝑟𝑖𝑠𝑡𝑎
18
1.4. Método de análisis matemático
Durante los últimos 100 años se ha acumulado un inmerso número de
soluciones analíticas para los problemas de transferencia de calor. Aun
así, en numerosas situaciones prácticas, la geometría o condiciones de
frontera son tales que no se ha podido obtener una solución analítica, o si
se ha desarrollado una solución esta involucra una solución en serie tan
compleja que la evaluación numérica se vuelve en extremo difícil. Para
tales situaciones el enfoque más fructífero es el que se basa en técnicas
de diferencia finita, cuyos principios básicos presentaremos brevemente
en esta sección. Claro está que el rápido desarrollo de las computadoras
de alta velocidad ha permitido al especialista práctico en transferencia de
calor, obtener soluciones numéricas a muchos problemas que hasta el
momento se pensaban imposibles.
Considérese un cuerpo bidimensional que se ha de dividir en incrementos
iguales tanto en la dirección x como en la y, tal como se muestra en la
figura 1.5. Los puntos nodales están designados como se muestran; las
posiciones m indican el incremento x y las n indican el incremento y.
Deseamos establecer las temperaturas en cualquiera de los puntos
nodales dentro del cuerpo, utilizando la ecuación (2) como condición
dominante. Se usan diferencias finitas para aproximar incrementos
diferenciales en las coordenadas de temperatura y espacio; y entre más
pequeños escojamos estos incrementos finitos, más exacta será la
aproximación de la verdadera distribución de temperatura.
19
𝜕𝑇 𝑇𝑚+1,𝑛 − 𝑇𝑚,𝑛
] ≈
𝜕𝑥 𝑚+1/2,𝑛 ∆𝑥
𝜕𝑇 𝑇𝑚,𝑛 − 𝑇𝑚−1,𝑛
] ≈
𝜕𝑥 𝑚−1/2,𝑛 ∆𝑥
𝜕𝑇 𝑇𝑚,𝑛+1 − 𝑇𝑚,𝑛
] ≈
𝜕𝑥 𝑚,𝑛+1/2 ∆𝑦
𝜕𝑇 𝑇𝑚,𝑛 − 𝑇𝑚,𝑛−1
] ≈
𝜕𝑦 𝑚,,𝑛−1/2 ∆𝑦
𝜕𝑇 𝑇𝑚+1,𝑛 − 𝑇𝑚,𝑛
] ≈
𝜕𝑦 𝑚+1/2,𝑛 ∆𝑦
𝜕𝑇 𝜕𝑇
] − ]
2
𝜕 𝑇 𝜕𝑥 𝑚+1/2,𝑛 𝜕𝑥 𝑚−1/2,𝑛 𝑇𝑚+1,𝑛 + 𝑇𝑚−1,𝑛 − 2𝑇𝑚,𝑛
] ≈ =
𝜕𝑥 2 𝑚,𝑛 ∆𝑥 (∆𝑥)2
𝜕𝑇 𝜕𝑇
] − ]
2
𝜕 𝑇 𝜕𝑦 𝑚,𝑛+1/2 𝜕𝑦 𝑚,𝑛−1/2 𝑇𝑚,𝑛+1 + 𝑇𝑚,𝑛−1 − 2𝑇𝑚,𝑛
] ≈ =
𝜕𝑦 2 𝑚,𝑛 ∆𝑦 (∆𝑦)2
20
También puede diseñarse un esquema de diferencias finitas para tener en
cuenta la generación de calor. Simplemente se añade el término en la
ecuación general, obteniendo
𝑇𝑚+1,𝑛 + 𝑇𝑚−1,𝑛 − 2𝑇𝑚,𝑛 𝑇𝑚,𝑛+1 + 𝑇𝑚,𝑛−1 − 2𝑇𝑚,𝑛 𝑞̇
+ + =0
(∆𝑥)2 (∆𝑦)2 𝑘
Entonces para una malla cuadrada en la que ∆𝑥 = ∆𝑦,
𝑞̇ (∆𝑥)2
𝑻𝒎+𝟏,𝒏 + 𝑻𝒎−𝟏,𝒏 + 𝑻𝒎,𝒏+𝟏 + +𝑻𝒎,𝒏−𝟏 + − 𝟒𝑻𝒎,𝒏 = 𝟎 …(14)
𝑘
21
Cuando el sólido está sometido a condiciones de contorno conectivas, la
temperatura de la superficie debe calcularse de forma distinta al método
dado anteriormente. Considérese el contorno mostrado en la Figura 1.7.
El balance energético en ese nodo (𝑚, 𝑛) es
𝑇𝑚,𝑛 − 𝑇𝑚−1,𝑛 ∆𝑥 𝑇𝑚,𝑛 − 𝑇𝑚,𝑛+1 ∆𝑥 𝑇𝑚,𝑛 − 𝑇𝑚,𝑛−1
−𝑘∆𝑦 −𝑘 −𝑘 = ℎ∆𝑦(𝑇𝑚,𝑛 − 𝑇∞ )
∆𝑥 2 ∆𝑦 2 ∆𝑦
𝒉∆𝒙 𝒉∆𝒙 𝟏
𝑻𝒎,𝒏 ( + 𝟐) − 𝑻∞ − (𝟐𝑻𝒎−𝟏,𝒏 + 𝑻𝒎,𝒏+𝟏 + 𝑻𝒎,𝒏−𝟏 ) = 𝟎 …(15)
𝒌 𝒌 𝟐
22
∆𝑦 𝑇𝑚,𝑛 − 𝑇𝑚−1,𝑛 ∆𝑥 𝑇𝑚,𝑛 − 𝑇𝑚,𝑛−1 ∆𝑥 ∆𝑦
−𝑘 −𝑘 = ℎ (𝑇𝑚,𝑛 − 𝑇∞ ) + ℎ (𝑇𝑚,𝑛 − 𝑇𝑥 )
2 ∆𝑥 2 ∆𝑦 2 2
Si ∆𝑥 = ∆𝑦,
𝒉∆𝒙 𝒉∆𝒙
𝟐𝑻𝒎,𝒏 ( + 𝟏) − 𝟐 𝑻 − (𝑻𝒎−𝟏,𝒏 + 𝑻𝒎,𝒏−𝟏 ) = 𝟎 …(16)
𝒌 𝒌 𝒙
23
TABLA 2
ℎ∆𝑥
𝐵𝑖 =
𝑘
(c) Esquina exterior con contorno
convectivo ℎ∆𝑥 ℎ∆𝑥
0=2 𝑇∞ + (𝑇𝑚−1,𝑛 + 𝑇𝑚,𝑛−1 ) − 2 ( + 1) 𝑇𝑚,𝑛
𝑘 𝑘
ℎ∆𝑥
𝐵𝑖 =
𝑘
(d) Esquina interior con contorno
convectivo ℎ∆𝑥 ℎ∆𝑥
0=2 𝑇∞ + 2𝑇𝑚−1,𝑛 + 𝑇𝑚,𝑛+1 + 𝑇𝑚+1,𝑛 + 𝑇𝑚,𝑛−1 − 2 (3 + ) 𝑇𝑚,𝑛
𝑘 𝑘
ℎ∆𝑥
𝐵𝑖 =
𝑘
25
(f) Nodo interior cercano a un
contorno curvo
2 2 2 2
0= 𝑇2 + 𝑇𝑚+1,𝑛 + 𝑇𝑚,𝑛−1 + 𝑇
𝑏(𝑏 + 1) 𝑎+1 𝑏+1 𝑎(𝑎 + 1) 1
1 1
− 2 ( + ) 𝑇𝑚,𝑛
𝑎 𝑏
𝑏 𝑏 𝑎+1 ℎ∆𝑥
(g) Nodo de contorno con 0= 𝑇1 + 𝑇3 + 𝑇𝑚,𝑛 + (√𝑐 2 + 1 + √𝑎2 + 𝑏 2 ) 𝑇∞
2
√𝑎 + 𝑏 2 2
√𝑐 + 1 𝑏 𝑘
convección de un contorno 𝑏 𝑏 𝑎+1 ℎ∆𝑥
curvo-nodo 2 de (f) anterior −[ + + + (√𝑐 2 + 1 + √𝑎2 + 𝑏 2 ) ] 𝑇2
√𝑎2 + 𝑏 2 √𝑐 2 + 1 𝑏 𝑘
26
1.5. Formulación numérica en
términos de elementos resistivos
Hasta aquí hemos mostrado cómo los problemas de conducción pueden ser
resueltos por diferencias finitas aproximaciones a las ecuaciones
diferenciales. Una ecuación se formula para cada nodo y el conjunto de
ecuaciones que resuelve para las temperaturas en todo el cuerpo. En la
formulación de las ecuaciones que podrían muy bien haber utilizado un
concepto de resistencia para la escritura de la transferencia de calor entre
los nodos. La designación de nuestro nodo de interés con el subíndice i y
los nodos adyacentes con el subíndice j, tenemos la situación de conducción
de nodos general que se muestra en la Figura 1.9.
𝑻𝒋 − 𝑻𝒊
𝒒𝒊 + ∑ =𝟎 …(17)
𝑹𝒊𝒋
𝒋
A estado estacionario la entrada de calor neto para el nodo i debe ser cero
o donde 𝑞𝑖 se entrega el calor al nodo i por la generación de calor, radiación,
etc. El 𝑅𝑖𝑗 puede tomar la forma de los límites de convección, conducción
interna, etc., y la ecuación (17) puede igualarse a un residuo de una solución
de relajación o a cero para el tratamiento con matriz y métodos iterativos.
TABLA 3
Resistencias para los nodos
Configuración física 𝑹𝒎+ 𝑹𝒎− 𝑹𝒏+ 𝑹𝒎− ∆𝑽
1 1 1 1
(a) Nodo interior (∆𝑥)2
𝑘 𝑘 𝑘 𝑘
1 1 2 2 (∆𝑥)2
(b) Contorno convectivo
ℎ∆𝑥 𝑘 𝑘 𝑘 2
(c) Esquina exterior, 2 2 2 2 (∆𝑥)2
convección ℎ∆𝑥 𝑘 ℎ∆𝑥 𝑘 4
(d) Esquina interior, 2 1 1 2 3(∆𝑥)2
convección 𝑘 𝑘 𝑘 𝑘 4
1 2 2 (∆𝑥)2
(e) Contorno aislado ∞
𝑘 𝑘 𝑘 2
2 2
2𝑎 2𝑏
(f) Nodo interior cercano (𝑏 + 1)𝑘 (𝑎 + 1)𝑘
a un contorno curvo (𝑏 + 1)𝑘 (𝑎 + 1)𝑘 0.25(1 + 𝑎)(1 + 𝑏)(∆𝑥)2
𝑎𝑙 𝑛𝑜𝑑𝑜 𝑎𝑙 𝑛𝑜𝑑𝑜
𝑎𝑙 𝑛𝑜𝑑𝑜 1 𝑎𝑙 𝑛𝑜𝑑𝑜 2
(𝑚 + 1, 𝑛) (𝑚, 𝑛 − 1)
28
2√𝑐 2 + 1
𝑅23 =
𝑏𝑘
2√𝑎2 + 𝑏 2
𝑅23 =
(g) Nodo de contorno con 𝑏𝑘
contorno curva nodo 2 ∆𝑉 = 0.125[(2 + 𝑎) + 𝑐](∆𝑥)2
de (f) anterior 2
𝑅2−∞ =
ℎ∆𝑥(√𝑐 2 + 1 + √𝑎2 + 𝑏 2 )
2𝑏
𝑅𝑛− = 𝑎𝑙 𝑛𝑜𝑑𝑜 (𝑚, 𝑛)
𝑘(𝑎 + 1)
𝒒𝒊 + ∑𝒋(𝑻𝒋 /𝑹𝒊𝒋 )
𝑻𝒊 = …(18)
∑𝒋(𝟏/𝑹𝒊𝒋 )
29
3. El proceso se repite hasta que los cálculos sucesivos difieren por una
cantidad suficientemente pequeña. En términos de un programa de
ordenador, esto significa que una prueba puede ser insertado para
detener la cálculos cuando
𝑇𝑖𝑛+1 − 𝑇𝑖𝑛
𝜖≥| |
𝑇𝑖𝑛
Obviamente, cuanto menor es el valor de 𝛿, mayor es el tiempo de
cálculo necesario para obtener el resultado deseado. El lector debe
observar, sin embargo, que la precisión de la solución del problema
físico no depende del valor de 𝛿 solo. Esta constante gobierna la
precisión de la solución para el conjunto de ecuaciones en diferencias.
La solución a la física problema también depende de la selección del
incremento ∆𝑥.
30
TABLA 4
Resistencias internas de los nodos para distintas sistemas de coordenadas
Cartesianas Cilíndricas Esféricas
𝑥, 𝑚 𝑟, 𝑚 𝑟, 𝑚
Nomenclatura para los incrementos 𝑦, 𝑛 ∅, 𝑚 ∅, 𝑚
𝑧, 𝑘 𝑧, 𝑚 𝜃, 𝑚
∆𝑥 ∆𝑟 ∆𝑟
𝑅𝑚+
∆𝑦∆𝑧𝑘 (𝑟𝑚 + ∆𝑟/2)∆∅∆𝑧𝑘 (𝑟𝑚 + ∆𝑟/2)2 sin 𝜃 ∆∅∆𝜃𝑘
∆𝑥 ∆𝑟 ∆𝑟
𝑅𝑚−
∆𝑦∆𝑧𝑘 (𝑟𝑚 − ∆𝑟/2)∆∅∆𝑧𝑘 (𝑟𝑚 − ∆𝑟/2)2 sin 𝜃 ∆∅∆𝜃𝑘
∆𝑦 𝑟𝑚 ∆∅ ∆∅ sin 𝜃
𝑅𝑛+
∆𝑥∆𝑧𝑘 ∆𝑟∆𝑧𝑘 ∆𝑟∆𝜃𝑘
∆𝑦 𝑟𝑚 ∆∅ ∆∅ sin 𝜃
𝑅𝑛−
∆𝑥∆𝑧𝑘 ∆𝑟∆𝑧𝑘 ∆𝑟∆𝜃𝑘
∆𝑧 ∆𝑧 ∆𝜃
𝑅𝑘+
∆𝑥∆𝑦𝑘 𝑟𝑚 ∆𝑟∆∅𝑘 sin(𝜃 + 𝜃/2)∆𝑟∆∅𝑘
∆𝑧 ∆𝑧 ∆𝜃
𝑅𝑘−
∆𝑥∆𝑦𝑘 𝑟𝑚 ∆𝑟∆∅𝑘 sin(𝜃 − 𝜃/2)∆𝑟∆∅𝑘
CONCLUSIONES
Existen diversas formas de solución para paredes multidimensionales,
sin importar su forma cuya que pueden ser cartesianas, cilíndricas o
esféricas, y para poder aplicarlos tenemos que realizar el método
analítico como grafico adecuadamente.
BIBLIOGRAFÍA
Jack P. Holman. Transferencia de Calor. Décima edición. Capitulo 3.