The purpose of this paper is to propose a new algorithm for obtaining approximate solutions to the Burgers' equation (BE). Integration in time by a quadratic B-spline collocation method is shown. To the best of our knowledge, B-splines have not previously been used to integrate partial differential equations in both time and space. First, the BE is integrated using quadratic B-spline functions in time, and then the time-integrated BE is further solved in space via the cubic B-spline collocation method. The resulting recursive algebraic equation is used to obtain both shock wave and front propagation solutions of the BE, demonstrating the effectiveness of the space--time B-spline collocation method.
- 论文ID: 2510.09408
- 标题: A time-space B-spline integrator for the Burgers' equation
- 作者: Idris Dag, Serkan Uğurluoğlu, Nihat Adar (Eskisehir Osmangazi University, Turkey)
- 分类: math.NA cs.NA
- 发表时间: 2025年10月10日 (arXiv预印本)
- 论文链接: https://arxiv.org/abs/2510.09408
本文提出了一种求解Burgers方程近似解的新算法。该方法采用二次B样条配置方法进行时间积分。据作者所知，这是首次同时在时间和空间两个维度上使用B样条来积分偏微分方程。首先使用二次B样条函数在时间方向上积分Burgers方程，然后通过三次B样条配置方法在空间方向上进一步求解时间积分后的方程。所得到的递归代数方程用于获得Burgers方程的激波和前沿传播解，证明了时空B样条配置方法的有效性。
- Burgers方程的重要性：Burgers方程是一个准线性抛物型偏微分方程，在工程科学中广泛应用于建模湍流、气体动力学、热传导、弹性、交通流和声波等现象。它提供了描述对流和扩散相互作用的数学表述。
- 数值求解的必要性：非线性偏微分方程通常不存在解析解，特别是当粘性常数υ较小时，级数解收敛缓慢，实用性差。小的υ值会导致激波形成，使得解析和数值方法都难以提供合理的解。
- 现有方法的局限：传统的时间积分器包括Crank-Nicolson格式、Runge-Kutta方法、几何积分器、指数积分器等，但这些方法在处理Burgers方程的激波和前沿传播问题时存在精度和稳定性方面的挑战。
本文的核心动机是开发一种新的数值方法，能够：
- 同时在时间和空间维度上使用B样条进行积分
- 有效处理小粘性值情况下的激波问题
- 提供更高精度的数值解
- 首创性方法：首次提出同时使用B样条作为时间和空间积分器求解Burgers方程的方法
- 双重B样条策略：在时间方向使用二次B样条，在空间方向使用三次B样条的创新组合
- 激波和前沿传播模拟：成功模拟了Burgers方程的激波解和行波前沿解
- 递归算法框架：建立了基于B样条配置的递归代数方程求解体系
求解一维时间相关的Burgers方程：
ut+uux−νuxx=0,x∈[a,b],t∈(0,T]
配合初始条件：u(x,0)=f(x)
边界条件：u(a,t)=u0,ux(a,t)=u0′,u(b,t)=um,ux(b,t)=um′
其中ν是粘性常数。
- 网格划分：将时间域[c,d]等分为网格点c=t0<t1<⋯<tk=d，其中tj=t0+jΔt
- 二次B样条基函数：定义Bj(t)∈C1[a,b]，j=−1,…,n：
(t_{j+2}-t)^2 - 3(t_{j+1}-t)^2 + 3(t_j-t)^2, & t_{j-1} \leq t \leq t_j \\
(t_{j+2}-t)^2 - 3(t_{j+1}-t)^2, & t_j \leq t \leq t_{j+1} \\
(t_{j+2}-t)^2, & t_{j+1} \leq t \leq t_{j+2} \\
0, & \text{其他}
\end{cases}$$
3. **近似解表示**：
$$U(x,t) = \sum_{j=-1}^n \delta(x,t_j)B_j(t)$$
#### 空间离散化（三次B样条）
1. **空间网格**：将空间域$[a,b]$等分，$x_i = x_0 + ih$，$i = 0,\ldots,m$
2. **三次B样条基函数**：$B_i(x) \in C^2[a,b]$，$i = -1,\ldots,m$
3. **参数展开**：
$$\delta(x,t^j) = \sum_{i=-1}^{m+1} \sigma_i^j B_i(x)$$
### 技术创新点
1. **分层积分策略**：先在时间方向用二次B样条积分，再在空间方向用三次B样条求解
2. **配置点法**：在网格点处强制满足微分方程
3. **线性化处理**：通过选择前一时间步的值来线性化非线性项
4. **迭代求精**：采用迭代过程提高每个时间步的精度
## 实验设置
### 测试问题
#### 问题1：激波解
解析解：
$$u(x,t) = \frac{x/t}{1 + \sqrt{t/t_0}\exp(x^2/(4\nu t))}, \quad t \geq 1, \quad 0 \leq x \leq 1$$
其中$t_0 = \exp(1/(8\nu))$
#### 问题2：行波前沿解
解析解：
$$u(x,t) = \frac{\alpha + \mu + (\mu - \alpha)\exp(\eta)}{1 + \exp(\eta)}$$
其中$\eta = \frac{\alpha(x - \mu t - \gamma)}{\nu}$，$\alpha = 0.4$，$\mu = 0.6$，$\gamma = 0.125$
### 评价指标
- **$L_\infty$范数**：$|U - U_N|_\infty = \max|U_i^j - (U_N)_i^j|$
- **$L_2$范数**：$|U - U_N|_2 = \sqrt{\sum_i |U_i^j - (U_N)_i^j|^2}$
- **峰值**：激波的最大值
### 实现细节
- 时空步长：$h = \Delta t = 0.01, 0.005, 0.001$
- 粘性参数：$\nu = 0.01, 0.005$
- 时间范围：$t = 3.25$（激波问题），$t = 1.2$（前沿波问题）
## 实验结果
### 主要结果
#### 激波问题结果
| $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ | Peak value |
|-------------|--------------|------------|-------|------------|
| | 0.01 | 0.0498 | 0.0211| 0.2233 |
| | 0.005 | 0.0508 | 0.0207| 0.2225 |
| | 0.001 | 0.0516 | 0.0206| 0.2203 |
| $\nu=0.005$ | $L_\infty$ | $L_2$ | Peak value |
|--------------|------------|-------|------------|
| | 0.0582 | 0.0166| 0.2503 |
| | 0.0571 | 0.0159| 0.2480 |
| | 0.0560 | 0.0155| 0.2467 |
#### 前沿波问题结果
| $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ |
|-------------|--------------|------------|-------|
| | 0.01 | 0.4875 | 0.1138|
| | 0.005 | 0.2964 | 0.0638|
| | 0.001 | 0.1514 | 0.0354|
| $\nu=0.005$ | $L_\infty$ | $L_2$ |
|--------------|------------|-------|
| | 0.8841 | 0.2532|
| | 0.6906 | 0.1329|
| | 0.2427 | 0.0367|
### 实验发现
1. **精度与网格关系**：使用更小的时空步长能获得更光滑的解和更好的精度
2. **粘性影响**：较小的粘性常数导致更尖锐的波形和更高的误差
3. **误差分布**：最大误差通常出现在激波峰值附近或波前区域
4. **稳定性**：方法在不同参数设置下都保持了良好的数值稳定性
## 相关工作
该领域的主要研究方向包括：
1. **传统数值方法**：有限差分、有限元、谱方法等
2. **B样条方法**：主要用于空间离散，很少用于时间积分
3. **配置方法**：Galerkin方法、配置方法、无网格方法等
4. **专用求解器**：针对Burgers方程的特殊数值方法
本文的创新在于首次将B样条同时应用于时间和空间积分，填补了该领域的空白。
## 结论与讨论
### 主要结论
1. 提出的时空B样条配置方法能够有效求解Burgers方程
2. 方法在激波和前沿传播问题上都表现出良好的精度和稳定性
3. 通过调节时空步长可以控制解的精度和光滑性
### 局限性
1. **计算复杂度**：需要在每个时间步求解线性方程组
2. **非线性处理**：需要线性化和迭代过程
3. **参数敏感性**：对粘性参数和网格参数较为敏感
4. **适用范围**：主要针对一维Burgers方程
### 未来方向
作者提出未来工作将重点关注：
- 使用更高阶B样条提高精度
- 扩展到多维Burgers方程
- 优化算法效率
## 深度评价
### 优点
1. **方法创新性**：首次提出时空B样条积分方法，具有原创性
2. **理论完整性**：从数学推导到数值实现都较为完整
3. **实验充分性**：测试了激波和前沿传播两类典型问题
4. **精度分析**：提供了详细的误差分析和收敛性研究
### 不足
1. **理论分析不足**：缺乏收敛性和稳定性的严格理论证明
2. **对比实验有限**：未与其他先进方法进行充分对比
3. **计算效率分析缺失**：未提供计算时间和内存使用分析
4. **参数选择指导不足**：缺乏B样条阶数选择的理论指导
### 影响力
1. **学术贡献**：为PDE数值求解提供了新的思路
2. **实用价值**：方法简单易实现，适合工程应用
3. **可扩展性**：框架可推广到其他非线性PDE
### 适用场景
该方法特别适用于：
- 需要高精度求解Burgers方程的科学计算
- 激波和前沿传播现象的数值模拟
- 对计算精度要求较高但对效率要求适中的应用
## 参考文献
论文引用了16篇相关文献，涵盖了Burgers方程的数值方法、B样条理论、配置方法等相关领域的重要工作，为本研究提供了坚实的理论基础。