Астрофизику сегодня невозможно представить без компьютерного моделирования: ученые воссоздают на ЭВМ космические процессы, не доступные для наблюдения, чтобы ставить эксперименты и подтверждать теории. Математики из новосибирского Академгородка пишут программы, которые решают такие задачи в несколько раз эффективнее, чем зарубежные аналоги. Недавно статья об этом вышла в сборнике Supercomputing издательства Springer.

«Задачами астрофизики я занимаюсь с самого начала своей научной жизни, уже примерно 15 лет. За это время мы моделировали самые разные объекты: от звезд и галактик до космической паутины крупномасштабной структуры Вселенной», — рассказывает старший научный сотрудник Института вычислительной математики и математической геофизики СО РАН доктор физико-математических наук Игорь Михайлович Куликов. Сейчас он руководит проектом по суперкомпьютерному моделированию сверхновой типа 1а.
 
Постановка задачи
 
Сверхновая — это заключительный этап в жизни некоторых звезд; взрыв, самый сильный, какой только можно наблюдать в космосе. При этом звезда резко становится ярче на 4—8 порядков, а потом постепенно гаснет. Во время взрыва сверхновой выделяется колоссальное количество энергии, и вещество из внешней оболочки звезды распространяется на дальние расстояния, а из остатков образуются или нейтронная звезда, или черная дыра.
 
При взрыве сверхновой разлетаются продукты термоядерного синтеза, происходившего на протяжении существования звезды, — все жизненно важные химические элементы тяжелее гелия. Это причина эволюции химических веществ: из остатков сверхновых состоит практически все во Вселенной, включая людей. 
 
Ученые исследуют излучение, которое достигает Земли через много лет после образования сверхновой, остатки от взрыва, а также потенциальные звезды-предшественники. На основании этих данных пишут уравнения и строят компьютерные модели, позволяющие наблюдать и исследовать процессы появления сверхновых в режиме реального времени.
 
«Сверхновая типа 1а рождается в результате взрыва белого карлика — старой звезды, от которой осталось только ядро. У белого карлика есть определенный предел массы, при котором он может существовать, — предел Чандрасекара, равный примерно 1,44 солнечных масс. Если масса белого карлика становится больше, его ядро начинает сжиматься, температура растет, запускается термоядерная реакция. В итоге — взрыв сверхновой, — говорит Игорь Куликов. — Это происходит в том случае, если у белого карлика есть звезда-компаньон, с которой он может слиться, и в результате увеличить свою массу».
 
Один из сценариев образования сверхновой типа 1а — это так называемый сценарий не центрального взрыва. «Предположительно, белый карлик может взорваться до того, как полностью срастется со своим компаньоном. В таком случае взрыв происходит на периферии, и остатки звезды будут распространяться преимущественно в одном направлении», — поясняет ученый. Эту идею высказал главный научный сотрудник Института астрономии РАН доктор физико-математических наук Александр Васильевич Тутуков, а коллектив ИВМиМГ СО РАН создал компьютерную модель такого сценария. Недавно она нашла подтверждение: в ноябре 2018 года в остатках взрыва сверхновой обнаружили остатки звезды-компаньона
 
Ресурсы
 
Космические вычисления требуют мощной техники. Игорь Куликов, сидя в своем кабинете, использует ноутбук для доступа к суперкомпьютерам, расположенным несколькими этажами ниже. «Современные технологии позволяют работать на расстоянии. Многие научные сотрудники удаленно заходят на наш кластер, передают данные, задача считается, и полученный результат они уже анализируют на своих персональных компьютерах. Можно отправлять задачи в наш центр, находясь как на рабочем месте, так и в командировке или на отдыхе», — рассказывает руководитель Сибирского суперкомпьютерного центра ИВМиМГ СО РАН кандидат физико-математических наук Игорь Геннадьевич Черных. «Это Игорь нас имеет в виду, — улыбается его тезка. — Мы оба даже в отпуск берем рабочие ноутбуки».
 
Основных вычислительных кластеров в суперкомпьютерном центре два. Новый и пока еще небольшой НКС-1П — современная вычислительная система, мощностью 100 терафлопсов. 
 

Новый кластер потребляет совсем не много электроэнергии: вычислительная часть — 15 кВт, а система охлаждения — всего 500 Вт (в четыре раза меньше, чем электрический чайник). НКС-1П построен на базе серверных процессоров Intel Xeon и специализированных процессоров Intel Xeon Phi. Каждый процессор содержит 72 ядра. 

 
«Старичок» НКС-30Т 2012 года не настолько энергоэффективен. С кондиционерами этот компьютер потребляет порядка 100 кВт. Он производит довольно много шума, а из системы охлаждения вырывается мощный поток воздуха. Но сегодня и без этой машины не обойтись.
 
