авторефераты диссертаций БЕСПЛАТНАЯ РОССИЙСКАЯ БИБЛИОТЕКА - WWW.DISLIB.RU

АВТОРЕФЕРАТЫ, ДИССЕРТАЦИИ, МОНОГРАФИИ, НАУЧНЫЕ СТАТЬИ, КНИГИ

 
<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

Pages:   || 2 | 3 |

Исследование гидродинамической неустойчивости и турбулентного перемешивания в задачах лазерного термоядерного синтеза

-- [ Страница 1 ] --

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

ЯХИН РАФАЭЛЬ АСХАТОВИЧ

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

01.04.21 – Лазерная физика

Автореферат

диссертации на соискание ученой степени

кандидата физико-математических наук

Автор:

Москва 2009

Работа выполнена в Национальном исследовательском ядерном университете "МИФИ" и Учреждении Российской академии наук Физическом институте им. П.Н. Лебедева РАН

Научный руководитель: доктор физико-математических наук, профессор Розанов Владислав Борисович

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

ведущий научный сотрудник Бельков Сергей Аркадьевич

доктор физико-математических наук,

профессор Лебо Иван Германович

Ведущая организация: Институт Математического Моделирования РАН, Москва

Защита состоится 30 сентября 2009 года в 15 час 00 мин. на заседании диссертационного совета Д 212.130.05 в аудитории K-608 Национального исследовательского ядерного университета "МИФИ" по адресу: 115409, г. Москва, Каширское шоссе, дом 31, тел. (495)324-87-66.

С диссертацией можно ознакомиться в библиотеке НИЯУ МИФИ.

Автореферат разослан «30» августа 2009 года.

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

Ученый секретарь диссертационного совета Евсеев И.В.

Общая характеристика работы

Актуальность работы.

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

В основе подхода лежит концепция сверхвысокого сжатия и нагрева термоядерных мишеней мощными лазерными импульсами. Независимо от выбранного типа мишени для достижения условий термоядерного зажигания необходимо сжимать горючее до огромных плотностей масштаба 100-2000 г/см3. Однако в процессе облучения и сжатия развиваются гидродинамические неустойчивости, препятствующие достижению оптимальных параметров термоядерного горючего. Например, в современных экспериментах, различие в наблюдаемом выходе нейтронов от предсказаний одномерных расчетов, не учитывающих влияние несимметрии, превышает два порядка при объемных сжатиях горючего более 104 раз [4,5]. Поэтому, на сегодняшний день, проблема устойчивого сжатия мишеней является одной из ключевых в лазерном термоядерном синтезе (ЛТС), а решением вопросов, связанных с гидродинамическими неустойчивостями, из которых основную роль играют неустойчивости Рэлея-Тейлора и Рихтмайера-Мешкова занимаются ведущие лазерные лаборатории мира.



На сегодняшний день не существует сравнительно простой теоретической модели для описания процессов турбулентного перемешивания, включающей зависимость от начальных условий. Следует отметить, что при общем вполне оправданном, стремлении к осуществлению высокосимметричного сжатия капсул ЛТС, представляется интересным исследование сжатия при сильных отклонениях от симметрии, т.к. стадия коллапса практически всегда оказывается несимметричной. Серьезная работа проводится по изучению развития крупномасштабного гидродинамического перемешивания и его влияния на нейтронный выход реакции [6-10].

Большое внимание, в том числе в связи со сложностью технического эксперимента и его дороговизной, уделяется численному моделированию процессов, происходящих в мишенях в процессе сжатия и горения, его сопоставлению с экспериментом и предсказаниями теоретических работ [6, 11 - 17].

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

Цели и задачи диссертационной работы.

Целями диссертационной работы являются:

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

2) исследование процессов сжатия и горения мишеней для лазерных установок "Искра V" [18] и HiPER [19-21].

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

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

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

3) проанализировать результаты эксперимента по прямому сжатию лазерных термоядерных мишеней на установке "Искра V";

4) для лазерной установки "Искра V" представить более эффективные мишени для дальнейших экспериментов, обеспечивающие более симметричное сжатие и высокий нейтронный выход;

