Построение параллельных вычислительных алгоритмов высокого порядка точности для уравнений газовой динамики тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Жалнин, Руслан Викторович

• Жалнин, Руслан Викторович
• кандидат физико-математических наук
• 2008

Оглавление диссертации кандидат физико-математических наук Жалнин, Руслан Викторович

Введение

1 Алгоритм расчета газодинамических течений

1.1 Разностная схема.

1.2 Вычисление дискретных потоков.

1.3 Алгоритм реконструкции газодинамических параметров

1.3.1 ЕШ' алгоритм.

1.3.2 алгоритм.

1.4 Дискретизация по времени.

1.5 Алгоритм расчета для газов с различными показателями адиабаты

1.6 Решение тестовой одномерной задачи о распаде произвольного разрыва.

2 Программный комплекс СВЛУЕГ^ОРАНЗБ для параллельных вычислений

2.1 Параллельный вычислительный алгоритм.

2.2 Описание программного комплекса.

2.2.1 Общий алгоритм работы.

2.2.2 Подпрограмма инициализации StartUp

2.2.3 Подпрограмма реализации одного вычислительного шага.

2.2.4 Подпрограмма сохранения результатов расчетов Бауе-Data.

2.2.5 Подпрограммы реализации межпроцессорного обмена и граничных условий (ExchangeBoundCond, ВоипсЮопс!)

2.2.6 Модуль декомпозиции и «склеивания» расчетных подобластей.

2.2.7 Модули интерпретации и визуализации результатов расчетов.

2.3 Тестирование работы программного комплекса.

3 Прямое численное моделирование развития неустойчивости Рихтмайера-Мешкова

3.1 Постановка задачи.

3.2 Постановка расчетной задачи.

3.3 Сравнение результатов ЗБ расчетов с экспериментальными данными.

3.3.1 Расчеты Р1 и Р2.

3.3.2 Расчет РЗ.

3.4 Анализ течения на основе результатов ЗБ расчетов. Сравнение результатов.

Введение диссертации (часть автореферата) на тему «Построение параллельных вычислительных алгоритмов высокого порядка точности для уравнений газовой динамики»

Уравнения газовой динамики встречаются при решении задач во многих областях современной науки. Нелинейность уравнений газовой динамики выдвигает численные методы их решения с использованием ЭВМ на первый план по отношению к другим методам. И на сегодняшний день численные методы являются универсальным методом решения этих уравнений.

Математическое' моделирование реальных газодинамических течений сводится к решению многомерных уравнений газовой динамики [1, 2], что в еще большей мере усложняет теоретические исследования. Поэтому в данной области ведущая роль отводится прямому численному моделированию (DNS - Direct Numeric Simulation) газодинамических течений.

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

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

Основными методами численного решения нелинейных уравнений газовой динамики являются разностные методы [3, 4]. Большинство используемых на сегодняшний день разностных схем для решения систем уравнений газовой динамики [5] конструируются так называемым методом конечных объемов (finite volume method). Суть этого метода в том, что область интегрирования разбивается на множество ячеек и интегралы, выражающие законы сохранения в ячейках, заменяются конечными суммами. При этом предполагается, что сеточные значения в ячейках есть интегральные средние соответствующей величины по объему и потоки на границах между ячейками (дискретные потоки) определяются тем или иным способом. От интегро-интерполяционного метода [3] данный метод отличается тем, что здесь для всех величин используется одна и та же разностная сетка.

Наиболее полное описание существующих на сегодняшний день методов можно найти в работах [5, б].

Множество существующих методов вычисления дискретных потоков основываются на идеях метода С. К. Годунова, предложенного в работе [7]. Здесь описываются разностные схемы для решения нестационарных задач, и предложен способ вычисления потоков на гранях между ячейками дискретной сетки как решений задачи о распаде произвольного разрыва. Решения задачи о распаде разрыва могут находятся или из линейных приближений - при дозвуковом течении, или итерационными методами - при сверхзвуковых течениях. Многие методы вычисления дискретных потоков такие, как методы Куранта-Изаксона-Риса [9], Лакса-Фридрихса [10], Ошера [13, 13], Роу [12] и Хартена-Лакса-Ван Лира [15], по сути являются упрощенными вариантами метода С. К. Годунова. В литературе эти методы носят название методов типа Годунова.

Но метод Годунова имеет только первый порядок точности. В работе [7] доказано, что схемы порядка выше первого не являются монотонными, что приводит к появлению нефизических осцилляций вблизи разрывов (в частности на сильных ударных волнах).

