【CFD理论】对流项-06高分格式

high resolution schemes

• high order (HO)
• high resolution (HR)
• 对流有界评价标准 convection boundedness criterion
• Normalized Variable Formulation, NVF
• 总变差变小，total variation diminishing
• 两种方法，downwind weighting factor-DWF, normalized weighting factor-NWF

$对应OPENFOAM/numerics/div/NVD$

the normalized variable formulation (NVF)

• Leonard最早提出
• Gaskell, Lau 1988简化提出 对流有界评价标准
• the normalized variable diagram-NVD可以用来分析高阶和高分辨率格式

NVF
$$\begin{gathered} \tilde{\varphi }=\frac{\varphi -{\varphi }_U}{{\varphi }_D-{\varphi }_U} \\ {\tilde{\varphi }}_f=f({\varphi }_U,{\varphi }_C,{\varphi }_D) \end{gathered}$$
TVD
$r=\frac{{\varphi }_P-{\varphi }_U}{{\varphi }_D-{\varphi }_P}$

分情况讨论$\widetilde{{\varphi }_C}$

• a. & b. 表示C出现极值
• c & d 表示C出现梯度突变
• e 表示单调变化

HO schemes into linear relations between ${\tilde{\varphi }}_f$ and ${\tilde{\varphi }}_C$

$$\begin{gathered} {\tilde{\varphi }}_U=0 \\ {\tilde{\varphi }}_D=1 \end{gathered}$$

所有基于三节点以内插值的高阶格式，${\tilde{\varphi }}_f$可以表示${\tilde{\varphi }}_C$的线性函数：
${\tilde{\varphi }}_f=\iota {\tilde{\varphi }}_C+k$
因此，可以把上述关系表现为NVD图

the convection boundedness criterion (CBC)

$$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& f({\tilde{\varphi }}_C)& ,& continuous,\\ & f({\tilde{\varphi }}_C)=1& ,& if{\tilde{\varphi }}_C=1,\\ & f({\tilde{\varphi }}_C)with& ,& if0<{\tilde{\varphi }}_C<1,\\ & f({\tilde{\varphi }}_C)=0& ,& if{\tilde{\varphi }}_C=0,\\ & f({\tilde{\varphi }}_C)={\tilde{\varphi }}_C& ,& if{\tilde{\varphi }}_C<0or{\tilde{\varphi }}_C>1,\end{array}\end{array}$$

• 高阶格式，对于${\tilde{\varphi }}_C$必须经过(0,0),(1,1)两个点
• 在${\tilde{\varphi }}_C<0,{\tilde{\varphi }}_C>1$非单调区间，采用迎风格式

high resolution (hr) schemes