«Когда мы ввели в строй новый кластер, то думали, что на год-два немного снизим нагрузку старого компьютера, но уже через две недели поняли, что надо еще увеличивать вычислительную мощность, потому что сразу захотелось считать более сложные задачи. Современные исследования уровня мегасайнс должны производиться на оборудовании класса мегасайнс, как минимум в пять раз мощнее того, что есть у нас», — говорит Игорь Черных. 
 
Сейчас совокупная мощность суперкомпьютеров центра — порядка 0,2 петафлопса. Для сравнения: пиковая вычислительная мощность самого производительного на сегодняшний день суперкомпьютера Summit (США) составляет 200 петафлопсов.
 
 
Технологии
 
Однако преимущество сибирских ученых в первую очередь не в технике, а в головах: они пишут код, который позволяет максимально использовать возможности процессоров.
 
«Мы применяем векторизацию. Это увеличивает производительность кода примерно в 10—20 раз, — объясняет Игорь Куликов. — Грубо говоря, если считать задачу лучшими зарубежными или российскими программами, то это займет примерно неделю или десять дней, а такие же вычисления с использованием нашего кода — около одного дня. На данный момент в мире мы единственные: для задач астрофизики нет других программ, которые используют векторизацию».
 
Векторизация — это вид параллельных вычислений, то есть таких, когда программа выполняет несколько вычислительных процессов одновременно. Суперкомпьютеры оснащены многоядерными процессорами, поэтому, чтобы использовать их эффективно, последовательной программы недостаточно: нужно распределить работу между ядрами процессора. Приложения, которые выполняют операции последовательно, при векторизации модифицируются таким образом, чтобы выполнять одновременно несколько однотипных операций. 
 
«Представьте, что есть два набора из четырех элементов: четыре цифры с одной стороны и четыре с другой. Нам нужно поэлементно их сложить. В обычных, последовательных вычислениях мы сначала складываем и записываем первый результат в первую ячейку, потом второй во вторую, третий в третью, четвертый в четвертую. У нас получается четыре инструкции вычисления (сложения) и четыре инструкции записи, то есть всего восемь инструкций. Векторизация позволяет это сделать за две инструкции. За одну — сложить и за одну — записать. А векторы, с которыми мы работаем на процессоре Intel Xeon Phi, — это уже даже не четыре, а восемь элементов двойной точности», — рассказывает Игорь Куликов.
 
Ученые начинают писать программу с уровня математической реализации, с уравнений для той или иной задачи. Игорь Куликов поясняет: «Еще на этапе постановки задачи, когда мы только записываем уравнение в тетрадке, мы уже должны понимать, на чем будем считать. Разработка конечного программного обеспечения сразу должна учитывать архитектуру, то есть общее устройство кода. Мы называем это содизайном (co-design)».
 
Для расчетов используются уравнения гидродинамики — раздела физики, который изучает движение идеальных жидкости и газа. Это может показаться странным, ведь межзвездное пространство характеризуется очень низкой плотностью и давлением и приближено к абсолютному вакууму. «В космосе тоже есть газ. В среднем плотность галактики — примерно один атом водорода на 1 см3. Это, конечно, крайне мало, но пока атом из одной части галактики перелетит в другую, он столкнется с другими атомами много раз. Если спроектировать эти столкновения на меньший объем, объем обычного воздуха, то количество атомов и количество столкновений будет сопоставимо с такими параметрами у идеального газа. Поэтому в астрофизических масштабах представление всей видимой материи как идеального газа — это вполне адекватная модель», — объясняет Игорь Куликов. 
 
Реализация
 
Решить уравнения для такой большой области, как в задачах астрофизики, на сегодняшний день невозможно обычным, аналитическим способом. Поэтому ученые применяют численный метод, дискретизируя интересующую их область — представляют ее как сетку и решают уравнения для каждой из ее ячеек.
 
«Как происходят вычисления: мы берем какую-то область пространства, делим ее по всем направлениям на кубики, и каждый кубик начинаем обсчитывать по нашим уравнениям. Известно, что если производительность, память, энергоэффективность, программное обеспечение компьютеров будут развиваться так же, как в последние десятилетия, то компьютер, достаточно мощный для того, чтобы при моделировании всей Вселенной достичь разрешения, когда один кубик по размеру равен такой звезде, как Солнце, появится примерно через 60 лет. И даже такого разрешения не всегда достаточно. Например, в нашей задаче взрыва сверхновой типа 1а фронт горения углерода очень мал по сравнению с самой звездой», — говорит Игорь Куликов.
 
Чтобы преодолеть этот барьер, ученые, во-первых, используют вложенные сетки — получают как бы кубик в кубике: дробят участки на более мелкие в тех местах, где необходимы более точные вычисления (в областях горения углерода и фронта ударной волны). А во-вторых, разработали метод с рабочим названием «спутниковые вычисления».
 
