автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Прецизионные модели ударных адиабат и база ТЕФИС

кандидата физико-математических наук
Иванченко, Евгений Сергеевич
город
Москва
год
2009
специальность ВАК РФ
05.13.18
Диссертация по информатике, вычислительной технике и управлению на тему «Прецизионные модели ударных адиабат и база ТЕФИС»

Автореферат диссертации по теме "Прецизионные модели ударных адиабат и база ТЕФИС"

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

Иванченко Евгений Сергеевич

Прецизионные модели ударных адиабат и база ТЕФИС.

Специальность - 05.13.18

математическое моделирование, численные методы и комплексы

программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

□0348Э480

Москва-2009

003489488

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

Научный руководитель

•Доктор физико-математических наук, член-корреспондент РАН Калиткин Николай Николаевич

Официальные оппоненты

■ Доктор физико-математических наук, Якобовский Михаил Владимирович.

-Доктор технических наук, Фунтиков Александр Иосифович.

Ведущая организация ФГУП РФЯЦ ВНИИЭФ.

Защита состоится «_»_2010 года в _часов_минут

на заседании диссертационного совета Д 002.058.01 при институте математического моделирования РАН по адресу: 125047, Москва, Миусская пл., 4а.

С диссертацией можно ознакомиться в библиотеке института математического моделирования РАН.

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

Ученый секретарь диссертационного совета доктор физико-математических наук

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ.

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

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

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

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

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

Первые эксперименты были статическими, в них исследовалась изотермическая сжимаемость веществ, и первый в этой области был Бриджмен [1]. Он исследовал сжимаемость нескольких десятков веществ до давлений в 50-100 кБар. На сегодняшний день наиболее эффективным направлением экспериментального подхода являются эксперименты по ударным сжатиям веществ: взрывчатые вещества создают давления до 0,5 Мбар. Разгоняя металлическую пластину

3 4

\

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

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

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

Для решения этой проблемы может быть предложено два основных направления:

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

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

Однако гораздо более выгодно использовать объединение этих подходов.

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

Собственно отсюда вытекает постановка задачи:

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

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

3. Естественно необходимо создать хранилище (базу данных), где могли бы храниться результаты экспериментов.

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

Исторически впервые задачи по исследованию поведений веществ в экстремальных условиях были поставлены в США [2], там начиная с первой половины 1940ых годо, в лабораториях Лос Аламоса [3] начались проводиться эксперименты по ударным сжатиям.

Си

и км/с

Рис 1 Участок ударной адиабаты меди.

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

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

Пример подобной ситуации приведен на рис. 1. В этом случае фактически невозможно, что либо сказать о поведении кривой.

Для преодоления этой проблемы в лабораториях Ливермора [4] в 1949 под руководством Эдварда Теллера была заложена база данных по теплофизическим свойствам веществ, получившая название СЕЗАМ. В эту базу вошли многолетние результаты работ Лос Аламоса [5], Ливермора [6], и лабораторий военно-морского флота САНТИЯ.

На сегодняшний день это наверняка крупнейшая база данных, но есть одна загвоздка: большая часть данных сосредоточена в закрытой части этой базы, к которой естественно мы доступа не имеем. Помимо закрытой части СЕЗАМ содержит так же и открытую часть, в которой содержится некоторое количество экспериментальных данных, но в основном содержимое открытой части базы сезам - это огромное количество модельных предположений. Причем в большинстве случаев ни о точности расчетов, ни об используемых моделях, ни о диапазоне адекватности этих моделей информации нет. Кроме этого в базе содержится достаточно большое количество неточностей, и подчас даже противоречивых данных (Возможно внесенных туда специально). Все это несколько затрудняет использование открытой части базы СЕЗАМ в практических целях.

Параллельно с исследованиями в США, в СССР начиная 1946 года [7] в ядерном центре в Сарове [8] при содействии таких выдающихся ученых как Н.Г. Харитонова, Я.Б. Зельдовича, Л.В. Альтшулера, H.A. Дмитриева, К.А.Семендяева, А.Д. Сахарова, Д.А. Ф ранка-Каменецкого, Е.И. Забабахина, в своем роде лидеров новой научной дисциплины -физики высоких плотностей, были также основаны лаборатории по исследованиям поведений веществ в экстремальных состояниях. За более чем полувековую историю у нас было проведено огромное количество экспериментов. Результаты этих экспериментов и стали явились фундаментом создаваемой в Сарове, под руководством Л.В. Альтшулера, собственного объединенного компендиума.

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

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

