STOCHASTIC MODEL OF SEPARATION PROCESS OF COMPLEX HETEROGENEOUS DISPERSION SYSTEMS IN CENTRIFUGES WITH RADIAL-AXIAL FLOW
Authors:
1.
Maritime State University named after G. I. Nevelskoy
(Doctor of Sciences in Technology, Professor; Head of the Department of Ship Internal Combustion Engines)
2.
Maritime State University named after admiral G. I. Nevelskoy
(Doctor of Technical Sciences, Professor; Professor of the Department of Marine Internal Combustion Engines)
Russian Federation
Russian Federation
3.
Maritime State University named after G. I. Nevelskoy
(Postgraduate Student of the Department of Ship Internal Combustion Engines)
Type:
Article
Pages:
from 34
to 45
Status:
Published
Received:
20.02.2019
Accepted:
25.02.2019
Published:
25.02.2019
Subject area:
UDK 621.43.013:629.3
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
Language:
Russian
Keywords:
stochastic model of centrifugation, Markov processes, separation of heterogeneous systems, hydrodynamics of centrifugal apparatuses, efficiency of fuel and oil cleaners, attrition rate
Abstract (English):
The article describes separation processes in apparatus with complex hydrodynamic conditions. The new approach to the theoretical study consists in taking into account the action of the centrifugal field and the field of random forces, as it was implemented earlier for the simplest structures, but also in a detailed examination of the radial and axial flows. The movement of particles is investigated in the cramped conditions of precipitation during their interaction with each other. The processes of separation in centrifugal devices of various types have been identified, subject to stochasticity, constrained movement of a set of real particles and the possibility of their dispersion and coagulation. The combined deterministic and stochastic effects on the dispersed phase in a centrifugal field are described by means of a differential equation, the solution of which makes it possible to calculate the fractional and total centrifugation efficiency with high accuracy. A semi-analytical method has been developed and implemented for one-and two-dimensional stochastic models for the approximate solution of non-stationary centrifugal boundary value problems. Compared to difference schemes Faedo-Galerkin method gives the solution of equation systems with a smaller number of unknowns defined in all values of the arguments. The use of the method is possible only in the case of homogeneous boundary conditions and is limited by the possibility of selecting a complete system of functions that would satisfy them. The proposed approach to identification of refining fuels and lubricants is preferable since it fully reflects the essence of the processes occurring in centrifuges with a radial-axial flow. Generalization of stochastic effects on centrifugation, taking into account the hydrodynamic features of the flow of shared media in a centrifugal field, allows calculating all indicators of cleaning fuels and oils from mechanical impurities without model experiments, predicting the total and fractional efficiency of fuel oil filters and oil filters of the internal combustion engines. The developed identification method of separating complex dispersed systems is aimed at the synthesis of systems and devices for oil purification and fuel oil preparation of increased efficiency with predetermined parameters of purification quality in operating conditions.
The article describes separation processes in apparatus with complex hydrodynamic conditions. The new approach to the theoretical study consists in taking into account the action of the centrifugal field and the field of random forces, as it was implemented earlier for the simplest structures, but also in a detailed examination of the radial and axial flows. The movement of particles is investigated in the cramped conditions of precipitation during their interaction with each other. The processes of separation in centrifugal devices of various types have been identified, subject to stochasticity, constrained movement of a set of real particles and the possibility of their dispersion and coagulation. The combined deterministic and stochastic effects on the dispersed phase in a centrifugal field are described by means of a differential equation, the solution of which makes it possible to calculate the fractional and total centrifugation efficiency with high accuracy. A semi-analytical method has been developed and implemented for one-and two-dimensional stochastic models for the approximate solution of non-stationary centrifugal boundary value problems. Compared to difference schemes Faedo-Galerkin method gives the solution of equation systems with a smaller number of unknowns defined in all values of the arguments. The use of the method is possible only in the case of homogeneous boundary conditions and is limited by the possibility of selecting a complete system of functions that would satisfy them. The proposed approach to identification of refining fuels and lubricants is preferable since it fully reflects the essence of the processes occurring in centrifuges with a radial-axial flow. Generalization of stochastic effects on centrifugation, taking into account the hydrodynamic features of the flow of shared media in a centrifugal field, allows calculating all indicators of cleaning fuels and oils from mechanical impurities without model experiments, predicting the total and fractional efficiency of fuel oil filters and oil filters of the internal combustion engines. The developed identification method of separating complex dispersed systems is aimed at the synthesis of systems and devices for oil purification and fuel oil preparation of increased efficiency with predetermined parameters of purification quality in operating conditions.
Keywords:
stochastic model of centrifugation, Markov processes, separation of heterogeneous systems, hydrodynamics of centrifugal apparatuses, efficiency of fuel and oil cleaners, attrition rate
stochastic model of centrifugation, Markov processes, separation of heterogeneous systems, hydrodynamics of centrifugal apparatuses, efficiency of fuel and oil cleaners, attrition rate
Введение Процесс разделения в центробежных очистителях (ЦО) с радиально-осевым сливом усложнён наличием радиальной скорости потока. Кроме того, в этих аппаратах из-за путевого отбора фугата осевая скорость потока зависит от координаты в направлении оси ротора. Ряд центрифуг имеет очень сложную гидродинамическую обстановку, обусловленную особенностями подвода и отвода жидкости в ротор [1, 2], конструкцией питающих каналов, наличием свободной поверхности центрифугируемой жидкости [3]. Гидродинамика таких объектов изучена недостаточно полно. Учёт её специфики в расчёте процесса разделения, особенно в условиях стохастичности, должного развития не получил. С позиций оценки эффективности центрифугирования ЦО указанного типа не исследованы. Весьма трудна для них идентификация осаждения при изменении агрегатного состояния дисперсной фазы (ДФ) [4]. Новый подход к идентификации процесса центрифугирования горюче-смазочных материалов (ГСМ) заключается в комплексном учёте полидисперсности и многофазности загрязнителя, взаимодействия частиц ДФ, стохастичности разделения и сложной гидродинамики потоков в аппаратах очистки. Разработка теории разделения сложных дисперсионных сред (ДС) в центробежном поле с исследованием специфики твёрдой фазы и гидродинамической обстановки в ЦО главным образом нацелена на адекватную оценку эффективности центрифугирования нефтепродуктов в системах топливоподготовки и маслоочистки дизелей при очистке их от вредных нерастворимых примесей. Особенностью подхода является использование дифференциального уравнения движения частицы в радиальном и осевом направлениях в поле действия центробежных сил инерции [5]. Уравнение, описывающее детерминированное воздействие на ДФ, дополняется членом, учитывающим стохастичность процесса разделения. Постановка и решение многомерных задач центрифугирования в аппаратах со сложной гидродинамикой Случайная составляющая скорости может быть вызвана турбулентными пульсациями, стесненностью движения ДФ, ограниченностью пространства, приводящими к флуктуации. Стохастичность обусловлена также объединением частиц при столкновениях и разрушением вторичных образований вследствие действия гидродинамических факторов [4]. Ранее [6] было доказано, что случайное воздействие на ДС в центробежном поле является дельта-коррелированной функцией времени с неравным средним значением и заданной интенсивности. Центрифугирование можно считать простым марковским процессом и выражать его эффективность через многомерную плотность вероятности. Перечисленные представления позволяют учитывать при оценке качества очистки ГСМ совместные действия на осаждение частиц в центробежном поле детерминированных и случайных факторов. Особенностью очистителей с комбинированным радиально-осевым сливом (рис. 1), которые нашли широкое применение на судах, является наличие напорного осевого потока Qz и безнапорного радиального слива Qr (первый, как правило, представляет собой фугат, поступающий потребителю; второй - поток, идущий на гидрореактивный привод ротора при истечении через сопла). Возможны и другие схемы реализации Qz и Qr в системах топливоподготовки и маслоочистки на судах [7]. Переносная скорость движения частиц в роторе центрифуги по z и r принята равной скорости жидкости в этом же направлении. Скорость потока центрифугируемой жидкости uz и ur, если считать её идеальной, можно определить решением уравнения Лапласа для потенциала скорости y [8]: (1) Приняв осевую скорость uz на входе в центрифугу равной и радиальную ur0 на радиусе получим решение уравнения (1) в виде (2) Рис. 1. Расчётная схема центрифуги с комбинированным сливом: ω - угловая скорость вращения ротора, м/с; u(r, z) - осевая и радиальная скорости движения потока в роторе, м/с; u(r, z) - скорость движения частицы в координатах r и z, м/с; ε - коэффициент выноса частиц из ротора; br - интенсивность случайного воздействия на частицу в направлении координаты r, Дж/(кг·с); F0(d) - дифференциальная функция распределения массы частиц по диаметру, мкм-1; Rц - радиус ротора, м; z, r - координаты положения частиц в роторе, м; Hц - высота ротора, м Вследствие путевого радиального отбора жидкости по высоте центрифуги выражение для осевой и радиальной скорости потоков внутри ротора, согласно (1), имеет вид [6]: Если пренебречь изменением осевой скорости потока по высоте ротора и принять её равной то уравнение, описывающее воздействия внешних сил на радиальное перемещение по радиусу r частицы диаметром d в центробежном поле с частотой вращения ω, преобразуется следующим образом: (3) где ξ(τ) - совокупность случайных воздействий на ДС; μм - вязкость центрифугируемой жидкости, Па·с; ρэф = ρd - ρ - эффективная плотность частиц, кг/м2; kψ, kc - коэфициенты, учитывающие отклонения размера частицы от стохастивного диаметра (учитывает влияние формы частиц и их концентрацию) [5]; τ - время процесса, с. Для рассматриваемой конструкции очистителя концентрация частиц W(r, d, t) может быть определена решением дифференциального уравнения (4) с соответствующими выражению (3) коэффициентами Граничные условия, указывающие на отсутствие перемещения частиц вдоль радиуса при достижении ими стенки аппарата (r = Rц), a также зоны противотока (r = r0) и определяющие унос частиц из центрифугируемой жидкости и с радиальным сливом соответственно, могут быть записаны в виде . За начальное при t = 0 принимали условие: характеризующее равномерное распределение частиц одного диаметра по площади поперечного сечения ротора. После введения безразмерных переменных дифференциальное уравнение (4), граничные и начальные условия имеют вид: (5) где . Решение системы (5) искали по методу Фаэдо - Галёркина [9, 10] по следующему алгоритму: - находятся собственные функции дифференциального оператора, соответствующего переменной d; - для заданных переменных определяются собственные функции дифференциального оператора и по ним составляется набор задач Штурма - Лиувилля; - после решения трансцендентного уравнения находятся собственные числа и по каждому из них записываются соответствующие функции; - решается система обыкновенных дифференциальных уравнений для определения Сd(t) и фиксируется приближенное решение для вычисления плотности вероятности; - исходя из граничных условий составляются интегральные выражения для расчёта эффективности центрифугирования. (6) где - ортонормированные собственные функции (2) задачи: (7) Коэффициенты Сk(t) определялись из соотношения: (8) Решение задачи (7) представлено в виде Подставляя найденные в (6) и (7), получили систему линейных обыкновенных дифференциальных уравнений первого порядка с постоянными коэффициентами, которую использовали для определения : (9) Интегрируя входящее в (9) выражение по частям, получили: Для вычисления Iik применили одну из квадратурных формул [11, 12]. Задача (9) в матричной форме записывается в виде где Окончательно приближенное решение задачи (6) приняло вид: (10) Относительная концентрация частиц диаметра d в поперечном сечении ротора центрифуги в момент составила: (11) По физическому смыслу эквивалентно коэффициенту пропуска edz в направлении оси z, т. е. характеризует относительное количество частиц, уходящих с потоком Qz. Доля их, уносимая потоком Qr, находится, согласно (5), как поток вероятности через внутреннюю обечайку ротора: (12) Аналогично рассчитываем поток частиц через наружную обечайку ротора, определяющий полный фракционный коэффициент отсева центрифуги: (13) Переход к локальным задерживающим характеристикам центрифуги можно осуществить через безразмерные потоки очищаемости жидкости (14) (15) где jdz и jdr - фракционные коэффициенты отсева для потоков Qz и Qr - отображают качество очистки жидкости в напорном и свободном сливах. Полный и локальные коэффициенты отсева связаны посредством относительных показателей эффективности очистки и : (16) Если состав загрязнителя на входе в ЦО задан дифференциальной функцией распределения F(d), то полноту отсева jц можно получить на основе (16) интегрированием выражения [6] (17) где F(d) - задаваемая дифференциальная функция распределения размера частиц по d. Интенсивность случайных воздействий b0r, Дж/(кг·с), является квазидиффузионным параметром и отвечает за подвижность частиц в роторе центрифуги. По своему физическому смыслу она может быть интерпретирована как количество энергии, передаваемое в единицу времени ДФ массой в 1 кг в результате взаимного столкновения частиц при стесненном их движении. Степень подвижности частиц в различных аппаратах неодинакова и зависит от конструктивных особенностей очистителя, его расходной характеристики и общей концентрации твёрдой фазы в суспензии С на входе в очиститель [13, 14]. Стохастичность процесса центрифугирования предлагается учитывать по выражению1, найденному методом множественной корреляции: (18) Анализ (18) доказывает рациональность конструкций с радиально-осевым сливом для управления фракционной задерживающей характеристикой очистителя. Увеличение радиального слива способствует уменьшению подвижности и флуктуации скорости частиц в радиальном направлении, что выгодно для полнопоточных центрифуг, когда необходимо понизить полноту отсева ГСМ от общих загрязнений при полном удалении опасных частиц механических примесей. Продолжительность работы ЦО между чистками ротора в таком случае повышается. Расчётно-экспериментальная оценка эффективности центрифуги с радиально-осе-вым потоком Возможности разработанной модели (5) проверены на очистителе МЦН-5П (КС) с комбинированным сливом. Из графических зависимостей (рис. 22) видно, что фракционные коэффициенты пропускания и для мелких частиц и коэффициенты отсева , для d ³ 30 мкм определяются в основном соотношением потоков Qz и Qr. Рис. 2. Общая (а) и локальная (б) фракционная эффективность центрифуги с комбинированным сливом: Qz = 0,67·10-3м3/с; Qr = 0,33·10-3м3/с; Vp = 1 560 см3; w = 760 рад/с; Hц = 152 мм; r0 = 18 мм; Rц = 60 мм; b0 = 2,5·10-13 Дж/(кг·с); m = 0,03 Па·с; rd = 2 500 кг/м3; r = 870 кг/м3; - расчётные данные; - - - экспериментальные результаты Адекватность расчётной модели (5) и соответствие её экспериментальным данным проверены по локальным фракционным коэффициентам отсева (рис. 3). Рис. 3. Влияние безразмерных параметров (а) и (б) на унос и отфуговывание ДФ при очистке моторного масла Сходимость результатов расчёта и эксперимента высока. Теоретические зависимости полностью располагались в доверительных границах экспериментальных значений jdz и jdr. Коэффициенты отсева, полученные в опытах, отличались от расчётных максимально на 12 %. Среднее квадратическое отклонение экспериментальных точек от теоретической и аппроксимирующей кривых примерно одинаково и составляло 4-7 %. Погрешность эксперимента при определении фракционных коэффициентов отсева находилась в пределах 3-6 %. Для тех же условий центрифугирования выполнено решение при разном времени = adrt пребывания частиц в центробежном поле. Расчётом и экспериментально определено влияние на jdц не только , но и безразмерного параметра (см. рис. 3, б), представляющего отношение детерминированного воздействия к стохастическому. Сходимость данных расчёта и эксперимента для функций edi, jdi() и jdц() - хорошая. Сопоставление расчётных и экспериментальных зависимостей по критерию Вилькоксона [14] подтверждает, что они принадлежат к одной генеральной совокупности. Теоретические зависимости почти полностью располагаются в поле рассеяния экспериментальных данных, полученных с доверительной вероятностью 95 % по результатам трёх параллельных проб. При понижении вследствие роста интенсивности случайных воздействий на ДФ эффективность отсева частиц увеличивается. Если > 8, стохастичностью процесса для частиц d > 10 мкм можно пренебречь. У большинства конструкций ЦО при очистке ГСМ от крупнодисперсных загрязнителей находится в диапазоне 2-25. Для некоторых ДС при разделении их в центрифуге с цилиндрическим ротором радиальные составляющие силы инерции и силы Кориолиса имеют один и тот же порядок, и пренебречь одной из них (что сделано во многих работах) недопустимо. Достаточно высокую точность расчёта отфуговывания тяжёлых частиц ДФ (ρэф ˃ 3 000 кг/см3) можно получить, учитывая в одной из составляющих ard в формуле (4) действие сил Кориолиса [5]. Адекватность детерминированной модели (ДМ) при определении полноты отсева видна из сопоставления результатов расчёта и эксперимента (табл.). Эффективность очистителя МЦН-5П (КС)*по полноте отсева Загрязнитель Полнота отсева, % Погрешность расчёта, % Расчётная Экспериментальная Кварцевый с Sуд = 1,05 м2/г: dm = 5,5 мкм; vd = 0,7 16,2** 18,4 16,8 ± 1,5 16,8 ± 1,5 3,6 9,5 Кварцевый с Sуд = 0,56 м2/г: dm = 11,7 мкм; vd = 0,85 46,3 44,8 47,2 ± 3,1 47,2 ± 3,1 1,9 5,1 Продукты регенерации саморегенерирующегося фильтра с параметрами: bd = 8,6 мкм; pd = 1,84 25,8 29,2 27,1 ± 1,9 27,1 ± 1,9 4,8 7,7 bd = 2,3 мкм; pd = 2,08 2,5 3,1 2,4 ± 0,2 2,4 ± 0,2 4,2 29,2 * Параметры ЦO приведены на рис. 2. ** Расчёт выполнен по стохастической модели (СМ). Для крупнодисперсных загрязнителей сходимость теоретических и опытных данных по полноте отсева - хорошая. Для мелкодисперсных нерастворимых продуктов погрешность расчёта превышает 15 %, что требует использования СМ. Формулы расчёта эффективности центрифугирования по ДМ [3] рекомендованы для условий центрифугирования: ³ 5. При < 5 удовлетворительные результаты расчёта jdц могут быть получены только по СМ. Баланс ДФ при прохождении ДС через ротор зависит от распределения потоков (см. рис. 3, а). Увеличение доли Qr способствует уменьшению коэффициента уноса частиц в главную магистраль дизеля. При этом общая эффективность очистки jdц и jц как моно-, так и полидисперсной фазы загрязнения практически не понижается. Указанная особенность важна для полнопоточных очистителей. Через поток Qr можно управлять балансом загрязнений в роторе и значительно повысить эффективность центрифугирования поступающего к потребителю масла [15]. Разработанная СМ показывает способы реализации в условиях стохастичности возможностей и преимуществ ЦО с радиально-осевыми потоками [7]. Расчётные общие характеристики по коэффициенту и полноте отсева очистителя с оптимизированными потоками Qz и Qr хорошо согласуются с результатами эксперимента (см. рис. 3, б; табл.). Расхождение между теоретическими и экспериментальными данными по jц значительно больше, чем по jdц. Теоретические зависимости располагаются в границах доверительного интервала почти всех экспериментальных точек. Разработанные теоретические основы центрифугирования с учётом стохастичности в 1,5-4 раза повысили точность расчёта эффективности ЦО. Расчётные методики универсальны и могут быть использованы для создания аппаратов различных конструкций и назначений [15]. Предложенная теория применима в системе автоматизированного проектирования высокоэффективных очистителей, а также при разработке новых научно-технических решений и направлений инженерного оформления аппаратов и систем топливо- и маслоочистки судовой дизельной энергетической установки. Выводы 1. Разработаны с использованием представлений и аппарата случайных марковских процессов принципиально новые СМ центробежной очистки многофазных гетерогенных систем с полидисперсным загрязнителем, характерные для работающих в двигателях внутреннего сгорания моторного масла и топлив. Модели позволяют синтезировать очистительные комплексы высокой эффективности для условий случайных интенсивных воздействий на ДФ рабочих сред судовой дизельной энергетической установки в эксплуатации. 2. Обобщение опыта центрифугирования и близких ему разделительных процессов подтвердило: - случайные воздействия на ДС в центробежном поле являются дельта-коррелированной функцией времени с нулевым средним значением и заданной интенсивностью; - центрифугирование можно считать простым марковским процессом и выражать его эффективность через многомерную плотность вероятности. Перечисленные представления позволяют учитывать при оценке качества очистки совместные действия на осаждение частиц в центробежном поле детерминированных и случайных факторов. 3. Приведена общая схема представления и реализации многомерных задач центрифугирования при любом числе действующих на ДФ детерминированных и случайных факторов. Проиллюстрированы способы формирования начальных и граничных условий для ЦО разного конструктивного исполнения. Разработаны расчётные зависимости для локальной и общей оценки эффективности центрифугирования по фракционному коэффициенту и полноте отсева (уноса). Идентификация процесса разделения сложных многофазных систем в центробежном поле осуществлена по многоступенчатой схеме: постановка и решение многомерных стохастических задач центрифугирования, реализация упрощенных задач в детерминированной постановке с коррекцией на стохастичность и без неё. Получен критерий, связывающий детерминированное и стохастическое воздействия на частицы, по которому определяется правомерность использования ДМ очистки. Стохастическая модель обеспечивает высокую точность расчёта эффективности центрифугирования как с разрушением, так и без разрушения ДФ. 4. Адекватность разработанных СМ доказана по критерию Вилькоксона. Согласно данному критерию, все теоретические кривые jdц(d) и edb(d) согласуются с экспериментальными зависимостями фракционного коэффициента отсева (уноса) от стоксовского диаметра частиц. Они принадлежат к одной генеральной совокупности. Разница между расчётными и экспериментальными jdц(d) кривыми статистически незначима. Теоретические зависимости, найденные полуаналитическим методом Фаэдо - Галёркина на основе усложненных одно- и многомерных моделей, располагаются в основном в доверительных границах экспериментальных значений jdц.
1. Shepel'skiy Yu. L. Razvitie konstrukciy masloochistiteley dlya sudovyh dizeley // Dvigatelestroenie. 1985. № 7. S. 20-23.
2. Kicha G. P. Ochistka masla v dvigatelyah vnutrennego sgoraniya // Himiya i tehnologiya topliv i masel. 1985. № 2. S. 28-30.
3. Sokolov V. I. Sovremennye promyshlennye centrifugi. M.: Mashinostroenie, 1967. 523 s.
4. Kicha G. P., Lapin A. M. Identifikaciya processa centrifugirovaniya sistem s neustoychivoy dispersnoy fazoy // Povyshenie urovnya tehnicheskoy ekspluatacii dizeley rechnogo flota: sb. nauch. tr. Novosibirsk: NIIVT, 1988. S. 140-153.
5. Kicha G. P. Novye stohasticheskie modeli processa ochistki goryuche-smazochnyh materialov v DVS // Dvigatelestroenie. 1989. № 11. S. 18-23.
6. Kicha G. P., Semenyuk L. A., Taraschan N. N. Optimizaciya parametrov kombinirovannogo maslyanogo fil'tra, funkcioniruyuschego v sostave kompleksa «dizel'-toplivo-maslo» // Transportnoe delo Rossii. 2017. № 4. S. 96-102.
7. Kicha G. P., Perminov B. N., Nadezhkin A. V. Resursosberegayuschee masloispol'zovanie v sudovyh dizelyah: monogr. Vladivostok: Izd-vo MGU im. adm. G. I. Nevel'skogo, 2011. 372 s.
8. Loycyanskiy L. G. Mehanika zhidkostey i gazov. M.: Vyssh. shk., 1982. 685 s.
9. Fletcher K. Chislennye metody na osnove metoda Galerkina. M.: Mir, 1988. 352 s.
10. Kicha G. P., Zagorodnikov Yu. I., Osipov O. V. Stohasticheskaya model' processa centrifugirovaniya motornogo masla v SEU // Povyshenie urovnya tehnicheskoy ekspluatacii dizeley rechnogo flota: sb. nauch. tr. Novosibirsk: NIIVT, 1988. S. 126-139.
11. Mihlin S. G. Variacionnye metody v matematicheskoy fizike. M.: Nauka, 1970. 512 s.
12. Mihlin S. G. Lineynye uravneniya v chastnyh proizvodnyh. M.: Vyssh. shk., 1977. 431 s.
13. Vorob'ev B. N., Nadezhkin A. V., Kicha G. P. Stohasticheskoe modelirovanie razdeleniya slozhnyh geterogennyh sistem sudovyh ustroystv na osnove predstavleniy i apparata sluchaynyh markovskih processov // Mor. intellektual. tehnologii. 2017. № 3. T. 2. S. 112-121.
14. Kicha G. P., Nadezhkin A. V., Semenyuk L. A. Novye stohasticheskie modeli ochistki topliv i masel sudovymi centrobezhnymi apparatami so slozhnoy gidrodinamicheskoy obstanovkoy // Mor. intellektual. tehnologii. 2018. № 4. T. 5. S. 77-90.
15. Kicha G. P. Ekspluatacionnaya effektivnost' novyh masloochistitel'nyh kompleksov v forsirovannyh dizelyah // Dvigatelestroenie. 1987. № 6. S. 25-29.
