欧拉网格中基于VOF方法的表面张力项处理方法.pdfVIP

第30卷第3期 北京理工大学学报 V01．30No．3 Transactionsof Instituteof Mar．2010 2010年3月 Beijing Technology 欧拉网格中基于VOF方法的表面张力项处理方法 姚雯， 冯其京， 郝鹏程， 洪滔 (北京应用物理与计算数学研究所，北京 100088) 摘 要：为了考察表面张力对液滴破碎的抑制作用，采用Brackbill的CSF模型，研究基于VOF界面处理方法的表 面张力项离散模式．将其用于欧拉网格中，模拟了激波和拉伸流场中液滴变形与破碎过程．算例结果表明， Bussmann模型能够很好地解决VOF函数的表面张力离散问题． 关键词：表面张力；CSF模型；VOF方法；欧拉网格；液滴破碎 241 中图分类号：0351；O 文献标志码：A 文章编号：1001—0645(2010)03—0361—04 ofSurfaceTensionSolverBasedon Investigation VoFMethodinEulerMeshes YAO FENG HAo HONGTao Wen， Qi—jing， Peng—cheng， (Instituteof and 100088，China) AppliedPhysicsComputationalMathematics，Beijing ordertoattestthatsurfacetensioncanembarrassdeformationand Abstract：In drop breakup． item CSFmodelisusedto surfacetensionbasedonVOFmethodinEuler disperse meshes．Drop deformationand aresimulatedduetoshockwaveandanextensionalflow．Theresult breakup indicatedthatBussmannmodelcouldsolve ofsurfacetension． dispersionproblem words：surfacetension；CSFmodel；VOFmethod；Euler Key meshes；dropbreakup We数和0矗数对液滴变形破碎起主要作用[1]． VOF界面是非连续斜线段，但能保证质量守 低0^下(Oh0．1)，根据We不同可分为袋型破恒．VOF函数只在界面存在梯度，二阶曲率项不能 碎、剥落破碎、爆炸破碎． 中心差分，张力项处理较困难． 表面层分子受两种分子作用．由于两种分子性 质不同，表面层分子受力的球对称性遭到破坏，因而 围扩至6×6网格，得到基本准确的张力值．作者还 产生合力即表面张力． 考察了表面张力对空气中液滴变形破碎的影 若界面两侧分子性质近似，表面张力可忽略；差 响[4剖，并与实验进行了比较． 别较大时就必须考虑张力影响．由于张力只在界面