«Мы вырезаем тот маленький кубик, который выделили внутри кубика, расширяем его фактически на отдельную задачу, на новую сетку, и снова запускаем суперкомпьютер. А затем полученный результат возвращаем в основную сетку. Это позволило смоделировать ядерное горение углерода, а также достаточно хорошо оценить выделившуюся энергию. И мы достаточно адекватно смогли воспроизвести сам взрыв сверхновой типа 1а», — говорит Игорь Куликов.
 
Математики планируют моделирование и других типов сверхновых. В частности, они уже начали заниматься коллапсом молекулярного облака, который в некоторых случаях приводит к гравитационному коллапсу и рождению новых звезд. По словам ученых, эта задача предшествует моделированию образования сверхновых второго типа и сверхновых типа 1b и 1c, основанному на коллапсе ядра.
 
«Задачи астрофизики, как мне кажется, — одни из самых интересных научных задач. В конечном итоге, они призваны ответить на вопрос, который всегда волновал человечество: откуда мы появились и как, а также что будет со Вселенной дальше. Поэтому нам интересна эволюция химических веществ: хочется понять те факторы, которые способствуют возникновению простейшей органики, основы жизни. Благодаря суперкомпьютерным вычислениям появилась возможность обрабатывать большие данные, полученные с современных радиотелескопов, с отдаленных галактик, искать закономерности и использовать эти данные для моделирования», — говорит Игорь Черных.
 
«Мы не претендуем на открытия в астрофизике, у нас другое направление — высокопроизводительные вычисления, но все расчеты мы делаем качественно: от записи уравнения до решения конкретных задач на компьютере», — добавляет Игорь Куликов.
 
Сибирский суперкомпьютерный центр ИВМиМГ СО РАН для научной работы используют около 200 человек из 30 организаций, главным образом Сибирского отделения РАН. Все институты новосибирского Академгородка подсоединены к суперкомпьютерам 10-гигабитной сетью. Возможна удаленная работа из других городов и стран.
 
Основные пользователи центра — химические институты. Почти половину вычислительного времени потребляет ФИЦ «Институт катализа им. Г.К. Борескова СО РАН». Здесь используют квантово-химические вычисления в разработке новых катализаторов, а также применяют численное моделирование химических реакторов, которые планируется построить.
 
Много ресурсов требуют задачи в области нефтегазовой геологии и геофизики, гидродинамики, теоретической и прикладной механики, прогнозирования различных природных явлений (например, цунами), а также генетики.
 
В рамках проекта «Академгородок 2.0» на базе ССЦ планируется создать Сибирский национальный центр высокопроизводительных вычислений и обработки данных (СНЦ ВВОД) с мощностью минимум 10 петафлопсов. Это совместный проект ИВМиМГ СО РАН, Института вычислительных технологий СО РАН и Новосибирского государственного университета. Конец вреза. 
 
Работа выполняется при поддержке гранта РНФ 18-11-0044.
 
Александра Федосеева
 

Источники

Вычислить путь звезды
Академгородок (academcity.org), 12/04/2019
Ученые впервые в истории получили снимки черной дыры
Красноярский научный центр Сибирского отделения Российской академии наук (ksc.krasn.ru), 12/04/2019
Вычислить путь звезды, или что программисты делают в астрофизике
Наука в Сибири (sbras.info), 12/04/2019
Ученые из РФ создали алгоритм по расчету моделей рождения и смерти звезд
INTERFAX.RU, 12/04/2019
Ученые из РФ создали алгоритм по расчету моделей рождения и смерти звезд
The world news (theworldnews.net), 12/04/2019
Ученые из РФ создали алгоритм по расчету моделей рождения и смерти звезд
NewsRbk.ru, 12/04/2019
"Алтайские студенты участвуют в космическом эксперименте"
Ivest.kz, 12/04/2019
Ученые из РФ создали алгоритм по расчету моделей рождения и смерти звезд
События дня (inforu.news), 12/04/2019
Ученые из РФ создали алгоритм по расчету моделей рождения и смерти звезд
Спутник Новости (news.sputnik.ru), 12/04/2019
Новосибирские ученые создали алгоритм, позволяющий быстро рассчитывать модели рождения и смерти звезд
Новосибирские новости (nscn.ru), 12/04/2019
В Новосибирске нашли быстрый алгоритм для моделирования взрывов звезд
Российская газета (rg.ru), 12/04/2019
Вычислить путь звезды: что программисты делают в астрофизике
Pcnews.ru, 12/04/2019
Вычислить путь звезды: что программисты делают в астрофизике
Популярная механика (popmech.ru), 12/04/2019
В Новосибирске нашли быстрый алгоритм для моделирования взрывов звезд
Novosibirsk.4geo.ru, 12/04/2019
Вычислить путь звезды: что программисты делают в астрофизике
Все о космосе (aboutspacejornal.net), 13/04/2019
Вычислить путь звезды: что программисты делают в астрофизике
Новости для гиков (supreme2.ru), 15/04/2019
Вычислить путь звезды, или что программисты делают в астрофизике
Российский научный фонд (рнф.рф), 15/04/2019
Вычислить путь звезды, или что программисты делают в астрофизике
Российский научный фонд (rscf.ru), 15/04/2019

