Internal wave radiation by a turbulent fountain in a stratified fluid

Large eddy simulation is applied to model a fountain in a density-stratified fluid. The fountain is formed, as a vertical turbulent jet penetrates through a pycnocline. The jet flow is initiated by the formulation of a boundary condition in the form of an upward neutral-buoyancy fluid flow with the Gaussian axisymmetric mean-velocity profile and a given fluctuation level. It is shown that at a Froude number Fr higher than a certain critical value the fountain executes self-oscillations accompanied by internal wave generation within the pycnocline. The predominant self-oscillation mode is axisymmetric, when the fountain top periodically breaks down generating internal wave packets traveling toward the periphery of the computation domain. The characteristic frequency of the internal waves coincides with that of the fountain top oscillations and monotonically decreases with increase in Fr. The Fr-dependence of the fountain top oscillation amplitude obtained in the numerical solution is in good agreement with the predictions of the theoretical Landau model for the instability mode in the soft self-excitation regime.     

Surface characters of internal waves generated by Rankine ovoid

A linear theory on the internal waves generated in the stratified fluid with a pycnocline is presented in this paper. The internal wave fields such as the velocity fields in the stratified fluid and velocity gradient fields at the free surface are also investigated by means of the theoretical and numerical method. From the numerical results, it is shown that the internal wave generated by horizontally moving Rankine ovoid is a sort of trapped wave which propagates in a wave guide, and its waveform is a kind of Mach front-type internal wave in the pycnocline. Influence of the internal wave on the flow fields at the free surface is represented by the velocity gradient fields resulted from the internal waves generated by motion of the Rankine ovoid. At the same time, it is also shown that under the hypothesis of inviscid fluid, the synchronism between the surface velocity gradient fields at the free surface and the internal wave fields in the fluid is retained. This theory opens a possibility to study further the modulated spectrum of the Bragg waves at the free surface.The project supported by the National Natural Science Foundation of China (40576010). The English text was polished by Keren Wang.     

Laboratory experiments on internal wave interactions with a pycnocline

Laboratory experiments have been performed to investigate the interaction of internal waves with a pycnocline. An oscillating cylinder generated internal wave beams, which were observed using the synthetic schlieren technique. Internal waves incident on the pycnocline layer excited higher-frequency modes. In the absence of shear, a discrete spectrum of harmonic modes was generated due to nonlinear effects. These harmonic modes might play a role in the formation of internal solitary waves which have been observed in ocean pycnoclines. With shear, a continuous spectrum of excited modes was found.     

Experimental Investigation of the Generation of Internal Waves by a Vertical Cylinder in a Near-Surface Pycnocline

A bench study of the amplitudes, mode composition, and phase structure of the internal waves generated by a vertical cylinder in the presence of a near-surface pycnocline has been performed; the pycnocline took the form of a stratified fluid layer located between two quasi-homogeneous layers of thicknesses h ₁ and h ₂=2h ₁. In the experiments, the cylinder traveled at velocities critical with respect to internal wave generation. Different cases of model submergence relative to the pycnocline are considered. The dependence of the mode structure and the amplitude-phase characteristics of the forced internal waves on the body velocity and its relative submergence is analyzed. The parameters of both steady and unsteady wave systems are studied.The data obtained make it possible to predict the forced wave parameters and the critical body velocities for given model dimensions and pycnocline parameters.     

Diffraction of internal waves by a circular cylinder near the pycnocline

Propagation of internal waves over a circular cylinder under the conditions of a continuous stratification characterized by the presence of a high-gradient density layer (the pycnocline) of finite thickness is studied. The dependences of the coefficent of wave propagation on the wavelength of the first-mode incident wave for various thicknesses of the pycnocline are obtained. In the diffraction of internal waves, substantial nonlinear effects are shown to occur, which result in the appearance of waves of double oscillation frequency compared to the frequency of the incident waves. The generation coefficient for these waves is found. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 79–85, March–April, 1999.     

Self-oscillatory spouting regimes of plane vertical submerged heavy-fluid jets in setups with near-bottom discharge

New results of the series of experimental and numerical investigations of the self-oscillatory regimes of plane vertical jet spouting from beneath the free surface of a heavy incompressible fluid in reservoirs of limited dimensions are presented. The experiments were performed on a rectangular-in-plan setup with near-bottom regime of fluid discharge. For several widths of the near-bottom orifices in the end walls of the setup the dependence of the self-oscillation period on the jet flow rate is studied. A considerable difference of these dependences from those earlier obtained in the case of the spouting with fluid discharge over a weir is found to exist. On certain ranges of control parameters it is established that the fountain self-oscillation periods are similar in value to those of natural oscillations of standing waves in the setup. A fairly narrow jet velocity range is revealed on which a hysteresis effect, that consists in the difference between the flow-rate-dependences of the self-oscillation period with gradual increase or decrease in the flow rate, is observable. The results of numerical calculations carried out using the STAR-CD software package are in good agreement with the experimental data.     

Solitary waves and conjugate flows in a three-layer fluid