5) провести расчеты сжатия и горения мишеней для будущей лазерной установки HiPER.

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

  1. предложена теоретическая модель, включающая информацию о начальных условиях и позволяющая описать зону перемешивания в широком диапазоне начальных условий;
  2. определена слабая зависимость развития зоны перемешивания от вклада высокочастотных гармоник;
  3. проведен анализ влияния начальных возмущений формы поверхности капсулы термоядерного синтеза на результаты сжатия, такие как нейтронный выход, средняя и максимальная плотность и температура горючего, время сжатия, минимальный радиус оболочки;
  4. для систем, обладающих осевой симметрией, предложена методика проведения численных расчетов несимметричного по полярному углу сжатия и горения мишеней лазерного термоядерного синтеза; впервые дано объяснение наблюдаемого в эксперименте низкого (в 100-1000 раз) по сравнению с данными одномерных расчетов нейтронного выхода.

Практическая ценность результатов диссертационной работы:

  1. представленная модель позволяет быстро и с хорошей точностью проводить исследования процессов в лазерных термоядерных мишенях;
  2. проведен анализ экспериментов на установке "Искра V";
  3. для планируемых экспериментов на установке HiPER проведены расчеты сжатия и горения, основанные на картах и гистограммах освещенности мишени в зависимости от числа лазерных пучков;
  4. предложены перспективные мишени для возможных экспериментов на установке "Искра V", обладающие большей степенью симметрии при сжатии и позволяющие достигать больший нейтронный выход.

На защиту выносятся следующие положения:

  1. предложенная теоретическая модель, которая описывает ширину зоны перемешивания в широком диапазоне начальных условий;
  2. методика описания процессов несимметричного по полярному углу сжатия и горения мишеней лазерного термоядерного синтеза, позволяющая быстро и с хорошей точностью описывать процессы в лазерных термоядерных мишенях;
  3. слабое влияние высокочастотных гармоник на развитие области перемешивания и снижение нейтронного выхода при сжатии мишеней лазерного термоядерного синтеза;
  4. для условий лазерной установки "Искра V" предложенные мишени, имеющие более высокую степень симметрии сжатия и позволяющие достигать больший нейтронный выход по сравнению с экспериментальными.

