Numerical investigation of turbulence near a sheared air–water interface. Part 1: Turbulence statistics beneath a flat water surface sheared by wind

Direct numerical simulation (DNS) and large eddy simulation (LES) of turbulent shear flow beneath a flat water surface with imposed wind shear stress are presented. The results of DNS indicate that there are clear differences, and also similarities, between wind-driven flow and the flow near a solid wall. The qualitative structures of turbulence are similar in both types of flow. Low-speed streaks are also present in wind-driven flow. On the other hand, profiles of the mean velocity and turbulence intensities in wind-driven flow are significantly different from those in flow near a solid wall. The differences are attributed to the lack of a viscous sublayer, and to the boundary condition which allows fluctuations of the tangential velocity components at the boundary. LES of the same flow was also carried out to evaluate subgrid-scale models. It was shown that the features of the flow observed in DNS are well simulated by all models tested, and that the discrepancies between DNS and LES with the dynamic mixed model are very small. Received: August 17, 2000 / Accepted: December 22, 2000     

大涡模拟在河道平面二维水流模拟中的初步应用

通过丁坝对水流进行控制是航道整治中常用的工程措施,丁坝附近水流呈强紊动特性,大涡模拟相对于基于雷诺方程的时均模型对涡旋有较强的捕捉能力。在直角坐标系的基础上,将大涡模拟中的亚格子应力模型(SGS模型)引入河道水流平面二维数学模型中,通过盒式滤波函数将控制方程进行滤波,大尺度量通过控制方程直接求解,小尺度量借助Smagorinsky提出的亚格子尺度应力模型进行求解。采用空间等步长的正方形网格和交替方向隐式差分法对其控制方程进行离散求解。将模型初步应用于非淹没丁坝绕流的数值模拟中,计算结果表明,该模型对湍流旋涡的捕捉和水深计算较为合理,计算精度可满足工程实践要求。     

航道清礁对西江鱼类保护区紊流体的影响

河底裸露的礁石增加水体紊动,有助于鱼卵正常形成,大量存在礁石的河段往往是鱼类栖息和产卵的重要场所。礁石也存在碍航的可能,航道清礁一定程度上影响水体紊动强度,改变鱼类栖息环境。为此,在重要鱼类栖息河段进行航道清礁需充分论证。以紊动动能作为评价水体紊动强度指标,开展西江中游珍稀鱼类保护区河段三维紊流数值模拟,计算并统计航道清礁工程前后水体紊动动能及不同强度的紊流水体体积,分析其变化规律及原因。研究结果表明,断面中紊动动能强度与水深成正相关关系;清礁工程导致河段中弱紊动动能的紊流水体体积增加,强紊动动能的紊流体积减小;底层水体紊动动能的变化量高于表层;强紊动动能水体体积减少的变化率大于弱紊动动能水体体积增加的变化率。     

壁面展向周期振动减阻机理及其二维波谱分析

利用直接数值模拟(DNS)对壁面作展向周期运动的槽道湍流进行研究,建立了槽道湍流数据库.通过改变振幅大小和振动周期,可以使壁面摩擦阻力明显减少.随着平均减阻率的增加,壁面阻力随时间的变化也更加稳定并呈现出较大周期的变化.将一个典型的阻力变化周期分成三个不同区域进行讨论,对湍流脉动能的二维波谱进行了定量分析,进一步揭示出壁面展向周期振动抑制湍流、实现减阻的内在机理.     

熵格子Boltzmann方法的重整化群代数湍流模型

重整化群理论所建立的湍流模型能够最大程度地减小模型经验性,因此文章尝试将重整化群代数湍流模型引入到熵格子Boltzmann方法中,建立新型的计算模型以对高雷诺数湍流进行模拟研究。同时为了进行比较研究,还建立了熵格子Boltzmann方法的标准大涡模拟模型。完成了对高雷诺数湍流绕流场的模拟计算。结果表明:所建立的熵格子Boltzmann方法重整化群代数湍流模型能够有效地模拟高雷诺数湍流流动问题;其对紧贴壁面处较小尺度湍涡的模拟结果趋近于大涡模拟的结果;重整化群代数湍流模型在对高雷诺数湍流的模拟中表现出耗散模型的特征。     

Large eddy simulations for quasi-2D turbulence in shallow flows: A comparison between different subgrid scale models

