动态显示在非恒定流场中的应用.pdfVIP

动态显示在非恒定流场中的应用 包芸 刘欢 (中山大学近岸海洋研究中心，广州，510275) 摘要t动态显示数值模拟非恒定流场的计算结果，能更清晰地看到整个流动过程中的连续变化情况。便 于了解复杂的流动机制．并可以发现一些特殊的流动现象。 关键词t动态显示数值模拟非恒定流 O引言 1987年，McCormick等人根据美国国家科学基金会召开的“科学计算可视化研讨会’’ inScientific 的内容撰写的一份报告中正式提出了科学计算可视化(Visualization Computing —VisC，或称科学可视化--ScientificVisualization)这一概念【l】，此后至今的十数年间，科 学计算可视化已被成功的应用到流体力学、气象、航空等众多领域12]。概括起来，计算流体 力学可视化的主要内容有以下几点：(1)计算区域与计算网格的显示；(2)计算过程及流体 结构的显示；(3)计算结果的显示与分析，包括速度、压力、水位、温度等【3】。 根据可视化的不同方式，大致可以划分为两类：静态显示和动态显示。静态显示主要是 用于显示流场典型时刻(如涨急、落急等)的静止状态，动态显示则主要用于表示矢量场的 连续变化。由于矢量场的可视化对于理解复杂的流体机制有着十分重要的作用，因此，动态 显示在计算流体力学中的应用也越来越广泛。 l数学模型 数值模拟是动态显示的基础，针对河口岸线曲折、滩槽交错，水动力特征较为复杂 的特点，中山大学近岸海洋研究中心自主开发了二维水动力模式CCOST--2D。该模式具有 以下几个特点：(1)采用平面正交曲线网格；(2)采用平面有限体积法求解水动力方程组； (3)通过旋转矩阵克服科氏力产生的计算弱不稳定问题；(4)水位计算类似于ECOM采用 隐式技术；(5)采用冻结网格法模拟滩地出没水面的动边界技术；(6)对水点记数，只计算 求解水点的物理量。 1．1水动力模块 连续性方程 ”z鲁+音(叩，。)+啬(％uz。)+％等=。 缈=w一去[叩。●詈+鲁)+叩：(盯鼍+老)]-(仃詈+鲁) 动量方程 迎掣+毒仳。圻)+亳(啊。u㈦+啊心掣+。％(’u：专川袭一啊也，) =噌。嚏曙+期一譬赡一丢詈冲一。每鲁+亳(2～鲁。刳 LLaa¨a) +亳卜鲁。薏)+啬(钆。蚴a4,J+--D-2-丽a 24l 盟等型+素咖暖)+亳瓴。u％)+^吃掣+。u(_u。薏+u：詈一^吃，] =一g砒暖+翔一墨字f慝一面o-砺ao磊aar],盯一。参薏+素(2～鲁。刳 +音(～鲁。刳+未(～。割 +警未心刳 其中U-和Uz分别是曲线网格￡，和毫。方向上的速度，h-和hz是对应方向上的网格步长，7／ 为水位，D为总水深，D=H+77；垂向采用的是随水深变化的口坐标。将以上两公式垂向积分， 就可用于进行二维水动力数值计算求解。 1．2泥沙模块 泥沙采用挟沙力模型。悬沙输运方程为 _+弋一+弋一+删J—s)2__)+_【洲_)TbHs+罢坐+下bvHs+咧J—s‘)：_b4／e(sa_倒【a_+)sb_ ) df ck dy c墉cI％dy dy 式中， U、v为水流速度，H=h+T1为水深，h为静水深，q为潮位，g为水中含沙量，￡为悬 沙扩散系数，s，-c、0【、∞分别为水流挟沙能力、恢复饱和系数和悬沙沉降速度。 2算例一：动态场显示整个潮周期中泥沙运动 伶仃洋河口湾是珠江最大的河口湾，是珠江主要的入海口，汇集了珠江八大口门中的东 四口门的来水来沙。珠江河口岛屿众多，地形复杂，不同口门的来水来沙条件各不相同；同 时伴有大量的人类活动诸如抛石促淤、围垦、航道疏浚等影响在内的环境演变过程。因此， 伶仃洋的泥沙沉积与运动非常复杂，对伶仃洋河口湾进行水沙数值模拟一直是众多学者关注 的问题。利用CCOST--2D模型模拟了伶仃洋的水沙分布情况，模拟时段为1999年7月16 日0点到1999年7月24日0点共192小时。