автореферат диссертации по технологии продовольственных продуктов, 05.18.12, диссертация на тему:Процессы переноса массы и тепла системы "жидкость-жидкость" в аппаратах с объемом, структурированном мембранными перегородками

На правах рукописи

СХАЛЯХОВ Анзаур Адамович

ПРОЦЕССЫ ПЕРЕНОСА МАССЫ И ТЕПЛА СИСТЕМЫ «ЖИДКОСТЬ-ЖИДКОСТЬ» В АППАРАТАХ С ОБЪЕМОМ, СТРУКТУРИРОВАННОМ МЕМБРАННЫМИ ПЕРЕГОРОДКАМИ

Специальность 05.18.12 - Процессы и аппараты пищевых производств

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук

2 3 СЕН 2010

Краснодар - 2010

004608451

Работа выполнена в ГОУ ВПО «Майкопский государственный технологический университет»

Научный консультант: доктор технических наук, профессор

Блягоз Хазрет Рамазанович

Официальные оппоненты: доктор технических наук, профессор

Шаззо Аслан Юсуфович

доктор технических наук, профессор Антипов Сергей Тихонович

доктор технических наук, профессор Дворецкий Станислав Иванович

Ведущая организация: Северо-Кавказский филиал

Всероссийского научно-исследовательского института жиров Россельхозакадемии

Защита состоится "5" октября 2010 года в 13е0 часов на заседании диссертационного совета Д 212.100.03 Кубанского государственного технологического университета по адресу: 350072, г. Краснодар, ул. Московская, 2

С диссертацией можно ознакомиться в библиотеке Кубанского государственного технологического университета

Автореферат диссертации разослан «3» сентября 2010 года

Ученый секретарь

диссертационного совета, /У

канд. техн. наук, доцент щ/ Филенкова М.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Диссертационная работа выполнена в Майкопском государственном технологическом университете в рамках подпрограммы: «Раз- -работка инновационных технологий и техники переработки сельско-

хозяйственного сырья и производства продовольствия на основе мембранных и экстракционных процессов» Республиканской целевой программы «Развитие сельского хозяйства и регулирование рынков сельскохозяйственной продукции, сырья и продовольствия», утвержденной постановлением Государственного Совета-Хасэ Республики Адыгея от 30 января 2008 года №689-ГС.

Цель и задачи работы. Целью работы является научное обоснование, разработка процессов и создание оборудования для переноса массы и тепла системы «жидкость-жидкость» в аппаратах с объемом, структурированном мембранными перегородками, создание на этой основе технологии комплексной переработки растительных масел.

Для достижения поставленной цели решались следующие основные задачи:

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

- поставить и решить сопряженные задачи переноса массы и тепла диффузией и теплопроводностью через трубчатую мембрану;

- провести моделирование и обоснование обобщенных кинетических зависимостей переноса массы и тепла в трубчатых мембранах;

- получить математическое описание и промоделировать перенос потоков в мембранах под действием трансмембранного перепада давления;

- разработать методику и определить проницаемости половоло-конных и трубчатых мембран;

- исследовать свойства компонентов, экстрагируемых из растительных масел и провести математическое моделирование процесса мембранной экстракции смеси триацилглицеринов и жирных кислот;

- исследовать процесс массопереноса в мембранном контакторе;

- провести математическое моделирование кинетики и свойств реакционной смеси при производстве биодизельного топлива из масел;

- построить математическую модель мембранного реактора для производства биотоплива и на этой основе провести оптимизацию процесса;

- обосновать конструкцию и математическую модель конденсатора с непористыми полипропиленовыми половолоконными мембранами для парогазовых смесей экстракционного производства;

- провести исследования теплообмена с непористыми полипропиленовыми половолоконными мембранами;

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

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

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

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

Изложенная концепция открывает научное направление - создание новых процессов и аппаратов тепло- и, массообмена в системе «жидкость-жидкость» со структурированием объема мембранными перегородками.

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

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

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

Развита теория переноса потоков в мембранных контакторах под действием трансмембранного перепада давления.

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

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

Разработана методика идентификации вязкости компонентов реакционной смеси в процессе переэтерификации растительного масла с метанолом на основе применения математической модели кинетики реакции и модели вязкости смеси, при этом температурная зависимость вязкости компонентов определена по усовершенствованному методу Льюиса и Сквайрса.

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

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

Новизна технических решений подтверждена 3 патентами и свидетельствами РФ на полезные модели и 2 решениями о выдаче патентов РФ на изобретения.

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

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

Разработаны методики расчетов мембранного контактора для экстракции свободных жирных кислот из растительных масел дву-

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

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

Разработанные технологическая схема переработки растительных масел и конструкции мембранных аппаратов приняты к внедрению в условиях филиала «Лабинский МЭЗ» ООО «МЭЗ Юг Руси», а также переданы Министерству сельского хозяйства Республики Адыгея для проектирования завода по переработке семян подсолнечника. Ожидаемый экономический эффект составит более 17 млн. рублей при переработке 15000 тонн нерафинированного подсолнечного масла в год.

Апробация работы. Результаты работы представлены на: СевероКавказской региональной научной конференции аспирантов и молодых ученых «Перспекггива-98», г. Нальчик, 1998 г.; П Международной научной конференции «Продовольственный рынок и проблемы здорового питания», г. Орел, 1999 г.; Международной научно-практической конференции «Индустрия продуктов здорового питания — третье тысячелетие (человек, наука, технология, экономика)», г. Москва, 1999 г.; Международной научно-практической конференции «Молодые ученые -пищевым и перерабатывающим отраслям АПК (технологические аспекты производства)», г. Москва, 2000 г.; Международной научной конференции «Прогрессивные пищевые технологии - третьему тысячелетию = Advanced food technologies - to the third millennium», г. Краснодар, 2000 г.; Международной научно-практической конференции «Продовольственная индустрия юга России. Экологически безопасные энергосберегающие технологии хранения и переработки сырья растительного и животного происхождения», г. Краснодар, 2000 г.; Международной научно-технической конференции, посвященной 70-летию Санкт-Петербургского государственного университета низкотемпературных и пищевых технологий «Низкотемпературные и пищевые технологии в XXI веке», г. С.-Петербург, 2001 г.; Международной научно-практической конференции «Научные основы процессов, аппаратов и машин пищевых производств», г. Краснодар, 2002 г.; V Всероссийской научно-практической конференции студентов, аспирантов, докторантов и молодых ученых «Наука XXI веку» (II сессия) «IX Неделя науки МГТУ», г. Майкоп, 2004 г.; Всероссийской научно-практической конференции с международным участием «Инновационные технологии в создании продуктов питания нового поколения», г. Краснодар, 2005 г.; III Международной научно-технической конференции «Наука, техника и тех-

нология XXI века (НТТ-2007)», г. Нальчик, 2007 г.; Десятой международной научно-практической конференции «Современные проблемы техники и технологии пищевых производств», г. Барнаул, 2007 г.; IX Международной конференции молодых ученых «Пищевые технологии и биотехнологии», г. Казань, 2008 г.

Публикации. По результатам исследований опубликовано 55 работ, из них 2 монографии, 37 научных статей, в том числе 24 научных статьи в журналах, рекомендуемых ВАК, 11 работ в материалах конференций, получено 3 свидетельства и патента РФ на полезные модели и 2 решения о выдаче патента РФ на изобретения.

Объем работы. Диссертация состоит из введения, литературного обзора, девяти глав, объединенных в три части, основных результатов и выводов, списка литературы, включающего 339 наименований, и приложений. Работа изложена на 347 страницах, включает 97 рисунков и 41 таблицы.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Во второй главе проведен анализ процессов переноса массы и тепла в аппарате при контакте через половолоконную и трубчатую мембранную перегородку и примыкающим к ним внутренним и внешним областям аппаратов [1, 7-9].

Рассмотрена мембранная экстракция компонента из растительного масла при условии переноса массы через трубчатую липофоб-ную или трубчатую липофильную мембрану.

Рисунок 1 - Структурная схема исследования

На рисунке 2 приведены схемы противоточного массопереноса через трубчатую липофобную мембрану, а на рисунке 3 - схемы противоточного массопереноса через трубчатую липофильную мембрану.

НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ

ЛИПОФОБНАЯ НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ , МЕМБРАНА 1

Рисунок 2 - Схемы противоточного массопереноса через трубчатую липофобную мембрану: а - масло снаружи липофобной мембраны; б — масло внутри липофобной мембраны

НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ

ЛИПОФИЛЬНАЯ НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ, МЕМБРАНА ........................ ■

Рисунок 3 - Схемы противоточного массопереноса через трубчатую липофильную мембрану: а — масло снаружи липофильной мембраны; б - масло внутри липофильной мембраны

На рисунке 4 приведены схемы прямоточного массопереноса через трубчатую липофобную мембрану, а на рисунке 5 - схемы прямоточного массопереноса через трубчатую липофильную мембрану.

ЛИПОФОБНАЯ МЕМБРАНА

. НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ ЛИПОФОБНАЯ НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ . ■ МЕМБРАНА ■

___Гчч- '

Рисунок 4 - Схемы прямоточного массопереноса через трубчатую липофобную мембрану: а - масло снаружи липофобной мембраны; б - масло внутри липофобной мембраны

ЛИПОФИЛЬНАЯ МЕМБРАНА

I НАПРАВЛЕНИЕ МАССОПЕРЕДАЧИ

V г» " ЧШЯ—гЗ"

ЛИПОФИЛЬНАЯ МЕМБРАНА

а б

Рисунок 5 - Схемы прямоточного массопереноса через трубчатую липофобную мембрану: а - масло снаружи липофильной мембраны; б - масло внутри липофильной мембраны

На рисунке 6 приведена схема теплопереноса в половолокон-ной непористой мембране при различном направлении теплового потока.

НАПРАВЛЕНИЕ ТЕПЛОПЕРЕДАЧИ

Рисунок 6 - Схема теплообмена в половолоконной мембране при различном направлении теплового потока

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

+ (1) в случае прямотока в межтрубном пространстве:

Му{юх-соу+\)={с10у1(14), (2)

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

+ 1 )=-(*»,/</£). (3)

Распределение потенциалов в фазе трубного и межтрубного пространства мембран решается: для прямотока с помощью системы уравнений П) и (2):

= -(ых/ЛГ>, = К/{Мх + Л^И"'^ -1], (4)

для противотока с помощью системы уравнений (1) и (3):

со,

= [Л""'* - 1\/[\ - (му,

(5)

'У 1Г ~ ~ J. (6)

Общее решение с учетом перемешивания по обеим фазам для прямотока и противотока имеет вид:

/=1

(7)

где а, определяется из решения характеристического уравнения для прямотока (верхний знак) и для противотока (нижний знак):

1 -а1 ±

РехРеу

1 -±- 1

кРе,Му РеуНх

кРе, ЫхМу

Ре.

\ / \ 1 1

а -

Их ЛГ„ V * > /

= 0 (8)

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

чм"

Л2(х) '

Л(*1

ФМ] И-Ф(х)-Ф{х)+1

К-Ф(Х)

м,

N. N.

1-Ф(х)

N.

О

//•фоо

лг.

N.

Рисунок 7 — Матричное уравнение для нахождения коэффициентов уравнения (8)

Проведен анализ эффективности (б) для массо- и теплообмена по зависимости:

1 - ехр -

1--¿-ехр \ >

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

а б

Рисунок 8 - Изменение эффективности переноса массы и тепла при различных показателях чисел единиц переноса по фазе*, при отсутствии перемешивания и при различном Мх/Ыу: а) режим противотока; б) режим прямотока

Был проведен анализ изменения эффективности от факторов перемешивания в полном факторном эксперименте (таблица 1), где факторами явились параметры перемешивания по фазе х и по фазе у (Рех и Реу), имеющие значения 10 и 50, при значении числа единиц переноса >1Х = 2, при отношении чисел единиц переноса Му/Ых= 0,15 и при режиме противотока.

Таблица 1 - Влияние параметров перемешивания в факторном эксперименте

Значение параметра перемешивания Значение эффективности е

Рех Реу

-10 -10 0,77

+50 -10 0,82

-10 +50 0,77

+50 +50 0,82

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

В третьей главе рассмотрен сопряженный конвективный перенос тепла и массы от внутреннего потока с теплопроводностью и диффузией через мембранную перегородку [23-25].

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

Лцр СарТцф

3

Рж /'* Сж ХжТжад

I

'тпр

Рисунок 9 - Основные параметры, характеризующие сопряженный процесс переноса тепла потока жидкости внутри мембраны

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

УУЖ (г )--—-- =

дх

81Тх(г,х)1 1 дТж(г,х)| д2Тж (г, х)

дгг г дг дх1

(10)

д2ТтЛг,х) _ 1 дт (г,х) , дгтт(г,х)

тр\' 7 + ----"" У у +-тр I = 0 (11)

дг г дг дх2 и и

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

Для решения данной задачи методом конечных разностей определяем одномерную сетку (рисунок 10) по текущему радиусу г. Для упрощения последующих выкладок, узлы сетки определяем с равномерным шагом АЬ, значение которого определяется из условий сходимости и устойчивости явной схемы решения задачи.

i=Oj=MaxR _ Ы^МахЯ

/=1,/=)/?+!

ЫМЙ-1

1=0,1=1 ¡=1,1=1 $1=0,1=0 О М^О

е-

О

о

ЬМя*У.1 ¡-¡ПИ

(=махХ,рМахЯ1

/=Ma*X¿=/R+í

/=MaxXj=^K

ШахХ^М

рМахХ-и=1 I ¡=МахХ,]=1 01=ШхХ-1^0 ( ) /=МэхХ>0

'О

Рисунок 10 - Определение узлов сетки для области решения сопряженной задачи

Для расчета температурного поля ламинарного потока жидкости использовалось итерационное матричное уравнение следующего вида:

«0.0 «0.1 0 0 0 '/-1.0

«1.0 «1,1 «1.2 0 0 '/-и

0 ти 0 ■ 'и =

0 0 '/-1../Я-1

0 0

(12)

где коэффициенты матрицы т,^ определяются коэффициентами уравнений:

1м,Г

^ 4'Р.кон - Чнач + ^ ЧтрЧнач " ^ ЧтрЧкон + «^»ж^.ж^.тр

^^(Рнач-Р.ко«) VI

;.ДЬ2.(Клр2^2.ДЬ2).(Ркон-Рнач)

(13)

1-1,0"

В-Ах-а^ц^Ьгр

ЛЬ 2 ЯГр2-(Рнач - Ркон)

+1

1,0

8-Дх-а.^ц^Ц.р

.ДьЧр^Ркон-Рнач).

\1 (14)

и, ¿л1

Дн'/лЬ-Яг -Р

'Р -Р мК~2

+ ^^Ь^уРнач-2^3-Укон-*- ^У^Чр

