【CFD理论】对流项-03

最新推荐文章于 2023-07-13 17:20:27 发布

咔咔酷酷冲冲

最新推荐文章于 2023-07-13 17:20:27 发布

阅读量584

点赞数

分类专栏： of初阶文章标签： python 开发语言

本文链接：https://blog.csdn.net/zq93538196/article/details/126337655

版权

of初阶专栏收录该内容

16 篇文章 7 订阅

订阅专栏

【CFD理论】对流项-03

first order upwind

Divergence schemes
一维稳态定常（稳态）对流扩散方程
$F_e\phi_e-F_w\phi_w=D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W)$
对流项（面心）=梯度（体心）

high resolution schemes

QUICK

Quadratic Upstream Interpolation for Convective Kinematics
Second order
Unbounded
Leonard
在网格规则的情况下是三阶

$\phi_{face}=\frac{6}{8}\phi_{i-1}+\frac{3}{8}\phi_i-\frac{1}{8}\phi_{i-2}$
$\phi_e=\frac{6}{8}\phi_{P}+\frac{3}{8}\phi_E-\frac{1}{8}\phi_{W}$

一维定常对流扩散方程

$F_e\phi_e-F_w\phi_w=D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W)$
$F_e(\frac{6}{8}\phi_{P}+\frac{3}{8}\phi_E-\frac{1}{8}\phi_{W})-F_w(\frac{6}{8}\phi_{P}+\frac{3}{8}\phi_W-\frac{1}{8}\phi_{WW})=D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W)$
$\Rightarrow$
$[D_w-\frac{3}{8}F_w+D_e+\frac{6}{8}F_e]\phi_P=[D_w+\frac{6}{8}F_w+\frac{1}{8}F_e]\phi_W+[D_e-\frac{3}{8}F_e]\phi_E-\frac{1}{8}F_w\phi_{WW}$

非对称矩阵

对流项存在
采样点3个
不稳定容易发散

example

在这里插入图片描述

TVD schemes

TVD (total variation diminishing)，总变差变小
upwind scheme 无条件稳定，有界。但是会带来 false diffusion
QUICK等高阶格式，当Pe较大时，可能spurious oscillation or wiggles。当计算一些物理量时候，诱发计算不稳定，TVD格式就是用来处理该问题的。
增加人工扩散或者增加上游的权重，基于此的想法，叫做flux corrected transport (FCT)
OpenFOAM TVD参考: HIGH RESOLUTION NVD DIFFERENCING SCHEME FOR
ARBITRARILY UNSTRUCTURED MESHES

upwind differencing (UD)
$\phi_e=\phi_p$
linear upwind differencing (LUD)
$\phi_e=\phi_P+\frac{1}{2}(\phi_P-\phi_W)$
$\psi=\frac{\phi_P-\phi_W}{\phi_E-\phi_P}$
CD格式
$\phi_e=\phi_P+\frac{1}{2}(\phi_P-\phi_W)$
QUICK格式
$\phi_e=\phi_P+\frac{1}{8}[3\phi_E-2\phi_P-\phi_W]$
$\psi=(3+\frac{\phi_P-\phi_W}{\phi_E-\phi_P})\frac{1}{4}$
广义高阶格式
$\phi_e=\phi_P+\frac{1}{2}\psi(r)(\phi_E-\phi_P)$