Advances in Applied Mathematics
Vol.
12
No.
01
(
2023
), Article ID:
60952
,
9
pages
10.12677/AAM.2023.121029
基于通量限制器的时间分数阶Burgers方程数值解法
魏英岚,白慧冉
长沙理工大学数学与计算科学学院,湖南 长沙
收稿日期:2022年12月28日;录用日期:2023年1月24日;发布日期:2023年1月31日
![](http://html.hanspub.org/file/29-2623091x1_hanspub.png?20230202083205945)
摘要
对于求解可能具有激波等不连续点的时间分数阶守恒律,当分数阶γ接近于0时,目前还没有有效的方法求解时间分数阶非线性离散系统,本文以时间分数阶Burgers方程为例,运用多重网格迭代方法进行求解,对于对流项,采用通量限制器,使得新的数值通量在光滑区域变为高阶通量而在间断附近变为低阶通量,从而使问题的解达到更高阶精度,并在不同的γ取值以及不同的初边值条件下进行了有效的数值实验。
关键词
时间分数阶Burgers方程,多重网格法,通量限制器
![](http://html.hanspub.org/file/29-2623091x2_hanspub.png?20230202083205945)
Numerical Solution of Time-Fractional Burgers Equation Based on Flux Limiter
Yinglan Wei, Huiran Bai
College of Mathematics and Computational Science, Changsha University of Science and Technology, Changsha Hunan
Received: Dec. 28th, 2022; accepted: Jan. 24th, 2023; published: Jan. 31st, 2023
![](http://html.hanspub.org/file/29-2623091x3_hanspub.png?20230202083205945)
ABSTRACT
For solving the fractional time conservation law that may have discontinuities such as shock waves, when the fractional order γ is close to 0, there is no effective method to solve the fractional nonlinear discrete system. Taking the time fractional Burgers equation as an example, this paper applies the multi-grid iterative method to solve the problem. For the convection term, a flux limiter is adopted to make the new numerical flux become a higher order flux in the smooth area and a lower order flux near the discontinuity, thus achieving a higher order precision. The effective numerical experiments are carried out with different values of γ and different initial boundary values.
Keywords:Time Fractional Burgers Equation, Multiple Grid Method, Flux Limiter
![](http://html.hanspub.org/file/29-2623091x4_hanspub.png?20230202083205945)
Copyright © 2023 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
1. 引言
时间分数阶标量守恒律方程
(1)
当 时,我们称其为时间分数阶Burgers方程。
本文考虑下述时间分数阶Burgers方程:
(2)
其中 ,。
满足初值条件
以及合适的边界条件,其中 , 为给定函数,
为 阶Caputo分数阶导数, 为Gamma函数。
近十年来,许多学者运用不同的数值方法对含黏时间分数阶Burgers方程进行了求解,例如,Duangpan运用了有限积分法结合移位Chebyshev多项式 [1],Akram等提出了有限差分格式 [2],Esen和Tasbozan提出了三次b样条有限元搭配方法 [3] 等等。
2. 离散
定义空间步长 ,空间网格节点 ,其中 ;定义空间网格半节点 ,得到 个网格单元 ,其中 。定义时间步长 ,时间网格节点 ,其中 。
对于函数 ,使用如下简写形式:
在 处计算式(1.1)并在单元 上对方程进行积分,得到
(3)
其中 表示 在区间 上的平均值。
定义 ,对式(3)的左端第一项运用 公式逼近,取 [4]:
(4)
通量近似为
(5)
(6)
结合式(3)、(4)、(5)及(6),得到最终的离散形式为
(7)
3. 通量限制器
为了使问题的解达到更高阶精度,我们首先选择一个在光滑区域表现很好的高阶通量和一个在间断附近表现很好的低阶通量,然后将数值通量重新定义为低阶通量和高阶通量的线性组合,使得新的数值通量在光滑区域变为高阶通量而在间断附近变为低阶通量:
式中 为低阶通量, 为高阶通量; 为通量限制器,是变差比 的函数,即
的取值与流动方向相关,对于均匀网格,定义
可以看到,当 时,通量为低阶格式;当 时,通量为高阶格式。
为保证格式稳定性, 是必要的;同时,格式需要满足二阶精度,即对于限制器 有如下约束条件(若 ,则令 ):
我们把Lax-Friedrichs通量作为低阶通量:
如果 取常数 ,令 ,则称Lax-Friedrichs通量为全局Lax-Friedrichs通量;如果 ,则称其为局部Lax-Friedrichs通量,其中 ,是一个常数标量。
使用Lax-Wendroff通量为高阶通量:
下面列举出一些常见的通量限制器 [5]:
1) Chakravarthy-Osher限制器:
当 时,称其为minmod限制器;
2) Koren限制器:
3) Sweby限制器:
当 时,称其为Superbee限制器;
4) OSPRE限制器:
5) van Leer限制器:
限制器 一般是对称的,即满足以下等式:
4. 多重网格法
本文使用多重网格迭代方法,由于空间系统是强非线性的,我们的多重网格方法将以线性多重网格方法的非线性推广——FAS多重网格方法为基础 [6]。
假设三个正整数 满足式:
构建空间网格层次结构: ,假设最细网格 的网格大小为 ,对于 , 的网格大小为 ,其中 表示 的网格数量, ,下文中与 有关的量以上标 的形式书写,例如 。
在多重网格法中,我们通过求解粗网格 上的离散方程对迭代解进行反复修正:
其中 为非线性算子。
对式(7)采用Richardson迭代方法求解:
(8)
其中 , 和 分别为一次迭代过程中的新近似解和旧近似解, 为一个能够保证迭代方法收敛性的合适的正数。
我们定义延拓和限制算子来实现相邻网格之间的数据传输,其中,延拓算子 实现 的数据传输,即
定义为:
限制算子 实现 的数据传输,即
定义为:
以下提出求解方程的多重网格迭代方法步骤:
定义 :
步1:在最网格上计算非线性离散方程 ,得到方程近似解 ;
步2:在粗网格上计算非线性残差方程 ;
步3:得到新的右端项 ;
步4:在粗网格上对近似解进行修正 。
5. 数值实验
例5.1设计算域 ,初值条件为 ,边界条件为 ,取 ,,分别使用OSPRE limiter (见图1)和Koren limiter (见图2)求解在 时方程(2)的解:
Figure 1. OSPRE limiter
图1. OSPRE限制器
Figure 2. Koren limiter
图2. Koren限制器
例5.2设计算域 ,初值条件为 ,边界条件为 ,取 ,,分别使用Koren limiter (见图3)和van Leer limiter (见图4)求解在 时方程(2)的解:
Figure 3. Koren limiter
图3. Koren限制器
Figure 4. Van Leer limiter
图4. Van Leer限制器
例5.3设计算域 ,初值条件为 ,边界条件为 ,,取 ,,分别使用OSPRE limiter (见图5)和van Leer limiter (见图6)求解在 时方程(2)的解:
Figure 5. OSPRE limiter
图5. OSPRE限制器
Figure 6. Van Leer limiter
图6. Van Leer限制器
6. 总结
本文运用多重网格迭代方法求解时间分数阶Burgers方程,对于对流项,引入了通量限制器,在三种不同的初值和边界条件下,选用了两种通量限制器,分别在 时方程的解,得到了较好的结果。
文章引用
魏英岚,白慧冉. 基于通量限制器的时间分数阶Burgers方程数值解法
Numerical Solution of Time-Fractional Burgers Equation Based on Flux Limiter[J]. 应用数学进展, 2023, 12(01): 254-262. https://doi.org/10.12677/AAM.2023.121029
参考文献
- 1. Duangpan, A., Boonklurb, R. and Treeyaprasert, T. (2019) Finite Integration Method with Shifted Chebyshev Polynomi-als for Solving Time-Fractional Burgers’ Equations. Mathematics, 7, 1201. https://doi.org/10.3390/math7121201
- 2. Akram, T., Abbas, M., Riaz, M.B., et al. (2020) An Efficient Numerical Technique for Solving Time Fractional Burgers Equation. Alexandria Engineering Journal, 59, 2201-2220. https://doi.org/10.1016/j.aej.2020.01.048
- 3. Esen, A., Bulut, F. and Oruç, Ö. (2016) A Unified Approach for the Numerical Solution of Time Fractional Burgers’ Type Equations. The European Physical Journal Plus, 131, 1-13. https://doi.org/10.1140/epjp/i2016-16116-5
- 4. Langlands, T.A.M. and Henry, B.I. (2005) The Accuracy and Sta-bility of an Implicit Solution Method for the Fractional Diffusion Equation. Journal of Computational Physics, 205, 719-736. https://doi.org/10.1016/j.jcp.2004.11.025
- 5. Hesthaven, J.S. (2017) Numerical Methods for Conserva-tion Laws: From Analysis to Algorithms. In: Computational Science & Engineering, xvi + 555. https://doi.org/10.1137/1.9781611975109
- 6. Briggs, W.L., Henson, V.E. and McCormick, S.F. (2000) A Multi-grid Tutorial. In: Other Titles in Applied Mathematics, xii + 187. https://doi.org/10.1137/1.9780898719505