Разработка расчетных моделей и численный анализ неоднородных режимов течения в каналах индукционных МГД-машин тема диссертации и автореферата по ВАК РФ 01.04.13, кандидат физико-математических наук Обухов, Денис Михайлович

• Обухов, Денис Михайлович
• кандидат физико-математических наук
• 2006

Оглавление диссертации кандидат физико-математических наук Обухов, Денис Михайлович

Введение.

Глава 1. Состояние вопроса и задачи исследования

1.1 Особенности индукционных МГД-машин.

1.2 Конструкция цилиндрического линейного индукционного МГД-насоса

1.3 Уравнения магнитной гидродинамики. Критерии подобия.

1.4 МГД-неустойчивость течения в каналах ЭМН.

1.5 Экспериментальные исследования МГД-неустойчивости.

1.6 Теоретические исследования МГД-неустойчивости.

1.6.1 Модель двух насосов.

1.6.2 Одномерная струйная модель течения.

1.6.3 Двумерные модели течения.

1.7 Выводы по главе 1.

Глава 2. Описание расчетных моделей

2.1 Введение.

2.2 Одномерная струйная модель течения с двумерным распределением магнитного поля EMP2D-JET.

2.2.1 Описание модели.

2.2.2 Расчет индуцированного магнитного поля.

2.2.3 Расчет электромагнитного давления.

2.2.4 Решение гидродинамической задачи.

2.2.5 Расчет амплитуды пульсаций давления с удвоенной частотой источника питания.

2.2.6 Расчет энергетических характеристик.

2.2.7 Выбор некоторых параметров модели.

2.3 Полностью двумерная модель EMP-MHD2D.

2.3.1 Описание модели.

2.3.2 Решение электромагнитной задачи.

2.3.3 Решение гидродинамической задачи.

2.3.4 Численные методы решения электромагнитной и гидродинамической задач.

2.4 Обмотка индуктора и построение приложенного магнитного поля.

2.4.1 Параметры обмотки индуктора.

2.4.2 Расчет внешнего магнитного поля.

2.5 Выводы по главе 2.

Глава 3. Результаты расчетов по разработанным моделям

3.1 Введение.

3.2 Результаты расчета по модели EMP2D-JET.

3.2.1 Распределение магнитного поля по длине канала.

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

3.2.3 Интегральные и пульсационные характеристики насосов.

3.3 Результаты расчета по модели EMP-MHD2D.

3.3.1 Распределение магнитного поля по длине канала.

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

3.3.3 Интегральные и пульсационные характеристики насоса ЦЛИН-В

3.4 Выводы по главе 3.

Введение диссертации (часть автореферата) на тему «Разработка расчетных моделей и численный анализ неоднородных режимов течения в каналах индукционных МГД-машин»

Актульность проблемы. Магнитогидродинамические (МГД) машины различного типа с жидкометаллическим рабочим телом нашли широкое применение в таких областях техники как ядерная энергетика, металлургия и химическая промышленность [1,2,3,4]- Перспективы развития ядерной энергетики во многом связаны с реакторами на быстрых нейтронах, обладающими возможностью воспроизводить ядерное топливо и позволяющими более полно использовать природный уран [5]. Теплоноситель в реакторе на быстрых нейтронах должен быть слабым замедлителем и обеспечивать отвод тепла при высокой плотности энерговыделения. Жидкие металлы отвечают этим требованиям. Наибольшее применение в качестве жидкометаллического теплоносителя в экспериментальных и промышленных реакторах нашел натрий, имеются попытки использовать также сплав натрий-калий, ртуть, свинец.

В настоящее время индукционные МГД-насосы или электромагнитные насосы (ЭМН) используются во вспомогательных системах реакторов на быстрых нейтронах и сопутствующих им экспериментальных и о технологических установках. ЭМН с расходом натрия до 500 м /ч будут применены в системах расхолаживания основных контуров строящегося реактора БН-800 и проектируемого БН-1800. Кроме того, рассматривается возможность применения электромагнитных насосов с расходом 25000 м /ч в качестве главных циркуляционных насосов реактора БН-1800.

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