Но сложность решаемых задач требует высокой точности получаемых решений. И первой [5] работой, в которой был предложен метод получения схем с порядком точности большим первого является работа [16]. Здесь и в последующих работах описываются способы перехода от разностных схем высокого порядка аппроксимации к монотонным схемам первого порядка около особенностей решения. Среди таких работ можно выделить работы Ван Лира [17] - описываются «монотонизированные» схемы, Бориса и Бука [18, 19, 20] - описывается алгоритм расчета переноса с коррекцией потоков (РСТ-метод).

В работах Хартена [21, 22] был предложен метод, получивший название ТУБ-метод или метод невозрастания полной вариации решения. Вместо монотонности этот метод обеспечивает невозрастание полной вариации и более точно передает характер поведения разрывных решений. Суть данных методов сводится к использованию разнообразных «монотонизирующих» ограничителей потоков (Нт^егз) с переключателями, зависящими от свойств решений [23, 24].

Также широкое распространение получил метод реконструкции численных решений, сводящийся к замене определенной в ячейке сеточной величины кусочно-постоянной функцией, имеющей разрывы на границах между ячейками. Тогда для повышения точности необходимо более точно проинтерполировать сеточные величины так, чтобы сохранить интегральные средние значения и получить распределение этой величины в ячейке с минимальными скачками на границах. А для обеспечения монотонности необходимо обеспечить отсутствие новых максимумов и минимумов при реконструкции и чтобы на границах скачок не менял свой знак. Наиболее популярным среди этих методов является метод кусочно-параболической реконструкции (метод PPM - pieccwise parabolic method).

В работе [25] предложен новый способ реконструкции сеточных значений с автоматическим анализом гладкости решения, получивший название ENO схемы (essentially non-oscillatory scheme). В данном методе вместо фиксированного шаблона для интерполяции значений в ячейках сетки используется несколько шаблонов и выбирается тот, на котором решение является наиболее гладким. Этот метод получил дальнейшее развитие в работах [26, 27].

Далее вместо использования при реконструкции значения, полученного на одном шаблоне, в работе [28] было предложено использовать выпуклую линейную комбинацию значений полученных на всех возможных шаблонах, где весовые коэффициенты в линейной комбинации подбираются в зависимости от гладкости решения на каждом шаблоне, используя так называемые индикаторы гладкости. Данный метод получил название WENO (weighted ENO). В работах [29, 30, 31] были предложены новые способы вычисления весовых коэффициентов, а точнее новые способы анализа гладкости решения на шаблоне.

В работе [32] показано, что вблизи гладкого экстремума в ШЕ1Ч0 схемах понижается порядок точности и предложен новый масштабированный метод ^/УЕМО. Где найденные весовые коэффициенты дополнительно масшабируются с помощью специально выбранных функций, что позволяет добиться высокого порядка точности вблизи экстремума.

Бурное развитие высокопроизводительных вычислительных систем и потребности современной науки в моделировании сложных газодинамических течений способствуют развитию параллельных алгоритмов численного решения уравнений газовой динамики [33, 34]. В настоящее время технологии параллельных вычислений посвящено достаточно много работ [35, 36, 37, 38, 39], также много и электронных ресурсов освещающих- актуальные вопросы данного направления [40, 41, 42].

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

• построение существенно неосциллирующих разностных схем высокого порядка точности для решения многомерной системы уравнений газовой динамики в Эйлеровых переменных;

• тестирование построенных схем на известных одномерных и двумерных модельных задачах;

• разработка параллельного вычислительного алгоритма для построенных схем и создание программного комплекса на их основе;

• проведение прямого численного моделирования для одной задачи о развитии неустойчивости Рихтмайера-Мешкова и сравнение результатов с экспериментальными данными и результатами, полученными с использованием других схем.

Научная новизна работы заключается в следующем:

• в работе предложены новые существенно неосциллирующие схемы высокого порядка точности с масштабированием весовых коэффициентов;

• построен параллельный вычислительный алгоритм для моделирования газодинамических течений;

• разработан программный комплекс для прямого численного моделирования газодинамических течений на высокопроизводительных вычислительных системах;

• проведено прямое численное моделирование развитой турбулентной стадии неустойчивости Рихтмайера-Мешкова с использованием разработанного программного комплекса на кластере, составленном из персональных ЭВМ;

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

Достоверность научных выводов и результатов подтверждается следующим:

• построенные вычислительные алгоритмы апробированы на тестовых задачах и показывают работоспособность для различных типов задач;

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

• используемые математические модели базируются на фундаментальных законах сохранения массы, импульса и энергии.

Автор данной работы выносит на защиту:

• построенную существенно неосциллирующую разностную схему высокого порядка точности для решения многомерной системы уравнений газовой динамики в переменных Эйлера;

