¿Cómo puedo determinar el coeficiente kkk en dTdt=−k(T−100∘C)dTdt=−k(T−100∘C) \dfrac{dT}{dt} = -k(T - 100 \mathrm{^\ círculo C}) ?

Hace poco dediqué un tiempo a cocinar y tengo curiosidad por la evolución temporal de la temperatura del agua. Hice un experimento y la temperatura es de la forma

T = 100 C + ( T ( t 0 ) 100 C ) mi k ( t t 0 )
para alguna constante positiva k de alguna manera (por ejemplo, vea mi publicación anterior y mi respuesta ). En otras palabras,
d T d t = k ( T 100 C ) .

Mi pregunta: ¿Cómo puedo determinar el coeficiente k ? Por supuesto, puedo determinarlo mediante experimentos como se muestra a continuación. Pero lo que quiero hacer es entender esto como una función del volumen de agua. V , potencia del calentador PAG , y algunas otras variables relacionadas con un bote. Entonces, mi respuesta esperada es algo así como "Si el volumen (potencia) es norte veces mayor que el coeficiente es norte veces ( 1 / norte veces) más grandes".

Aquí están los datos experimentales y los gráficos de ajuste . ( V el agua está en la misma olla con un radio de 9 cm y calentada por IH corresponde al poder PAG . Mido su temperatura cada 30 segundos.)

Mi objetivo es determinar a priori la evolución temporal de la temperatura del agua. Entonces, si otro enfoque es mejor, por favor dígame. Agradecería si me ayudas. Gracias.

Estoy tratando de responder, pero me preocupa que pueda haber alguna iteración en esta pregunta. En primer lugar, ¿qué están haciendo los 100 grados allí? Está utilizando la ley de calentamiento/enfriamiento de Newton, y las dos temperaturas relevantes son la del agua y la del aire. Además, necesitará una nueva forma de esa ecuación que incluya la entrada de calor adicional. Una complicación es que realmente no puedes creer la cantidad de vatios, porque no todos llegan al bote. Eso da 2 variables independientes, pero deberíamos poder resolverlas 1 por 1.
Lo que dijo AlanSE y también: ¿Cómo determinaste β ? Tus números parecen incorrectos. más grande β corresponde a un aumento de temperatura más rápido, pero sus números van en sentido contrario. Además, solo por el aspecto de sus datos, no puede reclamar ni cerca de tres cifras significativas. :) Una vez que identifique los parámetros importantes, probablemente podría obtener una estimación razonable para β simplemente por análisis dimensional.
@AlanSE Pensé que el gráfico es lineal. Pero si es así, la temperatura supera los 100 C . Así que probé alguna otra función creciente acotada arriba por 100. Así es como se me ocurre. Luego, lo reescribo de forma de la ley de enfriamiento de Newton. Para la segunda parte, también lo pienso. Pero no estoy seguro de qué valores de calentamiento son primarios: uno escapó de un lado a través de una olla y el otro escapó de la superficie del agua. Así que traté de resolver el más simple pero me quedé atascado.
@MichaelBrown Determino el β con temperatura de 0 segundos y 180 segundos. Pero es probable que mis experimentos estén equivocados. Solo lo hago con objetos comunes de cocina. No es exacto. Ya probé el análisis dimensional pero los parámetros importantes son desconocidos para mí (yo me especializo en matemáticas y no estoy familiarizado con este tipo de problemas).
Revisé sus datos y ajustes, y en mi opinión, su fórmula no se ajusta bien (sin barras de error en algunos lugares). Dado que el sistema está experimentando una transición obvia de un estado a otro (el agua comienza a hervir), tal vez podría lograr mejores ajustes si intentara dividir los datos en dos dominios (antes de hervir y dentro de la ebullición) y ajustarlos de forma independiente. , o con alguna condición coincidente en el borde de los dominios.

Respuestas (3)

La ley relevante es la conservación de la energía.

cambio en la energía térmica del agua = energía puesta energía perdida en el aire

Aplique la ley de enfriamiento de Newton para escribir la energía perdida en el aire. La entrada de energía del ojo es constante. Entonces el cambio de energía del agua es una cuestión de notación que estoy introduciendo.

d q d t = q ˙ = W ˙ A h ( T ( t ) T a i r )

También tenga en cuenta que a partir de la definición de calor específico, la energía en la olla es q = C pag metro T , que ignora cualquier compensación, lo cual no es un problema porque todas las temperaturas pueden ser relativas.

d q d t = d T d t metro C pag = W ˙ + A h T a i r A h T ( t )

d T d t = W ˙ metro C pag + A h metro C pag T a i r A h metro C pag T ( t )

Esta es la forma que yo quería, para poder sustituir en tu β .

