INTERNAL SOLITARY WAVE LOADS EXPERIMENTS AND ITS THEORETICAL MODEL FOR A CYLINDRICAL STRUCTURE

在大型重力式密度分层水槽中, 对内孤立波与圆柱型结构的相互作用特性开展了系列实验. 基于两层流体中 内孤立波的KdV,eKdV和MCC理论, 建立了圆柱型结构内孤立波载荷的理论预报模型, 给出了该载荷理论预报模型中3类内孤立波理论的适用性条件.研究表明, 圆柱型结构内孤立波水平载荷包括水平Froude-Krylov力、附加质量力和拖曳力3个部分, 可以由Morison公式计算, 而内孤立波垂向载荷主要为垂向Froude-Krylov力, 可以由内孤立波诱导动压力计算.系列实验结果表明, 附加质量系数可以取为常数1.0, 拖曳力系数与内孤立波诱导速度场的雷诺数之间为指数函数关系, 而且基于理论预报模型的数值结果与系列实验结果吻合.     

孤立波与淹没平板相互作用的三维波面和水动力实验研究

在波浪水池中进行了孤立波作用下有限长度和有限宽度淹没平板的三维模型水池实验. 首次应用多目视觉立体重构技术测量局部三维自由表面变形, 该系统的有效测量水平范围为1.7 m$\times $1.6 m. 用4个三分力测力传感器组成水下测力系统, 在不影响波面的情况下测量孤立波对平板的作用力和力矩. 针对波浪不破碎的情况, 选择0.4 m水深和0.16 m波高的来波条件, 平板淹没深度为0.1 m. 实验结果表明, 孤立波经过淹没平板时自由面有明显的三维变形, 导致孤立波波幅的时空变化. 波幅在平板尾缘中心线处达到最大值, 并沿展向逐渐减小. 利用多目视觉立体重构系统得到的波面变化过程与浪高仪给出定点波面时间序列相互印证, 表明建立标识码波面测量方法是有效的. 孤立波对淹没平板作用的水动力载荷变化分为6个典型阶段, 并与利用波面三维重构得到的波面测量标识码并讨论. 基于多目视觉立体重构技术得到了垂向力和俯仰力矩极值点出现时的三维波面形态. 建立的多目视觉立体重构系统将为海洋工程结构物的水池物理模型实验提供新的波面测量手段.      

孤立波与淹没平板相互作用的三维波面和水动力实验研究

在波浪水池中进行了孤立波作用下有限长度和有限宽度淹没平板的三维模型水池实验.首次应用多目视觉立体重构技术测量局部三维自由表面变形,该系统的有效测量水平范围为1.7 m×1.6 m.用4个三分力测力传感器组成水下测力系统,在不影响波面的情况下测量孤立波对平板的作用力和力矩.针对波浪不破碎的情况,选择0.4 m水深和0.16 m波高的来波条件,平板淹没深度为0.1 m.实验结果表明,孤立波经过淹没平板时自由面有明显的三维变形,导致孤立波波幅的时空变化.波幅在平板尾缘中心线处达到最大值,并沿展向逐渐减小.利用多目视觉立体重构系统得到的波面变化过程与浪高仪给出定点波面时间序列相互印证,表明建立的三维波面测量方法是有效的.孤立波对淹没平板作用的水动力载荷变化分为6个典型阶段,并与利用波面三维重构得到的波面测量结果一并讨论.基于多目视觉立体重构技术得到了垂向力和俯仰力矩极值点出现时的三维波面形态.建立的多目视觉立体重构系统将为海洋工程结构物的水池物理模型实验提供新的波面测量手段.     

静水槽孤立波的实验研究

本文在静水槽中复现出单孤立波,将实验得到的孤立波的波形和波速与理论值比较,极为接近。同时,考察了孤立波的碰撞行为,它表明了孤立波的粒子性。还验证了孤立波的非线性性质。     

细观结构固体层中孤立波的传播及相互作用特性