Создание подобной базы - это очень трудоемкий, но весьма важный проект. Одной из таких реализаций является база данных Института Теплофизики Экстремальных Состояний [9]. В эту базу уже включена большая часть данных из большинства упомянутых ранее источников. База ИТЭС насчитывает порядка 17000 точек. Однако при проведении сверок с последними редакциями оригинальных компендиумов выяснилось наличие пробелов в данных, ошибок и просто опечаток. Подобного рода недочеты могут влиять на результаты последующих обработок, в особенности, если речь идет об обработке пористых веществ.

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

Цель работы.

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

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

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

У данной работы было две основных цели:

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

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

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

1) Модель Томаса-Ферми с квантовой и обменной поправками.

2) Модель Ионизационного Равновесия.

3) Модель Электронного Переноса.

В дополнение ко всему вышеперечисленному, база ТЕФИС содержит в себе таблицы термодинамических величин, рассчитанных по разным моделям на достаточно подробных сетках. Что при наличии хороших методов аппроксимации, (в том числе и включенных в комплекс) позволит избежать временных затрат на расчет термодинамических таблиц для решения целого ряда прикладных задач при использовании базы ТЕФИС.

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

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

Практическая ценность работы.

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

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

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

выполнена как единый, интегрированный комплекс. Кроме того она рассчитана на дальнейшее расширение.

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

Результаты работы докладывались на международных конференциях:

3rd Moscow Workshop on Targets and Applications, October, 2007 Moscow, Russia.

