自由面对潜艇尾流场流动特性影响研究

潜艇周围绕流场流动特性会影响潜艇的机动性能, 特别是近水面航行时, 自由面的存在会增大潜艇尾流场的复杂程度.为探究潜艇在近水面航行时自由液面对潜艇尾流场流动特性的影响机理,借助大型水下三维粒子图像测速技术开展潜艇尾流场流动特性研究.首先通过美国泰勒水池标准模型实验结果对试验方法准确性进行验证; 随后,用验证后的模型试验方法对潜艇尾流场进行测量,得到不同潜深工况、不同速度下的桨盘面轴向速度以及脉动速度,同时辅以数值模拟对试验无法测得的兴波波系及中纵剖面速度场加以补充,从兴波角度阐述了自由面对潜艇尾流场流动特性影响机理. 研究结果表明:潜艇在近水面航行时, 随着Fr增大,桨盘面处轴向速度云图中上方等值线整体趋于扁平化, 较4D潜深工况,1.5D潜深工况出现局部脉动速度极大值, 且脉动速度结构整体下移; 自由面存在时,艇体与自由面间流场速度明显增大, 特别在桨盘面区域, 流场速度明显提升.随着Fr增大, 桨盘面处的自由液面高度逐渐降低,这就导致了桨盘面位置出现更大的流体速度, 即造成了桨盘面伴流场挤压现象.      

浅水舰船引起水底附近的扰动速度研究

基于浅水波浪理论,建立了亚临界和超临界航速条件下舰船运动引起水底附近扰动速度的理论计算模型。通过数值计算,揭示了扰动速度分量的变化规律和分布特点。建立了舰船扰动速度的测试系统,对船模运动引起的扰动速度进行了测量。横向扰动速度分量计算结果与实验结果基本吻合,纵向扰动速度分量计算结果与实验结果定性一致。理论与实验研究表明,...     

通气空泡尾部气泡流仿真和特性分析

基于欧拉-欧拉双流体模型对通气空泡尾部气泡流进行数值仿真, 并采用基于多尺寸分组模型的总体平衡方法预估气泡尺寸分布. 应用改进后湍流耗散系数计算模型, 考虑了气泡体积含量对在湍流作用下气泡扩散现象的影响. 基于上述模型对两种试验工况下流场进行了数值仿真.结果表明模型对空泡尾部回流区特性进行了准确预示, 在回流区高湍流度作用下气泡迅速破碎成小气泡. 并进一步得到试验体尾流区空泡体积分数和速度分布. 尾流区水流速度分布保持了流体经过非流线型对称体时产生的尾流分布规律. 仿真结果与试验数据相一致, 模型适用性得到验证.      

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

基于Fluent平衡大气边界层流场模拟及应用

采用Fluent进行平衡大气边界层流场模拟,在地面边界引入气动粗糙壁面函数,上边界分别采用剪力驱动和速度驱动。首先对平坦区域进行流场模拟,然后基于模拟的流场对立方体建筑物表面风压分布进行了计算。结果表明:基于本文的边界条件设置方法,三种k-ε湍流模型均能获得良好的自保持性平衡大气边界层流场;在地面边界引入气动粗糙壁面函数,可实现任意类型的平衡大气边界层流场模拟,且对近地单元的高度无限制,有利于精细化网格计算;在上边界采用剪力驱动和速度驱动是等效的,采用速度驱动边界条件有利于提高计算效率;三种k-ε湍流模型计算的立方体中心面表面压力分布基本在测量值变化范围内,说明基于本文的平衡大气边界层流场可应用于结构风荷载数值模拟研究中。     

定床弯道内水沙两相运动的数值模拟

在适体同位网格中采用非正交曲线坐标系下的三维k-ε-kp固液两相双流体湍流模型研究弯道内水流和悬浮泥沙运动,主要计算了试验室S型水槽内清水流动的三维流场、120°弯道内水沙两相流动中底沙与底流的运动轨迹以及S型水槽内水沙两相流动的两相流场和泥沙浓度场. 对于S型水槽内清水流动,数值结果与试验结果吻合良好. 120°弯道内水沙两相流动中固液两相的运动轨迹在弯道直线段基本重合,在弯道内泥沙轨迹逐步偏离水体轨迹,其偏离程度随泥沙粒径增大而增大. 从S型水槽内水沙两相流动计算结果中发现泥沙纵向流速在壁面附近比水流纵向速度大,在远离壁面区域比水流纵向速度小;弯道内泥沙横向流速比水流横向流速小;垂向流速在直线段和泥沙沉速相当,在弯道内受螺旋水流影响而变化;两相流速差别随泥沙粒径增大而变大;泥沙浓度呈现下浓上稀的分布,在弯道内横向断面上呈现凸岸大凹岸小的分布,泥沙浓度随泥沙粒径增大而减小.      

