Отчет ИВМ СО РАН за 2011 годПрограммы фундаментальных исследований Сибирского отделения РАН
IV.29. Системы автоматизации, CALS-технологии, математические модели и методы исследования сложных управляющих систем и процессовПрограмма IV.29.1. Теоретические основы и методы информационных и вычислительных технологий проектирования и принятия решений Координаторы программы: академик РАН Ю. И. Шокин, член-корреспондент РАН В. В. Шайдуров Проект IV.29.1.2 «Информационно-вычислительные технологии для анализа и проектирования сложных систем в ближнем космосе»№ гос. регистрации 01201056404 Научный руководитель проекта: член-корреспондент РАН В. В. Шайдуров Блок 1. Моделирование отклика магнитосферного хвоста на резкое падение динамического давления солнечного ветра, связанное с приходом к границе магнитосферы Земли фронта магнитного облака, движущегося от Солнца. Аналитико-численное исследование диффузионной области магнитного пересоединения в рамках двухжидкостной магнитогидродинамической модели (Н. В. Еркаев, В. В. Денисенко, А. В. Китаев). Разработана аналитическая холловская магнитогидродинамическая модель стационарного магнитного пересоединения в сжимаемой плазме. Показано, что задача сводится к решению уравнения Града-Шафранова относительно магнитного потенциала. Получено решение, на основе которого исследовано формирование двойных электрических слоев, а также слоев пониженной плотности плазмы. Показано, что двойные электрические слои располагаются вдоль сепаратрис, отделяющих пересоединившиеся магнитные силовые от непересоединившихся. Вне двойных слоев условие квазинейтральности плазмы выполняется с большой точностью. Исследованы эффекты, связанные с прохождением фронта падения давления вдоль магнитосферного хвоста. Показано, что резкий перепад давления на фронте приводит к резкому нарушению равновесия плазменного слоя и образованию перетяжек, приводящих к неустойчивости. Усовершенствована математическая модель проникновения электрических полей и токов в ионосферу с поверхности Земли. Проанализирована наиболее популярная в настоящее время модель такого проникновения, включающая сторонние токи. Показано, что она не применима в земных условиях. Блок 2. Разработка и отладка вычислительных алгоритмов, позволяющих менять параметры конструкции и режима работы элементов приборного блока космического аппарата во время проектирования (В. А. Деревянко, Е. Н. Васильев, А. В. Макуха, С. В. Кукушкин). Разработан вычислительный алгоритм, интерфейс и программный комплекс «Panel Emulator» для расчетов теплового режима приборного блока космического аппарата негерметичного исполнения. В «Panel Emulator» расчетная часть и интерфейс полностью разделены и реализованы в разных программных компонентах, связанных между собой с помощью XML-схемы. Такая реализация обеспечивает простоту адаптации пакета для решения других сходных вычислительных задач. Математическая модель, положенная в основу «Panel Emulator», учитывает циклограммы включения приборов и нагревателей, динамику изменения внешних радиационных потоков при движении по орбите. Модель позволяет задавать сложную геометрию приборной сотовой панели с отверстиями и выступами. Разработанный программный комплекс предоставляет конструктору удобный графический интерфейс, с помощью которого можно задавать конфигурацию приборов и элементов терморегулирования на поверхностях сотовой панели. «Panel Emulator» позволяет проводить для заданной конфигурации приборного блока моделирование теплового режима и наблюдать динамику изменения тепловых полей. Предусмотрена возможность отслеживания динамики изменения температуры в заданных точках поверхности и индикации перегрева приборов. Работа выполнена в интересах ОАО «Информационные спутниковые системы им. ак. М. Ф. Решетнева», г. Железногорск. Блок 3. Донное течение за летательным аппаратом в сверхзвуковом потоке (В. В. Шайдуров, Г. И. Щепановская, М. В. Якубович). Изучено донное течение за летательным аппаратом при сверхзвуковом обтекании вязким теплопроводным газом. Предложена математическая модель на основе полных уравнений Навье-Стокса для двумерного случая. Дискретизация по пространству осуществляется методом конечных элементов. Метод траекторий применяется при аппроксимации по времени, впервые для всей системы уравнений Навье-Стокса. Из тестовых расчётов следует, что совокупность метода конечных элементов и метода траекторий позволяет построить экономичный алгоритм с вычислительной точки зрения. На рис. III.15 показано распределение поля скоростей за летательным аппаратом. Донное течение, получено по расчётным значениям за летательным аппаратом, обтекаемым сверхзвуковым потоком вязкого теплопроводного газа. Числа Маха и Рейнольдса равны соответственно M=2 и Re=2000. За телом образуется вихревая зона с отрицательными значениями скорости. ![]() Рис. III.15. Распределение поля скоростей за летательным аппаратом Блок 4. Построение кинетической квазиклассической теории трехмерного охлаждения и устойчивого удержания атомов и ионов с триподной системой квантовых переходов (типа F=1 > F=0) в частично когерентном световом поле. I. Математическая модель локализации атомов в полихроматической оптической сверхрешетке (И. В. Краснов). Построена квазиклассическая модель локализации атомов и ионов с триподной конфигурацией квантовых переходов в трехмерной (3D) диссипативной оптической сверхрешетке нового типа. Ее действие основано на эффекте выпрямления градиентной оптической силы в полихроматическом поле: специальной суперпозиции трехчастотного когерентного поля и частично-когерентного поля. Показано, что движение атомов в такой оптической сверхрешетке описывается кинетическим уравнением типа Фоккера-Планка для вигнеровской функции распределения атомов в фазовом пространстве или эквивалентной ему системой стохастических дифференциальных уравнений, учитывающих квантовые флуктуации оптических сил. Получено приближенное стационарное решение кинетического уравнения для невзаимодействующих частиц (движущихся в полихроматической решетке) в виде разложения по малому параметру — числу Кнудсена в данной задаче: отношению эффективной длины свободного пробега атомов в «вязком газе фотонов» к периоду решетки. Это решение демонстрирует способность полихроматической сверхрешетки осуществлять эффективное охлаждение и сильную трехмерную пространственную локализацию триподных атомов и создавать кристаллоподобную систему локализованных сгустков атомов с хорошо управляемыми характеристиками (Рис. III.16, Рис. III.17). Другое, практически важное достоинство оптической сверхрешетки (вытекающее из решения кинетического уравнения): существование режимов охлаждения и локализации атомов, сочетающие в себе гигантские критические скорости (capture velocity) атомов, захватываемых в решетку с возможностью достижения субдоплеровских температур. Эти режимы локализации достигаются при медленном согласованном изменении интенсивностей полей (Рис. III.17). II. Моделирование формирования ионных кулоновских кластеров в полихроматической оптической решетке (Л. П. Каменщиков, И. В. Краснов). Ионные кластеры применяются в спектроскопии высокого разрешения, физике взаимодействующих кулоновских систем, квантовой информатике (хранение информации на квантовом уровне) и т.п. Цель работы этапа: обоснование с помощью численного моделирования возможности формирования ионных кластеров в оптическом (лазерном) поле. При этом следует решить две проблемы: 1) охлаждение ионов до достаточно низких температур; 2) добиться пространственной локализации частиц (с долговременным удержанием) с учетом силы кулоновского отталкивания. ![]() Рис. III.16. Типы пространственного распределения атомов в ячейках полихроматической сверхрешетки: (a) и (c) — пространственные распределения атомов в плоскости z=0; (b), (d) и (e) — области локализации атомов вблизи узлов сверхрешетки; (a), (b) — $\textit{L${}_{x}/\lambda$}$= $\textit{L${}_{y}/\lambda$ = L${}_{z}/\lambda$ = }$2.7$\times$10${}^{4}$${}_{ }$; (c), (d), (e) — $\textit{L${}_{x}/\lambda$}$ = 2.7$\times$10${}^{4}$; (d) — $\textit{L${}_{x}$/L${}_{y}$${}_{ }$ = L${}_{x}$/L${}_{z}$}$${}_{ }$ = 20; (e) — $\textit{L${}_{x}$/L${}_{y}$}$${}_{ }$ = 20, $\textit{L${}_{x}$/L${}_{z}$}$${}_{ }$ = 1; $\textit{s${}_{0}$/L ${}_{x}$} \approx$ 0.05, $\textit{s${}_{1}$/L ${}_{x}$} \approx$ 0.01. Низкие температуры получаются методом лазерного охлаждения ионов. Для пространственной локализации частиц ранее обычно использовались радиочастотные ловушки и статические электромагнитные ловушки (их недостаток: заметные паразитные движения). В 2010 г. была выдвинута идея локализации ионных частиц с помощью только оптических ловушек. В работах [1-2] показана перспективность освоения чисто оптических методов локализации ионов. В экспериментах Schneider'а и др. использовалась градиентная оптическая сила, порождаемая неоднородным монохроматическим лазерным излучением. В настоящей работе используем так называемые «выпрямленные градиентные силы» (rectified gradient forces) в полихроматическом лазерном поле с получением специальной пространственной конфигурации — полихроматической оптической сверхрешетки, предложенной работах [3-5]. Уравнения светоиндуцированного броуновского движения ансамбля взаимодействующих ионов в полихроматической оптической сверхрешетке выражают ньютоновский закон движения материальных частиц: \[\; \; \; \; \; \; \; \; \; \; \; \frac{d{\rm r}_{i} }{dt} ={\rm v}_{i} {\kern 1pt} ,\; \; \; \; \; m\frac{d{\rm v}_{i} }{dt} ={\rm F}_{i}^{trap} +{\rm F}_{i}^{laser} +{\rm F}_{i}^{Coul} +{\rm F}_{i}^{stoch} ,\; \; \; \; i=1,{\kern 1pt} {\kern 1pt} .{\kern 1pt} {\kern 1pt} .{\kern 1pt} {\kern 1pt} .{\kern 1pt} {\kern 1pt} ,N,\] где ${\rm r}_{i} $ и ${\rm v}_{i} $ — радиус-вектор и скорость $\textit{i}$-й частицы, ${\rm F}_{i}^{trap}, \; \; {\rm F}_{i}^{laser} $ — консервативная и диссипативная (трение) составляющие выпрямленной градиентной силы, ${\rm F}_{i}^{Coul} $ — кулоновская сила, действующая на $\textit{i-}$й ион со стороны всех других ионов, ${\rm F}_{i}^{stoch} $ — случайная сила (типа $\delta $-коррелированного гауссовского случайного процесса), обусловленная квантовыми флуктуациями и приводящая к диффузии в пространстве скоростей. ![]() Рис. III.17. Зависимость относительной температуры атомов $\textit{T'}$ от управляющих параметров — безразмерных интенсивностей полей $\textit{c}$ и $\textit{G}$; $\textit{A}$, $\textit{B}$ — характерные области, соответствующие супердоплеровским температурам ($\textit{T'}$>>1), $\textit{C}$ — область субдоплеровских температур ($\textit{T'}$<1); $\textit{G}$ и $\textit{G'}$ — траектории устойчивого (без потери локализованных частиц) и неустойчивого перехода $\textit{A$\rightarrow$C}$. Область $\textit{A}$ характеризуется гигантскими значениями критической скорости частиц (velocity capture), при которой они могут захватываться в оптическую сверхрешетку Численный пример. Проведен расчет движения ионов химического элемента иттербия в оптической решетке $3\times 3\times 3$ ячеек. В начальный момент в каждой ячейке среднего слоя (в плоскости $z=0$) были случайно разбросаны по шесть частиц с нулевыми начальными скоростями (рис. III.18 (а)), т.е. рассматривался ансамбль ионов из $N=54$ частиц. На рис. III.18 (б) показано положение частиц через $1.4{\kern 1pt} \cdot {\kern 1pt} 10^{-5} $c, а на рис. III.18 (в) — через $5.32{\kern 1pt} \cdot {\kern 1pt} 10^{-5} $с. При дальнейшем счете координаты частиц заметно не менялись, т.е. частицы достигали устойчивого состояния равновесия. Для решения системы дифференциальных уравнений применялся численный метод Штёрмера-Верле. ![]() Рис III.18. Положение ионов в решетке в разные моменты времени Построена вычислительная модель, описывающая формирование ионных кулоновских кластеров в полихроматической оптической сверхрешетке. Предварительные численные расчеты продемонстрировали принципиально новую возможность создания упорядоченного массива ионных кулоновских кластеров с помощью оптических сил в лазерном поле. В отличии от традиционных схем создания кластеров в радиочастотных ловушках и ловушках Пеннинга в предложенном методе отсутствуют паразитные движения типа микроосцилляций (с заметной энергией) или циклотронные и магнетронные вращательные движения. Отметим также, что построенная математическая модель позволяет рассчитывать реальные режимы формирования двумерных и трехмерных ионных кулоновских кластеров. Основные публикации:
III. Численный анализ воздействия лазерного излучения на агрегат наночастиц (А. П. Гаврилюк). На основе ранее построенной математической модели (см. A. P. Gavrilyuk, S. V. Karpov. Applied Physics B: v.97, issue 1, pp. 1 1) при интенсивностях, при которых начинают проявляться ПИФ (поляризационно-индуцированные силы), существенен нагрев электронов и кристаллической решетки частиц. Это означает, что необходимо в этих условиях учитывать термоэмиссию электронов из частиц димера; 2) фотомодификация димеров при умеренных плотностях энергии излучения, в основном, связана со смещением частиц внутри домена относительно друг друга. Это смещение обусловлено действием ПИФ и (или) плавлением полимерного слоя, происходящим при нагреве частиц; 3) чем меньше межчастичная щель димера, тем значительнее изменение его спектра при воздействии квазирезонансного лазерного излучения, то есть, меньше энергетический порог фотомодификации. Основные публикации: К блоку 1:
К блоку 2:
К блоку 3:
К блокам 4:
(Отделы Вычислительной математики, Вычислительной физики)
|
Webmaster |