Dynamic parameters in models of atmospheric vortex structures

Vortex simulation and the computation of fields of dynamic parameters of vortex structures (velocity, rotor velocity, and helicity) are carried out with the use of exact hydrodynamic equations in a cylindrical coordinate system. Components of centripetal and Coriolis accelerations are taken into account in the initial equations. Internal and external solutions are defined. Internal solutions ignore the disturbances of the pressure field, but they are considered in external solutions. The simulation is carried out so that the effect of accounting for spatial coordinates on the structure of the above fields is pronounced. It is shown that the initial kinetic energy of rotating motion transforms into the kinetic energy of radial and vertical velocity components in models with centripetal acceleration. In models with Coriolis acceleration, the Rossby effect is clearly pronounced. The method of an “inverse problem” is used for finding external solutions, i.e., reconstruction of the pressure field at specified velocity components. Computations have shown that tangential components mainly contribute to the velocity and helicity vortex moduli at the initial stage.     

Scour and Burial Mechanics of Objects in the Nearshore

A process-based, numerical, hydrodynamic vortex lattice mine scour/burial model (VORTEX) is presented that simulates scour and burial of objects of arbitrary shape resting on a granular bed in the nearshore. There are two domains in the model formulation: a far-field where burial and exposure occur due to changes in the elevation of the seabed and a near-field involving scour and transport of sediment by the vortices shed from the object. The far-field burial mechanisms are based on changes in the equilibrium bottom profiles in response to seasonal changes in wave climate and accretion/erosion waves spawned by fluxes of sediment into the littoral cell. The near-field domain consists of one grid cell extracted from the far-field that is subdivided into a rectangular lattice of panels having sufficient resolution to define the shape of the object. The vortex field induced by the object is constructed from an assemblage of horseshoe vortices excited by local pressure gradients and shear over the lattice panels. The horseshoe vortices of each lattice panel release a pair of vortex filaments into the neighboring flow. The induced velocity of these trailing vortex filaments causes scour of the neighboring seabed and induces hydrodynamic forces on the object. Scour around the object and its subsequent movement into the scour depression contribute to burial, while far-field changes in local sand level may increase burial depth or expose the object. Scour and burial predictions of mines and mine-like objects were tested in field experiments conducted in the nearshore waters off the Pacific coast of California at Scripps Pier, the Gulf Coast of Florida at Indian Rocks, and off the Atlantic coast of Massachusetts at Martha's Vineyard. Model predictions of mine scour and burial are in reasonable agreement with field measurements and underwater photographs.     

基于客观方法的海洋三维涡旋核心线提取

涡旋核心线是海洋中尺度涡旋结构的重要组成要素,涡旋核心线提取和可视化对于切入中尺度涡三维结构研究、开展海洋物质能量垂直运输分析具有重要意义。本文基于客观参考框架和准则,提出了使用客观化的流场参数得到区域,并提取其山谷线作为涡旋核心线的方法,实现了对海洋三维结构中尺度涡旋核心线的提取和可视化。首先,引入了最优局部参考系,使速度、速度梯度等测度转换为在运动的参考系下保持不变的客观量,提升了在海洋科学实践中的可靠性和实用性。其次,针对含有垂向速度的海洋三维流场数据,计算其空间雅可比矩阵,展示了涡旋核心区域的三维结构,实现了海洋涡旋研究从二维到三维的提升。最后,分别在多个半径大小播撒流线种子点,分析不同旋转方向的涡旋,对已提取的涡旋核心线实验结果进行验证,证明了客观海洋三维涡旋核心线提取方法的有效性及可行性。     

幂律流体偏心环空螺旋流轴向速度分布规律的PIV实验研究

