Ley de enfriamiento de Newton

Cruz Reyes Jonathan

ESIME Azcapotzalco, Instituto Politécnico Nacional,
Av. De las granjas 682, Col. Sta. Catarina,
Azcapotzalco 02250 CDMX, Fecha: 07/Noviembre/22

El siguiente documento de trabajo tiene como objetivo, aplicar los

temas aprendidos en ecuaciones diferenciales, como métodos de
resolución de ecuaciones diferenciales, ejemplo de aplicación de
ecuaciones diferenciales, ley de enfriamiento de Newton, así como
un experimento para comprobar la ley.

Palabras clave: ecuaciones diferenciales, métodos de resolución de

ecuaciones diferenciales, experimento, ley de enfriamiento de Newton.

I. INTRODUCCIÓN La forma general de las ecuaciones

diferenciales ordinarias es:
A. ECUACIÓN DIFERENCIAL
𝐹(𝑥, 𝑓(𝑥), 𝑓′(𝑥), 𝑓′′(𝑥), . . . , 𝑓 𝑛 (𝑥)) = 0
Definimos una Ecuación diferencial
como una ecuación que relaciona una
función y su variable (o sus variables) El orden de una ecuación diferencial
con sus derivadas. está dado por el grado mayor de las
Si una ecuación diferencial contiene derivadas que aparecen en la
derivadas respecto a una sola variable ecuación. Así pues, de los ejemplos
independiente se llama ecuación mostrados anteriormente, la ecuación
diferencial ordinaria. Si contiene (1) es una ecuación diferencial de
derivadas parciales respecto a dos o primer orden, la ecuación (2) es una
más variables independientes, se ecuación diferencial de tercer orden, y
llama ecuación en derivadas. Las finalmente las ecuaciones (3) y (4) son
siguientes ecuaciones son ejemplos ecuaciones diferenciales de segundo
de ecuaciones diferenciales ordinarias orden.
que llamaremos simplemente
ecuaciones diferenciales:
Las ecuaciones diferenciales
𝑑𝑦 aparecen de forma natural en las
1. = 3𝑦 2 sin(𝑡 + 𝑦)
𝑑𝑡 ciencias y en las humanidades. En
𝑑3 𝑦 (−𝑦) 𝑑2 𝑦 esta unidad estudiaremos las
2. = 𝑒 + 𝑡 + ecuaciones diferenciales de primer
𝑑𝑡 3 𝑑𝑡 2
orden, tomando en cuenta que la
3. 𝑓 ′′ (𝑥) − 5 = 0
ecuación está reducida y puede
4. 𝑓 ′′ (𝑥) + 6𝑓 ′ (𝑥) + 3𝑓(𝑥) + 2 = 0 llevarse a la forma:
 Donde 𝑇0 es la temperatura del medio
ambiente y 𝑘 es la constante de
𝑑𝑦 proporcionalidad positiva.
= 𝑓(𝑥, 𝑦)
𝑑𝑥
Dado que la incógnita es 𝑇 y es el
argumento de una derivada, se trata
En forma general se tiene: de una ecuación diferencial. Esto
significa que se debe buscar una
función 𝑇(𝑡) tal que, al ser derivada
una vez, resulte en la expresión
𝐹(𝑥, 𝑓(𝑥), 𝑓′(𝑥)) = 0
𝐾(𝑇 − 𝑇0 )

B. LEY DE ENFRIAMIENTO DE C. SOLUCIÓN DE LA E.D.

NEWTON
Una forma de resolver la ecuación
La ley de enfriamiento de Newton diferencial que representa la ley de
asegura que la tasa de cambio en el enfriamiento de Newton es mediante
tiempo de la temperatura de un el método de separación de variables.
cuerpo resulta directamente Consiste en manipular
proporcional a la diferencia entre la algebraicamente la ecuación, para
temperatura medida en el entorno y en dejar a un lado de la igualdad todo lo
el cuerpo. que contenga una de las variables, y
al otro lado de la igualdad, quedaría la
Tal como dicta la experiencia otra variable:
cotidiana, un objeto a cierta
temperatura mayor o menor que la del 𝑑𝑇
= 𝑘(𝑇 − 𝑇0 )
medio ambiente, tiende con el tiempo 𝑑𝑡
a establecer un equilibrio térmico con
este último. Por ejemplo, una taza de Resolviendo tenemos que:
café caliente que se queda un buen
rato sobre la mesa, finalmente se 𝑑𝑇
[ = 𝑘(𝑇 − 𝑇0 )] [𝑑𝑡]
enfría hasta que su temperatura es la 𝑑𝑡
misma que la de la habitación
𝑑𝑇 = 𝑘(𝑇 − 𝑇0 )𝑑𝑡
Este fenómeno se modela
matemáticamente a través de la ley de 1
[𝑑𝑇 = 𝑘(𝑇 − 𝑇0 )𝑑𝑡] [ ]
enfriamiento de Newton, recordando (𝑇 − 𝑇0 )
que una tasa de cambio se puede
representar como una derivada 𝑑𝑇
= 𝑘𝑑𝑡
primera. En tal caso, la ley queda (𝑇 − 𝑇0 )
expresada así:
𝑑𝑇
𝑑𝑇 ∫ = 𝑘 ∫ 𝑑𝑡
= 𝑘(𝑇 − 𝑇0 ) (𝑇 − 𝑇0 )
𝑑𝑡
𝑙𝑛|𝑇 − 𝑇0 | = 𝑘 ∙ 𝑡 + 𝐶
Despejando a 𝑇 fuimos tomando las lecturas en
intervalos de tiempos iguales,
𝑒 𝑙𝑛|𝑇−𝑇0 | = 𝑒 𝑘∙𝑡 + 𝐶
A continuación, los datos recopilados:

