Direct numerical simulation of a turbulent stably stratified air flow above a wavy water surface

The influence of the roughness of the underlaying water surface on turbulence is studied in a stably stratified boundary layer (SSBL). Direct numerical simulation (DNS) is conducted at various Reynolds (Re) and Richardson (Ri) numbers and the wave steepness ka. It is shown that, at constant Re, the stationary turbulent regime is set in at Ri below the threshold value Ri_c depending on Re. At Ri > Ri_c, in the absence of turbulent fluctuations near the wave water surface, three-dimensional quasiperiodical structures are identified and their threshold of origin depends on the steepness of the surface wave on the water surface. This regime is called a wave pumping regime. The formation of three-dimensional structures is explained by the development of parametric instability of the disturbances induced by the surface water in the air flow. The DNS results are quite consistent with prediction of the theoretical model of the SSBL flow, in which solutions for the disturbances of the fields of velocity and temperature in the wave pumping regime are found to be a solution of a two-dimensional linearized system with the heterogeneous boundary condition, which is caused by the presence of the surface wave. In addition to the turbulent fluctuations, the three-dimensional structures in the wave pumping regime provide for the transfer of impulse and heat, i.e., the increase in the roughness of the water–air boundary caused by the presence of waves intensifies the exchange in the SSBL.     

明渠恒定紊流的总流控制方程与能量损失特性

为进一步深化对明渠总流控制方程的认识,探究总流能量损失的构成与分布特点,通过理论分析,在黏性流体力学理论的框架下构建了描述明渠恒定紊流总体特性的总流积分模型与总流微分模型的控制方程,其模型参数能直接通过流动统计特征量在过流断面上的分布来获取,由此实现明渠流的流场特性描述与总流描述的统一,并得到了总流能量损失的显示表达式。同时指出,总流能量损失由黏性耗散与维持紊动两部分构成,在矩形明渠恒定均匀流紊流的分类能量损失构成中,壁面上以黏性耗散部分为主,维持紊动部分的能量损失密度随着离开壁面距离的增加而快速增加,且随着雷诺数的增大,其增加速度也越快。     

Statistical properties of the atmospheric turbulent boundary layer over steep surface waves

A new method of digital optical anemometry (Particle Image Velocimetry, PIV) of turbulent flows is suggested and implemented in the laboratory; it is based on the use of continuous laser radiation and high-speed video photography, providing continuous statistical ensembles of flow velocity fields. Application of the method to the study of wind field over waves has allowed us to perform, for the first time, direct measurements of velocity fields, averaged over turbulent pulsations induced by waves in the air flow. The experiments demonstrated that the velocity fields, averaged over the turbulent pulsations, are nonseparated even in the case of steep and breaking waves, when separation of the flow from the wave crests in the instantaneous fields is observed. Based on comparison with the experimental data, it is shown that the average wind fields over waves are described well quantitatively in the framework of semiempirical closure models of turbulence.     

初始裂隙对岩溶水紊流形成的影响

岩溶地区地下发育着大量的溶洞和地下河管道,地下水流状态既有层流也有紊流,而紊流是溶洞管道形成的重要条件。紊流的形成受到岩石初始裂隙的影响,初始裂隙的张开度、分布、走向、迹长、密度等因素都影响着裂隙发育过程中水流状态的变化。通过对不同统计特征的初始裂隙网络进行水流和溶蚀的数值模拟发现,以张开度标准差反映的裂隙网络非均匀性越强,模拟紊流出现的时间就越早;主要裂隙的存在使裂隙网络的非均性增强,主要裂隙与水力梯度总方向的角度越小,紊流出现的时间就越早;当裂隙平均迹长过小时会导致裂隙连通性较差,影响裂隙水流和溶蚀作用;裂隙密度,尤其是主要裂隙密度,对岩溶发育的影响较大。相对于次要裂隙,如果主要裂隙密度偏小,紊流形成时间会大大增加,甚至很难形成紊流。当初始裂隙张开度小于0.001 cm,增大水力梯度仍没有紊流发生,岩溶几乎不发育。      