• MINMOD
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{3}{2}{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{2}\\ & \frac{1}{2}{\tilde{\varphi }}_C+\frac{1}{2}& ,& \frac{1}{2}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$

divSchemes
{
    default         none;
    div(phi,U)      Gauss Minmod;
}

Minmod divergence scheme

• Bounded CD
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{1}{2}{\tilde{\varphi }}_C+\frac{1}{2}& ,& 0\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$

• OSHER
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{3}{2}{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{2}{3}\\ & 1& ,& \frac{2}{3}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$

• SMART
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{3}{4}{\tilde{\varphi }}_C+\frac{3}{8}& ,& 0\le {\tilde{\varphi }}_C\le \frac{5}{6}\\ & 1& ,& \frac{5}{6}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$

divSchemes
{
    default         none;
    div(phi,U)      Gauss QUICK;
}

SMART(QUICK)

• STOIC
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{1}{2}{\tilde{\varphi }}_C+\frac{1}{2}& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{2}\\ & \frac{3}{4}{\tilde{\varphi }}_C+\frac{3}{8}& ,& \frac{1}{2}\le {\tilde{\varphi }}_C\le \frac{5}{6}\\ & 1& ,& \frac{5}{6}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$
• MUSCL
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& 2{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{4}\\ & {\tilde{\varphi }}_C+\frac{3}{4}& ,& \frac{1}{4}\le {\tilde{\varphi }}_C\le \frac{3}{4}\\ & 1& ,& \frac{3}{4}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$

divSchemes
{
    default         none;
    div(phi,U)      Gauss MUSCL;
}

MUSCL

• SUPERBEE
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& \frac{1}{2}+\frac{1}{2}{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{2}\\ & \frac{3}{2}{\tilde{\varphi }}_C& ,& \frac{1}{2}\le {\tilde{\varphi }}_C\le \frac{2}{3}\\ & 1& ,& \frac{2}{3}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$
• Modified SMART
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& 3{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{6}\\ & \frac{3}{4}{\tilde{\varphi }}_C+\frac{3}{8}& ,& \frac{1}{6}\le {\tilde{\varphi }}_C\le \frac{7}{10}\\ & \frac{1}{3}{\tilde{\varphi }}_C+\frac{2}{3}& ,& \frac{7}{10}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$
• Modified STOIC
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& 3{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{5}\\ & \frac{1}{2}{\tilde{\varphi }}_C+\frac{1}{2}& ,& \frac{1}{5}\le {\tilde{\varphi }}_C\le \frac{1}{2}\\ & \frac{3}{4}{\tilde{\varphi }}_C+\frac{3}{8}& ,& \frac{1}{2}\le {\tilde{\varphi }}_C\le \frac{7}{10}\\ & \frac{1}{3}{\tilde{\varphi }}_C+\frac{2}{3}& ,& \frac{7}{10}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$
• Modified SUPERBEE
  $$\begin{array}{c}{\tilde{\varphi }}_f=\{\begin{array}{rlrl}& 2{\tilde{\varphi }}_C& ,& 0\le {\tilde{\varphi }}_C\le \frac{1}{3}\\ & \frac{1}{2}{\tilde{\varphi }}_C+\frac{1}{2}& ,& \frac{1}{3}\le {\tilde{\varphi }}_C\le \frac{1}{2}\\ & \frac{3}{2}{\tilde{\varphi }}_C& ,& \frac{1}{2}\le {\tilde{\varphi }}_C\le \frac{2}{3}\\ & 1& ,& \frac{2}{3}\le {\tilde{\varphi }}_C\le 1\\ & {\tilde{\varphi }}_C& ,& elsewhere\end{array}\end{array}$$
  

the TVD framework

$${\varphi }_f={\varphi }_C+\frac{1}{2}\psi (r_f)({\varphi }_D-{\varphi }_C),with{r}_f=\frac{{\varphi }_C-{\varphi }_U}{{\varphi }_D-{\varphi }_C}$$
$$\begin{array}{c}\{\begin{array}{rlrl}& Upwind& ,& \psi (r_f)=0\\ & Downwind& ,& \psi (r_f)=2\\ & FROMM& ,& \psi (r_f)=\frac{1+r_f}{2}\\ & SOU& ,& \psi (r_f)=r_f\\ & CD& ,& \psi (r_f)=1\\ & QUICK& ,& \psi (r_f)=\frac{3+r_f}{4}\end{array}\end{array}$$

The NVF-TVD relation

$$\begin{gathered} r_f=\frac{{\varphi }_C-{\varphi }_U}{{\varphi }_D-{\varphi }_C}=\frac{({\varphi }_C-{\varphi }_U)/({\varphi }_D-{\varphi }_U)}{({\varphi }_D+{\varphi }_U+{\varphi }_U-{\varphi }_C)/({\varphi }_D-{\varphi }_U)}=\frac{{\tilde{\varphi }}_C}{1-{\tilde{\varphi }}_C} \\ \Rightarrow {\tilde{\varphi }}_C=\frac{r_f}{1+r_f} \end{gathered}$$

example11.1

• 控制体面积
  $\Delta A_i=2\times 10^{-2}{m}^2$
• 节点之间的距离
  $\Delta \xi =2/\sqrt{3}\times 10^{-2}m$
• node1
  face AB:$k\times \frac{(T_2-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_2-T_P)$
  face BC:$k\times \frac{(T_2-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_8-T_P)$
  face AC:$k\times \frac{(T_2-T_P)}{1/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=2k\sqrt{3}(T_{AC}-T_P)$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_2-T_P)+k\sqrt{3}(T_8-T_P)+2k\sqrt{3}(T_{AC}-T_P)=0$
  $4T_1=T_2+T_8+2\times T_{AC}$
  $4T_1=T_2+T_8+1000$
• node2
  face AB:$k\times \frac{(T_1-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_1-T_P)$
  face BD:$k\times \frac{(T_3-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_3-T_P)$
  face AD:$0$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_1-T_P)+k\sqrt{3}(T_3-T_P)+0=0$
  $2T_2=T_1+T_3$
• node3
  face BD:$k\times \frac{(T_2-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_2-T_P)$
  face ED:$k\times \frac{(T_4-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_4-T_P)$
  face BE:$k\times \frac{(T_{BE}-T_P)}{1/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=2k\sqrt{3}(T_{BE}-T_P)$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_2-T_P)+k\sqrt{3}(T_4-T_P)+2k\sqrt{3}(T_{BE}-T_P)=0$
  $4T_1=T_2+T_4+2\times T_{BE}$
  $4T_1=T_2+T_4+800$
• node4
  face ED:$k\times \frac{(T_3-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_3-T_P)$
  face EF:$k\times \frac{(T_5-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_5-T_P)$
  face DE:$0$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_3-T_P)+k\sqrt{3}(T_5-T_P)+0=0$
  $2T_4=T_5+T_3$
• node6&8
  face EH:$k\times \frac{(T_{EH}-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=2k\sqrt{3}(T_{EH}-T_P)$
  face EG:$k\times \frac{(T_5-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_5-T_P)$
  face HG:$k\times \frac{(T_7-T_P)}{1/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_7-T_P)$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_5-T_P)+k\sqrt{3}(T_7-T_P)+2k\sqrt{3}(T_{EH}-T_P)=0$
  $4T_6=T_5+T_7+2\times T_{EH}$
  $4T_6=T_5+T_7+400$
  $4T_8=T_1+T_9+400$
• node7&9
  face IG:$k\times \frac{(T_{IG}-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=2k\sqrt{3}(T_{IG}-T_P)$
  face HG:$k\times \frac{(T_6-T_P)}{2/\sqrt{3}\times 10^{-2}}\times 2\times 10^{-2}=k\sqrt{3}(T_6-T_P)$
  face IH:$0$
  ${\sum }_{allsurfaces}{D}_i(T_{nb}-T_P)=0$
  $k\sqrt{3}(T_6-T_P)+0+2k\sqrt{3}(T_{IG}-T_P)=0$
  $3T_7=T_6+T_7+2\times T_{IG}$
  $3T_7=T_6+1000$
  $3T_9=T_8+1000$
  MATLAB

clc;
clear;
Su=[1000,0,800,0,1000,400,1000,400,1000]';
T=zeros(9,1);
a=zeros(9,9);
a(1,:)=[4,-1,0,0,0,0,0,-1,0];
a(2,:)=[-1,2,-1,0,0,0,0,0,0];
a(3,:)=[0,-1,4,-1,0,0,0,0,0];
a(4,:)=[0,0,-1,2,-1,0,0,0,0];
a(5,:)=[0,0,0,-1,4,-1,0,0,0];
a(6,:)=[0,0,0,0,-1,4,-1,0,0];
a(7,:)=[0,0,0,0,0,-1,3,0,0];
a(8,:)=[-1,0,0,0,0,0,0,4,-1];
a(9,:)=[0,0,0,0,0,0,0,-1,3];
T=a\Su;