Численное моделирование взаимодействия фемтосекундных лазерных импульсов с металлами тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

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

Шепелев Вадим Владимирович

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВЗАИМОДЕЙСТВИЯ ФЕМТОСЕКУНДНЫХ ЛАЗЕРНЫХ ИМПУЛЬСОВ С МЕТАЛЛАМИ

Специальность 01.02.05 - Механика жидкости, газа и плазмы

Автореферат

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

1 2 ДПР ¿012

Москва-2012

005019886

Работа выполнена в Институте Автоматизации Проектирования РАН

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

Олег Валентинович Трошкин

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

Владимир Иванович Мажукин

кандидат физико-математических наук, Андрей Викторович Конюхов

Ведущая организация: Институт теоретической физики им. Л. Д. Ландау РАЛ

Защита состоится «_»___ 2012 г. в_часов на заседании

Диссертационного совета Д 002.024.03 при Институте Прикладной математики им. М.В. Келдыша РАН по адресу: 125047, г. Москва, Миусская пл. д 4.

С диссертацией можно ознакомиться в библиотеке Института Прикладной математики им. М.В. Келдыша РАН.

Автореферат разослан « »_2012 г.

Ученый секретарь диссертационного совета Д 002.024.03

доктор физико-математических наук •*/ Н.В. Змитренко

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

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

1.1. Актуальность темы

В диссертационной работе рассматриваются ультракороткие лазерные импульсы (УКЛИ), порядок длительности которых находится в промежутке от единиц до сотен фемтосекунд (1 фс = 10"с). Лазеры с фемтосекундными импульсами появились около 30 лет назад, и они имеют большое количество приложений в науке и технике.

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

1.2. Цели и задачи работы

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

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

1.3. Научная новизна

Научная новизна настоящей работы состоит в следующем:

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

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

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

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

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

5. В представленной модели используется современная модель широкодиапазонной двухтемпературной электронной теплопроводности, хорошо согласующаяся с экспериментом и известными моделями теплопроводности.

6. Выполнены расчета в широком диапазоне температур: от порога плавления металла до плазменных температур в десятки эВ.

1.4. Научная и практическая ценность

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

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

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

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

1.5. Личный вклад автора

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

1.6. Апробация работы

Материалы настоящей работы докладывались на семинарах и конференциях:

• Семинар ИММ РАН 20.10.2009 г.;

• Japan-Russia Symposium on "Numerical Experiment in Hydrodynamical Instability and Turbulence with High-Performance Computing", Москва, 2009 г.

• 7th International Seminar on Mathematical Models and Modeling in Laser-Plasma Processes.

1.7. Объем и структура диссерт ации

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

1.8. Публикации

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

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

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

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

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

—+ — .0(2)

д! 8т

Ё^+р^^фЩ-^Ъ-Т,)*^ (3)

Э/ дгп дт\ 8т) р

гъ+рЬ'Ьсг.-т,)«)

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