Апробация результатов исследования.

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

    1. Змитренко Н.В., Розанов В.Б., Степанов Р.В., Яхин Р.А. “Моделирование перемешивания и роль несимметрии облучения при сжатии лазерных термоядерных мишеней”. XXXVI Международная (Звенигородская) конференция по физике плазмы и управляемому термоядерному синтезу, (г. Звенигород, 9-13 февраля 2009 г.).
    2. The fifth International Conference on Inertial Fusion Sciences and Applications (IFSA2007). V. Rozanov, R. Yakhin, N. Zmitrenko et al. “Analysis of the Conditions of Experiments on the Compression of Capsules With a Foam Absorber at Iskra V Facility”. Journal of Physics: Conference Series 112 (2008) 022006.
    3. The fifth International Conference on Inertial Fusion Sciences and Applications (IFSA2007). V. Rozanov, R. Yakhin, N. Zmitrenko et al. “The model of turbulent mixing zone evolution, which accounts for the deviation from spherical symmetry of laser thermonuclear target compression”. Journal of Physics: Conference Series 112 (2008) 022011.
    4. Яхин Р.А., Розанов В.Б., Степанов Р.В., Змитренко Н.В. “Анализ влияния несимметрии облучения в экспериментах по сжатию капсул на лазерной установке “ИСКРА-5”. Сборник трудов научной сессии МИФИ 2009, в печати.
    5. Gus'kov S.Yu., Demchenko N.N., Doskoch I.Ya., Rozanov V.B., Stepanov R.V., Tishkin V.F., Vergunova G.A., Yakhin R.A., Zmitrenko N.V. “Mathematical modeling of the fast ignition laser targets and the HiPER project”. XXIV International Conference on Interaction of Intense Energy Fluxes with Matter. March 1-6, 2009, Elbrus, Kabardino-Balkaria, Russia.
    6. Гуськов С.Ю., Демченко Н.Н., Доскоч И.Я., Змитренко Н.В., Розанов В.Б., Степанов Р.В., Тишкин В.Ф., Яхин Р.А., “Математическое моделирование сжатия лазерных мишеней для энергии лазера 1-2~кДж”. XXIII Международная конференция “Уравнения состояния вещества”. 1-6 марта 2008 г., Кабардино-Балкарская республика, Россия.
    7. V.B. Rozanov, R. Yakhin, N. Zmitrenko et al. “Compression and Thermonuclear Yield of Laser Targets under Conditions of HiPER Facility with due to Regards for Irradiation Asymmetry”. Proc. of the XXX ECLIM, Darmshtadt, 2008, be in print.
    8. V.B. Rozanov, R. Yakhin, N. Zmitrenko et al. “The mixing zone growth laws based on the evolution theory and evolution of neutron yield decreasing”. Proc. of the XXX ECLIM, Darmshtadt, 2008, be in print.
    9. Аристова Е.Н., Розанов В.Б., Змитренко Н.В., Яхин Р.А., Баришпольцев Д.В., Лимпоух И. “Перенос энергии в малоплотных пористых средах легированных тяжёлыми элементами: анализ и моделирование экспериментальных результатов, полученных на установке PALS”. XXXV Международная (Звенигородская) конференция по физике плазмы и управляемому термоядерному синтезу, (г. Звенигород, 11-15 февраля 2008 г.).
    10. Яхин Р.А., Розанов В.Б., Змитренко Н.В., Прончева Н.Г. “Модель перемешивания оболочек термоядерной лазерной мишени при сферическом сжатии”. XXXIV Международная (Звенигородская) конференция по физике плазмы и управляемому термоядерному синтезу, (г. Звенигород, 12-16 февраля 2007 г.).
    11. I. Doskoch, N. Proncheva, V. Rozanov, R. Stepanov, R. Yakhin, N. Zmitrenko. “The Evolution Model of Mixing Zone Growth in the Case of a Spherical Shell Compression”. The 10th IWPCTM, Paris, France 17-21 July 2006. Proceedings, pр. 337-342.

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

Структура и объем диссертации.

Диссертация состоит из введения, четырех глав и заключения. Работа содержит 135 страниц машинописного текста, 42 рисунка, 14 таблиц. Список литературы включает 105 наименований.

Содержание работы.

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

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

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





В главе 2 представлены двумерные расчеты перемешивания двух газов разной плотности (постановка соответствовала экспериментам С.Зайцева на ударных трубах [35]), изучались процессы развития гидродинамических неустойчивостей Рэлея-Тейлора (РТ) и Рихтмайера-Мешкова (РМ).

Проводилось исследование влияния разнообразных начальных условий на развитие турбулентного перемешивания; для плоской геометрии была проведена серия 2D расчетов развития неустойчивостей РТ и РМ при различных условиях по амплитуде и спектру мод начальных возмущений, числам Атвуда. Моделировалась переходная и нелинейная стадии развития гидродинамической неустойчивости. Для однотипных возмущений проведены расчеты при различном случайном выборе фаз возмущений. Геометрия задач – 2D и 3D (около 10%). Результатами расчетов явились поля плотности, компонент скорости, кинетической энергии, энстрофии, размеры зоны перемешивания, масса перемешенного вещества. Большинство расчетов были проведены по численному коду NUT [12] и МАХ [14].

В расчетах рассматривалась следующая плоская задача, соответствующая реальным экспериментам, - в двумерной области {х1 x х2, z1 z z2} c размерами х2-х1 = 72 мм и z2-z1 = 150 мм, z1 = -100 мм, z2 = 50 мм. В начальный момент времени задавались постоянная по пространству температура и соответствующие равновесные распределения плотностей и давления двух газов, находящихся в поле постоянной внешней силы (рис.1).

 Геометрия задачи. В отсутствие-1

Рис.1. Геометрия задачи.

В отсутствие возмущений поверхность раздела газов имела координату z=0, выше нее, в области 0 z z2 располагался тяжелый газ, ниже, в области z1 z 0 - легкий. Ускорение внешней силы мм/мкс2 ( м/с2 – ускорение свободного падения) было направлено вниз.

