Концентрационная конвекция в процессе обратного осмоса при ламинарном течении в плоских каналах тема диссертации и автореферата по ВАК РФ 05.17.08, кандидат технических наук Кириченко, Николай Борисович

• Кириченко, Николай Борисович
• кандидат технических наук
• 2006

Оглавление диссертации кандидат технических наук Кириченко, Николай Борисович

Введение.стр.

Глава 1. Анализ решений задачи устойчивости в тепло-и массообмене.стр.

1.1 Описание основного стационарного течения в задачах конвективного тепло-и массообмена.стр.

1.2 Теоретическое описание массообмена в баромембранных процессах.стр.

1.3 Решение задачи устойчивости в процессах тепло- и массообмена.стр.

1.4 Массообмен в смешанноконвективном режиме.стр.

1.5 Выводы. Цель работы.стр.

Глава 2. Теоретическое исследование стационарного массообмена при наличии вынужденноконвективного продольного течения.стр.

2.1 Упрощение уравнений Навье - Стокса для рассматриваемой задачи.стр.

2.2 Решение уравнения стационарной диффузии операционным методом.стр.

2.3 Расчет профиля концентраций методом численного обращения преобразования Лапласа.стр.

2.4 Обобщение полученных соотношений.стр.

2.5 Выводы.стр.

Глава 3. Теоретическое исследование концентрационной устойчивости при одномерном нестационарном массообмене в окрестности селективно-проницаемой стенки.стр.

3.1 Решение уравнения нестационарной одномерной диффузии в примембранном пространстве операционным методом.стр.

-83.1.1 Формулировка задачи. Описание стационарного состояния.стр.

3.1.2 Решение уравнения одномерного нестационарного массообмена.стр.

3.1.3 Анализ решения уравнения нестационарного массообмена.стр.

3.1.4 Описание задачи устойчивости в случае квазиэкспоненциального профиля концентраций.стр.

3.2 Математическое исследование задачи концентрационной устойчивости в плоском мембранном канале.стр.

3.2.1 Формулировка задачи.стр.

3.2.2 Разработка уравнений для анализа устойчивости.стр.

3.2.3 Алгоритм численного решения задачи устойчивости.стр.

3.2.4 Анализ случая Pev—>0.стр.

3.2.5 Численные результаты. Примеры расчета.стр.

3.3 Выводы.стр.

Глава 4. Экспериментальное исследование концентрационной устойчивости процесса обратного осмоса.стр.

4.1 Выбор объектов и методов исследования.стр.

4.2 Лабораторная установка для исследования одномерного нестационарного массообмена.стр.

4.2.1 Установка для поточного измерения концентрации проникающего раствора (ионометр с самописцем).стр.

4.2.2 Установка для визуализации распределения концентрации раствора над мембраной лазерный интерферометр) и узел фиксации интерферограмм.стр.

4.2.3 Установка для исследования процесса обратного осмоса в статических условиях (технологическая часть).стр.

4.3 Подготовка и проведение эксперимента.стр.

4.4 Методика обработки опытных данных.стр.

-94.5 Оценка погрешности расчетов и измерений.стр.

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

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

5.2 Основные экспериментальные результаты и их сравнение с теоретическим предсказанием.стр.

5.3 Выводы.стр.

Глава 6. Прогноз эффективности использования смешанноконвективного режима массообмена в проточных плоскокамерных баромембранных аппаратах.стр.

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

6.2 Предсказание массообменных характеристик проточных установок обратного осмоса в смешанноконвективном режиме.стр.

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

6.4 Выводы.стр.

Основные результаты работы.стр.

Введение диссертации (часть автореферата) на тему «Концентрационная конвекция в процессе обратного осмоса при ламинарном течении в плоских каналах»

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