Interfacial symmetric solitary waves propagating horizontally in a three-layer fluid with constant density of each layer are investigated. A fully nonlinear numerical scheme based on integral equations is presented. The method allows for steep and overhanging waves. Equations for three-layer conjugate flows and integral properties like mass, momentum and kinetic energy are derived in parallel. In three-layer fluids the wave amplitude becomes larger than in corresponding two-layer fluids where the thickness of a pycnocline is neglected, while the opposite is true for the propagation velocity. Waves of limiting form are particularly investigated. Extreme overhanging solitary waves of elevation are found in three-layer fluids with large density differences and a thick upper layer. Surprisingly we find that the limiting waves of depression are always broad and flat, satisfying the conjugate flow equations. Mode-two waves, obtained with a periodic version of the numerical method, are accompanied by a train of small mode-one waves. Large amplitude mode-two waves, obtained with the full method, are close to one of the conjugate flow solutions.     

Internal gravity waves in critical generation modes and in the vicinity of trajectories of motion of perturbation sources

The field of internal gravity waves in a layer of an arbitrary stratified fluid is studied for critical generation modes and in the vicinity of trajectories of motion of the perturbation sources. The exact solutions describing the structure of a separate mode of the wave field in the vicinity of the perturbation source in the critical generation modes are investigated, and expressions for the total field representing the sum of all wave modes are obtained. In the vicinity of the trajectories of the perturbation sources, asymptotic representations of the eigenfunctions and eigenvalues of the basic vertical spectral problem of internal waves are constructed in the approximation of large wave numbers and asymptotic expressions for a separate mode of the wave field are obtained that describe the spatial structure and features of the fields of internal gravity waves. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 70–79, September–October, 2008.     

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

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

Generation of internal waves by a turbulent jet in a stratified fluid

The process of generation of internal waves by an initially cylindrical, turbulent jet with a Gaussian profile of the average horizontal velocity component in a fluid with stable linear density stratification is investigated by direct numerical simulation. It is shown that on time intervals Nt < 30, where N is the buoyancy frequency, the vertical velocity pulsations collapse, which is accompanied by the generation of internal waves whose spatial period is close to the wavelength of the spiral mode of jet instability in a homogeneous fluid. The wave dynamics and kinematics can be satisfactorily described by the linear theory for a pulsed source and their parameters are in good agreement with the parameters of the “coherent” internal waves generated by a stratified wake in a laboratory experiment. At large times the wave generation ceases and the variations of the fluid density are localized in the neighborhood of the centers of large-scale vortices formed in the horizontal plane in the neighborhood of the jet.     

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

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

侧加热腔体内重力波演化过程的数值模拟

利用二维数值模拟的方法研究了侧加热腔体内的自然对流．基于数值模拟结果,描述了水平热入侵流(intrusion)的整个演化过程,并对该过程的物理机制进行了讨论．结果表明：当热入侵流抵达腔体冷壁后,由于冷壁无法卷入所有的热入侵流,热入侵流在冷上角堆积并产生一个反向流动,在冷壁边界层附近形成一个顺时针涡,该涡在浮力效应驱动下可返回热壁,并在腔体的冷热壁之间形成了腔体尺度的流体振荡,即内重力波．     

Generation of mode-2 internal waves in a two-dimensional stratification by a mode-1 internal wave

The generation of mode-2 nonlinear internal waves (IWs) by the evolution of a mode-1 IW in a two-dimensional stratification is investigated. A generation model accounting for intermodal interaction is derived based on a multi-modal approach in a weakly nonlinear and non-hydrostatic configuration. The generation model is numerically solved to simulate the evolution of mode-1 and mode-2 IWs in an inhomogeneous pycnocline. The numerical experiments confirm that mode-2 IWs are generated due to linear and nonlinear intermodal interaction. The mode-2 IW continues growing and gradually separates with the mode-1 IW during the generation process. The numerical results suggest that the pycnocline strength or thickness prominently affects the generation of mode-2 IWs, followed by pycnocline depth. A weakening or thinning pycnocline favors the generation of mode-2 IWs by evidently enhancing linear and nonlinear intermodal interaction, whereas a shoaling pycnocline favors a rapid growth rate mainly by enhancing linear intermodal interaction. The wave amplitude of an initial mode-1 IW strongly affects the generation of mode-2 IWs and increasing it can noticeably enlarge mode-2 IWs.     

The propagation of a passive admixture from a local instantaneous source in a turbulent mixing zone

The propagation of a passive admixture from a local instantaneous source in a turbulent mixing zone in a stable stratified fluid (two-dimensional problem) is simulated numerically. The location of the source does not coincide with the center of the turbulent zone. The calculation results indicate that the average admixture concentration distribution depends significantly on the initial data. The location of the maximum concentration in homogeneous and linearly stratified fluid is fairly slowly displaced toward the center of the zone. Calculations in a pycnocline show that situations are possible when the propagation of a passive admixture is largely determined by a convective flow generated by the turbulent mixing zone.     

Single cavitation bubble generation and observation of the bubble collapse flow induced by a pressure wave