Для описания термодинамических свойств газов применялось уравнение состояния идеального газа в виде:

, ,

- показатель Пуассона, который для всех газов принимался равным ; - универсальная газовая постоянная; - молярная масса газа (для ксенона - , гелия - , аргона - ).

В состоянии гидростатического равновесия плотность и давление изотермического газа имели экспоненциальное распределение:

pi(z)=poexp(-),

i=oiexp(-),

где индексы i=1, 2 означали соответственно легкий и тяжёлый газ, oi и po- плотности и давление газов на поверхности раздела, mi- масса молекулы газа, Т0=300 0К (температура газов выбрана соответственно условиям эксперимента на ударных трубах), k - постоянная Больцмана. Для двух выбранных пар газов Xе-Аr и Хе-Не на контактной границе: p0=0.5 aтм, 01=8.07·10-4 г/см3 (аргон), o1=8.07·10-5г/см3 (гелий), 02=2.63·10-3 г/см3 (ксенон).

В приведенных исследованиях возмущения задавались в соответствии со следующей формулой для формы контактной поверхности:

z(x) = - a0i cos (kix+i). (1)

Здесь i пробегает 6,8 или 10 значений по числу генерируемых мод возмущений, ki = 2/i - волновое число гармоники с длиной волны i, a0i - начальная амплитуда отдельного возмущения. Фазы i = 2/8·mi, где целое число mi из диапазона 0 mi 7, выбирались случайным образом. Длины волн выбирались из условия отсутствия в наборе гармоник кратных значений, для чего было принято следующее определение:

i=(х2-х1)/ n (2)

где ni - i-тое простое целое число (ni = 2,3,..29), х2-х1 = 72 мм. Тем самым, диапазон значений длин волн определялся как 2.531 i 36 мм.

Рассматривались две принципиально различные задачи. В первом случае изучался процесс перемешивания вследствие развития неустойчивости Рэлея-Тейлора, вызванного внешней постоянной силой. Во втором случае изучалась эволюция неустойчивости Рихтмайера-Мешкова. Для этого в начальный момент через вещество из легкого газа в тяжелый и наоборот вдоль оси z проходила ударная волна с числом Маха 2.5-4.5. В качестве примера на рисунке 2 представлен пример расчета развития неустойчивости Рэлея-Тейлора.

 Эволюция неустойчивости-16

Рис.2. Эволюция неустойчивости Рэлея-Тейлора. Поле плотности (представлено положение границы между тяжёлым и лёгким газом) – результат расчета для пары He-Xe (плоский случай, число гармоник возмущения n=6).

На основе численных расчетов была создана обширная база данных по развитию неустойчивости и перемешивания двух газов в плоской геометрии при постоянном ускорении для режимов, отличающихся типом неустойчивости (РТ и РМ), числом Атвуда (0.941 и 0.532), различным числом учитываемых мод и, соответственно, максимальным номером наиболее коротковолновой моды (n6=13, n8=19 или n10=29), различными начальными амплитудами а0i и их зависимостью от номера моды (а0i =const/k и а0i =const);

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

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

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

- на нелинейной стадии скорость рост пузыря легкого газа практически останавливался, ширина зоны росла за счет проникновения струй тяжелого газа;

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

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

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

В главе 3 на основании анализа проведенных численных расчетов и существующих теоретических моделей описания турбулентного слоя в процессе перемешивания двух разноплотных веществ была разработана теоретическая модель для описания ширины и скорости роста зоны турбулентного перемешивания для широкого диапазона плотностей и ускорений, учитывающая влияние начальных условий. Следует иметь в виду, что значения плотностей и ускорений в экспериментах на ударной трубе и мишенях лазерного сжатия отличаются на много порядков: плотности 10-3 г/см3 и 103 г/см3, ускорения 107 см/с2 и 1015 см/с2 соответственно. Поэтому предложенная нами модель будет опираться на фундаментальные свойства газодинамического подобия.



Pages:   || 2 | 3 |
 

Похожие работы:










 
© 2013 www.dislib.ru - «Авторефераты диссертаций - бесплатно»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.