不可壓縮流體定常流動(dòng)的動(dòng)量方程

第五章不可壓縮流體定常流動(dòng)的動(dòng)量方程是否已經(jīng)可解？HSJ第五章流動(dòng)計(jì)算第一節(jié)困難與辦法連續(xù)性方程問題一二三問題一：壓力如何與速度相關(guān)聯(lián)HSJp=10p=1問題二：壓力梯度如何表達(dá)PEWNS問題三：速度梯度如何表達(dá)交錯(cuò)網(wǎng)格法HSJ解決辦法之一：交錯(cuò)網(wǎng)格法uvpu、v、p各在不同的點(diǎn)上計(jì)算第二節(jié)動(dòng)量方程HSJ第二節(jié)動(dòng)量方程令：即得離散的動(dòng)量方程動(dòng)量方程-1將壓力項(xiàng)單獨(dú)列出，便于計(jì)算：HSJSIMPLE算法Simi-ImplicitMethdforPressure-LinkedEquations先假定壓力場(chǎng)求速度場(chǎng)，再用連續(xù)性方程修正壓力HSJ第三節(jié)SIMPLE算法一、基本思路二、推導(dǎo)過程Simple-1HSJ將不精確的值表示為：u=u*+uv=v*+v將修正關(guān)系表示為：對(duì)修正值也有：enPwsENWS連續(xù)性方程離散化化簡(jiǎn)HSJ連續(xù)性方程的離散化（定常不可壓）壓力修正方程整理，得壓力修正方程HSJ注意b的表達(dá)式！注意事項(xiàng)對(duì)u、v和p進(jìn)行計(jì)算時(shí)，其單元體是不同的，因此方程中各系數(shù)的表達(dá)式中的x、y和x、y可能是不同的。有些值需要進(jìn)行插值計(jì)算。實(shí)施步驟三、SIMPLE算法的實(shí)施步驟1、估計(jì)壓力場(chǎng)p*；2、求解動(dòng)量方程，得u*、v*、w*；3、求解壓力修正方程得p'；4、計(jì)算p＝p*＋p'；5、修正速度u＝u*＋u'、v＝v*＋v'、w＝w*＋w'；6、計(jì)算其它物理量（如溫度、濃度、湍流度等）；7、將p

視為新的p*，返回第二步。HSJ邊界條件四、邊界條件1、動(dòng)量方程的邊界條件HSJ速度邊界條件固壁邊界速度進(jìn)口邊界出口邊界單向坐標(biāo)法速度已知壓力修正邊界條件壓力邊界條件壓力修正方程中，aE=0HSJ2、壓力修正方程的邊界條件p’=0速度（流量）邊界條件ueue壓力的相對(duì)性五、壓力的相對(duì)性絕對(duì)壓力不確定壓力差計(jì)算的相對(duì)誤差增大當(dāng)邊界條件不包含壓力的時(shí)候，給定一點(diǎn)的壓力為零HSJ計(jì)算實(shí)例六、計(jì)算實(shí)例：方腔驅(qū)動(dòng)問題HSJu＝1v＝0固體邊界1×1流體域＝1.0＝1／Re交錯(cuò)網(wǎng)格的實(shí)施交錯(cuò)網(wǎng)格的實(shí)施HSJv2u1v1u2速度分布速度分布HSJ壓力分布?jí)毫Ψ植糎SJ壓力三維壓力分布HSJ速度star速度分布：Star－CD（單元）HSJ速度