Logo
Logo

Динамика волны кристаллизации

Оборудование:

  1. Ноутбук
  2. Электрическая схема с датчиком температуры
  3. Термометр
  4. Весы
  5. Секундомер
  6. Капельницы
  7. Линейка
  8. Маркер
  9. Ножницы
  10. Дистиллированная вода
  11. Ведёрки (кювета для измерений, ведро для слива, ведро с холодной водой)
  12. Салфетки для поддержания чистоты на рабочем месте
  13. Шприцы 10 мл
  14. Пластиковые пробирки 50 мл
  15. Пакет с ацетатом натрия
  16. Чайник

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

 

Переохлажденная жидкость — метастабильное (неустойчивое) состояние вещества, при котором вещество остаётся в жидкой фазе даже при температуре ниже точки кристаллизации (замерзания). Оно достигается либо при быстром охлаждении, либо в отсутствии центров кристаллизации (этот случай и рассматривается в рамках данной задачи). Как было сказано, термодинамическое равновесие такой жидкости неустойчиво, поэтому даже при слабом механическом воздействии или контакте с твёрдыми частицами вещества, метастабильное состояние рушится и жидкость кристаллизуется, переходя в устойчивое равновесие

На фото переохлажденную воду из бутылки наливают в бокал. Из-за соприкосновения с уже кристаллизовавшейся жидкостью, струя воды "на лету" переходит в кристаллическое состояние.

Раствор стоит изготавливать в пробирке с крышкой. При этом рекомендуется готовить объёмом около $25~мл$. Мы готовим перенасыщенный при комнатной температуре раствор, поэтому порошок не будет растворяться в воде полностью. Чтобы он растворился, смесь необходимо нагреть методом водяной бани — положив пробирку в чайник, заполненный водой.

Часть A. Разминка (1.2 балла)

Выше приведена фазовая диаграмма раствора ацетата натрия в координатах $T$ температуры от массовой концентрации $c=m_{ац}/(m_{ац}+m_{вод})$ ацетата в растворе. В зоне Liquid равновесное состояние – водный раствор, а в зоне Liquid+SAT равновесное состояние – насыщенный раствор с осадком.

Явление переохлаждения заключается в том, что раствор из состояния Liquid охлаждают ниже температуры перехода Liquid $\to$ Liquid+SAT и при этом не образуется осадок, т.е. система находится не в равновесии.

Температуру этого перехода мы будем называть температурой кристаллизации. Фазу Liquid будем называть жидкой, а фазу Liquid+SAT твердой.

A1  0.20 Опираясь на приведённый график, определите минимальную концентрацию $c_\min$, при которой возможно наблюдать переохлаждённое состояние при комнатной температуре.

На деле диапазон концентраций, в котором переохлажденное состояние оказывается достаточно устойчивыми: $c \in [0.40, 0.55]$. Именно в границах этого диапазона мы будем работать в ходе данной задачи.

Измерим температуру кристаллизации и плотность раствора с концентрацией $c_0=0.50$. 

Для измерения температуры фазового перехода нужно привести раствор в переохлаждённое состояние. Для этого необходимо:

  1. Приготовить смесь дистиллированной воды и порошка ацетата натрия в пробирке с крышкой. Запишите значения масс составляющих. Имейте ввиду, ацетат не полностью растворится в воде при комнатной температуре
  2. Растопить смесь в пробирке в горячей водяной бане (чайнике с кипятком). Имейте ввиду, что могут образовываться комки твёрдого вещества, которые очень долго растворяются сами, поэтому их стоит разбить стеклянной палочкой, либо встряхнув пробирку.
  3. Достать пробирку из бани. Положите в пробирку термометр.
  4. Охладить раствор до температуры $25^\circ\text{C}-35^\circ\text{C}$, поместив пробирку в ведёрко с холодной водой
  5. Запустить кристаллизацию, насыпав сверху в открытую пробирку немного порошка ацетата натрия

A2  0.30 Как можно точнее определите плотность раствора $\rho_l$, когда он находится в жидком состоянии.

Ответ запишите в единицах СИ.

