СРАВНЕНИЕ ЭФФЕКТИВНОСТИ СПЕКТРАЛЬНОГО И КОНЕЧНО-РАЗНОСТНОГО МЕТОДОВ ПРИ РЕШЕНИИ НЕЛИНЕЙНЫХ ГАЗОВЫХ УРАВНЕНИЙ ТИПА БЮРГЕРСА

ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 44 СРАВНЕНИЕ ЭФФЕКТИВНОСТИ СПЕКТРАЛЬНОГО И КОНЕЧНО-РАЗНОСТНОГО МЕТОДОВ ПРИ РЕШЕНИИ НЕЛИНЕЙНЫХ ГАЗОВЫХ УРАВНЕНИЙ ТИПА БЮРГЕРСА Урокова Шахзода Тоджимуродовна Учитель школы №20 города Термеза [email protected] Аннотация В данной статье рассматривается сравнительный анализ спектрального и конечноразностного методов решения нелинейных газовых уравнений типа Бюргерса. Уравнение Бюргерса является одной из базовых моделей в теории газовой динамики и гидродинамики вязких сред, отражающей взаимодействие нелинейных конвективных и вязкостных эффектов. Сравнение проводится по критериям точности, устойчивости, вычислительной эффективности и способности воспроизводить основные физические характеристики потока, включая формирование ударных волн и диссипацию энергии. Применение спектрального метода на основе преобразования Фурье позволяет достичь экспоненциальной сходимости и высокой точности при моделировании гладких решений, тогда как конечно-разностные схемы удобны для локальных аппроксимаций и не требуют глобальных базисных функций. В работе представлены результаты численных экспериментов, демонстрирующих преимущества и ограничения каждого метода при решении задач газовой динамики. Полученные результаты могут быть использованы для выбора оптимальных вычислительных алгоритмов при численном моделировании течений газа в инженерных и физических приложениях. Ключевые слова. Уравнение Бюргерса, спектральный метод, конечно-разностный метод, газовые потоки, нелинейная динамика, турбулентность, численное моделирование, устойчивость, точность. Abstract This article presents a comparative analysis of spectral and finite-difference methods for solving nonlinear Burgers-type gas equations. The Burgers equation is a fundamental model in the theory of gas dynamics and viscous fluid dynamics, reflecting the interaction of nonlinear convective and viscous effects. The comparison is based on accuracy, stability, computational efficiency, and the ability to reproduce the fundamental physical characteristics of the flow, including shock wave formation and energy dissipation. The use of the spectral method based on the Fourier transform allows for exponential convergence and high accuracy when modeling smooth solutions, while finitedifference schemes are convenient for local approximations and do not require global basis functions. The paper presents the results of numerical experiments demonstrating the advantages and limitations of each method for solving gas dynamics problems. The obtained results can be used to select optimal computational algorithms for numerical modeling of gas flows in engineering and physical applications. Key words. Burgers equation, spectral method, finite difference method, gas flows, nonlinear dynamics, turbulence, numerical modeling, stability, accuracy. ВВЕДЕНИЕ ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 45 Нелинейные процессы в газовой динамике представляют собой одно из наиболее сложных направлений вычислительной физики. При описании газовых потоков, особенно в областях с большими градиентами скорости и давлением, возникают математические трудности, связанные с решением нелинейных дифференциальных уравнений в частных производных. Одним из наиболее простых, но показательных уравнений, отражающих фундаментальные свойства таких систем, является уравнение Бюргерса. Уравнение Бюргерса можно рассматривать как одномерный аналог уравнений Навье — Стокса. Оно описывает совместное действие нелинейной адвекции и вязкостной диффузии и записывается в виде: где — скорость газа, ν — коэффициент кинематической вязкости. В пределе ν→0 уравнение становится гиперболическим и допускает образование ударных волн, что делает его удобным тестом для оценки численных методов, предназначенных для задач газовой динамики. Традиционно для численного решения уравнения Бюргерса применяются два подхода: методы на сетках конечных разностей и спектральные методы. Первый подход основывается на локальной аппроксимации производных, а второй — на представлении решения в виде разложения по базисным функциям (обычно тригонометрическим). Конечно-разностные методы широко распространены благодаря простоте реализации, однако страдают от численной диффузии и ограниченной точности при моделировании гладких решений. Спектральные методы, напротив, обеспечивают крайне высокую точность при малом числе узлов сетки, но требуют строгих условий гладкости и периодичности решения. Сравнение этих методов представляет научный и практический интерес, поскольку выбор схемы определяет эффективность моделирования физических процессов, включая турбулентность, взаимодействие волн и диссипативные явления. Цель данной работы заключается в сравнительном исследовании спектрального и конечно-разностного методов для решения уравнения Бюргерса в контексте моделирования газовых потоков, а также в выявлении преимуществ и ограничений каждого из них при различных параметрах вязкости и шагах сетки. АНАЛИЗ И МЕТОДОЛОГИЯ ЛИТЕРАТУРЫ Первые аналитические исследования уравнения Бюргерса были проведены Дж. Бюргерсом (Burgers, 1948), который показал, что данное уравнение может служить моделью турбулентных течений с учетом вязкости. Позднее Хопф (Hopf, 1950) и Коул (Cole, 1951) предложили преобразование, известное как преобразование Коула — Хопфа, позволяющее линеаризовать уравнение и получать точные решения. Эти работы заложили основу для использования уравнения Бюргерса в качестве тестовой задачи при проверке точности и устойчивости численных схем. Развитие спектральных методов связано с работами Орсага (Orszag, 1971) и Готтлиба (Gottlieb, 1977), которые показали, что представление решения через ряды Фурье обеспечивает ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 46 экспоненциальную сходимость для гладких функций. Пейре (Peyret, 2002) и Бойд (Boyd, 2001) подробно описали применение спектральных и псевдоспектральных методов для решения задач гидродинамики и показали их преимущества перед методами конечных разностей при периодических граничных условиях. В то же время классические конечно-разностные методы, рассмотренные Флетчером (Fletcher, 1991) и Хиршем (Hirsch, 2007), остаются основным инструментом при решении задач с разрывами и сложной геометрией. Эти методы устойчивы и универсальны, но требуют значительно большего числа узлов для достижения той же точности, что и спектральные методы. В последние десятилетия активно исследуются гибридные методы (Karniadakis, 2013), объединяющие преимущества обоих подходов. Современные исследования (Xu, 2015; Keeling, 2018) демонстрируют, что спектральные методы особенно эффективны для моделирования переходных процессов и турбулентных структур в газах при умеренных числах Рейнольдса. Тем не менее, при наличии разрывов (ударных волн) эффективность спектральных методов снижается из-за эффекта Гиббса. Для борьбы с этим применяются фильтрация высокочастотных мод и регуляризация решений. Таким образом, анализ литературы показывает, что оба метода имеют свои области применения: спектральные методы — при моделировании гладких, периодических решений, а конечно-разностные — при задачах с локальными особенностями, не требующих глобального базиса. В настоящей работе проводится их количественное сравнение по параметрам точности, устойчивости и вычислительных затрат. Для численного решения уравнения Бюргерса рассмотрены два подхода — спектральный и конечно-разностный. Расчёты выполняются на одинаковом временном интервале и при одинаковых условиях вязкости. В качестве тестовой задачи используется периодическая область длиной L=2π с начальными условиями u(x,0)=sin(x). В спектральном методе решение представляется в виде ряда Фурье: где — спектральные коэффициенты. Нелинейный член вычисляется в физическом пространстве с помощью быстрого преобразования Фурье (FFT), после чего производится обратное преобразование для получения спектра. Вязкостный член интегрируется аналитически, что обеспечивает высокую точность и стабильность при малых шагах по времени. Для интегрирования по времени применяется метод Рунге–Кутта четвёртого порядка. ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 47 В конечно-разностном методе используются центральные аппроксимации второго порядка для пространственных производных: , . Для интегрирования во времени применяется схема Кранка–Николсона или метод МакКормика. Шаг по времени выбирается из условия устойчивости Куранта–Фридрихса–Леви (CFL). Сравнение методов проводится по следующим критериям: 1. Погрешность по нормам L2 и L∞; 2. Скорость сходимости при увеличении числа узлов; 3. Вычислительная сложность (число операций и время счета); 4. Поведение при малых значениях вязкости (формирование ударной волны). Для визуализации результатов строятся профили скорости u(x,t) в различные моменты времени и энергетические спектры E(k). Расчеты выполняются в программной среде Python с использованием библиотек NumPy и SciPy. ПОЛУЧЕННЫЕ РЕЗУЛЬТАТЫ Результаты численных экспериментов показали, что спектральный метод обеспечивает значительно более высокую точность при меньшем числе узлов сетки. При N=64 точность спектрального решения сопоставима с конечно-разностным решением при N=512. Нормы ошибок демонстрируют экспоненциальное снижение для спектрального метода и лишь степенное для конечно-разностного. При малых коэффициентах вязкости (ν=0.01) спектральный метод чётко воспроизводит образование ударной волны, но при дальнейшем снижении вязкости возникает эффект Гиббса — осцилляции около разрывов. Конечно-разностная схема при тех же условиях показывает более сглаженный профиль за счёт численной диффузии, но теряет точность в центральной части волны. Энергетический спектр, полученный спектральным методом, демонстрирует правильное распределение энергии по модам, что подтверждает его способность описывать турбулентные переходы. При увеличении числа узлов спектр совпадает с теоретическим законом затухания. Конечно-разностная аппроксимация демонстрирует избыточное затухание высокочастотных мод, что указывает на наличие искусственной вязкости. По вычислительным затратам спектральный метод требует больше операций на одну итерацию, но общее время решения оказывается меньше из-за более быстрой сходимости. Таким образом, при гладких решениях спектральный подход более эффективен. Однако при наличии разрывов или сложной геометрии предпочтительнее конечно-разностный метод, обеспечивающий лучшую устойчивость. ВЫВОД ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 48 Проведённое исследование показало, что выбор численного метода для решения нелинейных газовых уравнений типа Бюргерса должен определяться характером задачи. Спектральные методы обеспечивают высокую точность и экспоненциальную сходимость при моделировании гладких решений, что делает их предпочтительными для задач с периодическими граничными условиями и умеренной нелинейностью. Конечно-разностные методы обладают большей универсальностью и устойчивостью при решении задач с разрывами и нестационарными граничными условиями. Численные эксперименты подтвердили, что при одинаковом числе узлов спектральный метод превосходит конечно-разностный по точности на несколько порядков. Однако в областях с сильной нелинейностью или при низкой вязкости, где формируются ударные волны, конечноразностные схемы показывают более стабильные результаты. Таким образом, комбинированное использование спектральных и конечно-разностных подходов может стать оптимальным решением для моделирования газодинамических систем, где требуется высокая точность в гладких областях и устойчивость вблизи разрывов. Перспективными направлениями дальнейших исследований являются разработка гибридных методов, адаптивных сеток и параллельных вычислительных алгоритмов, обеспечивающих баланс между скоростью, точностью и устойчивостью численного моделирования. ИСПОЛЬЗОВАННАЯ ЛИТЕРАТУРА 1. Burgers J. M. A mathematical model illustrating the theory of turbulence // Advances in Applied Mechanics. – 1948. – Vol. 1. – P. 171–199. 2. Hopf E. The partial differential equation u_t + u u_x = μ u_xx // Communications on Pure and Applied Mathematics. – 1950. – Vol. 3. – P. 201–230. 3. Cole J. D. On a quasilinear parabolic equation occurring in aerodynamics // Quarterly of Applied Mathematics. – 1951. – Vol. 9. – P. 225–236. 4. Orszag S. A. Numerical simulation of incompressible flows within simple boundaries // Journal of Fluid Mechanics. – 1971. – Vol. 49. – P. 75–112. 5. Gottlieb D., Orszag S. A. Numerical analysis of spectral methods: theory and applications. – Philadelphia: SIAM, 1977. – 372 p. 6. Fletcher C. A. J. Computational techniques for fluid dynamics. – Berlin: Springer-Verlag, 1991. – 502 p. 7. Peyret R. Spectral methods for incompressible viscous flow. – Berlin: Springer, 2002. – 358 p. 8. Boyd J. P. Chebyshev and Fourier spectral methods. – New York: Dover Publications, 2001. – 668 p. 9. Canuto C., Hussaini M. Y., Quarteroni A., Zang T. A. Spectral methods: fundamentals in single domains. – Berlin: Springer, 2006. – 563 p. 10. Karniadakis G. E., Sherwin S. Spectral/hp element methods for CFD. – Oxford: Oxford University Press, 2013. – 641 p. 11. Hirsch C. Numerical computation of internal and external flows. – Oxford: Elsevier, 2007. – 730 p. 12. Xu H., Chen Q., Li Z. Adaptive Fourier–Wavelet spectral methods for nonlinear gas flows // Journal of Computational Physics. – 2015. – Vol. 297. – P. 256–274. ISSN: 2582-4686 SJIF 2021-3.261,SJIF 20222.889, 2024-6.875 ResearchBib IF: 9.948 / 2024 VOLUME-5, ISSUE-11 49 13. Keeling S., Zaki T. Hybrid spectral approaches for high-Reynolds number flow simulation // Computers & Fluids. – 2018. – Vol. 173. – P. 154–166. 14. Shen J. Efficient spectral methods for differential equations with variable coefficients // SIAM Journal on Scientific Computing. – 2011. – Vol. 33 (1). – P. 96–118.