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

Об учете рефракции лазерных лучей при моделировании двумерно-неоднородного сжатия мишеней Текст научной статьи по специальности «Физика»

CC BY
99
18
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по физике, автор научной работы — А В. Искаков, И Г. Лебо, И В. Попов, В Б. Розанов, В Ф. Тишкин

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — А В. Искаков, И Г. Лебо, И В. Попов, В Б. Розанов, В Ф. Тишкин

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Об учете рефракции лазерных лучей при моделировании двумерно-неоднородного сжатия мишеней»

УДК 533.951

ОБ УЧЕТЕ РЕФРАКЦИИ ЛАЗЕРНЫХ ЛУЧЕЙ ПРИ МОДЕЛИРОВАНИИ ДВУМЕРНО-НЕОДНОРОДНОГО

СЖАТИЯ МИШЕНЕЙ

А. Б. Искаков1, И. Г. Лебо, И. В. Попов1, В. Б. Розанов, В. Ф. Тишкин1

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

В работах [1, 2] исследовалась рефракция лазерного излучения с учетом одномерной гидродинамики мишени с усреднением по телесному углу. Для двумерных лагранжевы.х расчетов сжатия мишеней был создан комплекс программ "АТЛАНТ" [3].

В настоящей работе развита физико-математическая модель и создана подпро грамма расчета распространения лазерных лучей в двумерно-неоднородной лазерной плазме сферических мишеней. Эта подпрограмма используется в комплексе программ "АТЛАНТ" для моделирования нагрева и сжатия термоядерных микромишеней с учетом эффектов рефракции и отражения лазерных лучей. Физическая модель подобна той. что описана в [4]. Однако она применяется к треугольной лаграижевой сетке в сферических осесимметричных координатах (г, (¿>). Коэффициенты поглощения и преломления взяты из работ [5, 6].

1 Институт математического моделирования РАН.

Постановка задачи и основные уравнения. В комплексе программ "А ГЛЛН I рассматривается многослойная лазерная мишень, близкая к сферической. Общим прел-полол<ением двумерного расчета является наличие устойчивой осевой симметрии II ¡»а-метры не зависят от долготы, скорость вещества и поток энергии вдоль соответствую щих координатных линий отсутствуют. На ми1пень падает осесимметричное монохрп ма гическое лазерное излучение (ЛИ) с заданным пространственным и временным ра< пределением мощности. Под действием ЛИ мишень частично испаряется, образуя пла менную корону, которая предполагается полностью ионизованной. Эта задача репы 11 в предположении трехтемпературной гидродинамики с учетом электронной, ионноп и радиационной теплопроводности, обмена энергией между электронной и ионной компонентами и между электронной и радиационной компонентами, объемных потерь на собственное излучение плазмы. В данной работе описывается алгоритм поглощения ЛИ в рамках приближения геометрической оптики (ГО). Более подробное- описание задачи и остальных частей алгоритма можно найти в [3, 7].

На поверхности мишени S (или на выпуклой поверхности, окружающей мишень задаются функции P®{ts), s"(rs), Ts = (R$. Zs) G S, определяющие мощность по к. ка (Вт/м2) и единичный вектор в направлении падающего лазерного излучения, где (Я-Z) радиус и высота в двумерных цилиндрических координатах (либо декартов!. координаты в случае плоской геометрии), г номер луча (лазерные пучки разбиваю i на дискретные лучи, каждый из которых распространяется в двумерно-неоднороднч плазме согласно законам геометрической оптики). Требуется найти распределение м<>т ноет и потока ЛИ в области мишени P{r) - Y1 Д'(г> s«)j гДе функции P,(r. s,) задаются параметрически через г(г, Г5), s, (r, r.s), Р, (т, r.s) системой дифференциальных уравнен и и ГО [4] и граничных условий

% = nst

< = i V г(0) = Ts, s,-(0) = s°(r5), Р,(0) = P?(rs), 0 < г < +оо.

f- =

где г радиус-вектор луча, т.е. траектории, вдоль которой распространяется эле к гро магнитная волна; s, - единичный вектор в направлении г-го луча; Рг - мощность / -го луча лазера; т длина оптического пути вдоль луча; n(r),/x(r) - известные функции, определяющие коэффициенты преломления и поглощения внутри мишени. Тогда объем ная мощность, выделяемая в точке мишени с координатами г, определяется выражент (г(г)Р(г).

При таком подходе используются следующие предположения:

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

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

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

Выражения для коэффиицентов //, и п определяются обратным тормозным эффек том и взяты из [5, 6]. При этом предполагается, что плазма полностью ионизована, электронная компонента квазиравновесна и имеет температуру Г/?. интенсивность ноля мала настолько, что нелинейных эффектов нет, междуэлектронными столкновениями можно пренебречь по сравнению с электронно-ионными (лоренцова плазма). Тогда ком плекспая диэлектрическая проницаемость определяется уравнениями:

V us и; pf ¡у ei е-ег+ ге2; ег = 1 - , е2 = , где г; = ——, л = -,

1 + s I + s и и

где введены обозначения vEi = эффективная частота столкновений эле к

Vm/iT£

тронов с ионами; ш частота лазерного излучения; пе, п/ концентрации электро нов и ионов в плазме; тв масса электрона; т/ = Лтр средняя масса иона: 1

средний атомный вес ионов; е - заряд электрона; Ze - средний заряд ионов: У'/.; электронная температура; A ei кулоновский логарифм; шр/, = плазменная

частота. Для значения кулоновского логарифма используется асимптотическое выра жение при ш >> шрь, но оно является достаточно точным даже при экстраполяции до частот и> ~ LtJpi [8]. На значение кулоновского логарифма влияет также- условие квантового рассмотрения стокновений: при < ^ используется борновское прибли

жение. при > 1 квазиклассическое выражение для сечения рассеяния, т.е. при

Т < 0,077г>[кэВ): А к/ = при Т > 0,077¿*[кэВ): Л Е, = ^

7 = 1,78.

Коэффициенты преломления п и поглощения /i вводятся по определению: у/<

Отсюда п = Ц =

В случае £] >> б2 имеем п = л/^ъ А* — ZJtt-

Дискретизация модели и расчет луча в ячейке. На границе двумерно!! р,и ностной сетки задается дискретное число лучей с заданными начальными условиях Каждый луч последовательно проходит набор ячеек, в которых определяется ero i [>«. ектория и выделенная им энергия. Произвольная выпуклая четырехугольная яченк в декартовом пространстве разбивается на 4 треугольника с общей вершиной в цеп тре ячейки. Поглощение и рефракция луча рассчитываются в полученных треугольных ячейках. Внутри каждой из них, следуя [4], предполагаем: у??.2 = const : у(n/i) - const В этом приближении из уравнений геометрической оптики в параметричес ком виде получаем (здесь /. - путь в плазме) уравнение направления луча:

ns(t) = (ns)0 + ^ V n2t,

уравнение траектории луча:

г(<) = го + (ns)0t + j(V"2)<2, (-)

уравнение поглощения мощности луча:

í _ t2 - - ¿3 1 P{t) = Роехр I -{nfi)0t - (y(n/i), (ns)0)- - (V(n/x), уд2) J2 Í '

уравнение для произведения коэффициентов поглощения и преломления:

пц{1) = (пц)о + (v(n^í), (■ns)o)t + yn2)i2. (4)

В точке входа луча в треугольник предполагаются известными значения параметров (»So), {пц) о, Го, Ро- Из уравнения (2) определяется точка пересечения луча с од поп из граней треугольника. Так как параметр t растет вдоль луча монотонно, то выбирается точка пересечения с минимальным значением <mm. Если прямая, содержащая грань задается точками A(x-i,í/i) и В{х2,у2)-, а входная точка С(х0,Уо), то возможная гочк пересечения грани и луча (х, у) определяется условием

■г2 — Х\ У г ~ V1 1 , 2\ 2 ^ , 2 2

-=-, где а; = х0 + nsxt + -(Ул-га )t , у = у 0 + nsyt. + -(УугГ)г, (•)>

х -xi у - yi 4 4

т.е. решением квадратного уравнения относительно параметра /..

Выбирается минимальный положительный действительный корень этого уравнены», если он вообще существует. Так как х2+у2 ос при I —* оо, то траектория луча, уходич на бесконечность и должна при этом пересечь хотя бы одну грань конечной ячейки (в случае п = 0 и Vп — 0 предполагается, что энергия луча полностью поглощаеч ся). Д. I» каждой из трех граней решается уравнение (5) и выбирается (обязательно существующий) минимальный положительный корень £ТО1П, по которому согласно уравнениям (2.1) (2.1) определяются выходные параметры ж(*т,п), у(1т4п), пз(/т,п), (пр)(1.„пп). /...„,„ I и следующая треугольная ячейка, в которую входит луч.

Проведенные расчеты и сравнения. Для тестирования была рассмотрена задача о падении плоской электромагнитной волны иод углом на слой плазмы с диэлек I ричс ской проницаемостью, зависящей линейно от расстояния. Получено хорошее соогвеч ел вне с аналитическим решением [5]. Выло проведено две серии расчетов для двух ч инов мишеней: I ) гонкое генные оболочечные миптни, 2) мишени типа "лазерный парник

[9, 10].

I. В расчетах исследовалась мишень прямого сжатия, состоящая из ОО-тъг в сте клянной оболочке радиусом 0,25 мм и толщиной 2,4 мкм. Строго в ее ценчр был сфокусирован импульс /V(/-лазера (длина волны 1,06 мкм) треугольной формы (с вершиной в момент 3,6 не) общей длительностью 4,6 и с и мощностью 250 Дж. Угловое распределение мощности импульса задавалось в виде

= 9о(1 + 0,25соз(4<р))Дж/с.- ср.

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

На рис. 1 показана зависимость коэффициента сжатия £ = ъ-Щ^ъ-, где До- Нмах-

IIм/д/ начальное, текущие максимальное и минимальное значения радиуса границы "оболочка горючее", от времени; на рис. 2 - рост относительной амплитуды возмущения: Ь/И = Пмлх-Ншн ох коэффициента сжатия Результаты расчетов показали.

Г*М ЛХ 1 ''МIN

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

8.00 п

5/Я

2.00

4.00

0.00

6.00

0.00

/

0.00 1.00 2.00 3.00 4.00 5.00 6.00 0.00

2.00 4.00

6.00

8.00

I, не

Рис. I. Зависимость коэффициента сжатия мишени (£) от времени (/) для расчетов сжатия оболочечной мшиепи с учетом (1) и без учета (2) рефракции.

Рис. 2. Зависимость относительной амплитуды возмущения 6/11 контактной границы "оболочка горючее" от коэффициента сжатия £ для случаев (1) и (2).

2. С целью увеличения равномерности сжатия в работе [9] была предложена но вая конструкция мишени типа "лазерный парник". Такая мишень представляет собой внешнюю камеру из вещества с большим атомным весом, имеющую отверстия д и ввода лазерного излучения. Внутри камеры содержится малоплотный поглотитель. центре мишени располагается рабочая мишень с ОТ-горючим. В мишени такого ти па энергия лазерного излучения переходит в тепловую энергию плазмы с малым /. а затем капсула с .ОТ-горючим сжимается за счет теплового давления. В [10] на основа нии двумерных численных расчетов, выполненных с помощью программы "АТЛАНТ без учета рефракции лазерных лучей, было показано, что с помощью теплопроводного выравнивания в малоплотном поглотителе удается подавить лишь мелкомасштабны! возмущения. Неоднородность нагрева, связанная с малым количеством лазерных пуч ков, приводит к сильному отклонению формы мишени от сферической. В настоящей статье мы привели результаты расчетов нагрева такого типа мишеней по программе "АТЛАНТ" с учетом рефракции лазерных лучей. Выбранная для расчета мишень сое го-яла из ОТ-газа в полимерной оболочке (СН) (/? = 0,35 мм), малоплотного поглотителя

(Я, = 0,60 м,м) и Си-оболочки с двумя отверстиями. Из центра мишени отверстия бы ли видны под углом в 60°. В отверстия симметрично с двух сторон вводился лазерный импульс (длина волны 0,353 мкм) треугольной временной формы (вершина в момен I 1,3 не, общая длительность 2,6 не) с энергией 6 кДж с угловым распределением

{ 0,|у>1 > 7Г/6

На рис. 3 показана форма мишени к моменту 1,05 не (ось симметрии расположена горизонтально).

Внутренняя I = 1.48 не

оболочка

Рис. 3. Форма мишени типа "лазерный парник" на момент времени ¿=1,05 не для случая лучей, сфокусированных в центр мишени. На рисунке показана лагранжевая сетка мтшни. Жирные линии - лазерные лучи.

Рис. 4. Форма мишени типа "лазерный парник" на момент времени ¿=1,48 не для случая распространения лазерных лучей параллельно оси.

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

этому моменту нейтронный выход составил 7,6 • 109. В момент 1,5 1,6 нс рабочая мишень разрушается.

Целью настоящей статьи было дать описание модели учета рефракции шзерпых iy чей в двумерно-неоднородной сферической мишени и продемонстрировать возможное : программы "АТЛАНТ". Были проведены двумерные расчеты нагрева и сжатия двух типов мишеней и показано, что учет рефракции приводит к более симметричному на. г реву короны мишеней прямого сжатия. В случае мишеней типа "лазерный парник' из-за сильной неодномерности засветки мишени (моделировалось облучение сферич< ской мишени двумя пучками) мишень разрушилась до окончания лазерного импульса. Видно, что в случае, когда лазерные лучи распространяются параллельно оси. мишень сжимается более равномерно, чем в случае фокусировки лучей строго по радиусу. К моменту разрушения мишени нейтронный выход составил 7,6 • 109. Работа поддержана МНТЦ (грант N 029-94).

ЛИТЕРАТУРА

[1] Афанасьев Ю. В., Гамалий Е. Г., Демченко Н. Н.. Розанов В. Б. Труды ФИАН, 134, 32 (1982).

[2] Р о з а и о в В. Б., Демченко Н. Н. Квантовая электроника, 12. N 9, 1895 (1985).

[3] Гамалий Е. Г., Демченко Н. Н., Л е б о И. Г. и др. Квантовая электроника, 15, 1622 (1988). См. также Волкова Р. А., Т и ш к и н В. Ф. и др. Сафра. Функциональное наполнение. Программа АТЛАН'1 . Решение двумерных задач управляемого термоядерного синтеза. Инструкция. НИМ им.

М. В. Келдыша, М., 1985.

[4] Б о л ь ш о в Л. А. и др. Численное моделирование распространения лазерного излучения в неоднородной поглощающей среде, ИАЭ-4732/16, ЦНИИатоминформ. 1988.

[о] Гинзбург В. Л. Распространение электромагнитных волн в плазме, Наука. М., 1967.

[6] Л и ф ш и ц Е. М., Питаевский Л. П., Теоретическая физика. 10. Физическая кинетика, 1979.

[7] Л е б о И. Г., Попов И. В., Розанов В. Б., Т и ш к и н В. Ф. Journal of Russian Laser Research, 15, 136 (1994).

[8] Б е к е ф и Дж. Радиационные процессы в плазме, Мир, М., 1971.

[9] Г у с ь к о в С. К)., 3 м и т р е н к о Н. В., Розанов В. Б. ЖЭТФ. 108. 548 (1995).

[10] JI е б о И. Г., Попов И. В., Розанов В. Б., Т и ш к и н В. Ф. Квантовая электроника, 22, 1257 (1995).

Поступила в редакцию 18 ноября 1996 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.