Среди ЭМН большой мощности наибольшее распространение получили трехфазные цилиндрические линейные индукционные насосы (ЦЛИН). В нашей стране в НИИЭФА им. Д.В. Ефремова были разработаны и успешно испытаны насосы ЦЛИН-5/700, ЦЛИН-5/850, ЦЛИН-8/1200, ЦЛИН-3/3500 с расходами 700, 850, 1200, 3600 м /ч и развиваемыми давлениями 0.5, 0.5, 0.8 и 0.3 МПа соответственно. Насос ЦЛИН-5/700 был в 1977 г, установлен во второй контур реактора БОР-бО, где безаварийно проработал 34000 ч. Самым мощным ЭМН, разработанным в нашей стране, является ЦЛИН-3/3500, предназначенный для одной из петель второго контура реактора БН-350 [6]. Насос развивает давление 0.3 МПа, расход натрия 3600 м /ч при его температуре 300 °С, потребляемая мощность 980 кВт, напряжение 650 В при частоте тока в сети 50 Гц, КПД 30%, масса 17000 кг. Самым мощным в мире на сегодняшний день ЭМН является индукционный насос с каналом цилиндрического типа с расходом натрия 9600 м3/ч и давлением 0.25 МПа [7], разработанный и испытанный в США.

Для мощных ЭМН характерны большие значения расхода, развиваемого давления и крупные размеры, следствием чего является увеличение характерных параметров - магнитного числа Рейнольдса Rm, числа Гартмана На, параметра МГД-взаимодействия N. Увеличение магнитного числа Рейнольдса свидетельствует об усилении индуцированного магнитного поля, что приводит к существенной зависимости характеристик ЭМН от характера распределения токов. При превышении некоторого критического значения (Rm)crit возникает так называемая МГД-неустойчивость течения рабочего тела, связанная с образованием сильно неоднородного профиля скорости. Это отражается на рабочей характеристике - развиваемое насосом давление при наличии МГД-неустойчивости имеет характерный провал. Следствием МГД-неустойчивости являются низкочастотные колебания давления, расхода и других параметров насоса и сопутствующие им вибрации насоса и контура, недопустимые на практике. Кроме того, для индукционных ЭМН при любых Rm характерны колебания параметров с удвоенной частотой питающего обмотку тока. Все перечисленные факторы отрицательно влияют на эффективность ЭМН большой мощности.

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

Существует ряд работ по экспериментальному исследованию ЭМН, работающих при больших параметрах электромагнитного взаимодействия. Так в НИИЭФА им. Д.В. Ефремова были проведены экспериментальные исследования нескольких мощных цилиндрических линейных индукционных насосов, работающих при расходах металла до 3600 м /ч, в том числе в условиях существования МГД-неустойчивости течения [6,8,9,10,11]. Эти эксперименты включали в себя исследование структуры магнитного поля в зоне индуктора, интегральных характеристик насоса, влияния конструктивных особенностей ЦЛИН на проявление МГД-неустойчивости, взаимосвязь между явлением МГД-неустойчивости и низкочастотными колебаниями параметров насоса.

Специфический характер напор-расходных характеристик ЭМН с плоским каналом при конечных магнитных числах Рейнольдса также обнаружен в ряде работ, например в [12,13].

Однако, теоретические аспекты явления МГД-неустойчивости были исследованы недостаточно. Существует ряд теоретических работ на основе одномерной струйной модели течения [14,15,16,17,18]. Одномерные модели позволяют получать интегральные и часть локальных характеристик ЭМН, но ввиду заложенных в них ограничений не позволяют в полном объеме описать явление МГД-неустойчивости. Для получения более полной картины неоднородного течения металла в канале ЭМН необходимо совместное решение двумерной электромагнитной и гидродинамической задач. Из-за наличия конвективных слагаемых задача является нелинейной и ее аналитическое решение в общем случае невозможно. Таким образом, построение более совершенных расчетных моделей, основанных на совместном численном решении двумерных электромагнитной и гидродинамической задач, и анализ с их помощью режимов течения в каналах мощных ЭМН является актуальной задачей. Целью работы являлось:

1. Построение расчетных моделей, способных описывать неоднородные течения жидкого металла в каналах индукционных ЭМН с учетом индуцированных в металле токов;

2. Анализ возникновения и развития неоднородного режима течения металла по длине и ширине канала ЭМН;

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

4. Оценка критических значений параметров МГД-взаимодействия, определяющих границу области устойчивой работы насоса;

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

Научная новизна и результаты, выносимые на защиту

1) Разработаны и реализованы в виде вычислительных программ (EMP2D-JET и EMP-MHD2D) квази-двумерная и полностью двумерная модели течения металла в канале электромагнитного насоса.

2) С использованием разработанных моделей проанализирован процесс потери устойчивости однородного течения в канале ЭМН; получены критические значения характерных параметров процесса, позволяющие определять границы перехода к неоднородному режиму течения.