A3  0.70 В ходе кристаллизации из переохлажденного состояния выделяется тепло и температура раствора растет до температуры $T_0$, которая немного меньше чем температура кристаллизации.

Определите величину $T_0$. Ответ запишите в $^\circ\text{C}$.

Часть B. Зависимость от температуры (7.5 балла)

Как вы могли заметить при выполнении эксперимента из части А, что кристаллизация идет не мгновенно по всему объему, а распространяется фронтом. В этой части мы измеряем скорость $v$ распространения фронта волны кристаллизации в зависимости от температуры, в которой находилась переохлажденная жидкость до начала кристаллизации.

 

Для проведения эксперимента по измерению скорости волны кристаллизации, необходимо следовать данной инструкции:

  1.  Засыпать в шприц порошок ацетата натрия (его мы будем использовать как стартер волны кристаллизации)
  2. Разрезать, затем сделать на ней пометки каждые несколько сантиметров
  3. Нагреть в кипятке один из концов отрезка трубки (чтобы можно было вставить шприц в трубку)
  4. Набрать в трубку раствор, оставив расстояние до шприца с порошком (чтобы случайно не запустить кристаллизацию)
  5. Закрыть свободный конец трубки роликовым зажимом
  6. Положить трубку с раствором в ведёрко с водой температуры $T_\text{out}$ и оставить на некоторое время, периодически слегка помешивая воду в ведёрке (чтобы раствор принял температуру среды)
  7. Запустить кристаллизацию, всосав раствор в трубку до шприца с порошком и включить секундомер

Известно, что при температурах близких к температуре плавления (обычно в районе $15-20^\circ C$ ниже $T_0$), характер фронта волны кристаллизации заметно отличается от меньших температур, из-за чего отличается вид зависимости $v(T_\text{out})$ на низких и на высоких температурах. Пусть характерная температура, отвечающая этому изменению — $T_\text{c}$.

B1  6.00 Пользуясь описанной выше методикой измерений для раствора с $c = c_0$ снимите зависимость $v(T_\text{out})$ в интервале температур $T_\text{out} \in [0^\circ C, T_0]$. В диапазонах $T_\text{out} < T_\text{c}$ и $T_\text{out} > T_\text{c}$ должно быть как минимум по 4 экспериментальные точки.

Загрузите обработанные данные в виде Excel таблицы (после автоматического сглаживания в выбранном вами диапазоне).

B2  0.50 Постройте график полученной зависимости. Сделайте скриншот и загрузите в систему.

B3  0.70 Определите значение $T_\text{c}$, соответствующее вашему эксперименту. Ответ запишите в $^\circ\text{C}$.

B4  0.30 Определите значение коэффициента наклона графика скорости $v$ от температуры $T_\text{out}$ в диапазоне $T_\text{out} > T_\text{c}$. Ответ запишите в единицах СИ.

Часть C. Определение теплоёмкостей (4.5 балла)

На рабочем столе вашего компьютера находятся две программы.

Первая из них: « recorder.exe » – программа для измерений с помощью автоматического термометра. Вторая из них: « analyzer.exe » — программа для начальной обработки полученных данных, полученный с помощью программы « recorder.exe ».

На выходе из программы « recored.exe » вы получаете Excel таблицу с предобработанными данными, с которыми вам предстоит работать.

Тепловая мощность $P$, переходящая от через стенки пробирки от окружающей воды к содержимому пробирки, обусловлена законом Ньютона-Рихмана:
\[
P = K \Delta T \tag{1},
\]
где $K$ — коэффициент теплопередачи.

Измерим коэффициент теплопередачи пробирки $K_\text{tube}$. Для этого снимем процесс остывания пробирки с водой, поставленной в ведерко с водой комнатной температуры. Уровни жидкости в ведерки и пробирке должны практически совпадать!

Объём воды в пробирке должен быть равным $25~мл$.

C1  0.50 Снимите зависимость $T$ температуры воды в пробирке от времени $t$  в течение 10 минут.

Загрузите обработанные данные в виде Excel таблицы (после автоматического сглаживания в выбранном вами диапазоне).

С помощью уравнения теплового баланса для воды внутри пробирки, можно определить её коэффициент теплопередачи $K_\text{tube}$.