为研究幂律流体偏心环空螺旋流轴向速度的分布规律,建立了幂律流体的偏心垂直环空螺旋流动实验模型及实验方案。利用PIV系统拍摄各种工况下偏心环空螺旋流场中粒子的图像,对实验数据进行处理分析,确定了幂律流体偏心环空螺旋流轴向速度分布的主要影响因素。着重分析了压力梯度、黏滞性、内管旋转角速度、偏心度这些影响因素对偏心环空宽、窄间隙处轴向速度分布的影响规律,为解决石油工业中相关问题提供了有益的参考。     

偏心环空螺旋流的PIV测试研究

首次将PIV技术应用于测试偏心环空幂律流体紊流螺旋流速度场，设计了一套可调偏心度的垂直环空管道实验装置。实验为不同浓度的聚丙烯酰胺水溶液在偏心度分别为40％和80％的垂直环空管道内做螺旋流动的PIV实验。得到了轴向速度影响规律：压力梯度一定时，黏滞性的减小或内管转速的增加将使轴向速度增大；流量一定时，黏滞性的减小或内管转速的增加都将使宽间隙处紊流核心区的轴向速度减小；偏心度的增大可以使紊流核心区轴向速度减小；当其它条件相同时，轴向速度随着压力梯度或流量的增加而增大。PIV实验结果与PHOENICS数模结果吻合良好。证明了PIV技术对幂律流体偏心环空螺旋流速度场进行测试是可行的、有效的。     

On the onset of the instability of a three-dimensional jet in a stratified fluid

The onset of a three-dimensional jet flow in a stratified fluid is studied with the aid of a direct numerical simulation. An initially cylindrical jet with a Gaussian velocity profile is considered in a fluid with stable linear density stratification. The results indicate that, if an initial small perturbation of the velocity field has a wide spectrum, an exponential growth of the isolated quasi-two-dimensional mode occurs and its spectral maximum is shifted toward smaller wave numbers in comparison with the maximum of the helical mode of the instability of a nonstratified jet. The growth rate is proportional to Ri^0.5, where Ri is the global Richardson number. The onset of the instability leads to the formation of the flow’s vortex structure, which consists of a collection of different-polarity quasi-two-dimensional vortices located in a horizontal plane near the longitudinal axis of the jet. At sufficiently long times (Nt > 100, where N is the buoyancy frequency and t is time), the growth of instability reaches the saturation stage and further fluctuations in velocity and density decay under the effect of viscous diffusion. At this stage, the flow becomes self-similar and the time dependences of the transverse and vertical widths of the jet are consistent with the asymptotic behaviors of integral parameters of the flow that are observed experimentally in the far stratified wake. The results suggest that the onset of the instability of a quasitwo-dimensional mode can play the determining role in the dynamics of flow in the far stratified wake.     

Mechanism for enhanced diffusivity in the deep-sea perpetual salt fountain

The mechanism of enhanced diffusivity occurring in the deep-sea perpetual salt fountain has been investigated experimentally and numerically. Some factors which possibly contribute to the enhanced diffusivity were found to be the pipe oscillation with ocean waves and its baffled wall surface. Field experiments in the ocean (Onagawa Bay of Miyagi, Japan) and numerical simulations were performed to study and confirm the dynamics of the flow and heat transport with enhanced diffusivity occurring in upwelling deep-sea water. The agreement between the field experimental data and the numerical solutions of an oscillating-wall boundary condition imposed on the baffled pipe is encouraging, and it indicates the baffled pipe surface subject to the oscillatory motion leads to the enhanced diffusivity. The buoyancy force and then upwelling velocity can be greatly increased by the enhanced diffusivity. The dominant mechanism is the occurrence of complicated vortices and vortex shedding leading to efficient mixing and enhanced diffusion.     

Modelling of flow around a near-bed pipeline with a spoiler

