一、项目背景详细介绍
在计算流体力学（CFD）、数值天气预报、海洋模拟以及输运问题中，
平流（Advection / Convection）方程是最基础、最核心的模型之一。
它描述的是：
某种物理量在给定速度场作用下，被“整体搬运”的过程
典型例子包括：
污染物随风场移动
温度随流体整体迁移
浓度、密度的输运
数值方法中误差传播的“放大镜模型”
为什么平流方程如此重要？
因为它具有以下特征：
方程形式极其简单
数值解却极其敏感
能直接暴露：
数值耗散
数值色散
相位误差
稳定性问题
因此，在 CFD 教学体系中：
“平流方程 = 数值格式的试金石”
本项目研究对象
即：
初始波形以速度 ccc 不变形地向右传播
为什么选择 Lax–Wendroff 方法？
在众多平流格式中：
一阶迎风：稳定但耗散大
中心差分：不稳定
Leapfrog：无耗散但振荡
Lax–Wendroff：
显式
二阶精度
兼顾稳定性与精度
因此：
Lax–Wendroff 是线性平流问题中“最经典的二阶格式”
二、项目需求详细介绍
2.3 功能需求
使用有限差分法离散空间导数
使用Lax–Wendroff 方法推进时间
处理一维恒速平流问题
正确实现周期边界条件
满足 CFL 稳定性条件
输出数值解用于对比解析解
三、相关技术详细介绍
3.1 一维平流方程的数学本质
3.1.1 双曲型偏微分方程
平流方程属于：
线性双曲型偏微分方程
其特点是：
信息沿特征线传播
特征线为：
3.1.2 精确解特性
解析解为平移：
这意味着：
无耗散
无变形
波形只发生相位移动
任何数值偏差：
都来自数值格式本身
该条件是：
线性平流方程数值稳定的核心约束
3.6 Lax–Wendroff 方法的数值特性
| 特性 | 表现 |
|---|---|
| 精度 | 时间、空间二阶 |
| 稳定性 | 条件稳定 |
| 数值耗散 | 较小 |
| 数值色散 | 存在 |
| 工程用途 | 教学 / 基准 |
四、实现思路详细介绍
4.1 整体算法流程
对空间区间进行均匀网格划分
初始化初始条件
选择满足 CFL 的时间步长
使用 Lax–Wendroff 格式推进时间
每个时间步施加周期边界条件
输出最终或中间结果
4.2 周期边界条件实现
左端点：引用右端相邻节点
右端点：引用左端相邻节点
通过取模索引实现。
4.3 数值行为预期
波形整体向右平移
振幅基本保持
波峰附近可能出现轻微振荡（色散）
五、完整实现代码
/**************************************************** * 文件名：Advection1D_LaxWendroff.cpp * 描述：C++ 使用有限差分法 + Lax–Wendroff 方法 * 求解一维恒速平流方程 ****************************************************/ #include <iostream> #include <vector> #include <cmath> using namespace std; /**************************************************** * 主函数 ****************************************************/ int main() { // 空间参数 int Nx = 200; double a = 0.0, b = 1.0; double dx = (b - a) / Nx; // 时间参数 double c = 1.0; // 平流速度 double dt = 0.002; // 时间步长 double T = 1.0; // 总时间 // CFL 条件检查 double CFL = fabs(c) * dt / dx; if (CFL > 1.0) { cout << "CFL 条件不满足，程序终止" << endl; return -1; } int Nt = static_cast<int>(T / dt); // 空间网格 vector<double> x(Nx); for (int i = 0; i < Nx; ++i) x[i] = a + i * dx; // 数值解 vector<double> u_curr(Nx, 0.0); vector<double> u_next(Nx, 0.0); // 初始条件 u(x,0) = sin(2πx) for (int i = 0; i < Nx; ++i) u_curr[i] = sin(2.0 * M_PI * x[i]); // 时间推进 for (int n = 0; n < Nt; ++n) { for (int i = 0; i < Nx; ++i) { int ip = (i + 1) % Nx; int im = (i - 1 + Nx) % Nx; u_next[i] = u_curr[i] - 0.5 * CFL * (u_curr[ip] - u_curr[im]) + 0.5 * CFL * CFL * (u_curr[ip] - 2.0 * u_curr[i] + u_curr[im]); } u_curr = u_next; } // 输出结果 cout << "x u(x,T)" << endl; for (int i = 0; i < Nx; ++i) cout << x[i] << " " << u_curr[i] << endl; return 0; }六、代码详细解读（仅解读方法作用）
u_curr/u_next：当前与下一时间层解中心差分：空间一阶、二阶导数近似
Lax–Wendroff 更新公式：二阶时间推进
周期边界条件：通过取模索引实现
CFL 检查：保证稳定性
七、项目详细总结
通过该项目，你已经系统掌握：
一维平流方程的数学与物理本质
Lax–Wendroff 方法的完整推导思想
二阶显式格式在平流问题中的表现
数值色散与相位误差的来源
平流问题作为“数值格式放大镜”的作用
这是从：
一阶迎风 → 二阶经典格式
的关键过渡，也是深入TVD / WENO / 高分辨率格式的坚实基础。
八、项目常见问题及解答
Q1：为什么会出现振荡？
A：Lax–Wendroff 是非 TVD 格式，存在数值色散。
Q2：可以完全无失真吗？
A：只能通过更高阶或限制器方法逼近。
Q3：工程 CFD 会用它吗？
A：更多用于教学与基准测试。
九、扩展方向与性能优化
一阶迎风 vs Lax–Wendroff 对比
TVD 平流格式（Minmod / Superbee）
MUSCL 二阶迎风
WENO 平流方程
非恒定速度平流方程