In this study, the performance of the horizontal large eddy simulation module, developed at the University of Leuven (HLES-KULeuven module) is assessed. A comparison between different subgrid scale models has been carried out. The study is concerned with the non-rotating and unstratified flows. The results of the simulation for an oscillatory backward facing (BFS) flow are presented in case of an expanding flume based on a one-length scale approach and a two-length scale approach. Three subgrid scale (SGS) models have been tested: Smagorinsky SGS model (Smagorinsky, J., (1963). General circulation experiments with the primitive equations, I. the basic experiments. Monthly Weather Review, 91(3), 99–164), Uittenbogaard SGS model (Uittenbogaard, R.E., and van Vossen, B., (2004). Subgrid-scale model for quasi-2D turbulence in shallow water. Shallow Flows. Jirka and Uijttewaal (Eds.), Taylor & Francis Group, London, ISBN 90 5809 700 5) and a proposed two-length scale approach. The first two models are considered to be a one-length scale models. A simulation without a subgrid scale model for the horizontal mixing has also been conducted. In all simulations, a quadratic friction model parameterizes the dissipation produced by the 3D-subdepth scale turbulence. The two-length scale concept uses a newly mixing length formulation for the quasi-2D turbulence and doesn't depend on the filter width in contrast to the one-length scale approach, in which the mixing length is function of the filter width. The outputs of the HLES-KULeuven module have been compared with the experimental data taken from Stelling, G.S., and Wang, L.X., (1984). Experiments and computations on separating flow in an expanding flume. Dept. Civil Engineering, Delft University of Technology, Report 2–84.). The two-length scale approach has been validated with experimental data from SERC Flood Channel Facility at HR Wallingford. In general, there is a qualitative agreement with the experimental data. It has also been found that the two-length scale approach produces more elongated and less isotropic vortex than the one-length scale models.     

非均衡湍流边界层脉动压力测试研究

通过建立非均衡湍流边界层脉动压力测试技术及可提供较高测试信噪比的低噪声风洞试验装置,采用表面贴装式硅微传声器线列阵,测量获得了不同压力梯度模型的湍流边界层脉动压力特性,压力梯度模型表面湍流边界层脉动压力低频平台区幅值增加3-5 dB;回归得到了非均衡湍流边界层流脉动压力频谱与Strauhal数及Re数和压力梯度相关的拟合模型,具有一定的适用性;压力梯度对湍流边界层脉动压力频率—波数谱的影响主要集中在300 Hz以下的低频段,传输波数附近峰脊区谱级增加3-5 dB。     

Effect of bed vicinity on vortex shedding and force coefficients of fluid flow on an offshore pipeline

The effect of rigid bed proximity on flow parameters and hydrodynamic loads in offshore pipelines exposed to turbulent flow is investigated numerically. The Galerkin finite volume method is employed to solve the unsteady incompressible 2D Navier–Stokes equations. The large eddy simulation turbulence model is solved using the artificial compressibility method and dual time-stepping approach. The proposed algorithm is developed for a wide range of turbulent flows with Reynolds numbers of 9500 to 1.5×10⁴. Evaluation of the developed numerical model shows that the proposed technique is capable of properly predicting hydrodynamic forces and simulating the flow pattern. The obtained results show that the lift and drag coefficients are strongly affected by the gap ratio. The mean drag coefficient slightly increases as the gap ratio increases, although the mean lift coefficient rapidly decreases. The vortex shedding suppression happen at the gap ratio of less than 0.2.     

Using high sampling-rate ADCP for observing vigorous processes above sloping [deep] ocean bottoms

Above sloping bottoms in the ocean mixing processes are not predominantly generated by shear-induced turbulence via bottom friction. Instead, the restratifying buoyancy forces and internal waves create a highly non-linearly varying environment including ‘stratified turbulence’. Most of the resulting vigorous mixing processes that dominate sediment resuspension occur during the passage of frontal bores or solitary boluses, ‘solibores’. Here, the observed evolution of different forms of highly non-linear strictly upslope moving ‘waves’, bores or boluses are reviewed from various NIOZ projects at deep sloping bottom sites ranging from 500 to 3000 m.Such fronts pass a fixed site within a few minutes, extending some 60 ± 30 m above the bottom and occurring over much larger periods at once per subinertial or meso-scale period or approximately, but not exactly, once per tidal harmonic period. In order to observe the details of such solibore one needs specific, high-sampling rate equipment. A suitable piece of equipment is a bottom-mounted 4-beam 300 kHz acoustic Doppler current profiler (ADCP), provided it samples at a rate of about once per second over a period of at least several weeks. Not just the three components of current velocity [u,v,w] are monitored over a range of some 80 m at 1 m intervals, but also the relative ‘echo intensity’ dI, which is a measure for suspended matter and stratified turbulence. Such ADCP-observations are combined and compared with high-resolution temperature measurements. Fine details show a turbulent inner core with more or less laminar streamlines outside it. Whether a front or a bolus, the bore is never observed as a completely closed contour, as swept up turbulent material is sucked into the core at the rear end.     