Flow around a pipeline with and without a spoiler near a smooth wall is simulated by solving the Navier–Stokes equations. Finite-difference formulation with a second-order upwind scheme in a curvilinear coordinate system is employed. The influences of the spoiler on hydrodynamic forces, pressure distribution, vortex shedding frequency, velocity profile under the pipe, as well as shear stress on the wall are investigated. The attachment of a spoiler significantly increases drag, root-mean-square (RMS) lift, flow through the gap between the pipe and the wall and shear stress on the seabed around the pipe. The spoiler also generates a non-zero mean downward force on the pipeline, which may enhance the self-burial of the pipeline.     

Inleakage of a vortex-ring bunch onto a plane screen

A method is suggested for simulating axisymmetric laminar or turbulent flows formed during the motion of a vortex-ring bunch of given geometry and circulation toward a plane screen. Earlier, similar problems were simulated with the numerical solution of the Navier-Stokes equations for laminar flows. Turbulent flows have remained unconsidered until now. When a vortex ring approaches the screen, the secondary nonstationary flow is induced near the screen’s surface and this secondary flow causes the formation of the radial boundary layer (provided that air viscosity is taken into account). First, the medium spreads out from the critical point at the screen’s center with the negative pressure gradient along the radial coordinate and then detaches in the region of the positive pressure gradient. This radial wall flow and the corresponding boundary layer are considered in the quasi-stationary approximation. When the boundary layer detaches at successive instances, the flow is replenished with the radially moving secondary vortex rings whose circulations have the sign opposite to that of the circulation of the primary vortex ring. It is the interaction of the primary and secondary vortices that governs process dynamics, which differs substantially from that in the case when the formation of secondary vortices is disregarded. The suggested method is based on the method of discrete vortices (a perfect liquid) and the boundary-layer (laminar or turbulent) theory. During the development of the flow under investigation, the nonstationary ascending flow in the direction perpendicular to the screen’s plane is formed and then this flow decays and dissipates. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer show that the velocity of ascending vortices in the plane of the initial vortex bunch is less than one-tenth of the initial velocity of the descending vortex ring. The boundary layer is introduced into calculations with the sole goal of determining the parameters of the secondary vortex rings formed during boundary-layer detachments. The interaction of the primary and secondary vortices is then considered within the framework of a perfect medium. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer on the screen were correlated with the available data obtained in laboratory experiments for small Reynolds numbers. Qualitative agreement between the simulations and experiments is fairly satisfactory. The simulation for one combination of the circulation and vortex-ring geometry takes at most 10–15 min with the use of an average PC.     

动波浪壁圆柱绕流三维数值模拟

利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。     

Degradation of some mechanical current meter measurements byhigh-frequency mooring or wave motion

Comparison of current measurements from moored (paddle-wheel rotor) Aanderaa current meters and acoustic Doppler current profilers in a strong tidal flow on Georges Bank indicates rate under-reading by the Aanderaa meters at some vertical positions. The under-reading may arise from mooring-line vibrations induced by vortex shedding from spherical buoyancy packages, and shielding of the paddle-wheel rotors due to the meters' inability to remain aligned with the fluctuating relative water velocity. Field tests and a simple model are used to investigate this explanation. The authors briefly review the Georges Bank observations that instigated this investigation. An explanation for the degradation of the Aanderaa measurements is proposed, observational information from field tests is discussed, and the model and its implications are presented. The model results support the hypothesis that high-frequency mooring-line vibration causes the degradation of the Aanderaa measurements     

Vortex dynamics of a fluid near a boundary with a circular cavity

The motion of a point vortex along a rectilinear boundary with a circle cavity, which models the coastline of a bay, and associated fluid particle advection are studied within a model of barotropic inviscid fluid. Using an analytical expression for the complex potential through which the velocity field is determined, we show that fluid particles start moving irregularly when the vortex is passing the cavity due to the nonstationarity of the velocity field generated by the vortex. Some of the fluid particles which were initially inside the vortex atmosphere leave it due to the irregularity and remain within the cavity vicinity. Depending on the initial position of the vortex and a parameter that determines the cavity size, the fraction of these fluid particles can differ significantly from fluid particles initially uniformly distributed within the vortex atmosphere. The escape of fluid particles from the vortex atmosphere is shown to be most efficient in the case of a relatively closed cavity under the condition that the initial vortex atmosphere area should be significantly smaller than the cavity area.     