В настоящее время развитие мембранной технологии приняло особенно динамичный характер, благодаря разработке новых высокопроизводительных композитных мембран, например мембран типа ОПМН (НПО "Полимерсинтез", г. Владимир), которые иногда называют "нанофильтрационными". На самом деле мембраны типа ОПМН легко модифицируются в заводских условиях как для процессов ультра- и нанофильтрации так и для процесса обратного осмоса. Естественно, что с нанесением каждого последующего слоя полимера на композитную мембрану типа ОПМН вместе с появлением селективных свойств по отношению к частицам всё меньших размеров, наблюдается и снижение производительности мембраны. Однако указанное снижение происходит всего лишь в 2 - 4 раза при переходе в ряду: ультрафильтрация - нанофильтрация - обратный осмос, что нельзя не признать выдающимся достижением.

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

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

Однако именно увеличение производительности мембран при заданном уровне селективности ведёт к резкому (практически экспоненциальному) росту внешнедиффузионного сопротивления массопереносу, которое в мембранной литературе принято называть "концентрационной поляризацией" [1]. Явление концентрационной поляризации всегда приводит к снижению эффективности мембранного аппарата по крайней мере в двух отношениях: а) снижается эффективность разделения, т. е. увеличивается средняя концентрация пермеата б) снижается производительность мембранного аппарата - т. е. количество пермеата с единицы объёма аппарата. Кроме того, увеличение концентрации растворённых веществ может привести к осадкообразованию на мембране и снижению её характеристик и срока службы.

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

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

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

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

Повышение скорости течения, согласно простейшим расчётам, которые можно провести на основании плёночной модели (1.7) - (1.8), довольно слабо влияет на величину концентрационной поляризации; вместе с тем увеличение скорости сопряжено с повышенными затратами энергии; также при высокой скорости в напорном канале возможна деформация мембраны взвешенными частицами и её порча.

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

В настоящей работе будут разработаны теоретические основы метода снижения концентрационной поляризации в баромембранных процессах, основанного на явлении естественной конвекции. Данный метод не требует каких-либо изменений в существующих типах мембранных модулей, не предполагает увеличения эксплуатационных затрат и не требует установки дополнительных механизмов. На основании теоретических и экспериментальных исследований процесса возникновения естественной конвекции в плоских мембранных каналах будет показано, что режим естественной конвекции должен возникать в большинстве имеющихся типов мембранных модулей. Для возникновения данного режима необходимо сочетание трёх условий: высокой (>5 г/л) концентрации исходного раствора, определённой ориентации модуля относительно поверхности Земли и низкой скорости течения раствора в напорном канале.

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

Заключение диссертации по теме «Процессы и аппараты химической технологии», Кириченко, Николай Борисович

Основные результаты работы

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

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

2. Подтверждено экспериментальным путем теоретическое число (модифицированное число Рэлея) потери устойчивости в каналах плоских обратноосмотических модулей Rac4 ~ 100.

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

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

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

2. С использованием аппарата преобразований Лапласа получено аналитическое решение задачи о нестационарном массообмене вблизи мембраны.

3. Методом Галёркина решена задача о концентрационно-гравитационной устойчивости горизонтального диффузионного пограничного слоя с экспоненциальным распределением концентрации. При решении данной задачи разработан математический аппарат численного решения задачи устойчивости концентрационных (температурных) пограничных слоев с произвольным распределением концентрации (единственное ограничение - данное распределение должно в пределе Z—>оо переходить в экспоненциальное). В процессе решения получены уравнения возмущений и амплитудные уравнения для анализа устойчивости концентрационных пограничных слоев (при наличии постоянной в пределах слоя скорости отсоса).

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

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

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

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

Список литературы диссертационного исследования кандидат технических наук Кириченко, Николай Борисович, 2006 год

1. Дытнерский Ю. И. Обратный осмос и ультрафильтрация, М., Химия, 1978. -352 с.

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

3. Шлихтинг Г. Теория пограничного слоя, М., Наука, 1974. 712 с.

4. Эккерт Э. Р., Дрейк Р. М. Теория тепло- и массообмена, М., Госэнергоиздат, 1961.-679 с.

5. Кем В. А. Массопередача в плоских каналах с селективно-проницаемыми стенками применительно к процессам газоразделения. Дисс. канд. тех. наук, М., МХТИ им. Д. И. Менделеева, 1989. 222 с.

6. Кастро А. X. X. Массообмен при обтекании проницаемого горизонтального цилиндра. Дисс канд. тех. наук, М., МХТИ им. Д. И. Менделеева, 1992. 189 с.