不同坝型丁坝坝体周围水流紊动特性试验研究*

通过水槽试验,结合理论分析,研究丁坝附近水流紊动强度沿水流方向、水深方向和横向的分布规律。绘制5种不同坝体结构形式其周围水流脉动动能和紊动强度的分布等值线图并分析分布规律,进一步探讨不同坝体挑角对丁坝附近水流紊动强度和脉动动能的影响。     

船坞竖井闸阀输水廊道水力特性的数值模拟研究

近年来,竖井闸阀输水廊道已广泛应用于新建大型修造船坞中。结合某竖井闸阀输水廊道的模型试验,运用FLUENT软件对廊道的水力特性进行三维数值模拟,研究廊道中水流的流速、压强、湍流动能和湍流耗散率等参数的分布规律和输水廊道的安全性。通过对消能段湍流参数的分析发现,在进水井中水流形成的漩涡和水流与壁面的撞击是消耗水流动能主因。不同水位条件下廊道内水流压力的数值模拟结果与模型试验结果吻合性较好,进一步验证了数值模拟计算的可靠性和采用数值模拟方法辅助优化输水廊道设计的可行性。     

Swimming for survival: A role of phytoplankton motility in a stratified turbulent environment

We investigate a role for vertical migration in stratified coastal water, where the swimming speed is generally significantly less than the typical turbulent fluctuations in a tidally-mixed bottom layer. In our modelling approach we use a k- turbulence model to describe the physical forcing, a Lagrangian random walk model to describe the vertical displacement of individual cells in response to turbulence and due to cell motility, and a phytoplankton growth model to direct the swimming behaviour of the phytoplankton according to their light and nutrient requirements. The model results show how the cells form a stable subsurface chlorophyll maximum (SCM) at the base of the thermocline where episodic tidal turbulence causes erosion of part of the SCM biomass into the bottom mixed layer (BML). We then focus on the question of whether an ability to swim (weakly, compared to typical bottom layer turbulent intensities) provides any advantage by allowing return to the SCM. Our results show that tidal turbulence in the BML helps both motile and neutrally-buoyant cells by periodically pushing them into the base of the thermocline. Motile cells then have the advantage that they can swim further into the thermocline towards higher light which also reduces the likelihood of being re-mixed back into the BML.     

湍流中螺旋桨激振力宽频谱及参数影响研究

针对螺旋桨在湍流中的非定常激振力问题,文章介绍并推导了相关性法,在相关性法的基础上,重新定义了湍流速度相关性公式,并考虑了螺旋桨抽吸作用导致的湍流涡拉伸变形的影响,引入了拉伸系数ax,ay和az,推导了新的螺旋桨激振力宽频公式。同时,使用了十桨叶螺旋桨模型,分别采用相关性法和改进后的相关性法对其进行了计算并与试验值进行对比,发现改进后的相关性法能较好地预报螺旋桨激振力宽频谱,特别是在叶频处能与试验值符合较好。最后,在此基础上,研究了螺旋桨转速n、湍流积分尺度Λ和螺旋桨桨叶数N等参数对螺旋桨激振力宽频谱的影响。     

一种河道紊流三维数值模拟离散格式的应用探讨

通过引用气体动力学中的Osher格式计算控制体交界面通量,建立了一种基于有限体积法的河道紊流三维数学模型。该格式是对准确的非线性黎曼问题的近似解法,具有很好的守恒性、逆风性、单调性和高分辨率捕捉间断的能力,解决了河床变化剧烈处水流计算易发散的难题,在紊流模型的选取上,考虑到紊动能量和紊动尺度在水平方向和垂直方向上的差异,平面上采用大涡模型,垂向上采用k-ε模型。通过桥墩和丁坝局部水流实测资料校核模型模拟局部复杂水流结构的能力,模型计算得到的马蹄涡、回流结构及流速分布与实测基本一致。     