深海悬垂取水管振动特性研究

深海悬垂取水管的设计趋向于大直径、高内流流速,内流对取水管振动特性的影响不可忽略,分析内部流动引起的管道动态失稳行为以及评估失稳临界流速具有重要工程意义。基于小尺度模型试验,研究内流对深海悬垂取水管道振动特性的影响规律,结果表明：随着内流流速增加,取水管模型会发生动态失稳行为,表现为一种间歇性的近周期运动,该运动主要由管道一阶弯曲模态引起的不稳定性诱导。管道材料、顶端连接方式及底部配重块均影响其动态稳定性,顶端固定连接时质量比小的管道更易发生动态失稳;顶端变为铰接时,质量比大的管道失稳临界流速变小,质量比小的管道失稳临界流速变大。相比顶端边界条件,底部配重块对管道失稳临界流速影响不显著但可以减小其振动幅值。     

基于重叠网格方法浮筒涡激运动特性数值分析

浮筒被广泛应用在海洋工程中,研究浮筒的涡激运动对于减少其对海洋平台构件的疲劳损坏具有指导意义。传统的动网格方法在处理浮筒转动运动时会因网格变形过大导致计算不收敛,采用了重叠网格方法以解决这一问题。数值试验采用了基于开源工具包Open FOAM自主开发的naoe-FOAM-SJTU求解器。分别进行了自由衰减数值试验和涡激运动数值试验。研究表明,随着折合速度的增加,浮筒的顺流向、横流向、垂荡和艏摇运动频率增加,且顺流向与垂荡频率相近,横流向与艏摇频率相近;其次,根据涡量场分布,表明浮筒前一时刻的泻涡会影响到浮筒之后时刻的周向涡量分布;最后,研究发现自由端对于浮筒尾流场泻涡有着显著影响,为将来探究减少浮筒涡激运动的方法提供指导。     

2012年“苏拉”和“达维”双台风影响的近海风暴潮过程

作为影响我国沿海的主要自然灾害之一,台风风暴潮的产生和影响机制与防灾减灾息息相关。双台风引起的风暴潮因台风强度、路径等相对关系复杂多变,目前双台风相互作用下的风暴潮研究还不充分。采用参数化台风模型对2012年典型双台风"苏拉"和"达维"的风场、气压场过程进行了模拟与融合,并采用ELCIRC模型对双台风作用下的风暴潮过程进行了模拟。引入单台风单独作用的假设算例,探讨了双台风之间对增水、流场的相互影响和影响区域。研究结果表明,虽然两个台风登陆强度相当,但台风"达维"在海州湾海域引起的增水要远大于"苏拉"在台湾、福建海域引起的增水。风暴潮引起的增水及流速变化与台风在海表的风应力密切相关,较大变化幅值分布在台风行进路径的右侧。与台风单独作用时相比,台风"苏拉"与"达维"引起的风暴潮增水与流速变化在两者相互作用下均有所削弱,其中"苏拉"引起的风暴潮受到的影响更大。双台风风暴潮之间的非线性效应在不同区域的强度存在差异,在台风"苏拉"主要影响区域内非线性效应较强,其他区域则相对较弱。以上结果表明产生风暴潮较弱的台风一方对气象环境敏感性更高,风暴潮的响应更显著。     

Hydrodynamic coefficients for a thick-walled bottomless cylindrical body floating in water of finite depth

The paper deals with the linearized hydrodynamic forces acting on a thick-walled, bottomless cylindrical body having vertical symmetry axis and oscillating in water of finite depth. For the solution of the radiation problem, the flow field around the structure is subdivided into ring-shaped fluid regions, in each of which an axisymmetric eigenfunction expansion for the velocity potential is made. By implementing Galerkin's method the various potential solutions are then matched and numerical results concerning the hydrodynamic coefficients for heave, surge and pitch motions, as well as the coupling terms between the last two modes are obtained.     