• параллельный вычислительный алгоритм для построенных схем;

• разработанный программный комплекс для параллельных вычислений;

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

По главам содержание работы распределено следующим образом. В первой главе описан класс существенно неосциллирующих разностных схем высокого порядка точности для численного решения Эйлеровой системы уравнений газовой динамики. Представлены результаты численного решения с помощью построенных схем одномерной задачи о распаде разрыва с начальными данными Сода [49] и с начальными данными Лакса [50].

Вторая глава включает в себя описание параллельного вычислительного алгоритма для построенных разностных схем. Описан разработанный программный комплекс ОВ"У^Е1ЮРА11ЗВ для вычислений па многопроцессорных системах, реализующий построенные алгоритмы.

Проведен анализ эффективности распараллеливания для кластера, построенного на базе персональных ЭВМ.

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

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

• семинарах Средневолжского математического общества под руководством профессора Е. В. Воскресенского;

• XXXIV, XXXV, XXXVI Огаревских чтениях в Мордовском государственном университете им. Н. П. Огарева (г. Саранск, 2005, 2006, 2007 гг.);

• XIV, XV, XVI научных конференциях молодых ученых в Мордовском государственном университете им. Н. П. Огарева (г. Саранск, 2005, 2006, 2007, 2008 гг.);

• Второй Международной научной школе «Математическое моделирование, численные методы и комплексы программ» (г. Саранск, 2005 г.);

• VII Международной конференции «Дифференциальные уравнения и их приложения» (г. Саранск, 2006 г.);

• на семинаре отдела №4 Института математического моделирования РАН (г. Москва, 2007 г.);

• Третьей Международной научной школе «Математическое моделирование, численные методы и комплексы программ» (г. Саранск, 2007 г.);

• Восьмой Международной конференции «Диффернциальные уравнения и их приложения» (г. Саранск, 2008 г.).

Основные результаты опубликованы в работах [52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63].

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Жалнин, Руслан Викторович

Заключение

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

1. Построена класс квазимонотонных (существенно неосциллирующих) консервативных разностных схем высокого порядка точности (5-го порядка) с реконструкцией значений на гранях между ячейками, представляющей из себя модификацию хорошо зарекомендовавшей себя известной WENO схемы. Модификация заключается в особом способе вычисления весовых коэффициентов при построении интерполяционного полинома на шаблоне ячеек при реконструкции значений на гранях между ячейками, что позволяет избежать понижения порядка точности вблизи экстремальных точек. Построенная схема протестирована на модельной одномерной задаче Римана о распаде разрыва с начальными данными Сода и Лакса. Результаты расчетов хорошо согласуются с известными решениями данных задач.

2. Разработан параллельный программный комплекс GDWEN0PAR3D для прямого численного моделирования газодинамических течений с использованием построенных разностных схем. Построен параллельный вычислительный алгоритм на основе технологии MPI, для выполнения расчетов по построенным разностным схемам, оптимизированный под вычисления на кластере, составленном из объединенных сетью Ethernet персональных ЭВМ. Рассмотрены различные способы оценки эффективности параллельного алгоритма. Проведен анализ эффективности алгоритма при различных способах декомпозиции расчетной области. Выполнен анализ объемов межпроцессорных обменов при различных способах декомпозиции.

3. Выполнено прямое численное моделирование развития неустойчивости Рихтмайера-Мешкова с использованием программного комплекса СО\¥ЕГЮРА11ЗВ. Выполнено сравнение результатов расчетов с экспериментальными данными на основе значений пульсаций продольной компоненты скорости. Показано, что предлагаемый вычислительный алгоритм дает более приближенные к экспериментальным результаты по сравнению с вычислительными экспериментами других авторов.

Список литературы диссертационного исследования кандидат физико-математических наук Жалнин, Руслан Викторович, 2008 год

1. Ландау Л. Д., Лифшиц Е.М., Механика сплошных сред, М., Гостехиздат, 1954.

2. Седов Л. И., Механика сплошной среды, в 2 т. Т. 1-2, СПб.: Издательство «Лань», 2004

3. А.А.Самарский, Теория разностных схем, М., Наука, 1977.

4. A.A.Самарский, Ю.П.Попов, Разностные методы решения задач газовой динамики, М., Наука, 1980.

5. Бондаренко Ю.А., Башуров В.В., Янилкин Ю.В. Математические модели и численные методы для решения задач нестационарной газовой динамики. Обзор зарубежной литературы. Препринт №88, Саров: РФЯЦ-ВНИИЭФ, 2003, 53 с.