Р,,РС,Р - ионное, электронное и полное давление, /* = /; +1\;

,г.г,С - ионная, электронная и полная энергия единицы массы, с = с, + с,; Т„Те - ионная и электронная температуры.

X

т - лагранжева массовая координата, йт = рсЬ., т = ; г - время;

к -коэффициент широкодиапазонной электронной теплопроводности в двухфазной двухтемпературной среде [3];

ас, =36-1016 •——^Ж - коэффициент электронно-ионной релаксации алюминия, ра с-К-м

значение взято из работы Ю.В. Петрова [4] и согласуется с оценками Ж. Лина и Л.В. Жигилея [5];

р0 = 2700 -К-- - плотность алюминия при нормальных условиях. м

р

ехр|--—-- | - энергия лазерного излучения,

поглощенная единицей лагранжевой координаты за единицу времени;

РаЬ, ~ энергия лазерного излучения, поглощенная единицей облученной поверхности

(поглощенный флюенс);

<5 = 10(ш - тол1цина скин-слоя;

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

Электронная теплопроводность. В работе используется широкодиапазонная модель электронной теплопроводности алюминия [3]. В разделе приводится ее краткое описание и сравнение с другими известными моделями. Используемая модель отдельно учитывает электрон-электронные и электрон-электронные столкновения для плазмы и конденсированной среды. Отдельно рассчитывается также теплопроводность твердой и жидкой фаз алюминия.

Уравнения состояния. В расчетах используются отдельные уравнения состояния для электронной и ионной подсистем металла. Для ионной подсистемы берется табличное полуэмпирическое многофазное уравнение состояния алюминия [6 - 8], которое включает в себя области фазовой диаграммы со смешанным существованием фаз и области метастабильных состояний алюминия. Для энергии и давления электронной подсистемы используется интерполяции между значениями из уравнений состояния классического идеального и вырожденного фермиевского электронного газа [9].

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

Г0,,=30(Ж,р0=2700^,у0=0.

см

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

справа ставятся граничные условия на давление Р1 = Р„ = 0. а также условия на нулевой дТ,

поток тепла через границы:

dm

97^ dm

= 0.

№1

№2

__.................__]

Рис. 1 Три конфигурации облучаемых образцов для трех различных классов экспериментов: 1. «Вакуум-алюминий-вакуум». 2. «Вакуум-алюминий-стекло». 3. «Стекло-алюминий-вакуум».

Во втором случае «вакуум-алюминий-стекло» рассматривается задача о расширении в вакуум алюминиевой пленки, напыленной на стеклянную подложку. Наблюдение методом pump-probe ведется за задней поверхност ью пленки, то есть, за границей между алюминием и стеклом. Вычислительная область разбивается на две подобласти: алюминий и стекло. Для стекла используется то же уравнение состояния, что и для алюминия, отражением от границы «стекло- алюминий» пренебрегаем из-за близости акустических характеристик стекла и алюминия. Отличие расчетов в подобластях с алюминием и стеклом состоит в том, что в стекле отсутствует электронная теплопроводность. Толщина алюминия dM = 400 нм, толщина стекла dgtal = 400 им. Как и в первой постановке, слева и справа ставятся граничные условия на давление PL= PR- 0. Теплопроводность в стекле отсутствует, поэтому справа условие на нулевой

поток тепла ставим на границе «алюминий-стекло»:

81 dm

ат;

am

0, где

ЭТ.

Al-gh,

- градиент температуры через границу между алюминием и стеклом.

В третьем случае «стекло-алюминий-вакуум» лазерный луч перед тем, как попасть на пленку, проходит сквозь стекло. Наблюдение методом pump-probe ведется за правой границей. Расчеты производились для конфигураций в диапазоне толщин алюминиевой пленки dv от 350 до 1200 нм. Толщина стеклянной подложки на практике составляет

около 150 микрон. Для численного моделирования достаточно взять dghss в 2-3 раза больше, чем dtl, чтобы за интересующий нас в численном эксперименте промежуток времени возмущение, вызванное поглощением импульса на границе «стекло-алюминий»,

9

не успело отразиться от другой границы стекла и вернуться. Слева и справа ставятся граничные условия на давление Р, =ГК= 0, а условия на нулевой поток тепла через

57'

дТ.

границы выглядят так:

8т

л ,-цш 5т

- О. Из-за отсутствия теплопроводности в

я

стекле левое граничное условие ставится на границе «алюминий-стекло».

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

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

Этап гидродинамики. На первом вычислительном этапе решаются уравнения гидродинамики с учетом поглощения металлом лазерного излучения. Для решения уравнений первого этапа используется неявная конечно-разностная полностью консервативная схема Самарского [10], адаптированная для двухтемпературного случая. Порядок аппроксимации схемы второй по координате и первый по времени

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

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

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

В третьей главе представлены результаты работы.

«Сверхзвуковое плавление» в алюминии. Рассмотрим экспериментальную

конфигурацию «стекло-алюминий-вакуум». Лазерное излучение проходит сквозь стекло и

поглощается в поверхностном слое пленки около границы «стекло-алюминий». Весь

происходящий процесс можно условно разделить на две стадии: двухтемпературную (2Т)

и следующую за ней однотемпературную (1Т). Лазерное излучение поглощается только

электронами и не нагревает ионы, поэтому поглощение ультракороткого лазерного

импульса металлом приводит к образованию двухтемпературного (2Т) слоя,

существующего в течение нескольких пикосекунд. Ширина двухтемпературного слоя

составляет несколько десятков нанометров. В начале электронная температура в 2Т-слое

10

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

Состояние металла в 2Т-слое в конце 2Т-стадии характеризуется большим давлением (пик превышает 10 ГПа) и плотностью, примерно равной плотности металла при комнатной температуре. Это сильное возмущение давления, сформированное на 2Т-стадии, впоследствии распространяется вглубь пленки. О формировании этого возмущения говорят как о «нагружении» 2Т-слоя. Распространение возмущения вглубь металла или «разгрузка» 2Т-слоя определяет все дальнейшие гидродинамические процессы, происходящие в пленке. Помимо электронно-ионной релаксации, существенное влияние на перенос энергии на 2Т-стадии оказывает электронная теплопроводность, которая приводит к распространению тепловой волны вглубь металла. Поскольку ионы успевают нагреться до температур, превышающих температуру плавления, фрмгг плавления движется вглубь металла, следуя за тепловой волной. Скорость тепловой волны значительно превышает звуковую (порядок звуковой скорости в алюминии 6000 м/с).

Поэтому об этом эффекте говорят как о неравновесном «сверхзвуковом плавлении» (supersonic melting, [11, 12]) или гомогенном плавлении. Плавление происходит за счет образования зародышей жидкой фазы в двухфазной области внутри 2Т-слоя. Быстрое продвижение фронта плавления возможно за счет большей температуропроводности электронов по сравнению с ионами, в работе приводятся соответствующие оценки. На рис. 2 показаны траектории движения фронтов плавления в пленке от границы «стекло-алюминий» (х = 0), полученные с помощью настоящей модели (кривые 1 и 2, Fab = 1300 Дж /мг) и полученные с помощью прямого молекулярно-динамического моделирования [12] (кривая 3, Fabl = 820 Дж!м2).

О 10 20 30 40

Рис. 2 Движение фронтов плавления внутри пленки. /,. время релаксации 2Т-слоя, штриховой линией показана звуковая характеристика.

Произвести исследование «неравновесного» эффекта «равновесными» методами прямого численного моделирования без выделения фронтов плавления возможно здесь благодаря использованию широкодиапазонного табличного уравнения состояния [6-8].

Исследование профиля волны сжатия в алюминии. Основным итогом 2Т-стадии является формирование в граничном слое алюминиевой пленки сложного профиля сжатия. Существенное влияние на форму и распространение этого профиля оказывает процесс плавления металла. На рис. 3-4 показана эволюция профилей давления в пленке. На 2Т-стадии процесса на профиле давления появляются изломы, показанные стрелками на рис. 3. Изломы появляются в связи с наличием зоны плавления в веществе и расположены на границах этой зоны: левый излом соответствует значению относительной массовой доли жидкой фазы £ = 1 (полностью жидкая фаза), правый излом соответствует значению £ = О (полностью твердая фаза). На рис. 5 показана изохора алюминия при нормальной плотности. Происходит «впечатывание» изломов на изохоре плавления в профиль волны сжатия.

40

*,пш

Рис. 3 Профили ионного давления, ионной и электронной температур при значениях времени 1 пс, 2 пс, 5 пс.

Стрелкой «1-2» на рис. 3 показан сдвиг зоны плавления за время с 1 пс до 2 пс, скорость этого сдвига превышает звуковую более чем в 4 раза. Далее скорость падает и становится близкой к звуковой. На ГГ-стадии профиль с изломами обгоняет зону плавления, и распространяется в твердом алюминии. Наклон линии между изломами постепенно становится круче и, в конце концов, опрокидывается, приводя к образованию ударной волны в металле.

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

ю

X, пш

Рис. 4 Профили ионного давления при значениях времени 10 пс и 30 пс.

(9 8

О

О 1000 М00 3000 4000

т,к

Рис. 5 Изломы, ограничивающие зону плавления на изохоре алюминия.

i. Предвестник перед нижней границей акустического «слсда» от зоны

плавления. Он относится к «хвосту» тепловой волны, проникшей в металл на 2Т-стадии, и постепенно поглощается догоняющим следом зоны плавления. . ii. След зоны плавления, граничными точками которой являются изломы (см. рис. 3). Наклон этого практически прямого участка профиля становится все круче, затем происходит его опрокидывание с образованием ударной волны.

iii. Максимум на границе с зоной ii - след верхнего излома на изохоре плавления (рис. 5) и расплав. Второй максимум, по-видимому, имеет численное происхождение и связан с осцилляциями за участком ii.

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

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

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

отражении от задней стенки пленки соответствующей ему по номеру части сложного профиля давления на рис.4.

i. Предвестник. Это след тепловой волны, которая возникает па 2Т-стадии. Экспериментально этот этап практически неразличим, для его исследования нужны более тонкие гшенки.

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

iii. Когда этап іі заканчивается, от задней стенки оіражается следующая за ударной волной волна разрежения.

¡V. Скорость задней границы падает на этапе отражения волны разрежения.

V. На этом этапе скорость становится практически постоянной из-за наличия плато на профиле давления на рис. 4. Плато образуется из-за расширения расплава.

80 —........................................,----

ЛХ«».Ц«{І), пт |

АІ * дЫ» і ео ~га>і> -<зо пи/ыкпает І —г;-

40 Т! »..І^Яі--уі.--------

2ц---:--------------,:лг£.. 1.А1.1.1.

I .я ОеЬу Бта

I :' Шз п1ах(Я1ит о; !Ье:рип!1р, рз

О <--—--'-'-1-

0 « 80 В9 *И 200

Рис. 6 Траектория задней границы пленки. Сплошной кривой показана расчетная траектория; точками, соединенными ломаной - экспериментальные данные.

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

ДХ«».й«{і), піт

AI * glass

"ТаЪа = {зи mo/ciir dgtBjvSSOten

\ т..іа>!*»___

"^-ЗГ-^Эб1

-Л,

--Qlai

Oelay ßme tx^th« maxi:

¿pump,, p-s

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

ПЙА Н}

1 \ \ \ __________А

¡1 V \

»у I / е-Бол* Л» я 2*:ПХ2!П С* с е ригтр, р4

а <3 80 (20 КО 200

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

Исследование абляции алюминия. При облучении металла лазерным импульсом субпикосекундного диапазона или более коротким импульсом возникает особый режим лазерной абляции, так называемая откольная абляция. Сам по себе процесс лазерной абляции, т.е., уноса, испарения вещества с поверхности образца благодаря действию лазера, происходит при любой мощности входного импульса. Специфика откольного режима по сравнению с обычным заключается в двух вещах. Во-первых, при отколе уносится гораздо больше вещества: в эксперименте при переходе в соответствующий диапазон поглощенной энергии масса унесенного вещества скачком возрастает на несколько порядков. А во-вторых, при откольной абляции в определенный момент времени от поверхности мишени отделяются не маленькие порции испаренного вещества, а целое жидкое образование, называемое откольной пластиной. Явление откольной абляции было впервые открыто в оптических экспериментах с металлами и диэлектриками, описанных в [13], а впоследствии детально описано и объяснено в [14], [15] и [8]. Откол может существовать при условии на поглощенную энергию > где величина Ра называется порогом абляции. При более низких энергиях абляция происходит в обычном, безоткольном режиме. Основной физической причиной откольной абляции является кавитация в расплавленной части алюминиевой пленки. Кавитацию, в свою очередь, вызывают растягивающие напряжения в пленке. Вещество в зоне будущего

откола находится в двухфазном состоянии «жидкость + газ», на фазовой диаграмме эти

16

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

-ЮС 0 100 200 300 400 500 x,nm

Рис.8 Профиль давления при откольной абляции алюминиевой пленки в момент

откола.

-100 0 100 200 300 «0 500

К, ГШ!

Рнс.9 Профиль плотности при откольной абляции алюминиевой пленки в момент

откола.

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

17

с меньшим значением флюенса откол не наблюдается. На рис. 8-9 показаны профили давления и плотности для абляции в отколыюм режиме с флюенсом, равным пороговому. Откол происходит в момент времени 39 пс, глубина кратера составляет 33 нм. Согласно [16], молекулярно-динамическое моделирование дает пороговый флюенс 600-650 Дж/м2 в зависимости от выбранного потенциала взаимодействия и глубину кратера 48 геи, а эксперимент - флюенс 640 - 700 Дж1мг и глубину кратера 50 нм, то есть, согласование удовлетворительное. Довольно большая погрешность в пороговом флюенсе и времени вызвана используемым в модели критерием обнаружения откола. В настоящей модели используется «чистая» двухтемпературная гидродинамика без каких-либо специальных добавок, поэтому «откол» на рис. 8-9 регистрируется в момент, когда под действием растягивающих напряжений вещество попадает в область абсолютно неустойчивых состояний, в спинодаль или ниже спинодали. В задачах по отколу модель может быть использована для оценочных расчетов. Более точное и аккуратное исследование процесса отеола требует учета образования в жидкости зародышей паровой фазы. Такой подход используется, например, в работе [8], где с помощью феноменологической модели учитывается образование и рост пузырьков.

Следующий численный эксперимент по исследованию абляции алюминия проводился по материалах! экспериментов Эванса [17], где поглощенный флюенс

существенно превышал порог абляции: F. = 25610 . Справа от алюминиевой

м

пленки располагалась стеклянная подложка, и экспериментально велось наблюдение за задней (соседней со стеклом) границей пленки с помощью метода pump-probe. Измерялся сдвиг фазы в луче, отраженном от граничной точки. Толщина пленки составляла 400 нм. На рис. 10 показана траектория границы пленки со стеклом, рассчитанная с помощью настоящей модели, в сравнении с экспериментальной кривой. Значения сдвига границы пленки со стеклом были получены из данных эксперимента [17] с помощью известной X ¡jj)

формулы: h= рп—, где Л„гок =800 нм - длинаprobe-импульса, i// - сдвигфазы, Nni =1.45332 — показатель преломления стекла.

Рис.10 Рассчитанная с помощью вычислительного эксперимента (плавная линия) и

экспериментальная (точки, соединенные ломаной) траектории граничной точки для пленки толщиной 400 нм.

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

Заключение

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

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

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

4. В представленной модели используется многофазное широкодиапазонное полуэмпирическое уравнение состояния, что позволяет учитывать в численном

19

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

5. Выполнены расчеты в широком диапазоне температур: от температуры плавления

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

Библиография

1. С.И. Анисимов, Б.Л. Капелиович, Т.Л. Перельман, Электронная эмиссия с поверхности металлов под действием ультракоротких лазерных импульсов, ЖЭТФ, том 66, вып. 2 (1974), с. 776-781.

2. М.И. Каганов, И.М. Лифшнц, Л.В. Танатаров, Релаксация между электронами и решеткой, ЖЭТФ, том 31, выи. 2(8) (1956), с. 232-237.

3. H.A. Иногамов, Ю.В. Петров, Теплопроводность металлов с горячими электронами, ЖЭТФ, т. 137, вып. 3 (2010), с. 505-529.

4. Yu. V. Petrov, Energy exchange between the lattice and electrons in a metal under femtosecond laser irradiation, Laser and Particle Beams (2005), 23,283-289.

5. Zh. Lin, L.V. Zhigilei, and V. Celli, Phys. Rev. В 77,075133 (2008).

6. Бушман A.B., Ломоносов И.В., Фортов В.Е., Уравнения состояния металлов при высоких плотностях энергии. Черноголовка: ИХФЧ РАН, 1992.

7. М.Е. Поварницын, Т.Е. Итина, П.Р. Левашов, К.В. Хищенко "Моделирование абляции металлических мишеней фемтосекундными лазерными импульсами"// Сборник трудов "Физика экстремальных состояний вещества 2007 Черноголовка: ИПХФ РАН, С. 16 (2007).

8. М.Е. Povarnitsyn, Т.Е. Itina, M. Sentis, K.V. Khishchenko, and P.R. Levashov, Material decomposition mechanisms in femtosecond laser interactions with metals, Phys. Rev. В 75,235414(2007).

9. Ландау Л.Д., Лифшнц E.M., Теоретическая физика: учеб. пособие в 10 т., T. V, Статистическая физика, Москва, Наука, 1988 г.

10. Самарский Л.Л., Попов Ю.П., Разностные методы решения задач газовой динамики, Наука, Москва, 1980.

11. N.A. Inogamov, V.V. Zhakhovsky, S.I. Ashitkov, V.A. Khokhlov, V.V. Shepelev, P.S. Komarov, A.Y. Ovchinnikov, D.S. Sitnikov, Yu.V. Petrov, M.B. Agranat, S.I. Anisimov, and V.E. Fortov, Laser Acoustic Probing of Two-Temperature Zone Created by Femtosecond Pulse, Contrib. Plasma Phys. 50, No.l, 151-158 (2010).

12. N. A. Inogamov, S. I. Ashitkov, V. V. Zhakhovsky, V. V. Shepelev, V. A. Khokhlov, P.S. Komarov, M. B. Agranat, S. I. Anisimov and V. E. Fortov, Acoustic probing of two-temperature relaxation initiated by action of ultrashort laser pulse, Applied Physics A, 2010, DOI: 10.1007/s00339-010-5765-2.

13. Sokolowski-Tinten K., Bialkowski J., Cavalleri A., Von der Linde D., Oparin A., Meyer-ter-Vehn J., Anisimov S.I., Transient States of Matter during Short Pulse Laser Ablation, Phys. Rev. Lett.,1998, V. 81. P. 224-227.

14. H.A. Иногамов, A.M. Опарин, Ю.В. Петров, H.B. Шапошников, С.И. Анисимов,

Д. фон дер Линде, Ю. Майер-тер-Фен, Разлет вещества, нагретого ультракоротким лазерным импульсом, Письма в ЖЭТФ, том 69, вып. 4, стр. 284-289.

15. S.I. Anisimov, N.A. Inogamov, Y.V. Petrov, V.A. Khokhlov, V.V. Zhakhovskii, K. Nishihara, M.B. Agranat, S.I. Ashitkov, P.S. Komarov, Tresholds for front-side ablation and rear-side spallation of metal foil irradiated by femtosecond laser pulse, Appl. Phys. A. (2008), v.92, pp. 797-801.

16. V.V. Zhakhovski, N.A.Inogamov, Yu.V. Petrov, S.I. Ashitkov, K. Nishihara, Molecular dynamics simulation of femtosecond ablation and spallation with different interatomic potentials, Appl. Surf. Sci. (2009), DOI: 10.1016/j.apsusc.2009.04.082.

17. R. Evans, A.D. Badger, F. Fallies, M. Mahdieh, T.A. Hall, P. Audebert, J.-P. Gcindrc, J.-C. Gauthier, A. Mysyrowicz, G. Grillon, A. Antonetti, Time- and Space-Resolved Optical Probing of Femtosecond-Laser-Driven Shock Waves in Aluminum, Phys. Rev. Lett. (1996), V. 77, N. 16, P.3359-3362.

Публикации по теме диссертации

1. N.A. Inogamov, A.Ya Faenov, V.V. Zhakhovsky, T.A. Pikuz, I.Yu. Skobelev,

Yu.V. Petrov, V.A. Khokhlov, V.V. Shepelev, S.I. Anisimov, V.E. Fortov, Y. Fukuda, M. Kando, T. Kawachi, M. Nagasono, H. Ohashi, M. Yabashi, K. Tono, Y. Senda,

Т. Togashi, and Т. Ishikawa, Two-Temperature Warm Dense Matter Produced by Ultrashort Extreme Vacuum Ultraviolet-Free Electron Laser (EUV-FEL) Pulse, Contrib. Plasma Phys. 51, No. 5,419 - 426 (2011).

2. Н А. Иногамов, B.B. Жаховский, B.A. Хохлов, ВВ. Шепелев, Сверхупругость и распространение ударных волн в кристаллах, Письма в ЖЭТФ, т. 93, вып. 4, с. 245-251,(2011).

3. Агранат М.Б., Аниеимов С.И., Ашитков С И., Жаховский В.В., Иногамов Н.А., Комаров П.С., Овчинников А.В., Фортов В.Е., Хохлов В.А., Шепелев В.В., Прочностные свойства расплава алюминия в условиях экстремально высоких темпов растяжения при воздействии фемтосекундных лазерных импульсов, Письма в ЖЭТФ, том 91, вып. 9, с. 517-523, (2010).

4. N.A. Inogamov, V.V. Zhakhovsky, S.I. Ashitkov, V.A. Khokhlov, V.V. Shepelev, P.S. Komarov, A.V. Ovchinnikov, D.S. Sitnikov, Yu.V. Petrov, M.B. Agranat, S.I. Anisimov, and V.E. Fortov, Laser Acoustic Probing of Two-Temperature Zone Created by Femtosecond Pulse, Contrib. Plasma Phys. 50, No.l, 151-158 (2010).

5. N. A. Inogamov, S. I. Ashitkov, V. V. Zhakhovsky, V. V. Shepelev, V. A. Khokhlov, P. S. Komarov, M. B. Agranat, S. I. Anisimov and V. E. Fortov, Acoustic probing of two-temperature relaxation initiated by action of ultrashort laser pulse, Applied Physics A, 2010, DOI: 10.1007/s00339-010-5765-2.

6. N.A. Inogamov, V.V. Zhakhovsky, A.Y. Faenov, V.A. Khokhlov,

V.V. Shepelev,I.Y. Skobelev, Y. Kato, M. Tanaka, T.A. Pikuz, M. Kishimoto, M. Ishino, M. Nishikino, Y. Fukuda, S.V. Bulanov, T. Kawachi, Y.V. Petrov, S.I. Anisimov, V.E. Fortov, Spallative ablation of dielectrics by X-ray laser, Applied Physics A, 2010, DOI: 10.1007/s00339-010-5764-3.

7. N.A. Inogamov, V.V. Zhakhovskii, A.Ya. Faenov, V.E. Fortov, Y. Kato, V.V. Shepelev, Y. Fukuda, M. Kishimoto, M. Nagasono, N. Ohashi, Y. Senda, M. Yabashi, K. Tano, M. Tanaka, I.Yu. Skobelev, T.A. Pikuz, S.I. Anisimov, Yu.V. Petrov, M. Ishino,

T. Ishikawa, T. Togashi, M. Nishikino, V.A. Khokhlov, M. Kando, and T. Kawachi, Ablation of insulators under the action of short pulses of X-ray plasma lasers and free-electron lasers, Journal of Optical Technology, Vol. 78, Issue 8, pp. 473-480 (2011).

8. S.I. Anisimov, N.A. Inogamov, V.A. Khokhlov, Yu.V. Petrov, A.Ya. Faenov,

V.V. Zhakhovsky, I.Yu. Skobelev, T.A. Pikuz, V.E. Fortov, M. Tanaka, M. Kishimoto, M. Ishino, M. Nishikino, Y. Fukuda, S.V. Bulanov, T. Kawachi, K. Nishihara, Y. Kato, V.V. Shepelev, Spallative ablation of metals and dielectrics by ultrashort X-ray lasers,

Book of abstracts, International Symposium of High Power Laser Ablation, Eldorado Hotel, Santa Fe, 19 April -22 April 2010, p. 57.

9. S.I. Anisimov, N.A. Inogamov, V.A. Khokhlov, Yu.V. Petrov, V. E. Fortov,

M.B. Agranat, V.V. Zhakhovsky, S.I. Ashitkov, V.P. Komarov, V.V. Shepelev, Pumpprobe Exploration of Ultrafast Electron Processes in Metal Irradiated by a Femtosecond Laser Pulse, Book of abstracts, International Symposium of High Power Laser Ablation, Eldorado Hotel, Santa Fc, 19 April - 22 April 2010, p. 58.

10. S.I. Anisimov, N.A. Inogamov, V.A. Khokhlov, Yu.V. Petrov, V.E. Fortov,

M.B. Agranat, V.V. Zhakhovskii, S.I. Ashitkov, V.P. Komarov, V.V. Shepelev, Acoustic probing of two-temperature relaxation initiated by action of ultrashort laser pulse, COLA 2009,10th International Conference on Laser Ablation, Singapore, 22-27 November 2009, Chairs: Boris Luk'yanchuk, Hong MingHui, Co-Chairs (International): Richard F. Haglund, Mikio Takai, Alfred Vogel, Book of Abstracts, p. 122.

11. S.I. Anisimov, N.A. Inogamov, V.A. Khokhlov, Yu.V. Petrov, A.Ya. Faenov,

V.V. Zhakhovskii, I.Yu. Skobelev, T.A. Pikuz, V.E. Fortov, M. Tanaka, M. Kishimoto, M. Ishino, M. Nishikino, Y. Fukuda, S.V. Bulanov, T. Kawachi, K. Nishihara, Y. Kato, V.V. Shepelev, Spallative ablation of dielectrics by X-ray laser, COLA 2009,10th International Conference on Laser Ablation, Singapore, 22-27 November 2009, Chairs: Boris Luk'yanchuk, Hong MingHui, Co-Chairs (International): Richard F. Haglund, Mikio Takai, Alfred Vogel, Book of Abstracts, p. 174.

12. S.I. Anisimov, N.A. Inogamov, V.A. Khokhlov, Yu.V. Petrov, V.E. Fortov,

M.B. Agranat, V.V. Zhakhovskii, S.I. Ashitkov, V.P. Komarov, V.V. Shepelev. Pumpprobe «seismology» of ultrafast two-temperature processes inside a film irradiated by a femtosecond laser pulse, XIII International Conference on Physics of Non-Ideal Plasmas, September 13-18,2009, Chemogolovka, Russia, Books of Abstracts, p.76-77.

13. InogamovN.A., Anisimov S.I., Petrov Yu.V., Khokhlov V.A., Zhakhovskii V.V., Agranat M.B., Ashitkov S.I., Sitnikov D.S., Ovchinnikov A.V., Nishihara K., Oparin A.M., Shepelev V.V., Femtosecond ablation (FSA): calculated and measured ablation thresholds agree for a first time, in: 8-th International workshop on Subsecond Thermophysics, Moscow, Russia, September 26-2S, 2007, Book of Abstracts, p. 30-31.

Подписано в печать:

27.03.2012

Заказ № 6891 Тираж - 90 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

Введение

1 Модель физических процессов

1.1 Постановка задачи.

1.2 Электронная теплопроводность.

1.3 Уравнения состояния.

1.4 Начальные и граничные условия.

2 Вычислительный алгоритм

2.1 Метод расщепления.

2.2 Этап гидродинамики

2.3 Этап электронной теплопроводности.

2.4 Этап электронно-ионного обмена

2.5 Замечания по реализации.

3 Результаты

3.1 «Сверхзвуковое плавление» в алюминии.

3.2 Исследование профиля волны сжатия в алюминии.

3.3 Исследование абляции алюминия.

Настоящая работа посвящена численному исследованию физических процессов, происходящих при облучении металлов ультракороткими лазерными импульсами (У К ЛИ). Ультракороткими считаются, вообще говоря, все лазерные импульсы длительностью порядка или менее одной пикосекунды (1 пс = Ю-12 с). В работе рассматриваются в основном лазерные импульсы, порядок длительности которых находится в промежутке от единиц до сотен фемтосекунд (1 фс = Ю-15 с). Такие импульсы также называются субпикосекундными и фемтосекундными. Существуют экспериментальные установки с аттосекундными импульсами (1 ас = Ю-18 с), но они еще не получили достаточно широкого распространения.

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

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

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

Особый интерес к физике взаимодействия УКЛИ с веществом возникает из-за существования так называемого двухтемпературного или 2Т-состояния, в котором облученное вещество пребывает в течение нескольких пикосекунд после облучения. Основы теории 2Т-СОСТОЯНИЯ были заложены в середине прошлого века в работах [1, 2]. Спустя несколько десятков лет работы [3, 4], опирающиеся на теорию 2Т-состояния, позволили объяснить результаты экспериментов [5], по облучению металлов и диэлектриков с помощью УКЛИ. Дальнейшее развитие теория получила в работах [6, 8, 7, 9, 10, 11]. Также с ее помощью были объяснены интересные эксперименты по «окрашиванию» металлов фемтосекундными лазерными импульсами А.Я. Воробьева и Ч. Гуо [12]: описание механизма образования рельефа, приводящего к изменению цвета металлов приведено в работе [13].

Оптическая диагностика в упомянутых выше экспериментах по облучению тонких пленок производилась с помощью методики pump-probe, детальное описание которой можно найти в работе В.В. Темнова [14]. От исходного лазерного импульса с помощью оптической системы отделяется менее мощный зондирующий импульс, попадающий на облучаемый материал с некоторой временной задержкой. Использование этой методики позволяет с высокой точностью измерять оптические и кинематические характеристики либо самой облучаемой поверхности, либо, наоборот, задней поверхности пленки. Вместе с масс-спектрометрией методика pump-probe составляет наиболее распространенный арсенал экспериментальной диагностики в опытах с УКЛИ. Однако, несмотря на высокую точность измерений, эти методы дают лишь косвенную информацию о внутренней «кухне» процессов. Аналитические же подходы, наоборот, не обеспечивают достаточной точности.

Поэтому необходимыми являются методы исследования, основанные на математическом моделировании процессов, происходящих при облучении металлов и диэлектриков УКЛИ. В настоящее время распространены и признаны два подхода: молекулярно-динамическое моделирование и моделирование на основе двухтемпературной гидродинамики.

Описания наиболее известных в последнее десятилетие моделей и кодов по молекулярно-динамическому численному моделированию взаимодействия УКЛИ с металлами и диэлектриками опубликованы в работах Д.С. Иванова и Л.В. Жигклея [15], В.В. Жаховского [4, 16, 17], С. Амору-зо [18], Дж.К. Чена [19], Г.Э. Нормана и С.В. Старикова [20]. Эти методы требовательны к вычислительным ресурсам и развиваются параллельно с развитием супер-ЭВМ и параллельных вычислений.

2Т-гидродинамические модели и коды, в отличие от молекулярно-динамических, могут работать и на обычном настольном компьютере, поскольку для получения результатов достаточно решать одномерные уравнения гидродинамики. Эти методы эволюционировали из однотем-пературных моделей, использовавшихся еще в 70-е годы прошлого века. Как правило, разные подходы уделяют большое внимание тому или иному аспекту процесса поглощения УК ЛИ. Пионерской в этом направлении можно назвать модель М.Ф. Иванова [21], активно применявшуюся до конца прошлого века. Последние результаты, полученные с помощью этой модели, используются, например, в работе [3]. В 1984 году появился хорошо известный сейчас код МШИ [22], изначально ориен--тированный на учет переноса излучения в плазме и режимы облучения наносекундными импульсами с максимальными температурами порядка десятков эВ и более с образованием горячей плазменной короны. В модели М.Е. Поварницына [23] для более точного учета течения в двухфазной области использован феноменологический подход, учитывающий кинетику образования кавитационных пузырьков при растягивающих напряжениях в облучаемой пленке. В модели [10], используемой французским физиком Ж.П. Коломбье, с помощью решения уравнения Гельмголь-ца более точно учитывается процесс поглощения ритр-луча в пленке. Также существуют подходы, опирающиеся на явное выделение фронтов фазовых переходов, такие, как, например, подход В.И. Мажукина [24], [25], [26], использующий решение задачи Стефана. Особенностью подхода, предлагаемого в настоящей работе, является использование современной широкодиапазонной двухтемпературной модели электронной теплопроводности [27, 28] и многофазного полуэмпирического табличного уравнения состояния [29, 30, 23].

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

Диссертация состоит из трех глав, введения и заключения. Первая глава диссертации посвящена детальному описанию модели физических процессов. Раздел 1.1. содержит физическую постановку задачи и уравнения двухтемпературной гидродинамики, которые являются ядром модели. В разделе 1.2 кратко описывается подход к описанию электронной теплопроводноЬти и его сравнение с другими известными моделям^. Приведен алгоритм по быстрому расчету коэффициента теплопроводности. Раздел 1.3 содержит описание фазовой диаграммы алюминия и описание уравнений состояния: табличного уравнения для ионной подсистемы и аналитической аппроксимации между классичесским и фермиевским газом для электронной подсистемы. В разделе 1.4 приведены возможные конфигурации экспериментальных установок и зависящие от них формальные начальные и граничные условия для уравнений двухтемпературной гидродинамики из раздела 1.1.

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

Третья глава диссертации полностью посвящена изложению результатов, полученных с помощью настоящей модели, их анализу и сравнению с экспериментом. В разделах 3.1-3.2 приводятся результаты численного моделирования эксперимента по облучению алюминиевой пленки единичным фемтосекундным лазерным импульсом^ проходящим первоначально через толстую стеклянную подложку. Описывается обнаруженный при помощи численного моделирования эффект «сверхзвукового плавления» металла. Производится анализ формы профиля нелинейной волны сжатия, возникающей в пленке после короткого периода релаксации двухтемпературного слоя и анализ влияния процесса плавления на движение задней границы пленки. В разделе 3.3 приведены результаты исследования фемтосекундной лазерной абляции алюминия с помощью настоящей модели и сравнение полученных результатов с данными известных экспериментов [5, 31] 90-х годов прошлого века.

Результаты работы опубликованы в 7 статьях в рецензируемых на-учйых журналах (см. Список литературы, [32] - [38]) и были включены в 6 докладов на международных конференциях, посвященных физике лазерных и плазменных процессов (см. Список литературы, [39] - [44]).

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

Автор выражает благодарность А.Г. Аксенову, H.A. Иногамову, A.B. Конюхову, П. И. Лебедеву, П. Р. Левашову, В. И. Мажукину,

A.M. Опарину, Ю.В. Петрову, M.B. Поварницыну, O.B. Трошкину, K.B. Хищенко, В.А. Хохлову за ценные замечания и комментарии, полезные обсуждения, дискуссии, предоставление полезных материалов и помощь в подготовке настоящей работы, в проектировании и написании кода и в получении результатов.

Заключение

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

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

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

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

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

1. M.И. Каганов, И.M. Лифшиц, Л.В. Танатаров, Релаксация между электронами и решеткой, ЖЭТФ, том 31, вып. 2(8) (1956), с. 232237.

2. С.И. Анисимов, Б.Л. Капелиович, Т.Л. Перельман, Электронная эмиссия с поверхности металлов под действием ультракоротких лазерных импульсов, ЖЭТФ, том 66, вып. 2 (1974), с. 776-781.

3. Н.А. Иногамов, A.M. Опарин, Ю.В. Петров, Н.В. Шапошников, С.И. Анисимов, Д. фон дер Линде, Ю. Майер-тер-Фен, Разлет вещества, нагретого ультракоротким лазерным импульсом, Письма в ЖЭТФ, 1999, том 69, вып. 4, с. 284-289.

4. С.И. Анисимов, В.В. Жаховский, Н.А. Иногамов, К. Нишихара, A.M. Опарин, Ю.В. Петров, Разрушение твердой пленки в результате действия ультракороткого лазерного импульса, Письма в ЖЭТФ, 2003, 77(11), с. 731-736.

5. Sokolowski-Tinten К., Bialkowski J., Cavalleri A., Von der Linde D., Oparin A., Meyer-ter-Vehn J., Anisimov S.I., Transient States of Matter during Short Pulse Laser Ablation, Phys. Rev. Lett.,1998, V. 81., p. 224227.

6. Yu.V. Petrov, Energy exchange between the lattice and electrons in a metal under femtosecond laser irradiation, Laser and Particle Beams (2005), 23, p. 283-289.

7. Zh. Lin, L.V. Zhigilei, and V. Celli, Phys. Rev. В 77, 075133 (2008).

8. N.A. Inogamov, V.V. Zhakhovskii, S.I. Ashitkov, V.A. Khokhlov, Yu.V. Petrov, P.S. Komarov, M.B. Agranat, S.I. Anisimov and K. Nishihara, Two-temperature relaxation and melting after absorption of femtosecond laser pulse, Appl. Surf. Sci., 255, 24, 2009.

9. J.P. Colombier, P. Combis, E. Audoard, and R. Stoian, Transient optical response of ultrafast nonequilibrium excited: Effects of electron-electron contribution to collisional absorption, Phys. Rev. E 77, 036409 (2008).

10. D. Fisher, M. Fraenkel, Z. Henis, E. Moshe, and S. Eliezer, Interband and intraband (Drude) contributions to femtosecond laser absorption in aluminum, Phys. Rev. E 65, 016409 (2001).

11. A.Y. Vorobyev and Chunlei Guo, Colorizing metals with femtosecond laser pulses, Appl. Phys. Lett. 92, 041914 (2008).

12. B.B. Жаховский, H.А. Иногамов, К. Нишихара, Новый механизм формирования нанорельефа поверхности, облученной фемтосекунд-ным лазерным импульсом, Письма в ЖЭТФ, т. 87, с. 491-496 (2008).

13. Vasiliy V. Temnov, Klaus Sokolowskii-Tinten, Ping Zhou, and Dietrich von der Linde, Ultrafast imaging interferometry at femtosecond-laser-excited surfaces, J. Opt. Soc. Am. B, Vol. 23, No. 9, (2006).

14. Dmitry S. Ivanov and Leonid V. Zhigilei, Combined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films, Phys. Rev. B, 68, 064114 (2003).

15. Анисимов С.И., Жаховский В.В., Иногамов H.А., Нишихара К., Петров Ю.В., Хохлов В.А., Разлет вещества и формирование кратера под действием ультракороткого лазерного импульса, ЖЭТФ, 2006, Том 130, Вып. 2(8), стр. 212-227.

16. V.V. Zhakhovski, N.A. Inogamov, Yu.V. Petrov, S.I. Ashitkov, К. Nishihara, Molecular dynamics simulation of femtosecond ablation and spallation with different interatomic potentials, Appl. Surf. Sci. (2009), DOI: 10.1016/j.apsusc.2009.04.082.

17. S. Amoruso, R. Bruzzese, X. Wang, N.N. Nedialkov and P.A. Atanasov, Femtosecond laser ablation of nickel in vacuum, J. Phys. D: Appl. Phys. 40 (2007) 331-340.

18. Yong Gan and J.K. Chen, Integrated continuum-atomistic modeling of nonthermal ablation of gold nanofilms by femtosecond lasers, Appl. Phys. Lett. 94, 201116 (2009).

19. С.И. Анисимов, М.Ф. Иванов, Н.А. Иногамов, П.П. Пашинин, М.А. Прохоров, Численное моделирование лазерного нагревания и сжатия простых оболочечных мишеней, Физика плазмы, 3(4), 723732 (1977).

20. К. Eidmann, J. Meyer-ter-Vehn, Th. Schlegel, S. Huiler, Hydro-code simulation of subpicosecond laser interaction with solid matter, Phys. Rev. E 62, 1202 (2000).

21. M.E. Povarnitsyn, Т.Е. Itina, M. Sentis, K.V. Khishchenko, and P.R. Levashov, Material decomposition mechanisms in femtosecond laser interactions with metals, Phys. Rev. В 75, 235414 (2007).

22. В.И. Мажукин, А.В. Мажукин, М.Г. Лобок, Динамика фазовых переходов и перегретых метастабильных состояний при нано-фемтосекундном лазерном воздействии на металлические мишени, Математичекое моделирование, 2009, т.21, н.11, стр. 99-112.

23. В.И. Мажукин, А.В. Мажукин, М.М. Демин, А.В. Шапранов, Эффекты неравновесности при воздействии импульсного лазерного излучения на металлы, Оптический журнал 78, 8, 2011.

24. A.V. Mazhukin, V.I.Mazhukin, М.М. Demin, Modeling of femtosecond ablation of aluminium film with single laser pulses, Appl. Surf. Sci., 257(2011) 5443-5446.

25. Н.А. Иногамов, Ю.В. Петров, Теплопроводность металлов с горячими электронами, ЖЭТФ, т. 137, вып. 3 (2010), с. 505-529.

26. Бушман A.B., Ломоносов И.В., Фортов В.Е., Уравнения состояния металлов при высоких плотностях энергии. Черноголовка: ИХФЧ РАН, 1992.

27. М.Е. Поварницын, Т.Е. Итина, П.Р. Левашов, К.В. Хищенко "Моделирование абляции металлических мишеней фемтосекундными лазерными импульсами"// Сборник трудов "Физика экстремальных состояний вещества 2007Черноголовка: ИПХФ РАН, С. 16 (2007).

28. H.A. Иногамов, В.В. Жаховский, В.А. Хохлов, В.В. Шепелев, Сверхупругость и распространение ударных волн в кристаллах, Письма в ЖЭТФ, т. 93, вып. 4, с. 245-251, (2011).

29. Y.T. Lee and R.M. More, An electron conductivity model for dense plasmas, Phys. Fluids 27, 1273 (1984).

30. S.I. Anisimov, B. Rethfeld, Theory of ultrashort laser pulse interaction with a metal, Proc. SPIE 3093, 192 (1997).

31. Ландау JI.Д., Лифшиц Е.М., Теоретическая физика: учеб. пособие в 10 т., Т. V, Статистическая физика, Москва, Наука, 1988 г.

32. Самарский А.А., Попов Ю.П., Разностные методы решения задач газовой динамики, Наука, Москва, 1980.