Excitation of turbulence in accretion disks of binary stars by non-linear perturbations

Accretion disks in binary systems can experience hydrodynamical influences at both their inner and outer edges. The former is typical for protoplanetary disks around young T Tauri stars, while the latter is typical for circumstellar disks in close binaries. This influence excites perturbations with various scales and amplitudes in the disk. The nonlinear evolution of perturbations with a finite, but small amplitude against the background of a sub-Keplerian flow is investigated. Nonlinear effects at the fronts of perturbation waves lead to the formation of discontinuities in the density and radial velocity; i.e., to formation of shocks. The tangential flow in the neighborhood of the shock becomes equivalent to a flow in a boundary layer. Due to an instability of the tangential flow, the disk becomes turbulent. The characteristics of the turbulence depend on the parameters of the perturbations, but the Shakura–Syunyaev α parameter does not exceed ~0.1.     

Toward modeling of CO2 multi-phase flow patterns using a stochastic multi-scale approach

Although pipe transport for storage of CO₂ captured from combustion exhaust is most efficient under supercritical conditions, subcritical multi-phase transport might sometimes occur intentionally or unintentionally. To adequately assess the consequences of subcriticality, the fidelity of subcritical multi-phase turbulent transport modeling, including confinement and buoyancy effects, must be improved. For this reason, an extension of one-dimensional turbulence, a stochastic turbulence modeling strategy, for application to this regime is underway. As a step toward this extension, a turbulent liquid jet issuing from a planar channel is simulated and results are compared with liquid jet measurements and channel flow numerical simulations. A previously noted turbulence decay scaling is reproduced, suggesting that the scaling is more robust than might be supposed, given the complicating factors. This work provides a basis for extension to the primary breakup regime of liquid jets and hence to subcritical multi-phase turbulence more generally.     

近岸波生流运动三维数值模拟及验证

开发建立了近岸波生流运动三维数值计算模式。模式中,引入了三维时均剩余动量、破波表面水滚、波浪水平与垂向紊动作为主要驱动力,同时考虑了波流共同作用的底部剪切力。推导了可综合反映底坡、能量传递率和密度影响的水滚能量传输方程;将Larson-Kraus的二维波浪水平紊动系数表达式拓展至三维。采用大量实测数据和文献资料测试验证了所建模式,表明所建模式可有效模拟波浪增减水、底部离岸流、沿岸流、裂流、堤后环流等不同维度的波生流现象。此外,研究也表明破波水滚效应可解释波生流峰值向岸推移的物理现象,从而在模拟中不能忽略;破波带内沿岸流速垂向较为均匀的现象与波浪附加垂向紊动有关。     

The dynamics of turbulent, transitional and laminar clay-laden flow over a fixed current ripple

Most aqueous sedimentary environments contain varying concentrations of fine‐grained, often clay‐rich, sediment that is transported in suspension and may modify the properties of the flow and underlying mobile bed. This paper presents results from a series of laboratory experiments examining the mean and turbulent properties of clay‐laden (kaolinite) flows, of various volumetric sediment concentrations between 0·046% and 12·7%, moving over a fixed, idealized current ripple. As the kaolinite concentration was raised, with flow velocity and depth constant, four flow types were observed to occur: (i) turbulent flow, in which flow separation is dominant in the leeside of the ripple; (ii) turbulence‐enhanced transitional flow, in which turbulence in the leeside separation zone region is enhanced; (iii) turbulence‐attenuated transitional flow, in which turbulence along the separation zone shear layer and in the free flow above it becomes damped, eventually leading to a reduction in the size of the separation zone wake region; and (iv) laminar plug flow, in which turbulence is damped and flow is almost stagnant in the lee of the ripple. Such modulation of turbulence by increasing clay concentrations suggests that many paradigms of flow and bedform dynamics, which have been based on extensive past work in clear water flows, require revision. The present results highlight a need to fully characterize the boundary conditions for turbulence modulation as a function of clay type and applied flow conditions, and the effects of such flows on fully mobile cohesionless beds.     