This study utilizes a U-shape platform device to generate a single cavitation bubble for a detailed analysis of the flow field characteristics and the cause of the counter jet during the process of bubble collapse caused by sending a pressure wave. A high speed camera is used to record the flow field of the bubble collapse at different distances from a solid boundary. It is found that a Kelvin–Helmholtz vortex is formed when a liquid jet penetrates the bubble surface after the bubble is compressed and deformed. If the bubble center to the solid boundary is within one to three times the bubble’s radius, a stagnation ring will form on the boundary when impinged by the liquid jet. The fluid inside the stagnation ring will be squeezed toward the center of the ring to form a counter jet after the bubble collapses. At the critical position, where the bubble center from the solid boundary is about three times the bubble’s radius, the bubble collapse flow will vary. Depending on the strengths of the pressure waves applied, the collapse can produce a Kelvin–Helmholtz vortex, the Richtmyer–Meshkov instability, or the generation of a counter jet flow. If the bubble surface is in contact with the solid boundary, the liquid jet can only move inside-out without producing the stagnation ring and the counter jet; thus, the bubble collapses along the radial direction. The complex phenomenon of cavitation bubble collapse flows is clearly manifested in this study.     

逆向喷流流场模态分析及减阻特性研究

逆向喷流减阻的基本原理是利用逆向高速喷流与飞行器绕流的相互作用，使飞行器周围的流场结构发生变化，致使飞行器的气动特性发生改变，从而改善飞行器的气动性能。利用数值模拟方法对轴对称球头、截锥的逆向喷流流场开展了研究，考虑了高温非平衡化学反应对流场的影响。模拟了球头和截锥在不同总压比时流场不同的模态：长穿透流模态（LPM）和短穿透流模态（SPM），得到了不同模态下钝体表面压力、气动力系数和不同模态之间转换的瞬态效应．简单分析了喷流在减阻方面的应用，给出了几个喷口参数与减阻效率之间的关系，提出了喷流减阻工程应用时应考虑的主要因素。     

Effect of fluid boundaries on flow around a low-aspect ratio wing

Large eddy simulation is applied to investigate flow around a finite-aspect-ratio wing in motion in a stratified fluid in different positions. It is shown that the wing position relative to a pycnocline has a considerable effect on the hydrodynamic wake development and the internal wave structure.     

Viscous fluid flow induced by rotational-oscillatory motion of a porous sphere

Viscous fluid flow induced by rotational-oscillatorymotion of a porous sphere submerged in the fluid is determined. The Darcy formula for the viscous medium drag is supplementedwith a term that allows for the medium motion. The medium motion is also included in the boundary conditions. Exact analytical solutions are obtained for the time-dependent Brinkman equation in the region inside the sphere and for the Navier–Stokes equations outside the body. The existence of internal transverse waves in the fluid is shown; in these waves the velocity is perpendicular to the wave propagation direction. The waves are standing inside the sphere and traveling outside of it. The particular cases of low and high oscillation frequencies are considered.     

Flow pattern behind a submerged body near a pycnocline

A flow past an ellipsoid immersed in a flow of a viscous stratified fluid is studied using a Large Eddy Simulation (LES) method for different locations of the body relative to a density discontinuity. It is shown that, even when the internal Froude number is large, for small angles of attack of the body the stratification affects its drag force. When the body is located above the pycnocline, the presence of the discontinuity results in the ascending of a vortex filament. A spectral analysis showed that the internal waves in the body wake have a multimode structure.     

Unsteady phenomena of an oscillating turbulent jet flow inside a cavity: Effect of aspect ratio

Self-sustained oscillatory phenomena in confined flow may occur when a turbulent plane jet is discharging into a rectangular cavity. An experimental set-up was developed and the flow analysis has been made using mainly hot-wire measurements, which were complemented by visualisation data. Previous studies confirmed that periodic oscillations may occur, depending on the location of the jet exit nozzle inside the cavity, and also the distance between the side-walls. The present study deals with the symmetrical interaction between a turbulent plane jet and a rectangular cavity and the influence of the geometrical characteristics of the cavity on the oscillatory motion. The size and aspect ratio of the cavity were varied together with the jet width compared to that of the cavity. The study is carried out both numerically and experimentally. The numerical method solves the unsteady Reynolds averaged Navier–Stokes equations (URANS) together with the continuity equation for an incompressible fluid. The closure of the flow equations system is achieved using a two-scale energy-flux model at high Reynolds number in the core flow coupled with a wall function treatment in the vicinity of the wall boundaries. The fundamental frequency of the oscillatory flow was found to be practically independent of the cavity length. Moreover, the oscillations are attenuated as the cavity width increases, until they disappear for a critical value of the cavity width. Contour maps of the instantaneous flow field are drawn to show the flow pattern evolution at the main phases of oscillation. They are given for several aspect ratios of the cavity, keeping constant values for the cavity width and the jet thickness. The proposed approach may help to investigate further the oscillation mechanisms and the entrainment process occurring in pressure driven jet–cavity interactions.