International Workshop «Numerical Analysis and Scientific Computing» (NASCom'08) 13-17 октября 2008, Ростов-на-Дону.

Physics of high-energy density in matter. November 24-25 2008. Москва.

XI Харитоновские научные чтения. Экстремальные состояния вещества. Детонация. Ударные волны. 2009 г. Саров.

Работа выполнена при поддержке грантов:

• РФФИ:08-01-00171

• РФФИ:05-01-00152

• РФФИ 05-01-08006

• НШ-1918.2003.1

• НШ-5772.2006.1

Структура и объем работы.

Диссертация состоит из введения, 4х глав, заключения и приложений. Итого - пять глав и приложение.

Общий объем -100 страниц.

Публикации.

По теме работы сделано 5 статей:

Е.С.Иванченко, Н.Н.Калиткин, Л.В.Кузьмина. Главные ударные адиабаты в базе теплофизических свойств веществ (ТЕФИС). Математическое моделирование. Том 20, №б, 2008г. Стр. 99-110.

Е.С.Иванченко, H.H.Калиткин, Л.В.Кузьмина. Квантово-статистический расчет термодинамики и ударных адиабат (программы КСМ). Математическое моделирование. Том 20, N27, 2008г. Стр. 30-44.

TEFIS - Database on Thermophysical Properties of Substances under Extreme conditions E.S. Ivanchenko, N.N. Kalitkin, I.A. Kozlitin, L.V.Kuzmina, K.I.Lutskiy, I.A. Panin

E.C. Иванченко, H.H. Калиткин. Глобальные ударные адиабаты пористых веществ и база ТЕФИС. с. 142-143. XI Харитоновские чтения. Сборник тезисов докладов. 2009 г. Саров.

E.S. Ivanchenko, N.N. Kalitkin Wide-range Hugoniots of Porous Substances and the TEFIS database, p. 143-144. XI Харитоновские чтения. Сборник тезисов докладов. 2009 г. Саров.

Иванченко Е.С. Калиткин H.H. Интерпретация ударно волновых данных в базе ТЕФИС. Математическое Моделирование. Том 22, №1, 2010 г. Москва.

Краткое содержание работы. Первая глава.

Содержит введение, обоснование актуальности темы, формулировку основных целей, начальные сведения, краткий обзор комплекса ТЕФИС и краткое содержание остальных глав.

Вторая глава.

Описывает используемые в базе ТЕФИС методы расчета термодинамических функций, необходимые для расчета асимптотик ударных адиабат. Глава содержит обоснование выбора именно этих методов и демонстрирует их эффективность. Схематически диапазоны применения моделей изображены на (Рис. 2).

В этой главе описывается модель Томаса-Ферми с квантовыми и обменными поправками, как модель, принятая для расчета сверхплотной плазмы в базе ТЕФИС.

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

Рис 2 Зоны адекватности моделей расчета термодинамики в базе ТЕФИС.

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

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

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

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

Третья глава.

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

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

Е-Е0=0.5(Р + Р0)(1/р0-1/р) p = p0D/(D-U),P = P0+pDU

Описываются требования к точности решения, и приводятся решения, опирающиеся как на классические методы [10], [11] так и на оригинальные методы с повышением точности. В совокупности предложенные методики позволяют получить точность порядка 0.1%.

В следующем разделе этой главы выводятся асимптотики ударных адиабат в диапазоне адекватности модели Томаса-Ферми с поправками. Эти асимптотики имеют следующий вид:

dkcu (U) ~ С + BU + ¿U\C > 0, В > 1,1» А > О Приводятся оценки точности данной аппроксимации.

Выводятся зависимости коэффициентов А,В,С от пористости материала, имеющие следующий вид:

lg С{т) и C + y\gm,

lg А(т) ^ a + a\gm,

lg[В(т) + 2yjA(m)C(m) -1] * b + /? lg m

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

DE№(U)*c0 + b0U

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

Copper (Си)

U km/s

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

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

Классически это решение основывается на построении следующей кривой:

В(и) = Вксм(и)- С~С\ =

1 + /ии + ци С-с0

Где ¿0,с0 это параметры линейной аппроксимации, рассчитанные на

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

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

В работе предложен несколько другой подход. Предложено искать широкодиапазонную ударную адиабату в следующем виде:

Причем, а и /3 вычисляются с учетом всех экспериментальных данных, а Окш {и) - это параболическая аппроксимация на основе квантово статистической модели.

Приводятся оценки эффективностей различных функций Р(х) и оценки эффективности наиболее успешной из них:

/?(*) =-1-

1 + х(1 + х)

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

И в заключение приводится метод расчета широкодиапазонных ударных адиабат для пористых веществ, показывается эффективность этого метода для расчетов пористых веществ с пористостью т$10 на примере веществ с наиболее представительной выборкой экспериментов с пористыми образцами (см. рис. 5).

Copper (Си)

Рис 4 Широкодиапазонная ударная адиабата меди. Звездочки - результаты расчета на основе модели ТФП. Линия - широкодиапазонная ударная адиобата. Другие маркеры - результаты экспериментов различных исследовательских групп.

Lead (РЬ)

ol£_1_1_1_1_L

0 1 2 3 4 5

UknVi

Рис 5 Ударные адиабаты для свинца с различными пористостями.

Четвертая глава.

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

Пятая глава.

Пятая глава полностью посвящена описанию архитектуры комплекса ТЕФИС, его схематичная структура изображена на Рис б.

Рис 6 Архитектура комплекса ТЕФИС.

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

Для того чтобы удовлетворить всем этим требованиям в работе был выбран достаточно нестандартный подход. В качестве формата данных для экспорта базы данных из приложения был выбран стандарт XML, с жестким определением схемы на языке XSD. В качестве основного хранилища был выбран ' Microsoft SQL Server Compact Edition -бесплатно распространяемая версия известного SQL сервера, отличающаяся тем, что нет необходимости устанавливать дополнительные компоненты. Это позволило сделать базу предельно простой и легко расширяемой. Кроме того, Microsoft SQL Server Compact Edition отличается весьма высокой производительностью, что позволяет проводить различные выборки из базы данных с огромной скоростью.

Так же описывается архитектура управляющей системы комплекса ТЕФИС, описывается дизайн пользовательского интерфейса. И наконец, описывается методика интеграций различных модулей, выполненных в разных языковых средах (С++, Fortran, MATLAB), в том числе и модулей, выполняющих расчеты на основе разработанных моделей.

Для реализации основной управляющей среды (менеджера) используется популярная среда Microsoft .NET Framework 3.5. Использование этой системы позволило реализовать удобный доступ, как к базе данных, так и к реализациям моделей, которые реализовывались на С++, Fortran и MATLAB. Столь тесное переплетение различных языковых платформ позволило более направлено использовать их инструментарии. Так, пользовательский интерфейс и доступ к данным реализовывался на языке С#, который позволяет быстро и эффективно решать подобные задачи, однако он не достаточно производителен при серьезных математических расчетах.

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

10

Name

ChemicalAbbriviatute Net« Phase

И Ю

s Г .. v ■ й

[fe_Kll,GetData 0

' J .......

И ю

| MateriallD

^ MethodID

.. >>, ...... ... ISlsI

■У* FiU,GetD>ta о ]

/ ^ | • ID

MateriallD

I x

i ElementlD

S FilluetOalaO

■-1 ' .

id

Name

Notes

Ж

H!j FiltGetData 0

: -J

ij* ID

lgT

TablelD

в

НМпШМк

hat Fill.GetData 0

I Level | Valut I ElementlD

p FiHGetDataO^

Filt.GctDHi 0

1

НШМ ШШШ

Vf 2

Name

ChemicalAbbriviature

ШШЛШш/х-'- ййейя ;

i^ii Fill.GetData 0

■1

pt; Fill,GetData ij

Рис 7 Схема базы данных.

В заключительной части этой главы приводится сравнение комплекса ТЕФИС с существующими базами данных Сезам и ИТЭС. Выявляются и анализируются основные отличия, достоинства и недостатки комплекса.

Заключение.

В заключении сформулированы основные результаты работы. Приложение.

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

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

И наконец, приведены исходные коды программ и модулей комплекса ТЕФИС.

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

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

^ Разработка прецизионных алгоритмов расчета ударных адиабат.

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

^ Создание методов итерационного повышения точности широкодиапазонных ударных адиабат.

^ Создание пополняемой базы данных ТЕФИС, объединившей в себе все доступные на сегодняшний день компендиумы ударно-волновых экспериментальных исследований.

^ Создание комплекса ТЕФИС объединившего в себе как базу данных ТЕФИС, так и целый ряд моделей и программ для расчета теплофизических свойств веществ в экстремальных состояниях.

Литература.

1. Бриджмен П.В. // Твердые тела под высоким давлением. / Под ред. В. Пола и Д. Варшауэра. М.: Мир, 1966.

2. R.W. Goranson, D. Bancroft, B.L. Burton et al. // Ibid. 1955. Vol. 26, N 12. P. 1472-1479.

3. http://en.wikipedia.org/wiki/Los Alamos National Laboratory

4. http://en.wikipedia.org/wiki/Lawrence Livermore National La boratory

5. LASL shock Hugoniot data /Ed. S.P. Marsh. Berkley etc.: Univ. Calif, press. 1980.

6. Compendium of shock wave data / Ed. M. van Thiel. Livermore: Univ. Calif, press. 1977.

7. Л.В. Альтшулер, K.K. Крупников // История советского атомного проекта: (40-е и 50-е годы): Международный научный симпозиум ИСАП-96. М.: ИздАТ., 1997. Т. 1 Стр. 184-191.

8. http://ru.wikipedia.org/wiki/CapoB

9. http://www.ihed.ras.ru/rusbank/

10. H.H. Калиткин, Л.В. Кузьмина, Широкодиапазонные ударные адиабаты, стр. 107-120. //Ударные волны и Экстремальные состояния вещества, под редакцией В.Е. Фортова, Л.В. Альтшулера, Р.Ф. Трунина, А.И. Фунтикова. Москва "Наука" 2000г.

11. H.H. Калиткин, Л.В. Кузьмина, Метод расчета ударных адиабат: Препр. ИПМат. АН СССР №159. М., 1988, 24 С.

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

Формат 60x84 1/16. Уч.-изд.л. ^Тираж/00 экз. Заказ 7М

Отпечатано в типографии КПК МИЭТ.

124498, Москва, г.Зеленоград, проезд4806, д.5, МИЭТ.

Оглавление автор диссертации — кандидата физико-математических наук Иванченко, Евгений Сергеевич

Оглавление.

Глава 1. Введение.!.

Параграф 1. Возникновение задачи.

Параграф 2. Методики проведения экспериментов.

Смещение поршня.

Алмазные наковальни (АН).

Изоэнтропическое сжатие.

Атомные пучки (АП).

Ударные сжатия.

Метод торможения.

Метод у-репера [16].

Метод отражения.

Сравнительная сжимаемость.1В

Двукратное ударное сжатие.

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

Параграф 3. Проблема анализа точности.

Параграф 4. Знакомство с базой ТЕФИС.

Глава 2. Термодинамика в базе ТЕФИС.

Параграф 1. Используемые модели.

Параграф 2. Модель Томаса-Ферми с квантовой и обменной поправками.

Модель Томаса-Ферми.

Квантово-статистическая модель (КСМ).

Введение 2009 год, диссертация по информатике, вычислительной технике и управлению, Иванченко, Евгений Сергеевич

Смеси элементов.30

Параграф 3. Модель ионизационного равновесия.31

Флуктуирующее микрополе.31

Модель Хольцмарка и ее уточнения.'.32

Обобщенные уравнения Саха с учетом вырождения.34

Термодинамические функции.36

Параграф 4. Квазизонная интерполяция.38

Глава 3. Ударные адиабаты.42

Параграф 1. Ударные адиабаты.42

Параграф 2. Квантово-статистические ударные адиабаты.52

Параграф 3. Расчет квантово-статистических ударных адиабат в пористых веществах. 53

Параграф 4. Оценка экспериментальных данных.54

Параграф 5. Широкодиапазонные ударные адиабаты.55

Параграф 6. Широкодиапазонные ударные адиабаты в пористых веществах.58

Глава 4. Оценка экспериментальной точности.61

Параграф 1. Предпосылки.61

Параграф 2. Статистическая обработка.62

Параграф 3. Определение «плохих» точек.66

Параграф 4. Метод повышения точности.67

Глава 5. Архитектура комплекса.71

Параграф 1. Общая архитектура комплекса.71

Параграф 2. База данных.72

Выбор базы.73

Разработка базы данных.74

Схема базы данных.76

Параграф 3. Пользовательский интерфейс.80

Оконный интерфейс.80

Программный интерфейс.83

Параграф 4. Сценарий использования комплекса ТЕФИС.84

Параграф 5. Реализации моделей.85

Параграф 6. Сравнение с существующими проектами.86

Open Sesame.86

База данных Института Теплофизики Экстремальных Состояний.87

Заключение.88

Заключение.90

Опубликованная литература.91

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

Глава 6. Приложение.92

Параграф 1. Таблицы.92

Параграф 2. Графики широкодиапазонных ударных адиабат.100

Параграф 3. Графики начальных участков ударных адиабат.109

Список литературы.141

Заключение диссертация на тему "Прецизионные модели ударных адиабат и база ТЕФИС"

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

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

• Разработаны методы и реализованы алгоритмы прецизионного расчета широкодиапазонных ударных адиабат на основе заранее рассчитанных теоретических ударных адиабат и совокупности экспериментальных данных.

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

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

• Создана пополняемая база данных ТЕФИС, объединившая в себя все доступные на сегодняшний день компендиумы. С открытым форматом и с жестко определенной схемой.

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

Опубликованная литература.

Al] Е. С. Иванченко, Н. Н. Калиткин, Л.В. Кузьмина "Главные ударные адиабаты в базе теллофизических свойств веществ (ТЕФИС)", Москва, Математическое моделирование том 20. Номер б, стр. 99-110. 2008. [А2] Е. С. Иванченко, Н. Н. Калиткин, Л.В. Кузьмина "Квантово-статистический расчет термодинамики и ударных адиабат", Москва, Математическое моделирование том 20. Номер 7, стр. 30-44. 2008. [A3] TEFIS - Database on Thermophysical Properties of Substances under Extreme conditions. E.S. Ivanchenko, N.N. Kalitkin, I.A. Kozlitin, L.V.Kuzmina, K.I.Lutskiy, l.A. Panin