On the relationship between flow and suspended sediment transport over the crest of a sand dune, Río Paraná, Argentina

The links between large‐scale turbulence and the suspension of sediment over alluvial bedforms have generated considerable interest in the last few decades, with past studies illustrating the origin of such turbulence and its influence on flow resistance, sediment transport and bedform morphology. In this study of turbulence and sediment suspension over large sand dunes in the Río Paraná, Argentina, time series of three‐dimensional velocity, and at‐a‐point suspended sediment concentration and particle‐size, were measured with an acoustic Doppler current profiler and laser in situ scattering transmissometer, respectively. These time series were decomposed using wavelet analysis to investigate the scales of covariation of flow velocity and suspended sediment. The analysis reveals an inverse relationship between streamwise and vertical velocities over the dune crest, where streamwise flow deceleration is linked to the vertical flux of fluid towards the water surface in the form of large turbulent fluid ejections. Regions of high suspended sediment concentration are found to correlate well with such events. The frequencies of these turbulent events have been assessed from wavelet analysis and found to concentrate in two zones that closely match predictions from empirical equations. Such a finding suggests that a combination and interaction of vortex shedding and wake flapping/changing length of the lee‐side separation zone are the principal contributors to the turbulent flow field associated with such large alluvial sand dunes. Wavelet analysis provides insight upon the temporal and spatial evolution of these coherent flow structures, including information on the topology of dune‐related turbulent flow structures. At the flow stage investigated, the turbulent flow events, and their associated high suspended sediment concentrations, are seen to grow with height above the bed until a threshold height (ca 0·45 flow depth) is reached, above which they begin to decay and dissipate.     

明渠恒定紊流的总流控制方程与能量损失特性

为进一步深化对明渠总流控制方程的认识,探究总流能量损失的构成与分布特点,通过理论分析,在黏性流体力学理论的框架下构建了描述明渠恒定紊流总体特性的总流积分模型与总流微分模型的控制方程,其模型参数能直接通过流动统计特征量在过流断面上的分布来获取,由此实现明渠流的流场特性描述与总流描述的统一,并得到了总流能量损失的显示表达式。同时指出,总流能量损失由黏性耗散与维持紊动两部分构成,在矩形明渠恒定均匀流紊流的分类能量损失构成中,壁面上以黏性耗散部分为主,维持紊动部分的能量损失密度随着离开壁面距离的增加而快速增加,且随着雷诺数的增大,其增加速度也越快。     

浅水非线性波作用下沙纹床面底层流动特性试验研究

通过水槽试验研究浅水非线性波作用下沙纹床面底层流动特性,利用CCD图像技术观测分析非对称沙纹的形成和演化规律。利用声学多普勒测速仪(ADV)测量非对称沙纹底床上的流场,得到了不同波高、周期、水深条件下的沙纹峰顶和谷底断面的瞬时速度。试验结果分析表明,浅水非线性波作用下床面上形成非对称沙纹,其近底流速具有较强紊动特性,随着距床面距离的增大紊动强度逐渐减弱。在水流方向改变时,沙纹背部具有明显漩涡运动。沙纹背后形成的漩涡能起到维持沙纹的作用。浅水非线性波作用下,沙纹的形成原因主要是床面泥沙颗粒在非对称流动和床面近壁粘性底层中漩涡结构动力作用下,作受迫摆动、推移所致。     

Formation of accretion disks in close-binary systems with magnetic fields

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from 10⁵ to 10⁸ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.     

The vertical turbulence structure of experimental turbidity currents encountering basal obstructions: implications for vertical suspended sediment distribution in non‐equilibrium currents