绕船体自由面周围三维粘性流场的数值模拟

本文采用有限体积法通过求解不可压缩的雷诺平均(RANS)方程数值模拟了包括兴波的三维船体周围的粘性流场,湍流模式使用了子网格尺度模式(SGS)和Baldwin-Lomax模式相结合的混合模式.对于自由表面的处理,采用了任意拉格朗日-欧拉方法,网格为不仅与物体表面贴体,而且与自由表面贴体的动网格,即随着自由表面的变化要不断地重新划分网格.虽然此方法需要很长的计算时间,但能较好地描述船体的兴波情况.本文计算了系列60船模在Fn=0.316,Re=1.9×106时带自由面的粘性流场,计算结果与试验结果吻合较好.，An finite - volume method solving of incompressible RANS Equations is developed for the numeri-cal simulation of three - dimensional viscous flow with free - surface about a ship. A hybrid turbulence mod-el by combining the sub - grid scale (SGS) model and the Baldwin- Lomax model is used in this pa-per. The arbitrary - Lagrange - Euler formulation is devised for the treatment of free surface boundary condi-tion. The boundary - fitted coordinate system is fitted not only to the ship hull surface but also to the freesurface, so the computational grid is regenerated to follow the free surface deformation. Although this methodneeds spend much longer computational time, but it can fairly display wave pattern made by a ship. The nu-merical simulation has been carried out for viscous flow with free surface past Series 60 ship hull at Fn =0. 316, Re = 1.9 × 106. The computed results are in reasonable agreement with experimental measure-ments.     

An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data

The air–sea CO₂ exchange is primarily determined by the boundary-layer processes in the near-surface layer of the ocean since it is a water-side limited gas. As a consequence, the interfacial component of the CO₂ transfer velocity can be linked to parameters of turbulence in the near-surface layer of the ocean. The development of remote sensing techniques provides a possibility to quantify the dissipation of the turbulent kinetic energy in the near-surface layer of the ocean and the air–sea CO₂ transfer velocity on a global scale. In this work, the dissipation rate of the turbulent kinetic energy in the near-surface layer of the ocean and its patchiness has been linked to the air–sea CO₂ transfer velocity with a boundary-layer type model. Field observations of upper ocean turbulence, laboratory studies, and the direct CO₂ flux measurements are used to validate the model. The model is then forced with the TOPEX POSEIDON wind speed and significant wave height to demonstrate its applicability for estimating the distribution of the near-surface turbulence dissipation rate and gas transfer velocity for an extended (decadal) time period. A future version of this remote sensing algorithm will incorporate directional wind/wave data being available from QUIKSCAT, a now-cast wave model, and satellite heat fluxes. The inclusion of microwave imagery from the Special Sensor Microwave Imager (SSM/I) and the Synthetic Aperture Radar (SAR) will provide additional information on the fractional whitecap coverage and sea surface turbulence patchiness.     

弯曲河道方形桥墩紊流宽度数值模拟

弯曲河道流态复杂,修建桥梁后桥墩对河道流场改变较大,对船舶航行安全非常不利。采用格子Boltzmann方法模拟桥墩三维紊流流场,并结合大涡模拟和运动边界处理方法,分析不同来流、弯曲半径等条件下的方形桥墩三维紊流宽度,由此建立了相对紊流宽度与弗劳德数之间耦合关系。利用上述成果,以四川新西林大桥为例进行验算,并与长江航道规划设计研究院的研究成果比较。结果发现,该方法得出的桥墩紊流宽度略大,但数值趋于合理。     

Using a random displacement model to simulate turbulent particle motion in a baroclinic frontal zone: A new implementation scheme and model performance tests