d T d t = W ˙ metro C pag + β T a i r β T ( t )

La sustitución que acabo de hacer define β . En cuanto a la unidad, puede describir esto más a fondo. Pero en resumen, h es el coeficiente de transferencia de calor específico del área.

Para crear experimentos útiles, veamos la solución de estado estacionario de esta ecuación.

β T ( ) = W ˙ metro C pag + β T a i r

Supongo que tienes una medida directa de la temperatura del aire y de la olla. Eso nos deja con dos incógnitas. no sabemos β , y no estoy de acuerdo en que sepas W ˙ . Me veo obligado a tratarlo como un desconocido. En realidad, introduciría otra variable α = W ˙ metro C pag .

β T ( ) = α + β T a i r

T ( ) = α β + T a i r

Aunque no lo sepamos α , podría hacer uso de la linealidad que viene con él. En otras palabras, si usa 200 W y luego 400 W, estoy bastante seguro de que la entrada de calor aumenta en un factor de 2. Para β tenemos que ser cuidadosos. No puede estar (muy) cerca de hervir o afectará las corrientes convectivas. Y el nivel del agua también afectará mucho su valor (principalmente debido a metro , pero también A ).

Propuesta de experimento

paso 1

Para cualquier volumen de agua arbitrario, obtenga valores de estado estable para diferentes niveles de potencia y utilícelos para determinar α / β , que es específico para ese nivel de potencia y volumen. Matemáticamente, eso da:

α β = T ( ) T a i r

Es posible que desee hacer esto para varios niveles de potencia. Simplemente aumente un poco el nivel de potencia, espere un poco y luego registre la temperatura. Eso debería ser lineal, pero el grado en que es lineal es una buena prueba de qué tan precisa es la imagen mental.

paso 2

Hemos estado cambiando el nivel de potencia, así que ahora detengámonos, mantengamos el mismo volumen de agua y el nivel de potencia, y observemos el cambio de temperatura a lo largo del tiempo. La ecuación diferencial que escribí tiene una solución que básicamente resulta en lo siguiente, si la olla comienza a temperatura ambiente y luego se calienta.

T ( t ) = T a i r + ( 1 mi β t ) ( T ( ) T a i r )

Debe saber del paso anterior cuál es la temperatura final, por lo que es conocida.

β = en ( T ( t ) T a i r T ( ) T a i r + 1 ) / t

β = en ( T ( ) T a i r T ( ) T ( t ) ) / t

Siempre que haya comenzado correctamente con la temperatura del aire y haya encontrado correctamente la temperatura final, solo necesita un punto de datos en el medio para resolver la ecuación anterior.

Esto brinda una forma de encontrar la constante de tiempo de calentamiento, así como formas de interpretar los coeficientes. Recuerda que la constante de tiempo es específica a la cantidad de agua en la olla. Si quisiera obtener una comprensión más completa, usaría los datos para resolver α , y luego correlacione eso con lo que la estufa dice que es el calor, para obtener un factor de corrección de la cantidad de calor perdido directamente en el aire desde el ojo.

Como sé poco de física, es difícil seguir tu explicación. Me tomaré tiempo para entenderlo. Déjame preguntar sobre la notación. En la primera ecuación, W significa trabajo que hace el calentador y por lo tanto W ˙ = 700 W , 1000 W en este caso no? Y qué A y h ¿representa?
@Taro Sí, lo entiendes. Los puntos indican que son por unidad de tiempo. Q y W y medidas de calor, y ponerles un punto los convierte en una medida de flujo de calor. A es el área y h es el coeficiente de transferencia de calor, que solo tienes que buscar en Google. Sin embargo, ninguno de estos es necesario para su problema. Matemáticamente, hice solo 2 predicciones que importan: la temperatura va linealmente con la entrada de calor y como exponencial decreciente con el tiempo. Eso significa que tienes dos parámetros para resolver. yo también debí haber escrito T ( ) = T a i r + γ W ˙ como alternativa al coeficiente de calefacción
me confundo un poco con la notación T ( ) . Creo que tiene sentido porque la temperatura está aumentando y acotada arriba; por lo tanto existe T = límite t T ( t ) . Así que para suficiente grande t , su derivada d T / d t es cero y conduce a α / β = T T aire . Sin embargo, en este caso, T = 100 C , ¿no es así? En ese caso, α / β ¿Solo depende solo de la temperatura ambiente? ¿Eso significa h T aire ? y obtengo T = T aire + α β + mi β t . (Probablemente esto podría ser solo un pequeño error tipográfico).
@Taro La física que hemos discutido hasta ahora es el balance de calor y la ley de enfriamiento de Newton. Si la temperatura final será de 100 C, supongo que quiere decir que hervirá. Eso es extremadamente problemático. Todas las suposiciones que usamos ya no son válidas. El comportamiento es por partes (en el sentido más simple). Si hierves, no puedes encontrar T del mismo modo. Todavía podría aplicar estas ecuaciones para el calentamiento hasta la ebullición, tal vez ajustando un ajuste exponencial a los datos en Excel. Podrías hacer eso con los datos que publicaste. Pero con el inicio de la ebullición, la física cambia.
Me temo que no entiendo bien. Si las ecuaciones se vuelven inválidas en el punto de ebullición, ¿qué quiere decir con "solución de estado estacionario"? ¿Podría explicarlo, por favor?
@Taro Inicialmente no sabía que ibas a hervir. Eso cambia el problema y lo hace más difícil. La solución de estado estacionario no funciona cuando se trata de ebullición porque conocemos la temperatura final. Un modelo simplista es que la temperatura aumenta de acuerdo con la exponencial hasta el punto de ebullición y luego permanece constante a medida que el calor del ojo se transfiere a la energía del vapor. Dado que va a estar hirviendo, no tengo los objetivos lo suficientemente claros como para sugerir un curso de acción conciso.

