有限差分法求解不可压NS方程

网上关于有限差分法解NS方程的程序实现不尽完备，这里是一些补充注解

现有的优秀资料

理论向

【1】如何从物理意义上理解NS方程？ - 知乎

【2】NS方程数值解法：投影法的简单应用 - 知乎

【3】[计算流体力学] NS 方程的速度压力法差分格式_「已注销」的博客-CSDN博客

【4】有限差分法（Finite Difference Method）解方程：边界和内部结点的控制方程 - 简书

编程向

【5】【游戏流体力学基础及Unity代码（五）】用欧拉方程模拟无粘性染料之代码实现 - 知乎

【6】（CFD）投影法求解二维不可压缩N-S方程_陈天文的博客-CSDN博客_二维ns方程公式

【7】CFD101: 2D Lid Driven Cavity Flow - File Exchange - MATLAB Central

有限差分法求解流程按照【2】，网格离散方法按照【7】，压力泊松方程迭代按照【5】。

注1：速度定义在网格线上，压力定义在网格中心。对于整个流体域，为了编程方便起见，在四周多加一层网格（虚拟网格）。但并不是所有虚拟网格上都定义了速度，见下图

x方向速度定义在网格左边缘，而我们计算中不会用到左侧虚拟网格的速度，所以不需要定义第1列（x=1）的x方向速度。而对y方向速度，用第1列（x=1/2）y方向速度的原因是要算x=2处的y方向速度（x=1/2与x=3/2做平均）。

注2：方腔顶盖流似乎是某种标准的仿真验证工况，但是如果只验证这一种可能会发现不了一些问题。例如上面的【6】，无法计算左侧边界流入的情形，因为（对流项的）后向差分方法没法将左侧边界纳入计算。

对流项使用【7】中的散度形式似乎能解决这个问题，它使用了更多的平均，因此可以把周边的节点都考虑进来。

注3：方腔顶盖流只涉及了压力梯度为0和速度为0的边界条件。在每次压力迭代和速度迭代之前更新边界条件。对于边界法向的流动入口和出口，实测压力梯度为0可以保留，入口边界直接定义速度，出口边界速度梯度为0。但是还不清楚对于压力入口和滞止出口如何定义。

有限差分法求解不可压NS方程相关推荐

C++：应用有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上，与恒定速度，使用FTCS方法，正向时间差，居中空间差(附完整源码)
C++:应用有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上,与恒定速度,使用FTCS方法,正向时间差, 居中空间差 # include <cmath> ...
C++：有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上，与恒定速度，使用Lax-Wendroff方法作为时间导数(附完整源码)
C++:有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上,与恒定速度,使用Lax-Wendroff方法作为时间导数 # include <cstdlib> ...
C++：有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上，与恒定速度，使用Lax方法作为时间导数(附完整源码)
C++:有限差分法求解随时间变化平流方程 ut = - c * ux 在一个空间维度上,与恒定速度,使用Lax方法作为时间导数 # include <cstdlib> # includ ...
工程流体力学笔记暂记24 (不可压缩粘性流体的运动微分方程**N-S方程**)
粘性流体牛顿第二动律分析x表面的受力情况为例当气体各流层之间有相对运动时,存在一种阻滞这种运动的力,一旦相对运动消失,这种力也不复存在,气体的这种性质称为粘性.其作用力称为粘性力,粘性力时和它的 ...
有限差分法求解一维非稳态对流扩散方程
此文转载自本人在简书平台上发布的文章 1. 源起最近想用有限差分法计算二维的顶盖驱动流,在推导过程中遇到了许多问题,例如对流项的离散."线性化"这一在有限体积法中很常见的操作. ...
二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维 Poisson 方程的 MATLAB 实现陈莲a ,郭元辉b ,邹叶童a ( 西华师范大学 a．数学与信息学院; b．教育信息技术中心,四川南充 6437009) 摘要: ...
十二步解N-S方程之第四步
2019独角兽企业重金招聘Python工程师标准>>> 十二步解N-S方程之第四步通过前面的学习,我们已经知道如何去解线性和非线性的一维对流方程,而且也知道CFL是如何影响数值求解 ...
12步解N-S方程之第三步
2019独角兽企业重金招聘Python工程师标准>>> 12步解N-S方程之第三步在前2步中我们已经学会如何使用python编写程序解决简单的一维线性/非线性对流问题.也通过调节参 ...
12步解N-S方程之第五步（1）
2019独角兽企业重金招聘Python工程师标准>>> Burgers方程 Burgers方程是流体动力学的基本偏微分方程,在大量的应用领域都可以看到它的身影,如气动动力学模型.交通 ...

有限差分法求解不可压NS方程

现有的优秀资料

有限差分法求解不可压NS方程相关推荐

最新文章

热门文章