//Росатом/ ВНИИЭФ
 
Главная / Издания /"ВАНТ" Серия: Математическое моделирование физических процессов /2010 год /Выпуск 3/2010 /

Выпуск 3/2010

АННОТАЦИИ:


УДК 517.958:536.2
ПРИМЕНЕНИЕ TVD-ПОДХОДА К РЕШЕНИЮ УРАВНЕНИЯ ПЕРЕНОСА ТЕПЛОВОГО ИЗЛУЧЕНИЯ КВАЗИДИФФУЗИОННЫМ МЕТОДОМ / А. Д. Гаджиев, С. А. Грабовенская, В. В. Завьялов, А. А. Шестаков // С. 3-14.
Рассматривается нелинейная схема типа TVD второго порядка аппроксимации для решения уравнения переноса теплового излучения квазидиффузионным методом. Приводятся результаты численных расчетов тестовых задач.

Ключевые слова:  схема TVD, квазидиффузионный метод.


УДК 517.958:536.2

МЕТОДЫ ВЫДЕЛЕНИЯ ДИАГОНАЛЬНОГО ЭЛЕМЕНТА   ДЛЯ УСКОРЕНИЯ ИТЕРАЦИЙ ПРИ ЧИСЛЕННОМ РЕШЕНИИ УРАВНЕНИЯ ПЕРЕНОСА ТЕПЛОВОГО ИЗЛУЧЕНИЯ В ДИФФУЗИОННЫХ ПРИБЛИЖЕНИЯХА. Д. Гаджиев, А. А. Шестаков // С. 15-27.
Рассматриваются численные методы выделения диагонального элемента для ускорения итераций при решении уравнения переноса теплового излучения в диффузионных приближениях. Основную идею методов можно реализовать в три этапа:
1. Для каждой счетной ячейки сначала решается уравнение энергии, из которого по методам выделения диагонального элемента находится температура через плотность излучения.
2. По найденной температуре вычисляется функция Планка.
3. Новые значения плотности и потока излучения находятся прогонкой при известной правой части.

Ключевые слова: ускорение итераций, перенос теплового излучения.


УДК 519.6

DSn-СХЕМА ХАРАКТЕРИСТИЧЕСКИХ ТРУБОК ДЛЯ РЕШЕНИЯ УРАВНЕНИЯ ПЕРЕНОСА, ЗАПИСАННОГО В ЦИЛИНДРИЧЕСКОЙ СИСТЕМЕ КООРДИНАТА. И. Бочков // С. 28-37.
Для уравнения переноса, записанного в цилиндрической системе координат, построена характеристическая схема. Схема включает в себя два этапа. На первом этапе в счетной ячейке строится обыкновенное дифференциальное уравнение баланса частиц, с помощью которого находится полный поток частиц на неосвещенной поверхности ячейки. На втором этапе по схеме типа DSn-метода происходит распределение выходящего потока частиц по неосвещенным граням ячейки.

Ключевые слова: уравнение переноса, разностная схема, DSn-метод, DSn-схема характеристических трубок.


УДК 517+519+533

ОСОБЕННОСТИ СЧЕТА ДЕТОНАЦИИ НА ЭЙЛЕРОВЫХ  СЕТКАХ С НАВЯЗАННОЙ СКОРОСТЬЮ ФРОНТАЮ. А. Бондаренко // С. 38-45.
Анализируются результаты одномерных газодинамических расчетов нормальной детонации на эйлеровых сетках. Показано существенное влияние небольшого (~1%) несохранения полной энергии на точность счета детонации. Построено автомодельное решение, учитывающее наличие отрицательного дисбаланса полной энергии при навязанной скорости фронта волны энерговыделения. Результаты расчетов по неконсервативным методам согласуются со свойствами этого автомодельного решения. Продемонстрировано, что использование более точных методов пересчета газодинамических величин с лагранжевой сетки на эйлерову позволяет существенно экономить вычисления при сохранении точности.

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


УДК 519.6