Tienes tres procesos en marcha:

  1. el calor inyectado calienta el agua a un ritmo constante

  2. la sartén se enfría de acuerdo con la ley de Newton, es decir, la tasa de pérdida de calor es proporcional a la diferencia de temperatura entre la sartén y el ambiente

  3. pierdes calor por la evaporación del agua

Es fácil escribir una ecuación diferencial para los mecanismos 1 y 2, pero creo que el mecanismo 3 será problemático. Hay algo de información al respecto en el sitio web de Engineering Toolbox , pero esto requiere parámetros como el contenido efectivo de humedad del aire alrededor de la bandeja. Sospecho que necesitará medir esto en lugar de intentar calcularlo a partir de principios básicos.

Tenga en cuenta también que cuando la temperatura alcanza los 100ºC, el mecanismo cambia ya que la evaporación se vuelve efectivamente infinita, es decir, incluso con una gran entrada de energía, no puede elevar la temperatura del agua (mucho) por encima de los 100ºC. Esto significa que obtendrá una discontinuidad en la curva de calentamiento a 100ºC.

Puede hacerse una idea de la contribución del enfriamiento por evaporación (por debajo de 100 ºC) calentando la sartén con y sin una tapa ajustada.

Estoy mirando su enlace, tratando de encontrar una ocurrencia de la temperatura del agua en la expresión de la tasa de evaporación, y es difícil. Estoy pensando que está horneado en la presión parcial del vapor de agua en condiciones de saturación. Supongo que probablemente no sea lineal...
@AlanSE: sí, así es como interpreté la información en esa página. Hay ajustes empíricos para la presión de vapor frente a la temperatura, pero todo se va a complicar un poco.
Recientemente resolví un problema de mi pregunta anterior y dije que estaba pensando en este problema nuevamente. Allí, un término T 100 C aparece Su sugerencia no parece contener este término. (Y no puedo derivar este término de ningún primer principio en este momento). ¿Qué piensas al respecto? ¿Es casualidad que el término encaje?
La temperatura no puede subir por encima de los 100 °C porque el agua comienza a hervir, pero creo que hay un valor distinto de cero de dT/dt justo antes de llegar a los 100 °C y hay una discontinuidad a los 100 °C. En la práctica, esto sería muy difícil de ver porque el flujo de calor a través del agua toma un tiempo finito, por lo que no toda el agua estará a la misma temperatura. Mi conjetura es que esto elimina la discontinuidad, por lo que una ecuación que se aproxima a 100C asintóticamente da un buen ajuste aunque no sea estrictamente cierto.

mira que usamos integración simple ..

d T / d t = k ( T 100 )
esto está de acuerdo con la Ley de Enfriamiento de Newton, que dice que la tasa negativa de cambio de temperatura es directamente proporcional a la diferencia de la temperatura inicial de la sustancia y la temperatura circundante del medio ambiente.

d T / d t = k ( T 100 )
d T / ( T 100 ) = k d t
d T / ( T 100 ) = k d t
en ( T 100 ) = k t + C (1)

Poniendo t = 0 en el momento de origen y el correspondiente valor de T , obtenemos el valor de C .

Poner el valor de C en (1) podemos encontrar k en términos de T y t .

Podemos poner los respectivos valores para la respuesta final.

Gracias por su respuesta. Pero eso es básicamente lo que hice en mi publicación anterior como escribí en la pregunta... El problema es cómo determinar el k con respecto a los alrededores algo como V o PAG . Y de donde sale el DE si realmente gobierna el sistema.