flux unity 流体_【游戏流体力学基础及Unity代码（八）】激波捕捉法和RiemannSolver...

本篇介绍各种解决burgers方程的数值方法，相应的代码可在clatterrr/CFDcodepython找到。

上一章介绍了最简单的非线性方程，即一维无粘性burgers方程，但我们将它稍微改写一下

微分项用泰勒展开微分实际上得到了

上式是精确解。但是把右侧的值全算出来不现实。如果只取等式右侧第一项，就是有限差分(finite difference)形式的迎风格式(Upwind Scheme)。

一阶upwind格式的缺点就是舍弃了太多泰勒展开项，所以很不精确。而有限差分法规定每个网格由一个离散的值表示，在处理有关激波的不连续问题时并不那么自然，方便。所以有限体积(finite volume)法便被提了出来，这种方法将一个格子里的值看成连续的，从前到后积分后再除以格子长度就是格子的值，这种方法在处理不连续问题时有很大优势。比如最简单的有限体积法下的Lax-Friedrichs方法：

然而这种方法也只有一阶精度，并且有一些人工粘性(artifial viscosity)所导致的耗散(disspative)。它看起来像是这样的。虚线是准确的解，而其它颜色的线是不同情况下的Lax-Friedrhs。人工粘性就像高斯模糊样，让尖锐的部分变得平滑了，然而我们要模拟的激波也消失了，这显然不是一种很好的方法。

所以我们继续尝试新方法，即Lax-Wendroff方法: