Гидродинамическая модель измельчения кальцийсодержащих пищевых добавок

УДК 621.047:621.926

Гидродинамическая модель измельчения кальцийсодержащих пищевых добавок

асп. Даниленко Е.А., проф. Алексеев Г.В.

Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий, gva2003@rambler.ru

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

Ключевые слова: измельчение, пищевая добавка, функциональное питание, яичная скорлупа, модель кавитационно-истирающего процесса.

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

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

Таким образом, представляется актуальным изучение вопроса получения более мелких фракций яичной скорлупы при использовании аналитических моделей, например, гидродинамического разрушения ее фрагментов, предварительно подвергнутых механическому воздействию [1].

При предварительном измельчении фрагменты яичной скорлупы под действием механических воздействий приобретают множественные трещины-

1

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

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

Y _ q

f x2 у2 ^

— + Т2 ^ a b j

где ¥ — функция тока;

х, у — текущие координаты жидкой частицы в декартовой системе; a, b — большая и малая, соответственно, полуоси эллипсов, по которым движется жидкая частица; q — постоянный параметр.

Составляющую окружной скорости вдоль оси Ох запишем в виде:

UC

у,

dy _ 2q

а составляющую вдоль оси Оу определим как:

U __dy__2lx

dx a

В этом случае проекции силы равномерного движения жидкости по эллиптическим траекториям на оси координат можно представить так:

2q

b

TLT

Nx _k • ux _k • — у _ cxy,

д т 2q

Ny _k2 • Uy _ k ~x _ C2x ,

a

где k, k2, Cj, c2 — постоянные коэффициенты.

На самом деле, поскольку рассматриваемое движение жидкости в дефекте является неравномерным в течение периода её обращения, за счёт возникающих при этом сил сопротивления, следует учесть внутреннюю силу сопротив-

2

ления движению жидкости в каверне. Она возникает от влияния вязкости, благодаря которой главным образом происходит диссипация энергии, накопленной жидкостью в каверне от транзитного потока.

Эту силу сопротивления логично принять пропорциональной скорости движения жидкости в каверне и представить в виде:

T = -АУп,

где 1 — некоторый коэффициент пропорциональности;

V — линейная скорость неравномерного движения жидкой частицы вдоль эллиптических, но уже неконцентрических траекторий; п — показатель степени, который может изменяться в пределах от 1 до 2.

Составляющие силы сопротивления Т1 вдоль координатных осей теперь запишутся следующим образом:

- вдоль оси Ох

dx

А1Х

-1

dt

ушу

- вдоль оси Оу

т„ = -1

dyл

dt

ушу

Знак «минус» указывает на то, что силы сопротивления направлены в противоположном направлении скорости движения жидкости.

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

dx

движения жидкости в дефекте вдоль оси Ох —.

dt

Учитывая эти соображения имеем:

{

T2 =j

V

dx U ± dt

j

dx

Знак «±» перед производной — показывает, что в общем случае скорость

dt

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

п

3

тенденция к отставанию скорости движения жидкости в дефекте против скорости транзитного потока.

Составим динамические уравнения движения жидкой частицы в выемке, предположив, что весом частицы можно пренебречь.

Тогда имеем:

К

X

К

-T\X + NX + T2

Ty + Ny

где через FX и Fy обозначены проекции равнодействующей силы на координатные оси.

Перепишем систему (1), подставляя значения найденных ранее сил и их проекций, а также учитывая, что F = mw, где m — масса жидкой частицы, а w — её ускорение. В этом случае получим систему уравнений:

d2 x

т-

dt2

+ 1

( dx 4 n f dx 4

v dt у - С1У = j U ~~F

V dt J

d2 У + о m—— +1 dt2

dy

v dt j

+ c2 x = 0

Полученная система представляет собой неоднородную нелинейную систему двух обыкновенных дифференциальных уравнений второго порядка.

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

Учитывая, что в общем случае составляющие силы сопротивления T вдоль координатных осей могут зависеть от степени механического воздействия на первоначальные фрагменты, следует дополнительно предположить зависимость коэффициента X в уравнениях системы (2) соответственно в первом уравнении от х, а во втором от у.

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

В связи с этим для решения системы записанных дифференциальных уравнений можно рекомендовать метод Булирша-Штера, который не только более

4

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

Полученные результаты отражены на рисунке.

Статистическая обработка полученных численных значений методом наименьших квадратов позволила записать для приведенных решений следующие уравнения регрессии:

dx/dt = 1.77 - 15.01? + 2.16t2, dy/dt = 0.27 - 2.85? + 2.06?2.

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

77 mV2

F =---.

r

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

Рис. Зависимость dx/dt (Y2) и dy/dt (Y4) от времени (Yo).

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

Экспериментальное значение содержания кальция в готовом полуфабрикате (на примере бисквита основного) в 12-15 раз выше расчетного. А это означает, что смесь муки со скорлупой яиц в процессе приготовления бисквитного

5

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

Список литературы

1. Даниленко Е.А., Алексеев Г.В. Совершенствование использования минеральных добавок в производстве бисквита, Тезисы докладов Х Международной конференции молодых ученых «Пищевые технологии и биотехнологии» (13-16 апреля 2010 года), Казань, 2010

2. Алексеев Г.В., Кондратов А.В., Верболоз Е.И. Модель кавитационного разрушения пищевого сырья. Хранение и переработка сельхозсырья, №11, 2007, с. 56-59

Hydrodynamic model of pounding calcium containing

food additives

Danilenko E.A., Alexeyev G.V.

Saint-Petersburg State University of Refrigeration and Food Engineering

Within the frames of works carried out by various investigators positive effect of calcium containing food additives on consumer properties of functional foods became manifest.

Keywords: milling, food additive, functional foods, eggshell, model of cavita-tional abrasion.

6