7. Калика Г. 3. Массообмен в каналах с селективнопроницаемыми стенками в условиях концентрационной неустойчивости (применительно к процессам мембранного газоразделения). Дисс. канд. тех. наук, М., МХТИ им. Д. И. Менделеева, 1993. 210 с.

8. Sourirajan S. Reverse Osmosis, London, Logos, 1970. 578 p.

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

10. Тарарышкин М. В. Внешний массоперенос в процессе обратного осмоса при ламинарном течении в плоских каналах. Дисс. канд. тех. наук, М., 1995. 151 с.

11. Berman A. S. Laminar Flow in Channels with Porous Walls// J. Appl. Phys. 1953, v. 24, №9, p. 1232-1235

12. Лыков А. В. Тепломассообмен. Справочник, M., Энергия, 1971. С. 210

13. Bhattacharyga D. and oth., Prediction of Concentration Polarization.// J. of Membrane Science, 1990, v. 48, № 2-3, p. 231-262

14. Sherwood Т.К., Brian P. L. Т., Fisher R. E., Dresner L. Salt Concentration at the Phase Boundaries in Desalination by Reverse Osmosis// Ind. Eng. Chem. Fund., 1965, v. 4, №2, p. 113-118

15. Principles of Desalination. Ed. By K.S. Spegler, New York, Academic Press,1966.-566 p.

16. Brian P. L. T. Concentration Polarization in Reverse Osmosis.// Ind. Eng. Chem. Fund., 1965, v. 4, №4, p. 439-445

17. Shah Y. Mass Transport in Reverse Osmosis in Base of Variable Diffusivity// Int. J. Heat and Mass Transfer, 1971, v. 14, № 7, p. 921-930

18. Gill W. N., Tien C., Zeh D. Concentration Polarization Effects in a Reverse Osmosis System// Ind. Eng. Chem. Fund., 1965, v. 4, № 4, p. 433-439

19. Srinivasan S., Tien C. Simultaneous Development of Velocity and Concentration Profiles in Reverse Osmosis Systems// Chem. Eng. Science, 1967, v. 22, № 3, p. 417 -433

20. Kimura S., Sourirajan S.// Amer. Interface of Chem. Eng. J., 1967, v. 19, № 13, p. 497 502

21. Kimura S., Sourirajan S.// Ind. Eng. Chem. Pr. Des. and Develop., 1968, v. 7, № 1, p. 41 -48

22. Kimura S., Sourirajan S.// Ind. Eng. Chem. Pr. Des. and Develop., 1968, v. 7, № 4, p. 539

23. Matsura Т., Baxler A. G., Sourirajan S.// Ind. Eng. Chem. Pr. Des. and Develop., 1977, v. 16, № l,p. 82-89

24. Поляков С. В., Волгин В. Д., Максимов Е. Д., Синяк Ю. Е. Расчет концентрационной поляризации в аппаратах обратного осмоса с плоскокамерными фильтрующими элементами // Химия и технология воды, 1982, т. 4, № 3, С. 299 303

25. Рид Т., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей, JL, Химия, 1982. 592 с.

26. Strathmann Н., Keilin В. Control of Concentation Polarization in Reverse Osmosis Desalination of Water// Desalination, 1969, v. 6, № 2, p. 179-201

27. Дмитриев E. А. Исследование явления концентрационной поляризации и его учет в процессах разделения растворов обратным осмосом. Дисс. канд. тех. наук, М., МХТИ им. Д. И. Менделеева, 1980. 179 с.

28. Gilron J., Hasson D. Analysis of Laminar Flow Precipitation Fouling on Reverse

29. Osmosis Membranes// Desalination, 1973, v. 12, № 2, p. 127 129

30. Ramanandha K., Gill W. NM Amer. Interface of Chem. Eng. J., 1969, v. 19, №15, p. 872-874

31. Sloan E. D., Harshman R. CM Amer. Interface of Chem. Eng. J., 1973, v. 19, №3, p. 618-623

32. Jonsson A. R., Acrivos A. Concentration Polarization in Reverse Osmosis under Natural Convection // Ind. Eng. Chem. Fund., 1969, v. 8, № 2, p. 359

