Как определить коэффициент kkk в dTdt=−k(T−100∘C)dTdt=−k(T−100∘C) \dfrac{dT}{dt} = -k(T - 100 \mathrm{^\ круг С}) ?

Недавно я трачу некоторое время на приготовление пищи, и мне любопытно, как изменяется во времени температура воды. Я провел некоторый эксперимент, и температура имеет вид

Т "=" 100 С + ( Т ( т 0 ) 100 С ) е к ( т т 0 )
для некоторой положительной постоянной к как-нибудь (например, см. мой предыдущий пост и мой ответ ). Другими словами,
г Т г т "=" к ( Т 100 С ) .

Мой вопрос: Как я могу определить коэффициент к ? Конечно, я могу определить это опытным путем, как показано ниже. Но то, что я хочу сделать, это понять это как функцию объема воды. В , мощность нагревателя п и некоторые другие переменные, связанные с горшком. Итак, мой ожидаемый ответ примерно такой: «Если громкость (мощность) н раз больше, чем коэффициент н раз ( 1 / н раз) больше».

Вот экспериментальные данные и подгоночные графики . ( В вода находится в той же кастрюле радиусом 9 см и нагрета IH соответствует мощности п . Я измеряю его температуру каждые 30 секунд.)

Моя цель - определить временную эволюцию температуры воды априори. Так что, если другой подход лучше, пожалуйста, скажите мне. Буду признателен, если вы мне поможете. Спасибо.

Я пытаюсь ответить, но боюсь, что на этот вопрос может быть какая-то итерация. Во-первых, что там делают 100 градусов? Вы используете закон Ньютона о нагревании/охлаждении, и две соответствующие температуры — это температура воды и температура воздуха. Кроме того, вам понадобится новая форма этого уравнения, которая включает дополнительное тепловложение. Сложность заключается в том, что вы не можете на самом деле поверить в количество Ватт, потому что не все они попадают в банк. Это дает 2 независимые переменные, но мы должны иметь возможность решать их 1 на 1.
Что сказал AlanSE, а также: Как вы определили β ? Ваши цифры кажутся неправильными. Больше β соответствует более быстрому повышению температуры, но ваши цифры идут в другую сторону. Также, судя по вашим данным, вы не можете претендовать на три значащие цифры. :) Как только вы определите важные параметры, вы, вероятно, сможете получить разумную оценку для β просто размерный анализ.
@AlanSE Я думал, что график линейный. Но если так, то температура становится выше 100 С . Поэтому я попробовал другую возрастающую функцию, ограниченную сверху 100. Вот как я ее придумал. Вскоре после этого я переписал его в форме закона охлаждения Ньютона. Что касается второй части, я тоже об этом думаю. Но я не уверен, какая теплотворная способность первична: одна ушла сбоку через котел, а другая ушла с поверхности воды. Итак, я попытался решить простейшую, но застрял.
@MichaelBrown Я определяю β с температурой 0 секунд и 180 секунд. Но вполне вероятно, что мои эксперименты ошибочны. Я просто делаю это с помощью обычных кухонных предметов. Это не точно. Я уже пробовал размерный анализ, но важные параметры мне неизвестны (я специализируюсь на математике и не знаком с такими задачами).
Я посмотрел на ваши данные и подгонки, и, на мой взгляд, ваша формула довольно плохо подходит (в некоторых местах нет погрешностей). Поскольку система совершает очевидный переход из одного состояния в другое (вода начинает кипеть), возможно, вы могли бы достичь лучшего соответствия, если бы попытались разделить данные на две области (до кипения и в кипящем состоянии) и сопоставить их независимо друг от друга. , либо с некоторым условием совпадения на границе доменов.

Ответы (3)

Соответствующий закон - сохранение энергии.

изменение тепловой энергии воды "=" вложенная энергия энергия теряется в воздухе

Примените закон охлаждения Ньютона, чтобы записать энергию, теряемую воздухом. Поступление энергии от глаза постоянно. Тогда изменение энергии воды — это вопрос обозначений, которые я ввожу.

г Вопрос г т "=" Вопрос ˙ "=" Вт ˙ А час ( Т ( т ) Т а я р )

Также обратите внимание, что из определения удельной теплоемкости энергия в горшке равна Вопрос "=" С п м Т , что игнорирует любое смещение, что не является проблемой, поскольку все температуры могут быть относительными.

г Вопрос г т "=" г Т г т м С п "=" Вт ˙ + А час Т а я р А час Т ( т )

г Т г т "=" Вт ˙ м С п + А час м С п Т а я р А час м С п Т ( т )

Это форма, которую я хотел, чтобы я мог заменить в вашем β .

