Идеальный эспрессо просчитали математически: физики вывели точное уравнение проницаемости кофейной таблетки

Каждое утро миллионы людей запускают сложный физический эксперимент. В небольшой металлической корзине — портафильтре (или, по-простому, рожке) — спрессовывается мелко смолотый кофе. Затем через эту кофейную таблетку просачивается вода, нагретая почти до температуры кипения, под давлением в 9 бар. Это в три-четыре раза превышает давление в автомобильной шине.

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

От скорости прохождения жидкости, напрямую зависит время контакта воды с твердым веществом. Именно это время определяет физические и химические процессы экстракции кофеина, масел и соединений, формирующих итоговый вкус напитка.

Для управления скоростью потока специалисты используют несколько переменных. Главные из них — это масса сухого кофе, степень уплотнения (темперовка) и размер помола. Изменение размера помола влияет на геометрию микроскопических каналов внутри кофейной таблетки. Более мелкий помол сужает каналы, увеличивая гидродинамическое сопротивление и замедляя поток жидкости. Однако до недавнего времени не существовало математической модели, которая могла бы точно связать физические характеристики молотого кофе с результирующей динамикой потока. Инженеры и геологи давно используют подобные модели для расчета проницаемости почв, горных пород или пористой керамики, но к процессам приготовления кофе такие методы системно не применялись.

Международная группа исследователей из университетов Германии, Великобритании, Франции и Канады решила восполнить этот пробел. Они поставили перед собой задачу сформулировать и экспериментально подтвердить уравнения, которые позволили бы предсказывать проницаемость кофейной таблетки на основе данных о размере частиц и плотности их упаковки.

Проблема лабораторных измерений и цифровой подход

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

Чтобы обойти это физическое ограничение, исследователи перенесли эксперимент в цифровую среду. В качестве материалов были выбраны два коммерческих сорта кофе с немного разной степенью обжарки: руандийский и колумбийский. Зерна измельчали на профессиональной кофемолке, используя 11 различных настроек — от самой мелкой фракции до крупной.

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

Затем образцы сухого молотого кофе помещали в специальные цилиндры и сканировали с помощью рентгеновской компьютерной микротомографии (XCT). Этот метод позволил создать высокодетальные трехмерные цифровые модели кофейных структур с разрешением до нескольких микрометров на воксель.

Получив 3D-модели, физики применили алгоритм решетчатых уравнений Больцмана. Это вычислительный метод гидродинамики, который моделирует поведение жидкости не как сплошной среды, а как совокупности микроскопических пакетов массы, перемещающихся по узлам дискретной решетки и сталкивающихся друг с другом. Запустив виртуальный поток жидкости через оцифрованные кофейные структуры, ученые смогли рассчитать перепады давления и точную скорость движения жидкости в каждой точке пространства.

Геометрия каналов: почему классические уравнения требуют корректировки

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

Исторически для расчета проницаемости пористых сред применяется уравнение Козени — Кармана. Эта формула устанавливает зависимость между объемом пустот, размером частиц и сопротивлением потоку. Уравнение отлично работает для материалов, состоящих из гладких сферических элементов. Однако кофейные гранулы после помола не являются сферами. Они имеют угловатую, фрагментированную форму.

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

Анализ показал, что классическое уравнение Козени — Кармана предсказывает слишком высокую проницаемость для кофе, поскольку не в полной мере учитывает специфику его изломов. Гораздо более точные результаты показала теория перколяции. Перколяционная модель математически описывает процесс формирования и исчезновения сквозных путей в среде. Она включает в себя понятие «порога протекания» — такого состояния упаковки, при котором частицы смыкаются настолько плотно из-за своей сложной геометрии, что проницаемость материала падает до нуля. Применив теорию перколяции и адаптировав расчеты удельной площади поверхности под несферические объекты, ученые получили уравнение, которое идеально совпало с результатами цифрового жидкостного моделирования.

Ламинарный и инерционный режимы потока

Еще один вопрос, который решила группа исследователей, касался характера самого потока жидкости при приготовлении эспрессо.

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

Используя полученные данные о проницаемости кофе, ученые рассчитали число Форхгеймера для типичных условий экстракции (порция около 40 граммов напитка за 30 секунд при давлении помпы 9 бар). Число Форхгеймера в пористых средах позволяет определить момент, когда силы инерции начинают доминировать над силами вязкого трения.

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

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

Практическое применение математической модели

Разработанная и валидированная физиками математическая модель доказывает, что проницаемость кофейной таблетки — это не строго детерминированная система, зависящая преимущественно от двух переменных: среднего размера частиц и связанной пористости.

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

Источник: RSOS