Похожие новости

  • 21/05/2019

    По итогам сочинского форума «Наука будущего — наука молодых»

    ​В Сочи завершились III Международная конференция «Наука будущего» и IV Всероссийский форум «Наука будущего — наука молодых». Мы попросили сибирских ученых, в них участвующих, рассказать, какие проекты они представляли на мероприятиях форума и с какими целями приехали сюда.
    437
  • 16/04/2019

    Восемь ответов на частые вопросы о СНЦ ВВОД

    Зачем нужен Сибирский национальный центр высокопроизводительных вычислений, обработки и хранения данных — СНЦ ВВОД? Откуда придут деньги на его создание? Как этот проект связан с синхротроном СКИФ? С другими проектами «Академгородка 2.
    511
  • 21/01/2019

    «Сотканные» в НМИЦ им. академика Е.Н. Мешалкина тканеинженерные протезы меняют идеологию операций

    В школьные годы я зачитывалась книгой кардиохирурга Николая Амосова “Мысли и сердце”. Воображение легко рисовало бригаду врачей над неподвижным телом пациента со вскрытой грудной клеткой. Но в ближайшем будущем операции по имплантации, например, аортального клапана, пораженного стенозом, будут выглядеть совсем по-другому.
    1026
  • 28/02/2019

    СКИФ обретает очертания

    В Новосибирске полным ходом идет проектирование уникального синхротрона четвертого поколения, который должны построить в рамках реализации проекта «Академгородок 2.0» к 2024 году. Ученые разработали эскиз первых шести пользовательских станций СКИФаНапомним, центр коллективного пользования СКИФ будет включать в себя, помимо собственно источника фотонов, пользовательское оборудование экспериментальных станций и лабораторного комплекса.
    540
  • 20/06/2017

    Международная выставка «НТИ ЭКСПО» в Новосибирске

    ​​​Уникальная международная выставка достижений технологического развития "НТИ ЭКСПО" пройдет в рамках V Международного форума технологического развития "Технопром-2017" 20-22 июня в Новосибирске при поддержке правительства РФ, коллегии ВПК, Минпромторга России, Минэкономразвития России, МИДа РФ, правительства Новосибирской области.
    2918
  • 14/03/2017

    Новосибирские ученые оценят опасность цунами на российских побережьях

    Берега России, как известно, омываются несколькими морями. Если за северные можно быть спокойными, то восточные и южное Черное способны преподнести неприятные и сверхопасные сюрпризы. Сибирские ученые создали карту цунами-опасности в регионах РФ и выяснили: даже «самое синее в мире» время от времени напоминает — оно не только приятное место отдыха, но и непредсказуемая стихия.
    1501
  • 23/01/2019

    Новосибирские физики смоделировали атмосферу экзопланет

    ​Сотрудники Института лазерной физики СО РАН в лабораторных условиях моделируют плазменный ветер, аналогичный тому, что испускают объекты в сотнях световых лет от Земли. Эти исследования имеют большое значение для изучения состава и динамики верхней атмосферы разных классов экзопланет, в том числе потенциально пригодных для жизни.
    810
  • 12/02/2019

    Ученые выяснили, в каких регионах России выгодны тепловые установки на энергии солнца

    Ученые из России и Италии рассчитали, в каких регионах РФ и для каких нужд выгодно использовать преобразователи тепла, работающие на солнечной энергии. Оказалось, что летом такие установки могут нагревать воду для душа, стирки и других бытовых потребностей по всей России, даже в Оймяконе, сообщила во вторник пресс-служба Российского научного фонда (РНФ), при поддержке которого проводилось исследование.
    510
  • 04/10/2017

    Ученые СО РАН ищут новые виды генерации

    ​​​Спрос на альтернативную энергетику в регионах страны неизменно растет. Успешные примеры внедрения инновационных проектов появляются и в Новосибирской области. ​В столице Сибири активно обсуждается идея создания энергетического технопарка на основе инновационных разработок и продукции научно-образовательных учреждений и энергетических компаний региона.
    891
  • 03/02/2018

    Ученые новосибирского Академгородка представили новейшие достижения СО РАН

    ​​Перед Днем российской науки-2018 три крупнейших института СО РАН – Институт ядерной физики им. Будкера, Институт химической биологии и фундаментальной медицины и Институт гидродинамики им. Лаврентьева  – открыли свои двери для посетителей.
    2439