Structurization and restructurization of a homogeneous suspension in a standing wave field

The process of formation of a spatially ordered fine structure in an originally uniform suspension of small-size particles in the field of two-dimensional standing surface gravity waves in a rectangular reservoir oscillating in the vertical direction is for the first time recorded in experiments. The parameters of structured regions are determined in reservoirs with both separate vortical ripples and their groups and with rough sandy and smooth regions of the bottom for the first and second wave modes. It is established that isolated structured regions are formed near the bottom or the free surface and gradually fill the fluid throughout its entire thickness. Minimum scales of the structures are on the Stokesian range, while maximum scales are determined by the geometry and the energy of the main flow. In all experiments performed the structurization was conserved during the entire interval of wave motion of the suspension and vanished during a few minutes after the parametric excitation of the waves has been stopped.     

Migrating sea ripples

Ripple formation under sea waves is investigated by means of a linear stability analysis of a flat sandy bottom subject to the viscous flow which is present in the boundary layer at the bottom of propagating sea waves. Nonlinear terms in the momentum equation are retained to account for the presence of a steady drift. Hence the work by Blondeaux is extended by considering steeper waves and/or less deep waters. Second order effects in the sea wave steepness are found to cause neither destabilizing nor stabilizing effects on the process of ripple formation. However, because of the presence of a steady velocity component in the direction of wave propagation, ripples are found to migrate at a constant rate which is predicted as function of sediment and wave characteristics. The analysis assumes the flow regime in the bottom boundary layer to be laminar and the results are significant for ripples at the initial stage of their formation or for mature ripples of small amplitude (rolling-grain ripples). A comparison of the theoretical findings with laboratory experiments supports the reliability of the approach and of the theoretical results.     

黄茅海波-流共同作用下的三维悬沙数值模拟

针对黄茅海河口区波流相互作用显著、三维空间结构明显及泥沙运动复杂等特点,联合第3代海浪模式SWAN,建立了波流共同作用下的三维悬沙数学模型.在水流计算中,考虑了波浪产生的辐射应力影响,波浪依赖的表面风应力影响及波浪影响下的底部剪切应力,并采用k-kl紊流闭合模型提供垂向涡粘系数和垂向扩散系数.悬沙扩散方程中的源函数利用切应力法确定,泥沙沉降速度考虑絮凝作用,从而提高黄茅海悬沙场数值模拟精度.通过含沙量实测资料验证,模拟值与实测值符合较好.在模型验证的基础上,讨论了不考虑波浪与考虑波浪两种情况下黄茅海的悬沙分布特征,为工程实际研究提供了依据.     

Stability and capillary dynamics of circular vortex sheets

We investigate the motion of circular vortex sheets with surface tension. A linear stability analysis shows that high modes of the circular vortex sheet are stabilized by surface tension, and the sheet is stable if surface tension is larger than a critical value. The modes of perturbations, n = 1 and 2, are always stable, regardless of surface tension, and the mode n = 3 is also stable for large surface tension. The numerical results show that a stable vortex sheet rotates and oscillates weakly. The oscillations of each mode of the interface mainly consist of two travelling waves of different frequencies in time. The amplitude and the period of the oscillation depend on the mode of the perturbation and surface tension. We also perform long-time computations for the unstable evolution of circular sheets. For a high Weber number, ripples are produced on the sheets, as well as pinching and self-intersection. It is found that the appearance of ripples is associated with the growth of noise. For an intermediate Weber number, the sheet evolves to an exotic structure with small spikes on the fingers, while for a low Weber number, it is nonlinearly stable.     

Waves on vortex cores

Here we discuss the types of waves which can be supported on compact regions of vorticity. This is a subject first studied by Lord Kelvin for waves propagating along the vortex lines of the Rankine and hollow-core vorticity distributions. Kelvin's major interest was in the stability of vortex rings and the numerous resurgencies in interest have usually been driven by practically important phenomena, e.g., the observations of vortex breakdown and waves on tornado and aircraft vortices in the early 1960's and more recently in technologically and geophysically significant flows and on the quantised vortices of super fluid HeII.

The major wave-types of interest are of varicose, helicoidal and fluted form and represent a periodic swelling and contraction, a bending and a “krinkling” of the core, respectively. The first two propagate along and the third around the vortex lines. All have been studied theoretically, experimentally and numerically in the limit of small wave amplitude and their major characteristics are now clear. Of particular interest is the extension of these results to the non-linear regime in which case the two first types are known to exhibit solitary wave or soliton characteristics in certain parameter ranges. It is these non-linear waves which often dominate observations of vortex flows both in nature and in technological applications and which have caused much controversy in the interpretation of results found under complex circumstances of flow and apparatus geometry.     

