引問題都伴隨著復(fù)雜的幾何固壁。然而對于這個(gè)問題的數(shù)值模擬遇到的最大就是彎本文基于網(wǎng)格結(jié)合了浸入邊界法(Immersedboundarymethod)、擴(kuò)散界面法(DiffuseInterfacemetho)以及曲面上的特征移動接觸線模從而提出一種可以處理包含復(fù)雜幾何邊界并帶有移動接觸線的多相流動問題的數(shù)值為了驗(yàn)證使用了液滴在圓柱上的鋪展和和多個(gè)圓柱半浸入液體之中的算例對算法進(jìn)行驗(yàn)證過液滴在圓柱上的鋪展的算例驗(yàn)證了算法能夠準(zhǔn)確的得到設(shè)定的接觸角；通過對于多個(gè)圓柱半浸入液體之中的算例驗(yàn)證了算法能夠處理復(fù)雜時(shí)對于移動接觸線的動態(tài)行為也可以進(jìn)行很好的模擬。

u 解，使用基于均勻交錯(cuò)網(wǎng)格上的步方法。其中對流項(xiàng)使用Adam-使用網(wǎng)格來簡單有效的模擬復(fù)雜幾何邊界而受到廣泛使用，本文正是使用一種改進(jìn)的浸入邊界法[6]。這個(gè)方法可以有效地保證固體壁面不會被流線同時(shí)能夠很好最后，關(guān)于移動接觸線模型的選擇，我們使用了彎曲壁面上的特征移動接觸線模型[7]。彎曲壁面上的特征移動接觸線模型是一種接觸線模型幾何形式的體現(xiàn)，它1）Theprojectwassupportedby 將一個(gè)二維液滴放置在一個(gè)固定圓圖1邊界為對稱邊圓關(guān)無量綱參數(shù)為：1Oh0005液密度比為0.001為0.025h0.01是圓們通個(gè)算例可以測試算法對于平衡接觸角計(jì)算的準(zhǔn)確性。當(dāng)初始接觸角大于平衡接觸角時(shí)生鋪展；當(dāng)初始接觸角小于平衡接觸角時(shí)發(fā)生回縮慮終平衡形狀是一個(gè)圓質(zhì)量守恒和設(shè)定的平衡接觸角就可以求出平衡時(shí)液滴的理論預(yù)測結(jié)果。圖2給出平衡接觸角設(shè)置為30°、60°、90°和120的界面形狀結(jié)果與理論預(yù)測結(jié)果的對比。從結(jié)果中可以明顯看出，我們的數(shù)值結(jié)果與理論預(yù)測定性上對比非常好。為了得到定量的對們需要測量出數(shù)值結(jié)們使用對于數(shù)值結(jié)果中界面形狀擬合最好的圓來進(jìn)行計(jì)圓結(jié)合圓柱的半徑可以得到數(shù)值結(jié)平衡接 平衡接 角為了驗(yàn)證本文算法對于復(fù)雜幾何邊界的適應(yīng)性以及對移動接觸線的動態(tài)行為的描慮維圓每個(gè)圓柱的平衡接觸角在圖中相應(yīng)位置已經(jīng)關(guān)無量綱參數(shù)為：1Oh0005為.001液粘性比為0.025h.01是圓 基于網(wǎng)格，提出一種可以處理包含復(fù)雜幾何邊界并帶有移動接觸線的多相流動HuDL,BushJWM,Thehydrodynamicsofwater otion,Nature424(2003)AristoffJM,BushJWM,Waterentryofsmallhydrophobicspheres,J.FluidMech.619(2009)ThoreauV,MalkiB,etal,Physico-chemicalanddynamicstudyofoil-dropremovalfrombareandcoatedstainless-steelsurfaces,J.Adhes.Sci.Technol.20(2006)1819–1831DingH,SpeltPDM,ShuC,Diffuseinterfacemodelfor pressibletwo-phase?owswithlargedensityratios,J.Comput.Phys.226(2007)KimJ,Phase?eldcomputationsforternary?uid?ows,Comput.MethodsAppl.Mech.Eng.196(2007)RenWW,ShuC,WuJ,etal.Boundarycondition-endimmersedboundarymethodforthermalflowproblemswithDirichlettemperatureconditionanditsapplications.Computers&Fluids.57(2012)40-51LiuHR,DingH,Adiffuse-interfaceimmersed-boundarymethodfortwo-dimensionalsimulationof?owswithmovingcontactlinesoncurvedsubstrates,J.Comput.Phys.294(2015)484–502ASIMULATIONOFMULTI-PHASEFLOWWITHCOMPLEXGEOMETRYLIU DINGUniversityofScienceandTechnologyof,DepartmentofEngineeringMechanics，Hefei:Weproposeanapproachtosimulatemulti-phase?owswithmovingcontactlines(MCLs)onthecompexgeometryonaCartesianmesh.Theapproachcombinesanimmersedboundarymethodwith ponentdiffuseinterfacemodelandacharacteristicMCLmodel.Theimmersedboundarymethodisabletoaccurayen theno-slipboundaryconditionatthesolidsurfaceandeffectivelyrepresentthecomplexgeometry.Ontheotherhand,usingthe ponentdiffuseinterfacemodelsimulatetheliquid,gasandsolidphases.ThecharacteristicMCLmodelnotonlyallowsthecontactlinestomoveonthecurvedboundaries,butensurestheprescribedcontactangle.Weexaminetheperformanceoftheapproachthroughthenumericalsimulationsofadropspreadingonacirclec