液面斜撞击的数值模拟.pdfVIP

　第 22 卷 　第 1 期 爆 炸 与 冲 击 Vol . 22 , No . 1 　 　2002 年 1 月 EXPLOSION AND SHOCK WAVES J an . ,2002 　 ( ) 文章编号 : 2002 Bingham 流体与固壁 、液面斜撞击的数值模拟 1 1 2 2 徐胜利 , 岳朋涛 , 刘大有 , 段新平 ( 1. 中国科学技术大学力学和机械工程系 ,安徽 合肥 　230026 ; 2 . 中国科学院力学研究所 ,北京 　100080) 　　摘要 : 在交错网格中用 MAC 方法求解二维不可压 NS 方程 ,对圆截面液柱与固壁 、液面斜撞击问题进行 了数值模拟 ,得到了自由面随时间演化的图像 。主要考察了: ( 1) 固壁和液面 ; (2) 牛顿 、非牛顿流体 ; (3) 碰 θ( ) 撞入射角 ; 4 Bingham 流体近似本构式中参数 q0 、K 对计算结果的影响 。结果表明:液柱与液面碰撞形成 的自由面更复杂 。碰撞初期 ,Bingham 流体和水的自由面相似 ;但碰撞后期 ,Bingham 流体的自由面相对简单 。 θ 对液柱与液面碰撞 自由面的影响较大 。本文条件下 ,当 K 20 时 , K 对 自由面的影响不大 ;当 q0 增大时 , 自 由面变得相对简单 。 　　关键词 : 非牛顿流体 ; 斜撞击 ; MAC 方法 ; 数值模拟 　　中图分类号 : O353 . 4 　　　文献标识码 : A 1 　引　言 　　生活中的许多现象 ,如雨水降落 、农作物的施肥和药物喷洒等 ,都涉及到液滴 、液柱与固壁 、液面的 [ 1] 撞击问题 ,因此 ,研究牛顿/ 非牛顿流体与固壁或液面的斜撞击是有意义的。F. H. Harlow 等 对圆柱形 水柱和液面正撞击过程作过数值模拟 。研究碰撞问题时 , 自由面的模拟十分重要 , 目前有 SL IC[2 ] 、Level set [3 ] 、VOF[4 ] 等方法 。分析表明:MAC 方法是模拟复杂自由面 、特别是大变形自由面的好方法 。因此 ,在 水动力学和爆炸力学研究中得到了广泛的应用 。MAC 方法由F. H. Harlow 和J . E. Welch 提出[5 ] , 以当地 流体速度运动的标记点分布来确定自由面位置 。该方法具有简单直观 、程序编制容易等优点 ,后来推广 到隐式 ,并用于复杂区域自由面的模拟[6 ] 。由于 MAC 方法不能精确地确定自由面位置 ,因此 ,在 自由面 上也不能准确地提出应力边界条件 。当用 MAC 方法模拟大剪切变形时 ,若网格中标记点的数密度较 小 ,流体中会出现虚假的空格子 ,但可通过增大标记点的数密度来抑制 。 　　本文的研究目的是 :选择非牛顿流体 ———Bingham 流体作为介质 ,利用数值模拟研究液柱和固壁 、液 面的斜碰撞问题 ,并与水的结果作对比。通过了解 自由面随时间的演化 ,探索斜撞击过程中所伴随的物 理现象 。主要内容包括 : ( 1) 液柱与固壁间的斜撞击 ; (2) 液柱与液面间的斜撞击 ; (3) 近似本构式中主 要参数对斜撞击的影响 。 2 　控制方程和边界条件 2. 1 　Bingham 流体本构关系 　　Bingham 流体存在屈服应力 ,其本构式具有如下形式[7 ]