3) Показано, что при работе насоса в области МГД-неустойчивости течение в канале имеет вихревой характер.

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

Практическая ценность Разработанные модели используются для расчетов характеристик мощных индукционных ЭМН, в частности для реакторов БН-800 и БН-1800. Они позволяют:

- рассчитывать основные локальные и интегральные характеристики индукционных ЭМН с учетом развития МГД-неустойчивости течения;

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

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

Результаты диссертации опубликованы в 4 работах [19,20,21,22], докладывались на международной конференции по «Фундаментальной и прикладной магнитной гидродинамике» (Рига, 2005 г.)

Объем и структура диссертации. Диссертация состоит из введения, трех глав и заключения. Диссертация изложена на 100 машинописных листах, содержит 34 рисунка, 6 таблиц и список литературы из 40 наименований.

Заключение диссертации по теме «Электрофизика, электрофизические установки», Обухов, Денис Михайлович

Основные результаты работы сформулировать следующим образом.

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

2. Проведенный анализ показал, что превышение критических значений характерных параметров исследуемых процессов - магнитного числа Рейнольдса Rms и параметра МГД-взаимодействия N, приводит к формированию вихревого течения в канале ЭМН.

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

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

5. Разработанные модели позволяют качественно описать пульсации давления с удвоенной частотой источника питания и низкой частоты, что также является существенным достижением, так как пульсации давления низкой частоты впервые охарактеризованы количественно при помощи расчетных моделей, а исследуемое явление МГД-неустойчивости достаточно сложно для анализа.

Заключение

Список литературы диссертационного исследования кандидат физико-математических наук Обухов, Денис Михайлович, 2006 год

1. Глухих В.А., Таиаиаев А.В., Кириллов И.Р. Магнитная гидродинамика в ядерной энергетике. М.: Энергоатомиздат. 1987.

2. Верте JI.A. Магнитная гидродинамика в металлургии. М.Металлургия. 1975.

3. Баринберг А.Д., Капуста А.Б., Чекин Б.В. Обзор работ донецких исследователей по прикладной магнитной гидродинамике. Магнитная гидродинамика. 1975. №1. с. 111-117.

4. Полищук В.П. Промышленное использование магнитодинамических устройств. Магнитная гидродинамика. 1975. №1. с. 118-128.

5. Александров А.П. Технические аспекты ядерной энергетики на грани веков. Атомная энергия. 1984. Т.56 Вып. 6 С.339-342.

6. Андреев A.M., Безгачев Е.А., Карасев Б.Г., Кириллов И.Р., Огородников А.П., Преслицкий Г.В., Чвартацкий Р.В. Электромагнитный насос ЦЛИН-3/3500. Магнитная гидродинамика. 1988. №1. с.61-67.

7. Ota Н., Katsuki К., Funato М., Taguchi J., Fanning A., Doi Y., Nibe N., Uetaо

8. M., Inagaki T. Development of 160 m /min large capacity sodium-immersed self-cooled electromagnetic pump. Journal of nuclear science and technology. 2004. Vol.41, No.4,p.511-523

9. Андреев A.M., Карасев Б.Г., Кириллов И.Р., Огородников А.П., Остапенко В.П., Семиков Г.Т. Результаты экспериментального исследования электромагнитных насосов для установки БОР-бО. Магнитная гидродинамика, 1978, №4, с. 93-100.

10. Кириллов И.Р., Огородников А.П., Остапенко В.П. Экспериментальное исследование неоднородности течения в цилиндрическом линейном индукционном насосе. Магнитная гидродинамика. 1980. №2. с. 107-113.

11. Araseki Н., Kirillov I.R, Preslitsky G.V., Ogorodnikov А.Р. Double-supply-frequency pressure pulsation in annular linear induction pump, part I:measurement and numerical analysis. Nuclear Engineering and Design. 2000. 195, 85-100.

12. Araseki H., Kirillov I.R, Preslitsky G.V., Ogorodnikov A.P. Magnetohydrodynamic instability in annular linear induction pump, part I. Experiment and numerical analysis. Nuclear Engineering and Design. 2004. 227, 29-50.

13. Валдмане P.A., Кришберг P.P., Лиелпетер Я.Я. и др. Локальные характеристики течения в канале индукционной МГД-машины при больших параметрах электромагнитного взаимодействия. // Магнитная гидродинамика, 1977, № 3, с. 99-104.

