[计算流体力学] NS 方程的速度压力法差分格式

对于二维不可压流体，写出速度压力形式得连续性方程和运动方程；在如图所示的交错网格下，写出 p(1,j) 满足的差分格式

Latex

 重写为
u_{i+1/2,j}^{n+1}=u_{i+1/2,j}^n-\frac{\Delta t}{4\Delta x}\left[u_{i+\frac{3}{2},j}^n\left(u_{i+\frac{3}{2},j}^n+2u_{i+\frac{1}{2},j}^n\right)-u_{i-\frac{1}{2},j}^n\left(u_{i-\frac{1}{2},j}^n+2u_{i+\frac{1}{2},j}^n\right)\right]-\frac{\Delta t}{4\Delta y}\left[\left(u_{i+\frac{1}{2},j}^n+u_{i+\frac{1}{2},j+1}^n\right)\left(v_{i,j+\frac{1}{2}}^n+v_{i+1,j+\frac{1}{2}}^n\right)-\left(u_{i+\frac{1}{2},j-1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i,j-\frac{1}{2}}^n+v_{i+1,j-\frac{1}{2}}^n\right)\right]-\frac{1}{\rho}\frac{\Delta t}{\Delta x}\left(p_{I+1,J}^n-p_{I,J}^n\right)+\frac{\mu\Delta t}{\left(\Delta x\right)^2}\left(u_{i+\frac{3}{2},j}^n-2u_{i+\frac{1}{2},j}^n+u_{i-\frac{1}{2},j}^n\right)+\frac{\mu\Delta t}{\left(\Delta y\right)^2}\left(u_{i+\frac{1}{2},j+1}^n-2u_{i+\frac{1}{2},j}^n+u_{i+\frac{1}{2},j-1}^n\right)
v_{i,j+1/2}^{n+1}=v_{i,j+1/2}^n-\frac{\Delta t}{4\Delta x}\left[\left(u_{i+\frac{1}{2},j+1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i+1,j+\frac{1}{2}}^n+v_{i,j+\frac{1}{2}}^n\right)-\left(u_{i-\frac{1}{2},j+1}^n+u_{i-\frac{1}{2},j}^n\right)\left(v_{i,j+\frac{1}{2}}^n+v_{i-1,j-\frac{1}{2}}^n\right)\right]-\frac{\Delta t}{4\Delta y}\left[v_{i,j+\frac{3}{2}}^n\left(v_{i,j+\frac{3}{2}}^n+2v_{i,j+\frac{1}{2}}^n\right)-v_{i,j-\frac{1}{2}}^n\left(v_{i,j-\frac{1}{2}}^n+2v_{i,j+\frac{1}{2}}^n\right)\right]-\frac{1}{\rho}\frac{\Delta t}{\Delta y}\left(p_{I,J+1}^n-p_{I,J}^n\right)+\frac{\mu\Delta t}{\left(\Delta x\right)^2}\left(v_{i+1,j+\frac{1}{2}}^n-2v_{i,j+\frac{1}{2}}^n+v_{i-1,j+\frac{1}{2}}^n\right)+\frac{\mu\Delta t}{\left(\Delta y\right)^2}\left(v_{i,j+\frac{3}{2}}^n-2v_{i,j+\frac{1}{2}}^n+v_{i,j-\frac{1}{2}}^n\right)
\frac{1}{\left(\Delta x\right)^2}\left(p_{i+1,j}^n-2p_{i,j}^n+p_{i-1,j}^n\right)+\frac{1}{\left(\Delta y\right)^2}\left(p_{i,j+1}^n-2p_{i,j}^n+p_{i,j-1}^n\right)=\left(S_p\right)_{i,j}^n
\left(S_p\right)_{i,j}^n=\rho\left\{\frac{D_{i,j}^n}{\Delta t}-\frac{1}{{4\left(\Delta x\right)}^2}\left[\left(u_{i+\frac{3}{2},j}^n\right)^2-\left(u_{i+\frac{1}{2},j}^n\right)^2-\left(u_{i-\frac{1}{2},j}^n\right)^2+\left(u_{i-\frac{3}{2},j}^n\right)^2+2u_{i+\frac{1}{2},j}^n\left(u_{i+\frac{3}{2},j}^n-u_{i-\frac{1}{2},j}^n\right)+2u_{i-\frac{1}{2},j}^n\left(u_{i-\frac{3}{2},j}^n-u_{i+\frac{1}{2},j}^n\right)\right]-\frac{1}{2\Delta x\Delta y}\left[\left(u_{i+\frac{1}{2},j+1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i+1,j+\frac{1}{2}}^n-v_{i,j+\frac{1}{2}}^n\right)+\left(u_{i-\frac{1}{2},j}^n+u_{i-\frac{1}{2},j-1}^n\right)\left(v_{i,j-\frac{1}{2}}^n+v_{i-1,j-\frac{1}{2}}^n\right)\right]-\frac{1}{{4\left(\Delta y\right)}^2}\left[\left(v_{i,j+\frac{3}{2}}^n\right)^2-\left(v_{i,j+\frac{1}{2}}^n\right)^2-\left(v_{i,j-\frac{1}{2}}^n\right)^2+\left(v_{i,j-\frac{3}{2}}^n\right)^2+2v_{i,j+\frac{1}{2}}^n\left(v_{i,j+\frac{3}{2}}^n-v_{i,j-\frac{1}{2}}^n\right)+2v_{i,j-\frac{1}{2}}^n\left(v_{i,j-\frac{3}{2}}^n-v_{i,j+\frac{1}{2}}^n\right)\right]+\frac{\mu}{\left(\Delta x\right)^2}\left(D_{i+1,j}^n-2D_{i,j}^n+D_{i-1,j}^n\right)+\frac{\mu}{\left(\Delta y\right)^2}\left(D_{i,j+1}^n-2D_{i,j}^n+D_{i,j-1}^n\right)\right\}其中
D_{i,j}^n=\frac{1}{\Delta x}\left(u_{i+\frac{1}{2},j}^n-u_{i-\frac{1}{2},j}^n\right)+\frac{1}{\Delta y}\left(v_{i,j+\frac{1}{2}}^n-v_{i,j-\frac{1}{2}}^n\right)固壁条件为：
u_{i-1/2,j}^n=v_{i-1/2,j}^n=0进而 (i-1,j)，(i-3/2,j) 等 x<=i-1/2 得点上的物理量的值均为零，代入得：
u_{i+1/2,j}^{n+1}=u_{i+1/2,j}^n-\frac{\Delta t}{4\Delta x}\left[u_{i+\frac{3}{2},j}^n\left(u_{i+\frac{3}{2},j}^n+2u_{i+\frac{1}{2},j}^n\right)-0\right]-\frac{\Delta t}{4\Delta y}\left[\left(u_{i+\frac{1}{2},j}^n+u_{i+\frac{1}{2},j+1}^n\right)\left(v_{i,j+\frac{1}{2}}^n+v_{i+1,j+\frac{1}{2}}^n\right)-\left(u_{i+\frac{1}{2},j-1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i,j-\frac{1}{2}}^n+v_{i+1,j-\frac{1}{2}}^n\right)\right]-\frac{1}{\rho}\frac{\Delta t}{\Delta x}\left(p_{I+1,J}^n-p_{I,J}^n\right)+\frac{\mu\Delta t}{\left(\Delta x\right)^2}\left(u_{i+\frac{3}{2},j}^n-2u_{i+\frac{1}{2},j}^n+0\right)+\frac{\mu\Delta t}{\left(\Delta y\right)^2}\left(u_{i+\frac{1}{2},j+1}^n-2u_{i+\frac{1}{2},j}^n+u_{i+\frac{1}{2},j-1}^n\right)
v_{i,j+1/2}^{n+1}=v_{i,j+1/2}^n-\frac{\Delta t}{4\Delta x}\left[\left(u_{i+\frac{1}{2},j+1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i+1,j+\frac{1}{2}}^n+v_{i,j+\frac{1}{2}}^n\right)-0\right]-\frac{\Delta t}{4\Delta y}\left[v_{i,j+\frac{3}{2}}^n\left(v_{i,j+\frac{3}{2}}^n+2v_{i,j+\frac{1}{2}}^n\right)-v_{i,j-\frac{1}{2}}^n\left(v_{i,j-\frac{1}{2}}^n+2v_{i,j+\frac{1}{2}}^n\right)\right]-\frac{1}{\rho}\frac{\Delta t}{\Delta y}\left(p_{I,J+1}^n-p_{I,J}^n\right)+\frac{\mu\Delta t}{\left(\Delta x\right)^2}\left(v_{i+1,j+\frac{1}{2}}^n-2v_{i,j+\frac{1}{2}}^n+0\right)+\frac{\mu\Delta t}{\left(\Delta y\right)^2}\left(v_{i,j+\frac{3}{2}}^n-2v_{i,j+\frac{1}{2}}^n+v_{i,j-\frac{1}{2}}^n\right)
\frac{1}{\left(\Delta x\right)^2}\left(p_{i+1,j}^n-2p_{i,j}^n+0\right)+\frac{1}{\left(\Delta y\right)^2}\left(p_{i,j+1}^n-2p_{i,j}^n+p_{i,j-1}^n\right)=\left(S_p\right)_{i,j}^n
\left(S_p\right)_{i,j}^n=\rho\left\{\frac{D_{i,j}^n}{\Delta t}-\frac{1}{{4\left(\Delta x\right)}^2}\left[\left(u_{i+\frac{3}{2},j}^n\right)^2-\left(u_{i+\frac{1}{2},j}^n\right)^2-0+0+2u_{i+\frac{1}{2},j}^n\left(u_{i+\frac{3}{2},j}^n-0\right)+0\right]-\frac{1}{2\Delta x\Delta y}\left[\left(u_{i+\frac{1}{2},j+1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i+1,j+\frac{1}{2}}^n-v_{i,j+\frac{1}{2}}^n\right)+0\right]-\frac{1}{{4\left(\Delta y\right)}^2}\left[\left(v_{i,j+\frac{3}{2}}^n\right)^2-\left(v_{i,j+\frac{1}{2}}^n\right)^2-\left(v_{i,j-\frac{1}{2}}^n\right)^2+\left(v_{i,j-\frac{3}{2}}^n\right)^2+2v_{i,j+\frac{1}{2}}^n\left(v_{i,j+\frac{3}{2}}^n-v_{i,j-\frac{1}{2}}^n\right)+2v_{i,j-\frac{1}{2}}^n\left(v_{i,j-\frac{3}{2}}^n-v_{i,j+\frac{1}{2}}^n\right)\right]+\frac{\mu}{\left(\Delta x\right)^2}\left(D_{i+1,j}^n-2D_{i,j}^n+D_{i-1,j}^n\right)+\frac{\mu}{\left(\Delta y\right)^2}\left(D_{i,j+1}^n-2D_{i,j}^n+D_{i,j-1}^n\right)\right\}其中
D_{i,j}^n=\frac{1}{\Delta x}\left(u_{i+\frac{1}{2},j}^n-0\right)+\frac{1}{\Delta y}\left(v_{i,j+\frac{1}{2}}^n-v_{i,j-\frac{1}{2}}^n\right)