¿ь3-(йь-г.1д(рк0н-рн^

на'

АЬ3(^-2.1д.(рнач-РКО1).0Я-1)

- 2-0(ряяч- Р^-да-1)

Матричное уравнение (12), содержащее разряженную ленточную матрицу, решалось методом декомпозиции с использованием пакета Ма&САЛ.

Рассмотрена постановка и решение сопряженной задачи конвективного массопереноса при течении жидкости в трубчатой мембране с диффузией в стенке мембраны (рисунок 11).

Рисунок 11 - Схема сопряженной задачи конвективного массопереноса при течении жидкости в трубчатой мембране с диффузией в стенке мембраны

В зависимости от смачиваемости мембраны ее поры заполняются соответствующей фазой. Характерными размерами трубчатой мембраны являются внешний (!*„) и внутренний радиус (Й^) мембраны. В этом случае для решения сопряженной задачи конвективной диффузии важно учитывать условия сопряжения конвективного потока внутри трубы и диффузионного потока в стенке мембраны.

Применяя метод конечных разностей, аналогично решенной задачи сопряженного теплообмена, построена схема решения задачи.

В качестве параметров варьирования использовались: температура потока жидкости в трубе, начальное давление потока жидкости в

трубе, длина трубы, коэффициент сопротивления мембраны, внутренний и внешний радиус трубы.

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

В результате была построена степенная модель, определяющая влияние изучаемых параметров процесса на интенсивность массопереноса:

N , = О,1932 • & 1 ' • К 00098 • К ;о отз ■ Ре "°'7 (16)

В четвертой главе представлены результаты исследования теплообмена в аппарате с непроницаемыми половолоконными мембранами [17, 25, 26,39].

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

К,мм

Рисунок 12 - Двумерное температурное поле сопряженной задачи теплообмена в половолоконной мембране

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

Представленные данные обобщены в виде степенной зависимости следующего вида:

Ыимем = 3,66 - 8,08 ■ Ре~°-2т ■ вголт ■ К?™ • (П)

Для проведения экспериментальных исследований переноса тепла в теплообменнике с непроницаемыми половолоконными мембра-

нами была создана экспериментальная установка, представленная на рисунке 13.

V

Рисунок 13 - Экспериментальная установка: 1 — теплообменник с полимерными половолоконными мембранами; 2 - термостат; 3 — приемный бак для воды; 4 — сливная воронка;

5 - записывающее устройство; Т1, Т2, ТЗ, Т4 — прибор для измерения температуры; Р1, Р2 - прибор для измерения расхода жидкости

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

Результаты проведенных исследований на экспериментальной установке приведены в таблице 2.

Таблица 2 - Результаты исследований на экспериментальной установке

Расход Расход Начальная Конечная Начальная Конечная

воды в воды температура температура температура температура

мембранах, в кожухе, в мембра- в мембра- в кожухе, в кожухе,

кг/сек кг/сек нах, К нах, К К К

0,288 0,266 293,4 299,5 310,1 303,5

0,284 0,367 294,4 300,7 309,7 304,8

0,174 0,605 294,1 301,1 304,7 302,7

0,244 0,694 294,2 . 301,4 306,8 304,3

0,297 0,835 293,6 300,9 308,8 306,3

0,255 0,874 293,6 301,1 307,8 305,6

0,304 1,338 294,1 302,0 309,6 307,8

0,292 1,423 293,8 301,9 308,9 307,3

Мембранные модули были предоставлены фирмой Е1с1о8 (Чехия) и представляют собой мембранные пучки, состоящие из 1400 полипропиленовых половолоконных мембранных волокон, наружным диаметром 0,275 мм, внутренним диаметром 0,225 мм и длиной 0,75 м.

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

С учетом определенных значений общего коэффициента теплопередачи и оценки значений коэффициента внутренней теплоотдачи по зависимости (18) получили значения коэффициентов внешней теплоотдачи.

На рисунке 14 приведена графическая зависимость внешнего теплообмена от степени упаковки волокна в модуле, влияние которой особенно существенно при низких значениях числа 11еМОд.

Полученные результаты обобщены зависимостью:

Ш.,^ =

4,3486-0,1301

1-0 0

О- 30,758-1,1401

1-0 0

Рг3

(18)

Среднее отклонение расчетных значений от экспериментальных значений составляет 11%.

1000 2000

6000 6000 7000

3000 4000 Пвмод

Рисунок 14 - Определение вида функциональной зависимости экспериментальных данных в группах с различной степенью упаковки волокон в модуле: 1 - ф=0,042; 2 - <¡>=0,052; 3 - ф=0,2

В пятой главе получены зависимости поперечной фильтрации для определения проницаемости половолоконных и трубчатых мембран и проведено моделирование переноса потоков в мембранных контакторах под действием трансмембранного перепада давления [13,16,20,37].

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

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

\А,.Л(4.[\.Х) (19)

Рь = Р, + Л, •

С учетом найденных выражений А1 и А2 из начальных и граничных условий реализации потоков, уравнение (19) принимает вид:

Pl-P. + 0я

4 ■

sh

■\Р,-Р, + ch

l.

■Ср.- я.)

(20)

Если давление снаружи мембраны и на выходе внутри мембраны совпадают {РК=РЕ), то уравнение (20) принимает вид:

Pl ' Р. + (P« - О-

sh 4 о

sh №4

(21)

Характер зависимости изменения давления при поперечной фильтрации определен при следующих параметрах: длина волокна Lm = 0,215 м; радиус волокна RL= 1,15-Ю"4 м; проницаемость мембраны гр = 6,2-10' 5м; давление на входе в мембрану Рн = 2,026-105 Па (2 бар); давление на выходе из мембраны Рк= 1,52-105 Па (1,5 бар); давление снаружи мембраны Ps= 1,013-Ю5 Па (1 бар).

Как видно из представленных данных (рисунок 15) давление по длине мембраны меняется линейно как для случая при отличающихся значениях давления на выходе из мембраны и внешнего давления, так и при совпадающих значениях давления.

0.1 0.15 0.2 0,25

X, м

Рисунок 15 - Изменение давления при поперечной (1,2) и тупиковой (3) фильтрации в половолоконной мембране: 1 - Р* Ф Р5; 2 - Рк = Р5; 3 - тупиковая фильтрация

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

ей

ей

(22)

С учетом найденных постоянных интегрирования и уравнений объемных расходов поперечной фильтрации во входном и выходном сечении, получаем разность объемных расходов во входном и выходном сечениях (24), которую использовали для определения значения проницаемости полипропиленовых половолоконных мембран фирмы

ЕГООБ (Чехия) РР - М5 г„ = 4,4- 1(Г15м и РР - Мб гр = 1,33- 1(Г V

д е =

Г ГГ~ 1

ей -1

2 • ц ■ як

■I.

(23)

Используя методику тупиковой фильтрации с помощью уравнения (25), определили значение проницаемости трубчатых керамических мембран произведенных фирмой НПО «Керамикфильтр» (Москва) - гр = 7,43- 10"13м.

г-ц-ск

(24)

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

1 + у ■ сИ

( ЛхУ

1тг)

у - у - ей

У + 1 Я • х

!Г

В, +

X ■ х + у

у + 1

в,

я(г + 0

у Я х- у -УГ'^ А-(г + 1)

Я

17П

•В4 (25)

Я ■ х - •

Я-0- + 0

/■Я -х + ^[у -яИ

Эти уравнения содержат четыре постоянных интегрирования, что совпадает с числом точек контроля давлений на схеме (рисунок 16).

Рисунок 16 - Схема давлений в точках мембранного модуля

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

\ = Рш (27)

{(Рр ~ Рю,(Рив -

(28)

В,=РШ (29)

1 ру {{Рив -Р„н)-ску]-{РКИ -ркв)}

(30)

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

Используя полученные уравнения, произведен расчет давлений для керамической мембраны с проницаемостью, изменяющейся от гр= 7,4-Ю"13 м до гр= 7,4-Ю"9 м, с длиной мембраны 1^,= 800 мм, внутренним радиусом мембраны = 2 мм, внешним радиусом мембраны Ям = 5 мм, радиусом картриджа II = 16 мм, вязкостью жидкости р = 8,9-10"4 Па с, начальным давлением Рн= 0,2 МПа и конечным давлением Рк= 0,1 МПа.

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

200000*

180000

_Г О.

.160000

140000

120000

100000

г

¡6

д Ч

а а й 1 )

"3 ! \

Л О * 5

1 к с 1 ж

1 1 1 н

А i 1 к

* 1 Л й *

с ' А , > «

100

200

300

600

700

400 500 X, мм

Рисунок 17 - Изменение давления потока по оси мембраны Значения проницаемости: © - 7,4-Ю"13; ® - 7,4-Ю"12;® - 7,4-10"п;

@ - 7,4-Ю'10; ® - 7,4-Ю'9

800

26, 21 16

га

С11 <о

а. <в

1 Л , к Л 4 к

1 к А

А

Л к 1 Л

1 а с > а г > а с 1 а г > а 1 > а г 1 1 А

£ « ) Е 5 Е # 5 в Л 1 К Й! I с Й? % * ш \ < Я* I

100

200

300

400

500

600

700

X, мм

Рисунок 18 - Избыточное давление фильтрата по оси мембраны Значения проницаемости: © - 7,4-Ю"13; <8> - 7,4-Ю'12;® - 7,4-Ю-11;

© - 7,4-Ю"10; @ - 7,4-Ю"9

Графики изменения расхода потока и расхода фильтрата при различных величинах проницаемости представлены на рисунках 19 и 20.

зооо

2500

2000

3

3

-Г

о

1500

1000

500

100

ТЯн

«

д

Н В П

' • д Л

200

300

400 500 X, мм

600

700

800

Рисунок 19 - Изменение расхода потока по оси мембраны Значения проницаемости: © - 7,4-Ю"13; <8> - 7,4-Ю'12;® - 7,4-10"";

© - 7,4-Ю'10; ® - 7,4-Ю"9

л Л ' к V & '

* к

Л X

1 >

1 г в « < 5 1 , а с 1 а ( 1 а г 1 □ | ) о I I а и

(1 1 Ю 2( 0 30 ) 40 ) 50 1 60 ) 70 > 80(>

Рисунок 20 - Изменение расхода фильтрата по оси мембраны Значения проницаемости: ® - 7,4-Ю"13; <8> - 7,4-Ю'12;® - 7,4-Ю'11;

©- 7,4-10"10;® -7,4-10"9

В шестой главе определены данные по растворимости и коэффициентам диффузии в системе «триацилглицерины - жирные кислоты - двуокись углерода», необходимые для определения параметров массообмена при экстракции в массообменнике с объемом, структурированном мембранными перегородками [1,4-6, 22,29, 33,40, 41].

На основании уравнения состояния Пенга-Робинсона и уравнения растворимости Крестила была создана модель растворимости растительного масла в двуокиси углерода.

Результаты моделирования представлены на рисунке 21.

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

Для предсказания бинарных коэффициентов диффузии использовалась зависимость:

А, = 1,4974-10"9 [(44010/ РсоУ ~ 9,7243^/(77 М2), (31)

где к = (2 / 3)[1 - 0,28 ехр (- 4,25 рСОг)]

Рисунок 21 - Результаты моделирования растворимости нерафинированного подсолнечного масла в двуокиси углерода: 1-Т = 302К; 2-Т = 311 К; 3-Т = 320К

В результате использования уравнения состояния Пенга-Робинсона и уравнения (31) проведено моделирование поведения коэффициентов диффузии в сверхкритической области давлений и температур.

На рассматриваемой диаграмме (рисунок 22) выделяется гребень, так называемая «квазиспинодаль», ось которого проходит через критическую точку и может быть представлен на рисунке линейным уравнением.

Температура, К

Рисунок 22 - Изменения коэффициентов диффузии жирных кислот при бесконечном разбавлении в С02: область А - значения коэффициентов диффузии увеличиваются с ростом

температуры и уменьшаются с уменьшением давления; область Б - значения коэффициентов диффузии не зависят от температуры и значительно уменьшаются с уменьшением давления

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

С учетом специфики растворяемого вещества в уравнении (31), определяемое молекулярным весом, предложено универсальное для жирных кислот уравнение для расчета коэффициентов диффузии:

£>° 10 6 = (0,0696 -Т + 18 ,7)/ у/М 2 (32)

В работе получены значения коэффициентов распределения для жирных кислот в двуокиси углерода, которые представлены в таблице 3 для режимов (Р и Т), выбранных вдоль оси гребня на диаграмме (рисунок 19).

Таблица 3 - Значения коэффициента распределения в двуокиси углерода

(у*») для различных жирных кислот

Наименование жирной кислоты Режим 1 Режим 2 Режим 3 Режим 4

(Р = 80бар, Т = 308К) (Р = 120бар, Т = 329К) (Р = 160бар, Т = 351К) (Р = 200бар, Т = 373К)

Стеариновая кислота 2,59-10"4 9,84-103 0,130 1,12

Олеиновая кислота 0,326 2,172 8,253 25,7

Линолевая кислота 0,529 2,75 8,94 26,6

Получены уравнения для анализа массопереноса в мембранном контакторе: по фазе х:

О, д9,

1 дгр,

дг!2

по фазе у\

12*91.

1 д2р}

8Ъ

(33)

(34)

Общее уравнение массопередачи представлено в виде обобщенных переменных при базовых параметрах переноса массы в фазе х:

1 = 1 + 1 + 1 .

Бк а, + 5Й.Г,0 + ' ( )

Для нахождения зависимости числа Шервуда по фазе х и по фазе у от числа единиц переноса по фазе л: и отношения чисел единиц переноса фазы у к числу единиц переноса фазы х решали краевые задачи (33) и (34), определяли среднее число Шервуда вдоль продольной координаты, усредняли и проводили итерационные вычисления. В результате получили зависимости, приведенные на рисунках 23 и 24.

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

■Н

гнср,'

з.оо: 2,00. 1,00 0,00

8П,

8,00., 7,00.1 6.оо1

5.00| 4,0« ■ср» 3,00 2,00 1,00 0,00

мшщяр -

0.27 о,39

0,51 0,63 0,75 '

Му/Ых

Рисунок 23 - Зависимость 8Ьорх от числа единиц переноса по фазе х и отношения числа единиц переноса по фазе х к числу единиц переноса по фазе .у

' 0.27 0.39 Ыу/М;

0,51 0,63

0,75

Рисунок 24 - Зависимость 8Ьсру от числа единиц переноса по фазе х и отношения числа единиц переноса по фазе х к числу единиц переноса по фазе у

В качестве сырья использовали нерафинированное прессовое подсолнечное масло со следующими исходными показателями: кислотное число - 3,79 мг КОН/г (содержание свободных жирных ки-

слот 1,89 %); перекисное число - 7,85 ммоль активного кислорода/кг; массовая доля фосфолипидов - 0,65 %, анизидиновое число - 0,63.

В таблице 4 приведены значения показателей, полученных в результате экспериментов, а в таблице 5 представлена обработка экспериментальных данных.

Таблица 4 - Результат проведенных экспериментальных исследований экстракции жирных кислот из нерафинированного подсолнечного масла двуокисью углерода

№ опыта Характеристики режима экстракции Расход масла, л/ч Расход С02, л/ч КЧ масла на выходе, мг КОН/г

Т,К Р, бар

1 308 80 2,0 20 0,10

2 308 80 2,5 20 0,63

3 308 80 3,0 20 1,46

4 328 120 14,0 20 0,10

5 328 120 16,0 20 0,32

6 328 120 . 18,0 20 0,74

Таблица 5 -Результаты обработки экспериментальных данных

№№ опытов е N. 5ЬЭКСП ^Цзасч

1 0,975 0,307 4,797 7,472 7,487

2 0,834 0,383 2,287 4,442 4,443

3 0,614 0,460 1,148 • 2,679 2,673

4 0,975 0,322 4,874 7,747 7,516

5 0,917 0,368 3,286 5,787 5,759

6 0,805 0,414 2,098 4,155 4,138

Отмечается близость результатов экспериментов и расчетов по предложенной модели.

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

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

В седьмой главе рассмотрена непрерывная схема работы мембранного реактора для производства биотоплива и ее обоснование на основе математического моделирования [2,10,15, 18, 21, 43].

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

Схема циркуляционного контура с мембранным реактором представлена на рисунке 25.

Рисунок 25 - Схема циркуляционного контура с мембранным реактором

Непрореагировавшая часть реагента остается в капельном состоянии и, предварительно подогретая в теплообменнике, рециркули-рует на вход с помощью циркуляционного насоса.

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

Математическая модель мембранного реактора представляется тремя взаимосвязанными блоками: модель кинетики реакции, модель вязкости, модель гидравлики.

Масляные капельки

Смесь нейтрального —рНр?— I? "т° "V){"? масла, метанола и (I \\

катализатора Ь ¿ а> д> о> >рЩ

Непрореагированная часть реагента

а-

Продукты реакции, содержащие биотопливо

Рисунок 26 - Схема мембранного реактора для производства биотоплива

Из модели кинетики получаем информацию о составе реакционной смеси, что является основой определения вязкости. Имея информацию о вязкости, получаем оценку разделения потоков на мембранах и уточняем время процесса.

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

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

Для пересчета вязкости по температурной шкале известной вязкости чистого компонента использовали метод Льюиса и Сквайрса, который представлен аналитической зависимостью (рисунок 27).

Зависимость вязкости компонента чй от температуры, которая определяется как сумма опорной температуры 90 и разности между расчетной температурой I и температурой известной вязкости определяется по уравнению:

1](() = ехр[а • [1п(0о + (*-0)]21п(б>0 + (*~0)+с] (36)

6 4 2. О

-о" 4

•5 -2 -5

-4

-6 -8 -10

Рисунок 27 - Логарифмированный график Льюиса и Сквайрса с аналитической аппроксимацией 1п 7]0 = а •

Для проверки метода получения температурных зависимостей вязкости сравнили известные значения вязкости при различных температурах с расчетными значениями, что представлено на рисунке 28.

Рисунок 28- Зависимости экспериментальных и расчетных данных вязкости от температуры для триацилглицеринов (ТАГ), глицерина (Г) и метанола (М): 1 - глицерин (♦ опыт; - расчет); 2 - триацилглицерины (• опыт; — расчет);

3 - метанол (А опыт; - расчет)

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

N

1п(0 = Х[М>1п(?7,)] (37)

(-1

ч-

а 4 - ¿V ( 0 5 5 С 0 6

Щво)

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

Проведено математическое моделирование мембранного реактора для производства биодизельного топлива из растительных масел.

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

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

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

При совпадении значения вязкости и времени с точностью до 0,1% итерации прекращались, и результат численного моделирования считался окончательным.

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

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

Давление * модуле (Ры), бар Рисунок 29 - Выход целевого продукта в зависимости от числа фильтрующих элементов: 1 - пМ = 4; 2 - пМ = 5; 3 - пМ = 6; 4 - пМ = 7

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

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

-^ = 6,569-Ю-4 -ф'}лш -Л°'И6 'Кт-°'ш* .К/-т\ (39)

^м

где безразмерные критерии:

относительный перепад давления по

Ф — - степень заполнения поперечного сечения мембра-

нами;

. 1пКьп1

Ау ~ —— - удельная площадь мембран в единице объема К

модуля;

% _ п&1п(Рцв ~ Ркв) _

длине мембраны;

К, =--———-г-- - относительный перепад

давления поперек мембраны.

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

В восьмой главе рассмотрен теплообмен при конденсации парогазовых смесей и обоснована схема конденсатора с полимерными половолоконными мембранами в линии комплексной переработки растительных масел [11,12,14,32,36,38,42].

На рисунке 30 представлен принцип конденсации парогазовых смесей с полимерными половолоконными мембранами с внутренним источником тепла.

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

т.

Т.У

т*.

Поверхность теплообмена

1 . парогазовая смесь

2 - холодильный агент

б)

Рисунок 30 - Принцип конденсации парогазовых смесей: а) принципиальная схема; б) температурные профили потоков

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

Математическая модель процесса конденсации включает массовые и энергетические потоки в дифференциальном объеме мембранного конденсатора, представленные на рисунке 31 и поперечные профили температуры от смеси «газ-пар», представленные на рисунке 32.

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

9-йФ

П1 -

У-йу '

в

Парогаюмя смесь

пленка конденсата

Уи ^ Парка*««» - сиось

«дд Пленка •^конденсата ППНМ

хладагент И Стенка

Рисунок 31 - Массовые и энергетические потоки в дифференциальном объеме мембранного конденсатора

Рисунок 32 - Поперечные профили температуры от смеси «газ-пар»

На рисунке 33 приведена конструкция мембранного конденсатора

М5 ^-12 М8 '—19 20^-17 4 \—

Рисунок 33 - Мембранный конденсатор: 1 - цилиндрический корпус; 2 - фланец для подвода исходной смеси; 3 - фланец для отвода газовой фазы; 4 - фланец для отвода конденсата паровой фазы; 5,6 - крышка; 7, 8 - фланец для подвода исходной смеси; 9 - стержень; 10,11 - гайка; 12 - трубный пучок; 13,14 - распределительные

камеры; 15 - фланец для подвода хладагента; 16 - фланец для отвода хладагента; 17 - полипропиленовый непористый половолоконный мембранный пучок; 18,19,20 - поперечные перегородки; 21,22,23,24 - фиксаторы

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

В девятой главе обоснована технологическая схема и линия комплексной переработки растительных масел с применением мембранных аппаратов[43,44].

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

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

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

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

Полученное при обезжиривании фосфолипидов масло является сырьем для получения биотоплива.

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

Обезжиренные фосфолипиды

Теплоноситель N«.32 Т«плоно :ит«льза

Рисунок 34 - Технологическая схема комплексной переработки растительных масел: 1 - резервуар для нерафинированного масла; 2,5,8,10,16,18,24,25,27,30,33,36,39,40,42-насосы; 3-мембранный экстрактор; 4-резервуар для двуокиси углерода; 6,9,26,37,41 -выпарные аппараты; 7,38-конденсаторы; 11-теплообменник; 12-инжектор; 13-экспозитор; 14-отстойник для фосфолипидной эмульсии; 15 - вакуум - сушильный аппарат; 17-резервуар для фосфолипидной эмульсии; 19 - ротационно-пленочный сушильный аппарат; 20 - экстрактор; 21 - резервуар для ацетона; 22 — чанный испаритель; 23 - мембранный парогазовый конденсатор; 28 - сборник смеси; 29 - резервуар для метанола; 31 - мембранный реактор; 32 -датчик уровня; 34 - теплообменник; 35 - отстойник

Таблица 6 - Физико-химические показатели рафинированного подсолнечного масла, обезжиренных фосфолипидов, жирных кислот и биотоплива

Значение показателя

Наименование показателя Традиционная Разработанная

схема схема

Рафинированное масло

Кислотное число, мгКОН/г 0,25 0,10

Перекисное число, ммоль активного кислорода / кг 4,75 1,25

Массовая доля фосфолипидов, % отсутствие отсутствие

Выход рафинированного масла, % к исходному нерафинированному маслу 95,37 97,34

Обезжиренные фосфолипиды

Массовая доля фосфолипидов в продукте, % 94,50 98,70

Перекисное число, ммоль активного кислорода/ кг 3,78 0,25

Выход обезжиренных фосфолипидов, % к их содержанию в исходном нерафинированном масле Выход обезжиренных фосфолипидов, % к исходному нерафинированному маслу 53,84 0,35 84,61 0,55

Жирные кислоты

Массовая доля жирных кислот в продукте, % 95,30 99,15

Выход жирных кислот, % к их

содержанию в исходном нерафинированном масле 80,50 97,35

Выход жирных кислот, % к исходному нерафинированному маслу 1,48 1,85

Биотопливо

Кислотное число, мгКОН/г - 0,10

Температура вспышки, °С - 118

Плотность при 15 °С, г/см3 - 0,883

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

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

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

6. Зависимости для поперечной и тупиковой фильтрации позволяют определить по результатам экспериментов определения проницаемости мембран, которые в исследованном диапазоне (избыточное давление до 1 бар) практически не зависели от давления и соответственно составили для полипропиленовых половолоконных мембран фирмы ЕГООБ (Чехия): РР-М5 - 4,4-10'15 м; РР-М6 - 1,33-10'14м и керамических трубчатых мембран произведенных фирмой НПО «Кера-микфильтр» (Москва) - 7,43-10"13м.

7. Растворимость нерафинированного растительного масла в двуокиси углерода в сверхкритической области растет с ростом давления

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

8. Для коэффициентов диффузии жирных кислот в двуокиси углерода в диапазоне изменения давления от 80 до 400 бар и температуры от 293 до 373 К выделяются две области, разделенных линейной зависимостью, проходящей через критическую точку. Ниже этой линии значения коэффициентов диффузии практически не зависят от изменения значения температуры и значительно снижаются со снижением значения давления. Выше этой линии значения коэффициентов диффузии увеличиваются с ростом температуры и уменьшаются при снижении давления.

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

10. Вязкость реакционной смеси определяется вязкостью ее компонентов и составом. Температурная зависимость вязкости компонентов определяется по усовершенствованному методу Льюиса и Сквайрса.

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

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

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

14. Построенная математическая модель конденсации парогазовой смеси в конденсаторе с полимерными половолоконными мембранами может быть использована для моделирования и анализа работы конденсатора.

15. Обоснована комплексная технология и система процессов переработки нерафинированных растительных масел с новым процессом физической рафинации мембранной экстракцией и экстракционной очисткой фосфолипидов в схеме, в которой применен мембранный конденсатор парогазовой смеси, а отделенное масло перерабатывается в мембранном реакторе с получением биотоплива. Разработаны и переданы для реализации в производство филиала «Лабинский МЭЗ» ООО «МЭЗ Юг Руси» комплексная технология переработки нерафинированного подсолнечного масла.

ОБОЗНАЧЕНИЯ

т - т -

V* - 7--<РУ= 7-\ - безразмерные кон-

\Ртхо - УупУо) Ууг.Уо)

а с а - (*х~(хо) д _

центрации с обеих сторон мембраны; - т-и - г—_ \ -

\fjr0 ~ *у0/ \,'т0 *у<>)

безразмерная температура с обеих сторон мембраны; = {крц/хт)/<2х - числа единиц переноса массы по фазе в трубном и межтрубном пространстве; Нтх = (ктр)/ 1УХ и Л^ = (ктр)/1¥у -числа единиц переноса тепла по фазе в трубном и межтрубном пространстве; у/хт = х*/х„ и у/ущ = у*/у„ - коэффициенты распределения растворенного вещества; х*, х„ и у*, у„ - равновесные концентрации, относящиеся к условиям разных сторон поверхности мембраны между фазами; ?хи¥у- водяные эквиваленты обоих теплоносителей; ()х и - объемные расходы фазы в трубном в межтрубном пространстве; Ь - длина мембран; - площадь поверхности мембраны; 5 - сечение колонны; К — коэффициент массопередачи; <ах и соу - безразмерное значение потенциалов переноса массы и тепла в мембранных трубках и межтрубном пространстве; КТ - коэффициент теплопереда-

ту_ „ ^

чи; "ех = —и Реу = —- - безразмерные модифицированные

числа Пекле, характеризующие перемешивание по каждой из обеих фаз; Д, и коэффициенты продольного перемешивания по каждой из обеих фаз; Тж и Ттр - соответственно температура потока жидкости в трубе и непроницаемой стенки трубы, К; Ршч и Ркон - давления в начале и конце трубы; X«, Сж, рж, цж - соответственно, теплопроводность, теплоемкость, плотность и вязкость потока жидкости; ^ и Яст - радиусы трубы и стенки, м; Ьтр - длина трубы, м; С^, р^ -соответственно теплопроводность, теплоемкость, плотность трубы;

f P -P

иан кон m

Щ

mp

Pe = Pr-Re^/2 - число Пекле; Remax =' ,

Мж ' Pж

осевое число Рейнольдса потока жидкости; PL; Ps и QL; Qs - соответственно давления и расходы во внутреннем и внешнем объеме между мембранами; kXfS и kx¡L - компоненты тензора проницаемости Дарси; ц - динамическая вязкость, Па с; гр - мембранная проницаемость, м; Rl - внутренний радиус волокна, м; х - координата по оси мембраны, м; п - количество мембран; NuCT - критерий Нуссельта стенки, учитывающий объединенное сопротивление тепловому потоку в стенке и

„ г. w-dj-p „ _ „

внешней конвективной теплоотдачи; Re =-1— - критерии Реи-

м

нольдса; ф - степень упаковки половолоконных мембран в кожухе Ñ-u

модуля Рг = —---критерий Прандтля; г|0 - опорная вязкость чис-

А

того компонента; 90 - опорная температура чистого компонента, °С; ТАГ - триацилглицерины; М - метанол; Г - глицерины.

Основное содержание диссертации изложено в следующих работах:

Монографии:

1. Схаляхов, А.А. Мембранная экстракция двуокисью углерода в пищевой технологии/А.А. Схаляхов; ред. журн. «Известия ВУЗов. Пищевая технология». - Краснодар, 2007. - 159 е.: ил.- Библиогр.: с. 157-159 (141 назв.). - Деп. в ВИНИТИ 12.02.07, №128-В2007.

2. Схаляхов, А.А. Производство биотоплива из масел и жиров: монография / А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой; Федер. агентство по образованию, ГОУ ВПО МГТУ. - Майкоп: МГТУ, 2008. - 132 с.

Научные статьи в журналах, рекомендуемых ВАК:

3. Оценка развития научного направления «Экстракция двуокисью углерода»/А.А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю. Чундышко//Известия ВУЗов. Пищевая технология. -1999.-№1.-С. 8-11.

4. Поведение коэффициентов диффузии жирных кислот в СОг при сверхкритических условиях/А.А. Схаляхов, Е.П. Кошевой, Д.Э. Тихонов, Х.Р. Блягоз/УИзвестия ВУЗов. Пищевая технология. - 1999. -№2-3.-С. 72-73.

5. Универсальная установка для экстракции двуокисью угле-рода/А.А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю.

Чундышко//Известия ВУЗов. Пищевая технология. - 1999. - №4. -С. 67-69.

6. Схаляхов, A.A. Коэффициенты диффузии триглицеридов и жирных кислот в двуокиси углерода при сверхкритических усло-виях/АА. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз//Известия ВУЗов. Пищевая технология. - 2000. - №2-3. - С. 62-63.

7. Теоретический анализ экстракции в массообменнике с пористой перегородкой/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев, Х.Р. Блягоз// Известия ВУЗов. Пищевая технология. - 2001. - №5-6. - С. 66-68.

8. Схаляхов, A.A. Массообмен при экстракции в аппаратах с пористой перегородкой/А .А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз// Известия ВУЗов. Пищевая технология.- 2003.- №5-6.- С. 82-84.

9. Схаляхов, A.A. Влияние продольного перемешивания во взаимодействующих фазах при десорбции в пленочном трубчатом противоточном аппарате/А.А. Схаляхов, А.Г. Верещагин, Е.П. Кошевой // Известия ВУЗов. Пищевая технология.- 2007,- №2.- С. 71-73.

10. Схаляхов, A.A. Организация процесса переэтерификации в мембранном реакторе/А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой// Новые технологии. - 2008. - Вып. 5. - С. 34-36.

11. Конденсатор для парогазовых смесей с полимерными поло-волоконными мембранами/А.А. Схаляхов, ЕЛ.Кошевой, B.C. Косачев, А.Г. Верещагин//Известия ВУЗов. Пищевая технология-2009.-№2-3 С. 68-70.

12. Конденсация парогазовых смесей с полимерными половоло-конными мембранами/А.А. Схаляхов, Х.Р. Блягоз, А.Г. Верещагин, Е.П. Кошевой//Новые технологии. - 2009.- Вып. I. - С. 35-39.

13. Определение проницаемости половолоконных и трубчатых мембран/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев, Е.О. Никонов// Известия ВУЗов. Пищевая технология. - 2009. - №2-3 С. 96-98.

14. Разработка модели конденсации парогазовых смесей с полимерными половолоконными мембранами/А.А. Схаляхов, А.Г. Верещагин, B.C. Косачев, Е.П. Кошевой//Новые технологии. - 2009.-Вып. 1.-С. 39-43.

15. Схаляхов, A.A. Вязкость компонентов реакционной смеси при производстве биодизельного топлива из растительных масел/А.А. Схаляхов, B.C. Косачев, Е.П. Кошевой//Известия ВУЗов. Пищевая технология - 2009. - №1. - С. 113-115.

16. Схаляхов, A.A. Математическое моделирование процесса разделения жидких смесей в мембранном модуле с различной организацией потоков/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев// Известия ВУЗов. Пищевая технология.- 2009. - №2-3. - С. 71-74.

17. Теплообмен с полимерными поволоконными мембранами/ A.A. Схаляхов, Е.П. Кошевой, B.C. Косачев, А.Г. Верещагин // Известия ВУЗов. Пищевая технология.- 2009. - №2-3.- С. 79-81.

18. Схаляхов, A.A. Математическое моделирование кинетики и свойств реакционной смеси при производстве биодизельного топлива из масел / A.A. Схаляхов // Новые технологии. - 2009. - №3. - С, 29-35.

19. Схаляхов, A.A. Обоснование системы процессов получения фосфолипидных продуктов и биодизеля при переработке отходов гидратации масла/А.А. Схаляхов, А.Г. Верещагин, Е.А. Бутина//Новые технологии. - 2009. - №3. - С. 39-42.

20. Схаляхов, A.A. Математическая модель гидравлики мембранного реактора в линии производства биодизеля/А.А. Схаляхов// Новые технологии. - 2009. - №3. - С. 39-42.

21. Схаляхов, A.A. Результаты моделирования работы мембранного реактора для производства биодизельного топлива из масел/

A.A. Схаляхов//Новые технологии - 2009. - №3.- С. 42-48.

22. Схаляхов, A.A. Оценка эффективных режимов работы массо-обменника с пористой перешродкой/А.А. Схаляхов//Новые технологии. - 2009. - №3. - С. 48-54.

23. Схаляхов, A.A. Постановка сопряженной задачи переноса тепла через трубчатую мембрану от протекающего в ней потока/А.А. Схаляхов//Новые технологии. - 2009. - №4. - С. 53-56.

24. Схаляхов, A.A. Численное решение сопряженной задачи переноса тепла через трубчатую мембрану/А.А. Схаляхов//Новые технологии. - 2009. - №4. - С. 57-62.

25. Описание теплообмена в половолоконных мембранах для внутренней сопряженной задачи/А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой, B.C. Косачев//Новые технологии. - 2010. - №1. - С. 54-59.

26. Перенос тепла в теплообменнике с непроницаемыми поло-волоконными мембранами/А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой,

B.C. Косачев//Новые технологии. - 2010. - Ksl. - С. 59-64.

Научные статьи и отчеты:

27. Двуокись углерода при сверхкритических условиях как растворитель для экстракционной очистки жиров/А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой, Д.Э. Тихонов//Процессы, аппараты и машины пищевой технологии: межвуз. сб. науч. тр. / М-во образования Рос. Федерации, С.-Петерб. гос. акад. холода и пищевых технологий. -СПб., 1999.-С. 153-159.

28. Анализ технологии и техники экстрагирования двуокисью углерода при сверхкритических условиях/А.А. Схаляхов, Е.П. Коше-

вой, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю. Чундышко // Новые технологии: [сб. науч. тр. МГТУ] / М-во образования и науки Рос. Федерации, Майкоп, гос. технол. ун-т. - Майкоп, 2005. -С. 178-180.

29. Массоперенос в массообменнике с пористой перегородкой (мембраной) / A.A. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов // Новые технологии: [сб. науч. тр. МГТУ] / М-во образования и науки Рос. Федерации, Майкоп, гос. технол. ун-т. - Майкоп, 2005.-С. 175-177.

30. Технология и техника экстрагирования двуокисью углерода при сверхкритических условиях/А.А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю. Чундышко//Экономика и технологии. -2005. -№1.- С. 51-54.

31. Эффективность процессов экстракции и разделения в системе «триацилглицерины-двуокись углерода-этанол»/А.А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю. Чундышко//Новые технологии: [сб. науч. тр. МГТУ] I М-во образования и науки Рос. Федерации, Майкоп, гос. технол. ун-т. - Майкоп, 2005. - С. 181-183.

32. Разработка конструкции конденсатора для смесей растворителя с неконденсирующимся газом/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев, A.B. Гукасян, А.Г. Верещагин//НовЫе технологии/М-во образования и науки Рос. Федерации, Федер. Агентство по образованию, Майкоп, гос. технол. ун-т. - Майкоп, 2006. - Вып. 2. - С. 90-92.

33. Схаляхов, A.A. Коэффициенты распределения жирных кислот в системе с двуокисью углерода при сверхкритических условиях / A.A. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз//Научные основы процессов, аппаратов и машин пищевых производств: междунар. науч.-практ. конф., Россия, Краснодар, 24-26 сент. 2002 г.: [материалы]. - Краснодар, 2002.-С. 29-31.

34. Схаляхов, A.A. Моделирование растворимости триглице-ридов в смешанном растворителе при сверхкритических условиях/ A.A. Схаляхов, В.Ю. Чундышко, A.JI. Яковлев//1Х Неделя науки МГТУ: V Всерос. науч.-практ. конф. студентов, аспирантов, докторантов и молодых ученых «Наука XXI веку» (II сессия). - Майкоп, 2004.-С. 125-129.

35. Технология получения фракционированных фосфолипидных продуктов/Е.А. Бутина, A.A. Схаляхов, В.В. Сорокина, E.H. Пахо-мова // Инновационные технологии в создании продуктов питания нового поколения: материалы Всерос. науч.-практ. конф. с междунар. участием. - Краснодар, 2005. - С.183-185.

36. Конструкции конденсатора для смесей паров ацетона с неконденсирующимся газом азотом/А.А. Схаляхов, А.Г. Верещагин,

Е.П. Кошевой, B.C. Косачев, A.B. Гукасян//Наука, техника и технология XXI века (НТТ-2007): материалы III Междунар. науч.-техн. конф. (Нальчик, 4-6 октября 2007 г.) / [редкол.: Батыров У. Д. (отв. ред. и др.)]. - Нальчик, 2007. - Т. 1. - С. 221-224.

37. Схаляхов, A.A. Математическая модель гидравлики реактора с пористыми мембранами для производства биодизеля/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев//Современные проблемы техники и технологии пищевых производств: материалы Десятой междунар. науч.-практ. конф. (11-12 декабря 2007 г.). - Барнаул, 2007. - С. 347-349.

38. Моделирование процесса конденсации смеси паров ацетона с инертным газом в парогазовом конденсаторе/А.А Схаляхов, А.Г. Верещагин, Е.П. Кошевой, B.C. Косачев//Современные проблемы техники и технологии пищевых производств: материалы Десятой междунар. науч.-практ. конф. (11-12 декабря 2007 г.). - Барнаул, 2007.-С. 339-342.

39. Теплообменные характеристики половолоконных непористых мембран/А.А. Схаляхов, А.Г. Верещагин, Е.П. Кошевой, A.B. ГукасяшУНаука, техника и технология XXI века (НТТ-2007): материалы III Междунар. науч.-техн. конф. (Нальчик, 4-6 октября 2007 г.)/ [редкол.: Батыров У. Д. (отв. ред. и др.)]. - Нальчик, 2007. - Т. 1.-С. 216-221.

Патенты и свидетельства:

40. Универсальная установка для экстракции двуокисью углерода / A.A. Схаляхов, Е.П. Кошевой, Х.Р. Бляшз, Х.Р. Сиюхов, В.Ю. Чундышко // Свидетельство на полезную модель №16458. - 2001.

41. Экстрактор/А.А. Схаляхов, Е.П. Кошевой, Х.Р. Блягоз // Свидетельство на полезную модель №16503. - 2001.

42. Конденсатор/А.А. Схаляхов, Е.П. Кошевой, B.C. Косачев, А.Г. Верещагин, A.B. Гукасян/ЛТатент на полезную модель №61401.2006, БИ №6, 2007.

43. Линия комплексной переработки растительного масла. Решение о выдаче патента РФ на изобретение по заявке №2009141410 от 11.11.2009 г. Схаляхов A.A., Блягоз Х.Р., Кошевой Е.П., Бутана Е.А., Верещагин А.Г.

44. Линия получения биотоплива из растительного масла. Решение о выдаче патента РФ на изобретение по заявке №2009141411 от 11.11.2009 г. Схаляхов A.A., Блягоз Х.Р., Кошевой Е.П., Бутана Е.А., Верещагин А.Г.

Материалы конференций:

45. Схаляхов, A.A. Оценка коэффициентов диффузии различных веществ в двуокиси углерода при сверхкритических условиях / A.A.

Схаляхов, Х.Р. Блягоз, Е.П. Кошевой//Тез. докл. Сев.-Кавк. регион, науч. конф. студентов, аспирантов и молодых ученых «Перспектива-98», 24-27 алр. 1998г./[редкол.: Шаков Х.Х. (отв. ред.) и др.]. -Нальчик, 1998.-С. 10-11.

46. Разработка оборудования для экстракции газами при сверхкритических давлениях/Х.Р. Блягоз, Е.П. Кошевой, Х.Р. Сиюхов, A.A. Схаляхов, В.Ю. Чундышко//Материалы II Междунар. науч. конф. «Продовольственный рынок и проблемы здорового питания»: тез. докл. - Орел, 1999. - С. 262.

47. Энергообеспечение экстракции газами на основе теплового насоса/ЕЛ. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, A.A. Схаляхов, В.Ю. Чундышко//Индустрия продуктов здорового питания - третье тысячелетие (человек, наука, технология, экономика): междунар. науч.-практ. конф., 24-25 февр. 1999 г.: тез. докл.- М., 1999. - Ч. 1. - С. 93-94.

48. Блягоз, Х.Р. Разработка процесса экстракции двуокисью углерода/Х.Р. Блягоз, A.A. Схаляхов, Е.П. Кошевой//Технология, химия и экологические проблемы Северного Кавказа: сб. науч. тр. по материалам регион, науч.-практ. конф. «Неделя науки МГТИ»/М-во образования Рос. Федерации, Майкоп, гос. технол. ин-т. - Майкоп, 2000. - С. 34-35.- Библиогр. в конце ст.

49. Проблемы в развитии экстракции двуокисью углерода/ Х.Р. Блягоз, Е.П. Кошевой, Х.Р. Сиюхов, A.A. Схаляхов, В.Ю. Чун-дышко// Тез. докл. Междунар. науч.-практ. конф. «Молодые ученые -пищевым и перерабатывающим отраслям АПК (технологические аспекты производства)», г. Москва, 18-19 декабря 2000 г. - М., 2000.-С. 173-174.

50. Рафинация растительных масел экстракцией двуокисью углерода/А.А. Схаляхов, Е.П. Кошевой, А.Б. Боровский, Х.Р. Блягоз// Прогрессивные пищевые технологии - третьему тысячелетию = Advanced food technologies - to the third millenium: междунар. науч. конф. [19-22 сент. 2000 г.: тез. докл./редкол.: Петрик A.A. (гл. ред.) и др.]. - Краснодар, 2000. - С. 172-173.

51. Совершенствование техники для экстракции двуокисью углерода/Е.П. Кошевой, A.A. Схаляхов, Х.Р. Блягоз, Х.Р. Сиюхов, В.Ю. Чундышко//Материалы междунар. науч.-практ. конф. «Продовольственная индустрия юга России. Экологически безопасные энергосберегающие технологии хранения и переработки сырья растительного и животного происхождения»/Краснодар. НИИ хранения и перераб. с.-х. продукции; под ред. Сизенко Е.И. - Краснодар, 2000. -Ч. 2.-С. 157-158.

52. Физическая рафинация растительных масел экстракцией двуокисью углерода в сверхкритическом состоянии/А.А. Схаляхов,

Е.П. Кошевой, А.Б. Боровский, Х.Р. Блягоз//Материалы междунар. науч.-практ. конф. «Продовольственная индустрия юга России. Экологически безопасные энергосберегающие технологии хранения и переработки сырья растительного и животного происхождения»/ Краснод. НИИ хранения и перераб. с.-х. продукции; под ред. Сизенко Е.И.-Краснодар, 2000.-Ч. 1.-С. 137-138.

53. Энергообеспечение экстракционной установки с использованием теплового насоса / Е.П. Кошевой, Х.Р. Блягоз, Х.Р. Сиюхов, А.А. Схаляхов, В.Ю. Чундышко // Прогрессивные пищевые технологии -третьему тысячелетию = Advanced food technologies - to the third millenium: междунар. науч. конф. [19-22 сент. 2000 г.: тез. докл./ редкол.: Петрик А.А. (гл. ред.) и др.]. - Краснодар, 2000. - С. 187-188.

54. Мембранная экстракция жидких смесей двуокисью углерода/ А.А. Схаляхов, Х.Р. Блягоз, Е.П. Кошевой, B.C. Косачев// Низкотемпературные и пищевые технологии в XXI веке: междунар. науч.-техн. конф., посвящ. 70-летию С.-Петерб. гос. ун-та низкотемператур. и пищевых технологий, 6-7 июня 2001 г./[редкол.: Малышев А.А. и др.].-СПб., 2001.-С. 170.

55. Схаляхов, А.А. Состояние и перспективы производства биодизеля/А.А. Схаляхов, Е.П. Кошевой, А.И. Хомяков//Тезисы IX Междунар. конф. молодых ученых «Пищевые технологии и биотехнологии». - Казань, 2008. - С. 216.

СХАЛЯХОВ Анзаур Адамович

ПРОЦЕССЫ ПЕРЕНОСА МАССЫ И ТЕПЛА СИСТЕМЫ «ЖИДКОСТЬ-ЖИДКОСТЬ» В АППАРАТАХ С ОБЪЕМОМ, СТРУКТУРИРОВАННОМ МЕМБРАННЫМИ ПЕРЕГОРОДКАМИ

АВТОРЕФЕРАТ

Подписано в печать 03.09.2010 г. Формат бумаги 60х847)6. Бумага ксероксная. Гарнитура Тайме. Усл. печ. л. 3,0. Заказ №055. Тираж 100 экз.

Издательство МГТУ 385000, г. Майкоп, ул. Первомайская, 191

ВВЕДЕНИЕ

ЛИТЕРАТУРНЫЙ ОБЗОР

Глава 1 Состояние и направления развития процессов переноса массы и тепла в аппаратах с объемом, структурированном мембранными перегородками

Раздел 1.1 Общая характеристика состояния мембранных технологий 8 Раздел 1.2 Переработка растительных масел и жиров с использованием мембранной технологии

Раздел 1.3 Гидравлика, массообмен и теплообмен в мембранных процессах

Раздел 1.4 Заключения по обзору, цель, задачи, концепция и структура работы

ЧАСТЬ 1 РАЗВИТИЕ ТЕОРИИ ПЕРЕНОСА МАССЫ И ТЕПЛА В АППАРАТАХ С ОБЪЕМОМ, СТРУКТУРИРОВАННОМ МЕМБРАННЫМИ ПЕРЕГОРОДКАМИ

Глава 2 Анализ переноса массы и тепла в аппарате при контакте через мембранную перегородку

Раздел 2.1 Общая характеристика процессов переноса массы и тепла между фазами, разделенными мембранными перегородками

Раздел 2.2 Анализ эффективности переноса массы и тепла в мембранном контакторе

Глава 3 Сопряженный конвективный перенос тепла и массы от внутреннего потока с теплопроводностью и диффузией через мембранную перегородку

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

Раздел 3.2 Постановка и решение сопряженной задачи конвективного переноса массы с диффузией в трубчатой мембране от внутреннего^ потока

Глава 4 Конвективный перенос тепла от внутреннего потока с теплопроводностью в непроницаемых половолоконных мембранах

Раздел 4.1 Моделирование и обоснование обобщенной кинетической зависимости переноса тепла в трубчатой мембране от внутреннего потока

Раздел 4.2 Перенос тепла в теплообменнике с непроницаемыми половолоконными мембранами

Глава 5 Перенос потоков в мембранах под действием трансмембранного перепада давления

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

Раздел 5.2 Моделирование переноса потоков в мембранных контакторах под действием трансмембранного перепада давления

ЧАСТЬ 2 МОДЕЛИРОВАНИЕ ЭФФЕКТИВНЫХ РЕЖИМОВ МАССО-И ТЕПЛОПЕРЕНОСА В АППАРАТАХ С ОБЪЕМОМ, СТРУКТУРИРОВАННОМ МЕМБРАННЫМИ ПЕРЕГОРОДКАМИ

Глава 6 Определение параметров массообмена процесса экстракции в мембранном контакторе

Раздел 6.1 Равновесные и массообменные свойства в системе «жидкость-двуокись углерода» при экстракционной очистке растительных масел

Раздел 6.2 Процесс переноса массы в мембранном контакторе

Раздел 6.3 Экспериментальные исследования процесса экстракции свободных жирных кислот из растительных масел двуокисью углерода

Раздел 6.4 Оценка эффективных режимов работы мембранного контактора

Глава 7 Определение параметров массопереноса реакции переэтерификации при производстве биотоплива в мембранном реакторе

Раздел 7.1 Общая структура мембранного реактора и математической 238 модели процесса получения биотоплива

Раздел 7.2 Математическое моделирование кинетики и свойств 243 реакционной смеси при производстве биодизельного топлива из масел

Раздел 7.3 Математическое моделирование мембранного реактора для производства биодизельного топлива из масел

Глава 8 Определение параметров процесса конденсации парогазовых смесей в конденсаторе с полимерными половолоконными мембранами 282 Раздел 8.1 Обоснование схемы конденсатора для парогазовых смесей с полимерными половолоконными мембранами

Раздел 8.2 Моделирование конденсации парогазовых смесей с полимерными половолоконными мембранами

ЧАСТЬ 3 ТЕХНИЧЕСКИЕ И ТЕХНОЛОГИЧЕСКИЕ РЕШЕНИЯ ПЕРЕРАБОТКИ РАСТИТЕЛЬНЫХ МАСЕЛ С ИСПОЛЬЗОВАНИЕМ СТРУКТУРИРОВАНИЯ ОБЪЕМА АППАРАТОВ МЕМБРАННЫМИ ПЕРЕГОРОДКАМИ

Глава 9 Разработка технологической схемы и линии комплексной переработки растительных масел с применением мембранных аппаратов

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

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

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

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

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

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

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

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

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

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

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

Диссертационная работа выполнена в Майкопском государственном технологическом университете в рамках подпрограммы: «Разработка инновационных технологий и техники переработки сельскохозяйственного сырья и производства продовольствия на основе мембранных и экстракционных процессов» республиканской целевой программы «Развитие сельского хозяйства и регулирование рынков сельскохозяйственной продукции, сырья и продовольствия», утвержденной постановлением Государственного Совета-Хасэ Республики Адыгея от 30 января 2008 года № 689-ГС.

ЛИТЕРАТУРНЫЙ ОБЗОР

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

2. Анализ эффективных режимов переноса массы и тепла в контакторе с половолоконными мембранами на основе полученных зависимостей концентрационных и температурных распределений по длине аппарата для различного взаимного направления потоков и продольного перемешивания во взаимодействующих фазах определил, что условием достижения предельной эффективности является противоточный процесс с числом единиц переноса по фазе х до 4-5, превышение числа единиц переноса по фазе * по сравнению с числом единиц переноса по фазе у, перемешивание по фазе х должно быть минимальным, перемешивание по фазе у существенно на эффективность не влияет.

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

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

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

6. Зависимости для поперечной и тупиковой фильтрации позволяют определить по результатам экспериментов определения проницаемости мембран, которые в исследованном диапазоне (избыточное давление до 1 бар) практически не зависели от давления и соответственно составили для полипропиленовых половолоконных мембран фирмы ЕГО08 (Чехия): РР-М5

- 4,4-Ю-15 м; РР-М6 - 1,33-1СГ14м и керамических трубчатых мембран,

11 произведенных фирмой НПО «Керамикфильтр» (Москва) - 7,43- Ю- м.

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

8. Для коэффициентов диффузии жирных кислот в двуокиси углерода в диапазоне изменения давления от 80 до 400 бар и температур от 293 до 373 К выделяются две области, разделенных линейной зависимостью, проходящей через критическую точку. Ниже этой линии значения коэффициентов диффузии практически не зависят от изменения значения температуры и значительно снижаются со снижением значения давления. Выше этой линии значения коэффициентов диффузии увеличиваются с ростом температуры и уменьшаются при снижении давления.

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

10. Вязкость реакционной смеси определяется вязкостью ее компонентов и составом. Температурная зависимость вязкости компонентов определяется по усовершенствованному методу Льюиса и Сквайрса.

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

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

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

14. Построенная математическая модель конденсации парогазовой смеси в конденсаторе с полимерными половолоконными мембранами может быть использована для моделирования и анализа работы конденсатора.

15. Обоснована комплексная технология и система процессов переработки нерафинированных растительных масел с новым процессом физической рафинации мембранной экстракцией и экстракционной очисткой фосфолипидов в схеме, в которой применен мембранный конденсатор парогазовой смеси, а отделенное масло перерабатывается в мембранном реакторе с получением биотоплива. Разработаны и переданы для реализации в производство филиала «Лабинский МЭЗ» ООО «МЭЗ Юг Руси» комплексная технология переработки нерафинированного подсолнечного масла.

316

1. Алтунин B.B. Теплофизические свойства двуокиси углерода. М.: Издательство стандартов, 1975. 546 с.

2. Альперт Л.З Основы проектирования химических установок. М.: Высш. шк., 1989. 304 с.

3. Арутюнян, Н. С., Корнена Е.П., Янова А.И. и др. Технология переработки жиров. М.: Пищепромиздат, 1998. 452 с.

4. Блягоз Х.Р., Кошевой Е.П., Схаляхов A.A. Коэффициенты диффузии триглицеридов и жирных кислот в двуокиси углерода при сверхкритических условиях//Известия ВУЗов. Пищевая технология. 2000.- № 2-3. - С. 62-63.-Библиогр. в конце ст.

5. Блягоз Х.Р., Кошевой Е.П., Схаляхов A.A. и др. Массоперенос в массообменнике с пористой перегородкой (мембраной)//Новые технологии: сб. науч. тр. МГТУ. / М-во образования и науки Рос. Федерации, Майкоп, гос. технол. ун-т. Майкоп, 2005. - С. 175-177.

6. Блягоз, Х.Р. Развитие научных основ и разработка процессов и техники экстракции двуокисью углерода в пищевой промышленности: автореферат диссертации на соискание ученой степени доктора технических наук/Блягоз Х.Р.-.-Краснодар, 2002.- 49 с.

7. Блягоз Х.Р., Схаляхов A.A., Кошевой Е.П. Массообмен при экстракции в аппаратах с пористой перегородкой/УИзвестия ВУЗов. Пищевая технология.- 2003.- №5-6.- С. 82-84.- Библиогр. в конце ст.

8. Блягоз Х.Р., Схаляхов A.A., Кошевой Е.П. Технология и техника экстрагирования двуокисью углерода при сверхкритических условиях. Экономика и техника, 2005, №1, С.51-54.

9. Броунштейн Б.И., Щеголев В.В. Гидродинамика, массо- и теплообмен в колонных аппаратах.- Л.:Химия, 1988.-336 с.

10. Бутина Е.А. Научно-практическое обоснование технологии и оценка потребительских свойств фосфолипидных биологически активных добавок. Автореферат диссертации на соискание ученой степени доктора технических наук, КубГТУ, Краснодар, 2003. 53 с.

11. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. М. 1972. 720 с.

12. Верещагин, А.Г., Кошевой Е.П., Корнена Е.П., Бутина Е.А., Герасименко Е.О. Обоснование системы процессов получения фосфатидного концентрата в технологии производства БАД Витол//Известия ВУЗов «Пищевая технология», 2007. №1. с.108-109.

13. Верещагин А.Г., Схаляхов A.A., Кошевой Е.П. Влияние продольного перемешивания во взаимодействующих фазах при десорбции в пленочном трубчатом противоточном аппарате/ТИзвестия ВУЗов. Пищевая технология.-2007.- № 2,- С. 71-73.-Библиогр. в конце ст.

14. Гухман A.A. Применение теории подобия к исследованию процессов тепло-массообмена. М.:Высш. шк.,1967.303 с.

15. Дёч Г., Руководство к практическому применению преобразования Лапласа и Z-преобразования. М.: Изд. «Наука» ГРФМЛ, 1971. 288 с.

16. Диткин В.А., Прудников А.П. Операционное исчисление. М.:Высшая школа, 1975. 407 с.

17. Дытнерский Ю.И. Баромембранные процессы. Теория и расчет. М.: Химия, 1986. 272 с.

18. Касаткин А.Г. Основные процессы и аппараты химической технологии. М.: Химия, 1973. 752 с.

19. Кафаров В.В. Основы массопередачи. Системы газ-жидкость, пар-жидкость, жидкость-жидкость. М.: Высш. шк., 1972. 496 с.

20. Кафаров В.В., Глебов М.Б. Математическое моделирование основных процессов химических производств. М.: Высш. шк., 1991. 400 с.

21. Кафаров В.В., Дорохов И.Н. Системный анализ процессов химической технологии. Основы стратегии. М.: Наука, 1976. 500 с.

22. Коган В.Б. Теоретические основы типовых процессов химической технологии. Л.: Химия, 1977. 592 с.

23. Корнена Е.П. Химический состав, строение и свойства фосфолипидов подсолнечного и соевого масла. Диссертация доктора технических наук.-Краснодар, 1986.- 272 с.+ Прил. 47 с.

24. Кошевой Е.П., Блягоз Х.Р., Схаляхов А.А и др. Оценка развития научного направления «Экстракция двуокисью углерода//Известия ВУЗов. Пищевая технология. — 1999. № 1. - С. 8-11.- Библиогр. в конце ст.

25. Кошевой Е.П., Блягоз Х.Р., Схаляхов А.А и др. Поведение коэффициентов диффузии жирных кислот в С02 при сверхкритических условиях//Известия ВУЗов. Пищевая технология. 1999. - № 2-3. - С. 72-73.-Библиогр. в конце ст.

26. Кошевой Е.П., Блягоз Х.Р., Схаляхов А.А и др. Универсальная установка для экстракции двуокисью углерода//Известия ВУЗов. Пищевая технология. — 1999. № 4. — С. 67-69.- Библиогр. в конце ст.

27. Кошевой Е.П., Блягоз Х.Р., Схаляхов А.А и др. Свидетельство РФ на полезную модель №16458. Универсальная установка для экстракции двуокисью углерода. БИ №1, 2001.

28. Кошевой Е.П., Косачев B.C., Блягоз Х.Р., Схаляхов A.A. Теоретический анализ экстракции в массообменнике с пористой перегородкой//Известия ВУЗов. Пищевая технология. — 2001. № 5-6. - С. 66-68.- Библиогр. в конце ст.

29. Кошевой Е.П., Попова С.А., Масликов В.А. Использование обобщенных переменных для корреляции экстракционных свойств растворителя. Известия вузов. «Пищевая технология». 1973. №6. с. 116-119.

30. Крюков A.M. Разработка мембранно-компрессорной установки для хранилищ сельскохозяйственной продукции в регулируемой газовой среде. Автореферат кандидатской диссертации. Кр-р, 2004.

31. Кутателадзе С.С. Теплопередача при конденсации в кипении. 2-е изд. М.: Машгиз, 1952.232с.

32. Кутателадзе С.С. Теплопередача и гидродинамическое сопротивление. Справочное пособие. М.:Энергоатомиздат, 1990.367 с.

33. Кэйс В.М. Конвективный тепло- и массообмен. М.:Энергия, 1971.448с.

34. Кэйс В.М., Лондон А. Л. Компактные теплообменники. М.:Энергия, 1967.224 с.

35. Лейси Р. Технологические процессы с применением мембран. М.: Изд-во «Мир», 1976.370 с.

36. Лонцин М., Мерсон Р. Основные процессы пищевых производств. М.: Легкая и пищевая промышленность, 1983.384 с.

37. Лыков А.В: Теория теплопроводности. М.:Высшая школа, 1967.-600с.

38. Маркман A.JI. Химия липидов. Вып. 1,Ташкент,Изд-во АН УзССР, 1963. 176с.

39. Маркман А.Л. Химия липидов. Вып.2, Ташкент, Изд-во«ФАН»УзССР, 1970. 223с.

40. Меретуков М.А. Разработка процесса экструзионной агломерации обезжиренного фосфатидного концентрата при подготовке к отгонке растворителя. Автореферат диссертации на соискание ученой степени кандидата технических наук, КубГТУ, Краснодар, 2005. 21 с.

41. Орвис В. EXCEL для ученых, инженеров и студентов. К.: Юниор, 1999. 528с.

42. ПавловК.Ф., РоманковП.Г., НосковА.А. Примеры и задачи по курсу процессов и аппаратов химической технологии./Под ред. чл.-корр. АН СССР П.Г. Романкова. 9-е изд.,перераб. и доп. Л.:Химия, 1981. 560с.

43. ПлисА.И., СливинаН.А. Mathcad 2000. Математический практикум для экономистов и инженеров. М.гФинансы и статистика, 2000.-656с.!

44. Протодьяконов И.О., Глинский В.А. Экспериментальные методы исследования гидродинамики двухфазных систем в инженерной химии. Л.: Изд. ЛУ, 1982. 196 с.

45. Протодьяконов И.О., Марцулевич H.A., Марков A.B. Явления переноса в процессах химической технологии.- Л.: Химия, 1981. 264 с.

46. Реометрия пищевого сырья и продуктов: Справочник /Под ред. Ю.А. Мачихина. М.: Агропромиздат. 1990. 271 с.

47. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей: Справочное пособие. Л.: Химия, 1982. 592 с.

48. Романков П.Г., Курочкина М.И. Гидромеханические процессы химической технологии. Л.: ЛО «Химия», 1982.- 288 с.

49. Романков П.Г., Фролов В.Ф. Массообменные процессы химической технологии. Л.: Химия, 1990. 384 с.

50. Рудобашта С.П. Массоперенос в системах с твердой фазой. М.: Химия, 1980. 248 с.

51. Самарский A.A. Введение в численные методы. М.: ГРФМЛ "Наука", 1982. 272 с.

52. СамарскийА.А. Теория разностных схем. М.:ГРФМЛ "Наука",1983.616с.

53. Справочник по теплообменникам: в 2-х т. Т.1/ Пер. с англ., под ред. Б.С. Петухова, В.К. Шикова.- М.: Энергоатомиздат, 1987.-560 с.

54. Схаляхов, A.A. Мембранная экстракция двуокисью углерода в пищевой технологии / Схаляхов A.A.; ред. журн. «Известия ВУЗов. Пищевая технология».- Краснодар, 2007.- 159 е.: ил.- Библиогр.: с. 157-159 (141 назв.).-Деп. в ВИНИТИ 12.02.07, № 128-В2007.

55. Схаляхов, А.А Математическое моделирование кинетики и свойств реакционной смеси при производстве биодизельного топлива из масел / Схаляхов А.А.//Новые технологии. 2009. № 3, с. 29-35.

56. Схаляхов, А.А Математическая модель гидравлики мембранного реактора в линии производства биодизеля/Схаляхов А.А.//Новые технологии. 2009. №3.

57. Схаляхов, A.A. Результаты моделирования работы мембранного реактора для производства биодизельного топлива из масел//Новые технологии 2009.-№ З.-С. 42-48.

58. Схаляхов, A.A. Оценка эффективных режимов работы массообменника с пористой перегородкой//Новые технологии 2009.- № 3.-С. 48-54.

59. Схаляхов, A.A. Постановка сопряженной задачи переноса тепла через трубчатую мембрану от протекающего в ней потока//Новые технологии — 2009.-№4.- С.53-56.

60. Схаляхов, A.A. Численное решение сопряженной задачи переноса тепла через трубчатую мембрану//Новые технологии 2009. - № 4.- С. 57-62.

61. Схаляхов A.A., Блягоз Х.Р., Кошевой Е.П. Производство биотоплива из масел и жиров. Майкоп, Изд-во МГТУ, 2008. 132 с.

62. Схаляхов A.A., Блягоз Х.Р., Кошевой Е.П. Организация процесса переэтерификации в мембранном реакторе. Новые технологии. 2008, вып.5, с.34-36.

63. Схаляхов A.A., Блягоз Х.Р., Верещагин А.Г., Кошевой Е.П. Конденсация парогазовых смесей с полимерными половолоконными мембранами. Новые технологии.- 2009.- Вып. 1.- С. 35-39.

64. Схаляхов A.A., Верещагин А.Г., Бутина Е.А. Обоснование системы процессов получения фосфолипидных продуктов и биодизеля при переработке отходов гидратации масла // Новые технологии 2009. - № 3.-С. 39-42.

65. Схаляхов А.А, Верещагин А.Г., Косачев B.C., Кошевой Е.П. Конденсатор для парогазовых смесей с полимерными половолоконными мембранами. Известия ВУЗов «Пищевая технология». 2009. №2-3, с.68-70.

66. Схаляхов А.А, Верещагин А.Г., Косачев B.C., Кошевой Е.П. Теплообменники с полимерными половолоконными мембранами. Известия ВУЗов «Пищевая технология». 2009. №2-3, с.79-81.

67. Схаляхов A.A., Верещагин А.Г., Косачев B.C., Кошевой Е.П. Разработка модели конденсации парогазовых смесей с полимерными половолоконными мембранами // Новые технологии.- 2009.- Вып. 1.- С. 3943.

68. Схаляхов А.А, Верещагин А.Г., Кошевой Е.П., Косачев B.C., Гукасян A.B. Конденсатор. Патент на полезную модель № 61401 //Бюллетень изобретений, 2007. №6.

69. Схаляхов А.А, Косачев B.C., Кошевой Е.П. Математическое моделирование процесса разделения жидких смесей в мембранном модуле с различной организацией потоков. Известия ВУЗов «Пищевая технология». 2009. №2-3, с.71-74.

70. Схаляхов А.А, Косачев B.C., Кошевой Е.П., Никонов Е.О. Определение проницаемости половолоконных и трубчатых мембран. Известия ВУЗов «Пищевая технология». 2009. №2-3, с.96-98.

71. Схаляхов A.A., Кошевой Е.П., Косачев B.C. Вязкость компонентов реакционной смеси при производстве биодизельного топлива из растительных масел. Известия ВУЗов «Пищевая технология». 2009. №1, 113115.

72. Схаляхов A.A., Кошевой Е.П., Хомяков А.И. Состояние и перспективы производства биодизеля. Сборник тезисов докладов IX Международной конференции молодых ученых «Пищевые технологии и биотехнологии», Казань, 2008.-е. 216.

73. Теория тепломассообмена /С.И. Исаев и др.; Под ред. А.И. Леонтьева. М.: Высш. шк., 1979. 495 с.

74. Трейбал Р. Жидкостная экстракция. Пер. с англ./Под ред. С.З. Кагана. М.: Химия, 1966: 724 с.

75. Уэйлес С. Фазовые равновесия в химической технологии: В 2-х ч. 4.2. Пер. с англ. М.: Мир, 1989. 360 с.

76. Филиппов Л.П. Методы расчета и прогнозирования свойств веществ. М.: Изд-во МГУ, 1988. 252 с.

77. Флетчер К. Численные методы на основе метода Галеркина. М.: Мир, 1988. 352 с.

78. Хванг С.-Т., Каммермейер К. Мембранные процессы разделения. М.: Химия, 1981. 464 с.

79. Шервуд Т., Пигфорд Р., Уилки Ч. Массопередача.-М.:Химия, 1982-696с.

80. Ши Д. Численные методы в задачах теплообмена. М.: Мир, 1988.-544с.

81. Шлихтинг Г. Теория пограничного слоя. М.: ГРФМЛ "Наука", 1969. 742с.

82. Abdullah, N.S., Das D.B. Modelling nutrient transport in hollow fibre membrane bioreactor for growing bone tissue with consideration of multi-component interactions. Chemical Engineering Science 2007, 62, 5821-5839.

83. Abdullah, N.S., Jones, D.R., Das, D.B. Nutrient transport in bioreactors for bone tissue growth: Why do hollow fibre membrane bioreactors work? Chemical Engineering Science 2009, 64, 109-125.

84. Adrian, T., Wendland, M, Hasse, H., Maurer, G. High-pressure multiphase behavior of ternary systems carbon dioxide-water-polar solvent: review and modeling with the Peng-Robinson equation of state. Journal of Supercritical Fluids 1998, 12, 185-221.

85. Afrane, G., Chimowitz, E.H.,. Experimental investigation of a new supercritical fluid-inorganic membrane separation process. Journal of Membrane Science 1996, 116, 293-299.

86. Ahmed, T., Semmens, M. J. Use of sealed end hollow fibres for bubbeless membrane aeration: Experimental studies. Journal of Membrane Science, 1992, 69, 1-10.

87. Albaraka Z., Trebouet D., Tuna M., Loureiro J. M., Burgard M. SFGP 2007 -An Analysis of the Mass Transfer in a Non-Dispersive Solvent Extraction Closed System. Int. J. of Chem. Reactor Eng., 2008, vol. 6, A13.

88. Alexander P.R., Callahan R.W. Liquid-liquid extraction and stripping of gold with microporous hollow fibers. J. Membrane Sci., 1987, 35, 57.

89. Alonso, A. I., Urtiaga, A. M., Irabien, A., Ortiz, M. I. Extraction of Cr(VI) with Aliquat 336 in hollow fiber contactors: mass transfer analysis and modeling. Chem. Engng Sci. 1994, 49, 901-909.

90. Apelblat, A., Katzir-Katchalsky, A., Silberberg, A. A mathematical analysis of capillary-tissue fluid exchange. Biorheology 1974, 11, 1-49.

91. Bailey, A. F. G., Barbe, A. M., Hogan, P. A., Johnson, R. A., Sheng, J. The effect of ultrafiltration on the subsequent concentration of grape juice by osmotic distillation. J. Membr. Sci. 2000, 164, 195.

92. Baker, R. W.; Cussler, E. L.; Eykamp, W.; Koros, W. J.; Riley, R. L.; Strathmann, H. Membrane separation systems; Noyes Data Corp.: Park Ridge, NJ, 1991.

93. Bandelier, P.; Deronzier, J. C.; Lauro, F. Plastic heat exchangers. Mater. Tech. (Paris) 1992, 9-10, 67.

94. Bao, L.; Liu, B.; Lipscomb, G. Entry mass transfer in axial flows through randomly packed fiber bundles. AIChE J. 1999, 45 (11), 2346-2356.

95. Bao, L.; Lipscomb, G. Well-developed mass transfer in axial flows through randomly packed fiber bundles with constant wall flux. Chemical Engineering Science 2002, 57, 125-132.

96. Barbe, A. M., Bartley, J. P., Jacobs, A. L., Jonhson, R. A. Retention of volatile organic flavour/fragrance components in the concentration of liquid foods by osmotic distillation. J. Membr. Sci. 1998, 145, 67.

97. Basheer, S., Mogi, K., Nakajima, M. Development of a novel hollow-fiber membrane reactor for the interesterification of triglycerides and fatty acids using modified lipase. Proc. Biochem., 1995, 30, 531-536.

98. Basu R., Prasad R., Sirkar K.K. Nondispersive membrane solvent back extraction phenol. A.I.Ch.E.J., 1990, 36 (3), 450-460.

99. Basu R., Sirkar K.K. Hollow fiber contained liquid membrane separation of citric acid. A.I.Ch.E.J., 1991, 37, 383.

100. Beck R.E., Schultz J.S. Hindered diffusion in microporous membranes. Science, 1970, 170, 1302.

101. Beck R.E., Schultz J.S. Hindrance of solute diffusion within membranes as measured with microporous membranes of known pore geometry. Biochim. Biophys. Acta. 1972, 255, 273.

102. Belfort, G., Membranes and Bioreactors: A Technical Challenge in Biotechnology, Biotechnology and Bioengineering, 1989,33, 1047-1066.

103. Bhattachayya, D.K., De B.K., Das R., Dutta B.K. Membrane degumming and dewaxing of rice bran oil and it's refining. Fett/Lipid 100 Nr. 9, S. 416-421 (1998).

104. Bigg, D. M.; Stichford, G. H.; Talbert, S. G. Application of polymeric materials for condensing heat exchangers. Polym. Eng. Sci. 1989, 29 (16), 1111.

105. Bocquet, S., Torres, A., Sanchez, J., Rios, G.M., Romero, J. Modeling the mass transfer in solvent-extraction processes with hollow fiber membranes. A.I.Ch.E. Journal 2005, 51 (4), 1067-1079.

106. Bothun, G., Knutson, B., Strobel, H., Nokes, S., Brignole, E., Diaz, S.,. Compressed solvents for the extraction of fermentation products within a hollow fiber membrane contactor. Journal of Supercritical Fluids 2003a, 25, 119-134.

107. Biodiesel standards (http://www.biofuelsystems.com/ specification.htm).

108. Bourouni, K.; Martin, R.; Tadrist, L.; Tadrist, H. Experimental investigation of evaporation performances of a desalination prototype using the aero—evapo—condensation technique. Desalination 1997, 114, 111.

109. Bourouni, K.; Martin; R.; Tadrist, L.; Tadrist, H. Modeling of heat and mass transfer in a horizontal-tube falling-film evaporator for water desalination. Desalination 1998, 116, 165.

110. Bourouni, K.; Martin, R.; Tadrist, L.; Chaibi, M. T. Heat transfer and evaporation in geothermal desalination units. Appl. Energy 1999, 64, 129.

111. Bourouni, K.; Deronzier, J: C.; Tadrist, L. Experimentation and modeling of an innovative desalination unit. Desalination 1999, 125, 147.

112. Breslau, B. R., Testa, A. J., Milnes, B. A. and Medjanis, G. Advances in hollow fiber ultrafiltration technology. Polym. Sci. Technol. 1980, 13, 109-127.

113. Brewster, M. E.; Chung, K.-Y.; Belfort, G. Dean vortices with wall flux in a curved channel membrane system. A new approach to membrane module design. J. Membr. Sci. 1993, 81, 127.

114. Brinkman, H. C. A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles. Appl. Sci. Res. 1947, Al, 27-34.

115. Brotherton, J.D., Chau, P.C. Modelling of axial-flow hollow fibre cell culture bioreactors. Biotechnology Progress, 1996, 12, 575-590.

116. Brouwers, H. J. H.; Van der Geld, C. M. M. Heat transfer, condensation and fog formation in crossflow plastic heat exchangers. Int. J. Heat Mass Transfer 1996, 39 (2), 391.

117. Brudi K., Dachmen N., Schmieder H. Partition Coefficients of Organic Substances in Two-Phase Mixtures of Water and Carbon Dioxide at Pressures of 8 to 30 MPa and Temperatures of 313 to 333 K; J. Supercrit. Fluids, 1996, 9, 146151.

118. Bruining, W. J. A general description of flows and pressures in hollow fiber membrane modules. Chem. Engng Sci. 1989, 44, 1441-1447.

119. Budich, M., Brunner, G. Supercritical fluid extraction of ethanol from aqueous solutions. Journal of Supercritical Fluids, 2003, 25, 45-55.

120. Burns, J. R.; Jachuck, R. J. J. Condensation studies using cross-corrugated polymer film compact heat exchangers. Appl. Therm. Eng. 2001, 21, 495.

121. Cakl, J., Hruba, M., Jirankova, H., Dolecek, P. Combined membrane processes : Comparing of dead -end and cross-flow modes. 17th International Congress of Chemical and Process Engineering 27-31 August 2006. Prague, Czech Republic.

122. Calabro, V., Curcio, S., Iorio, G., A theoretical analysis of transport phenomena in a hollow fiber membrane bioreactor with immobilized biocatalyst, J. Membrane Sci., 2002,206,217-241.

123. Calabro, V., Jiao, B. L., Drioli, E. Theoretical and experimental study on membrane distillation in the concentration of orange juice. Ind. Eng. Chem. Res. 1994,33, 1803.

124. Can, R. Processing of oilseed crops. 1989, In Oil Crops of the World: Their Breeding and Utilization. Ch. 11, eds G. Robbelen, R. K. Downey & A. Ashfi. McGraw-Hill, New York.

125. Cao P., Tremblay A.Y., Dube M.A., Morse K. Effect of Membrane Pore Size on the Performance of a Membrane Reactor for Biodiesel Production Ind. Eng. Chem. Res. 46, 2007, 52-58.

126. Catchpole, O., King, M. Measurement and correlation of binary diffusion coefficients in near critical fluids. Industrial & Engineering Chemistry Research 1994, 33, 1828-1837.

127. Ceynowa, J., Adamczak, P., Staniszewski, Kinetics of olive oil hydrolysis in an enzyme membrane reactor bond graph network model of reactor performance (parti), ActaBiotechnologica., 1997, 17, 161-176.

128. Ceynowa, J., Adamczak, P., Analysis of the bond graph network model of the membrane reactor for olive oil hydrolysis, Sep. Purific. Tech., 2001, 22-23, 443-449.

129. Chantong A., Massoth F.E. Restrictive diffusion in aluminas. A.I.Ch.E.J., 1983,29, 725-731.

130. Chao-Hong He. Infinite-dilution diffusion coefficients in supercritical and high-temperature liquid solvents. Fluid Phase Equilibria, 1998, 147, 309-317.

131. Chao-Hong He, Yong-Sheng Yu, Estimation of Infinite-Dilution Diffusion Coefficients in Supercritical Fluids. Ind. Eng. Chem. Res., 1997, 36, 4430-4433.

132. Chao-Hong He; Yong-Sheng Yu; Wei-Ke Su. Tracer diffusion coefficients of solutes in supercritical solvents. Fluid Phase Equilibria, 1998, 142, 281-286.

133. Chen S.-H. A Rough-Hard-Sphere Theory for Diffusion in Supercritical Carbon Dioxide. Chem. Eng. Sci., 1983, 38, 655-660.

134. Chiu, Y., Tan, C. Regeneration of supercritical carbon dioxide by membrane at near critical conditions. Journal of Supercritical Fluids, 2001, 21, 81-89.

135. Chrestil, J., of solids and liquids in supercritical gases, J. Phys. Chem., 1982, vol.86, pp.3016 3021.

136. Coenen, H., Eggers, R., Kriegel, E., Die Trennung von Stoffgemischen durch Extraktion mit uberkritischen Gasen. Tech. Mitt. Krupp Forschungsber. 1982, vol.40,№l, pp. 1-11.

137. Colburn, A.P.; Edison, A.G.: Prevention of fog in cooler-condensers. Ind. Eng. Chem. 1941, 33 (4), 457-58.

138. Cooney, D., Jackson, C. Gas absorption in a hollow fibre device. Chem. Eng. Comm., 1989, 79, 153-163.

139. Costello M.J., Fane A.G., Hogan P.A., Schofield R.W. The effect of shell side hydrodynamics on the performance of axial flow hollow fiber modules. J. Membr. Sei., 1993, 80, 1-11.

140. Crowder, R., Cussler, E. Mass transfer resistances in hollow fiber pervaporation. Journal of Membrane Science 1998, 145, 173-184.

141. Cuperus, F. P., Nijhuis, H. H. Applications of membrane technology to food processing. Trends Food Sei. Technol, 1993, 4, 277-282.

142. Dahuron L., Cussler E.L. Protein extraction with hollow-fiber. A.I.Ch.E.J., 1988, 34, 130-136.

143. Darcovich K., Dal-Cin, M.M., Gros B. Membrane mass transport modeling with the periodic boundary condition. Computers and Chemical Engineering, 2009, 33, 213-224.

144. Das, D.B., Nassehi, V. Modelling of contaminants mobility in underground domains with multiple free/porous interfaces. Water Resources Research, 2003, 39 (3), 1072.

145. Das, D.B., Nassehi, V., Wakeman, R.J.,. A finite volume model for the hydrodynamics of combined free and porous flow in sub-surface regions. Advances in Environmental Research, 2002, 7 (1), 35-58.

146. Deen W.M. Hindered transport of large molecules in liquid-filled pores. A.I.Ch.E.J., 1987, 33, 1409.

147. D'Elia N.A., Dahuron L., Cussler E.L. Liquid-liquid extractions with microporous hollow fibers. J. Membrane Sci., 1986, 29, 309.

148. De Loach, J. R., Sheffield, C. L. and Spates, G E. A continuous-flow high-yield process for preparation of lipid-free hemoglobin. Anal. Biochem. 1986, 157, 191-198.

149. Ding H.B., Cussler E.L. Fractional extraction with hollow fibers with hydrogel-filled walls. A.I.Ch.E.Journal, 1991, 37, 855.

150. Ding H.B., Carr P.W., Cussler E.L. Racemic leucine separation by hollow-fiber extraction. A.I.Ch.E.Journal, 1992, 38, 1493.

151. Dohnal M., Foller B., Millionnyy A. Plastic capillary heat exchangers 18th International Congress of Chemical and Process Engineering 24-28 August 2008. Prague, Czech Republic.

152. Drioli E. Macedonio F., Curcio E., Di Profio G. Membrane-Based Systems for Seawater Desalination: Analysis and Comparison. 17th International Congress of Chemical and Process Engineering 27-31 August 2006. Prague, Czech Republic.

153. Drioli E., Romano M. Progress and New Perspectives on Integrated Membrane Operations for Sustainable Industrial Growth. Ind. Eng. Chem. Res. 2001, 40, 1277-1300.

154. Dymond J.H. Corrected Enskog Theory and the Transport Coefficients of Liquids. J. Chem. Phys., 1974, 60, 969-973.

155. El-Desssouky, H. T.; Ettouney, H. M. Plastic compact heat exchangers for single-effect desalination systems. Desalination 1999, 122, 271.

156. Eggers R.; Sievers U.; Stein W. High pressure extraction of oil seed. J. Am. Oil Chem. Soc., 1985, 62, 1228-1230.

157. Encinar J.M., Gonzalez J.F., Sabio E., Ramiro M.J. Preparation and Properties ofBiodisel from Cynara cardunculus L. Oil Ind. Eng. Chem. Res. 1999, 38,2927-2931.

158. Estay H., Bocquet S., Romero J., Sanchez J., Rios G.M., Valenzuela F. Modeling and simulation of mass transfer in near-critical extraction using a hollow fiber membrane contactor. Chemical Engineering Science 2007, 62, 5794-5808.

159. Ettouney, H. M.; El-Dessouky, H. T.; Alatiqi, I. Understand thermal desalination. Chem. Eng. Prog. 1999, 95 (9), 43.

160. Freedman, B., Butterfield, R.O., Pryde, E.H. Transesterification kinetics of soybean oil. JAOCS, 1986, 63, 1375-1380.

161. Frenander, U.; Jonsson, A. S. Cell harvesting by cross-flow microfiltration using shear-enhanced module. Biotechnol. Bioeng. 1996, 52, 397.

162. Friedrich J.P. Supercritical C02 extraction of lipids from lipid-containing materials. U.S. Patent 4466923, Apr. 1, 1982.

163. Funazukuri, T., Ishiwata, Y., Wakao, N. Predictive correlation for binary diffusion coefficients in dense carbon dioxide. A.I.Ch.E. Journal 1992, 38, 17611768.

164. Futselaar, H.; Reith, T.; Racz, I. G. The countercurrent transverse flow hollow fiber membrane module for the separation of liquid streams. Engineering of Membrane Processes, Gramish-Partenkirchent, Bavaria, Germany, May 1992.

165. Gabelman, A., Hwang, S. Hollow fiber membrane contactors. Journal of Membrane Science 1999, 159, 61-106.

166. Gabelman, A., Hwang, S. Experimental results versus model predictions for dense gas extraction using a hollow fiber membrane contactor. Journal of Supercritical Fluids 2005, 35, 26-39.

167. Gabelman, A., Hwang, S. A theoretical study of dense gas extraction using a hollow fiber membrane contactor. Journal of Supercritical Fluids, 2006, 157-172.

168. Gabelman, A., Hwang, S., Krantz, W. Dense gas extraction using a hollow fiber membrane contactor: experimental results versus model predictions. Journal of Membrane Science 2005, 257, 11-36.

169. Garcha A., Blvarez S., Riera F.A., Blvarez R., Coca J. Sunflower oil miscella degumming with polyethersulfone membranes effect of process conditions andjL

170. MWCO on fluxes and rejections. 16 International Congress of Chemical and Process Engineering 22-26 August 2004. Prague, Czech Republic.

171. Gawronski, R., Wrzesinska, B. Kinetics of solvent extraction in hollow fiber contactors. Journal of Membrane Science 2000, 168, 213-222.

172. Ghosh, R.; Cui, Z. F. Fractionation of BSA and lysozyme using ultrafiltration—Effect of gas sparging. AIChE J. 1998, 44, 61.

173. Giorno, L., Molinari, R., Natoli, M., Drioli, E. Hydrolysis and regioselective transesterification catalised by immobilized lipases in membrane bioreactors. J. Membr. Sci. 1997, 125, 177.

174. Giorno, L., Molinari, R., Drioli, E., Bianchi, D., Cesti, P. Performance of biphasic organic/aqueous hollow fibre reactor using immobilised lipase. J. Chem. Technol. Biotechnol. 1995, 64, 345.

175. Githens, R. E.; Minor, R. W.; Tomsic, V. J. Flexible tube heat exchangers. Chem. Eng. Prog. 1965, 61 (7), 55.

176. Godongwana, B., Sheldon, M.S., Solomons, D.M., Momentum transfer inside a vertically orientated capillary membrane bioreactor, J. Membrane Sci., 2007, 303, 86-99.

177. Grzywna Z.J. Scaling in diffusive transport through membranes. Chem. Eng. Sci., 1996, 51, No. 17, 4115-4125.

178. Guo J.-J., Ho C.-D. Theoretical and experimental studies of membrane extraction74.of Cu with D2EHPA through rectangular conduits. Chemical Engineering and Processing 2009, 48, 111-119.

179. Hall D.W., Scott K., Jachuck R.J.J., Determination of mass transfer coefficient of a cross-corrugated membrane reactor by the limiting-current technique. International Journal of Heat and Mass Transfer 2001,44,2201-2207.

180. Handbook of industrial membrane technology. Edited by Mark C. Porter. Reprint Edition. NOYES PUBLICATIONS, Westwood, New Jersey, U.S.A. 1988. 660 p.

181. Happel J. Viscous flow relative to arrays of cylinders. A.I.Ch.E. J., 1959, 5, 174-177.

182. Haraldsson, G. Degumming, dewaxing and refining. J Am. Oil Chem. Soc., 1983, 60(2), 203A-208A.

183. Hetsroni, G.; Mosyak, A. Heat transfer and pressure drop in a plastic heat exchanger with triangular channels. Chem. Eng. Process. 1994, 33, 91.

184. Hickman, H. J. An asymptotic study of the Nusselt-Graetz problem. Part I: Large x behavior. J. Heat Transfer 1974, 96, 354.

185. Hirohama, S., Takatsuka, T., Miyamoto, S., Muto, T. Measurement and correlation of phase equilibria for the carbon dioxide-ethanol-water system. Journal of Chemical Engineering of Japan, 1993, 26, 408-415.

186. Hogan, P. A., Canning, R. P., Peterson, P. A., Johnson, R. A., Michaelis, A. S. A New Option: Osmotic Distillation. Chem. Eng. Prog. 1998, 94, 49.*'

187. Hoq, M. M., Yamane, T., Shimizu, S. Continuous synthesis of glycerides by lipase in a microporous membrane bioreactor. J. Am. Oil Chem. Soc., 1984,61,776-781.

188. Hoq, M. M., Yamane, T., Shimizu, S., Funada, T., Ishida, S. Continuous hydrolysis of olive oil by lipase in microporous hydrophobic membrane bioreactor. J. Am. Oil Chem. Soc., 1985, 62, 1016-1021.

189. Iritani, E., Murakami. Y., Murase, T. Analysis of flow of power-law non-Newtonian fluids in hollow fiber membrane modules. Chem. Engng Sci. 1994; 49, 917-919.

190. Iwama, A. New process for purifying soybean oil by membrane separation and economic evaluation of the process. J.Am. Oil Chem. Soc., 1987, 64(9), 1258.

191. Jaakkola, H. Cost-effective evaporators for desalination. Desalination 1996,108,357.

192. Jachuck, R. J. J.; Ramshaw, C. Process intensification: Polymer film compact heat exchanger (PFCHE). Chem. Eng. Res. Des. 1994, 72A, 255.

193. Kaufmann S., Loretz Y., Hilfiker K. Prevention of fog in a condenser by simultaneous heating and cooling. Heat and Mass Transfer 1997, 32, 403-410.

194. Kaufman, J., Ruebusch, R. J. Oleochemicals a look at world trends. INFORM, 1990, 1, 1034-1048.

195. Kays, W. M.; London, A. L. Compact Heat Exchangers, 3rd ed.; McGraw-Hill: New York, 1984.

196. Kelsey, L. J., Pillarella, M. R., Zydney, A. L. Theoretical analysis of convective flow profiles in a hollow-fiber membrane bioreactor. Chem. Engnq Sei. 1990, 45,3211-3220.

197. Kertesz R., Schlosser S. Simulation of simultaneous membrane based solvent extraction and stripping of phenylalanine in hollow fiber contactors. 16th Int. Congress of Chem. and Proc. Eng. 22-26 August 2004. Prague, Czech Republic.

198. Keurentjes, J. T. F., Sluijs. J. T. M., Franssen, R. J. H., van't Riet, K. Extraction and fractionation of fatty acids from oil using an ultrafiltration membrane. Ind. Engng Chem. Res. 1992, 31, 581-587.

199. Kiani A., Bhave R.R., Sirkar K.K. Solvent extraction with immobilized interfaces in microporous membrane. J. Membrane Sei., 1984, 20, 125.

200. Kibboua R., Kechnit A, Azzi A. Laminar film condensation on an elliptical tube with porous coating //17th International congress of chemical and process engineering 27-31 August 2006 Prague-CZECH Republic.

201. Kim I.-C., Kim J.-H., Lee K.-H., Tak T.-M. Phospholipids separation (degumming) from crude vegetable oil by polyimide ultrafiltration membrane. J. Membr. Sei., 2002,205, 113-123.

202. Kimura, Y., A. Tanaka, K. Sonomoto, T. Nihira, and S. Fukui, Application of Immobilized Lipase to Hydrolysis of Triacylglyceride, Eur. J. Appl. Microbiol. Biotechnol. 1983, 17, 107-112.

203. Koike S., Subramanian R., Nabetani H., Nakajima M. Separation of oil constituents in organic solvents using polymeric membranes, JAOCS, 2002, 79, 937-942.

204. Koonaphapdeelert S., Wu Z., Li K. Carbon dioxide stripping in ceramic hollow fibre membrane contactors. Chemical Engineering Science, 2009, 64, 1-8.

205. Koris, A., Vatai G. Membrane filtration for vegetable oil degumming. 15th International Congress of Chemical and Process Engineering 25 29 August 2002 Praha, Czech Republic.

206. Koseoglu, S. S. Membrane technology for edible oil refining. Oils Fats Int., 1991,5, 16-21.

207. Koseoglu, S.S., Engelgau, D. E. Membrane applications and research in the edible oil industry: An assessment. J.Am. Oil Chem. Soc., 1990, 67(4), 239249.

208. Krogh, A. The number and distribution of capillaries in muscles with calculations of the oxygen pressure head necessary for supplying the tissue. J. Physiol. 1919, 52,409-415.

209. Kubaczka A., Burghardt A., Mokrosz T. Membrane-based solvent extraction s multicomponent systems. Chem. Eng. Sci., 1998, 53, №5, 899-917.

210. Kubaczka A., Burghardt A. Effect of mass transport resistances in multicomponent membrane extraction on the overall mass fluxes Chem. Eng. Sci. 2000, 55,2907-2916.

211. Labecki M., Piret J.M., Bowen B.D. Two-dimensional analysis of fluid flow in hollow-fibre modules. Chemical Engineering Science, 1995, Vol. 50, No. 21,3369-3384.

212. Labecki, M., Bowen, B.D., Piret, J.M.,. Two-dimensional analysis of protein transport in the extracapillary space of hollow-fibre bioreactors. Chemical Engineering Science, 1996, 51 (17), 4197-4213.

213. Laborie, S.; Cabassud, C.; Durand-Bourlier, L.; Laine, J. M. Flux enhancement by a continuous tangential gas flow in ultrafiltration hollow fibres for drinking water production: effects of slug flow on cake structure. Filtr. Sep. 1997, 34, 887.

214. Lauro, P. Projet d'usine de dessalement a multiple-effet a surface d1 échangé en matiere plastique. Desalination 1979, 31, 221.

215. Li, K., Tai, M. S. L., Teo, W. K. Design of a C02 scrubber for self-contained breathing systems using a microporous membrane. Journal of Membrane Science, 1994, 86, 119-125.

216. Li K., Kong Jianfeng, Tan Xiaoyao Design of hollow fibre membrane modules for soluble gas removal. Chem. Eng. Sci. 2000, 55, 5579-5588.

217. Liang T.-T., Long R.L., Corrections to Correlations for Shell-Side Mass-Transfer Coefficients in the Hollow-Fiber Membrane (HFM) Modules. Ind. Eng. Chem. Res. 2005, 44, 7835-7843.

218. Lin L., Rhee K.C., Koseoglu S.S. Bench-scale membrane degumming of crude vegetable oil: Process optimization, J. Membr. Sci., 1997, 134, 101-118.

219. Liu, W.; Davidson, J.; Mantell, S. Thermal analysis of polymer heat exchangers for solar water heating: A case study. J. Sol. Energy Eng. 2000, 122, 84.

220. Loeb S., Sourirajan S., Advan; Chem. Ser., 1963, 38, 117.

221. London, A L. Compact heat exchangers—design methodology. In Low Reynolds Number Flow Heat Exchangers; Kakac, S., Shah, R. K., Bergles, A. E., Eds.; Hemisphere Publishing Co.: New York, 1983; pp 815-844.

222. Luck, E., Marr, R., Estimation of The Process Parameter For High-Pressure Carbon Dioxide Extraction of Nature Products, Separ. Sci. Technology,1988, vol.23, № 1-3, pp.63-76.

223. Macroric, C.; Freoman, S. Design and operation of membrane filtration plants for water treatment. World Filtration Congress, Brighton, U.K., April 2000. In Proceedings, Vol. 1, pp 525-528.

224. Malcata, F.X., Hill, C.G., Amundson, C.H., Use of a lipase immobilized membrane in a membrane reactor to hydrolyze the glycerides of butteroil, Biotecnol. Bioeng., 1991, 38, 853-868.

225. Malcata, F.X., Reyes, H.R., Garcia, H.S., Hill, C.G., Amundson, C.H., Kinetics and mechanisms of reactions catalyzed by immobilized lipases, Enzyme Microb. Technol., 1992,14,426-446.

226. Malcata, F.X., Hill, C.G., Amundson, C.H., Hydrolysis of butteroil by immobilized lipase using a hollow-fiber reactor: Part III. Multiresponse kinetic studies, Biotecnol. Bioeng., 1992, 39, 1002-1012.

227. Manno, P.; Moulin, P.; Rouch, J. C.; Clifton, M.; Aptel, P. Mass transfer improvement in helically wound hollow fibre ultrafiltration modules. Yeast suspensions Sep. Purif. Technol. 1998, 14, 175.

228. McHugh M.A., Krukonis V.J. Supercritical Fluid Extraction: Principles and Practice. Published by Butterworth-Heinemann, 1994.

229. Meldrum, A., Hollow fibre membrane reactors, The Chemical Engineer 1987, 441,28-31.

230. Michaels, A. S. In Separation for Biotechnology 2; Pyle, D. L., Ed.; Elsevier Applied Science: Cambridge, U.K, 1990; p 3.

231. Michaels, S. L. Crossflow microfilters: the ins and outs. Chem. Engng1989, 96, 84-91.

232. Mikulasek, P., Hrdy, J. Permeate flux enhancement using a fluidized bed in microfiltration with ceramic membranes. Chem. Biochem. Eng. 1999, 13 (3), 133.

233. Miyatake O., Iwashita H. Laminar flow heat transfer to a fluid flowing axially between cylinders with a uniform wall heat flux. Int. J. Heat Mass Transfer, 1991,34, 322-327.

234. Moon, W. S.; Park, S. B. Design guide of a membrane for a membrane reactor in terms of permeability and selectivity. J. Membr. Sci. 2000, 170,43-51.

235. Myerson, A. S. Handbook of Industrial Crystallization, 2nd ed.; Butterworth-Heinemann: Boston, MA, 2002.xL

236. Nagy E. Mass transport through biocatalytic membrane reactors. 18 Int. Congress of Chem. and Proc. Eng. 24-28 August 2008. Prague, Czech Republic.

237. Nassehi, V. Modelling of combined Navier-Stokes and Darcy flows in crossflow membrane filtration. Chemical Engineering Science, 1998, 53 (6), 1253-1265.

238. Nunge, R.J., Gill, W.N. An analytical study of laminar counterflow double-pipe heat exchangers, AIChE J. 1966, 12, 279-289.

239. Nwuha V. Novel studies on membrane extraction of bioactive components of green tea in organic solvents: part I. J. Food Engineering, 2000, 44, 233-238.

240. Ochoa N., Pagliero C., Marchese J., Mattea M. Ultrafiltration of vegetable oils. Degumming by polymeric membranes. Sep. Purif. Technol., 2001, 22-23, 417-422.

241. Ohkuma, N., Shinoda, T., Aoi, T., Okaniwa, Y., Magara, Y. Performance of rotary disk modules in a collected human excreta treatment plant. Water Sci. Technol. 1994, Vol. 30, No. 4,141.

242. Padley, F. B., Gunstone, F. D. & Harwood, J. L. Occurrence and characteristics of oils and fats. In The Lipid Handbook. 2nd edn. Ch. 3. Chapman and Hall, London, 1994.

243. Pagliero C., Ochoa N., Márchese J., Mattea M. Degumming of crude soybean oil by ultrafiltration using polymeric membranes, JAOCS, 2001, 78, 793796.

244. Pagliero C., Ochoa N., Márchese J., Mattea M., Vegetable oil degumming with polyimide and polyvinilidenfluoride ultrafiltration membranes, J. Chem. Technol. Biotechnol., 2004, 79, 148-152.

245. Pangrle, B. J., Alexandrou, A. N., Dixon, A. G., DiBiasio, D. An analysis of laminar fluid flow in porous tube and shell systems. Chem. Engng Sci. 1991,46, 2847-2855.

246. Pangrle, B. J., Walsh, E. G., Moore, S., DiBiasio, D. Magnetic resonance imaging of laminar flow in porous tube and shell systems. Chem. Engng Sci. 1992, 47, 517-526.

247. Patel, A. B.; Brisson, J. G. Design, construction and performance of plastic heat exchangers for sub-Kelvin use. Cryogenics 2000, 40, 91.

248. Patil V. E. Membrane Technology for the Regeneration of Supercritical Carbon Dioxide. PROEFSCHRIFT, 2006.

249. Patkar, A. Y., Koska, J., Taylor, D. G., Bowen B. D., Piret, J. M. Protein transport in ultrafiltration hollow-fiber bioreactors. A.I.Ch.E. J. 1995, 41, 415-425.

250. Peng D.-Y., Robinson D. B., A new two-constant equation of state. Ind. Eng. Chem. Fundament., 1976, vol.15, pp.59 64.

251. Pierre, E, Souchon, I., Martin, M. Recovery of sulfur aroma compounds using membrane-based solvent extraction. Journal of Membrane Science, 2001, 187, 239-253.

252. Prasad R., Kiani A., Bhave R.R., Sirkar K.K. Further studies on solvent extraction with immobilized interfaces in microporous membrane. J. Membrane Sci., 1986, 26, 79.

253. Prasad R., Sirkar K.K. Dispersion-free solvent extraction with microporous hollow-fiber modules. A.I.Ch.E.J., 1988; 34, 177-188.

254. Prasad, R., Sirkar, K. K. Hollow fiber solvent extraction of pharmaceutical products: A case study, Journal of Membrane Science, 1989, Vol. 47,235-259.

255. Prasad R., Sirkar K.K. Hollow Fiber Solvent Extraction: Performances and Design. J. Memb. Sci., 1990, 50, 153-175.

256. Pronk, W., van der Burgt, M., van't Boswinkel, G., Riet, K. A hybrid membrane-emulsion reactor for the enzymatic hydrolysis of lipids. J. Am. Oil Chem. Soc., 1991, 68, 852-856.

257. Pryde, E. H., Rothfus, J. A. Industrial and non food use of vegetable oils. In Oil Crops of the World: Their Breeding and Utilization. Ch. 5. eds. G., 1989.

258. Pugazhenthi, G., Kumar, A., Enzyme membrane reactor for hydrolysis of olive oil using lipase immobilized on modified PMMA composite membrane, J. Membr. Sci., 2004, 228, 187-197.

259. Pugazhenthi, G., Kumar, A., Modeling of the Mass Transfer Effect in Biphaac Enzyme Membrane Reactor for Hydrolysis of Olive Oil. Int. J.ofFoodEng., 2008, Vol.4, № 5, A8.

260. Raman, R.; Mantell, S.; Davidson, J.; Wu, C.; Jorgensen, G. A review of polymer materials for solar water heating systems. J. Sol. Energy Eng. 2000, 122, 92.

261. Raman, L. P., Rajagopalan, N., Cheryan, M. Membrane technology. Oils Fats Int., 1994, 6(10), 28-36.

262. Raman, L. P., Cheryan, M., Rajagopalan, N. Deacidification of soybean oil by membrane technology. J. Am. Oil Chem. Soc., 1996, 73(2), 219224.

263. Ramirez, J. A.; Davis, R. H. Application of cross-flow microfiltration with rapid backpulsing to wastewater treatment. J. Hazard. Mater. 1998, 63, 179.

264. Reddy K., Subramanian R., Kawakatsu T., Nakajima M., Decolorization of vegetable oils by membrane processing, Eur. Food Res. Technoh, 2002, 213, 212-218.

265. Renkin E.M. Filtration, diffusion and molecular sieving through porous cellulose membranes. J. Gen. Physiol. 1954, 38, 225.

266. Robbelen, R. K. Downey & A. Ashri. McGraw-Hill, Publishing Company, New York.

267. Robinson J.R., Sims M., Method and System for Extracting a Solute from a Fluid Using Dense Gas and a Porous Membrane. U.S. Patent 5490884; Feb. 13, 1996.

268. Rossell, J. B. Vegetable oils and fats. 1991. In Analysis of Oilseeds, Fats and Fatty Foods. Ch. 7, eds. J. B. Rossell & J. L. R. Pritchard. Elsevier Applied Science, Amsterdam.

269. Rousse, D. R.; Martin, D. Y.; Theriault, R.; Levillee, F.; Boily, R. Heat recovery in greenhouses: A practical solution. Appl. Therm. Eng. 2000, 20, 687.

270. Sandvig, H. J. Total solvent plant safety. J.Am. Oil Chem. Soc., 1983, 60,243.

271. Sarrade, S., Rios, G.M., Carles, M. Nanofiltration membrane behavior in supercritical medium. Journal of Membrane Science 1996, 114, 81-91.

272. Sarrade, S., Guizard, C., Rios, G.M. Membrane technology and supercritical fluids: chemical engineering for coupled processes. Desalination 2002, 144, 137-142.

273. Satone, H. Comparison between MSF distillation and RO. Technol. Proc. 9th Annu. Conf. NWSIA 1981, Vol. I, Session II.

274. Schonberg, J.A., Belfort, G., Enhanced nutrient transport in hollow fiber perfusion bioreactors: a theoretical analysis, Biotechnology Progress, 1987, 3, 81-89.

275. Schoner P., Plucinski P., Nitsch W., Daiminger U. Mass transfer in the shell side of cross flow hollow fiber modules. Chemical Engineering Science, 1998, 53,2319-2326.

276. Schultz W.G. Process for extraction of flavors. U.S. Patent 3477856, Nov. 10,1965.

277. Schultze С., Donohue M.D. Prediction of Henryks constants for supercritical fluids using a van der Waals equation of state. Fluid Phase Equilibria, 1998, 142, 101-114.

278. Sen Gupta, A.K. Neuere Entwicklungen auf dem Gebiet der Raffination der Speiseole. Fette Seifen Anstrichm., 1986, 88(3), 79-86.

279. Seo G., Massoth F.E. Effect of pressure and temperature on restrictive diffusion of solutes in aluminas. A.I.Ch.E.J., 1985, 31, 494-496.

280. Sims, M. Porocritical fluid extraction from liquids using near-critical fluids. Membrane Technology 1998, 97, 11-12.

281. Sims M., McGovern W. E., Robinson J. R. Porocritical fluid extraction application: continuous pilot extraction of natural products from liquids with near critical fluids. PoroCrit LLC. (Информация из Интернета).

282. Sims M., Robinson J.R., Dennis AJ. Paper presented at the American Chemical Society National Meeting, New Orleans, March 24-26, 1996.

283. Sirkar К. K., Shanbhag P. V., Kovvali A. S. Membrane in a Reactor: A Functional Perspective.Ind. Eng. Chem. Res. 1999, 38, 3715-3737.

284. Snape J. В., Nakajima M. Processing of Agricultural Fats and Oils using Membrane Technology. Journal of Food Engineering 1996, 30, 1-41.

285. Stankiewicz, A. L; Moulijn, J. A. Process Intensification: Transforming Chemical Engineering. Chem. Eng. Prog. 2000, 96, 22.

286. Suarez J.J.; Medina I.; Bueno J.I. Diffusion coefficients in supercritical fluids: available data and graphical correlations. Fluid Phase Equilibria, 1998, 153, 167-212.

287. Subramanian R., Nakajima M. Membrane degumming of crude soybean and rapeseed oils. JAOCS, 1997, vol. 74, no 8.

288. Subramanian R., Nakajima M., Kawakatsuf T. Processing of Vegetable Oils Using Polymeric Composite Membranes. Journal of Food Engineering 1998, 38, 4156.

289. Subramanian R., Raghavarao K.S.M.S., Nabetani H., Nakajima M.5 Kimura T., Maekawa T., Differential permeation of oil constituents in nonporous denser polymeric membranes, J. Membr. Sci., 2001, 57-69.

290. Subramanian R., Raghavarao K.S.M.S., Nakajima M., Nabetani H., Yamaguchi T., Kimura T. Application of dense membrane theory for differential permeation of vegetable oil constituents, J. Food Eng., 2003, 60, 249-256.

291. Sun, S., Koseoglu, S. S. Membrane degumming of cottonseed miscella. Am. Oil Chem. Soc. Annual Meeting, Atlanta, GA, INFORM, 1994, 5(4), 481.

292. Takeuchi, H., Takahashi, K., Nakano, M. Mass transfer in single oil-containing microporous hollow fiber contactors. Industrial & Engineering Chemistry Research 1990, 29, 1471-1476.

293. Taylor, D.G., Piret, J.M., Bowen, B.D. Protein polarization in isotropic membrane hollow fibre bioreactors. A.I.Ch.E. Journal 1994, 40 (2), 321-333.

294. Ternan M. The diffusion of liquid in pores. Can. J. Chem. Engng., 1987, 65, 244.

295. Traub, P., Stephan, K. High-pressure phase equilibria of the system C02-water-acetone measured with a new apparatus. Chemical Engineering Science, 1990, 45, 751-758.

296. Vaidya, A. M., Bell, G., Hailing, P. J. Surfactant-induced breakthrough effects during the operation of two-phase biocatalytic membrane reactors. Biotechnol. Bioengng., 1994, 44, 765-771.

297. Vaidya, A. M., Hailing, P. J., Bell, G. Aqueous-organic membrane bioreactors. Part 1. A guide to membrane selection. J. Membrane Sci., 1992, 71, 139-149.

298. Van der Padt, A., Edema, M. J., Sewalt, J. J. W., van't Riet, K. Enzymatic acylglycerol synthesis in a membrane bioreactor. J. Am. Oil Chem. Soc., 1990, 67, 347-352.

299. Vicente, G.; Martinez, M.; Aracil, J., Esteban A., Kinetics of Sunflower Oil Methanolysis. M Eng. Chem. Res.,44,2005, 5447-5454.

300. Wang, Y. J., Wang, Y., Chen, F., Luo, G. S., Dai, Y. Y. Mass transfer characteristics of cadmium (II) extraction in hollow fiber modules. Chemical Engineering Science, 2003, Vol. 58, 3223-3231.

301. Weaver, J. M. Plastic film heat exchangers. Chem. Eng. Prog. 1960, 56 (7), 49.

302. Werdelmann, B. W., Schmid, R. D. The biotechnology of fats-a challenge and an opportunity. Fette Seifen Anstrichm., 1982, 84(11), 436-443.

303. Whitley, D. M. Plastic heat exchangers gain in severe service. Chem. Eng. 1957, 64 (9), 308.

304. Wickramasinghe, S. R., Semmens, M. J., Cussler, E. L. Better hollow fiber contactors. J. Membrane Sei. 1991, 62,371-388.

305. Wickramasinghe S.R., Semmens M. J., Cussler E.L. Mass transfer in various hollow-fiber geometries. J. Membr. Sei., 1992, 69, 235.

306. Wickramasinghe S.R., Semmens M. J., Cussler E.L. Hollow-fiber modules made with hollow-fiber fabric. J. Membr. Sei., 1993, 84, 1.

307. Wu, J., Chen, V, Shell-side mass transfer performance of randomly packed hollow fiber modules. Journal of Membrane Science 2000, 172, 59-74.

308. Yang M.C., Cussler E.L. Designing hollow-fiber contactors. A.I.Ch.E.J., 1986, 32, 1910-1916.

309. Yang, W., Cicek, N., Ilg, J., State-of-the-art of membrane bioreactors: Worldwide research and commercial applications in North America, J. Membrane Sei., 2006, 270, 201-211.

310. Yeh H.-M, Chen Y.-K. The effect of multipass arrangement on performance in membrane extractor of fixed configuration. Chem. Eng. Sei., 2000, 55, 5873-5880.

311. Yeh H.M., Huang C.M. Solvent extraction in multipass parallel-flow mass exchangers of microporous hollow-fiber modules. J. Membrane Sci., 1995, 103, 135-150.

312. Yeh H.M., Hsu Y.S. Analysis of membrane extraction through rectangular mass exchangers. Chemical Engineering Science, 1999, 54, 897-908.

313. Yeh H.-M, Peng Y.Y., Chen Y.-K. Solvent extraction through a doublepass parallel-plate membrane cannel with recycle. Journal of Membrane Science, 1999, 163, 177.

314. Yoshikawa, S., Ogawa, J.C. Minegishi, S. Distributions of pressure and flow rate in a hollow-fiber membrane module for plasma separation. J. Chem. Engng Japan 1994,27,385-390.

315. Zabaloy M.S.; Vera J.H. Cubic Equation of State for Pure Compound Vapor Pressures from the Triple Point to the Critical Point. Ind. Eng. Chem. Res., 1996,35, 829-836.

316. Zarkadas D. M., Sirkar K. K. Polymeric Hollow Fiber Heat Exchangers: An Alternative for Lower Temperature Applications. Ind. Eng. Chem. Res. 2004, 43, 8093-8106.

317. Zhang L.-Z. Heat and mass transfer in a cross-flow membrane-based enthalpy exchanger under naturally formed boundary conditions. International Journal of Heat and Mass Transfer 2007, 50, 151 -162. : ;

318. Zhang, Q., Cussler, E. L. Microporous hollow membranes for gas absorption I: Mass transfer in the liquid. Journal of Membrane Science, 1985, 23, 321-332.

319. Zhang, Q., Cussler, E. L. Microporous hollow membranes for gas absorption II: Mass transfer across the membrane. Journal of Membrane Science, 1985, 23, 333-345.

320. Zizovic, I., Skala, D., Calculations of the solubility of vegetable oils based on cubic equations of state, 12th International Congress of Chemical and Process Engineering (CHISA^). Praha, Czech Republic, 25-30 August 1996.

321. Zosel K. Process for deodorizing fats and oils. U.S. Patent 4156688, Jul. 11, 1977.

322. Zosel K. Process for simultaneous hydrogenation and deodorisation of fats and/or oils. U.S. Patent 3969382, Aug. 29,1974.

323. Zaheed L., Jachuck R J. J., Review of polymer compact heat exchangers, with special emphasis on a polymer film unit, Applied Thermal Engineering, Volume 24, Issue 16, November 2004, Pages 2323-2358.