МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА КРИСТАЛЛИЗАЦИИ С ЭЛЕМЕНТАМИ АНАЛИЗА НАДЕЖНОСТИ ОСНОВНОГО ПРОЦЕССА
Применение нами стохастического системного подхода к анализу химико-технологического процесса и оборудования, в целом как физической системы, позволяет определить все внутренние связи между элементарными процессами, выявить определенные закономерности их протекания или же доказать полное отрицание их детерминированности. Этот метод позволяет дать математическое описание надежной работы процесса и оборудования. Получающаяся система уравнений включает обычно нелинейные интегро-дифференциальные и дифференциальные уравнения в частных производных.
Решение полученной системы требует задания начальных и граничных условий, а это возможно лишь при изучении элементарных актов процесса. Наличие математической модели с элементами теории надежности в общем случае позволяет с достаточной достоверностью перенести результаты лабораторных исследований вплоть до промышленных условий. Все вышеизложенное дает возможность своевременно получить исходный материал для расчетов с учетом внутренних связей и определить главные условия в технической реализации процесса.
За основу построения модели нами было взято уравнение скорости выхода годной продукции (производительности). В общем виде эта зависимость описывается следующим уравнением:
(1)
где J - скорость зародышеобразования;
- скорость роста кристаллов;
λагл - скорость роста кристаллов за счет агломерации;
λдробл - скорость дробления кристаллов.
В понятие надежности процесса кристаллизации и работы кристаллизационного оборудования кроме стабильной производительности необходимо включить также и такие показатели: 1. определенное процентное содержание требуемого грансостава кристаллического продукта, 2. определенная чистота продукта, иначе минимальный процент примесей, 3. оптимальные (и даже min) расходы электроэнергии, тепла, холода, растворителя и т.д.
Скорость зародышеобразования - величина весьма сложная и зависящая более чем от десятка факторов. На нее оказывают влияние физико-химические свойства вещества и растворителя, параметры процесса, режим работы установки и ее конструктивные особенности.
Скорость зародышеобразования можно описать двумя зависимостями аналитической (2) и функциональной (3)
J = K (C-C0)m (2)
где К, m - постоянные. зависящие от физико-химических свойств растворенного вещества и растворителя, например, таких как произведение валентностей ионов вещества, число молекул кристаллогидратной воды, тип кристаллической решетки и т.д.;
С0 - абсолютная растворимость вещества, или концентрация пересыщенного раствора для данных условий;
С - действительная концентрация вещества в пересыщенном растворе.
(3)
где Пс = С - Со = f (C, Со) - величина пересыщения,
Пt = (tнас - tд) = f (tнас ,tд) -величина переохлаждения,
tнас - температура насыщения раствора;
tд - действительная температура при данном пересыщении (температура раствора);
Пc (φ) = f (tр, Со, Пвал., Мквg ,КР) - устойчивость пересыщения;
tp - температурный коэффициент растворимости;
Пвал - произведение валентностей ионов составляющих соль;
Мквg - число молекул кристаллогидратной воды;
КР - тип кристаллической решетки;
t0C = tд - температура раствора;
Θпр - количество и физико-химические свойства примесей, снижающих качество продукта, а, следовательно, влияющих на надежность процесса;
Vпер - скорость перемешивания.
Скорость роста кристаллов является функцией устойчивости пересыщения, скорости перемешивания раствора и его температуры, скорости инкрустации, если она имеет место, а также количества и свойств примесей.
Скорость роста кристаллов описывается в общем виде следующей зависимостью:
(4)
От величины устойчивости пересыщения будет зависеть форма образующихся кристаллов Фк, которая является одним из показателей качества получаемого продукта.
Скорость перемешивания Vпер раствора как один из параметров процесса кристаллизации, влияет на скорость диффузии вещества из раствора к граням кристалла Д, а это, в свою очередь, влияет как на скорость роста, так и на форму кристаллов.
Л.А. Разумовский использует в своей работе [2] для описания кинетики роста среднего кристалла следующее уравнение:
(5)
где А - коэффициент времени, определяемый режимом работы кристаллизатора. Следовательно, можно считать, что коэффициент А в определенной мере хаpaктеризует надежность процесса кристаллизации.
Взаимосвязь определяет основное условие надежности процесса по показателю «размер кристалла».
Согласно выводам Л.А. Разумовского [2] уравнение 5 в координатах дает прямую линию, k = 3 и не зависит от режима работы кристаллизатора. Итак, один из показателей надежности процесса кристаллизации и работы кристаллизационного оборудования можно всегда поддерживать на должном уровне, зная время роста кристаллов до требующегося размера.
В непосредственной близости от стенки аппарата (кристаллизатора, кристаллорастителя, трубопроводов) основной процесс идет несколько иначе (медленнее), а это позволяет побочным процессам развиваться быстрее. Этим и обусловливается процесс инкрустации, который снижает надежность основного процесса кристаллизации, т.е. замедляет скорость роста кристаллов, изменяет их форму, увеличивает наличие примесей в готовом продукте (товарных кристаллах), вследствие перехода полезной части раствора в инкрустационный слой и т.п.
Уравнение, описывающее процесс инкрустации, запишется в виде функциональной зависимости:
Fинкр = f (Пt ,Vпер ,Qпр ,Ст.ф ) (6)
где Ст.ф - комплекс величин, которые определяют пристеночный эффект.
Ст.ф = φ (Δt, x, ξ) (7)
Величина Ст.ф должна непосредственно отражать теплопроводность через слой l (пристеночаая зона), в котором проходит инкрустация.
Δt - градиент температур в пристеночной зоне;
х - координата фронта охлаждения;
ξ - положение границы раздела фаз;
δ - толщина стенки;
tд > tвоздуха; t0С = tж = t2; tст = tт = t1
Класс задач, связанный с фазовыми превращениями с выделением тепла или его поглощением (это хаpaктерно для процесса ивкрустации), относится к типу задач Стефана. Такая задача формулируется как задача о сопряженности температурных полей двух фаз (1 - твердой, 2 - жидкой) при наличии особого граничного условия на движущейся границе раздела фаз (уравнения 12, 13). Решение приведенной задачи (8 - 13) в данной математической постановке предложено В.Т. Меламедом [1], который учитывал теплоту фазовых превращений в виде непрерывно распределенных на фронте кристаллизации источников тепла. Это особый метод решения, т.к. условие на границе фаз (8, 9) относит задачу Стефана к классу нелинейных и решить её обычный путем не представляется возможным. Уравнение теплопроводности:
Дополнительные условия:
где а1, а2 - температуропроводность соответственно твердого и жидкого слоев;
λ- теплопроводность;
Q - количество теплоты кристаллизации в 1 м3;
т - время кристаллизации.
В функциональной зависимости (1) одним из аргументов роста кристаллов является λагл, т.е. параллельно основному процессу - росту кристаллов, увеличение их размеров может происходить еще за счет агломерации кристаллов.
Агрегированию могут подвергаться как зародыши, так и кристаллы любых размеров, разница лишь в том, что скорость агрегирования будет зависеть, в частности, от размеров контактируемых элементов.
В общем виде скорость процесса агломерации можно выразить следующим уравнением:
(14)
где Ипер - интенсивность перемешивания;
т - время кристаллизации (среднее время пребывания кристалла в слое);
Vt - скорость охлаждения раствора;
dо, V - размер (средний диаметр, объем) соответственно зародыша и кристалла;
σ - поверхностные силы натяжения;
Скр - концентрация кристаллов в слое.
На скорость роста кристаллов будет оказывать нежелательное воздействие дробление кристаллов, которое будет зависеть в общем виде от следующих величин:
λдробл = f (Vпер , Vt , V, Скр ,T) (15)
Рассматривать дробление по отношению к качеству продукта, как процесс уравновешивающий агломерацию, нельзя, вследствие того, что образующиеся осколки, как и образующиеся агломераты, снижают качество готового продукта, увеличивая неоднородность размеров.
Итак, анализируя уравнение скорости выхода годной продукции (1) можно установить следующее:
- Не все величины, входящие в это уравнение, оказывают положительное воздействие на процесс кристаллоращения и степень их влияния не одинакова.
- Такие факторы, как наличие примесей Qпр, протекание процессов дробления λдр и инкрустации Fинкр оказывают полное отрицательное воздействие на основной процесс, то есть снижают его надежность пропорционально величинам Qпр, λдр, Fинкр.
- Такие параметры, как Мквд, t°С, Vпер, т, Скр и процесс агломерации Хагл, могут оказывать нежелательное воздействие на основной процесс, если не поддерживать их на оптимальном уровне.
Литература.
- Меламед В.Т. Сведение задачи Стефана к системе обыкновенных дифференциальных уравнений. Известия АН СССР. Серия геофизика, 1958. № 7.
- Разумовский Л.А. Исследование кинетики процесса кристаллизации солей из растворов во взвешенном слое. 1967. Кандидатская диссертация. г. Иваново. ИХТИ.
Статья в формате PDF 130 KB...
26 04 2024 17:29:27
Статья в формате PDF 187 KB...
25 04 2024 19:41:12
Статья в формате PDF 474 KB...
24 04 2024 10:12:34
В статье рассматриваются две разновидности оттепели изменение глубины путем восстановления этапов нарушенных ландшафтов вечной мерзлоты, которые функционируют на суглинистых и песчаных отложениях высоких террас на правом и левом берегах реки Лены. Качественные изменения в динамике глубины сезонного оттаивания был обнаружен в определенные промежутки времени сукцессии этапов: трава, кустарники, березы, лиственницы (сосна) – березы и лиственницы (сосна). ...
23 04 2024 13:40:48
Статья в формате PDF 106 KB...
22 04 2024 12:25:21
Статья в формате PDF 119 KB...
21 04 2024 18:36:41
Статья в формате PDF 141 KB...
20 04 2024 22:55:25
Статья в формате PDF 170 KB...
18 04 2024 18:53:53
Территориальные различия электopaльных предпочтений отличаются высокой устойчивостью в современной России. Этот феномен подтверждается методом корреляционного анализа. Выделяются шесть основных социальных факторов, влияющих на различия в электopaльной географии: 1) доля городского населения; 2) приближенность к центру; 3) этнический фактор; 4) доля молодежи в составе населения; 5) преобладающие виды деятельности населения; 6) структура социальных связей. Электopaльное поведение в России менее индивидуально, чем в западных странах, большее значение имеют объективные социальные факторы. ...
16 04 2024 21:22:47
Статья в формате PDF 279 KB...
15 04 2024 16:52:44
Статья в формате PDF 184 KB...
14 04 2024 13:48:26
В этой статье рассматриваются вопросы применения инновационных технологии формирования навыков здорового образа жизни учащихся сельских школ во внеурочное время. ...
12 04 2024 2:54:40
Статья в формате PDF 257 KB...
11 04 2024 7:49:34
Статья в формате PDF 123 KB...
09 04 2024 6:26:21
08 04 2024 0:47:44
При анализе количества видов гельминтов (возбудителей зоонозов) у человека в 1999–2012 гг. увеличилось с 7 до 10 видов (на 30 %), в том числе цестод с 3 до 5 видов (на 40 %) и нематод с 4 до 5 видов (на 20 %). У человека и собак прослеживается биологический прогресс возбудителей зоонозов. Количества видов гельминтов у собак увеличилось с 5 до 8 видов (на 37,5 %), в т.ч. цестод с 2 до 3 видов (на 33,3 %) и нематод с 3 до 5 видов (на 40 %). В составе гельминтофауны общих для человека и животных доминировали классы Nematoda (6 видов) и Cestoda (5 видов) над классом Trematoda (3 вида). ...
07 04 2024 11:44:51
Статья в формате PDF 257 KB...
06 04 2024 19:37:24
Статья в формате PDF 148 KB...
05 04 2024 7:59:51
Статья в формате PDF 206 KB...
04 04 2024 11:16:31
Статья в формате PDF 109 KB...
02 04 2024 10:11:40
Статья в формате PDF 110 KB...
01 04 2024 17:31:46
Статья в формате PDF 253 KB...
31 03 2024 20:38:46
Статья в формате PDF 120 KB...
30 03 2024 7:18:37
Статья в формате PDF 112 KB...
29 03 2024 0:12:41
Статья в формате PDF 142 KB...
28 03 2024 23:11:25
Статья в формате PDF 109 KB...
27 03 2024 8:50:27
Статья в формате PDF 118 KB...
26 03 2024 20:18:43
Статья в формате PDF 108 KB...
25 03 2024 7:18:26
Статья в формате PDF 112 KB...
22 03 2024 10:26:15
Статья в формате PDF 257 KB...
21 03 2024 6:39:56
Статья в формате PDF 109 KB...
20 03 2024 6:10:33
19 03 2024 14:38:45
Статья в формате PDF 248 KB...
18 03 2024 14:11:30
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::