波浪作用下水平柱体水动力系数研究.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 波浪作用下水平柱体水动力系数研究∗ 李勇 林缅 ( 中国科学院力学研究所，北京，100190) 摘要： 本研究通过建立数值波浪水槽，对波浪作用下的水平柱体的受力进行了计算， 并利用最小二乘法计算得到水动力系数。通过与实验数据的比较验证了模型及其计算方法的 正确性，探讨了线性波和非线性波对CD 以及CM 的影响。 关键词： 水平柱体；波浪力；水动力系数；数值计算 1 引言 水平柱体是海洋油气资源开发过程中经常被用到的构件之一。柱体所受的水动力是工程 设计中必须考虑的主要问题之一。近年来随着海洋油气开采的迅速发展，人们越来越重视对 这类构件承载力的研究。 关于水下柱体水动力系数的实验研究非常多。Bearman[1] 和Tanaka[2]在U 型振荡流水槽 中对长（正）方形、圆形和菱形等截面形状柱体的试验发现，不同形状柱体的CM 随着KC 的 变化表现出不同的变化趋势。然而Ikeda 等人[3]在波浪水槽中也做了类似的实验，发现KC 比 较小时波浪槽中测得的CD 和CM 值比U 型水槽中测得的值要低。这说明应该重新认识波浪的 垂向特征，特别是在水表面附近。本研究将通过数值实验的方法研究波浪作用下水平圆柱体 的水动力系数CD 以及CM 值。 2 数学模型 控制方程为不可压缩流体的N-S 方程和连续性方程： 抖 r 抖t + xj (r uj )= 0 (1) 抖 抖p 殒 ¶ uj u i 鳄 ( ) 犏 ç ÷ (2) r u + (r u u )= - + m ç + ÷ + r g i j i e ÷ i 抖 抖 犏 ç 抖 ÷ t x x x 桫x x j i j 臌 i j 其中， xj (j = 1, 2)表示坐标分量，uj (j = 1, 2)表示时均流速分量，p 为压强，r 为水的 ∗国家科技部“863”项目(2006AA09Z301)和国家自然科学基金项目资助. - 678 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 密度，g 为重力加速度，m = m + m 为有效黏性系数，m 为水的分子黏性系数，m 为紊动 e f f 黏性系数,根据标准的k -e 模型给定，即m = C r k2 e ，k 是紊动动能，e 是紊动耗散率， f m 分别满足以下方程： 抖 抖轾 m k 骣 f ÷ ç ( )