г Т г т "=" Вт ˙ м С п + β Т а я р β Т ( т )

Замена, которую я только что сделал, определяет β . По единицам вы можете описать это более подробно. Но вкратце, час - удельный коэффициент теплопередачи по площади.

Чтобы провести полезные эксперименты, давайте посмотрим на стационарное решение этого уравнения.

β Т ( ) "=" Вт ˙ м С п + β Т а я р

Я предполагаю, что у вас есть прямое измерение температуры воздуха, а также горшка. Это оставляет нас с двумя неизвестными. мы не знаем β , и я не согласен, что вы знаете Вт ˙ . Я вынужден относиться к этому как к неизвестному. На самом деле, я бы ввел еще одну переменную α "=" Вт ˙ м С п .

β Т ( ) "=" α + β Т а я р

Т ( ) "=" α β + Т а я р

Хотя мы не знаем α , я мог бы использовать связанную с этим линейность. Другими словами, если вы используете 200 Вт, а затем 400 Вт, то я уверен, что тепловложение увеличится в 2 раза. β мы должны быть осторожны. Вы не можете быть (очень) близко к кипению, иначе это повлияет на конвективные потоки. И уровень воды тоже сильно повлияет на его стоимость (в основном из-за м , но и А ).

Предложение эксперимента

шаг 1

Для любого произвольного объема воды получите установившиеся значения для различных уровней мощности и используйте их для определения α / β , который зависит от этого уровня мощности и громкости. Математически это дает:

α β "=" Т ( ) Т а я р

Возможно, вы захотите сделать это для нескольких уровней мощности. Просто немного увеличьте уровень мощности, немного подождите, затем запишите температуру. Это должно быть линейно, но степень линейности является хорошим показателем того, насколько точна ментальная картина.

шаг 2

Мы меняли уровень мощности, так что теперь давайте остановимся, сохраним объем воды и уровень мощности одинаковыми и посмотрим на изменение температуры с течением времени. Дифференциальное уравнение, которое я написал, имеет решение, которое в основном сводится к следующему, если кастрюля начинается с температуры окружающей среды, а затем нагревается.

Т ( т ) "=" Т а я р + ( 1 е β т ) ( Т ( ) Т а я р )

Вы должны знать из предыдущего шага, какова конечная температура, так что это известно.

β "=" п ( Т ( т ) Т а я р Т ( ) Т а я р + 1 ) / т

β "=" п ( Т ( ) Т а я р Т ( ) Т ( т ) ) / т

При условии, что вы правильно начали с температуры воздуха и правильно нашли конечную температуру, вам нужна только одна точка данных в середине, чтобы решить приведенное выше уравнение.

Это дает возможность найти постоянную времени нагрева, а также способы интерпретации коэффициентов. Помните, что постоянная времени зависит от количества воды в горшке. Если бы вы хотели получить более полное представление, я бы использовал данные для решения α , а затем сопоставьте это с тем, что печь говорит о тепле, чтобы получить поправочный коэффициент количества тепла, теряемого непосредственно в воздух из глаз.

Поскольку я мало разбираюсь в физике, мне трудно понять ваше объяснение. Я возьму время, чтобы понять это. Позвольте мне спросить об обозначениях. В первом уравнении Вт означает работу, которую выполняет нагреватель, и, следовательно, Вт ˙ "=" 700 Вт , 1000 Вт в этом случае, верно? И что А и час означает?
@ Таро Да, ты понял. Точками указано, что они в единицу времени. Q и W и меры тепла, а расставление точек делает их мерой теплового потока. A — это площадь, а h — коэффициент теплопередачи, который вам просто нужно найти в Google. Однако ни один из них не является необходимым для вашей проблемы. Математически я сделал только 2 предсказания, которые имеют значение: температура изменяется линейно с подводимой теплотой и как экспоненциально убывающая со временем. Это означает, что у вас есть два параметра для решения. Я также должен был написать Т ( ) "=" Т а я р + γ Вт ˙ как альтернатива коэффициенту нагрева
немного путаюсь в обозначениях Т ( ) . Я думаю, что это имеет смысл, потому что температура увеличивается и ограничена сверху; следовательно, существует Т "=" лим т Т ( т ) . Так что для достаточно больших т , его производная г Т / г т равен нулю, и это приводит к α / β "=" Т Т воздух . Однако в этом случае Т "=" 100 С , не так ли? Если так, α / β зависит только от комнатной температуры? Это значит час Т воздух ? И я получаю Т "=" Т воздух + α β + е β т . (Возможно, это просто небольшая опечатка.)
@Taro Физика, которую мы обсуждали до сих пор, - это тепловой баланс и закон охлаждения Ньютона. Если конечная температура будет 100 С, то, я так понимаю, вы имеете в виду, что она будет кипеть. Это крайне проблематично. Все предположения, которые мы использовали, больше не действительны. Поведение кусочное (в самом простом смысле). Если кипятить, то не найти Т таким же образом. Вы все еще можете применить эти уравнения для нагрева до кипения, возможно, подгоняя экспоненциальную подгонку к данным в Excel. Вы можете сделать это с данными, которые вы разместили. Но с началом кипения физика меняется.
Боюсь, я плохо понимаю. Если уравнения становятся недействительными при температуре кипения, что вы подразумеваете под «стационарным решением»? Не могли бы вы уточнить это, пожалуйста?
@ Таро, я изначально не знал, что ты собираешься закипеть. Это меняет проблему и усложняет ее. Стационарное решение не работает, когда речь идет о кипении, потому что мы знаем конечную температуру. Упрощенная модель состоит в том, что температура увеличивается по экспоненте до точки кипения, а затем остается постоянной по мере того, как тепло глаза передается энергии пара. Учитывая, что вы собираетесь кипеть, я недостаточно ясно понимаю ваши цели, чтобы предложить краткий план действий.

