Three-layer model for vertical velocity distribution in open channel flow with submerged rigid vegetation

Based on the detailed laboratory experiments and theoretical analysis, a new three-layer model is proposed to predict the vertical velocity distribution in an open channel flow with submerged vegetation. The time averaged velocity and turbulence behaviour of a steady uniform flow with fully submerged artificial rigid vegetation was measured using a 3D Micro ADV, and the vertical distribution of velocity and Reynolds shear stress at different vegetation height, vegetation density and measuring positions were obtained. The results show that the velocity profile consists of three hydrodynamic regimes (i.e. the upper non-vegetated layer, the outer and bottom layer within vegetation); accordingly different methods had been adopted to describe the vertical velocity distribution. For the upper non-vegetated layer, a modified mixing length theory combined with the concept of ‘the new vegetation boundary layer’ was adopted, and an analytical model was presented to predict the vertical velocity distribution in this region. For the bottom layer within vegetation, the depth average velocity was obtained by numerically solving the momentum equations. For the upper layer within vegetation, the analytical solution was presented by expressing the shear stress as a formula fitted to the experimental data. Finally, the analytical predictions of the vertical velocity over the whole flow depth were compared with the results obtained by other researchers, and the good agreement proved that the three-layer model can be used to predict the velocity distribution of the open channel flow with submerged rigid vegetation.     

Fluid flow dynamics under location uncertainty

We present a derivation of a stochastic model of Navier Stokes equations that relies on a decomposition of the velocity fields into a differentiable drift component and a time uncorrelated uncertainty random term. This type of decomposition is reminiscent in spirit to the classical Reynolds decomposition. However, the random velocity fluctuations considered here are not differentiable with respect to time, and they must be handled through stochastic calculus. The dynamics associated with the differentiable drift component is derived from a stochastic version of the Reynolds transport theorem. It includes in its general form an uncertainty dependent subgrid bulk formula that cannot be immediately related to the usual Boussinesq eddy viscosity assumption constructed from thermal molecular agitation analogy. This formulation, emerging from uncertainties on the fluid parcels location, explains with another viewpoint some subgrid eddy diffusion models currently used in computational fluid dynamics or in geophysical sciences and paves the way for new large-scales flow modeling. We finally describe an applications of our formalism to the derivation of stochastic versions of the Shallow water equations or to the definition of reduced order dynamical systems.     

横向沙丘表面气流湍流特征野外观测

沙丘形态动力学特征是风沙地貌的重要研究内容之一,受观测技术等的限制,现有的关于沙丘形态动力学特征的研究成果还不能满足沙丘形成与演化过程研究的需要。利用三维超声风速仪对腾格里沙漠的横向沙丘表面气流特征进行野外观测,结果表明,沙丘脊线和迎风坡中部的气流均为未偏向流,而沙丘背风坡不同高度的气流包括反向流,偏向流和未偏向流。湍流强度、雷诺应力和气流丰富度随气流方向变化而变化,其在沙丘脊线和中部的变化趋势比较简单且有大致相似的规律,而其在沙丘背风坡不同高度比较复杂。     

Numerical Simulation on Hydraulic Performances of Quarter Circular Breakwater

Quarter circular breakwater （QCB） is a new-type breakwater developed from senti-circular breakwater （SCB）. The superstructure of QCB is composed of a quarter circular front wall, a horizontal base slab and a vertical rear wall. The width of QCB＇ s base slab is about half that of SCB, which makes QCB suitable to be used on relatively finn soil foundation. The numerical wave flume based on the Reynolds averaged Navier-Stokes equations for impressible viscosity fluid is adopted in this paper to simulate the hydraulic performances of QCB. Since the geometry of both breakwaters is similar and SCB has been studied in depth, the hydraulic performances of QCB are given in comparison with those of SCB.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

On the structure of supercritical western boundary currents