A4] E.C. Иванченко, H.H. Калиткин. Глобальные ударные адиабаты пористых веществ и база ТЕФИС. с. 142-143. XI Харитоновские чтения. Сборник тезисов докладов. 2009 г. Саров. [А5] E.S. Ivanchenko, N.N. Kalitkin Wide-range Hugoniots of Porous Substances and the TEFIS database, p. 143-144. XI Харитоновские чтения. Сборник тезисов докладов. 2009 г. Саров. [А6] Иванченко Е.С. Калиткин Н.Н. Интерпретация ударно волновых данных в базе ТЕФИС. Математическое Моделирование. Том 22, №1, 2010 г. Москва.

Апробация работы:

• 3rd Moscow Workshop on Targets and Applications, October, 2007 Moscow, Russia.

• International Workshop «Numerical Analysis and Scientific Computing» (NASCom'08) 13-17 октября 2008, Ростов-на-Дону.

• Physics of high-energy density in matter. November 24-25 2008. Москва.

• XI Харитоновские научные чтения. Экстремальные состояния вещества. Детонация. Ударные волны. 2009 г. Саров.

Заключение.

Библиография Иванченко, Евгений Сергеевич, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Бриджмен, П.В. Твердые тела под высоким давлением., Д. Варшауэра, В. Пол, Eds. / Москва, СССР: Мир, 1966.

2. Н. Н. Калиткин, Л. В. Кузьмина, "Проверка моделей вещества в экспериментах по сжатию.," АН СССР, Москва, препринт №116 1985.

3. F. P. Bundy. Design and development of apparatuses to achieve the highest possible static pressure. //X AIRAPT International high pressures conference, 1986, p. 42-51.

4. Тезисы докладов. // Всесоюзное совещание "Координате чувствительные фотоприемники и оптико-электронные устройства на их основе", 1985.

5. Н. Н. Калиткин, В. Б. Леонас, И. Д. Родионов, "Сжимаемость водорода при высоких давлениях," АН СССР, Москва, препринт №127 1985.

6. I. F. Silvern. New phases of molecular and atomic hydrogen under extreme conditions. // Europ. Phys. News., 1982, т. 13, №. 8-9, с. 4-6.

7. Ф. В. Григорьев, С. Б. Кормер, О. Л. Михайлова, А. П. Толочко, В. Д. Урлин. Уравнение состояния молекулярного водорода; о фазовом переходе в металлическое состояние. // ЖЭТФ, 1978, т. 75, №. 5(11), с. 1683-1692.

8. В. Б. Леонас, И. Д. Родионов. Исследования высокоэнергетического рассеяния атомов и молекул. //УФН, 1985, т. 146, №. 1, с. 7-34.

9. М. Ross, F. Н. Ree, D. A. Young. The equations of state of molecular hydrogen at very high densith. //J. Chem. Phys., 1983, т. 79, №. 3, с. 1487-1494.

10. Жирков B.H. Динамическин исследования твердых тел при высоких давлениях. / Сборник Москва: Мир, 1965. - 236 с.

11. Я В. Алыпшулер. Приминение ударных волн в физике высоких давлений. // УФН, 1965, т. 85, №. 2, с. 197-258.12. 77. В. Альтшулер, Б. Н. Моисеев, Л. В. Попов, Г. В. Симонов, Р. Ф. Трунин.

12. Сравнительная сжимаемость железа и свинца при давлениях 31-34 Мбар. // ЖЭТФ, 1968, т. 34, №. 3, с. 785-789.

13. Р. Ф. Трунин, М. А. Подурец, Г. В. Симонов, Л. В. Попов, Б. Н. Моисеев. Экспериментальная проверка модели Томаса-Ферми для металлов при высоких давлениях. // ЖЭТФ, 1972, т. 62, №. 3, с. 1043-1048.

14. Е. Н. Аврорин, Б. К. Водолага, Л. П. Волков, А. С. Владимиров, В. А. Симоненио, Б. Т. Черноволюк. Ударная сжимаемость свинца, кварца, аллюминия и воды при давлении ~100Мбар. // ЖЭТФ, 1980, т. 31, №. 12, с. 727-729.

15. Е. Н. Аврорин, Б. К. Водолага, Н. П. Волошин, В. Ф. Куропатенко, Г. В. Коваленко, В. А. Симоненко, Б. Т. Черноволюк. Экспериментальное подтверждение оболочечных эффектов на ударных адиабатах аллюминия и свинца. // ЖЭТФ, 1986, т. 43, N9. 5, с. 241-244.

16. В. А. Симоненко, Н. П. Волошин, А. С. Владимиров, А. П. Нагибин, В. Н. Ногин, В. А. Попов, В. А. Сальников, Ю. А. Шойдин. Абсолютное изменение сжимаемости алюминия при давлениях ~1ТПа. // ЖЭТФ, 1985, т. 88, №. 4, с. 1452-1459.

17. М. van. Thiel, "Compendium of shock wave data.," Lawrence Livwermore National laboratory, CA, Report 50108,1977.

18. S.P. Marsh. LASL shock Hugoniot data. / Los Angeles London: University of California Press, Berkley, 1980.

19. J. P. Romain, F. Cottet, M. Hallouin, R. Fabbro, B. Faral, H. Pepin. Laser shock experiments at pressures above 100 Mbar. // X AIRAPT International high pressure conference, 1986, p. 595-598.

20. В. E. Фортов, И. Т. Якубов. Физика неидеальной плазмы. / Черноголовка, СССР: ОИХФ АН СССР, 1984. 263 с.

21. W. J. Nellis, А. С. Mitchell, М. van Thiel, G. J. Devine, R. J. Trainor. Equations-of-state data for molecular hydrogen and deuterium at shock pressure in range 2-76 Gpa (20-760 Kbar). //J. Chem. Phys., 1983, т. 79, №. 3, с. 1480-1486.

22. W. J. Nellis, N. C. Holmes, A. C. Mitchell, R. J. Trainor, G. K. Governo, M. Ross, D. A. Young. Shock compression of liquid helium to 56 Gpa (560 kbar). // Phys. Rev. Letters, 1984, т.53, №. 13, с. 1248-1251.

23. Et. al. //X AIRAPT International high pressure conference, 1986, v. 139-140.

24. Я. Б. Зельдович. Об исследовании уравнения состояния с помощью мезанических измерений. // ЖЭТФ, 1957, т. 32, №. 2, с. 1577-1578.

25. В. Е. Фортов. Уравнения состояния конденсированных сред. // ПМТФ, 1972, №. 6, с. 156-166.

26. Online]. http://en.wikipedia.orR/wiki/Los Alamos National Laboratory