ДВУМЕРНАЯ СКАЛЯРНАЯ ИСКУССТВЕННАЯ ВЯЗКОСТЬ МЕТОДИКИ ЭГАК В СФЕРИЧЕСКИХ СИСТЕМАХЮ. В. Янилкин, О. О. Топорова // С. 46-54.
Рассматривается задача о сферическом схождении оболочки,  внутри которой содержится легкий газ. Показано, что при численном решении этой задачи важное значение имеет форма искусственной вязкости. На сферических сетках классическая вязкость Неймана–Рихтмайера может давать большую погрешность, связанную с тем, что в ней не учитывается однородное сжатие вещества при сферическом схождении. Значительного уменьшения этой погрешности в данной работе удалось добиться применением ранее предложенной авторами вязкости.

Ключевые слова: газодинамика, искусственная вязкость, двумерные течения, разностная схема.


УДК 519.6

КИНЕТИКА ИНИЦИИРОВАНИЯ ВЗРЫВЧАТЫХ ПРЕВРАЩЕНИЙ ЗА СЧЕТ ТРЕНИЯ ПРИ НИЗКОСКОРОСТНЫХ ВОЗДЕЙСТВИЯХ В РАМКАХ МЕТОДИКИ Д / М. Е. Мартынова, А. Р. Гушанов, Ю. Д. Чернышев // С. 55-64.
Рассматривается один из возможных механизмов инициирования взрывчатых превращений при низкоскоростных воздействиях. Он основан на поверхностном нагреве взрывчатого вещества при его трении об окружающие материалы. Представлена реализация уравнений кинетики в рамках методики Д. Описан алгоритм учета сил трения, возникающих на контактных границах при взаимодействии взрывчатого вещества с окружающими веществами. Приведены результаты численного моделирования экспериментов в постановке Stiven Test с использованием описанной кинетики.

Ключевые слова: методика Д, низкоскоростные воздействия, Stiven Test, коэффициент трения, инициирование взрывчатых превращений.


УДК 519.6

ВЕРИФИКАЦИЯ МОДЕЛИ РАСЧЕТА КОНВЕКТИВНОГО  ПЕРЕНОСА ПРИМЕСИ НА ПРИМЕРЕ КРИСТАЛЛИЗАЦИИ  ДВУХКОМПОНЕНТНЫХ РАСПЛАВОВ Pb-Sn / С. М. Ганина, В. П. Гинкин // С. 65-72.
Разработан алгоритм расчета сегрегации примеси при кристаллизации двухкомпонентных расплавов, характеризующихся большими размерами двухфазной области. При расчете конвекции в пористой среде учитывается добавочное сопротивление течению по аналогии с законом Дарси. Температура ликвидус линейно зависит от концентрации примеси в жидкой фазе.
Предложенная модель реализована в программном комплексе GIGAN-2D. Для верификации модели просчитаны два численных бенчмарка: по кристаллизации расплава Sn c 10 %-ным содержанием Pb и расплава Pb с 18 %-ным содержанием Sn.

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


УДК 519.6

АЛГОРИТМ МОДЕРНИЗИРОВАННОГО МНОГОЧАСТИЧНОГО ПОТЕНЦИАЛА ДЛЯ МОЛЕКУЛЯРНО-ДИНАМИЧЕСКОГО МОДЕЛИРОВАНИЯ НА ГРАФИЧЕСКОМ АРИФМЕТИЧЕСКОМ УСКОРИТЕЛЕ / С. В. Копкин, И. А. Крючков // С. 73-82.
Представлены основные этапы адаптации алгоритма вычисления значений потенциала МЕАМ в комплексе МД для гибридной вычислительной системы. Показано ускорение, достигаемое на различном количестве вычислительных устройств. Проведены тестовые расчеты, подтвердившие точность вычислений.
Ключевые слова: арифметический ускоритель, молекулярная динамика, потенциал, MEAM, ускорение, гибридная вычислительная система.


Сведения об авторах

© 2011 ФГУП РФЯЦ-ВНИИЭФ