Известно, что $c_\text{w} = 4200~\dfrac{кДж}{кг\cdot^\circ\text{C}}$ — удельная теплоёмкость воды, $\rho_\text{w} = 1000~кг/м^3$ — плотность воды.

C2  0.40 Используя уравнение (1) и уравнение теплового баланса, выразите скорость изменения температуры $\dot T = dT/dt$ воды внутри пробирки через температуру воды в ведерке $T_\text{out}$, температуру $T$ и объем $V$ воды в пробирке и известные константы.

C3  0.30 Постройте график зависимости, измеренной в пункте C1 в тех координатах, в которых он будет линейным. Загрузите скриншот графика.

C4  0.50 Определите значение коэффициента теплопередачи пробирки $K_\text{tube}$. Ответ запишите в единицах СИ.

Теперь повторим эксперимент для ацетата ацетата натрия. Приготовьте $25~мл$ раствора с концентрацией $c = 0.45$ (можете приготовите больше и слить избыток).

C5  0.70 Снимите зависимость температуры $T$ раствора ацетата натрия от времени $t$ при остывании в ведерке с водой комнатной температурой в течение 10 минут.

Загрузите обработанные данные в виде Excel таблицы (после автоматического сглаживания в выбранном вами диапазоне).

C6  0.70 Снимите зависимость температуры $T$ кристаллизовавшегося ацетата натрия от времени $t$ при остывании в ведерке с водой комнатной температуры в течение 15 минут.

Загрузите обработанные данные в виде Excel таблицы (после автоматического сглаживания в выбранном вами диапазоне).

Далее используются следующие обозначения: $\rho_l, c_l$ — плотность и удельная теплоёмкость смеси в жидкой фазе, $\rho_s, c_s$ — в кристаллической.

C7  0.30 Постройте график зависимости измеренной в пункте C5 в тех координатах, в которых он будет линейным. Загрузите скриншот графика.

C8  0.30 Постройте график зависимости измеренной в пункте C6 в тех координатах, в которых он будет линейным. Загрузите скриншот графика.

C9  0.40 Определите значения произведения $c_l \rho_l$. Ответ запишите в единицах СИ.

C10  0.40 Определите значение произведения $c_s \rho_s$. Ответ запишите в единицах СИ.

Часть D. Определение теплопроводности раствора (6.8 балла)

В этой части мы различаем показания датчика температуру $T_\text{m}$ и реальную температуру среды вокруг него $T$.

D1  0.40 Нагрейте воду в чайнике до кипения. Снимите зависимость температуры датчика $T_\text{m}$ от времени $t$, начав измерения до погружения в горячую воду $\sim 50^\circ C$, а закончив через $\sim 10 ~с$ после погружения.

Загрузите обработанные данные в виде Excel таблицы после автоматического сглаживания в диапазоне $\sim 5~с$ от момента, когда датчик попал в воду.

Как можно увидеть из графика измеренной вами зависимости, датчик температуры сам по себе имеет инерцию, то есть не моментально принимает температуру вещества, в которое он погружен. Это происходит из-за ненулевой теплоёмкости $C_\text{m}$ и конечного коэффициента теплопередачи $K_\text{m}$ его головки. Таким образом, при быстро протекающих процессах, которые мы будем рассматривать далее, его показания $T_\text{m}$ могут заметно отличаться от измеряемой им температуры $T$ вещества вблизи него.

D2  0.40 Запишите уравнение теплового баланса для датчика. Используя уравнение (1), выразите скорость изменения его показаний $\dot T_\text{m} = \dfrac{dT_\text{m}}{dt}$. Ответ выразите через $T, T_\text{m}$ и параметры датчика: $K_\text{m}$ и $C_\text{m}$.

D3  0.30 Постройте график зависимости, измеренной в пункте D1 в тех координатах, в которых он будет линейным. Загрузите скриншот графика.

 

D4  0.50 Определите характерное время инерции датчика $\tau = \dfrac{C_\text{m}}{K_\text{m}}$. Ответ запишите в единицах СИ.

 

