Balance Diferencial de Energia
Balance Diferencial de Energia
Balance Diferencial de Energia
𝒒̇ 𝑮 es la razón volumétrica de generación de energía (W/m3), que es negativa cuando el sistema consume
𝑩.𝑬.
energía, y positiva cuando genera. Al dividir entre el volumen diferencial y tomar 𝐥𝐢𝐦 se tiene que:
∆𝑽→𝟎 ∆𝑽
𝑲 𝛁𝕋|𝒋 − 𝑲 𝛁𝕋|𝒋
∆𝒋 (𝝆 𝑪𝑷 ) ∆𝑽 𝝏𝕋 𝒒̇ 𝑮 ∆𝑽 𝝏𝕋
𝐥𝐢𝐦 = 𝐥𝐢𝐦 − 𝐥𝐢𝐦 ⟺ 𝛁 ∙ (𝑲 𝛁𝕋) = (𝝆 𝑪𝑷 ) − 𝒒̇ 𝑮
∆𝑽→𝟎 ∆𝒋 ∆𝑽→𝟎 ∆𝑽 𝝏𝒕 ∆𝑽→𝟎 ∆𝑽 𝝏𝒕
Al considerar que las propiedades no varían significativamente con la temperatura, y definir la difusividad
𝑲
térmica (α) como la relación entre la energía transferida y la energía acumulada 𝜶 = , se obtiene:
𝝆 𝑪𝑷
𝝏𝕋 𝝆 𝑪𝑷 𝝏𝕋 𝒒̇ 𝑮 𝟏 𝝏𝕋 𝒒̇ 𝑮
𝛁 ∙ (𝑲 𝛁𝕋) = 𝛁 ∙ 𝛁 (𝑲 𝕋) = 𝑲 𝛁 𝟐 𝕋 = (𝝆 𝑪𝑷 ) − 𝒒̇ 𝑮 ∴ 𝛁𝟐𝕋 = − = −
𝝏𝒕 𝑲 𝝏𝒕 𝑲 𝜶 𝝏𝒕 𝑲
𝟏 𝝏𝕋 𝒒̇ 𝑮
𝐵𝑎𝑙𝑎𝑛𝑐𝑒 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎𝑙 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔í𝑎: 𝛁𝟐𝕋 = −
𝜶 𝝏𝒕 𝑲
El balance diferencial para flujo estacionario unidireccional de energía térmica es una E.D.O.
Al resolver una ecuación diferencial, aparecen constantes de integración cuyo valor depende de las
condiciones de frontera. Con el fin de describir por completo un problema de intercambio térmico, se
deben establecer dos condiciones de frontera para cada dirección del sistema de coordenadas donde
ocurra un flujo significativo de energía térmica. Por lo tanto, se necesitan especificar dos condiciones de
frontera en problemas unidimensionales, cuatro para los bidimensionales y seis para los tridimensionales.
Esto es, se deben especificar tantas condiciones de frontera en una dirección como sea el orden de la
ecuación diferencial en esa dirección.
Los dispositivos que se colocan perpendiculares a una superficie, a fin de incrementar la velocidad del
intercambio térmico, se conocen como aletas, mientras que la superficie que posee estos dispositivos se
denomina superficie extendida. Las aletas se fabrican de materiales muy conductivos, y mejoran el
intercambio térmico por convección y radiación entre la superficie y el fluido adyacente, al exponer un
área más grande a la convección y a la radiación térmica.
Flujo bidireccional de energía a través de mecanismos diferentes
Se considera un sólido de longitud L y sección transversal constante A, perpendicular a una superficie que
está a una temperatura 𝑻𝟎 , y en contacto con un fluido de temperatura 𝑻 que posee un coeficiente
de película h. Al aplicar un balance de energía en un elemento diferencial de volumen ∆𝑉 = 𝐴 ∆𝑥 bajo
condiciones estacionarias 𝑬̇𝒂𝒄𝒖𝒎 = 𝟎 , y suponiendo que las propiedades no varían significativamente
con la temperatura (𝑲 = 𝒄𝒕𝒆), se tiene que:
𝒅𝑻(𝒙) 𝒅𝑻(𝒙)
𝑬̇𝒆 − 𝑬̇𝒔 = 𝑬̇𝒂𝒄𝒖𝒎 − 𝑬̇𝒈𝒆𝒏 ⟺ −𝑲𝑨 − −𝑲𝑨 + 𝒉𝑺 [𝑻(𝒙) − 𝑻 ] = −𝒒̇ 𝑮 ∆𝑽
𝒅𝒙 𝒙 𝒅𝒙 𝒙 ∆𝒙
El área de intercambio por conducción (A) es diferente a la de convección (𝑺 = 𝒑 ∆𝒙), donde p es el
𝑩.𝑬.
perímetro del elemento diferencial. Al dividir entre ∆𝒙 y tomar 𝐥𝐢𝐦 se tiene que:
∆𝒙→𝟎 ∆𝒙
𝒅𝑻(𝒙) 𝒅𝑻(𝒙)
𝑲𝑨 − 𝑲𝑨
𝒅𝒙 𝒙 ∆𝒙 𝒅𝒙 𝒙 𝒉 (𝒑 ∆𝒙) 𝒒̇ 𝑮 (𝑨 ∆𝒙)
𝐥𝐢𝐦 − 𝐥𝐢𝐦 [𝑻(𝒙) − 𝑻 ] = − 𝐥𝐢𝐦
∆𝒙→𝟎 ∆𝒙 ∆𝒙→𝟎 ∆𝒙 ∆𝑽→𝟎 ∆𝒙
𝒅 𝒅𝑻(𝒙) 𝒅𝟐 𝑻(𝒙) 𝒉 𝒑 𝒒̇ 𝑮
𝑲𝑨 − 𝒉 𝒑 [𝑻(𝒙) − 𝑻 ] = −𝒒̇ 𝑮 𝑨 ∴ 𝟐
− [𝑻(𝒙) − 𝑻 ] = −
𝒅𝒙 𝒅𝒙 𝒅𝒙 𝑲𝑨 𝑲
Al definir la variable de desviación 𝜽(𝒙) = 𝑻(𝒙) − 𝑻 , donde 𝒅[𝜽(𝒙)] = 𝒅[𝑻(𝒙)], se obtiene una
ecuación diferencial ordinaria lineal no homogénea de II orden, cuya solución viene dada por:
𝒉𝒑 𝒒̇ 𝑮
𝑫𝟐 − 𝜽(𝒙) = 𝑷(𝑫) 𝜽(𝒙) = − ⇒ 𝜽(𝒙) = 𝜽(𝒙)(𝒉𝒐𝒎𝒐𝒈é𝒏𝒆𝒂) + 𝜽(𝒙)(𝒑𝒂𝒓𝒕𝒊𝒄𝒖𝒍𝒂𝒓)
𝑲𝑨 𝑲
𝑷(𝑫)
𝒒̇ 𝑮 𝒉𝒑
𝜽(𝒙) = 𝑪𝟏 𝒆𝒎𝒙 + 𝑪𝟐 𝒆 𝒎𝒙
+ 𝑑𝑜𝑛𝑑𝑒 𝒎=
𝒎𝟐 𝑲 𝑲𝑨
𝒙=𝟎
Se aplica una condición de frontera en la base del sólido ⇒ 𝜽(𝟎) = 𝑻𝟎 − 𝑻 = 𝜽𝟎
𝑻(𝟎) = 𝑻𝟎
El error cometido al suponer que la temperatura en una aleta de espesor t sólo varía a lo largo de ella
(análisis unidireccional), y no en otras direcciones, es despreciable (menos de 1%) cuando (𝒉 𝒕⁄𝑲) < 𝟎, 𝟐.
La segunda condición de frontera se aplica en la punta del sólido (𝒙 = 𝑳), lo que da origen a tres casos
particulares, dependiendo de la consideración que se haga.
𝒙=𝑳→∞
Sólido infinitamente largo: ⇒ 𝜽(𝑳) = 𝑻 − 𝑻 = 𝟎
𝑻(𝑳) = 𝑻
𝒒̇ 𝑮
𝑃𝑒𝑟𝑓𝑖𝑙 𝑑𝑒 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑒𝑛 𝑒𝑙 𝑠ó𝑙𝑖𝑑𝑜: 𝑻(𝒙) = 𝑻 + (𝑻𝟎 − 𝑻 ) 𝒆 𝒎𝒙 + 𝟐 (𝟏 − 𝒆 𝒎𝒙 )
𝒎 𝑲
𝐼𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜 𝑡é𝑟𝑚𝑖𝑐𝑜 𝑒𝑛 𝑒𝑙 𝑠ó𝑙𝑖𝑑𝑜: ̇ (𝑻
𝑸𝒇 = 𝑲 𝑨 𝒎 𝟎 − 𝑻 ) 𝑠𝑖 𝒒̇ 𝑮 = 𝟎
𝒙=𝑳 𝒅𝜽(𝒙)
Sólido finito sin intercambio en la punta: 𝑸̇(𝑳) = −𝑲𝑨
𝒅𝑻(𝒙)
=𝟎 ⇒ =𝟎
𝒅𝒙 𝒙 𝑳
𝒅𝒙 𝒙 𝑳
𝒄𝒐𝒔𝒉[𝒎 (𝑳 − 𝒙)]
𝑃𝑒𝑟𝑓𝑖𝑙 𝑑𝑒 𝑇 𝑒𝑛 𝑒𝑙 𝑠ó𝑙𝑖𝑑𝑜: 𝑻(𝒙) = 𝑻 + (𝑻𝟎 − 𝑻 )
𝒄𝒐𝒔𝒉(𝒎𝑳) 𝑠𝑖 𝒒̇ 𝑮 = 𝟎
𝐼𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜 𝑡é𝑟𝑚𝑖𝑐𝑜: ̇
𝑸𝒇 = 𝑲 𝑨 𝒎 (𝑻𝟎 − 𝑻 ) 𝒕𝒈𝒉(𝒎𝑳)
𝒅𝑻(𝒙)
𝑸̇(𝑳) = −𝑲𝑨 = 𝒉𝑨 [𝑻(𝑳) − 𝑻 ]
𝒅𝒙 𝒙 𝑳
Sólido finito con intercambio en la punta: 𝒅𝜽(𝒙) 𝒉
= − 𝜽(𝑳)
𝒅𝒙 𝒙 𝑳 𝑲
𝒉
𝒄𝒐𝒔𝒉[𝒎 (𝑳 − 𝒙)] +
𝒔𝒆𝒏𝒉[𝒎 (𝑳 − 𝒙)]
𝑃𝑒𝑟𝑓𝑖𝑙 𝑑𝑒 𝑇 𝑒𝑛 𝑒𝑙 𝑠ó𝑙𝑖𝑑𝑜: 𝑻(𝒙) = 𝑻 + (𝑻𝟎 − 𝑻 ) 𝒎 𝑲
𝒉
𝒄𝒐𝒔𝒉(𝒎𝑳) + 𝒔𝒆𝒏𝒉(𝒎𝑳)
𝒎𝑲
𝒉
𝒔𝒆𝒏𝒉(𝒎𝑳) + 𝒄𝒐𝒔𝒉(𝒎𝑳)
𝐼𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜 𝑡é𝑟𝑚𝑖𝑐𝑜: 𝑸̇𝒇 = 𝑲 𝑨 𝒎 (𝑻𝟎 − 𝑻 ) 𝒎 𝑲 𝑠𝑖 𝒒̇ 𝑮 = 𝟎
𝒉
𝒄𝒐𝒔𝒉(𝒎𝑳) + 𝒔𝒆𝒏𝒉(𝒎𝑳)
𝒎𝑲
Al diseñar una aleta es importante establecer su longitud apropiada, una vez especificados el material y
su sección transversal. Se pudiera pensar que cuanto más larga sea la aleta, mayor es el área de contacto
con el fluido y, en consecuencia, el intercambio térmico aumenta; sin embargo, la variación exponencial
de la temperatura a lo largo de ella hace que alcance la del fluido a cierta longitud. Si 𝒎𝑳 ≈ 𝟓, la aleta
se puede considerar infinita, y su intercambio sólo disminuye 1% cuando la longitud se reduce a la mitad.
En la práctica, si 𝒎𝑳 ≈ 𝟏 la aleta intercambia más del 75% de la energía que transfiere una infinita.
Intercambio térmico en una superficie extendida
𝑸̇𝒑𝒂𝒓𝒆𝒅 = 𝒉 𝑨𝟎 (𝑻𝟎 − 𝑻 )
𝑸̇𝒈𝒍𝒐𝒃𝒂𝒍 = 𝑸̇𝒑𝒂𝒓𝒆𝒅 + 𝑸̇𝒂𝒍𝒆𝒕𝒂𝒔 𝑑𝑜𝑛𝑑𝑒
𝑸̇𝒂𝒍𝒆𝒕𝒂𝒔 = 𝑵𝒂 𝑸̇𝒇 = 𝒉 𝜼𝒇 𝑨𝒇 (𝑻𝟎 − 𝑻 )
𝑵𝒂 es la cantidad de aletas que posee la superficie, 𝑨𝟎 = 𝑨𝒑𝒂𝒓𝒆𝒅 − 𝑵𝒂 𝑨𝒃𝒂𝒔𝒆 es el área de la pared que
está en contacto directo con el fluido (área libre de aletas), y 𝑨𝒇 = 𝑵𝒂 𝑺 es el área de intercambio por
convección de todas las aletas (área total de las aletas). Al sustituir en la expresión global, se obtiene:
𝟏
𝑸̇𝒈𝒍𝒐𝒃𝒂𝒍 = 𝒉 𝑨𝟎 (𝑻𝟎 − 𝑻 ) + 𝒉 𝜼𝒇 𝑨𝒇 (𝑻𝟎 − 𝑻 ) = 𝒉 𝑨𝟎 + 𝜼𝒇 𝑨𝒇 (𝑻𝟎 − 𝑻 ) ∴ 𝑹𝒕 =
𝒉 𝑨𝒆𝒇
á𝒓𝒆𝒂 𝒆𝒇𝒆𝒄𝒕𝒊𝒗𝒂
Como puede observarse, la resistencia térmica convectiva del fluido que está en contacto con la superficie
extendida, disminuye al aumentar el área efectiva de contacto 𝑨𝒆𝒇 = 𝑨𝟎 + 𝜼𝒇 𝑨𝒇 , lo que incrementa
el intercambio térmico. Los parámetros de desempeño de la superficie extendida vienen dados por:
𝑸̇𝒈𝒍𝒐𝒃𝒂𝒍 𝒉 𝑨𝟎 + 𝜼𝒇 𝑨𝒇 (𝑻𝟎 − 𝑻 ) 𝑨𝟎 + 𝜼𝒇 𝑨𝒇
𝐸𝑓𝑒𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑 𝑔𝑙𝑜𝑏𝑎𝑙: 𝜺𝒈 = = = >𝟏
𝑸̇ (𝒔𝒊𝒏 𝒂𝒍𝒆𝒕𝒂𝒔) 𝒉 𝑨𝒑𝒂𝒓𝒆𝒅 (𝑻𝟎 − 𝑻 ) 𝑨𝒑𝒂𝒓𝒆𝒅
𝑸̇𝒈𝒍𝒐𝒃𝒂𝒍 𝒉 𝑨𝟎 + 𝜼𝒇 𝑨𝒇 (𝑻𝟎 − 𝑻 ) 𝑨𝟎 + 𝜼𝒇 𝑨𝒇
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑔𝑙𝑜𝑏𝑎𝑙: 𝜼𝒈 = = = <𝟏
𝑸̇𝑴Á𝑿 (𝒊𝒅𝒆𝒂𝒍) 𝒉 𝑨𝟎 + 𝑨𝒇 (𝑻𝟎 − 𝑻 ) 𝑨𝟎 + 𝑨𝒇
𝒉𝒑 𝒉 [𝟐 (𝒘 + 𝒕)] 𝒉 (𝟐 𝒘) 𝟐𝒉 𝟐𝒉
𝒎= = ≈ = ∴ 𝒎=
𝑲𝑨 𝑲 (𝒘 𝒕) 𝑲 (𝒘 𝒕) 𝑲𝒕 𝑲𝒕
Rectangular: tiene sección transversal y espesor constantes, y no termina en punta.
𝒕
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺 = 𝟐𝒘𝑳+𝟐𝒕𝑳+𝒘𝒕≈ 𝟐𝒘 𝑳+ = 𝟐 𝒘 𝑳𝑪
𝟐
𝒕𝒈𝒉(𝒎 𝑳𝑪 )
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 = (𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑎𝑙 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒)
𝒎 𝑳𝑪
Triangular: su espesor decrece linealmente con la longitud, hasta terminar en punta.
𝟐 𝟐
𝒕 𝒕𝑳 𝒕
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺 = 𝟐 𝒘 𝑳𝟐 + +𝟐 ≈ 𝟐 𝒘 𝑳𝟐 +
𝟐 𝟐 𝟐
𝟏 𝑰𝟏 (𝟐 𝒎 𝑳) 𝑰𝒏 (𝟐𝒎𝑳): 𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑚𝑜𝑑𝑖𝑓𝑖𝑐𝑎𝑑𝑎
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 =
𝒎 𝑳 𝑰𝟎 (𝟐 𝒎 𝑳) 𝑑𝑒 𝐵𝑒𝑠𝑠𝑒𝑙 𝑑𝑒 𝑝𝑟𝑖𝑚𝑒𝑟𝑎 𝑒𝑠𝑝𝑒𝑐𝑖𝑒
𝒙 𝟐
Parabólica: su espesor decrece según 𝒕(𝒙) = 𝒕 𝟏 − , hasta terminar en punta.
𝑳
𝟐 𝟐
𝒕 𝑳 𝒕 𝒕
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺=𝒘𝑳 𝟏+ + 𝒍𝒏 + 𝟏 +
𝑳 𝒕 𝑳 𝑳
𝟐
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 =
𝟏 + (𝟐 𝒎 𝑳)𝟐 + 𝟏
Aguja: su sección transversal es circular (𝑨𝒃𝒂𝒔𝒆 = 𝝅 𝑫), con diámetro D y longitud L; agrupa:
𝒉𝒑 𝒉 (𝝅 𝑫) 𝟒𝒉 𝟒𝒉
𝒎= = 𝝅 𝟐 = 𝑲𝑫 ∴ 𝒎= 𝑲𝑫
𝑲𝑨 𝑲 𝑫
𝟒
Cilíndrica: tiene sección transversal y diámetro constantes, y no termina en punta.
𝝅 𝑫
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺 = 𝝅 𝑫 𝑳 + 𝑫𝟐 = 𝝅 𝑫 𝑳 + = 𝝅 𝑫 𝑳𝑪
𝟒 𝟒
𝒕𝒈𝒉(𝒎 𝑳𝑪 )
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 = (𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑎𝑙 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒)
𝒎 𝑳𝑪
Cónica: su espesor decrece linealmente con la longitud, hasta terminar en punta.
𝟏 𝑫 𝟐
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺= 𝝅 𝑫 𝑳𝟐 +
𝟐 𝟐
𝟐 𝑰𝟐 (𝟐 𝒎 𝑳)
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 =
𝒎 𝑳 𝑰𝟏 (𝟐 𝒎 𝑳)
𝒙 𝟐
Aguja parabólica: su diámetro decrece según 𝑫(𝒙) = 𝑫 𝟏 − , hasta terminar en punta.
𝑳
𝟐 𝟐 𝟐
𝝅𝑳𝟑 𝑫 𝑫 𝑳 𝑫 𝟐𝑫 𝒕 𝟐
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜: 𝑺= 𝟏+𝟐 𝟏+ + 𝒍𝒏 𝟏 + 𝟐 + 𝟏+
𝟖𝑫 𝑳 𝑳 𝟐𝑫 𝑳 𝑳 𝑳
𝟐
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎: 𝜼𝒇 =
𝟐
𝟐
𝟏+ 𝒎𝑳 +𝟏
𝟑
Anular: tiene forma de arandela (𝑨𝒃𝒂𝒔𝒆 = 𝟐 𝝅 𝑹𝟎 𝒕), con radio interno 𝑹𝟎 , radio externo 𝑹𝑳 y
𝒉𝒑 𝒉 [𝟐(𝟐𝝅 𝒓)] 𝟐𝒉 𝟐𝒉
longitud 𝑳 = 𝑹𝑳 − 𝑹𝟎 ; no termina en punta. 𝒎= = = ∴ 𝒎=
𝑲𝑨 𝑲 (𝟐𝝅 𝒓 𝒕) 𝑲𝒕 𝑲𝒕
𝒕
Á𝑟𝑒𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑎𝑚𝑏𝑖𝑜 𝑺 = 𝟐𝝅 𝑹𝟐𝑳 − 𝑹𝟐𝟎 + 𝟐𝝅𝑹𝑳 𝒕 ≈ 𝟐𝝅 𝑹𝟐𝑳𝑪 − 𝑹𝟐𝟎 𝑐𝑜𝑛 𝑹𝑳𝑪 = 𝑹𝑳 +
𝟐
𝑻(𝒓) − 𝑻 𝑰𝟎 (𝒎𝒓) 𝑲𝟏 (𝒎𝑹𝑳 ) + 𝑲𝟎 (𝒎𝒓) 𝑰𝟏 (𝒎𝑹𝑳 )
𝑃𝑒𝑟𝑓𝑖𝑙 𝑑𝑒 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 =
𝑻𝟎 − 𝑻 𝑰𝟎 (𝒎𝑹𝟎 ) 𝑲𝟏 (𝒎𝑹𝑳 ) + 𝑲𝟎 (𝒎𝑹𝟎 ) 𝑰𝟏 (𝒎𝑹𝑳 )
𝟐 (𝑹𝟎 ⁄𝒎) 𝑲𝟏 (𝒎𝑹𝟎 ) 𝑰𝟏 (𝒎𝑹𝑳𝑪 ) − 𝑰𝟏 (𝒎𝑹𝟎 ) 𝑲𝟏 (𝒎𝑹𝑳𝑪 )
𝐸𝑓𝑖𝑐𝑖𝑒𝑛𝑐𝑖𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎 𝜼𝒇 =
𝑹𝟐𝑳𝑪 − 𝑹𝟐𝟎 𝑲𝟎 (𝒎𝑹𝟎 ) 𝑰𝟏 (𝒎𝑹𝑳𝑪 ) + 𝑰𝟎 (𝒎𝑹𝟎 ) 𝑲𝟏 (𝒎𝑹𝑳𝑪 )