6. Куликовский А.Г., Погорелов Н.В., Семенов А.Ю., Математические вопросы численного решения гиперболических систем уравнений, М.: ФИЗМАТЛИТ, 2001, 608 с.

7. Годунов С. К., Разностный метод численного расчета разрывных решений уравнений гидродинамики // Матем. сб., 1959, 47, вып. 3, 271-306.

8. Годунов С.К., Забродин М.Я., Иванов М.Я., Крайко А.Н., Прокопов Г.П. Численное решение многомерных задач газовой динамики, М., Наука, 1976

9. Courant R., Isaacson E., Rees M. On the solution of nonlinear hyperbolic differential equations by finite differences, Comm. Pure Appl. Math. 5, No. 3, 1952, p. 243-255

10. Weak solutions of nonlinear hyperbolic equations and their numerical computation, Comm. Pure Appl. Math. 7, No. 1, 1954, p. 159-193

11. Холодов А. С. О построении разностных схем с положительной аппроксимацией для уравнений гиперболического тина, Ж. выч. матем. и матем. физики 18, № б, 1978, с. 1476 1492

12. Roe P.L. Approximate Rieman problem solvers, parameter vectors, and diference schemes, Journal of computational physics 43, No. 2, 1981, p. 357 372

13. Osher S. Numerical solution of singular perturbation problems and hyperbolic systems of conservation laws, North Holland Mathematical Studies 47, 1981, p. 179 205

14. Osher S. Rieman solvers, the entropy conditions, and differenceapproxima-tions, SIAM J. Numer. Anal. 21, No. 2, 1984, p. 217 235

15. Harten A., Lax P.D., van Leer B. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Review 25, No. 1, p. 35 61

16. Колган В.П. Применение принципа минимальных значений производной к построению конечно-разностных схем для расчета разрывных решений газовой динамики, Ученые записки ЦАГИ, 1972, т. 3, № 6, с. 68 77.

17. Van Leer В. Flux-vector splitting for the Euler equations, Lecture notes in physics, 1982, v. 170, p. 507 512

18. Boris J.P., Book D.L. Flux-corrected transport, I. SHASTA, a fluid transport algorith that works, Journal of computational physics, 1973, v. 11, No. 1, p. 38 69

19. Boris J.P., Book D.L., Hain K. Flux-corrected transport, II. Generalizations of the method, Journal of computational physics, 1975, v. 18, No. 3, p. 248 283

20. Boris J.P., Book D.L. Flux-corrected transport, III. Minimal-error FCT algorithms, Journal of computational physics, 1976, v. 20, No. 4, p. 397 -431

21. Harten A. High resolution schemes for hyperbolic conservation laws, Journal of computational physics, 1983, v. 49, p. 357 393 1

22. Harten A. On a class of high resolution total-variation-stabe finite-difference schemes, SIAM J. Numer. Anal., 1984, v. 21, p. 1 23

23. Вязников К. В., Тишкин В. Ф., Фаворский А. П., Построение монотонных разностных схем повышенного порядка аппроксимации для систем уравнений гиперболического типа, Математическое моделирование, 1989, т. 1, № 5, с. 95 120.

24. Тишкин В. Ф., Никишин В. В., Попов И. В., Фаворский А. П., Разностные схемы трехмерной газовой динамики для задачи о развитии неустойчивости Рихтмаера-Мешкова, Математическое моделирование, 1995, т. 7, № 5, с. 15 25.

25. A. Harten, В. Engquist, S. Osher and S. Chakravarthy, Uniformly high order essentialy non-oscillatory schemes. Ill, Journal of computational physics, 71, p. 231 303 (1987).

26. Shu C.-W., Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of computational physics, 1988, v. 77, No. 2, p. 439 471

27. Shu C.-W., Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes. II, Journal of computational physics, 1989, v. 83, No. 1, p. 32 78

28. Liu X.-D., Osher S., Chan T. Weighted essentially non-oscillatory schemes, Journal of computational physics, 1994, v. 115, p. 200

29. Jiang G.-S., Shu C.-W. Efficient implementation of weighted ENO schemes, Journal of computational physics, 126 (1996). P. 202-228.

30. C.-W. Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, ICASE Report 97-65, (1997).

31. Z. J. Wang, R. F. Chen, Optimized weighted essentially non-oscillatory schemes for computational aeroacoustics, Published by AIAA Inc. with permission, (2001).

32. Henrick A. K., Aslam T. D., Powers J. M. Mapped weighted essentially non-oscillatory schemes: Archieving optimal order near critical points, Journal of Computational Physics, 207, 2005, p. 542 567.

