Теоретические основы химической технологии, 2022, T. 56, № 6, стр. 748-759
Математическое описание процесса гранулирования расплавов на охлаждаемых поверхностях
a МИРЭА – Российский технологический университет
(Институт тонких химических технологий имени М.В. Ломоносова)
Москва, Россия
b Московский политехнический университет
Москва, Россия
* E-mail: taran_yu@mirea.ru
Поступила в редакцию 11.08.2022
После доработки 15.08.2022
Принята к публикации 15.08.2022
- EDN: KVPUYI
- DOI: 10.31857/S0040357122060173
Аннотация
Для экологически безопасных, ресурсо-энергосберегающих процессов гранулирования кристаллизацией капель расплавов на охлаждаемых поверхностях предложено математическое описание. Оно включает процессы кристаллизации, энантиотропные полиморфные превращения, охлаждение, с учетом переменных по поверхности гранул условий отвода теплоты, усадочные явления, кинетические параметры и переменные коэффициенты модели, а также учитывает в отличие от ранее предложенных нами моделей не только радиальный, но и продольный перенос теплоты и динамику межфазных границ. Усовершенствован итерационный локально-одномерный алгоритм численного решения с подвижными узлами на межфазных границах на шести точечной, неявной, нелинейной, устойчивой консервативной конечно разностной схеме. Методом скоростной термографии подтверждена адекватность модели ходу реального процесса. В качестве научно-практической иллюстрации предложенного проведены вычислительные эксперименты и анализ работы ленточных, валковых, дисковых кристаллизаторов применительно к ряду промышленных продуктов: азотсодержащих минеральных удобрений, щелочей, серы, нитрата калия, ацетонанила, шлаков фосфорных печей и других. Предложенное математическое описание в рамках формальной аналогии процессов со структурной перестройкой исходной системы может быть применено как для “классических” фазовых превращений (кристаллизация–плавления, энантиотропные полиморфные превращения) так и для некоторых формально аналогичных им процессов таких, как переходы гель–золь.
ВВЕДЕНИЕ
Наиболее экологически безопасным процессом гранулирования расплавленных веществ является кристаллизация их капель на охлаждаемых поверхностях. Процесс можно организовать с замкнутой по хладоагенту схемой, утилизировав теплоту гранулирования с использованием тепловых насосов или других способов “съема” теплоты [1]. Включение в технологическую схему теплового насоса более эффективно при гранулировании из расплавов веществ, кристаллизующихся и претерпевающих энантиотропные полиморфные превращения при высоких (~200–400°С) температурах (нитрат калия, щелочи, аммиачная селитра, шлаки и другие) [1, 2]. Важными преимуществами данной технологии являются реализация замкнутой по потоку хладоагента схемы, утилизация теплоты превращения, стабильное качество гранулированного продукта (высокая статическая прочность, гладкая поверхность, монодисперсность, обеспечиваемая применением современного оборудования, например, ротационных диспергаторов фирмы “Sandvik”) [1–5]. Для токсичных, технологически опасных, высокотемпературных малотоннажных продуктов процесс безальтернативен [5–7]. Используемые в промышленности валковые (все), дисковые (в основном) и ленточные (иногда) аппараты выпускаются с замкнутой полостью для хладоагента, что позволяет легко организовать замкнутый цикл по хладоагенту. “Съем” теплоты, превращения и отсутствие контакта продукта с хладагентом обеспечивает экологическую и технологическую безопасность [1]. Эти аппараты являются многофункциональными, в промышленности они широко использовались, как кристаллизаторы для получения чешуированного продукта (при наличии “насечки”) в виде плиток [8, 9]. Так же они используются, как грануляторы для кристаллизации капель расплавов на охлаждаемых поверхностях [1, 2, 5–7, 9]. Известны способы использования этих аппаратов для глубокой очистки веществ по технологии непрерывной зонной плавки [10–12] или фракционирования совмещенного с гранулированием.
МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА
Реализация таких разработок требует детального рассмотрения процессов кристаллизации полиморфных превращений, охлаждения, усадочных явлений, динамики межфазных границ и температурных полей при переменных по поверхности и во времени условиях отвода теплоты от лежащей на охлаждаемой поверхности несферической гранулы. Тепло – и массообменные процессы в подобных системах и их математическое описание рассмотрены в [13–16]. Начиная расчет процесса, каплю расплава считали сферической без сегмента, объем которого определяли по краевому углу смачивания φсм расплавом охлаждающей поверхности. Полагали, что условия переноса теплоты между кристаллизующейся каплей расплава и охлаждающей поверхностью кондуктивные (граничные условия IV рода), а “свободная” поверхность капли обменивается теплотой с окружающей средой конвективно (в граничных условиях III рода). С учетом симметрии задачи:
(1)
$\begin{gathered} {{c}_{j}}\left( t \right){{\rho }_{j}}\left( t \right)\frac{{\partial t\left( {{{r}_{i}},{{\theta }_{{i,}}}\tau } \right)}}{{\partial \tau }} = \,\,~\frac{1}{{r_{i}^{2}}}\left[ {\frac{\partial }{{\partial {{r}_{1}}}}\left( {{{\lambda }_{i}}\left( t \right)r_{i}^{2}\frac{{\partial t\left( {{{r}_{{i,}}}{{\theta }_{{i,}}}\tau } \right)}}{{\partial {{r}_{i}}}}} \right)} \right] + \\ + \,\,\frac{1}{{r_{i}^{2}\sin {{\theta }_{i}}}}\left[ {\frac{\partial }{{\partial {{\theta }_{i}}}}\left( {{{\lambda }_{j}}\left( t \right)\sin {{\theta }_{i}}\frac{{\partial t\left( {{{r}_{i}},{{\theta }_{i}},\tau } \right)}}{{\partial {{\theta }_{i}}}}} \right)} \right], \\ 0 \leqslant \theta \leqslant \pi ;\,\,\,\,0 \leqslant {{r}_{1}} \leqslant {{R}_{i}}. \\ \end{gathered} $Уравнение переноса теплоты в плоской стенке аппарата:
(2)
$\begin{gathered} \frac{\partial }{{\partial x}}\left[ {{{\lambda }_{{{\text{ст}}}}}\left( t \right)\frac{{\partial t\left( {x,z,\tau } \right)}}{{\partial x}}} \right] + \frac{\partial }{{\partial z}}\left[ {{{\lambda }_{{{\text{ст}}}}}\left( t \right)\frac{{\partial t\left( {x,z,\tau } \right)}}{{\partial z}}} \right] - \\ - \,\,{{c}_{{{\text{ст}}}}}\left( t \right){{\rho }_{{{\text{ст}}}}}\left( t \right)\frac{{\partial t\left( {x,z,\tau } \right)}}{{\partial \tau }} = 0. \\ \end{gathered} $Начальные условия
(3)
$\begin{gathered} t({{r}_{i}}_{,}{{\theta }_{i}},0) = {{t}_{{{\text{жн}}}}} = {\text{const}};\,\,\,\,0 \leqslant {{r}_{i}}_{~} \leqslant {{R}_{i}}; \\ {\text{ }}0 \leqslant {{\theta }_{i}} \leqslant \pi ;\,\,\,\,{{r}_{\xi }}{{_{,}}_{j}}{{_{,}}_{i}}({{\theta }_{i}},0) \equiv {{R}_{i}} + 0;\,\,\,\,{{r}_{{{\text{П}},}}}_{i}({{\theta }_{i}},0) \equiv 0; \\ t(x,z,0) = t(x,z,{{\tau }_{{{\text{охл}}}}}){\text{ }} \vee {{t}_{c}}. \\ \end{gathered} $Условия начала фазового превращения (игнорируя не сферичность капли при контакте со стенкой):
(4)
$\begin{gathered} 2\pi \int\limits_{{{\tau }_{{j,i}}}}^{{{\tau }_{{{\text{инд}}j,i}}}} {\int\limits_{{{r}_{{j,i}}}}^{{{R}_{i}}} {\int\limits_0^{{{\theta }_{{j,i}}}} {{{w}_{{{\text{з}}j}}}\left( {t\left( {{{r}_{i}},{{\theta }_{i}},\tau } \right)} \right)r_{i}^{2}\sin {{\theta }_{i}}d{{\theta }_{i}}d{{r}_{r}}d\tau } } } = 1; \\ t\left( {{{r}_{{j,i}}},{{\theta }_{{j,i}}},{{\tau }_{{j,i}}}} \right) = {{t}_{{{\text{ф}}j}}}, \\ \end{gathered} $(5)
$t({{x}_{{k,i}}},{{\theta }_{{k,i,}}}{{\tau }_{{{\text{инд}}{\kern 1pt} }}}{{_{{\text{к}}}}_{{,i}}}) = {{t}_{{\text{к}}}}.$Условие на границе сосредоточенной усадочной поры радиусом rПj,i
(6)
$\begin{gathered} \overline {{\text{grad}}} {\kern 1pt} t({{r}_{{\text{П}}}}_{{j,i}},{{\theta }_{{j,i}}},\tau ) = 0; \\ {{V}_{{{\text{П}}i}}}\left( \tau \right) = {{V}_{{{\text{к}}i}}}\left( \tau \right)\frac{{{{\rho }_{{\text{к}}}}\left( {{{t}_{{{\text{фк}}}}}} \right) - {{\rho }_{{\text{ж}}}}\left( {{{t}_{{{\text{фк}}}}}} \right)}}{{{{\rho }_{{\text{ж}}}}\left( {{{t}_{{{\text{фк}}}}}} \right)}}; \\ \frac{{{{S}_{{{\text{П}}i}}}\left( \tau \right)\sigma \left( \tau \right)}}{{{{V}_{{{\text{П}}i}}}\left( \tau \right)}} \to {\text{мин}}. \\ \end{gathered} $Объем распределенной усадочной пористости VПР,j,i(τ) в виде трещин, возникающих при термической усадке:
(7)
$\begin{gathered} {{V}_{{{\text{ПР}}j,i}}}\left( \tau \right) = {{V}_{{j,i}}}\left( \tau \right)\frac{{\overline {{{\rho }_{j}}} - {{\rho }_{j}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}{{\overline {{{\rho }_{j}}} }}; \\ \overline {{{\rho }_{j}}} = \frac{1}{{{{V}_{{j,i}}}\left( \tau \right)}}2\pi \int\limits_{{{r}_{{\zeta j + 1,i}}}}^{{{r}_{{\zeta j,i}}}} {\int\limits_0^\pi {{{\rho }_{j}}\left( {{{r}_{i}},{{\theta }_{i}},\tau } \right)r_{i}^{2}\sin {{\theta }_{i}}d{{\theta }_{i}}d{{r}_{i}}} \, = \,f(\tau )} ; \\ {{\rho }_{j}} = f\left( t \right);\,\,\,\,t = \left( {{{r}_{i}},{{\theta }_{i}},\tau } \right). \\ \end{gathered} $Объем распределенной усадочной пористости, возникающей при полиморфном превращении j – j + 1 (суммируется по всем 1 – J превращениям, которые претерпевает кристаллическая фаза):
(8)
$V_{{{\text{ПР}}j + 1,i}}^{'}\left( \tau \right) = {{V}_{{j + 1,i}}}\left( \tau \right)\frac{{{{\rho }_{{j + 1}}}\left( {{{t}_{{{\text{ф}}j}}}} \right) - {{\rho }_{j}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}{{{{\rho }_{j}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}.$Условия на межфазных (кристаллизации и энантиотропных полиморфных превращений) границах:
(9)
${{t}_{j}}\left( {{{r}_{{\xi j,i}}}{{\theta }_{{j,i}}},\tau } \right) = {{t}_{{j + 1}}}\left( {{{r}_{{\xi j,i}}}{{\theta }_{{j,i}}},\tau } \right) < {{t}_{{{\text{ф}}j}}}.$Граница раздела фаз смещается по нормали к границе со скоростью:
(10)
${{v}_{{\Lambda j}}} = \frac{{d{{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}{{\bar {n}}_{{ji}}}$Зная вектор $\overline {{\text{grad}}{\kern 1pt} } t({{r}_{{j,i}}} = {{\xi }_{{j,i}}}(\theta ,\tau ),{{\theta }_{{j,i}}},\tau )$ в этой области фазовой границы:
(11)
$\begin{gathered} \overline {{\text{grad\;}}} {\kern 1pt} t({{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i}}},\tau ) = \\ = \frac{{\partial t\left( {{{r}_{{\xi j,i}}},{{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{r}_{i}}}}{{{\bar {i}}}_{r}} + \frac{1}{{{{r}_{i}}}}\frac{{\partial t\left( {{{r}_{{\xi j,i,}}}{{\theta }_{{j,\tau }}}} \right)}}{{\partial {{\theta }_{i}}}}{{{\bar {i}}}_{\theta }}, \\ \end{gathered} $(12)
${{\bar {n}}_{{j,i}}} = \frac{{\frac{{\partial t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,\tau }}}} \right)}}{{\partial {{r}_{i}}}}{{{\bar {i}}}_{r}} + \frac{1}{{{{r}_{i}}}}~\frac{{\partial t\left( {{{r}_{{i,j}}} = {{\xi }_{{i,j}}},{{\theta }_{{j,i}}},\tau } \right)}}{{\partial {{\theta }_{i}}}}{{{\bar {i}}}_{\theta }}}}{{{{{\left[ {{{{\left( {\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i}}}\tau } \right)}}{{\partial {{r}_{i}}}}} \right)}}^{2}} + {{{\left( {\frac{1}{{{{r}_{i}}}}~\frac{{\partial t\left( {{{r}_{{i,j}}} = {{\xi }_{{i,j}}}\left( {\theta ,\tau } \right),{{\theta }_{{i,j,}}}\tau } \right)}}{{\partial {{\theta }_{i}}}}} \right)}}^{2}}} \right]}}^{{1/2}}}}} = \frac{{\overline {{\text{grad}}{\kern 1pt} ~t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji}}},\tau } \right)~}}{{\left| {{\text{grad}}{\kern 1pt} ~t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji}}},\tau } \right)} \right|}}.$Косинус угла между вектором скорости перемещения межфазной границы ${{\overline {d\xi } }_{{j,i}}}$(θ, τ)/dτ: ${{\bar {v}}_{{\Lambda ,j}}}$ = $ = \frac{{d{{{\bar {\xi }}}_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }} = $ $\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial \tau }}{{\bar {i}}_{\xi }} + $ $\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial {{\theta }_{{j,i}}}}}~\frac{{\partial {{\theta }_{{j,i}}}}}{{\partial \tau }}{{\bar {i}}_{\theta }}$ и единичным вектором нормали nj,i к межфазной поверхности j – j + 1 равен:
(13)
$\begin{gathered} \cos \left( {{{{\bar {n}}}_{{j,i}}}\overline {\frac{{d{{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}} } \right) = \frac{1}{{{{{\left( {{{{\left( {\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i.}}}\tau } \right)}}{{\partial {{r}_{i}}}}} \right)}}^{2}} + {{{\left( {\frac{1}{{{{r}_{i}}}}\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i}}},\tau } \right)}}{{\partial {{\theta }_{i}}}}} \right)}}^{2}}} \right)}}^{{1/2}}}}} \times \\ \times \,\,\frac{1}{{{{{\left( {{{{\left( {\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial \tau }}~} \right)}}^{2}} + ~{{{\left( {\frac{{\partial \xi \left( {{{\theta }_{i}}\tau } \right)}}{{\partial {{\theta }_{i}}}}\frac{{\partial {{\theta }_{i}}}}{{\partial \tau }}} \right)}}^{2}}} \right)}}^{{1/2}}}}} \times \\ \times \,\,\left( {\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial \tau }}~\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{r}_{i}}}}{{{\bar {i}}}_{r}} + \frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial {{\theta }_{{j,i}}}}}~\frac{{\partial {{\theta }_{{j,i}}}}}{{\partial \tau }}~\frac{1}{{{{r}_{i}}}}~\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{\theta }_{i}}}}{{{\bar {i}}}_{\theta }}} \right). \\ \end{gathered} $Проекция ${{\bar {v}}_{{\Lambda ,j}}}$ линейной скорости движения фронта на единичный вектор нормали к этой поверхности ${{\bar {n}}_{{j,\tau }}}$ равна:
(14)
$\begin{gathered} {{v}_{{\Lambda j}}} = ~\,\,\frac{{d{{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}{{{\bar {n}}}_{{j,i}}}~\,\, = {{\left( {{{{\left( {\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial \tau }}} \right)}}^{2}} + {{{\left( {\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial {{\theta }_{{j,i}}}}}~\frac{{\partial {{\theta }_{{j,i}}}}}{{\partial \tau }}} \right)}}^{2}}} \right)}^{{\frac{1}{2}}}}~\cos \left( {\bar {n},~\frac{{d{{{\bar {\bar {\xi }}}}_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}} \right) \times \\ \times \,\,\frac{1}{{{{{\left( {{{{\left( {\frac{{\partial t\left( {{{r}_{{j,i = {{\xi }_{{j,i}}}}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{r}_{i}}}}} \right)}}^{2}} + {{{\left( {\frac{1}{{{{r}_{i}}}}\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{\theta }_{i}}}}} \right)}}^{2}}} \right)}}^{{\frac{1}{2}}}}}} \times \\ \times \,\,\left( {\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}~\frac{{\partial t\left( {{{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right){{\theta }_{{j,i,}}}\tau } \right)}}{{\partial {{r}_{i}}}}{{{\bar {i}}}_{r}} + ~\,\,\frac{{\partial {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{\partial {{\theta }_{{j,i}}}}}~\frac{{\partial {{\theta }_{{j,i}}}}}{{\partial \tau }}~\frac{1}{{{{r}_{i}}}}~\frac{{\partial t({{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{j,i,}}}\tau }}{{\partial {{\theta }_{i}}}}{{{\bar {i}}}_{\theta }}} \right) \times \\ \times \,\,\frac{1}{{~\left| {{\text{grad}}{\kern 1pt} ~t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right|}}~\left( {\frac{{\overline {d{{\xi }_{{j,i}}}} \left( {\theta ,\tau } \right)}}{{d\tau }}\overline {{\text{grad}}~\,t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right).~~ \\ \end{gathered} $Сам вектор перемещения границы раздела фаз по нормали к поверхности межфазной границы ${{\bar {v}}_{{\Lambda ,j}}}$ равен произведению проекции линейной скорости движения фронта ${{v}_{{\Lambda ,j}}}$ (скалярная величина) на единичный вектор нормали к межфазной границе ${{\bar {n}}_{{j,i}}}$:
(15)
$\begin{gathered} {{v}_{{\Lambda j}}} = \,\,~\frac{{d{{\xi }_{{j,i}}}\left( {\theta ,\tau } \right)}}{{d\tau }}~{{{\bar {n}}}_{{j,i}}}~{{{\bar {n}}}_{{j,i}}} = \frac{1}{{~\left| {{\text{grad}}{\kern 1pt} ~t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right|}}~\left( {\frac{{\overline {d{{\xi }_{{j,i}}}} \left( {\theta ,\tau } \right)}}{{d\tau }}\overline {{\text{grad}}~t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right) \times \\ \times \,\,\frac{{\overline {{\text{grad}}\,~t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)}}{{\left| {{\text{grad}}\,~t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right|}} = \frac{{\overline {{\text{grad}}\,~t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)}}{{{{{\left| {{\text{grad}}\,~t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right|}}^{2}}}}\left( {\frac{{\overline {d{{\xi }_{{j,i}}}} \left( {\theta ,\tau } \right)}}{{d\tau }}\overline {{\text{grad}}~\,t} \left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right). \\ \end{gathered} $Баланс теплоты на подвижных границах кристаллизации и полиморфных превращений:
(16)
${{\bar {q}}_{{j + 1,i}}} = {{\bar {q}}_{{j,i}}} + {{r}_{{j + 1~}}}{{\rho }_{{j + 1}}}{{\bar {v}}_{{\Lambda ,j}}},~$При кондуктивном переносе теплоты в фазах удельные потоки в них равны:
(17)
${{\bar {q}}_{{j + 1,i}}} = {{\lambda }_{{j + 1}}}\left( t \right)\overline {{\text{grad}}} ~{{t}_{{j + 1}}}({{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right){{\theta }_{{ji,}}}\tau ),$(18)
${{q}_{{j,i}}} = {{\lambda }_{j}}\left( t \right){\text{grad}}~{{t}_{j}}({{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau ).$Тогда условие Стефана на межфазных границах имеет вид:
(19)
$\begin{gathered} {{\lambda }_{{j + 1}}}\left( t \right)\overline {{\text{grad}}\,} ~{{t}_{{j + 1}}}\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),~{{\theta }_{{ji}}},\tau } \right) - {{\lambda }_{j}}\left( t \right)\overline {{\text{grad}}} ~\,{{t}_{j}}\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),~{{\theta }_{{ji}}},\tau } \right) = \\ = {{r}_{{{\text{ф}},j + 1}}}{{\rho }_{{j + 1}}}\left( {\frac{{\overline {d{{\xi }_{{j,i}}}} \left( {\theta ,\tau } \right)}}{{d\tau }}\overline {{\text{grad}}} ~{{t}_{{j + 1}}}\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),~{{\theta }_{{ji}}},\tau } \right)} \right)\frac{{\overline {{\text{grad}}} \,~{{t}_{j}}\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)~}}{{{{{\left| {{\text{grad}}~\,{{t}_{j}}\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)~} \right|}}^{2}}}} = {{v}_{{\Lambda j}}}{{r}_{{{\text{ф}},j + 1}}}{{\rho }_{{j + 1}}}; \\ \end{gathered} $(20)
${{v}_{{\Lambda j}}} = \left\{ \begin{gathered} 0,{{r}_{{\xi j,i}}}\left( {\theta ,\tau } \right) = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right) = {{R}_{i}} + 0;~\tau < {{\tau }_{{{\text{инд}}}}}; \hfill \\ {{f}_{\Lambda }}\left( {\Delta {{t}_{3}}} \right) = {{f}_{\Lambda }}\left[ {{{t}_{{{\text{ф}}j}}} - t\left( {{{r}_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji,}}}\tau } \right)} \right], \hfill \\ {{r}_{{j,i}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right) < {{R}_{i}};\tau \geqslant {{\tau }_{{{\text{инд}}}}}. \hfill \\ \end{gathered} \right.$Условие устойчивости плоской межфазной границы имеет вид:
(21)
$\begin{gathered} - {{\lambda }_{j}}\left( t \right)\overline {{\text{grad}}\,~} {{t}_{j}}\left( {{{\tau }_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji}}},\tau } \right) \leqslant \\ \leqslant C~{{\lambda }_{{j + 1}}}\left( t \right)~{{\lambda }_{{j + 1}}}\left( t \right)\overline {{\text{grad}}\,~} {{t}_{{j + 1}}}\left( {{{\tau }_{{ji}}} = {{\xi }_{{j,i}}}\left( {\theta ,\tau } \right),{{\theta }_{{ji}}},\tau } \right), \\ \end{gathered} $Условие не появления новых центров новообразований в метастабильной области перед фазовой границей имеет вид, если считать, что коэффициент теплоотдачи разный в области контакта гранулы со стенкой и на свободной поверхности гранулы:
(22)
$\begin{gathered} 2\pi ~\mathop \smallint \limits_0^{{{\tau }_{{vj}}}} \mathop \smallint \limits_{{{r}_{{j,i}}}}^{{{r}_{{\xi j,i}}}} \mathop \smallint \limits_0^{{{\theta }_{{j,i}}}} {{\omega }_{{з,j}}}\left[ {t\left( {{{r}_{i}},{{\theta }_{i}},\tau } \right)} \right]r_{i}^{2}\sin {{\theta }_{i}}d{{\theta }_{i}}d{{r}_{i}}d\tau < 1; \\ t({{r}_{{j,i,}}}{{\theta }_{{j,i,}}}\tau ) \leqslant {{t}_{{{\text{ф}},j}}};\,\,\,\,~{{\tau }_{{\upsilon ,j}}}\left( {{{\theta }_{i}}} \right)~ = \mathop \smallint \limits_{{{r}_{{j,i}}}}^{{{r}_{{\xi j,i}}}} \frac{{d{{r}_{i}}\left( {{{\theta }_{i}}} \right)}}{{{{v}_{{\Lambda j}}}}}. \\ \end{gathered} $Если считать, что поверхность контакта капли расплава со стенкой мала, а коэффициент теплоотдачи постоянен по поверхности гранулы и равен его величине для свободной поверхности, то условие не появления новых центров в метастабильной области имеет вид:
(23)
$\begin{gathered} 4\pi \mathop \smallint \limits_0^{{{\tau }_{{''j}}}} \mathop \smallint \limits_{{{r}_{{j,i}}}}^{{{r}_{{\xi j,i}}}} {{\omega }_{{з,j}}}\left[ {t\left( {{{\tau }_{i}},\tau } \right)} \right]r_{i}^{2}d{{r}_{i}}d\tau < 1; \\ t\left( {{{r}_{{j,i}}},\tau } \right) \leqslant {{t}_{{{\text{ф}},j}}};~\,\,\,\,{{\tau }_{{\upsilon ,j}}} = \mathop \smallint \limits_{{{r}_{{j,i}}}}^{{{r}_{{\xi j,i}}}} \frac{{dr}}{{{{v}_{{\Lambda ,j}}}~}}. \\ \end{gathered} $Условие на внешней границе капли имеет вид:
(24)
$\begin{gathered} {{\alpha }_{{c,i}}}\left( {{{\theta }_{{i,\tau }}},\tau } \right)\left[ {t\left( {{{R}_{{i,}}}{{\theta }_{i}},\tau } \right) - {{t}_{c}}\left( {{{r}_{{\text{Б}}}},\tau } \right)} \right] = \\ = - {{\lambda }_{j}}\left( t \right)\overline {{\text{grad}}\,} ~{{t}_{j}}\left( {{{R}_{i}},{{\theta }_{i}},\tau } \right);\,\,\,\,{{\varphi }_{{{\text{см}}}}} \leqslant {{\theta }_{{j,i}}} \leqslant \pi , \\ \end{gathered} $Условие на границе капли расплава и плоской стенки имеет вид:
(25)
$\begin{gathered} - {{\lambda }_{{{\text{ст}}}}}\left( t \right){\text{grad}}\,~{{t}_{{{\text{ст}}}}}\left( {х{\kern 1pt} ',\theta ,\tau } \right)~\,\, = ~ \\ = - {{\lambda }_{j}}\left( t \right){\text{grad}}\,~{{t}_{j}}\left( {R_{j}^{'},{{\theta }_{{j,I}}},\tau } \right);~0 \leqslant {{\theta }_{{j,i}}} \leqslant {{\varphi }_{{{\text{см}}}}}; \\ 0 \leqslant х{\kern 1pt} ' \leqslant {{R}_{i}}~\sin ~{{\varphi }_{{{\text{см}}}}};~\,\,\,\,R_{i}^{'} = {{R}_{i}}\cos {{\theta }_{{j,i}}}~\,. \\ \end{gathered} $Условия на охлаждаемой внутренней поверхности стенки аппарата имеют вид:
где δ – шаг между гранулами, “выкладываемыми” роторным диспергатором по шестиграннику (четырехугольнику).Условие на охлаждаемой поверхности, не контактирующей с гранулами, имеет вид:
Система уравнений (1)–(27) замкнута и может быть решена численными методами.
МЕТОД ЧИСЛЕННОГО РЕШЕНИЯ ПРЕДЛОЖЕННОГО МАТЕМАТИЧЕСКОГО ОПИСАНИЯ
В работе использовали консервативный, неявный, абсолютно устойчивый, нелинейный, решаемый лишь итерационно на n-временном слое разностный аналог системы (1)–(27):
(28)
$\begin{gathered} \left( {c\rho } \right)_{{j,n}}^{S}\frac{{T_{{j,x,y,n}}^{{S + 1}} - T_{{j,x,y,n - 1}}^{{}}}}{{\Delta \tau }} = \\ = \lambda _{{j,n}}^{S}\Lambda _{{{{r}_{i}}}}^{2}T_{{j,x,y,n}}^{{S + 1}} + \lambda _{{j,n}}^{S}\Lambda _{{{{\varphi }_{i}}}}^{2}T_{{j,x,y,n}}^{S}, \\ \end{gathered} $(29)
$\begin{gathered} \left( {{{c}_{{{\text{СТ}}}}}{{\rho }_{{{\text{СТ}}}}}} \right)_{n}^{S}\frac{{T_{{m,l,n}}^{{S + 1}} - T_{{m,l,n - 1}}^{{}}}}{{\Delta \tau }} = \\ = \lambda _{{{\text{СТ}},n}}^{S}\Lambda _{z}^{2}T_{{m,l,n}}^{{S + 1}} + \lambda _{{{\text{СТ}},n}}^{S}\Lambda _{{{{x}_{i}}}}^{2}T_{{m,l,n}}^{S}, \\ \end{gathered} $Начальные условия:
(30)
$\begin{gathered} {{T}_{{j,x,y,0}}} = {{t}_{{{\text{C}}{\kern 1pt} {\text{'}}}}};{{r}_{{\zeta ,j,,i}}}\left( {x,y,0} \right) = {{R}_{i}} + 0; \\ {{r}_{{{\text{П}}j,i}}}\left( {x,y,0} \right) = 0;\,\,\,\,x = 0,1,2...X,\,\,\,y = 0,1,2...Y, \\ {{T}_{{m,l,0}}} = {{T}_{{m,l,{{n}_{{{\text{охл}}}}}}}} \vee {{t}_{{\text{C}}}};\,\,\,\,m = 0,1,2...M; \\ l = 0,1,2...L;\,\,\,\,Mx = {\delta \mathord{\left/ {\vphantom {\delta 2}} \right. \kern-0em} 2};\,\,\,\,Lz = {{\delta }_{{{\text{СТ}}}}};{{n}_{{{\text{охл}}}}}\tau = {{\tau }_{{{\text{охл}}}}}. \\ \end{gathered} $Условие начала фазового превращения имеет вид:
(31)
$\begin{gathered} \frac{2}{3}\pi \sum\limits_{n = {{n}_{1}}}^{{{n}_{2}}} {\sum\limits_{y = 0}^Y {\sum\limits_{x = {{x}_{1}}}^X {\left[ {\cos \Delta {{\theta }_{i}}y - \cos \Delta {{\theta }_{i}}\left( {y + 1} \right)} \right]} } } \times \\ \times \,\,\left( {r_{{i,x + 1,y,n}}^{3} - r_{{i,x,y,n}}^{3}} \right){{w}_{{3j}}}\left( {T_{{j,x,y,n}}^{{S + 1}}} \right)\Delta \tau = 1, \\ {{n}_{1}} = {{\tau \left( {{{t}_{{{\text{ф}}j}}},{{R}_{i}}} \right)} \mathord{\left/ {\vphantom {{\tau \left( {{{t}_{{{\text{ф}}j}}},{{R}_{i}}} \right)} {\Delta \tau }}} \right. \kern-0em} {\Delta \tau }};\,\,\,\,{{n}_{2}} = {{{{\tau }_{{{\text{инд}}j,i}}}} \mathord{\left/ {\vphantom {{{{\tau }_{{{\text{инд}}j,i}}}} {\Delta \tau }}} \right. \kern-0em} {\Delta \tau }}; \\ Y = {\pi \mathord{\left/ {\vphantom {\pi {\Delta {{\theta }_{i}};\,\,\,{{x}_{1}} = {{{{{{{r}_{i}}\left( {{{t}_{{{\text{ф}}j}}},{{\theta }_{i}},\tau } \right)} \mathord{\left/ {\vphantom {{{{r}_{i}}\left( {{{t}_{{{\text{ф}}j}}},{{\theta }_{i}},\tau } \right)} {\Delta r}}} \right. \kern-0em} {\Delta r}}}}_{i}}}}} \right. \kern-0em} {\Delta {{\theta }_{i}};\,\,\,{{x}_{1}} = {{{{{{{r}_{i}}\left( {{{t}_{{{\text{ф}}j}}},{{\theta }_{i}},\tau } \right)} \mathord{\left/ {\vphantom {{{{r}_{i}}\left( {{{t}_{{{\text{ф}}j}}},{{\theta }_{i}},\tau } \right)} {\Delta r}}} \right. \kern-0em} {\Delta r}}}}_{i}}}};\,\,\,\,X = {{{{R}_{i}}} \mathord{\left/ {\vphantom {{{{R}_{i}}} {\Delta {{r}_{i}}}}} \right. \kern-0em} {\Delta {{r}_{i}}}}. \\ \end{gathered} $Условие на границе сосредоточенной усадочной полости имеет вид:
(32)
$\begin{gathered} {{\Lambda }_{{{{n}_{i}}}}}T_{{j,x,y,n}}^{{S + 1}} = 0;\,\,\,\,V_{{{\text{П}}i,n}}^{{S + 1}} = V_{{{\text{к}}i,n}}^{{S + 1}}\frac{{{{\rho }_{{\text{к}}}}\left( {{{t}_{{{\text{фк}}}}}} \right) - {{\rho }_{{\text{ж}}}}\left( {{{t}_{{{\text{фк}}}}}} \right)}}{{{{\rho }_{{\text{ж}}}}\left( {{{t}_{{{\text{фк}}}}}} \right)}}; \\ \frac{{{{S}_{{{\text{П}}i,n}}}{{\sigma }_{n}}}}{{{{V}_{{{\text{П}}i,n}}}}} \to {\text{мин}}. \\ \end{gathered} $Объем распределенной усадочной пористости за счет образования трещин при термической усадке кристаллической фазы определяется зависимостями:
(33)
$\begin{gathered} V_{{{\text{ПР}}~j,i,n}}^{{S + 1}} = V_{{~j,i,n}}^{{S + 1}}\frac{{{{{\bar {\rho }}}_{{j,i,n}}} - {{\rho }_{j}}\left( {{{t}_{{\varphi j}}}} \right)}}{{{{{\bar {\rho }}}_{{j,i,n}}}}}; \\ ~{{{\bar {\rho }}}_{{j,i,n}}} = \frac{1}{{V_{{j,i,n}}^{{S + 1}}}}\frac{4}{3}\pi \mathop \sum \limits_{x = {{x}_{2}}}^{{{x}_{1}}} \left( {r_{{j,i,x,n}}^{3} - r_{{j + 1,i,x,n}}^{3}} \right)\rho _{{j,x,n}}^{{S + 1}};~ \\ \rho _{{j,x,n}}^{{S + 1}} = f\left( {T_{{j,x,n}}^{{S + 1}}} \right). \\ \end{gathered} $Объем распределенной усадочной пористости за счет появления трещин при полиморфных превращениях определяется зависимостями:
(34)
$V_{{{\text{ПР}}~j + 1,i,n}}^{'} = {{V}_{{j + 1,i,n}}}\frac{{{{\rho }_{{j + 1}}}\left( {{{t}_{{{\text{ф}}j}}}} \right) - {{\rho }_{j}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}{{{{\rho }_{j}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}.~$Условия на межфазных границах (кристаллизации и полиморфных превращений)
(35)
$\begin{gathered} T_{{j,x,y,n}}^{{s + 1}} = T_{{j + 1,x,y,n}}^{{s + 1}} < {{t}_{{{\text{ф}}j}}};~\,\,\,r_{{j,i,y,n}}^{{s + 2}} = \Delta {{r}_{i}}x; \\ ~{{\theta }_{i}}_{{{{\xi }_{i}},i,y,n}} = \Delta {{\theta }_{i}}y. \\ \end{gathered} $Граница раздела фаз смещается по нормали к границе со скоростью:
– разностный аналог проекции вектора скорости движения межфазной границы в пространстве x, y на единичный вектор нормали ͞ni к поверхности раздела фаз $j~--{{j}_{{\text{H}}}}$, (т.е. к локально изотермической поверхности).Градиент сеточной функции температуры в конкретной области границы j, j+1 можно представить в виде:
(37)
$\begin{gathered} {\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}} = ~{{\Lambda }_{{{{r}_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{{{\bar {i}}}_{r}} + \\ + \,\,\frac{1}{{{{r}_{{j,i,{{x}_{\xi }},y,n}}}}}{{\Lambda }_{{{{\varphi }_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{{{\bar {i}}}_{\varphi }} = \bar {\Lambda }~{\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}, \\ \end{gathered} $Этот вектор ортогонален к локально-изотермическому, линейно интерполированному участку межфазной границы. Из него получили единичный вектор нормали ${{\bar {n}}_{{j,i}}}$ к элементу кусочно-аппроксимированной межфазной границы, деля вектор градиента сеточной функции температуры у межфазной границы на его модуль. Единичный вектор нормали к аппроксимированной межфазной границе ${{\bar {n}}_{{j,i}}}$ в конечно-разностной форме можно определить:
(38)
${{\bar {n}}_{{j,i}}} = \frac{{{{{\bar {\Lambda }}}_{{{{r}_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{{{\bar {i}}}_{r}} + \frac{1}{{{{r}_{{j,i,{{x}_{\xi }},y,n}}}}}{{{\bar {\Lambda }}}_{{{{\varphi }_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{{{\bar {i}}}_{\theta }}}}{{{{{\left[ {{{{\left( {{{\Lambda }_{{{{r}_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right)}}^{2}} + {{{\left( {\frac{1}{{{{r}_{{j,i,{{x}_{\xi }},y,n}}}}}{{\Lambda }_{{{{\varphi }_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right)}}^{2}}} \right]}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}} = \frac{{\bar {\Lambda }~{\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}}}{{\left| {~\Lambda ~{\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right|}}.$Косинус угла между вектором скорости смещения межфазной границы:
(39)
${{\bar {\Lambda }}_{\tau }}{{\xi }_{{j,i,x,y,n}}} = {{\Lambda }_{\tau }}{{\xi }_{{j,i,x,y,n}}}\bar {\iota }{{{\kern 1pt} }_{\xi }}\,\, + {{\Lambda }_{{{{\varphi }_{j}}}}}{{\xi }_{{j,i,x,y,n}}}{{\Lambda }_{\tau }}{{\varphi }_{{j,i}}}{{i}_{{{{\varphi }_{i}}}}};$(40)
$\begin{gathered} \cos \left( {\bar {n},{{{\bar {\Lambda }}}_{\tau }}{{\xi }_{{j,i,x,y,n}}}} \right) = \frac{1}{{{{{\left( {{{{\left[ {{{\Lambda }_{{{{r}_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right]}}^{2}} + {{{\left[ {\frac{1}{{{{r}_{{j,i,{{x}_{\xi }},y,n}}}}}{{\Lambda }_{{{{\varphi }_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right]}}^{2}}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}{{{\left[ {{{{\left( {{{\Lambda }_{\tau }}{{\xi }_{{j,i,x,y,n}}}} \right)}}^{2}} + {{{\left( {{{\Lambda }_{{{{\varphi }_{j}}}}}{{\xi }_{{j,i,x,y,n}}}{{\Lambda }_{\tau }}{{\varphi }_{{j,i}}}{{i}_{\varphi }}} \right)}}^{2}}} \right]}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}} \times \\ \times \,\,{{\Lambda }_{\tau }}{{\xi }_{{j,i,x,y,n}}}{{\Lambda }_{{{{r}_{i}}}}}{{T}_{{j,i,{{x}_{\xi }},y,n}}} + {{\Lambda }_{{{{\varphi }_{j}}}}}{{\xi }_{{j,i,x,y,n}}}{{\Lambda }_{\tau }}{{\varphi }_{{j,i}}}{{i}_{\varphi }}\frac{1}{{{{r}_{{j,i,{{x}_{\xi }},y,n}}}}}{{\Lambda }_{\varphi }}{{T}_{{j,i,{{x}_{\xi }},y,n}}}. \\ \end{gathered} $Проекция ${{v}_{{\Lambda ,j}}}$ линейной скорости движения фазового фронта на единичный вектор к поверхности раздела фаз ${{\bar {n}}_{{j,i}}}$, рассчитывается с использованием разностного аналога:
(41)
$\begin{gathered} {{v}_{{\Lambda ,J}}} = {{\Lambda }_{\tau }}{{\xi }_{j}}_{{,i,x,y,n}}~{{{\bar {n}}}_{{j,i}}} = ~{{[{{({{\Lambda }_{\tau }}{{\xi }_{j}}_{{,i,x,y,n}})}^{2}} + {{({{\Lambda }_{{\varphi ,}}}_{i}{{\xi }_{j}}_{{,i,x,y,n}}{{\Lambda }_{\tau }}{{\varphi }_{{j,i}}})}^{2}}]}^{{1/2}}}\cos (~{{{\bar {n}}}_{{j,i}}}{{\Lambda }_{\tau }}{{\xi }_{j}}_{{,i,x,y,n}}) = \\ = \,\,{{\left( {{{{({{\Lambda }_{{r,i}}}{{T}_{j}}_{{,i,x,\xi ,y,n}})}}^{2}}\, + \,{{{\left( {\frac{1}{{{{r}_{{j,i,x\xi ,y,n}}}}}{{\Lambda }_{{\varphi ,}}}_{i}{{T}_{j}}_{{,i,x\xi ,y,n}}} \right)}}^{2}}} \right)}^{{ - 1/2}}}\left( {{{\Lambda }_{\tau }}{{\xi }_{{,j,i,y,n}}}~{{\Lambda }_{{r,i}}}{{T}_{{j,i,x\xi ,y,n}}}\, + \,{{\Lambda }_{{\varphi i,\xi ,j,i,x,y,n}}}{{\Lambda }_{\tau }}_{{,\varphi ,j,}}\frac{1}{{{{r}_{{j,i,x,\xi ,y,n}}}}}{{\Lambda }_{{I\varphi }}}{{T}_{{j,i,x\xi ,y,n}}}} \right)\, = \\ = ~\,\,\frac{1}{{\left[ {\Lambda {\text{\;grad}}~{{T}_{{j,i,{{x}_{i}},\xi ,y,n}}}} \right]}}({{{\bar {\Lambda }}}_{\tau }}{{\xi }_{{j,i,y,n}}}{{{\bar {\Lambda }}}_{{{\text{grad}}}}}{{T}_{{j,i,x}}}_{{\xi ,y,n}})~. \\ \end{gathered} $Вектор скорости движения фазовой границы по нормали к ее поверхности ${{\bar {n}}_{{j,i}}}$ ${{\bar {v}}_{{\Lambda ,j}}}$ равен величине ее проекции на нормаль ${{v}_{{\Lambda ,j}}}$ (скалярная величина) на единичный вектор нормали к межфазной поверхности ${{\bar {n}}_{{j,i}}}$, записанный через разностные аналоги производных по ri, θi:
(42)
$\begin{gathered} {{{\bar {v}}}_{{\Lambda ,J}}} = {{v}_{{\Lambda ,}}}_{j}{{{\bar {n}}}_{{j,i}}} = {{\Lambda }_{\tau }}{{\xi }_{{j,i,y,n}}}{{{\bar {n}}}_{{j,i}}}{{{\bar {n}}}_{{j,i}}} = \\ = \frac{1}{{\left[ {\Lambda {\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}} \right]}}({{{\bar {\Lambda }}}_{\tau }}{{\xi }_{{j,i,y,n}}}\bar {\Lambda }{\text{grad}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{\text{\;}}) \times \\ \times \,\,\frac{{\Lambda {\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}~}}{{\left[ {\Lambda {\text{grad}}~{{T}_{{j,i,x,y,n}}}} \right]}} = \\ = \frac{{\Lambda {\text{grad}}~{{T}_{{j,i,{{x}_{\xi }},y,n}}}}}{{{{{\left[ {\Lambda {\text{grad}}~{{T}_{{j,i,x,y,n}}}} \right]}}^{2}}}}({{{\bar {\Lambda }}}_{\tau }}{{\xi }_{{j,i,y,n}}}\bar {\Lambda }{\text{grad}}{{T}_{{j,i,{{x}_{\xi }},y,n}}}{\text{\;}}). \\ \end{gathered} $Тепловой баланс удельных потоков теплоты у межфазной границы (условие Стефана), записанный в разностной форме, имеет вид:
(43)
$\begin{gathered} \lambda _{{j + 1,n}}^{s}\bar {\Lambda }{\text{grad}}~T_{{j + 1,i,x,y,n}}^{{s + 1}} - ~\,\,\lambda _{{j,n}}^{s}\Lambda {\text{grad}}T_{{j,i,x,y,n}}^{{s + 1}} = \\ = r_{{{\text{ф}},j + 1}}^{s}\rho _{{j + 1}}^{s}\left( {\frac{{\Delta \xi _{{ijyn}}^{{s + 2}}}}{{\Delta \tau }}\bar {\Lambda }{\text{grad}}~T_{{j,i,x,y,n}}^{{s + 1}}} \right) \times \\ \times \,\,\frac{{\bar {\Lambda }{\text{grad}}~{{T}_{{~j,i,x,y,n}}}^{{s + 1}}}}{{{{{\left| {\Lambda {\text{grad}}~T_{{~j,i,x,y,n~~}}^{{s + 1}}} \right|}}^{2}}}} = v_{{\Lambda ,j}}^{{s + 2}}~r_{{{\text{ф}},~j + 1}}^{s}\rho _{{j + 1}}^{s}, \\ \end{gathered} $(44)
$\begin{gathered} v_{{\Lambda ,j}}^{{s + 2}} = ~\,\,{{v}_{{\Lambda ,j~}}}\,\,~ = \\ = \left\{ \begin{gathered} 0,\,\,\,\,~~r_{{\xi jiyn~}}^{{s + 2}}~\,\, = ~{{R}_{i}}~\,\, + \,\,0;\,\,\,\,~n~\Delta \tau ~\,\, < ~\,\,{{\tau }_{{{\text{инд}}}}} \hfill \\ ~{{f}_{1}}~~\left[ {{{t}_{{{\text{ф}},j}}}~ - T_{{j,x,y,n}}^{{s + 1}}} \right],\,\,\,\,~r_{{j,x,y,n}}^{{s + 2}}~\,\, < {{R}_{i}};\,\,\,\,~n~\Delta \tau ~\,\, \geqslant {{\tau }_{{{\text{инд}}}}}, \hfill \\ \end{gathered} \right. \\ \end{gathered} $Разностный аналог условия устойчивости межфазной границы имеет вид:
(45)
$c\lambda _{{j + 1,n~}}^{s}\Lambda {\text{grad}}T~_{{J + 1,x,y,n~}}^{{S + 1}}\,\, \geqslant \,\,~ - ~\lambda _{{jn~}}^{s}\Lambda {\text{grad}}T_{{j,i,x,y,n}}^{{s + 1}}.$Условие возникновения новых центров фазового превращения в метастабильной области перед межфазной границей:
(46)
$\begin{gathered} \frac{2}{3}\pi \sum\limits_{n = 0}^{{{n}_{з}}} {\sum\limits_{y = 0}^Y {\sum\limits_{x = {{x}_{1}}}^{{{x}_{{\zeta ji,y}}}} {\left[ {\cos \Delta {{\theta }_{i}}y - \cos \Delta {{\theta }_{i}}\left( {1 + y} \right)} \right]} } } \times \\ \times \,\,\left[ {r_{{i,x + 1,y,n}}^{3} - r_{{i,x,y,n}}^{3}} \right]{{w}_{{зj}}}\left( {T_{{j,x,y,n}}^{{S + 1}}} \right)\Delta \tau < 1, \\ {{n}_{{\text{з}}}} = \frac{{r_{{\zeta j,i,y,n}}^{{S + 2}} - r_{{i,y,n}}^{{S + 2}}\left( {{{t}_{{{\text{ф}}j}}}} \right)}}{{{{v}_{{\Lambda j}}}}};\,\,\,\,{{x}_{{\zeta j,i,y}}} = {{r_{{\zeta j,i,y,n}}^{{S + 2}}} \mathord{\left/ {\vphantom {{r_{{\zeta j,i,y,n}}^{{S + 2}}} {\Delta {{r}_{i}}}}} \right. \kern-0em} {\Delta {{r}_{i}}}}. \\ \end{gathered} $Условие на внешней границе капли:
(47)
$\begin{gathered} {{\alpha }_{{c{\kern 1pt} '}}}\left( {y,n} \right)\left[ {T_{{j,X,y,n}}^{{S + 1}} - {{t}_{{c{\kern 1pt} '}}}} \right] = - \lambda _{{j,n}}^{S}{{\Lambda }_{{{{r}_{i}}}}}T_{{j,X,y,n}}^{{S + 1}}; \\ y{\kern 1pt} ' < y \leqslant Y;\,\,\,\,y{\kern 1pt} ' = \frac{{{{\varphi }_{{{\text{см}}}}}}}{{\Delta {{\theta }_{i}}}}. \\ \end{gathered} $Условие на границе капля – стенка аппарата (гранулятора):
(48)
$\begin{gathered} - \lambda _{{{\text{СТ}},n}}^{S}{{\Lambda }_{Z}}T_{{m{\kern 1pt} ',0,n}}^{{S + 1}} = - \lambda _{{j,n}}^{S}{{\Lambda }_{{{{r}_{i}}}}}T_{{j,X{\kern 1pt} ',y,n}}^{{S + 1}}; \\ 0 \leqslant y \leqslant y{\kern 1pt} ';\,\,\,\,0 \leqslant m{\kern 1pt} ' \leqslant M; \\ X{\kern 1pt} ' = \frac{{{{R}_{i}}\cos {{\varphi }_{{{\text{см}}}}}}}{{{{r}_{i}}}};\,\,\,M{\kern 1pt} ' = \frac{{{{R}_{i}}\sin {{\varphi }_{{{\text{см}}}}}}}{x}. \\ \end{gathered} $Условие на охлаждаемой поверхности, не контактирующей с гранулами:
(49)
$\begin{gathered} {{\alpha }_{{c{\kern 1pt} ',i}}}\left( n \right)\left[ {T_{{m{\kern 1pt} '',0,n}}^{{S + 1}} - {{t}_{{c{\kern 1pt} '}}}} \right] = - \lambda _{{{\text{СТ}},n}}^{S}{{\Lambda }_{z}}T_{{m{\kern 1pt} '',0,n}}^{{S + 1}}; \\ M{\kern 1pt} ' < m{\kern 1pt} '' \leqslant M. \\ \end{gathered} $Условие на внутренней охлаждаемой поверхности, омываемой хладоагентом:
АЛГОРИТМ РЕШЕНИЯ ПОСТАВЛЕННОЙ ЗАДАЧИ
Учитывая, что математическое описание процесса имеет нелинейности I, II и III рода [17] и носит достаточно сложный характер, его решали численно [18–20]. Неявный (использовали в зависимости от удобства решения шести точечный (Кранка–Николсона) и трех точечный шаблоны), консервативный, устойчивый, нелинейный конечно-разностный аналог, решаемый лишь итерационно [17, 18]. С учетом этого предложена и использовалась локально- одномерная нелинейная итерационная на каждом шаге по времени схема решения с “подвижными узлами сетки”, отслеживающими положение межфазных фронтов после вычислений градиентов значений сеточной функции у подвижных границ, направлений, скорости перемещения последних, с совпадением с допустимой погрешностью материальных балансов, образующейся с учетом теплового баланса на границе новой фазы и “отсекаемой” переместившейся границей области старой фазы, с учетом всех видов усадки, контролем устойчивости межфазной границы и с оценкой вероятности зарождения центров новой фазы в метастабильном объеме старый фазы. Локально-одномерный сектор, отслеживающий изменение сеточной функции температуры и положения межфазных границ по радиусу (в пределах шага по углу Δθ) перемещался по угловой координате θ шагами Δθ с учетом симметрии задачи в пределах 0 ≤ θ ≤ π. Т.е. на каждой итерации в пределах шага по времени можно учитывать изменения всех теплофизических коэффициентов, граничных условий и геометрии границ (от поверхности капли, межфазных границ до границы и формы сосредоточенной усадочной полости). Этот алгоритм, в отличие от предложенных и апробированных на практике ранее [1, 9, 21], кроме радиального переноса теплоты, учитывает продольный по углу θ тепловой поток и обеспечивает уточненный расчет положений межфазных границ. Он особенно значим в условиях резкого изменения условий отвода теплоты по поверхности, претерпевающей фазовые превращения, капли расплава. При оценочных расчетах отмечена экономичность алгоритма. Например, погрешность численного расчета, оцениваемая по методу Рунге [18–20], не превышала 5% ± 1% с вероятностью 95% при разбивке пространственно-временного континуума на 10 шагов по координатам при 3–5 итерациях на каждом шаге по времени с использованием простых ускоряющих алгоритмов [20, 22]. Известным методом [18, 20] аппроксимировали граничные условия. На каждой итерации в пределах шагов по времени синхронизировали параметры пространственно-временных сеток в капле и охлаждающей стенке.
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
Адекватность вычислительного эксперимента физическому подтверждали сравнением с данными скоростного термического анализа (СТА). Каплю расплава выдавливали на рамку с 3-мя расположенными друг под другом термопарами Au–Au–Pt диаметром 0.05 мм, (т.о. один королек d = 0.01 мм термопары оказывался в близи свободной поверхности капли, другой – в центре, третий – около границы с охлаждаемой поверхностью (проточная ячейка из нержавеющей стали 50 × 51 мм δ = 10 мм).
Положение корольков термопар контролировали с помощью оптико-вычислительной системы, описанной в [1]. Работу проводили с расплавами ряда крупнотоннажных промышленных продуктов, которые условно разделили на вещества с высокой температурой кристаллизации (KNO3, NaOH, KOH, NH4NO3, шлак доменных печей) и вещества с умеренной температурой кристаллизации (карбамид, сера, ацетонанил, ряд других органических веществ). Большинство расплавов были оптически прозрачными, поэтому динамику процесса наблюдали с помощью оптико-вычислительной системы [1, 12]. Фиксировали температуру в 3-х вышеназванных точках (рис. 1). Обращали внимание на продолжительность индукционного периода τинд, который фиксировался оптико-расчетной системой [1] и термограммой рис. 1 (по времени появления вершины пика на термограмме вблизи охлаждаемой поверхности). Процесс всегда начинался у охлаждаемой поверхности (на стальной полированной поверхности он был систематически на 15–20% короче, чем рассчитанный по уравнению (5) для гомогенного зародышеобразования. Такая близость ко времени гомогенного зародышеобразования наблюдалось в случае, если охлаждающая поверхность была полированная промывалась, сушилась, обезжиривалась. При учете в расчетах радиального и продольного потоков теплоты совпадение расчета и эксперимента по предложенному математическому описанию было лучшим, чем при игнорировании продольного переноса теплоты, как полагали в ранее предложенных нами моделях [1, 9]. Хорошее согласование, термограмм и расчета, а также полученных экспериментально и рассчитанных по математической модели (1)–(27) размеров и формы усадочной поры, распределенной усадочной пористости дало основание полагать, что оценки устойчивости межфазной границы, механизма превращения (последовательный), определения скоростей зарождения и роста кристаллов (кинетических параметров превращения), теплообмена (параметров переноса при превращении) и описания в целом, адекватны ходу реального процесса.
ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ
Применительно к вышеназванным веществам (в качестве примера на рис. 1–5 приедены данные по KОН) были проведены расчеты процесса, когда отвод тепла с внешней стороны близок к его интенсивности на охлаждаемой поверхности (обдув снаружи продукта, испарительное охлаждение, аппарат без кожуха и т.д.). Процесс производительный, но не всегда допустим и экологически опасен (рис. 2). Процесс, когда с внешней стороны гранул отвод теплоты минимален (аппарат закрыт кожухом), с позиций экологии и безопасности предпочтителен, но производительность ниже (рис. 3). Заметна разница в динамике параметров кристаллизующейся капли (положений межфазных границ, полей температур, структуры усадочной полости). При расчете процесса и при оценке качества гранул – это надо учитывать (рис. 2, рис. 3).
Схема процесса гранулирования отверждением капель расплава или кристаллизации (чешуирования) на охлаждаемых поверхностях с использованием теплоты гранулирования с применением теплового насоса включает валковые, ленточные, дисковые и другие кристаллизаторы, запитываемые расплавом с последующим срезанием ножом чешуированного продукта. В случае гранулирования расплав “выкладывается” в виде капель ротационным монодисперсным гранулятором на охлаждаемую распыляемой форсункой водой поверхность. Горячие, со степенью кристалличности η ~ 0.7, предотвращающей их комкование [1, 9, 23] гранулы поступают в охладитель с псевдоожиженным слоем, охлаждающий их воздухом до температуры затаривания ~ 50°С. Возможна схема капсулирования не успевших закристаллизоваться, но уже не комкующихся η = 0.6–0.7 гранул, полимерными оболочками, образующимися из растворов полимеров (мономеров с инициаторами полимеризации) [24, 25] или их “догранулирования” с удалением воды (растворителя) из связующего за счет аккумулированной в гранулах теплоты – один из вариантов предлагаемой технологии “fattening” [26]. Особенно интересен этот процесс для веществ, кристаллизующихся при высоких температурах. В нашем случае KNO3, NaOH, KOH, NH4NO3, шлаки фосфорных и металлургических печей. Нами показано [1], что процесс можно использовать для гранулирования расплавов и более низкотемпературных веществ: карбамид, сера, NPK, NP, NS, NЦ, NK-удобрения, средства защиты растений, коагулянты, тугоплавкие парафины и др.
ЗАКЛЮЧЕНИЕ
Для экологически безопасных ресурсосберегающих процессов гранулирования кристаллизации капель расплава на охлаждаемых поверхностях с энергосбережением за счет использования теплоты гранулообразования с применением тепловых насосов рассмотрены принципиальные схемы с замкнутыми потоками хладагента.
Поскольку для технологического расчета и выбора необходимой аппаратуры нужен детальный анализ стадии гранулирования и охлаждения продукта, математические описания этого процесса дополнены учетом продольного переноса теплоты и двухмерного движения межфазных и усадочных границ в условиях переменных по поверхности гранул условий охлаждения. Предложен и апробирован алгоритм решения этой двухмерной задачи локально-одномерным численным методом конечных разностей с “ловлей” в узел сетки положений межфазных и усадочных границ и итерационным уточнением нелинейности I, II, III рода на каждой итерации “внутри” каждого временного шага.
Адекватность предложенного математического описания ходу реального процесса подтверждена сравнением расчетов с данными скоростного термического анализа, положением и объемом усадочной полости, экспериментальным определением наиболее вероятной продолжительности индукционного периода.
Работа проведена с крупнотоннажными веществами, кристаллизующимися при высоких температурах (около 200°С и выше) KNO3, NaOH, KOH, NH4NO3, NPK, NP, NS, NЦ, NK-удобрения, шлаки и веществами с температурами фазового перехода менее 150°С – карбамид, сера, ацетонанил и др.
Статья подготовлена в рамках выполнения базовой части государственного задания ФГАОУ ВО “Московский политехнический университет” (проект АААА-А20-120092190052-9).
ОБОЗНАЧЕНИЯ
i | пространственная координата, м |
j | номера фаз (название фаз) |
q | удельный поток тепла, Вт/м2 |
R, r | радиус гранулы, м |
S | номер итерации |
t | температура, град |
V | объем, м3 |
α | коэффициент теплоотдачи, Вт/(м2 K) |
∆ | приращение физической величины |
δ | толщина (линейный размер), м |
δст | толщина стенки кристаллизатора, м |
ε= Vn/V | коэффициент объемной пористости |
η | степень превращения |
λ | коэффициент теплопроводности, Вт/(м K) |
с | теплоемкость, Дж/кг |
ξ | координата подвижной межфазной границы, м |
ρ | плотность, кг/м3 |
σ | поверхностное натяжение, дж/м2 |
τ | время, с |
${{v}_{\Lambda }}$ | линейная скорость, м/с |
φсм | краевой угол смачивания, ° |
wз | скорость зародышеобразования, [м3 с]–1 |
Λ2 | разностный аналог оператора второй производной по соответствующим координатам |
$ \vee $ | или |
$ \wedge $ | и |
ИНДЕКСЫ
г | газовая фаза |
ж | жидкость (расплав) |
з | зарождение |
инд | индукционный |
к | кристаллическая фаза, конечное значение параметра |
кр | кристаллизация |
н | начальное значение параметра |
охл | охлаждение |
П | порошок, полость |
ПР | распределенная пористость |
ПЗ | переходная зона |
С | среда, охлаждающая (нагревающая) среда |
СТ | стенка |
ф | фаза |
n | нормаль |
Список литературы
Таран Ю.А. Разработка и анализ процессов гранулирования расплавов с использованием экологически безопасных энергосберегающих схем: дис. … канд. техн. наук. М.: МИТХТ, 2011.
Taran Y.A. Mathematical description of phase transformations and their formally analogous processes for engineering design // Theor. Found. Chem. Eng. 2020. V. 54 I. 3. P. 431–445 [Таран Ю.А. Математическое описание фазовых и формально аналогичных им превращений для инженерных расчетов // Теор. основы хим. технол. 2020. Т. 54. № 3. С. 309–324.]
Бесков В.С., Гетманцев В.С., Гетманцев С.В. Гранулирование на ленточном формователе коагулянта для водоочистки // Хим. промышленность сегодня. 2008. № 9. С. 50–53.
Процессы и аппараты химической технологии. Общий курс // Под ред. Айнштейна В.Г. Санкт-Петербург: Лань, 2019.
Пестова Г.С., Рычков Г.С., Сырнев Б.В. Метод грануляции расплава фторида бериллия. Ч.1. Диспергирование // Цветные металлы. 2006. № 5. С. 54–57.
Пестова Г.С., Рычков Г.С., Сырнев Б.В. и др. Метод грануляции расплава фторида бериллия. Ч. 2. Охлаждение // Цветные металлы. 2006. № 6. С. 62–67.
Пестова Г.С., Рычков Г.С., Сырнев Б.В. и др. Метод грануляции расплава фторида бериллия. Ч. 3. Установка грануляции фторида бериллия // Цветные металлы. 2007. № 6. С. 69–71.
Классен П.В., Гришаев И.Г. Шомин И.П. Гранулирование. М.: Химия, 1991.
Таран А.Л. Теория и практика процессов гранулирования расплавов и порошков: Дис. … докт. техн. наук. М.: МИТХТ, 2001.
Холин А.Ю., Таран А.Л. Способ кристаллизационной очистки веществ полунепрерывной зонной плавкой и аппарат для его осуществления (варианты). Пат. 2308555 РФ. 2005.
Таран Ю.А., Холин А.Ю. Таран А.Л. Кристаллизатор полунепрерывной зонной плавки. Пат. 24392213 РФ. 2012.
Taran Yu.A., Kholin A.Yu. Crystal nucleation and growth rates in single-component and binary melts under conditions of growth in a polycrystalline layer and in the case of single crystals // Theor. Found. Chem. Eng. 2011. V. 45. I. 3. P. 272–283. [Таран Ю.А., Холин А.Ю. Исследование кинетики скоростей зарождения и роста кристаллов в однокомпонентных и бинарных расплавах при их росте в поликристаллическом слое и в случае одиночных кристаллов. // Теор. основы хим. технол. 2011. Т. 45. № 2. С. 1–12.]
Kosheleva M.K., Meshalkin V.P., Dornyak O.R. Mathematical modeling of heat and mass transfer in the drying of granules used as a support for a nickel catalyst // Theor. Found. Chem. Eng. 2021. V. 55. P. 410–416. [Кошелева М.К., Мешалкин В.П., Дорняк О.Р. Математическое моделирование тепло- и массопереноса при сушке гранул носителя для никелевого катализатора. // Теор. основы хим. технол. 2020. Т. 54. № 2. С. 163–175.]
Флисюк О.М., Марцулевич Н.А., Шининов Т.Н. Гранулирование порошкообразных материалов в скоростном грануляторе // Мат. методы в технике и технологиях. 2016. № 2(84). С. 69–70.
Pokusaev B.G., Vyazmin A.V., Zakharov N.S., Khramtsov D.P., Nekrasov D.A. Unsteady Mass Transfer of Nutrients in Gels with Channels of Different Spatial Structures // Theor. Found. Chem. Eng. 2020. V. 54. I. 2. P. 277–288 [Покусаев Б.Г., Вязьмин А.В., Захаров Н.С., Храмцов Д.П., Некрасов Д.А. Нестационарный массоперенос питательных веществ в гелях с каналами различной пространственной структуры // Теор. основы хим. технол., 2020. Т. 54. № 2. С. 163–175.]
Polyanin A.D., Sorokin V.G., Vyazmin A.V. Reaction-Diffusion Models with Delay: Some Properties, Equations, Problems, and Solutions // Theor. Found. Chem. Eng. 2018. V. 52. P. 334–348. [Полянин А.Д., Сорокин В.Г., Вязьмин А.В. Нелинейные реакционно-диффузионные модели с запаздываением: некоторые свойства, уравнения, проблемы и решения // Теорет. основы хим. технологии. 2018. Т. 52. С. 334–348.]
Тихонов А.Н., Арсенин В.Я. Методы решений некорректных задач. М.: Наука. Главная редакция физико-матем. литературы, 1986.
Самарский А.А. Введение в теорию разностных схем. М.: Наука, 1971.
Самарский А.А. Теория разностных схем. М.: Наука. 1983.
Самарский А.А. Введение в численные методы. М.: Наука, 1982.
Таран А.Л., Долгалев Е.В., Таран Ю.А. Решение многомерных задач с фазовыми превращениями на локально-одномерном разностном аналоге // Математические методы в технике и технологиях ММТТ–20: сборник трудов ХХ Международной научной конференции. Изд-во ЯПУ, 2007. Т. 1. С. 150–154.
Химмельблау Д.М. Прикладное нелинейное программирование. М.: Мир, 1975.
Протодьяконов И.О., Марцулевич Н.А., Марков А.В. Явления переноса в химической технологии. Л.: Химия, 1981.
Таран А.Л., Олевский В.М., Шмелев С.Л. и др. Способ получения водоустойчивых гранул. Пат. 2038346 РФ. 1995.
Taran, Yu. A. Taran A.L. Revamping Equipment and Technology of Existing Granulating Plants to Improve the Quality of Granulated Products // Chemical and Petroleum Engineering. 2016. V. 51. I. 9–10. P. 581–5891 [Таран Ю.А., Таран А.Л. Аппаратурно-технологическое оформление реконструкции существующих производств гранулированных продуктов с целью повышения их качественных показателей // Хим. и нефтегаз. машиностроение. 2015. № 9. С. 3–8.]
Таrаn Yu.А., Моrоzоv R.V., Таrаn А.L. Calculation of process of finish-granulation by fattening technology with random granule-wise distribution of finish-granulating flow // Chemical and Petroleum Engineering. 2017. V. 53 I. 5–6 P. 355–363 [Таран Ю.А., Морозов Р.В., Таран А.Л. Расчет процесса догранулирования по технологии “fattening” со случайным распределением по гранулам догранулирующего потока // Хим. и нефтегаз. машиностроение. 2017. № 6. С. 3–8.]
Дополнительные материалы отсутствуют.
Инструменты
Теоретические основы химической технологии