自由面对潜艇尾流场流动特性影响研究

潜艇周围绕流场流动特性会影响潜艇的机动性能,特别是近水面航行时,自由面的存在会增大潜艇尾流场的复杂程度.为探究潜艇在近水面航行时自由液面对潜艇尾流场流动特性的影响机理,借助大型水下三维粒子图像测速技术开展潜艇尾流场流动特性研究.首先通过美国泰勒水池标准模型实验结果对试验方法准确性进行验证;随后,用验证后的模型试验方法对...     

逆流式冷却塔填料及淋水分布的数值优化设计

在逆流式自然通风冷却塔流场、温度场和湿度场数值模拟工作的基础上 ,采用数值分析的方法研究了非均匀填料分布及非均匀淋水密度分布对提高冷却效率的影响 ,结合某电厂的工程设计实例 ,给出了最佳的填料层分布及淋水密度分布     

逆流式冷却塔填料及淋水分布的数值优化设计

在逆流式自然通风冷却塔流场、温度场和湿度场数值模拟工作的基础上,采用数值分析的方法研究了非均匀填料分布及非均匀淋水密度分布对提高冷却效率的影响,结合某电厂的工程设计实例,给出了最佳的填料层分布及淋水密度分布。     

方柱尾流流态对床面冲淤形态的影响

应用平面二维悬沙数学模型对方柱尾流区的泥沙输运及床面调整进行了数值模拟，并对重要物理参数做了分析说明。数值计算采用时间分裂一全隐式有限差分格式，流场的计算基于沿水深平均的RANS方程。通过系统的数值模拟，揭示了不同来流情况下，方柱尾流区的流态变化，及与之相应的床面变化规律，并对不同流态下的泥沙运动作了分析。计算结果表明在方柱尾流区，不同的流场流态对尾流区的床面调整有明显影响。流场较弱的情况下，尾流区中湍动强度相应较弱，此时悬浮泥沙由于流速减小而普遍落淤，床面应力的减小也致使床面冲刷量减小。随着流场强度的逐渐增大，尾流区中湍动强度相应增强，床面应力增大，同时湍流的强烈交换作用增强了对泥沙的输运作用，床面变化不再是简单的淤积状态，部分区域出现了冲刷。     

PIV analysis of flow around a container ship model with a rotating propeller

The flow characteristics of the propeller wake behind a container ship model with a rotating propeller were investigated using a two-frame PIV (Particle Image Velocimetry) technique. Ensemble-averaged mean velocity fields were measured at four different blade phases and ensemble-averaged to investigate the flow structure in the near-wake region. The mean velocity fields in longitudinal planes show that a velocity deficit is formed in the regions near the blade tips and hub. As the flow develops in the downstream direction, the trailing vortices formed behind the propeller hub move upward slightly due to the presence of the hull wake and free surface. Interaction between the bilge vortices and the incoming flow around the hull causes the flow structure to be asymmetric. Contour plots of the vorticity give information on the radial distribution of the loading on the blades. The radial velocity profiles fluctuate to a greater extent under the heavy (J=0.59) and light loading (J=0.88) conditions than under the design loading condition (J=0.72). The turbulence intensity has large values around the tip and trailing vortices. As the wake develops in the downstream direction, the strength of the vorticity diminishes and the turbulence intensity increases due to turbulent diffusion and active mixing between the tip vortices and the adjacent wake flow.     

The interaction of background ocean air bubbles with a surface ship

