伴随着对海洋环境保护日益重视和大规模海洋开发，各种无人海洋机器人如远程控制机器人(Remoted Operation Vehicle或ROV)、自治水下机器人(Automated Underwater Vehicle或AUV)、水面无人艇(Unmanned Surface Vessel或USV)等在海洋科学考察、海洋工程、海水养殖、军事等领域得到了广泛的应用。从能源供给角度来看以上传统的机动观测装置的能量供给大大限制了它们的长时间、大范围观测的能力。美国Liquid Robot公司开发了波浪能驱动的机动浮标产品Wave Glider，它利用波浪能垂荡运动的能量提供运动能量，该装置已经在海洋动物学、溢油泄漏检测、海平面升高等多领域得到应用，取得较好的效果[1]。

波浪能滑翔机(见图1)由水面浮体、水下推进机构及连接两者的缆索组成，它的优点是不需要外界能量的供给，只需要从波浪的上下运动中获取能量来驱动探测机构的前进(见图2)，因此该类型的机构能够长期、大范围、经济灵活、可靠自主的进行海洋观测任务，在民用及军用领域都具有广泛的应用前景[2]。它的前进动力主要由翼片在流体作用下的被动摆动产生。Eslam等[3]分析了在固定雷诺数下理想推进模式下的翼片在不同垂荡和纵荡运动规律下的涡泄脱落情况；Triantafyllou等[4]研究了扑翼推进性能和翼片参数的关系；Silva等[5]基于雷诺平均纳维-斯托克(Reynolds Average Navier-Stokes，即RANS)方程和动网格技术，模拟了其水动力学性能；Politis等[6]利用边界元方法对波浪垂荡驱动激励下的翼片水动力学进行了研究，分析了翼片产生的推进力和涡泄状态；B?ckmann等[7]通过控制翼片的纵荡运动获得了更好的推进性能；刘鹏、胡合文等[8-9]基于RANS方程研究了串列扑翼推进的水动力学，发现翼间涡系干扰是导致推进性能差异的主要原因。贾丽娟、郑炳欢、李小涛等[10-12]利对翼片做了水动力定常计算分析，贾丽娟计算了翼型在不同来流速度、不同水翼间隔分布情况下的升阻比值，郑炳欢计算得出水翼的旋转角度为45°时水翼产生的水平推力达到极大值，李小涛对水面母船不同来流方向和不同吃水做静水阻力计算并与经验公式对比发现趋势相同，并对固定水翼在不同来流攻角下做了定常阻力计算。

图1 波浪滑翔机Fig.1 Wave glider

图2 波浪滑翔机工作原理Fig.2 Operational mechanism

本文基于有限体积法和Fluent软件，用雷诺平均Navier-Stokes方程(RANS)，对开始运动阶段给定垂荡运动和摆动运动的推进装置翼片水动力学进行了分析和仿真，模拟了单个翼片、多翼片纵向阵列的压力分布、涡泄状况和推进动力，着重分析了翼片间距变化对推进力的影响作用。

1 波浪滑翔机推进装置翼片水动力学建模

1.1 波浪滑翔机推进装置翼片水动力学模型

波浪滑翔机是依靠推进机构随浮体上下运动对翼片产生水动力来实现运动的，波浪滑翔机在启动、稳定运转和回收等过程的动态性能的各不相同。本文对该设备的启动阶段进行了研究，启动阶段是典型的非稳态流动或湍流问题。现有的湍流数值模拟的方法有三种：直接数值模拟(Direct Numerical Simulation, DNS)，大涡模拟(Large Eddy Simulation, LES)和Reynolds平均方法(Reynolds Average Navier-Stokes, RANS)[13]。直接数值模拟可以获得湍流场的精确信息，是研究湍流机理的有效手段，需要较大的计算资源，限制了其应用范围。LES大涡模拟的计算量介于另外两者之间，相较于RANS方法，LES可以模拟更多的湍流大尺度运动，LES所用的湍流亚网格应力模型受边界的几何形状和流动类别的影响小。Reynolds平均方法只要设定合理的Reynolds应力模型，可以很容易地从RANS方程解出湍流的统计量。

对非稳态的N-S方程作时间演算，得到方程：

动量方程：

湍流运动的连续性方程：

其他变量输运方程：

式中附加应力可记为并称为雷诺应力。

在Fluent计算中采用RNGK-ε模型，使得附加项与时均值联系起来，采用SIMPLE算法进行压力速度耦合来求解运算。

1.2 边界条件及参数设置

计算中应用两方程的RNGK-ε模型，采用非平衡壁面函数法进行壁面修正，使用二维非结构划分整个流域，并采用动网格技术来保证翼片运动过程中的网格质量，翼片网格初始状态如图3所示。

对波浪滑翔机启动阶段进行分析，设定来流速度为0，对波浪滑翔机的工作原理和实验过程进行观察[8]，估计翼片摆动的时间占整个运动周期的1/4，翼片的摆动规律设为：