МОДЕЛИРОВАНИЕ СТРУКТУР АМОРФНЫХ СПЛАВОВ СРЕДСТВАМИ ГРАФИЧЕСКИХ ПРОЦЕССОРОВ (GPGPU)

Научная статья
Выпуск: № 11 (18), 2013
Опубликована:
2013/12/08
PDF

Полищук С.В.1, Смехун Я.А.2

1Магистрант, 2магистрант, Дальневосточный федеральный университет

МОДЕЛИРОВАНИЕ СТРУКТУР АМОРФНЫХ СПЛАВОВ СРЕДСТВАМИ ГРАФИЧЕСКИХ ПРОЦЕССОРОВ (GPGPU)

Аннотация

Данная работа посвящена моделированию атомных структур аморфных сплавов средствами GPU. Модель аморфного сплава представляется как набор атомов. Моделирование производится в двух вариантах: двумерном, как электронно-микроскопическое изображение, и трехмерном, как 3d-модель атомов сплава. В первом, случае высокая скорость расчета кадра достигается путем использования возможностей параллельных вычислений, с помощью программного интерфейса CUDA. Для визуализации объемных моделей используется Vertex Array Object (VAO) – расширение OpenGL.

Ключевые слова: моделирование структур, аморфные сплавы, графические процессоры.

Polischuk S.V.1, Smekhun Y.A.2

1Undergraduate, 2undergraduate, Far Eastern Federal University

MODELING THE STRUCTURE OF AMORPHOUS ALLOYS BY MEANS OF GRAPHICS PROCESSORS (GPGPU)

Abstract

This work is devoted to modeling the atomic structures of amorphous alloys by means of GPU. The model of the amorphous alloy is represented as a set of atoms. There are two variants of modeling: it is a two-dimensional variant, where images are as electron microscopic and three-dimensional images as the 3d- model atoms of the alloy. In the first case, the high speed calculation frame is achieved through the use of parallel computing, using the API of CUDA. For visualization of three-dimensional models used Vertex Array Object (VAO) - extension of OpenGL.

Keywords: modeling of structures, amorphous alloys, graphics processors.

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

Среди наиболее эффективных методов моделирования и визуализации выступает использование графического процессора (GPGPU) в качестве вычислителя.

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

Так как необходимо визуализировать статичные данные: координаты, цвета и индексы, которые в процессе моделирования не изменяются, то наиболее оптимальным решением является использование ARB_vertex_buffer_object (VBO) – расширения OpenGL.

Для упрощения процесса кодирования целесообразно использовать следующее расширение OpenGL – ARB_vertex_array_object (VAO), которое позволяет упростить процедуры активации буферов до одной команды.

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

Построение 2D моделей. Моделирование средствами GPGPU [2] производится в несколько этапов [3]. Производится вычисление проекций потенциалов отдельных атомов (1) и последующее формирование потенциала модели путем линейной суперпозиции потенциалов каждого атома в данной модели (2):

12-05-2021 11-36-23        (1)

где 12-05-2021 11-36-32 – положение атома в плоскости, перпендикулярной к оптической оси микроскопа; 12-05-2021 11-36-40 – модифицированная функция Бесселя нулевого порядка; 12-05-2021 11-36-50 – радиус орбиты Бора для электрона; 12-05-2021 11-37-03 – коэффициенты, значения которых зависят от атомного номера и приводятся в соответствующих таблицах [4].

12-05-2021 11-37-15        (2)

Далее необходимо вычислить функцию прохождения электронов через исследуемый объект (3):

12-05-2021 11-37-26      (3)

где σ – постоянная взаимодействия электронного пучка с образцом.

Получаем преобразование Фурье от функции прохождения (4).

12-05-2021 11-37-54        (4)

Умножаем функцию (4) на передаточную функцию объектива (5), чтобы получить волновую функцию в обратной фокальной плоскости (6).

12-05-2021 11-38-07        (5)

где 12-05-2021 11-41-46 - дефокусировка; 12-05-2021 11-41-55 - коэффициент сферической аберрации; 12-05-2021 11-42-03 - функция апертуры.

12-05-2021 11-41-28         (6)

Вычисляем квадрат модуля волновой функции изображения (в реальном пространстве), чтобы получить конечную интенсивность (7).

12-05-2021 11-41-37      (7)

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

12-05-2021 11-45-38

Рис. 1 - Электронно-микроскопические изображения: а) реальное; б) смоделированное

3D визуализация. Перед вычислением проекции потенциалов модели, необходимо убедиться в корректности дефектов внесенных в модель, так как процесс создания проекций очень затратный по времени и ресурсам. Визуализация реализуется средствами OpenGL 2.1 с помощью расширения ARB_vertex_array_object (VAO). Данные вершин, цветов и индексов, считываются, полностью загружаются в память видеоустройства, и удаляются из системной памяти. Расчеты производились на станции NVIDIA GTX 690. Примерное время обработки модели 51*51*51 нм составляет 2-3 минуты. Расчет проекций потенциалов для данной модели составляет около 5-6 минут.

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

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

  • Кристаллография: Лабораторный практикум / Под ред. Проф. Е.В. Чупрунова: Учеб. Пособие для вузов. – М.: Издательство физико-математической литературы, 2005. – 412 с. – ISBN 5-94052-103-7

  • David B. Kirk, Wen-mei W. Hwu. Programming Massively Parallel Processors, 2012. – p. 514.

  • Earl J. Kirkland. Advanced Computing in Electron Microscopy, 2010. – p. 289.

  • A. Doyle and P.S. Turner. Relativistic Hartree-Forck x-ray and electron scattering factors. ActaCryst., 1968, р. 390-397.