Ley de enfriamiento de Newton para la condición de frontera de la ecuación del calor

La ley de enfriamiento de Newton dice que la temperatura de un objeto satisface

(1) d T d t = k ( T ( t ) T 0 ) ,

donde T 0 es la temperatura ambiente. Consulte estas notas HTML, por ejemplo.

Ahora si tu ( X , t ) denota la temperatura de una barra aislada lateralmente en un punto X y tiempo t , entonces la ecuación del calor dice

tu t = α 2 tu X 2

Si el extremo izquierdo X = 0 se expone a un ambiente a temperatura T 0 , todos los libros de física matemática dicen que la ley de enfriamiento de Newton es

(2) C tu X ( 0 , t ) = k ( tu ( 0 , t ) T 0 ) ,

donde C es una constante Consulte el Problema 5 en esta vista previa del libro de Google de Problemas de valor límite: y Ecuaciones diferenciales parciales de David Powers . Véase también la página 131 del mismo libro. (Tres páginas antes del problema 5, donde aparece por primera vez la ley de enfriamiento de Newton).

Sin embargo, de acuerdo con la ley de enfriamiento de Newton (1), obtenemos

tu t ( 0 , t ) = k ( tu ( 0 , t ) T 0 ) ,

Por la ecuación del calor, tu t = α 2 tu X 2 , concluimos que

(3) α 2 tu X 2 ( 0 , t ) = k ( tu ( 0 , t ) T 0 ) .

Esta condición de contorno no es la misma que (2). ¿Por qué (3) no es correcto?

A diferencia de su ecuación (2), su ecuación (3) no es una condición límite.
@AmeyJoshi ¿Por qué? Por definición, una "condición de contorno" es una condición en el contorno requerido de la función. Según esta definición, (3) es ciertamente una "condición de contorno".
Parte del problema aquí es su fuente para (1). No está mal per se, pero utiliza un modelo de transferencia de calor muy simplificado. En ese caso, k no se usa conductividad en (2), es más como el coeficiente de transferencia de calor, C . Está tratando el objeto como agrupado (temperatura constante) en lugar de resolver la transferencia de calor dentro de él, como lo hacen las ecuaciones posteriores que ha utilizado. Por eso también han usado la notación para ODE ( d / d t ) en lugar de PDE ( / t ).
@ user3823992, debería elaborar esto un poco en una respuesta; será más preciso que la respuesta existente, creo.

Respuestas (1)

La ley de enfriamiento de Newton en realidad proviene de la ecuación más general para el calor q transferida entre un sistema (temperatura T ) y su entorno (temperatura T 0 ):

d q d t = h A ( T T 0 )
donde A es el área a través de la cual se produce la transferencia de calor (ver, por ejemplo, aquí ). Para un objeto macroscópico ordinario, donde d q = metro C   d T , obtenemos la ley convencional de enfriamiento de Newton en términos de temperatura:
metro C d T d t = h A ( T T 0 )
Sin embargo, para el caso de la barra conductora, de la ley de Fourier:
1 A d q d t = k tu X
La condición de contorno es por lo tanto:
k tu ( 0 , t ) X = h ( tu ( 0 , t ) T 0 )

Pero, esto todavía no responde "por qué" --- ¿por qué no puedes usar la fórmula? d q = metro C d T para el caso de la barra para concluir la ecuación (3)? Es decir, suponga que no conoce la ley de Fourier. ¿Qué le impediría usar la fórmula? d q = metro C d T para deducir la fórmula tu / t ( 0 , t ) = k ( tu ( 0 , t ) T 0 ) ? ¿No creerías que esta fórmula es válida si no te dijeran que usaras la ley de Fourier en su lugar?
En el límite, hay un calor entrante (por conducción) q 1 y la transferencia de calor al entorno q 2 por unidad de tiempo. Al considerar un segmento de longitud d X en el límite, vemos que el calor neto que ingresa a la región es de hecho d q = q 1 q 2 = q ( d X ) q ( 0 ) . Es esto lo que contribuye al cambio de temperatura, no solo la pérdida de calor hacia el entorno. Este cambio está dado por q X d X tiempo d q = metro C T ˙ = ρ A C d X T ˙ . Así, obtenemos T ˙ = 1 ρ A C q X . Por lo tanto, es equivalente a la ecuación del calor...
...con la pieza faltante de lo que q en realidad está dada por la ley de Fourier.