У вас происходит три процесса:

  1. подаваемое тепло нагревает воду с постоянной скоростью

  2. сковорода охлаждается в соответствии с законом Ньютона, т.е. скорость потери тепла пропорциональна разнице температур между сковородой и окружающей средой

  3. вы теряете тепло из-за испарения воды

Легко написать дифференциальное уравнение для механизмов 1 и 2, но я думаю, что механизм 3 будет проблематичным. На веб-сайте Engineering Toolbox есть некоторая информация об этом, но для этого требуются такие параметры, как эффективное содержание влаги в воздухе вокруг поддона. Я подозреваю, что вам нужно измерить это, а не пытаться вычислить из первых принципов.

Обратите также внимание, что когда температура достигает 100ºC, механизм меняется, так как испарение становится фактически бесконечным, т.е. даже при большом подводе энергии вы не можете поднять температуру воды (намного) выше 100ºC. Это означает, что вы получите разрыв кривой нагрева при 100ºC.

Вы можете получить некоторое представление о вкладе испарительного охлаждения (ниже 100ºC), нагревая кастрюлю с плотно закрывающейся крышкой и без нее.

Я смотрю по вашей ссылке, пытаюсь найти вхождение температуры воды в выражение для скорости испарения, и это сложно. Я думаю, что это связано с парциальным давлением водяного пара в условиях насыщения. Я думаю, что это, вероятно, нелинейно...
@AlanSE: да, именно так я интерпретировал информацию на этой странице. Есть эмпирические подгонки зависимости давления пара от температуры, но все это становится немного запутанным.
Недавно я решил проблему моего предыдущего вопроса и заявил, что снова думаю об этой проблеме. Там термин Т 100 С появляется. Ваше предложение, похоже, не содержит этого термина. (И я не могу сейчас вывести этот термин из каких-либо первых принципов.) Что вы думаете об этом? Случайно ли подходит этот термин?
Температура не может подняться выше 100°С, потому что вода начинает кипеть, но я думаю, что существует ненулевое значение dT/dt как раз перед тем, как вы достигнете 100°С, и есть разрыв при 100°С. На практике это было бы очень трудно увидеть, потому что поток тепла через воду занимает конечное время, поэтому вся вода не будет иметь одинаковую температуру. Я предполагаю, что это размывает разрыв, поэтому уравнение, которое асимптотически приближается к 100 ° C, дает хорошее соответствие, даже если оно не совсем верно.

см. мы используем простую интеграцию ..

г Т / г т "=" к ( Т 100 )
это соответствует закону охлаждения Ньютона, который гласит, что отрицательная скорость изменения температуры прямо пропорциональна разнице начальной температуры вещества и окружающей температуры окружающей среды.

г Т / г т "=" к ( Т 100 )
г Т / ( Т 100 ) "=" к г т
г Т / ( Т 100 ) "=" к г т
п ( Т 100 ) "=" к т + С (1)

положить т "=" 0 на момент возникновения и соответствующее значение Т , получаем значение С .

Помещение значения С в (1) мы можем найти к с точки зрения Т и т .

Мы можем поставить соответствующие значения для окончательного ответа.

Спасибо за ваш ответ. Но это в основном то, что я сделал в своем предыдущем посте, как я написал в вопросе... Проблема в том, как определить к по отношению к окружению что-то вроде В или п . И откуда берется DE, если он действительно управляет системой.