利用直接微扰方法.确定了孤立波的放大或衰减与孤立波的初始幅度以及介质的结构参数之间的关系.然后利用线性化技术构造出一种二阶精度的稳定差分格式,并对孤立波在细观结构固体层中传播特性进行了数值模拟,特别对细观结构固体层中传播的不同幅度的孤立波的相互作用进行了详细的数值模拟,从而得到在适当条件下细观结构固体层中孤立波传播时即可以衰减、放大又可以稳定传播,且相互作用不影响这种传播特性.     

等波高双孤立波直墙爬高的数值模拟

基于RANS方程、VOF方法以及修正的Goring造波方法建立了模拟活塞式推波板运动的二维数值波浪水槽,实现了双孤立波直墙爬高的数值模拟．利用动边界技术模拟造波机推波板的运动,有效地实现了不同波峰间距双孤立波的造波方法．在验证单孤立波直墙爬高的基础上,模拟了不同相对波高、相对波峰间距的等波高双孤立波的直墙爬高过程,给出了波面、速度场及波动能量的变化规律．数值模拟结果表明：对于等波高的双孤立波,当入射波波高较大及两个波峰间距相对较小时,跟随在后孤立波的爬高放大系数小于先导孤立波的爬高放大系数;双孤立波在直墙爬高过程中,波动场的势能时间过程线呈现三峰形态,其中居中的最大势能峰值出现在第二个孤立波与经直墙反射后反向传播的第一个孤立波完全对撞的时刻．     

大气中非线性波和孤立波——有热源影响的非线性重力惯性内波

本文考虑有外界热源情况下的非线性重力惯性内波,首先我们讨论了波解存在的条件及得到了解析解的表达式。然后我们利用拟能影响函数中的根与系数之间关系,导出一个无量纲的量M,利用M把周期解存在条件转化为M>2/3,当M→2/3时,就得到了区别于KdV方程的孤立波解。最后利用M我们建立了非线性波的波速公式,波速C与振幅,特征散度,M之间的关系,当波速公式向线性化波速公式退化时,我们发现当M≥l时系统呈线性效应;当2/3相似文献   

孤立波与多个潜水障碍物相互作用的数值研究

本文以完整二维Navier-Stokes方程为控制方程,采用VOF界面跟踪技术和差分方法,数值计算孤立波与多个淹没水下物体相互作用的问题。本文对潜水物体的高度接近水深时,孤立波通过水下孤立直立方柱、两个间距较大的水下直立方柱和两个间距较小的水下直立方柱等三种情况分别进行了计算,给出了波形随时间的演化图,可以看到反射波、前传波和跟随振荡型小波列的生成。对孤立波通过水下孤立直立方柱情形的计算结果,与实验结果和势流理论结果进行了比较。     

基于势流理论的内孤立波与立管作用数值研究

内孤立波是一种发生在水面以下的在世界各个海域广泛存在的大幅波浪, 其剧烈的波面起伏所携带的巨大能量对以海洋立管为代表的海洋结构物产生严重威胁, 分析其传播演化过程的流场特征及立管在内孤立波作用下的动力响应规律对于海洋立管的设计具有重要意义. 本文基于分层流体的非线性势流理论, 采用高效率的多域边界单元法, 建立了内孤立波流场分析计算的数值模型, 可以实时获得内孤立波的流场特征. 根据获得的流场信息, 采用莫里森方程计算内孤立波对海洋立管作用的载荷分布. 将内孤立波流场非线性势流计算模型与动力学有限元模型结合来求解内孤立波作用下海洋立管的动力响应特征, 讨论了内孤立波参数、顶张力大小以及内部流体密度对立管动力响应的影响. 发现随着内孤立波波幅的增大, 海洋立管的流向位移和应力明显增大. 由于上层流体速度明显大于下层, 且在所研究问题中拖曳力远大于惯性力, 因此管道顺流向的最大位移发生在上层区域. 顶张力通过改变几何刚度阵的值进而对立管的响应产生明显影响. 对于弱约束立管, 内部流体的密度对管道的流向位移影响较小.      

移动载荷作用下浮冰的非线性动力学响应

