CC BY
56
www.volsu.ru
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ
DOI: http://dx.doi.org/10.15688/jvolsu1.2015.1.4
УДК 519.6 ББК 22.193
РАЗРАБОТКА ПРОГРАММЫ ДЛЯ ЧИСЛЕННОГО ГАЗОДИНАМИЧЕСКОГО МОДЕЛИРОВАНИЯ НА ОСНОВЕ
____о____о ____л
ЛАГРАНЖЕВО-ЭИЛЕРОВОИ СХЕМЫ LES - ASG1
Белоусов Антон Владимирович
Студент Института математики и информационных технологий, Волгоградский государственный университет anton.belousov@mail.ru, infomod@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Храпов Сергей Сергеевич
Кандидат физико-математических наук, доцент кафедры информационных систем и компьютерного моделирования, Волгоградский государственный университет xss-ip@mail.ru, infomod@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Аннотация. Представлены результаты численного моделирования уравнения переноса с неоднородным распределением скорости для одномерного случая. Проведен анализ и сравнение численных схем LES и MUSCL для кусочно-линейной и кусочно-постоянной реконструкции.
ю
3 Ключевые слова: численные схемы, LES, ASG, MUSCL, лагранжево-
эйлерова схема.
и и
m
g Введение
a X
шг В данной работе проведено сравнение двух типов решения уравнения переноса с
< неоднородным распределением скорости, а именно LES (лагранжево-эйлерова схема) и
о MUSCL (Monotonic Upstream-Centered Scheme). По результатам сравнения можно опре-
о делить эффективность схем, а также точность по отношению к вычислительным затра-
рэ там для схем равного порядка точности. Все полученные данные будут способствовать @ дальнейшей реализации схемы LES — ASG уравнений газодинамики.
1. Численные схемы
Разработка новых и модификация уже существующих численных методов для повышения их эффективности всегда были актуальными задачами вычислительной математики. Это связано как с практической необходимостью нахождения решений новых сложных задач, так и самой логикой развития численных методов как раздела прикладной математики.
Для построения численного алгоритма необходимо прежде всего четко поставить математическую задачу, эквивалентную данной физической задаче. Кроме того, алгоритм, позволяющий корректно производить вычисления, должен удовлетворять ряду требований, чего не всегда просто добиться в каждой конкретной ситуации.
Известно, что решения различных задач математической физики, которые описываются гиперболическими системами уравнений, могут быть гладкими в одних областях и разрывными в других [1; 5; 6; 8; 9]. Разрывные решения могут при этом возникать из гладких начальных данных. Такие свойства решений налагают на алгоритмы численного решения достаточно жесткие требования. С одной стороны, численный метод должен уметь сохранять свойство монотонности в тех областях, где искомые решения имеют большие перепады значений. С другой стороны, тот же метод должен обладать высоким порядком точности в тех областях, где решение является гладким.
1.1. Схема MUSCL
Данная схема относится к методу конечных объемов, обеспечивает высокую точность решения. Идея состоит в том, чтобы заменить кусочно-постоянную аппроксимацию на реконструкцию ячеек и получить их промежуточные состояния на каждом временном шаге. Для каждой ячейки ограничен угол наклона, реконструкция левой и правой ячейки производится непосредственно через вычисление потоков на границах ячеек [3; 11].
1.2. Схема LES
Численные схемы, основанные на лагранжевом или эйлеровом подходах, обладают как преимуществами, так и недостатками. Комбинирование лагранжева и эйлерова подходов в рамках одного алгоритма позволяет объединить преимущества обоих методов. Для численного решения можно применить схему «предиктор — корректор». LES является обобщением новой численной схемы cSPH — TVD (combined Smoothed Particle Hydrodynamics — Total Variation Diminishing) [4; 7; 10; 12; 13]. Схема LES строится для уравнения переноса, не только для гидродинамического моделирования, но и прекрасно подходит для задач по переносу загрязненных примесей в приземном слое атмосферы [2].
1.3. Кусочно-постоянная реконструкция
Рассмотрим одномерное уравнение переноса с неоднородным распределением скорости:
f + ^ = 0, V = V (t) = const. (!)
Применив схему MUSCL к уравнению (1), получим
^п+1
где
Q-+1 = (1 -Ai(Ï))Q! + MÏ)QU, (2)
Ai(0 = c^i/a (1 - Cj+1/2 -Cj-1/2) , A2{Ï) = Cj—1/2 (l - C-1/2 - C-3/2) .
Следует учитывать условия устойчивости
Для реализации схемы LES уравнение (1) будет иметь вид
В2(г ¿+1
где
— Сг+1/2, В2(г) — Сг-1/2.
Погрешность для обеих схем рассчитываем через соотношение (7)
1.4. Кусочно-линейная реконструкция
Применив схему MUSCL к уравнению (1), получим
l2 + A4(W-1 -2 где В1,В2,В3 и В4 рассчитываем через соотношение (9)
Cj+1 /2 — Cj_ 1
A2(i) = Cj—1/2 (l - Cj-1/2 - Cj-3/2 )
Аз(г) = Cj+1/2 (l - Cj+1/2 + Cj-1/2) A4 (г) = Cj-1/2 (l - Cj-1/2 + Cj-3/2 )
(3)
ViT
Cj = V-, Ci < 1. (4)
h
Q7+1 =(1 - + В+%1, (5)
(6)
1 N
£ = ^E IQj-Qjl, Cj < 1. (7)
j=1
Q7+1 = (1 - A1(i))Q? + A2(Î)Q?_ 1 - А3(г)А?- + А4(г)Л™_^, (8)
(9)
Как и в кусочно-постоянной реконструкции, следует учитывать условия устойчивости (4).
Для реализации схемы LES применим шаг «предиктор — корректор» к уравнению (1) и в результате получим
= (1 - ^pjQ1! + -ЗДА?^ + f2, (10)
где B1,B2,B3,B4, 5+ и 5— рассчитываем через соотношения (11) и (12)
В lit) = Ci+1/2, B2(i) = Ci-1/2, Вз(г) = Ci+i/2(1 -Ci), B4(i) = Q-1/2(1 - C*_i).
(11)
= 1 + C
г+1/2
-Ci.
= 1 + Ci-1/2 — Ci-1.
(12)
2. Проведение вычислительных экспериментов
Проведем сравнение численных схем MUSCL и LES для кусочно-линейной и кусочно-постоянной реконструкции для решения уравнения переноса с неоднородным распределением скорости. Следует обратить внимание, что обе схемы совпадают при U = = const, а также схемы хороши для реализации переноса гладких профилей, поэтому рассмотрим, как будет вести себя профиль, представленный на рисунках 1 и 2.
Рис. 1. Исходный профиль Рис. 2. Профиль скорости
Проведем тестовые расчеты и получим решения на сетках с частотой дискретной пространственной дискретизации N: 300, 600, 1 200, 2 400, 4 800, 9 600. Будем использовать ограничитель minmod. Количество итераций для N = 300 будет 970, а для каждого последующего N количество итераций увеличивается в два раза, для удовлетворения условия (4).
A? = L
иг+1 -иг иг -U-
h
h
L(a, b) = — (sign(a) + sign(b))min(la\, |6|).
(13)
Рис. 3. Конечные профили для схемы LES при N = 300, где: сплошная линия — это кусочно-постоянная реконструкция; прерывистая линия — это кусочно-линейная реконструкция; пунктирная линия — это исходный профиль
Рис. 4. Конечные профили для схемы МШСЬ при N = 300, где: сплошная линия — это кусочно-постоянная реконструкция; прерывистая линия — это кусочно-линейная реконструкция; пунктирная линия — это исходный профиль
Обратим внимание на то, что обе схемы (рис. 3 и 4) почти полностью совпадают в конечных профилях. Более подробные результаты испытаний представлены на рисунках 5-8, а также в таблице 1.
Таблица 1
Результаты тестов, погрешность вычислений
N Кусочно-постоянная реконструкция Кусочно-линейная реконструкция
MUSCL LES MUSCL LES
300 7,5455 х 10-2 7,5451 х 10-2 2,0163 х 10-2 1,0191 х 10-2
600 4,2171 х 10-2 4,2169 х 10-2 5,7162 х 10-3 5, 7275 х 10-3
1200 2,2633 х 10-2 2,2633 х 10-2 2,6335 х 10-3 2, 6309 х 10-3
2 400 1,1818 х 10-2 1,1818 х 10-2 1,2402 х 10-3 1,2411 х 10-3
4 800 6, 0794 х 10-3 6, 0794 х 10-3 5, 9244 х 10-4 5,9223 х 10-4
9600 3,0903 х 10-3 3, 0903 х 10-3 2,7158 х 10-4 2, 7163 х 10-4
Составим формулу, по которой будем сравнивать кусочно-линейную схему MUSCL с LES.
N
F(x) = log1o -Q**, (14)
г=о
где Qг — решение по MUSCL, а Q* — решение по LES.
Рис. 5. Вырезка из конечного профиля схемы LES с кусочно-линейной реконструкцией на сетках N = 300 (сплошная линия), N = 2 400 (прерывистая линия), N = 9 600 (пунктирная
линия), исходный профиль (плюс)
Рис. 6. Профиль схемы LES с кусочно-линейной реконструкцией в моменты итераций 150 (квадрат), 300 (плюс), 450 (круг), 700 (крест), 850 (ромб), для N = 300
При проведении тестовых расчетов мы получили данные, на основе которых можем полагать, что численные схемы MUSCL и LES одинаково хорошо применимы для моделирования уравнений переноса.
По итогам проведенной работы мы удостоверились в правильности выбора схемы LES для последующей разработки программы для численного газодинамического моделирования на основе схемы LES — ASG.
Рис. 7. Отличие кусочно-линейной схемы MUSCL от LES по формуле (14)
Рис. 8. Погрешность вычислений кусочно-линейной схемы ЬББ
ПРИМЕЧАНИЕ
1 Работа выполнена при поддержке гранта РФФИ № 15-02-06204.
СПИСОК ЛИТЕРАТУРЫ
1. Алиев, А. В. Иерархический БРИ-метод для математического моделирования в гравитационной газовой динамике / А. В. Алиев, Г. А. Тарнавский // Сиб. электрон. мат. изв. - 2007. - С. 376-434.
2. Белослудцев, А. А. Информационно-компьютерный комплекс для моделирования динамики примесей от предприятий химической промышленности / А. А. Белослудцев, Д. В. Гусаров, М. А. Еремин, Н. М. Кузьмин, С. А. Хоперсков, С. С. Храпов // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2009. — № 12. — С. 95-102.
3. Еремин, М. А. Конечно-объемная схема интегрирования уравнений гидродинамики / М. А. Еремин, А. В. Хоперсков, С. А. Хоперсков // Изв. Волгогр. гос. техн. ун-та. — 2010. — № 6:8. — С. 24-27.
4. Жумалиев, А. Г. Численная схема сБРИ — ТУБ: моделирование фронта ударной волны / А. Г. Жумалиев, С. С. Храпов // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2012. — № 16. — С. 60-67.
5. Кузьмин, Н. М. Численное моделирование эволюции неустойчивых мод джетов, выходящих из молодых звездных объектов / Н. М. Кузьмин, В. В. Мусцевой, С. С. Храпов // Астрон. журн. — 2007. — № 84:12. — С. 1089-1098.
6. Куликовский, А. Г. Математические вопросы численного решения гиперболических систем уравнений / А. Г. Куликовский, Н. В. Погорелов, А. Ю. Семенов. — М. : ФИЗМАТ-ЛИТ, 2001. — 608 с.
7. Кузьмин, Н. М. Численная схема сБРИ — ТУБ: исследование влияния ограничителей наклонов / Н. М. Кузьмин, А. В. Белоусов, Т. С. Шушкевич, С. С. Храпов // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2014. — № 1. — С. 22-33.
8. Писарев, А. В. Особенности динамики затопления Волго-Ахтубинской поймы в зависимости от режимов испарения и инфильтрации / А. В. Писарев, С. С. Храпов,
А. А. Воронин, Т. А. Дьяконова, Е. А. Циркова // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2012. — № 1. — C. 36-41.
9. Писарев, А. В. Численная модель динамики поверхностных вод в русле Волги: оценка коэффициента шероховатости / А. В. Писарев, С. С. Храпов, Е. О. Агафонникова, А. В. Хоперсков // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. — 2013. — № 1. — C. 114-130.
10. Храпов, C. C. Численная схема для моделирования динамики поверхностных вод на основе комбинированного SPH-TVD-подхода / C. C. Храпов, А. В. Хоперсков, Н. М. Кузьмин, А. В. Писарев, И. А. Кобелев // Вычислительные методы и программирование. —
2011. — Т. 12, № 1. — C. 282-297.
11. Хоперсков, А. В. Динамика газового диска в неосесимметричном темном гало / А. В. Хоперсков, М. А. Еремин, С. А. Хоперсков, М. А. Бутенко, А. Г. Морозов // Астрон. журн. — 2012. — № 89:1. — C. 19-31.
12. Хоперсков, А. В. Задача управления гидрологическим режимом в эколого-экономической системе «Волжская ГЭС — Волго-Ахтубинская пойма». Ч. 1. Моделирование динамики поверхностных вод в период весеннего паводка / А. В. Хоперсков, С. С. Храпов, А. В. Писарев, А. А. Воронин, М. В. Елисеева, И. А. Кобелев // Проблемы управления. —
2012. — № 5. — C. 18-25.
13. Khrapov S. The Numerical Simulation of Shallow Water: Estimation of the Roughness Coefficient on the Flood Stage / S. Khrapov, A. Pisarev, I. Kobelev, A. Zhumaliev, E. Agafonnikova, A. Losev, A. Khoperskov // Advances in Mechanical Engineering. — 2013. — Vol. 2013. — Article ID 787016, 11 pages.
REFERENCES
1. A^v A.V., Tarnavskiy G.A. Iеrarkhichеskiy SPH-mеtod dlya matеmatichеskogo modеlirovaniya v gravitatsionnoy gazovoy dinamikе [The hierarchical SPH-method for mathematical modeling in gravitational gas dynamics]. Sib. e^ktron. mat. izv., 2007, pp. 376434.
2. Bеlosludtsеv A.A., Gusarov D.V., Eremin M.A., Kuzmin N.M., Khopеrskov S.A., Khrapov S.S. Informatsionno-kompyutеrnyy kompkks dlya modеlirovaniya dinamiki primеsеy ot predpriyatiy khimichеskoy promyshknnosti [Information and computer system for modeling the dynamics of impurities from the chemical industry]. Vеstnik Volgogradskogo gosudars^nnogo univеrsitеta. Sеriya 1, Matеmatika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2009, no. 12, pp. 95-102.
3. Eremin M.A., Khopеrskov A.V., Khopеrskov S.A. Konеchno-obyеmnaya skhеma intеgrirovaniya uravnеniy gidrodinamiki [A finite-volume scheme for integrating the equations of hydrodynamics]. Izv. Volgogr. gos. tеkhn. un-ta, 2010, no. 6:8, pp. 24-27.
4. Zhuma^v A.G., Khrapov S.S. Chistennaya skhеma cSPH — TVD: modеlirovaniе fronta udarnoy volny [Numerical scheme cSPH-TVD: simulation of the shock front]. Vеstnik Volgogradskogo gosudars^nnogo univеrsitеta. Sеriya 1, Matеmatika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2012, no. 16, pp. 60-67.
5. Kuzmin N.M., Mustsеvoy V.V., Khrapov S.S. Chislеnnoе modеlirovaniе evolyutsii nеustoychivykh mod dzhеtov, vykhodyashchikh iz molodykh zvеzdnykh obyеktov [Numerical Modeling of the Evolution of Unstable Modes of Jets from Young Stellar Objects]. Astron. zhurn. [Astronomy Reports], 2007, no. 84:12, pp. 1089-1098.
6. Kulikovskiy A.G., Pogorelov N.V., Sеmеnov A.Yu. Matеmatichеskiе voprosy chishnnogo rеshеniya gipеrbolichеskikh sistеm uravmniy [Mathematical problems in the numerical solution of hyperbolic systems]. Moscow, FIZMATLIT Publ., 2001. 608 p.
7. Kuzmin N.M., Bеlousov A.V., Shushkеvich T.S., Khrapov S.S. Chistennaya skhеma sSPH — TVD: isslеdovaniе vliyaniya ogrem^^^ naklonov [Numerical scheme CSPH - TVD: study of the impact limiters slopes]. Vеstnik Volgogradskogo gosudarstvеnnogo univеrsitеta. Sеriya 1, Matеmatika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2014, no. 1, pp. 22-33.
8. Pisarev A.V., Khrapov S.S., Voronin A.A., Dyakonova T.A., Tsirkova E.A. Osobennosti dinamiki zatopleniya Volgo-Akhtubinskoy poymy v zavisimosti ot rezhimov ispareniya i infiltratsii [Features of the dynamics of flooding of the Volga-Akhtuba floodplain, depending on the mode of evaporation and infiltration]. Vеstnik Volgogradskogo gosudars^nnogo univеrsitеta. Sеriya 1, Matеmatika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2012, no. 1, pp. 36-41.
9. Pisarev A.V., Khrapov S.S., Agafonnikova E.O., Khoperskov A.V. Chislennaya model dinamiki poverkhnostnykh vod v rusle Volgi: otsenka koeffitsienta sherokhovatosti [Numerical model of dynamics of surface waters in line with the Volga: the estimated coefficient of roughness]. Vеstnik Udmurtskogo univеrsitеta. Matеmatika. Mеkhanika. Kompyutеrnyе nauki, 2013, no. 1, pp. 114-130.
10. Khrapov C.C., Khoperskov A.V., Kuzmin N.M., Pisarev A.V., Kobelev I.A. Chislennaya skhema dlya modelirovaniya dinamiki poverkhnostnykh vod na osnove kombinirovannogo SPH-TVD-podkhoda [Numerical scheme for modeling the dynamics of surface water based on the combined SPH-TVD-approach]. Vychislitеlnyе mеtody i programmirovaniе, 2011, vol. 12, no. 1, pp. 282-297.
11. Khoperskov A.V., Eremin M.A., Khoperskov S.A., Butenko M.A., Morozov A.G. Dinamika gazovogo diska v neosesimmetrichnom temnom galo [Dynamics of Gaseous Disks in a Non-axisymmetric Dark Halo]. Astron. zhurn. [Astronomy Reports], 2012, no. 89:1, pp. 19-31.
12. Khoperskov A.V., Khrapov S.S., Pisarev A.V., Voronin A.A., Eliseeva M.V., Kobelev I.A. Zadacha upravleniya gidrologicheskim rezhimom v ekologo-ekonomicheskoy sisteme «Volzhskaya GES — Volgo-Akhtubinskaya poyma». Ch. 1. Modelirovanie dinamiki poverkhnostnykh vod v period vesennego pavodka [The task of managing the hydrological regime in the ecological and economic system «Volga Hydroelectric Power Station — Volga-Akhtuba floodplain». Part 1. Modeling the dynamics of surface water during spring flood]. Probhmy uprav^niya [Control Sciences], 2012, no. 5, pp. 18-25.
13. Khrapov S., Pisarev A., Kobelev I., Zhumaliev A., Agafonnikova E., Losev A., Khoperskov A. The Numerical Simulation of Shallow Water: Estimation of the Roughness Coefficient on the Flood Stage. Advances in Mechanical Engineering, 2013, vol. 2013. Article ID 787016, 11 p.
DEVELOPMENT OF THE PROGRAM FOR NUMERICAL GASDYNAMIC
MODELING ON THE BASIS OF LAGRANGE - EULER SCHEME THE LES - ASG
Belousov Anton Vladimirovich
Student, Institute of Mathematics and Information Technologies, Volgograd State University anton.belousov@mail.ru, infomod@volsu.ru Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Khrapov Sergey Sergeevich
Candidate of Physical and Mathematical Sciences,
Associate Professor, Department of Information Systems and Computer Simulation, Volgograd State University xss-ip@mail.ru, infomod@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Abstract. In this paper, the authors compare the two types of solving the transport equation with inhomogeneous distribution of velocity, namely LES (Lagrange — Euler scheme) and MUSCL (Monotonic Upstream-Centered Scheme)
Computational experiments LES schemes and MUSCL, using a piecewise linear and piecewise constant reconstruction, for solving the transport equation with inhomogeneous distribution of velocity. Note that both circuits coincide with U = const, as well as circuitry to implement good smooth transfer profiles. In the calculations limiter minmod is used. According to the results of comparison circuits to determine the effectiveness, as well as accuracy with respect to the computational cost for the accuracy of the order of equal circuits.
The authors study the results of calculations of numerical schemes for LES and MUSCL nets dimension 300, 600, 1 200, 2 400, 4800 and 9600. According to the research we can assume that the scheme LES and MUSCL equally well applicable for modeling the transport equation.
According to the results of the work we made sure that LES scheme and its subsequent development will continue on the basis of the scheme LES-ASG.
Key words: numerical schemes, LES, ASG, MUSCL, Lagrange — Euler scheme.