Numerical artifacts can limit accurate simulation of turbulent particle motion when Lagrangian particle-tracking models are implemented in hydrodynamic models with stratified conditions like fronts. Yet, modeling of individual particle motion in frontal regions is critical for understanding sediment dynamics as well as the transport and retention of planktonic organisms. The objective of this research was to develop a numerical technique to accurately simulate turbulent particle motions in a particle-tracking model embedded within a hydrodynamic model of a frontal zone. A new interpolation scheme, the ‘water column profile’ scheme, was developed and used to implement a random displacement model for turbulent particle motions. A new interpolation scheme was necessary because linear interpolation schemes caused artificial aggregation of particles where abrupt changes in vertical diffusivity occurred. The new ‘water column profile’ scheme was used to fit a continuous function (a tension spline) to a smoothed profile of vertical diffusivities at the x–y particle location. The new implementation scheme was checked for artifacts and compared with a standard random walk model using (1) Well Mixed Condition tests, and (2) dye-release experiments. The Well Mixed Condition tests confirmed that the use of the ‘water column profile’ interpolation scheme for implementing the random displacement model significantly reduced numerical artifacts. In dye-release experiments, high concentrations of Eulerian tracer and Lagrangian particles were released at the same location up-estuary of the salt front and tracked for 4 days. After small differences in initial dispersal rates, tracer and particle distributions remained highly correlated (r = 0.84 to 0.99) when a random displacement model was implemented in the particle-tracking model. In contrast, correlation coefficients were substantially lower (r = 0.07 to 0.58) when a random walk model was implemented. In general, model performance tests indicated that the ‘water column interpolation’ scheme was an effective technique for implementing a random displacement model within a hydrodynamic model, and both could be used to accurately simulate diffusion in a highly baroclinic frontal region. The new implementation scheme has the potential to be a useful tool for investigating the influence of hydrodynamic variability on the transport of sediment particles and planktonic organisms in frontal zones.     

Turbulent mixing and internal tides in Gaoping (Kaoping) Submarine Canyon, Taiwan

Turbulent overturning on scales greater than 10 m is observed near the bottom and in mid-depth layers within the Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) in southern Taiwan. Bursts of strong turbulence coexist with bursts of strong sediment concentrations in mid-depth layers. The turbulence kinetic energy dissipation rate in some turbulence bursts exceeds 10^− 4 W kg^− 1, and the eddy diffusivity exceeds 10^− 1 m² s^− 1. Within the canyon, the depth averaged turbulence kinetic energy dissipation rate is ~ 7 × 10^− 6 W kg^− 1, and the depth averaged eddy diffusivity is ~ 10^− 2 m² s^− 1. These are more than two orders of magnitude greater than typical values in the open ocean, and are much larger than those found in the Monterey Canyon where the strong turbulent mixing has also been. The interaction of tidal currents with the complex topography in Gaoping Submarine Canyon is presumably responsible for the observed turbulent overturning via shear instability and the breaking of internal tides and internal waves at critical frequencies. Strong 1st-mode internal tides exist in KPSC. The depth averaged internal tidal energy near the canyon mouth is ~ 0.17 m² s^− 2. The depth integrated internal tidal energy flux at the mouth of the canyon is ~ 14 kW m^− 1, propagating along the axis of the canyon toward the canyon head. The internal tidal energy flux in the canyon is 3–7 times greater than that found in Monterey Canyon, presumably due to the more than 10 times larger barotropic tide in the canyon. Simple energy budget calculations conclude that internal tides alone may provide energy sufficient to explain the turbulent mixing estimated within the canyon. Further experiments are needed in order to quantify the seasonal and geographical distributions of internal tides in Gaoping Submarine Canyon and their effects on the sediment flux in the canyon.     

孔腔脉动压力及其波数—频率谱的大涡模拟研究

孔腔流动中包含着流动分离和失稳以及涡旋相互干扰等复杂的流动现象。孔腔涡旋流动引起的流体振荡能够引起脉动压力的显著增加从而产生强烈的噪声,在工程实际中备受关注。湍流脉动压力是流激噪声的重要来源,也是湍流研究中的基础性问题,对其进行数值计算研究是流声耦合领域的重要内容,而湍流脉动压力波数—频率谱的构建更是该领域的技术难点。文章采用大涡模拟方法(LES)对孔腔脉动压力进行了数值模拟,考察了四套网格和四种亚格子应力模型对计算结果的影响,并与试验结果进行比较,验证数值计算方法的可靠性。首先采用大涡模拟方法计算了孔腔的脉动压力,并与中国船舶科学研究中心的空泡水筒试验结果进行对比分析。接着详细地分析孔腔脉动压力,研究亚格子应力模型和网格数量对计算结果的影响。最后,对数值计算得到的脉动压力多元阵列结果进行时间/空间Fourier变换,构建了三维脉动压力波数-频率谱。该文工作对今后流激结构振动噪声的预报和流动控制研究奠定了基础。