Non‐hydrostatic 3D free surface layer‐structured finite volume model for short wave propagation

In this paper a layer‐structured finite volume model for non‐hydrostatic 3D environmental free surface flow is presented and applied to several test cases, which involve the computation of gravity waves. The 3D unsteady momentum and mass conservation equations are solved in a collocated grid made of polyhedrons, which are built from a 2D horizontal unstructured mesh, by just adding several horizontal layers. The mesh built in such a way is unstructured in the horizontal plane, but structured in the vertical direction. This procedure simplifies the mesh generation and at the same time it produces a well‐oriented mesh for stratified flows, which are common in environmental problems. The model reduces to a 2D depth‐averaged shallow water model when one single layer is defined in the mesh. Pressure–velocity coupling is achieved by the Semi‐Implicit Method for Pressure‐Linked Equations algorithm, using Rhie–Chow interpolation to stabilize the pressure field. An attractive property of the model proposed is the ability to compute the propagation of short waves with a rather coarse vertical discretization. Several test cases are solved in order to show the capabilities and numerical stability of the model, including a rectangular free oscillating basin, a radially symmetric wave, short wave propagation over a 1D bar, solitary wave runup on a vertical wall, and short wave refraction over a 2D shoal. In all the cases the numerical results are compared either with analytical or with experimental data. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

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

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

Circulation flow patterns in the Faraday waves in the presence of damper plates

The effect of transverse damper plates on the restructurization of a suspension with originally homogeneous concentration in the field of standing surface waves in a vertically oscillating rectangular reservoir is for the first time experimentally investigated. The flow pattern evolution and the formation of the fine flow structure are traced in reservoirs with smooth walls, with one vertical deflecting plate, and with two plates mounted in parallel or antiparallel at the center of the larger side. The symmetry of the surface wave flow is determined by the plate configuration.     

The loss of stability of laminar flow in open channel and the mechanism of sand ripple formation

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｒｉｎｇｔｈｅｔｉｍｅｏｆｂｅｄｓｕｒｆａｃｅｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔａｔｉｏｎ,ｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｆｂｅｄｌｏａｄｍａｋｅｔｈｅｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｖａｒｉｏｕｓｄｉｆｆｅｒｅｎｔｆｏｒｍｓｏｎｔｈｅｒｉｖｅｒｂｅｄｓｕｒｆａｃｅ,ｔｈｕｓｔｈｅｆｏｒｍｏｆｂｅｄｓｕｒｆａｃｅｃｈａｎｇｅｓｃｏｎｓｔａｎｔｌｙ .Ｔｈｉｓｋｉｎｄｏｆｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｎｔｈｅｂｅｄｓｕｒｆ…     

Wave attenuation in a vessel equipped with damper plates

The effect of vertical damper plates mounted at the center of a rectangular vessel normal or at an angle to a wave flow is experimentally investigated for the first mode of the standing surface waves excited at the parametric resonance. The variation of the resonance curves and the wave attenuation degree are discussed. The fluid depth effect on the wave motion damping is evaluated.     

Reflection of a spherical blast wave from a planar surface

Numerous authors have carried out rather extensive studies in the last twenty to thirty years of the problem of the interaction of shock and blast waves with obstacles in their paths. Owing to the complexity of the problem, they assumed certain limiting cases for the shock wave interactions in which the parameters behind the shock wave were usually taken to be constants. The first wave diffraction studies involving variable parameters behind the front were presented in [1, 2], wherein a development of the theory of “short waves” (blast waves at a substantial distance from the center of an explosion) and their reflection from a planar surface was given. The theory of short waves assumes that the jump in pressure at the wave front and the region over which the parameters vary are small. The problem concerning reflection of a blast wave from a surface was also considered in [3, 4], wherein a solution in the region behind the reflected wave was obtained at initial times. The initial stage in the reflection of a blast wave from a planar, cylindrical, or spherical surface (the one-dimensional case) was studied in [5]. In this paper we investigate the interaction of a spherical blast wave, resulting from a point explosion, with a planar surface; we consider both regular and non-regular reflection stages. In solving this problem we use S. L. Godunov's finite-difference method. We obtain numerical solutions for various values of the shock strength at the instant of its encounter with the surface. We present the pressure fields in the flow regions, the pressure distribution over the surface at various instants of time, and the trajectories of the triple point. The parameter values at the front of the reflected wave are compared with results obtained from the theory of regular reflection of shock waves.     

细颗粒泥沙净冲刷和输移的大涡模拟研究