The structure of supercritical western boundary currents is investigated using a quasi-geostrophic numerical model. The basic flow is of meridional Munk balance, and the input boundary is perturbed by the most unstable wave solution obtained from linear spatial instability calculations. Self-preserving (or equilibrium) solutions are obtained for the model runs at Re=30, 60, 90, and 120, and their energy and vorticity budgets are analyzed. In an analogy with the laboratory turbulence of wall boundary layers, the western boundary layer is divided into inner and outer layers. In the inner layer, the mean energy is dissipated via direct viscous dissipation, while in the outer layer it is converted to the eddy energy via turbulence production. The main scenario is that the mean energy is produced in the inner layer via ageostrophic pressure work divergence, and it is partly removed due to viscous action within a narrow region near the wall, defined here as viscous sub-layer. The remaining portion is converted to the eddy energy via turbulence production in the outer layer, which is in turn transported to the inner layer, then again to the viscous sub-layer where it is ultimately dissipated. In the near-wall side, the vorticity balance of the mean flow is maintained by viscous effect and Reynolds flux divergence, while in the offshore side it is maintained by beta effect and Reynolds flux divergence. The length scale of the supercritical boundary current is roughly , where L_M is the Munk length, as observed from a dimensional analysis.     

Interaction of Streamwise and Wall-Normal Velocities in Combined Wave-Current Motion

The aim of this paper is to present an analytical expression for the streamwise velocity distribution in a non-uniform flow in the presence of waves; the correlation between the horizontal and vertical velocity components has been compreheusively examined. Different from previous researches which attributed the deviation of velocity from the classical log-law to the wave Reynolds stress, i.e. - ρ uv^- only, this study demonstrates that the momentum flux caused by mean velocities, i.e., u^- and v^-, is also responsible for the velocity deviation, and it is found that the streamwise velocity for a flow in the presence of non-zero wall-normal velocity does not follow the classical log-law, but the modified log-law proposed in this study based on simplified mixing-length theorem. The validity of the modified log-law has been verified by use of available experimental data from published sources for combined wave-current flows, and good agreement between the predicted and observed velocity profiles has been achieved.     

波浪要素对水平圆柱杆件水动力系数影响研究

通过试验研究了波浪要素对水平固定圆柱杆件水动力系数的影响。试验利用垂直杠杆原理,提取了杆件在波浪下的水平受力历时曲线,基于改造的莫里森方程,计算了不同波浪要素下的水动力系数C_D和C_M取值。定义分析了新的波高参数KH,周期参数KT和雷诺数Re对水动力系数的影响。研究表明,整体拟合法与四点拟合法获得的水动力系数在数值上差异不大。KH、KT及Re均对水动力系数的取值有较显著影响:1)随着波高KH的增大,C_D、C_M整体呈幂指数衰减,且C_M的衰减更为迅速;2)随着波周期KT的增大,C_D、C_M整体亦呈衰减趋势;3)相同波高条件下,C_M随Re数增大而增大,而C_D值相对稳定。最后给出了波浪条件下水动力系数C_D和C_M的经验计算公式。     

Large-eddy simulation of the flow past a circular cylinder at sub- to super-critical Reynolds numbers

Large-eddy simulation of turbulent flow past a circular cylinder at sub- to super-critical Reynolds numbers is performed using a high-fidelity orthogonal curvilinear grid solver. Verification studies investigate the effects of grid resolution, aspect ratio and convection scheme. Monotonic convergence is achieved in grid convergence studies. Validation studies use all available experimental benchmark data. Although the grids are relatively large and fine enough for sufficiently resolved turbulence near the cylinder, the grid uncertainties are large indicating the need for even finer grids. Large aspect ratio is required for sub-critical Reynolds number cases, whereas small aspect ratio is sufficient for critical and super-critical Reynolds number cases. All the experimental trends were predicted with reasonable accuracy, in consideration the large facility bias, age of most of the data, and differences between experimental and computational setup in particular free stream turbulence and roughness. The largest errors were for under prediction of turbulence separation.