27. Online]. http://en.wikipedia.org/wiki/Lawrence Livermore National Laboratory

28. Л. В. Альтшулер, К. К. Крупников. История советского атомного проекта: (40-е и 50-е годы). // Международный научный симпозиум ИСАП-96, 1997, р. 184-191.

29. Р. Ф. Трунин, Ed. Экспериментальные данные по ударно-волновому сжатию и адиабатическому расширению конденсированных веществ., 2nd изд. / Саров, Россия, 2006.

30. База данных Института Теплофизики Экстремальных Состояний. Online]. http://www.ihed.ras.ru/rusbank/

31. Н. Н. Калиткин, Л. В. Кузьмина. Ударные адиабаты 83 веществ. // Химическая физика., 1999, т. 18, №. 10, с. 74-78.

32. Н. Н. Калиткин, Л. В. Кузьмина, "Таблицы термодинамических функций вещества при высоких концентрациях энергий.," ИПМ, Москва, Препринт 1975.

33. S. P. Rothman, А. М. Evans, R. Т. Eagleton, L. В. Pearson. Material properties experiments using AWE high power laser. // Int. J. Impact Eng., 1999, т. 20, с. 803.

34. H. H. Калиткин. Модели вещества в экстремальном состоянии. // Математическое Моделирование; Физикохимические свойства веществ., 1989, с. 114-167.