在传统水沙输移数值模拟研究中一般采用雷诺时均模拟技术(Reynolds-averaged simulation,RANS).与RANS相比,大涡模拟技术(large eddy simulation,LES)能够更加精确反映细部流动结构,计算机的发展使得采用LES探讨水流和泥沙运动规律成为可能.本文尝试给出净冲刷条件下悬沙计算的边界条件,采用动态亚格子模式对循环槽道和长槽道中的水流运动和泥沙输移进行了三维大涡模拟研究.利用直接数值模拟(direct numerical simulation,DNS)结果对LES模型进行了率定,计算结果符合良好,在此基础上初步探讨了泥沙浓度、湍动强度和湍动通量等的分布特征.结果表明,净冲刷条件下输沙平衡时泥沙浓度符合Rouse公式分布,单向流动中泥沙浓度沿着流向逐渐增大.泥沙浓度湍动强度和湍动通量都在近底部达到最大值,沿着垂向迅速减小.湍动黏性系数和扩散系数基本上在水深中间处达到最大.湍动Schmidt数沿着水深方向不是常数,在近底部和自由水面附近较大,在水深中间处较小.     

Wavelength selection of ripples in a vertically vibrating dynamically thick granular layer due to density-wave refraction

A numerical study was made on the wavelength selection mechanism of the ripples observed on the surface of a granular layer that is oscillated vertically. Multiple collisions of the one-dimensional array of beads show a time-dependent particle distribution, which induces a density wave. The magnitude of the wave velocity is estimated by the theory of elasticity, which reveals the refraction of the density wave in a quasi-two-dimensional granular layer. Our theory explains how the vertical excitation of particles determines the horizontal characteristic scale on the surface, or the wavelength of the ripples, which is applicable even to the saturation regime of a dynamically thick granular layer, e.g., the one that is thick enough to allow an immobilized region in the lower part.     

Modelling solitons under the hydrostatic and Boussinesq approximations

An examination of solitary waves in 3D, time‐dependant hydrostatic and Boussinesq numerical models is presented. It is shown that waves in these models will deform and that only the acceleration term in the vertical momentum equation need be included to correct the wave propagation. Modelling of solitary waves propagating near the surface of a small to medium body of water, such as a lake, are used to illustrate the results. The results are also compared with experiments performed by other authors. Then as an improvement, an alternative numerical scheme is used which includes only the vertical acceleration term. Effects of horizontal and vertical diffusion on soliton wave structure is also discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Measuring the two-dimensional structure of a wavy water surface optically: A surface gradient detector

A new method of measuring the slopes of a water surface covered with short waves is developed. A camera is placed far above the water surface looking downward so that it receives only approximately vertical rays of light emerging from the water surface from a source below. A large lens is positioned horizontally underwater. A plane light source in the form of a translucent colored screen is placed horizontally in the focal plane below this lens. Corresponding to each value of water surface slope, regardless of observer position, there is one and only one point of origin on the color screen from which light rays can enter the camera. When the color screen has a suitable two-dimensional color pattern, we are able to detect the gradient of the surface elevation throughout the field of view of the camera. This refraction slope detector has been used to find statistical properties of short wind waves in a wind-wave channel where a broad angular beam width of capillary ripples and short gravity waves contribute to the surface slopes. In these experiments waves were generated by winds ranging from 5 m/s to 10 m/s at a fetch of 24 m. The wavenumber spectra of short wave slopes have two distinguishing features: a dip at the capillary-gravity transition and steep slopes in the capillary range. Surface shapes resembling the shape of solitary capillary-gravity waves have been found from profiles of wave elevation deduced by integration of the elevation gradient.We are especially grateful for the advice of Dr. M. Gharib on the use of the HSI color system. John Lyons provided expert help in the laboratory and materials for and advice on photography. We thank the staff of the SIO Hydraulics Laboratory for making the wind-wave channel available for our use, and the staff of UCSD library for enabling us to use the Barneyscanner photometer-digitizer. We thank an anonymous reviewer who pointed out a numerical error and improved the clarity of the text.     

Experimental investigation of secondary steady flows in Faraday surface waves

It is shown that in the two-dimensional Faraday surface waves excited in a vertically oscillating rectangular water-filled vessel there is a system of secondary circulatory flows that occupies the entire fluid volume between the vessel bottom and the free surface. In parallel with the oscillations at the wave frequency, the fluid particles are slowly displaced in accordance with these circulatory flows. The secondary flow velocity field is measured and the trajectories of individual fluid particles in the standing wave are determined. The experimental data are compared with the Longuet-Higgins model. It is shown that the initial stage of formation of regular structures on the surface of a sediment layer on the vessel bottom may be related with the presence of secondary circulatory flows.