某些不可压缩分层流的计算.pdfVIP

某些不可压缩分层流的计算 李谊乐刘应中张怀新翁剑锋 (上海交通大学，上海200030) 摘要本文数值求解NS方程，计算_r分层流体绕过方拄和在箝内晃荡时的流态、密度和压 力分布，用图形显示了所撙结果。形象地表刀i肘0们和常密度情形的异J叫． 关键词分层流，@s方程数值解≥≯桂绕流．流体晃荡 1前言 大气和海洋一般是分层的，在环境流体力兰中分层流也di据着十分重要的地位，冈 此多年来一直受到多方面的广+泛关注。}抟艇在水r航行，也受到海洋分层的影响；特别 是潜艇的隐身性能更与海洋的分层密切相关。近年来人们也注意到了海水分层对水面船 舶和海洋结构物运动的影响，只是多局限丁间断分层的特殊情形。连续分屡对物体运动 的影响则鲜有提及。本文!】!IJ尝试Hj数值求解NS方样的方法，米搦示物体在分层流中的绕 流和层结流体在容器中振荡的特点，以_l{f』获得一些感性认识，为避一步研究这类问题打 下基础。 计算程序是以原有的二维SOLA．VOF程序[1]为基础．经过修改以计入密度的变 化。SOLA-VOF是个老程序．数值格式的精度不太高，但描述一些人幅三薹动还是成功的【2】。 因此除了对原有程序的适当修改外，在算例的选抨上也充分注意到这一点。所选算例一 个是方拄绕流，一个是层结流体在容器内的人幅晃荡。经过对选定算例f}{】计算，看到了 一些耐人寻味的现象。 2控制方程与数值格式 计算中用到的主要控制方程有： 不可压缩流体3_矿20害mVp．o (1) 连续性方程 V·厅=0 (2) 动量方程j 丝+厅．v厅：户一里+旦甲2厅 (3) c3t p P 系统是封闭的，未知数与方程的数日相等。 在SOLA-VOF[1】中，动量方程(3)离散为 2-287． ∥“一∥ 。 v∥=一一芝+旦v：厅” —石一州‘ P p 实际计算时拆成 业+i一．v厅一：F‘“+EV：i”竺=竺+面一．V厅一=‘“+面” (4) At P 和 厅…：磊～Atv@“1 (5) p 其中 羊Ⅱ (6) p一“：p一+6p一 磊：i。一_AtVP” p 将式15)代入离散的连续性方程，得到压力增量的Poissoo方程。 方程，得到压力增量6口。然后可以算出一F--时刻的速度和压力，进入下一时间步。 j}{中心著分。压力Poisson方程刑共轭余量法求解。 密度的传输方程(1)采心守恒形式，离散中用算子分裂方法；时间上为一阶向前显 式格式，形式是 型二堕+塑＆!二塑盏!：o (7) At Ax， 型二型+鲨坠￡鲨坠垒：o (8) At 姆i 以式(7)为例，假定左边界上岳手{≤争霸／，通过左边界的p流量可写作 中i为限制因子．西，取0时，(9)为一阶精度的迎风格式：中。等于1时，(9)为二阶 在斜率变化较大处产生振荡，出现数值解的上溢或下溢。TVD格式中o，取为两边P的斜