A two-fluid model suitable for the calculation of the two-phase flow field around a naval surface ship is presented. This model couples the Reynolds-averaged Navier–Stokes (RANS) equations with equations for the evolution of the gas-phase momentum, volume fraction and bubble number density, thereby allowing the multidimensional calculation of the two-phase flow for monodisperse variable size bubbles. The bubble field modifies the liquid solution through changes in the liquid mass and momentum conservation equations. The model is applied to the case of the scavenging of wind-induced sea-background bubbles by an unpropelled US Navy frigate under non-zero Froude number boundary conditions at the free surface. This is an important test case, because it can be simulated experimentally with a model-scale ship in a towing tank. A significant modification of the background bubble field is predicted in the wake of the ship, where bubble depletion occurs along with a reduction in the bubble size due to dissolution. This effect is due to lateral phase distribution phenomena and the generation of an upwelling plume in the near wake that brings smaller bubbles up to the surface. © 1998 John Wiley & Sons, Ltd.     

横摆角下汽车尾部湍流特征量的PIV试验分析

汽车尾部湍流场是汽车压差阻力的主要来源,在HD-2汽车模型风洞中,首先使用测力天平和测压系统,对横摆角工况下汽车模型的气动六分力和纵对称截面48个测点的表面压力进行了测量,然后利用PIV测量技术对模型在横摆角分别为0°、15°的尾部湍流场进行了测量,获得该模型尾流场的速度场、涡量场和雷诺应力流场信息,通过计算得出尾流场区域空间相关系数和湍流积分尺度。结果表明:在横摆角工况下,汽车模型尾部涡流的结构呈现向上发展的趋势;尾流场拖拽涡的范围和强度的增大导致了模型气动力出现较大的增加;湍流积分尺度的变化表明,尾部涡流区的分离噪声与涡流分离位置有关,在汽车尾部造型设计中,要尽量推迟尾部涡流的分离。     

Measurement of tonal-noise characteristics and periodic flow structure around NACA0018 airfoil

The characteristics of tonal noise and the variations of flow structure around NACA0018 airfoil in a uniform flow are studied by means of simultaneous measurement of noise and velocity field by particle-image velocimetry to understand the generation mechanism of tonal noise. Measurements are made on the noise characteristics, the phase-averaged velocity field with respect to the noise signal, and the cross-correlation contour of velocity fluctuations and noise signal. These experimental results indicate that the tonal noise is generated from the periodic vortex structure on the pressure surface of the airfoil near the trailing edge of the airfoil. It is found that the vortex structure is highly correlated with the noise signal, which indicates the presence of noise-source distribution on the pressure surface. The vorticity distribution on the pressure surface breaks down near the trailing edge of the airfoil and forms a staggered vortex street in the wake of the airfoil.     

Experimental study of non-uniform inlet conditions and three-dimensional effects of vertical air–water two-phase flow in a narrow rectangular duct

To study the three-dimensional interfacial structure development in vertical two-phase flow, air–water upflow experiments were performed in a rectangular duct. Various non-uniform two-phase profiles were created by injecting air from individually controlled spargers at the duct inlet into uniformly injected water flow. A four-sensor conductivity probe was used to measure local void fraction, interfacial area concentration, bubble velocity and Sauter mean diameter at three axial locations to record the development of two-phase parameters. Experimental results showed that the lateral development across the wider dimension of the duct was significant with a non-uniform inlet profile when compared to a uniform inlet profile. It is postulated that lift, wall and turbulent forces are the major contributors to the lateral distribution of the two-phase interfacial structures making this an useful experiment for benchmarking three-dimensional two-fluid models. In examining the interfacial area, the shearing-off of group 1 bubbles (defined as the smaller spherical and distorted bubbles) from the skirt region of group 2 bubbles (defined as the bigger cap and churn bubbles), the coalescence of group 2 bubbles due to wake entrainment, and random collision are the major source and sink mechanisms of interfacial area concentration.     

Mean velocity and moments of turbulent velocity fluctuations in the wake of a model ship propulsor