[计算流体力学] NS 方程的速度压力法差分格式相关推荐

【计算流体力学】Python实现加权余量法求微分方程数值解比较伽辽金法（Galerkin法）、最小二乘法和矩法的求解精度分析误差随n增大的变化情况
一.简介微分方程的经典数值方法有Ritz法.加权余量法.有限元法.有限体积法等等.对于微分方程的边值问题,如果能找到与微分方程相对应的泛函,可以通过求取相应泛函的极小值,将微分方程转化为关于基函数待 ...
有限差分法求解不可压NS方程
网上关于有限差分法解NS方程的程序实现不尽完备,这里是一些补充注解现有的优秀资料理论向 [1]如何从物理意义上理解NS方程? - 知乎 [2]NS方程数值解法:投影法的简单应用 - 知乎 [3][ ...
计算流体力学简介(一)——一些基本概念
偏微分方程与常微分方程偏微分方程和常微分方程的区别主要就是体现在待求解函数是一元函数还是多元函数. 多元函数存在对不同自变量的偏导数,因此这种带有多元函数偏导数的方程就是偏微分方程.而对于只有一个自 ...
c语言温度转换csdn,NTC热敏电阻温度计算方法，Steinhart-Hart方程和B值法（转）
NTC热敏电阻随环境温度(T)升高,电阻值(R)会下降,反之,当温度(T)下降,电阻值(R)会上升.其对温度感应非常灵敏,NTC热敏电阻电路相对简单,价格低廉,组件精确,可以轻松获取项目的温度数据,因 ...
计算机模拟求解流体力学方程,基于CFD方法对圆盘空化器超空泡流动的数值模拟_计算流体力学-论文网...
论文摘要:为分析不同来流速度对超空泡形态及流体动力特性的影响,基于粘流理论和有限体积方法,对水下运动的圆盘空化器产生的超空泡流动进行了数值模拟.得到了空泡形态.阻力.阻力系数与航行体速度之间的关系.数 ...
Matlab 编程《计算流体力学基础及应用（约翰D安德森）》全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法
Matlab 编程 <计算流体力学基础及应用(约翰D安德森)> 全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法问题之全亚声速等熵喷管流动的CFD解法初始化参数迭代过 ...
流体力学-----动量方程的推导和理解---NS方程
文章目录前言推导动量方程的流动模型推导过程书中给的剪切力分析前提条件:速度的三个分量u.v.w的正增量和坐标轴一致前言可以参考之前的博客计算流体力学1-流体力学的控制方程推导动量方程的流 ...
工程流体力学笔记暂记24 (不可压缩粘性流体的运动微分方程**N-S方程**)
粘性流体牛顿第二动律分析x表面的受力情况为例当气体各流层之间有相对运动时,存在一种阻滞这种运动的力,一旦相对运动消失,这种力也不复存在,气体的这种性质称为粘性.其作用力称为粘性力,粘性力时和它的 ...
计算流体力学1-流体力学的控制方程
文章目录前言一.流动控制方程的推导及其物理含义的讨论.适用于CFD 计算的控制方程形式 1.1 引言 1.2 流动模型 1.3 实质导数 1.4 速度的散度及其物理意义 1.5 连续方程 1.5. ...
ftcs格式 matlab,计算流体力学_对流方程_有限差分法_Lax格式_迎风格式_FTCS格式（8页）-原创力文档...
有限差分法数值求解一维伯格斯方程作者:潭花林引言本文利用有限差分法计算了一维伯格斯方程的初边值问题.采用FTCS格式,并深入讨论了它的相容性.收敛性与稳定.有限差分法在计算流体力学.数值传热学中 ...

[计算流体力学] NS 方程的速度压力法差分格式

[计算流体力学] NS 方程的速度压力法差分格式相关推荐

最新文章

热门文章