33. Воеводин В. В., Воеводин Вл. В., Параллельные вычисления, С-Пб., БХВ-Петербург, 2002

34. Антонов А. С., Параллельное программирование с использованием технологии MPI, Изд-во МГУ, Москва, 2004

35. Якобовский М.В., Распределенные системы и сети, Изд-во Станкин, 2000.

36. Гергель В. П., Теория и практика параллельных вычислений: учебное пособие, М., БИНОМ. Лаборатория знаний, 2007, 423 с.40. www.parallel.ru41. www.mpi-forum.org42. www.cfd-online.com

37. Змитренко H. В., Ладонкина M. Е., Тишкин В. Ф., Численное исследование турбулентного перемешивания для одной задачи о развитии неустойчивости Рихтмаера-Мешкова, ВАНТ, сер. Мат. моделир. физ. процессов, 2004, вып. 1.

38. Poggi F., Thorembey М.-Н., Rodrigues G. Velocity measurements in turbulent gaseous mixtures induced by Richtmyer-Meslikov instability // Phisics of Fluids. 1998. Vol. 10, No 11. P. 2698-2700.

39. Левитан Ю. Л., Соболь И. M. О датчике псевдослучайных чисел для персональных компьютеров // Математическое моделирование. 1990. Т. 2, № 8. С. 119-126.

40. Ладонкина M. Е., Численное моделирование турбулентного перемешивания с использованием высокопроизводительных систем, Диссертация на соискание ученой степени канд. физ.-мат. наук, ИММ РАН, Москва, 200447. www.dislin.de

41. Woodward P., Colella P. The numerical simulation of two-dimensional fluid flow with strong shocks // Journal of Computational Physics. 1984. V.54, PP.115-173.

42. Sod A.G. Review. A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, J. Comput. Phys. 27, No. 1, p. 1 31

43. Lax P.D., Weak solutions of nonlinear hyperbolic equations and their numerical computation, Comm. Pure Appl. Math. 1954, No. 7, p. 159 — 193

44. Double Mach Reflection of a Strong Shock (http://www.math.ntnu.no/ andreas /fronttrack / gas / wedge / index, html)

45. Жалнин P.B. Построение неосциллирующих алгоритмов (ENO-схем) для гиперболических уравнений, (http: //svmo.mrsu.ru/lib / cmu05 / zhalnin.pdf)

46. Жалнин Р. В. Построение ENO-схем для одномерного случая, Труды СВМО, 2005, т. 7, № 1, сс. 407 408

47. Жалнин Р. В., Численная реализация неосциллирующих вычислительных алгоритмов, Материалы научной конференции XXXIV Огаревкие чтения, Саранск: СВМО, 2005, сс. 35 38

48. Жалнин Р.В. Сравнение методов минимизации вариации интерполяционного полинома для WENO схем, Материалы XI научной конференции Мордовского государственного университета молодых ученых, аспирантов и студентов, Саранск: СВМО, 2006, сс. 16 18

49. Тишкин В. Ф., Жалнин Р. В. Одномерные неосциллирующие схемы на равномерной сетке с минимизацией вариации интерполяционного полинома в ячейках, содержащих точку локального экстремума, Труды СВМО, 2006, т. 8, № 1, сс. 115 121

50. Жалнин Р. В., О реализации параллельных вычислительных алгоритмов на персональных ЭВМ объединенных локальной сетью, Материалы научной конференции XXXV Огаревкие чтения, Саранск: СВМО, 2006, сс. 20 25

51. Жалнин Р.В. О прямом численном моделировании развития неустойчивости Рихтмайера-Мешкова, Материалы XII научной конференции Мордовского государственного университета молодых ученых, аспирантов и студентов, Саранск: СВМО, 2007, сс. 17 19

52. Жалнин Р. В. О построении параллельного вычислительного алгоритма высокого порядка точности для гиперболических систем уравнений, Труды СВМО, 2007, т. 9, № 1, сс. 145 153

53. Жалнин Р. В., К вопросу о численном моделировании развития неустойчивости Рихтмайера-Мешкова, Материалы научной конференции XXXVI Огаревкие чтения, Саранск: СВМО, 2007, сс. 18 20

54. Жалнин Р.В., Змитренко Н.В., Ладонкина М.Е., Тпшкин В.Ф., Численное моделирование развития неустойчивости Рихтмайера-Мешкова с использованием схем высокого порядка точности // Мат. моделирование, 2007, том 19, №10, сс. 61-66

55. Жалнин Р. В. О построении параллельного вычислительного алгоритма высокого порядка точности для гиперболических систем уравнений, Труды СВМО, 2008, т. 10, № 1, 137 146