Pod drives are modern outboard ship propulsion systems with a motor encapsulated in a watertight pod, whose shaft is connected directly to one or two propellers. The whole unit hangs from the stern of the ship and rotates azimuthally, thus providing thrust and steering without the need of a rudder. Force/momentum and phase-resolved laser Doppler anemometry (LDA) measurements were performed for in line co-rotating and contra-rotating propellers pod drive models. The measurements permitted to characterize these ship propulsion systems in terms of their hydrodynamic characteristics. The torque delivered to the propellers and the thrust of the system were measured for different operation conditions of the propellers. These measurements lead to the hydrodynamic optimization of the ship propulsion system. The parameters under focus revealed the influence of distance between propeller planes, propeller frequency of rotation ratio and type of propellers (co- or contra-rotating) on the overall efficiency of the system. Two of the ship propulsion systems under consideration were chosen, based on their hydrodynamic characteristics, for a detailed study of the swirling wake flow by means of laser Doppler anemometry. A two-component laser Doppler system was employed for the velocity measurements. A light barrier mounted on the axle of the rear propeller motor supplied a TTL signal to mark the beginning of each period, thus providing angle information for the LDA measurements. Measurements were conducted for four axial positions in the slipstream of the pod drive models. The results show that the wake of contra-rotating propeller is more homogeneous than when they co-rotate. In agreement with the results of the force/momentum measurements and with hypotheses put forward in the literature (see e.g. Poehls in Entwurfsgrundlagen für Schraubenpropeller, 1984; Schneekluth in Hydromechanik zum Schiffsentwurf, 1988; Breslin and Andersen in Hydrodynamics of ship propellers, 1996; Schneekluth and Bertram in Ship design for efficiency and economy, 1998), the co-rotating propellers model showed a much stronger swirl in the wake of the propulsor. The anisotropy of turbulence was analyzed using the anisotropy tensor introduced by Lumley and Newman (J Fluid Mech 82(1):161–178, 1977). The invariants of the anisotropy tensor of the wake flow were computed and were plotted in the Lumley–Newman-diagram. These measurements revealed that the anisotropy tensor in the wake of ship propellers is located near to the borders of the invariant map, showing a large degree of anisotropy. They will be presented and will be discussed with respect to applications of turbulence models to predict swirling flows.     

Criteria of negative wake generation behind a cylinder

It was demonstrated by simulation in our previous study that both the normal stress and its gradient are responsible for the negative wake generation (overshoot in the axial velocity) and streamline shifting. Extensional properties of the fluids dominate the generation of the negative wake, while other factors strengthen or weaken the formation of velocity overshoot. In this study, the criteria for the negative wake generation are discussed in detail for various fluid models, including the PTT, the FENE-CR, the FENE-P, and the Giesekus models. With the FENE-CR fluid, it is easier to generate negative wake than with the FENE-P fluid. This confirms that the constant shear viscosity FENE-CR fluid enhances the velocity overshoot, and that the shear-thinning viscosity FENE-P fluid delays the negative wake generation. The Giesekus fluid has a similar behaviour to the PTT fluid with regarding to the critical conditions of negative wake generation when appropriate fluid parameters are selected. The mechanism of wall proximity in enhancing the negative wake generation is also demonstrated with the analysis for the first time.     

Measurements of a stator-induced circumferentially varying flow

The flow physics associated with the generation of both axisymmetric and non-axisymmetric swirl by various deflection patterns of a stator array was investigated experimentally through surface pressure and Stereoscopic Particle Image Velocimetry measurements. A three-dimensional rendering technique was developed to reconstruct the flow field around the model and in its wake. The three-dimensional fluid volume was reconstructed from multiple two-dimensional measurement planes. A cyclic distribution of the stators’ deflections resulted in non-axisymmetric distributions of the surface pressure and the flow field downstream of the stator array. The addition of a shroud had an amplifying effect: accelerating the flow through the stator array while reducing the non-uniform tangential velocity component generated by the stators. In the model near wake the flow field is associated with secondary flow patterns in the form of coherent streamwise vortical structures that can be described by potential flow mechanisms. The collective pitch distribution of the stators produces a flow field that resembles a potential Rankine vortex, whereas the cyclic pitch distribution generates a flow pattern that can be described by a potential vortex pair in a cross-flow.     

Calculation of flow around a lifting surface in the case of large deformations

Many studies have been made of the nonstationary flow of an ideal incompressible fluid around a lifting surface. The present state of the numerical methods of solution of this problem is reviewed in [1]. The present paper studies three-dimensional nonstationary flow around a lifting surface which undergoes deformation and behind which a wake vortex surface is formed. The lifting and wake vortex surfaces are represented in parametric form. The metrics of these surfaces are used, and the introduced vortex function is approximated by bicubic splines. For the convenient application of the theory developed here to the flapping flight of insects, for which it is sometimes difficult to distinguish the lateral and trailing edges of the wings, the following terminology is introduced. The part of the edge of the lifting surface from which the wake vortex surface is shed is called the trailing edge. The remaining part is called the leading edge. On the leading edge, the velocity has a singularity. Test calculations have demonstrated the effectiveness of the method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 72–79, July–August, 1980.