35. А. Киржниц. О границе применимости квазиклассического уравнения сотстояния вещества. // ЖЭТФ, 1958, т. 35, №. 6(12), с. 1545-1557.

36. А. Ф. Никифоров, В. Г. Новиков, В. Б. Уваров. Квантово-статистические моделивысокотемпературной п'лазмы. Методы расчета росселандовых пробегов и уравнений состояния. / Москва: Физматлит, 2000.

37. Д. А. Киржниц. Квантовые поправки к уравнению Томаса-Ферми. // ЖЭТФ, 1957, т. 32, №. 1, с. 115-123.

38. Н. Н. Калиткин, Л. В. Кузьмина, "Квантово-статистическое уравнение состояния и ударные адиабаты," ИПМат, Москва, препринт №14 1976.

39. Д. А. Киржниц, Г. В. Шпатаковская. // ЖЭТФ, 1995, т. 102, с. 1238.

40. Н. Н. Калиткин, Л. В. Кузьмина. Квантово статистическое уравнение состояния. // Физика плазмы, 1976, т. 2, №. 5, с. 858-868.

41. Н. Н. Калиткин, И. А. Козлитин. Компенсированная микрополевая модель неидеальности плазмы. // ДАН, 2006, т. 411, с. 36-40.

42. Н. Н. Калиткин, Козлитин И. А. Модель квазинезависимых частиц для плазменного микрополя. //ДАН, 2008, т. 418, №. 5, с. 614-618.