In-line force on a cylinder translating in oscillatory flow

Experiments were conducted with smooth and sand-roughened cylinders moving with constant velocity in a sinusoidally oscillating flow to determine the drag and inertia coefficients and to examine the effect of wake biasing on the modified Morison equation. The various flow parameters such as the relative cylinder velocity. Reynolds number, and the Keulegan-Carpenter number were varied systematically and the in-line force measured simultaneously. The principal results, equally valid for both smooth and rough cylinders, are as follows: the drag coefficient decreases with increasing relative current for a given Reynolds number and Keulegan-Carpenter number; the effect of wake biasing on the drag and inertia coefficients is most pronounced in the drag-inertia dominated regime; and the two-term Morison equation with force coefficients obtained under no-current conditions is not applicable to the prediction of wave and current induced loads on circular cylinders.     

A three-dimensional wake impingement model and applications on tandem oscillating foils

Potential flow based vortex numerical methods have been widely used in aerodynamics and hydrodynamics. In these methods, vortices shed from lifting bodies are traced by using vortex filaments or dipole panels. When the wake elements encounter a downstream body, such as a rudder behind a propeller or a stator behind a rotor, a treatment is necessary to divert the wake elements to pass by the body. This treatment is vital to make wake simulations realistic and to satisfy the non-penetration condition during wake body interaction. It also helps to avoid pure numerical disturbances such as when a vortex filament or an edge of a dipole panel passes through the collection point of a body element; this is a singularity for induced velocity and it will introduce a large numerical disturbance. This necessary treatment for three-dimensional problems with geometrical complexity has not been found to date. In this study, a wake impingement model was developed to divert wake elements to slip over the body surface, model the vortex/body interaction, and predict forces on fluctuating components. The model was also tested on configurations of oscillating foils in tandem with an existing panel method code. Simulation results with the wake impingement model are shown to be in closer agreement with limited published experimental data than those without the model. With the established wake impingement model, force fluctuations on the after body due to the wake vortex impingement were investigated based on a series of simulations. The series varied several parameters including distance between two foils, oscillating frequency, span, rear foil pitch angle, swap angle and vertical position.     

Cyclostrophic adjustment and nonlinear oscillations in the core of an intense atmospheric vortex

Intense atmospheric vortices are characterized by a regime of cyclostrophic balance, i.e., the balance between the pressure gradient and centrifugal force. To describe motions in the core of an axisymmetrical vortex, a class of exact solutions to the equations of gas dynamics with a linear dependence on radius is derived for the velocity components and with a quadratic dependence for temperature. It is shown that small deviations from the balance state give rise to oscillations of the hydrothermodynamic fields in the vortex core with a frequency proportional to the angular velocity of the rotation of the core. For fairly large initial deviations, oscillations are clearly anharmonic and, under the conditions of the prevailing centrifugal force, result in a significant temperature decrease on the vortex axis. The application of this class of solutions to describing the Ranque vortex effect (the intense cooling of gas during rapid rotations) and the acoustic radiation from tornadoes is discussed.     

Vortex shedding suppression for flow over a circular cylinder near a plane boundary

In this study, the Navier-Stokes equations and the pressure Poisson equation for two-dimensional time-dependent viscous flows are solved with a finite difference method in a curvilinear coordinate system. With this numerical procedure, the vortex shedding flow past a circular cylinder near a wall is investigated. The flow is calculated for a broad range of gap ratios for different Reynolds numbers ranging from 80 to 1000. Based on the numerical solutions, the vortex shedding is observed using various methods, and the mechanism for the vortex shedding suppression at small gap ratios is analyzed. The critical gap ratio at which the vortex shedding is suppressed is identified at different Reynolds numbers.