33. Jonsson A. R., Acrivos A. Experimental Investigation of Polarization Effect in Reverse Osmosis // Amer. Interface of Chem. Eng. J., 1974, v. 20, № 5, p. 966-974

34. Ляпунов A. M. Общая задача об устойчивости движения. Дисс. докт. мат. наук, Харьков, 1892

35. Гершуни Г. 3., Жуховицкий Е. М. Конвективная устойчивость несжимаемой жидкости, М., Наука, 1972. 392 с.

36. Platten J.-K. Convection in Liquids. Bruxelles, 1977. 697 p.

37. Le Gal P., Croquette V. Appearance of a square pattern in a Rayleigh-Benard experiment // Phys. Fluids, 1988, № 31 (11), p. 3440 3442

38. Норден П. А. Экспериментальное исследование условий возникновения конвекции в горизонтальных жидкостных слоях. Автореф. дисс. канд. тех. наук, Казань, КХТИ им. С. М. Кирова, 1969. 18 с.

39. Гленсдорф П., Пригожин И. Термодинамическая теория структуры, устойчивости и флуктуаций (2-е издание), М., УРСС, 2003. 280 с.

40. Busse F. Н., Non-linear properties of thermal convection // Rep. Prog. Phys., 1978, № 41 (12), p. 1929 1967

41. Демидович Б. П., Марон И. А., Шувалова Э. 3. Численные методы анализа, М., Наука, 1967-489 с.

42. Platten J.-K. These de doctrat. Chimie Physique, Bruxelles, Universite de Bruxelles, 1970.-493 p.

43. Петухов Б. С., Поляков А. Ф. Теплообмен при смешанной турбулентной конвекции, М., Наука, 1986. 192 с.

44. Шварцблат Д. Л. Стационарные конвективные движения в плоскомгоризонтальном слое жидкости с проницаемыми границами. Изв. АН СССР, МЖГ, 1969, № 5

45. Шварцблат Д.Л., Численное исследование стационарного конвективного движения в плоском горизонтальном слое жидкости. В Сб.: Гидродинамика (вып. 3). Издательство Пермского Государственного Университета, 2001

46. Костырев К. Г. Концентрационная конвекция в растворах сильных электролитов. Дисс. канд. физ.-мат. наук, Пермь, Издательство Пермского Государственного Университета, 1992. 158 с.

47. Ерошенко В. М., Зайчик Л. И., Рабовский В. Б. Об устойчивости течения жидкости в плоском канале с равномерным вдувом или отсосом через проницаемые стенки // Инж.-физ. журнал, 1981, т. 41, № 3. С. 436 - 440

48. Непомнящий А. А. Расчет коэффициентов в амплитудном уравнении для вторичных течений, С. 3 8 // Исследование тепловой конвекции и теплопередачи, Свердловск, изд. УНЦ АН СССР, 1981. - 96 с.

49. Ландау Л. Д., Лившиц Е. М. Теоретическая физика. Т. 6, Гидродинамика, М., Наука, 1986 (3-е изд.). 736 с.

50. Eckhaus W. Studies in Non-Linear Stability Theory. In: Springer Tracts in Natural Philosofy, v. 6, Berlin, Springer, 1965

51. Palm E. On the tendency towards hexagonal cells in steady convection // J. Fluid Mech., 1960, v. 8, № 2, p. 183 192

52. Oliver D. S., Booker J. R., Planform of convection with strongly temperature-dependent viscosity, Geophys. Astrophys. Fluid Dyn., 1983, v. 27, № 1-2, p. 73 85

53. Scanlon J. W., Segel L. A. Finite amplitude cellular convection induced by surface tension // J. Fluid Mech., 1967, v. 30, № 1, p. 149 162

54. Krishnamurti R. Finite amplitude convection with changing mean temperature (Part 1. Theory; Part 2. Experimental test of the theory) // J. Fluid Mech., v. 33, № 3, pp. 445-455, 457 463

55. Гершуни Г. 3., Жуховицкий E. M. Устойчивость конвективных течений, М., Наука, 1989.-318 с.