Для следующих экспериментов для измерения температуры непосредственно внутри канала с волной вместо силиконовой трубки необходимо использовать саму пробирку с раствором, при этом датчик температуры расположите на середине высоты раствора в ней. Концентрацию раствора нужна по-прежнему $c = 0.45$, а температура окружающей воды $T_\text{out} \approx$ комнатной.

D5  0.70 Снимите зависимость $T_\text{m}(t)$ для кристаллизации раствора при температуре воды примерно равной комнатной.

Загрузите обработанные данные в виде Excel таблицы (после автоматического сглаживания в выбранном вами диапазоне).

D6  0.50 По измеренной в D5 зависимости пересчитайте температуру вещества $T(t)$ и постройте её график. На том же графике (в тех же осях) должна быть исходная кривая $T_\text{m}(t)$. Загрузите скриншот графика.

Опираясь на измеренную зависимость, можно вычислить $\lambda$ удельной теплоты кристаллизации раствора. Считайте, что вся тепловая энергия, выделяемая при кристаллизации переохлаждённой жидкости равна разности энергий, $Q_l - Q_s$ необходимых для нагрева жидкой фазы до температуры фазового равновесия $T_\text{melt}$ и охлаждения твёрдой от $T_\text{melt}$ до устанавливающейся после окончания процесса температуры.

D7  0.20 Из графика кривой фазового равновесия (из части A), определите температуру $T_\text{melt}$ отвечающую концентрации исследуемого раствора. Ответ запишите в $^\circ C$.

D8  0.50 Определите значение $\lambda$, соответствующее вашему эксперименту. Ответ запишите в единицах СИ.

Для описания теплового равновесия для смести внутри пробирки (канала) введём следующие величины:

  • $P$ — мощность потока тепловой энергии через рассматриваемое сечение канала
  • $P' = \dfrac{\partial P}{\partial x}$ — градиент мощности вдоль канала
  • $T' = \dfrac{\partial T}{\partial x}$ — градиент температуры вещества вдоль канала
  • $\dot{T} = \dfrac{d T}{d t}$ — скорость измерения температуры вещества во времени

Рассмотрим кусочек канала сечением $S$, длиной $\Delta x$ со смесью плотностью $\rho_l$, удельной теплоёмкостью $с_l$ и теплопроводностью $\varkappa_l$. Считайте, что если через левое сечение в этот объём жидкости поступает тепловая мощность $P$, то через правое — выходит $P + P' \Delta x$ (см. рисунок).

D9  0.40 Используя уравнение теплового баланса для выделенного объёма жидкости, получите выражение для $\dot T$. Ответ выразите через $P'$, $c_l$, $\rho_l$ и $S$.

Закон теплопроводности Фурье является следствием из закона Ньютона-Рихмана и связывает мощность теплового потока $P$ с градиентом температуры $T'$ и может быть записан как: \[P = \varkappa S T'.\]

D10  0.20 Пользуясь законом Фурье и результатом пункта D9, получите связь между $\dot T$ и $T'' = (T')'$.

Опираясь на стационарность движения волны кристаллизации, можно получить связь между производными, не привязанными ко времени \[T'' = A T', \tag{2}\] где $A = \text{const}$.

D11  1.00 Чему равно $A$? Ответ выразите через скорость волны $v$ и параметры системы.

Проинтегрировав выражение (2) получим, что температура $T$ связана со своим градиентом вдоль канала $T'$ как: \[T' = A T + \text{const}. \tag{3}\]А значит мы можем вернуться к измеряемым величинам, приведя уравнение (3) к виду: \[\dot T = B T + \text{const}. \tag{4}\]

D12  0.30 В уравнении (4) выразите значение $B$ через $v$ и параметры системы.

D13  0.40 Постройте график зависимости, измеренной в пункте D6 в тех координатах, в которых он будет линейным. Загрузите скриншот графика.

D14  0.50 Для определения $\varkappa_l$ необходима скорость фронта при новой исследуемой концентрации. Определите скорость распространения волны кристаллизации $v$ тем же методом, который использовался в части B. Ответ запишите в единицах СИ.

D15  0.50 Определите коэффициент теплопроводности раствора $\varkappa_l$. Ответ запишите в единицах СИ.