43. Н. Н. Калиткин, И. А. Козлитин. Самосогласованная компенсированная микрополевая модель неидеальности плазмы. // Математическое моделирование, 2008, т. 20, №. 1, с. 61-76.

44. Н. Н. Калиткин, И. А. Козлитин. Сравнение детального состава плазмы в различных моделях. // Математическое моделирование, 2008, т. 20, №. 4, с. 69-77.

45. И. А. Козлитин. Моделирование распределения Хольцмарка методом Монте-Карло. // Математическое моделирование, 2008, принято в печать.

46. Н. Н. Калиткин, И. А. Козлитин. Микрополевое обрезание статистических сумм в плазме. // Труды VIII Харитоновских чтений, 2006, р. 252-258.

47. Н. Н. Калиткин, И. А. Козлитин. Микрополе и компенсированная модель неидеальности плазмы. //Труды IX Харитоновских чтений, 2008, р. 253-258.

48. И. А. Козлитин, "Микрополевая модель квазинезависимых частиц и неидеальная плазма.," ИММ, Москва, Диссертационная работа. 2008.

49. А. Ф. Никифоров, В. Г. Новиков, В. Б. Уваров. Модифицированная модель Хартри

50. Фока-Слетера для вещества с заданной температурой и плотностью. // ВАНТ, серия: методики и программы., 1979, т. 4, №. 6, с. 16-35.

