【CFD理论】对流项-04-TVD

咔咔酷酷冲冲

已于 2022-08-29 10:54:23 修改

阅读量929

点赞数

分类专栏： of初阶文章标签：算法

于 2022-08-18 16:00:56 首次发布

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

版权

of初阶专栏收录该内容

16 篇文章 7 订阅

订阅专栏

【CFD理论】对流项-04

total variation and TVD schemes

UD scheme is the most stable scheme and does not give any wiggles
whereas the CD and QUICK schemes have higher-order accuracy and give rise to wiggles under certain conditions.
our goal is to find a scheme with a higher-order of accuracy without wiggles
monotinicity preserving 保证单调性
total variation

criteria for TVD schemes

在这里插入图片描述
HIGH RESOLUTION SCHEMES USING FLUX LIMITERS FOR HYPERBOLIC CONSERVATION LAWS

Sweby (1984)给出了TVD格式的充要条件

0<r<1，上界是 $\psi(r)=2r$ ，所以TVD格式 $\psi(r)<2r$
r>1，上界是 $\psi(r)=2$ ，所以TVD格式 $\psi(r)<2$

$\psi (r)是个通量限制器函数$
$\psi(r)对于二阶精度格式，必须通过点(1,1)，所以UD不是二阶精度，CD，LUD,QUICK都是二阶精度$

flux limiter functions 通量限制器

Van Leer
Van Albada
Min-Mod
SUPERBEE
Sweby
QUICK
UMIST

TVD schemes for unstructured grids

$r_e=\frac{\phi_P-\phi_W}{\phi_E-\phi_P}\Rightarrow$
$r_{i+1/2}=\frac{\phi_i-\phi_{i-1}}{\phi_{i+1}-\phi_{i}}$
通量求解器
$\phi_{i+1/2}=\phi_i+\frac{1}{2}\psi(r_{i+1/2})(\phi_{i+1}-\phi_i)$

在这里插入图片描述

implementation of TVD schemes

one-dimensional convection-diffusion equation
$\frac{d}{dx}(\rho u \phi)=\frac{d}{dx}[\Gamma\frac{d\phi}{dx}]$
discretised form
$F_e\phi_e-F_w\phi_w=D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W)$
positive x-direction u>0 and the value of $\phi_e$ and $\phi_w$ using a TVD scheme:
$\phi_e=\phi_P+\frac{1}{2}\psi(r_e)(\phi_E-\phi_P)$
$\phi_w=\phi_W+\frac{1}{2}\psi(r_w)(\phi_P-\phi_W)$
$r_e=\frac{\phi_P-\phi_W}{\phi_E-\phi_P},r_w=\frac{\phi_W-\phi_{WW}}{\phi_P-\phi_W}$
substition of $\phi_e$ and $\phi_w$ into equation discretised form
$F_e[\phi_P+\frac{1}{2}\psi(r_e)(\phi_E-\phi_P)]-F_w[\phi_W+\frac{1}{2}\psi(r_w)(\phi_P-\phi_W)]$
$=D_e(\phi_E-\phi_P)-(\phi_P-\phi_W)$
$[D_e+F_e+D_w]\phi_P=[D_w+F_w]\phi_W+D_e\phi_E-F_e[\frac{1}{2}\psi(r_e)(\phi_E-\phi_P)]+F_w[\frac{1}{2}\psi(r_w)(\phi_P-\phi_W)]$
$a_P\phi_P=a_W\phi_W+a_E\phi_E+S_u^{DC}$
$a_W=D_w+F_w$
$a_E=D_e$
$a_P=a_W+a_E+(F_e-F_w)$
$S_u^{DC}=-F_e[\frac{1}{2}\psi(r_e)(\phi_E-\phi_P)]+ F_w[\frac{1}{2}\psi(r_w)(\phi_P-\phi_W)]$
TVD scheme

TVD neighbour coefficients
$a_W=D_w+max(F_w,0)$
$a_E=D_e+max(-F_e,0)$
TVD deferred correction source term
$S_u^{DC}=-F_e\frac{1}{2}[(1-\alpha_e)\psi(r_e^-)-\alpha_e\cdot \phi(r_e^+)](\phi_E-\phi_P)+ F_w\frac{1}{2}[\alpha_w\cdot\psi_w(r_w^)-(1-\alpha_w)\psi(r_w^-)](\phi_P-\phi_W)$
$\alpha_w=1$ for $F_w>0$ and $\alpha_e=1$ for $F_e>0$
$\alpha_e=1$ for $F_w>0$ and $\alpha_e=0$ for $F_e<0$