Large roughness features, caused by erosion of the sea floor, are commonly observed on the modern sea floor and beneath turbidite sandstone beds in outcrop. This paper aims to investigate the effect of such roughness elements on the turbulent velocity field and its consequences for the sediment carrying capacity of the flows. Experimental turbidity currents were run through a rectangular channel, with a single roughness element fixed to the bottom in some runs. The effect of this roughness element on the turbulent velocity field was determined by measuring vertical profiles of the vertical velocity component in the region downstream of the basal obstruction with the Ultrasonic Doppler Velocity Profiling technique. The experiments were set up to answer two research questions. (i) How does a single roughness element alter the distribution of vertical turbulence intensity? (ii) How does the altered profile evolve in the downstream direction? The results for runs over a plane substrate are similar to data presented previously and show a lower turbulence maximum near the channel floor, a turbulence minimum associated with the velocity maximum, and a turbulence maximum associated with the upper flow interface. In the runs in which the flows were perturbed by the single roughness element, the intensity of the lower turbulence maximum was increased between 41% to 81%. This excess turbulence dissipated upwards in the flow while it travelled further downstream, but was still observable at the most distal measurement location (at a distance ca 39 times the roughness height downstream of the element). All results point towards a similarity between the near bed turbulence structure of turbidity currents and free surface shear flows that has been proposed by previous authors, and this proposition is supported further by the apparent success of a shear velocity estimation method that is based on this similarity. Theory of turbulent dispersal of suspended sediment is used to discuss how the observed turbulent effects of a single large roughness element may impact on the suspended sediment distribution in real world turbidity currents. It is concluded that this impact may consist of a non‐equilibrium net‐upwards transport of suspended sediment, counteracting density stratification. Thus, erosive substrate topography created by frontal parts of natural turbidity flows may super‐elevate sediment concentrations in upper regions above equilibrium values in following flow stages, delay depletion of the flow via sedimentation and increase their run‐out distance.     

斜坡上异重流的三维数值模拟

针对异重流的流动特征,建立了适用于具有各向异性浮力紊动特征的三维异重流运动的数学模型,并模拟了异重流在15°斜坡底面上的潜行过程。计算结果准确地模拟了异重流的运动特征和形态,其前锋的潜行速度与实验结果相当吻合。该模型采用非结构同位网格上的SIMPLEC算法能适应复杂边界和地形,可应用于自然界实际环境中异重流的演进计算。     

温度和盐度分层流的数值模拟

以N-S方程和k-ε湍流模型为基础,针对水温和盐度分层流的流动特性,考虑了浮力对紊动的影响,建立了温度和盐度分层流的k-ε模型,对同时存在温度和盐度梯度的一类密度分层流进行了数值模拟计算,其中紊动平均量控制方程中的紊动项是由k-ε紊流模型计算,给出了速度、温度和盐度的分布。计算结果正确地反映了分层及紊动特征,与实测结果吻合较好。     

Amplification of fluctuations and current dynamo due to helical and chiral effects in geophysical and plasma-like media

Dynamo effect is considered in more general case than well known one when not only mean motion but also conducting components possess nonzero mean helicity. Dispersion equation for helical motions is studied, including inhomogeneous case. Criteria for development of instability are found. Interaction between large-scale internal wave and small-scale helical turbulence in plane Couette flow of fluid with statically stable uniform density gradient in gravitational field is studied basing on a set of equations for scalar and quasi-scalar wave fields and for helicity of turbulence. This system of equation describes acceleration of wave growth due to its backward action on turbulence, i.e. wave-turbulent instability. Chiral media are further considered and effects topologically close to helical ones are analyzed. It is shown that in such media ranges of scales exist, in which fluctuations of electric field can be anomalously amplified.     

Quantifying velocity and turbulence structure in depositing sustained turbidity currents across breaks in slope

