摘要

风波增长现象广泛存在于自然中.尽管通过实验研究已经提出了多种理论假设,但由于风场和波场的随机性和多尺度特性,未能通过实验获得可靠的原理解释其机理.随着计算流体力学的发展,对风波演化进行数值模拟以研究其生成机理更为便利.文章基于自主开发的直角网格两相流求解器,结合高阶数值格式和VOF(volume of fluid)界面捕捉方法,对阶段和指数增长阶段组成.为分析能量特征,计算了能量集中的主要波数,并结合波高、波斜率与波数的变化与实验相对比.在一维波数空间中统计能量分布,与Kolmogorov-ZakharovE (k)~k-15/4幂律进行对比.在二维波数和时域空间中,增长前期应用Phillips线性增长理论和共振加强后的增长模型对风波线性增长阶段进行验证;增长后期应用Miles理论对风波指数增长进行验证.揭示了风波增长不同阶段的机理特征,观察到在指数增长开始前,波浪的线性增长存在初始响应阶段、风压稳定阶段和风压受扰阶段,分别具有不同的机理特征.证明了波浪演化的初始线性阶段是由湍流波动产生的随机尾流叠加占主导,而后期的指数增长阶段是由增长率和摩擦速度正相关占主导.原始的Phillips理论仅在增长初始响应阶段准确,使用共振加强修正后则在初始响应后也符合实际增长,而Miles理论仅在单色波阶段时才准确.相较于对波陡有限制的动态曲线网格求解器,应用5阶WENO (weighted essentially non-oscillation)和动量-质量一致推进方法修正的VOF方法也可用于分析风波增长问题的求解,并能准确捕捉到更多风波增长的细节.