𝑇 − 𝑇0 = 𝑒𝑘∙𝑡 + 𝐶 Tabla de lecturas recopiladas a

temperatura ambiente de 24.5°𝑪
Tiempo en Temperatura en °C
𝑇 − 𝑇0 = 𝑒𝑘∙𝑡 ∙ 𝑒𝐶 minutos
0 86.5
𝑇 − 𝑇0 = 𝑒𝑘∙𝑡 ∙𝐶 5 80.4
𝑘∙𝑡 10 75
𝑇 − 𝑇0 = 𝐶𝑒 15 70.5
20 66.6
𝑘∙𝑡
𝑇 = 𝐶𝑒 + 𝑇0 25 63.2
30 60.4
Solución general: 35 57.6
40 55.4
𝑘∙𝑡
𝑇(𝑡) = 𝐶𝑒 + 𝑇0 45 53.2
50 51.3
𝑇(𝑡): 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑜𝑏𝑗𝑒𝑡𝑜 𝑒𝑛 𝑡 55 50
60 48.4
𝑇0 : 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎𝐴𝑚𝑏𝑖𝑒𝑛𝑡𝑒(𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒) 65 46.8
70 45.5
𝐶: 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
75 44.3
𝑘: 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑝𝑟𝑜𝑝𝑜𝑟𝑐𝑖𝑜𝑛𝑎𝑙𝑖𝑑𝑎𝑑 80 43.2

𝑡: 𝑡𝑖𝑒𝑚𝑝𝑜 Recopilamos estos datos ya que no

contábamos con mucho tiempo para
medir, pero al final después de horas
el agua llegó a los 24.5°𝐶
II. DESARROLLO
Con los siguientes datos recopilados
Comenzamos por colocar 1 𝐿𝑡 de podemos realizar la comprobación con
agua en un recipiente, una vez en el la ley de enfriamiento de Newton, para
recipiente la colocamos en la estufa obtener la solución particular de la
para tratar de llegar al punto de ecuación:
ebullición, una vez cuando el agua
𝑘∙𝑡
estaba hirviendo, apagamos la estufa 𝑇(𝑡) = 𝐶𝑒 + 𝑇0
y comenzamos a medir, la primera
lectura fue de 86.5°𝐶 y al final después
de horas el agua se enfrío hasta llegar
a la temperatura ambiente de 24.5°𝐶,
Tenemos entonces:
18.7
= 𝑒𝑘(80)
Primero calcularemos la constante 𝐶. 62
Para esto necesitamos establecer
condiciones iniciales para sustituir en Eliminamos la exponencial
nuestra ecuación general y obtener
así la constante particular. 18.7
[ = 𝑒 𝑘(80) ] [𝑙𝑛]
Tomando los datos iniciales de 62
nuestra tabla:
18.7
𝑡 = 0 , 𝑇(𝑡) = 86.5°𝐶 𝑦 𝑐𝑜𝑛 𝑇0 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 = 24.5°𝐶 𝑙𝑛 = 𝑘 ∙ 80
62
Sustituyendo:
18.7
86.5 = 𝐶𝑒 𝑘∙0
+ 24.5 𝑙𝑛
62
=𝑘
Despejamos a 𝐶
80

86.5 − 24.5 = 𝐶𝑒 𝑘∙0 𝑘 = −0.0149826

62 = 𝐶𝑒 0 Ahora que hemos obtenido las

constantes, podemos escribir nuestra
𝐶 = 62 ecuación particular para el
enfriamiento de nuestro sistema,
Ahora tomaremos otras condiciones entonces tenemos que:
de nuestra tabla y vamos a sustituir
igualmente en la ecuación, pero ahora Sustituyendo las constantes
también vamos a sustituir la constante obtenidas:
𝐶
ECUACIÓN PARTICULAR
Las segundas condiciones son:

𝑡 = 80 , 𝑇(𝑡) = 43.2°𝐶 𝑦 𝑐𝑜𝑛 𝑇0 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 = 24.5°𝐶 𝑇(𝑡) = 62𝑒 (−0.0149826)∙𝑡 + 24.5

𝐶 = 62
Sustituyendo: Vamos a comprobar con un dato al
azar para verificar que nuestra
43.2 = 62𝑒𝑘(80) + 24.5 ecuación se aproxima a los datos
obtenidos en la tabla de nuestra
Ahora despejando a la constante de experimentación.
proporcionalidad 𝑘, tenemos que:

43.2 − 24.5 = 62𝑒𝑘(80)

Para un tiempo de 75 min, ¿cuánta
18.7 = 62𝑒𝑘(80) temperatura tendrá nuestro sistema?
 CAMPO DIRECCIONAL
(−0.0149826)∙(75)
𝑇(75) = 62𝑒 + 24.5
Ahora calcularemos el campo
𝑇(75) = 44.65 direccional de nuestra ecuación
diferencial.
Y como observamos en nuestra tabla
para un tiempo de 75 min obtuvimos Partiendo de nuestra ecuación
una lectura de 44.3°𝐶 que se diferencial.
aproxima muchísimo a la temperatura
que obtuvimos mediante la ecuación 𝑑𝑇
= 𝑘(𝑇 − 𝑇0 )
particular. 𝑑𝑡

𝟒𝟒. 𝟔𝟓°𝑪~𝟒𝟒. 𝟑°𝑪 Para determinar el campo direccional

de la ecuación anterior
Ahora que hemos comprobado debemos verificar que tenga la forma
calcularemos las demás temperaturas
usando la ecuación y así comprobar la 𝑑𝑇
= 𝑓(𝑡, 𝑇)
aproximación y compararlas con los 𝑑𝑡
datos obtenidos mediante el
experimento. Entonces se verifica que:

𝑑𝑇
Tabla de datos calculados mediante = 𝑘(𝑇 − 𝑇0 )
la ecuación particular 𝑑𝑡
Tiempo en Temperatura en °C
minutos donde claramente
0 86.5
𝑑𝑇
5 82.02 𝑇′ = , 𝑓(𝑡, 𝑇) = 𝑘(𝑇 − 𝑇0 )
10 77.87 𝑑𝑡
15 74.02
20 70.44 Para encontrar el campo direccional
25 67.13 𝑑𝑇
se debe sustituir 𝑑𝑡 por 𝑚, para
30 64.05
35 61.19 encontrar el valor de la pendiente en
40 58.55 un punto (𝑡; 𝑇) dado por la función
𝑓(𝑡, 𝑇).
45 56.09
50 53.81
Por lo que la ecuación a trabajar será
55 51.69
60 49.73 𝑚 = 𝑘(𝑇 − 𝑇0 )
65 47.91
70 46.22 Sustituyendo 𝑘 y 𝑇0
75 44.65
80 43.2 𝑚 = −0.0149826(𝑇 − 24.5)
GRÁFICAS

Ahora vamos a demostrar

gráficamente tanto la gráfica del
campo direccional de la ecuación
diferencial y las gráficas de las tablas
de datos tanto las que obtuvimos
midiendo experimentalmente y los
datos obtenidos mediante la ecuación
diferencial particular.

GRÁFICA DE CAMPO Como podemos observar tanto la

DIRECCIONAL gráfica azul correspondiente a los
datos obtenidos mediante la ecuación
𝑚 = −0.0149826(𝑇 − 24.5) y la gráfica verde correspondiente a
los datos obtenidos
experimentalmente, se asemejan
demasiado y se aproximan y así
podemos comparar el fenómeno físico
y observar si la ley es verídica.

III. CONCLUSIONES

Como podemos ver la ecuación

diferencial de la ley de enfriamiento de
Newton, nos ayuda a encontrar la
variación de la temperatura respecto
al tiempo, podemos notar que dicha
ley dice que mientras más se vaya
enfriando el objeto o sistema este
tenderá a tardar más tiempo para
llegar hasta la temperatura ambiente,
entre más se acerque al límite más
tiempo tardará en enfriarse, hemos
comprobado la ley tanto
experimentalmente y analíticamente,
Podemos apreciar el campo observando las gráficas podemos
direccional de la ecuación diferencial. comparar que tanto los resultados
Ahora a continuación tenemos las analíticos como las mediciones
gráficas de las tablas de datos de la experimentales se aproximan o se
temperatura, de las medidas tomadas asemejan demasiado, también
experimentalmente y las calculadas observamos el campo direccional de
mediante la ecuación para nuestra ecuación, concluyendo así la
compararlas. comprobación de dicha ley y de
manera satisfactoria.
 IV. AGRADECIMIENTOS

El autor agradece al IPN y en especial

al profesor Josué Roberto Hernández,
por el asesoramiento y aclarar mis
dudas respecto a la asignatura de
ecuaciones diferenciales.

V. REFERENCIAS

ENCICLOPEDIA ASIGNA (2021)

Definición de La Ley de Enfriamiento
de Newton. enciclopedia.net.
Consultado el 07 de noviembre de
2022. https://enciclopedia.net/ley-de-
enfriamiento-de-newton/

UNAM (s/f) Definición de las

ecuaciones diferenciales.
prometeo.matem.unam.mx.
Consultado el 07 de noviembre de
2022.
http://prometeo.matem.unam.mx/recu
rsos/Licenciatura/UNADM_Adapt/scor
m_player/941/content/escenas/2_Inici
o_1.html