SIMPLEC 二维、结构网格单相流SIMPLE 算法的C 语言实现的说明(17)

的原因可以参考《传热与流体流动的数值计算》５．１￣５．２节。　

Sc与SP分别是源项的常数项与线性项，参阅《传热与流体流动的数值计算》４．２－５节。　

(Je Feφp)可以按照如下的形式进行模拟：　

Je Feφp=aE(φP φE)

Je Fwφp=aW(φW φP)

Je Fnφp=aN(φP φN)

Je Fsφp=aS(φS φP)

aE=DeA(|PE|) § Fe,0¨,

aW=DwA(|P|) §Fw,0¨W

aN=DnA(|PN|) § Fn,0¨

aE=DsA(|PS|) §Fs,0¨

（传热５．５７与５．５８）搭建了通量JeDe=Dw=Ds=Dn=Γe yδxeΓw yδxwΓs xδxs　　　　　　　　　　　　　　　　　　　　　　　　　　（传热５．５７与５．５８）　Γn xδxn,Jw,Js,Jn与变量φP,φE,φW,φS,φN之间的桥梁，从而构成了一个完整的差分方程。微分方程（传热２．１５）到差分方程（传热５．５６）的转变过程也就此完成。（传热５．５７与５．５８）中的Ａ（｜Ｐ｜）代表包括中心差分、上风、混合、幂函数以及指数形式不同的差分格式。具体见（传热表５．２）。ＳＩＭＰＬＥ＿Ｃ采用了幂函数格式。它的具体实现过程是ｄｉｆｌｏｗ（　）函数。　

应该注意的是，ＳＩＭＰＬＥ算法中并没有把动量方程的压力项计入源项。动量方程（速度ｕ与速度ｖ）的离散方程最后的形式是通用变量的离散方程（传热５．５６）上加入压力梯度项。　

aeue=∑anbunb+b+(PP PE)Ae

aeue=∑anbunb+b+(PP PE)Ae　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　（传热　６．６）　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　（传热　６．７）　式子（传热　６．６）与（传热　６．７）分别与（传热图６．７）与（传热图６．８）对应。（传热　６．６）与（传热　６．７）更详细的推导过程请参阅《传热与流体流动的数值计算》　６．４节。　

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库SIMPLEC 二维、结构网格单相流SIMPLE 算法的C 语言实现的说明(17)在线全文阅读。