51. В. Г. Новиков, "Модифицированная модель Хартри-Фока-Слетера и ее применение для расчетов уравнения состояния вещества.," ИПМат, Москва, Кандидатская диссертация. 1981.

52. Н. Н. Калиткин, Л. В. Кузьмина, И. В. Ритус, "Табличное и графическое представление теплофизических свойств веществ (проект стандарта)," ИПМат, Москва, Препринт 1987.

53. Н. Н. Калиткин, К. И. Луцкий. Прецизионная аппроксимация квантово-статистических кривых холодного сжатия. // Математическое, 2008.

54. Н. Н. Калиткин, К. И. Луцкий. Аппроксимация гладких поверхностей методом двойного периода. // Математическое моделирование, 2008.

55. Holtsmark. //J. Ann. Phys., 1919, т. 577 , №. 58.

56. Г. Грим. Уширение спектральных линий в плазме. / Москва, СССР: Мир, 1978. 491 с.

57. И. О. Голосной. Аналитические аппроксимации для уширения спектральных линий водородоподобных ионов. // Физика плазмы, 2001, т. 27, N9. 2, с. 526-535.

58. Ю. К. Куриленков, В. С. Филинов. К теории микрополя в неидеальной плазме. // Теплофизика высоких температур., 1976, т. 14, №. 4, с. 886-888.

59. С. A. Iglesias, J. Т. Lebowitz. Electric microfield distribution in multicomponent plasmas. // Phys.Rev.A, 1984, т. 30, №. 4, с. 2001-2004.

60. J. W. Dufty, D. B. Boercker, C. A. Iglesias. Electric field distributions in strongly coupled plasmas. 11 Phys.Rev.A., 1985, т. 31, №. 3, с. 1681-1686.

61. С. A. Iglesias, H. Е. DeWitt, J. Т. Lebowitz, D. McGowan, W. B. Hubbard. Low-frequency electric microfield distribution in plasmas. // Phys.Rev.A., 1985, т. 31, №. 3, с. 16981702.

62. A A. Broyles. Stark fields from ions in plasmas. I I Phys.Rev.A., 1955, т. 100, №. 11811190.

63. М. Barange, В. Mozer. Electric field distributions in an ionized gas. // Phys.Rev., 1959, т. 115, №. 3, с. 521-525.

64. H. H. Калиткин, И. В. Ритус, А. М. Миронов, "Ионизационное равновесие с учетом вырождения электронов.," ИПМат, Москва, Препринт 46, 1983.

65. М. М. Баско, "Уравнение состояния металлов в приближении среднего иона.," ИТЭФ, Москва, Препринт 57, 1982.

66. Н. Н. Калиткин. Квазизонное уравнение состояния. // Математическое Моделирование, 1989, т. 1, №. 2, с. 64-108.

67. Н. Н. Калиткин, Л. В. Кузьмина, И. В. Ритус, "Квазизонная интерполяция термодинамических функций.," ИПМат, Москва, Препринт 1988.

68. Н. Н. Калиткин, Л. В. Кузьмина, "Программы расчета теплофизических свойств веществ ТФП1-6," ИПМ, Москва, Препринт 153,1984.

69. Е. С. Иванченко, Н. Н. Иалиткин, Л. В. Кузьмина. Главные ударные адиабаты в базе теплофизических свойств веществ (ТЕФИС). // ическое моделирование., 2008, т. 20, №. 6, с. 99-110.

70. Н. Н. Калиткин, И. А. Панин. О сходимости приближений Чепмена-Энскога при расчете электронных транспортных коэффициентов. // Математическое моделирование, 2008, т. 20, №. 6, с. 86-98.

71. Н. Н. Калиткин, И. А. Панин. Вычисление электронного переноса в плазме в высоких приближениях Чепмена-Энскога. // Математическое моделирование, т. 20, №. 4, с. 103-116.

72. Н. Н. Калиткин, И. А. Панин. Электронный перенос в неидеальной плазме. // Сборник трудов международной конференции "IX Харитоновские тематические научные чтения", 2007.

73. Саров. Online]. http://ru.wikipedia.org/wiki/CapoB

74. Демо версия WEB ресурса ТЕФИС. Online], http://

75. R. W. Goranson, D. Bancroft, В. L Burton. // Ibid, 1955, т. 26, с. 1472-1479.

76. Е. С. Иванченко, Н. Н. Иалиткин, Л. В. Кузьмина. Квантово-статистический расчет термодинамики и ударных адиабат (программы КСМ). // Математическое моделирование, 2008, т. 20, №. 7, с. 30-44.