14. Валдмане Р.А., Кришберг P.P., Лиелпетер Я.Я. и др. Интегральные характеристики индукционной МГД-машины при больших параметрах электромагнитного взаимодействия.// Магнитная гидродинамика, 1977, №4, с. 107-109.

15. Гайлитис А., Лиелаусис О. О внутренней гидравлике МГД-машин при неоднородном распределении сил.// Магнитная гидродинамика, 1971, № 2 с. 123-127.

16. Валдмане Р.А., Кириллов И. Р., Огородников А.П., Остапенко В.П., Улманис Л.Я. Расчет характеристик индукционного МГД-насоса при Rms>l с учетом неоднородности распределения внешнего магнитного поля. Магнитная гидродинамика. 1982. №3. с.98-104

17. Вадмане Р., Лиелаусис О., Улманис Л. Модель неоднородного течения в канале индукционного насоса. // Магнитная гидродинамика, 1983,№ 2,с. 98-102.

18. Половко Ю.А., Романова Е.П., Тропп Э.А. Исследование устойчивости индукционных течений в МГД насосах и генераторах с широким каналом. Журнал технической физики. 1996. Том 66. №4. с. 36-44

19. Половко Ю.А., Романова Е.П., Тропп Э.А. Индукционная цилиндрическая МГД машина в режиме идеального источники давления. Журнал технической физики. 1997. Том 67. №6. с.5-8

20. Kirillov I.R., Obukhov D.M. Two dimensional model for analysis of cylindrical linear induction pump characteristics: model description and numerical analysis. Energy conversion and management. 2003. 44. p.2687-2697.

21. Кириллов И.Р., Обухов Д.М. Полностью двумерная модель для анализа характеристик линейного цилиндрического индукционного насоса. Журнал технической физики. 2005. Том 75. Вып. 8.

22. Igor R. Kirillov, Denis М. Obukhov, Anatoly P. Ogorodnikov, Hideo Araseki. Comparison of computer codes for evaluation of double-supply-frequency pulsations in linear induction pumps. Nuclear Engineering and Design. 2004.231,177-185.

23. Лиелпетер Я.Я. Жидкометаллические индукционные МГД-машины. Рига: Зинатне, 1969, 246 с.

24. Вольдек А.И. Индукционные магнитогидродинамические машины с жидкометаллическим рабочим телом. Л.: Энергия, 1970. 272 с.

25. Баранов Г.А., Глухих В.А. Кириллов И.Р. Расчёт и проектирование индукционных МГД-машин с жидкометаллическим рабочим телом. М.: Атомиздат, 1978,247 с.

26. Кириллов И.Р., Огородников А.П., Остапенко В.П. Экспериментальное исследование неоднородности течения в цилиндрическом линейном индукционном насосе. Магнитная гидродинамика. 1980. №2, с. 107-113.

27. Шерклиф Дж. Курс магнитной гидродинамики. М: Мир. 1967

28. Ватажин А.Б., Любимов Г.А., Регирер С.А. Магнитогидродинамические течения в каналах. М.: Наука. 1970

29. Кириллов И.Р., Огородников А.П., Остапенко В.П., Преслицкий Г.В. Исследования влияния конструктивных особенностей на течение в канале цилиндрического индукционного насоса. Магнитная гидродинамика. 1981. №2, с.79-84.

30. Кириллов И.Р., Остапенко В.П. Локальные характеристики цилиндрического индукционного насоса при Rm s>l. Магнитная гидродинамика. 1987. №2. с. 95-102.

31. Гайлитис А., Лиелаусис О. Неустойчивость однородного распределения скоростей в индукционной машине.// Магнитная гидродинамика, 1975, № 1, с. 87-101.

32. Кириллов И.Р. Учет проводимости стенок в задачах о продольном концевом эффекте в индукционных каналах. Магнитная гидродинамика. 1973. №4. 117-122.

33. Флетчер К. Вычислительные методы в динамике жидкостей. T.l М.: Мир. 1991. 504 с.

34. Андерсон Д., Таннехил Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. М.: Мир, 1990. Т.2. 728 с.

35. Флетчер К. Вычислительные методы в динамике жидкостей. Т.2 М.: Мир. 1991. 552 с.

36. Марчук Г.И. Методы расщепления. М.: Наука. 1988. 264 с.

37. Марчук Г.И. Методы вычислительной математики. М.: Наука. 1989. 608 с.

38. Идельчик И.Е. Справочник по гидравлическим сопротивлениям. M.-JI. Госэнергоиздат.1960.464 с.