基于Fluent的自由表面初扰动现象的数值模拟.pdfVIP

第12卷 第8期 中国水运 VoI．12 No．8 2012年 8月 OhnaWater．Transport August 2012 侯玲 (广东海洋大学工程学院，广东湛江524088) 面的襁始扰动，并通过软件的相关功能，对自由表面流动中有初始位移的重力波动问题进行数值模拟。将数值计算 结果与小振幅波重力波理论解进行比较分析，研究数值计算过程中的一些影响因素。研究表明，理论分析的小振幅 假设，只有当波动的振幅与波长比小于一定数值时，才可以得到比较好的线性结果。此外，通过比较分析不同工况 下的计算结果，发现在数值模拟时，存在着自由表面厚度的问题，对计算的精确性有一定影响。 关键词：Fluent；自由表面；数值计算；小振幅；C语言；UDF 中图分类号：18675 文献标识码：A 一、引言 液面以下是静水压力，所以韧始速度为零，：初始压力为一个 近年来随着计算机技术的高速发展，计算流体力学软件 大气压。 四、数值计算结果 Fluent在各个领域均得到了广泛应用，Fluent软件能够比 较真实的再现流体的运动情况，因为受到用户的普遍关注。 0．01 应用Fluent软件进行数值模拟的一个比较典型的问题就是 自由表面问题。目前处理自由表面的方法有很多，如MAC 0．00．5 法，VOF法，SPH法等。相关的已有大量的文献[103]等对此占0 ^ 做了综述。文献【4】将数值模拟与理论、试验结果相比较，验 一0．005 证了Fluent计算模型的正确性。 本文借助该软件模拟流体力学中一种常见的波动现象 —0．0l ——初始扰动引起的重力波问题。通过数值模拟，将数值计 O 2 4 6 8 算结果和理论值比较的同时，通过不同工况的计算分析影响 图2 自由表面质点波动的时间历程(理论值) 计算结果准确性的因素。 为了便于比较数值计算结果的合理性，图2给出了自由 二、几何模型 液面在初值扰动作用下，通过理论计算得到的自由表面某一 图1为数值模拟自由液面初始扰动的几何模型，整个计 点在重力作用下的波动的时间历程。由于理论分析时做了一 算域为一矩形区域。 个小振幅假设，即波动振幅相对于波长为小量(A，A1)。 广’。————_‘—————————————～——。‘r——…I 通过数值计算发现，A／2的比值变化对数值计算结果有一定 影响。所以在数值模拟时分别研究了AIA的比值大小和网格 变化对数值计算结果的影响。 y 1．小振幅假设的影响 ● I一一一X 为了确定,4／2的比值对计算结果的影响。表1中分别给 圈1 自由衷面运动的二维模型 出了6种工况，基于不同的．4／2进行数值模拟。 三、数值方法 表1计算工况(基于不同的彳，2) 本文采用计算流体软件ANSYS—FLUENT进