本文考虑非线性、惯性和阻尼的影响, 研究了任意深度二维理想流体顶部浮冰的振动. 对相关的拟微分算子进行展开并将非线性项保留至三阶后, 完全非线性问题被简化为仅与自由面上的变量相关的三阶截断模型. 为了验证简化模型的准确性, 重点关注了自由孤立波解. 在不考虑阻尼的情况下, 采用多重尺度方法推导了三阶非线性薛定谔方程(NLS), 利用该方程预测了任意水深下原始欧拉方程中自由波包型孤立波解的存在性及三阶截断模型的准确性. 相比于Dinvay等所提出的二阶模型, 三阶截断模型的优势在于其对应的三阶NLS具有准确的非线性项系数, 能够在最小相速度附近更好地模拟冰层的动力学响应. 进一步地对自由孤立波解进行数值计算, 数值结果表明三阶截断模型在分岔曲线和孤立波波形上均与完全欧拉方程吻合良好, 准确性高于二阶截断模型. 基于三阶截断模型, 探究了匀速局域化载荷作用下的浮冰非线性动力学响应并将时间依赖解与实验测量数据进行比较, 数值计算结果与实验记录吻合良好.      

Wave formation on vertical falling liquid films

The method of integral relations is used to derive a nonlinear “two-wave” structure equation for long waves on the surface of vertical falling liquid films. This equation is valid in a wide range of Reynolds numbers and reduces to the known equations for high and low Re. Theoretical data for the fastest growing waves are compared with the experimental results on velocities, wave numbers and growth rates of the waves in the inception region. The validity of theoretical assumptions is also confirmed by the direct measurements of the instantaneous velocity profiles in a wave liquid film.     

Local film thickness and temperature distribution measurement in wavy liquid films with a laser-induced luminescence technique

Heat transfer in falling liquid film systems is enhanced by waviness. Comprehension of the underlying kinetic phenomena requires experimental data of the temperature field with high spatiotemporal resolution. Therefore a non-invasive measuring method based on luminescence indicators is developed. It is used to determine the temperature distribution and the local film thickness simultaneously. Results are presented for the temperature distribution measurement in a laminar-wavy water film with a liquid side Reynolds number of 126 flowing down a heated plane with an inclination angle of 2° at two positions in flow direction. The measured temperature distributions are used to calculate the local heat transfer coefficient for solitary waves at two positions in flow direction.     

Long waves between a subsonic gas flow-film interface flowing down an inclined plate

The present work deals with temporal stability properties of a falling liquid film down an inclined plane in the presence of a parallel subsonic gas flow. The waves are described by evolution equation previously derived as a generalization of the model for the Newtonian liquid. We confirm linear stability results of the basic flow using the Orr–Sommerfeld analysis to that obtained by long wave approximation analysis. The non-linear stability criteria of the model are discussed analytically and stability branches are obtained. Finally, the solitary wave solutions at the liquid–gas interface are discussed, using specially envelope transform and direct ansatz approach to Ginzburg–Landau equation. The influence of different parameters governing the flow on the stability behavior of the system is discussed in detail.     

Bifurcation analysis of the travelling waves on a falling power-law fluid film

Non-linear waves on the surface of a falling film of power-law fluid on a vertical plane are investigated. The waves are described by evolution equations previously derived as a generalization of the model for the Newtonian liquid. This paper presents a numerical bifurcation analysis of the steady travelling waves on a falling film described by an equation containing three parameters. It is shown that the wave regimes are sensitive to the power-law index as well as to the film parameter and the wavenumber that is typical for the Newtonian liquid.     

Numerical analysis on flow dynamics and heat transfer of falling liquid films with interfacial waves

Flow dynamics and heat transfer of falling liquid films with interfacial waves flowing on a vertical plate have been studied with originally proposed numerical simulation method. To discretize basic equations a staggered grid fixed on a physical space is employed. A small amplitude disturbance generated at inflow boundary develops to a solitary wave which consists of a large amplitude roll wave and small amplitude capillary waves. Instantaneous streamwise velocity profiles at the wave crest and trough are very different from a laminar flow. A circulation flow occurs in the roll wave and it affects temperature distributions, especially the strong effect is observed for high Prandtl number liquids. The interfacial wave enhances the heat transfer by two kinds of effects which are a film thinning effect and a convection effect. The dominating effect depends on the Prandtl number. Received on 23 December 1998     

内孤立波作用下潜深对潜体运动响应和载荷特性的影响研究

内孤立波常发生于海洋密度跃层, 因其峰高谷深、携带能量巨大, 在传播过程中会导致跃层上下的海水流动呈现剪切状态, 并引起突发性的强流. 潜体在水下悬停时极有可能会遭遇内孤立波, 由于内孤立波的流场特性, 置于跃层上下的悬浮潜体所产生运动响应和水动力载荷变化不尽相同, 甚者会出现掉深现象. 为探究潜深对波体耦合作用的影响, 基于不可压缩N-S方程和mKdV理论, 采用速度入口造波, 结合重叠网格技术和流固耦合方法, 建立了分层流中内孤立波耦合水下潜体多自由度运动的数值模型, 通过该模型分析了不同潜深下悬浮潜体的运动响应和载荷特性. 结果表明: 在内孤立波作用下, 位于密度跃层上方和跃层中的潜体顺着波的前进方向运动, 先下沉后抬升, 位于跃层下方的潜体则会逆流持续下沉; 潜体与波面的垂向距离越小, 对其纵荡、垂荡和速度的影响越显著, 而位于密度跃层中的潜体在分界面处沿着波形运动, 其运动响应和载荷变化受影响较小; 潜体在跃层上、下流体中所受水平力的方向相反, 水平力峰值小于垂向力峰值, 且位于跃层下方的潜体一直受到低头力矩, 最终导致掉深.      

Wave regimes on power-law fluid film flowing down a porous plane

Non-linear waves on the surface of a falling film of power-law fluid on a vertical porous plane are investigated. The waves are described by evolution equations generalising equations previously derived in the case of solid plane. It is shown that the slip condition on the interface between pure liquid and the porous substrate drastically changes structure of the steady waves travelling in the film.     

Numerical simulation on vapor absorption by wavy lithium bromide aqueous solution films

Numerical simulation has been made on heat and mass transfer of vapor absorption by wavy lithium bromide aqueous solution films. The velocity fields and interface positions are obtained by VOF model. Solitary waves are generated by periodically disturbed inflow boundary. Based on these, the temperature and concentration fields are obtained with a stationary interface shape. The effect of solitary waves on the heat and mass transfer across the film is investigated. It is shown that due to the mixing of circulation and stretch of large film thickness, the gradient of concentration and absorption rate decrease for solitary wave region. The region of capillary waves shows a significant amount of absorption enhancement. The percentage of absorption for the different regions is quantified.     

Selection of solitary waves in vertically falling liquid films

Two-dimensional solitary waves at the surface of a film flow down a vertical plane are considered. When the system is subjected to inlet white noise, solitary waves are formed after an inception region and interact with each other. Using open-domain simulations of reduced equation models, we investigate numerically their late time process dynamics. Close to the instability threshold, the waves synchronize themselves into bound states. For higher values of the Reynolds number, the separation distance between the waves increases and the synchronization process at work is weaker. Performing statistics, we show that the mean characteristics of the waves correspond to the minimal value of the mean film thickness along the traveling-wave branch of solutions. In this regime, synchronization occurs through the waves tails which is associated with a change of scaling of the waves features. A similar behavior is observed performing simulations in periodic domains: the selected waves maximize the mean flow rate.     

Investigations of the Marangoni effect on the regular structures in heated wavy liquid films

Formation and development of quasi-regular metastable structures within laminar-wavy falling films were studied. These structures emerge within the residual layer between large waves and could be one reason for the break up of the falling film. The temperature field of the film surface was visualised using IR-thermography. The film thickness was obtained from point measurements with the chromatic confocal imaging method and converted into a film thickness field, based on a quasi-steady assumption and IR thermography images. The thermo-capillary nature (Marangoni effect) of the regular structures was proven experimentally.