-15-2014中文第15章有限元方法在流體力學(xué)中的應(yīng)用

Chapter15有限元方法在流體力學(xué)中的應(yīng)用ApplicationsofFEMinFluidMechanics(FundamentalsoffiniteelementanalysisbyDavidV.Hutton)1引言不可壓縮流體(incompressibleflow):密度不變可壓縮流體(compressibleflow)Newton’slawofviscosityμ為絕對粘度(absoluteviscosity).是流體的基本材料參數(shù),與其抗剪應(yīng)力性能直接相關(guān)。如果流體的粘度很小，則可忽略其剪應(yīng)力，理想化為無粘性流體。2不可壓縮流體的控制方程質(zhì)量守恒連續(xù)性方程u,v,w,是速度在x,y,z方向的分量質(zhì)量守恒要求一個體積里面的質(zhì)量變化率等于該體積的質(zhì)量凈流入速度。體積dV里的總質(zhì)量為ρdV,注意到dV為常數(shù),所以有：從x,y,z三個方向流入質(zhì)量造成的控制體積中的質(zhì)量變化率分別表示為：所以質(zhì)量的變化率為：注意到dV=dxdydz,于是得到連續(xù)性方程為：對于穩(wěn)態(tài)的不可壓縮流體(steadyflowofanincompressiblefluid),密度與時間和空間坐標(biāo)無關(guān)，于是有旋流和無旋流(RotationalandIrrotationalFlow)把流體流動分為：有旋流動(rotational)

---平動和轉(zhuǎn)動混合無旋流動(irrotational)---僅有平動。流體微元體積沒有純轉(zhuǎn)動。如果下列三式不能同時滿足，則為有旋流動:我們目前僅考慮無旋流動。二維流的流函數(shù)(steamfunction)對于2D,穩(wěn)態(tài),不可壓縮,無旋流動連續(xù)性方程為無旋條件退化為如果引入

ψ(x,y)(流函數(shù))，則連續(xù)性條件自動滿足:無旋條件變?yōu)?Laplace’sequation)流函數(shù)的物理意義流線(streamlines):x-y平面內(nèi)的曲線，其上的流函數(shù)為常數(shù)。流線上的任一點的切向量可以表示為

nt

=dxi+dyj，該點處的流體速度為V

=ui+vj.則

V×nt=(?vdx+udy)k=0因為：兩個非零向量的叉積為零意味著這兩個向量平行，所以：在流線上任一點，流體速度正切于流線。有限元列式具有M個節(jié)點的有限單元內(nèi)的流函數(shù)可以表示為:利用Galerkin方法,單元的殘差方程為：or應(yīng)用Green-Gauss定理，上式變?yōu)槠渲?/p>

S

為單元邊界(nx

,ny

)為邊界單位外法線向量。代入流函數(shù)表達式，有：寫為矩陣形式二維流動的流速勢函數(shù)(VelocityPotentialFunction)假定存在流速勢函數(shù)φ(x,y)使得則無旋條件能自動滿足:連續(xù)性條件變?yōu)?(again,weobtainLaplace’sequation)在流速勢函數(shù)為常數(shù)的曲線上，有所以可見，速度向量垂直于流速勢為常數(shù)的曲線。于是，流線和流速勢等值線形成了正交網(wǎng)格，稱為流網(wǎng)(flownet).有限元列式與流函數(shù)的情形類似：單元剛度陣與流函數(shù)法的相同，但是節(jié)點力完全不同。不可壓縮粘性流假定:可以看作二維問題不涉及熱3.密度和粘度為常數(shù)4.穩(wěn)態(tài)(不隨時間變化)表示動量守恒的Navier-Stokes方程為:u,v=x,y方向速度分量

ρ=密度p=壓力

μ=絕對粘度FBx

,FBy

=x,

y方向單位體積的體力斯托克斯流(StokesFlow)Stokesflow(orcreepingflow)：流體流動的速度很小，慣性項比粘性項小很多，可以忽略。高粘度流體，如融化的高分子材料。動量方程變?yōu)樯鲜胶瓦B續(xù)性方程一起構(gòu)成了三個未知數(shù)u(x,y),v(x,y),andp(x,y)的三個方