56. Newell А. С., Whitehead J. A. Finite bandwidth, finite amplitude convection, J.-206

57. Fluid Mech., 1969, v. 38, № 2, p. 279 303

58. Swift J., Hohenberg P. C. Hydrodynamic fluctuations at the convective instability // Phys. Rev., 1977, v. A15, № 1, p. 319 328

59. Jones W. P., Launder В. E. // Int. J. Heat Mass Transfer, 1972, v. 15, № 193

60. Whitehead J. A., Chen M. M. Thermal instability and convection of a thin fluid layer bounded by a stably stratified region // J. Fluid Mech., v. 40, № 3, 1970, p. 549 -576

61. Пасконов В. M., Полежаев В. И., Чудов JI. А. Численное моделирование процессов тепло- и массообмена, М., Наука, 1984. 285 с.

62. Clever R. М., Busse F. Н. Nonlinear oscillatory convection // j. Fluid Mech., v. 176, 1987.-p. 403-417

63. Келлер Д., Фридман А. Дифференциальные уравнения турбулентного движения сжимаемой жидкости // Методы расчета турбулентных течений. Пер. с англ. под ред. А.Д. Хонькина, М., Мир, 1984. 366 с.

64. Непомнящий А. А. Расчет коэффициентов в амплитудном уравнении для вторичных течений, С. 3 9 // Исследование тепловой конвекции и теплопередачи (Сборник статей), Свердловск, УНЦ АН СССР, 1981. - 96 с.

65. Getling А. V. Evolution of two-dimensional disturbances in the Rayleigh -Benard problem and their preffered wavenumbers // J. Fluid Mech., 1983, v. 130, p. 165-186

66. Дзюбенко В. Г., Бон А. И., Кадыкова Н. Е. Обратноосмотические мембраны среднего давления. Тез. докл. Всес. конф. Мембранные методы разделения смесей. Черкассы, 23 27 декабря 1991, С. 62 - 63

67. Зайцев В. Ф., Полянин А. Д. Справочник по линейным обыкновенным дифференциальным уравнениям, М., Факториал, 1997. 304 с.

68. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции (т. 1 Гипергеометрические функции, Функции Лежандра), М., Наука, 1973. 296 с.

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

70. Свешников А. Г., Тихонов А. Н. Теория функции комплексного переменного, М., Наука, 1999.-417 с.-20769. Диткин В. А., Прудников А. П. Интегральные преобразования и операционное исчисление, М., Физматгиз, 1961. 341 с.

71. Диткин В. А., Кузнецов П. И. Справочник по операционному исчислению, М., Гостехиздат, 1951. 255 с.

72. Крылов В. И., Шульгина JI. Т. Справочная книга по численному интегрированию, М., Наука, 1966. 514 с.

73. Jonsson A. R. Experimental Investigation of Polarization Effect in Reverse Osmosis. // AIChE J., 1974, v. 20 №5, p. 966 974

74. Себиси Т., Бредшоу П. Конвективный теплообмен, М., Мир, 1987. 590 с.

75. Федоренко Р. П. Вычислительная физика, М., Изд. МФТИ, 1999. 512 с.

76. Рыжик И. М., Градштейн И. С. Таблицы интегралов, рядов, сумм и произведений (3-е издание), М., Гостехиздат, 1951. 805 с.

77. Макаров Е.Д. Инженерне расчеты в MathCAD. Учебный курс, СПб., Издательский дом «Питер», 2005. 448 с.

78. Хауф В., Григуль У. Оптические методы в теплопередаче, пер. с англ. под ред. проф. Лихушина В.Я., М., Мир, 1973. 240 с.

79. Петров Ф.А. Голографическая интерференция, М., Химия, 1990. 442 с.

80. Сенявин М. М. Ионный обмен в технологии и анализе неорганических веществ, М., Химия, 1980. 358 с.

81. Иоффе Б. В. Рефрактометрические методы химии (2-е издание), Л., Химия, 1974.-400 с.

82. Дытнерский Ю.И., Брыков В. П., Каграманов Г. Г. Мембранное разделение газов, М., Химия, 1991. 344 с.

83. Aluminium Dusseldorf J., 2005, № 6, p. 22