Two-dimensional experiments investigating sediment transport and turbulence structure in sustained turbidity currents that cross breaks in slope are presented as analogue illustrations for natural flows. The results suggest that in natural flows, turbulence generation at slope breaks may account for increased sand transport into basins and that the formation of a hydraulic jump may not be necessary to explain features such as the occurrence of submarine plunge pools and the deposition of coarser-grained beds in the bottomsets of Gilbert-type fan deltas. Experimental flows were generated on 0°, 3°, 6° and 9° slopes of equal length which terminated abruptly on a horizontal bed. Two-component velocities were measured on the slope, at the slope break and downstream of the slope break. Flows were depositional and non-uniform, visibly slowing and thickening with distance downstream. One-dimensional continuous wavelet transforms of velocity time series were used to produce time-period variance maps. Peaks in variance were tested against a background red-noise spectrum at the 95% level; a significant period banding occurs in the cross-wavelet transform at the slope break, attributed to increased formation of coherent flow structures (Kelvin–Helmholtz billows). Variance becomes distributed at progressively longer periods and the shape of the bed-normal-velocity spectral energy distribution changes with distance downstream. This is attributed to a shift towards larger turbulent structures caused by wake stretching. Mean velocity, Reynolds shear stress and turbulent kinetic energy profiles illustrate the mean distribution of turbulence through the currents. A turbulent kinetic energy transfer balance shows that flow non-uniformity arises through the transfer of mean streamwise slowing to mean bed-normal motion through the action of Reynolds normal stresses. Net turbulence production through the action of normal stresses is achieved on steeper slopes as turbulence dissipation due to mean bed-normal motion is limited. At the slope break, an imbalance between the production and dissipation of turbulence occurs because of the contrasting nature of the wall and free-shear boundaries at the bottom and top of the flows, respectively. A rapid reduction in mean streamwise velocity predominately affects the base of the flows and steeper proximal slope flows have to slow more at the break in slope. The increased turbulent kinetic energy, limited bed-normal motion and strong mixing imposed by steep proximal slopes means rapid slowing enhances turbulence production at the break in slope by focusing energy into coherent flow structures at a characteristic period. Thus, mean streamwise slowing is transferred into turbulence production at the slope break that causes increased transport of sediment and a decrease in deposit mass downstream of the slope break. The internal effects of flow non-uniformity therefore can be separated from the external influence of the slope break.     

采样频率及样本容量对明渠紊流统计值的影响

为了研究采样频率和样本容量对明渠紊流统计特征值的影响,利用粒子图像测速技术,对明渠恒定均匀水流的二维流速分布进行测量,获得了包含340535个样本的流场时间序列.对该大容量时间序列进行抽样,得到不同采样频率和容量条件下的一系列子序列.对比分析各子序列的时均流速、紊动强度及雷诺应力的变化特点,发现存在临界的采样频率和样本容量,紊动统计特征参数在临界值以下波动较大,而在临界值以上基本趋于稳定;要满足同一测量精度,时均流速对采样频率和样本容量的要求最低,紊动强度次之,雷诺应力最高.     

弯型截止阀流场的PIV显示和数值模拟

为了解弯型进口截止阀的流场特性,用粒子成像流速仪(PIV)对弯型截止阀对称面模型进行了流场显示;此外,采用了RNGk ε紊流模型和贴体坐标对弯型进口和斜进口截止阀对称面流场进行了模拟。通过用粒子成像流速仪流场显示技术和数值计算,揭示了弯型进口和斜进口截止阀流场特性,数值计算的流态与实验结果也较为吻合。计算结果显示斜进口截止阀流态好于弯型进口截止阀的流态。     

Prediction of magnetite segregation in dense medium cyclone using computational fluid dynamics technique

Multi-phase simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted in Fluent, using the Algebraic Slip Mixture model to model the dispersed phases and the air-core, and both the Large Eddy Simulation turbulence model (LES) and Reynolds Stress Models (RSM) for turbulence closure. The predicted air-core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible to use the LES turbulence model with ASM multi-phase model to predict the air/slurry interface accurately. Multi-phase simulations (air/water/medium) show appropriate medium segregation effects but over-predict the level of segregation compared to that measured by gamma ray tomography near the wall. This is believed to be because of unaccounted back-mixing of the dispersed phase due to turbulence in the basic Algebraic Slip Mixture model. The predictions of accurate axial segregation of magnetite medium are investigated using the LES turbulence model in conjunction with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. At higher feed densities the agreement between the Dunglison [Dunglison, M.E., 1999, A general model of the dense medium cylone, PhD thesis, JKMRC, University of Queensland] correlations and experimental measurements and the CFD is reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. The effect of size distribution of the magnetite has been fully studied. As expected, the ultra-fine magnetite sizes (i.e. 2 and 7 microns) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall.