波浪环境中垂直射流紊动特性的实验研究

利用粒子图像测速技术PIV(particle imagevelocimetry)对有限水深规则波浪环境下垂直射流紊动特性进行了实验研究. 应用相位分析法从测量数据中分离出速度脉动项,用4种不同波高的波浪研究波高对射流紊动特性的影响,对紊动量的分布以及大小进行了分析. 结果表明波高对射流的紊动特性有显著影响,并且对流项对波高的变化较紊动扩散项更为敏感,紊动扩散项量值约是对流项的$1/8\sim1/3$, 在时均化的N-S方程中起的作用不可忽略.      

波浪破碎的一种混合湍流模拟模式

近岸波浪的变形与破碎,一方面影响水体与泥沙运输,另一方面对消波护岸具有指导意义.本文提出一种三维混合湍流模拟模式方法,将求解区域分为造波区、波浪传播区和消波区.造波区采用层流模式,通过基于Fluent的二次开发的UDF方法在边界进行速度造波.这种方法在给定入口速度的条件下,根据已知波高进行精准插值,从而控制水的体积分数.在波浪传播区域,采用大涡模拟进行模拟研究,在消波区,采用RANS模型并利用多孔介质消波法进行消波.模型通过VOF方法捕捉波浪破碎过程中的自由面变化.本文对波高为5.5 cm的规则波(M1)、波高为13.5cm的规则波(M3)、有效波高为7.75cm的TMA谱单向不规则波(U1)和有效波高为19cm的TMA谱多向不规则波(B5)展开了模拟研究,并与前人的相关实验结果作比较,各条件下模拟结果与实验结果吻合.模拟结果说明本文提出的模型能够准确模拟出波浪传播过程中的折射和绕射现象,并且能够捕捉的波浪破碎过程中的自由面变化,为三维波浪的传播与破碎的数值模拟提供一种模拟方法.     

波浪破碎的一种混合湍流模拟模式

近岸波浪的变形与破碎,一方面影响水体标识码运输,另一方面对消波护岸具有指导意义.本文提出一种三维混合湍流模拟模式方法,将求解区域分为造波区、波浪传播区和消波区. 造波区采用层流模式,通过基于 Fluent 的二次开发的 UDF 方法在边界进行速度造波.这种方法在给定入口速度的条件标识码据已知波高进行 精准插值,从而控制水的体积分数. 在波浪传播区域,采用大涡标识码行模拟研究,在消波区,采用 RANS 模型并利用多孔介质消波法进行消波. 模标识码 VOF 方法捕捉波浪破碎过程中的自由面变化.本文对波高为 5.5cm 的规则波 (M1)、波高为 13.5cm 的 规则波 (M3)、有效标识码 7.75cm 的 TMA 谱单向不规则波 (U1) 和有效波高为 19cm 的 TMA 谱多向不规则波 (B5) 展开标识码研究,并与前人的相关标识码果作比较,各条件下模拟结果与实验结果吻合.模拟结果说明本文提出的模型能够准确模拟出波浪传播过程中的折射和标识码象,并且能够捕捉的波浪破碎过程中的自由面变化,为三维波浪的传播与破碎的数值模拟提供一种模拟方法.      

波浪中舰船远程尾流区气泡运动研究

为准确描述波浪中舰船远程尾流区气泡运动规律,分别建立了波浪中远程尾流场气泡质心运动和径向运动模型;并结合二维深水波模型,构建了波浪中尾流气泡运动耦合模型。采用变步长龙格库塔法进行数值求解,模拟得出波浪中尾流气泡的两个运动特征规律:(1)气泡在波浪中的上浮轨迹呈螺旋状,且气泡尺度越小,受波浪牵引作用越显著;(2)波浪的波长越短,波高越高,气泡在波浪传播方向的位移越大,但气泡的存留时间受波浪的影响很小。该模型考虑了真实海况中波浪对气泡运动的影响,可作为研究舰船尾流气泡运动特性及尾流场演化规律的理论基础。     

基于P-M谱的波浪荷载数值模拟及其对跨海斜拉桥的冲击响应规律研究

针对深水桥梁在海洋中面对波浪的冲击问题,以某跨海斜拉桥为例,设定P-M谱为目标谱,依据Welch算法对海浪功率谱进行验证,基于以上理论采用fluent软件建立数值波浪水槽,根据重现期为10年、50年和100年分别计算对应的波高为6.4 m,7.2 m和9.6 m的波浪力,并将波浪力导入斜拉桥进行动力响应计算。综合分析大桥的动力响应, 结果表明,(1) 纵桥向波浪力是由横桥向波浪冲击桥墩出现绕射效应而产生,所以斜拉桥所受横桥向波浪力远大于纵桥向波浪力;(2) 随着波高的增大,斜拉桥的动力响应峰值逐渐增大,主跨跨中和塔顶动力响应峰值随波高的增长率最快,次边跨跨中动力响应峰值随波高的增长率次之,边跨跨中动力响应峰值随波高的增长率最慢;(3) 随着波浪荷载的增大,塔顶横向位移响应增幅显著,塔顶位移响应峰值最大达到0.0598 m。相关研究对大跨度深水桥梁的设计和动力分析有一定的参考意义。     

孤立波斜坡堤堤顶及后坡越浪流数值分析

孤立波往往被用来模拟海啸等浅水大波.运用基于RANS方程、VOF方法以及Goring造波方法建立了二维数值波浪水槽,进行了孤立波在简单斜坡堤上越浪过程的数值模拟.通过数值模拟,研究了不同相对堤顶宽及不同相对波高条件下对孤立波越堤流基本特征,包括堤顶前缘水层厚度及速度、堤顶末端水层厚度及速度、堤顶流厚度及速度沿程分布、后坡流厚度及速度沿程分布.依据数值实验结果,建立了堤顶流最大厚度分布及最大速度分布的经验公式.计算结果表明,堤顶流厚度沿程减小,而堤顶流速度则沿程增大;后坡流厚度先增大,随后再沿程缓慢减小;后坡流速度先减小,随后再沿程缓慢增大.     

浮泥对非线性波的运动响应简

将浮泥作为高黏度牛顿流体,与水体构成互不相混的两相流体,基于垂向二维RANS 方程(Reynolds-averaged Navier-Stokes)、SST 湍流模式(shear-stress transport),利用流体体积法(volume of fluid,VOF),结合数值造波、消波技术模拟波浪,建立了波泥耦合运动的数学模型,研究了线性波、Stokes波和椭余波作用下,浮泥界面的运动响应. 发现在各种波浪条件下,界面位移与表面位移间有相近的波形,两者间存在相位差,其值随周期和波高的变化不大;界面波高基本正比于表面波高,它随周期的增长呈单调增大的趋势.     

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型. 在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽. 分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟. 结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性. 本文工作 为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.      

风生短波谱的实验研究

根据在风水槽上波高仪测得的风生波波高时序,讨论了风生短波频谱与水波频率及摩擦风速的关系,通过实验及与他人实验结果的比较,发现风生短波频谱与摩擦风速有幂函数关系,摩擦风速的幂指数并非为常数,实验显示幂指数应是频率的函数（有近似的对数关系）．     

离岸式潜堤对海堤越浪量影响的实验研究

外海波浪作用下海堤越浪是一个复杂的水动力学过程,越浪量则是评估海堤安全的重要因素之一。离岸式潜堤作为一种典型的护岸工程结构物,已被应用于淤泥质海岸复合式海堤工程,其主要作用是促淤保滩和降低越浪量。本文针对不同的潜堤与海堤距离以及不同的潜堤高程,进行了规则波作用下越浪量系列物理模型实验,给出了各工况下的实测越浪量。实验结果表明:潜堤是否导致波浪破碎是控制海堤越浪量的主要因素之一;当潜堤使得波浪破碎时,潜堤位置对海堤越浪量有较大影响。潜堤与海堤距离接近1/2波长时越浪量偏大。     

波作用量守恒原理在波流相互作用中的应用

波和流相遇后，经过相互作用产生稳定的波流场．基于波作用量守恒原理，本文提出了解决Ｓｔｏｋｅｓ波与指数剖面流同向或反向相遇过程中波传播特性改变问题的五阶理论，给出了这种相互作用前后波参数间的关系和相应流场的变化，并与线性结果进行了比较．结果表明，当波流同向时，波变平缓；当波流反向时，波明显变陡．随着波流强度的增加，非线性的影响变得越来越明显，尤其当波同反向流相遇时，波幅显著增大，这时必须用非线性理论来考虑波流相互作用的影响     

Simulating the formation of spiral wave in the neuronal system

Some experimental results have confirmed that a spiral wave could be observed in the cortex of brain. The biological Hodgkin–Huxley neurons are used to construct a regular network with nearest-neighbor connection, artificial line defects are generated to block the traveling wave in the network, and the potential mechanism for formation of spiral wave is investigated. A target wave is generated in a local area by imposing two external forcing currents with diversity (I ₀?I ₁) in different areas of the network. It is confirmed that spiral wave could be induced by the defects even if no specific initial values are used. A single perfect spiral wave can occupy the network when the coupling intensity exceeds certain threshold; otherwise, a group of spiral waves emerges in the network. Certain channel noise can enhance the diversity (I ₀?I ₁) for generating target wave, and then spiral waves are induced by blocking the target wave with defects under no-flux and/or periodic boundary conditions in the network.     

An experimental investigation into non-linear wave loading on horizontal axis tidal turbines

Tidal turbines are subject to large hydrodynamic loads from combinations of currents and waves, which contribute significantly to fatigue, extreme loading and power flow requirements. Physical model testing enables these loads and power fluctuations to be assessed and understood in a controlled and repeatable environment. In this work, a 1:15 scale tidal turbine model is utilised to further the fundamental understanding of the influence of waves on tidal turbines. A wide range of regular waves are generated in both following-current and opposing-current conditions. Wave frequencies range from 0.31 Hz to 0.55 Hz & wave heights from 0.025 m to 0.37 m in a fixed 0.81 m/s current velocity. Waves are selected and programmed specifically to facilitate frequency domain analysis, and techniques are employed to isolate the effect of non-linear waves on turbine power and thrust.Results demonstrate that wave action induces large variations in turbine power and thrust compared to current only conditions. For the range of conditions tested, peak values of thrust and power exceed current-only values by between 7%–65% and 13%–160% respectively. These wave-induced fluctuations are shown to increase with wave amplitude and decrease with wave frequency. Following wave conditions exhibit greater variations than opposing for waves with the same wave height and frequency due to the lower associated wavenumbers.A model is developed and presented to aid the understanding of the high-order harmonic response of the turbine to waves, which is further demonstrated using steady state coefficients under assumptions of pseudo-stationarity. This approach is proven to be effective at estimating wave-induced power and thrust fluctuations for the combinations of waves, currents and turbine state tested. The outcome of which shows promise as a rapid design tool that can evaluate the effect of site-specific wave–current conditions on turbine performance.     

Formation of solitary waves on gas-sheared liquid layers

The origin of solitary waves on gas-liquid sheared layers is studied by comparing the behavior of the wave field at sufficiently low liquid Reynolds number, R_L, where solitary waves are observed to form, to measurements at higher R_L where solitary waves do not occur. Observations of the wave field with high-speed video imaging suggest that solitary waves, which appear as a secondary transition of the stratified gas-liquid interface, emanate from existing dominant waves, but that not all dominant waves are transformed. From measurements of interface tracings it is found that for low R_L, waves which have amplitude/substrate depth (a/h) ratios of 0.5–1 occur while for higher R_L, no such waves are observed. A comparison of amplitude/wavelength ratios shows no distinction for different R_L. Consequently, it is conjectured that solitary waves originate from waves with sufficiently large a/h ratios; this change of form being similar to wave breaking. The dimensionless wavenumber is found to be smaller at low R_L, where solitary waves are observed. This suggests that perhaps, larger precursor (to solitary wave) waves are possible because the degree of dispersion, which acts to break waves into separate modes, is lower.     

Extended Mild-Slope Equation

ＩｎｔｒｏｄｕｃｔｉｏｎＡｃｃｕｒａｔｅｍｏｄｅｌｌｉｎｇｏｆｓｕｒｆａｃｅｗａｖｅｄｙｎａｍｉｃｓｉｎｃｏａｓｔａｌｒｅｇｉｏｎｓｈａｓｂｅｅｎｔｈｅｇｏａｌｏｆｍｕｃｈｒｅｃｅｎｔｒｅｓｅａｒｃｈ ,ｗｈｉｃｈｈａｓｂｅｅｎｓｕｍｍａｒｉｚｅｄｉｎｓｕｒｖｅｙｓｂｙＤｉｎｇｅｍａｎｓ( 1 997) [1]ａｎｄＫｉｒｂｙ( 1 997) [2 ].Ｔｈｅｒｉｃｈｎｅｓｓｏｆｃｏａｓｔａｌｗａｖｅｄｙｎａｍｉｃｓａｒｉｓｅｓｆｒｏｍｔｈｅｓｔｒｏｎｇａｍｂｉｅｎｔｃｕｒｒｅｎｔｓａｎｄｔｈｅｗｉｄｅｖａｒｉａｔｉｏｎｓ…     

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     

The Patterns of Surface Capillary-gravity Short-crested Waves with Uniform Current Fields in Coastal Waters

A fully three-dimensional surface gravity-capillary short-crested wave system is studied as two progressive wave-trains of equal amplitude and frequency, which are collinear with uniform currents and doubly-periodic in the horizontal plane, are propagating at an angle to each other. The first- and second-order asymptotic analytical solutions of the short-crested wave system are obtained via a perturbation expansion in a small parameter associated with the wave steepness, therefore depicting a series of typical three-dimensional wave patterns involving currents, shallow and deep water, and surface capillary waves, and comparing them with each other.The project supported by the Foundation for the Author of National Excellent Doctoral Dissertation of China (200428), the National Natural Science Foundation of China (10272072 and 50424913), the Shanghai Natural Science Foundation (05ZR14048), and the Shanghai Leading Academic Discipline Project (Y0103). The English text was polished by Yunming Chen.     

Solitary waves in a peridynamic elastic solid

The propagation of large amplitude nonlinear waves in a peridynamic solid is analyzed. With an elastic material model that hardens in compression, sufficiently large wave pulses propagate as solitary waves whose velocity can far exceed the linear wave speed. In spite of their large velocity and amplitude, these waves leave the material they pass through with no net change in velocity and stress. They are nondissipative and nondispersive, and they travel unchanged over large distances. An approximate solution for solitary waves is derived that reproduces the main features of these waves observed in computational simulations. It is demonstrated by numerical studies that the waves interact only weakly with each other when they collide. Wavetrains composed of many non-interacting solitary waves are found to form and propagate under certain boundary and initial conditions.     

A special solution of wave dissipation by finite porous plates

The reflection and transmission of water waves caused by a small amplitude incident wave through finite fine porous plates with equal spacing and permeability in an infinitely long open channel of constant water depth and zero slope are studied. A special solution is obtained when the distance between the two neighbouring plates is an integral multiple of the half-wavelength of the incident wave. It is found, that when the dimensionless porous-effect parameter G_0 is equal to half the total plate number, the wave dissipation reaches a maximum, and only 50% of the incident wave energy remains in the reflected and transmitted waves. Meanwhile, the reflected and transmitted waves have the same amplitude.     

非均匀水流中非线性波传播的数值模拟

以一种考虑波流相互作用的新型{Boussinesq}型方程为控制方程组， 采用五阶{Runge}-{Kutta}-{England}格式离散时间积分，采用七点 差分格式离散空间导数，并通过采用恰当的出流边界条件，从而建立了非均匀水流中非线性 波传播的数值模拟模型. 通过对均匀水流与水深水域内和潜堤地形上存在弱流或强流时波浪 传播的数值模拟，说明模型能有效地反映水流对波浪传播的影响.