A VOF-based numerical model for breaking waves in surf zone

This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes (RANS) equations for a mean flow, and (he k-s equations for turbulence kinetic energy k and turbulence dissipation rate e. To track a free surface, the volume of fluid (VOF) function, satisfying the advection equation was introduced. In the numerical treatment, third-order upwind difference scheme was applied to the convection terms of the RANS equations in order to reduce the effect of numerical viscosity. The shoaling and breaking processes of a periodic wave train on gently sloping beaches were modeled. The computed wave heights of a sloping beach and the distribution of breaking wave pressure on a vertical wall were compared with laboratory data.     

Parameterization of ocean wave-induced mixing processes for finite water depth

Three dimensional wave-induced mixing plays an important role in shallow water area. A quite direct approach through the Reynolds average upon characteristic length scale is proposed to parameterize the horizontal and vertical shallow water mixing. Comparison of finite depth case with infinite depth results indicates that the difference of the wave-induced mixing strength is evident. In the shallow water condition, the infinite water depth approximation overestimates the mixing strength in the lower layers. The nonzero horizontal wave-induced mixing presents anisotropic property near the shore. The Prandtl’s mixing length theory underestimated the wave-induced mixing in the previous studies.     

Sediment transport under the presence and absence of emergent vegetation in a natural alluvial channel from Brazil

Studies regarding the influence of emergent vegetation on sediment transport are scarce and have mainly focused on flume conditions. To fill this gap and also meet the international need, we aimed to evaluate the influence of emergent vegetation (Echinodorus macrophyllus) on sediment transport of Capibaribe River, Brazil. Bedload and suspended sediment measurements were carried out using the US BLH 84 and US DH 48 samplers, respectively. Measurements of stem diameter, stem spacing and plant density were performed in conjunction with flow and sediment field measurements. Based on our results, 0.45 m s ? 1 was the threshold of mean flow velocity supported by E. macrophyllus under field conditions. This value can be helpful for other rivers with gravel-bed river to armoured layer ratio (AR ? D50-surface)/D50-subsurface ? 12.50) – natural conditions observed in Capibaribe River – or where the vegetation can provide positive effects, such as increase the bed stability, assist water restoration/rehabilitation and decrease water turbidity. Our results can hopefully be used in engineering practice and ecosystem management. In general, both the drag coefficient and drag force varied inversely and directly with the mean flow velocity and vegetation density, respectively. The vegetation resistance force was inversely proportional to the bedload transport owing to the resistance caused by emergent vegetation. This finding was supported by the clear decoupling between nonvegetated and emergent vegetated conditions indicated by cluster analysis. The study results provided a reasonable understanding of the interaction between emergent vegetation, water flow and sediment transport in the Capibaribe River.&2016 International Research and Training Centre on Erosion and Sedimentation/the World Association     

Wind-Driven Ocean Circulation in Shallow Water Lattice Boltzmann Model

A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximation for the collision operator, the model becomes fully explicit. In this case, any iterative technique is not needed. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretized accuracy of the LB equation. The numerical results show correct physics of the ocean circulation driven by the double-gyre wind stress with different Reynolds numbers and different spatial resolutions. An intrinsic low-frequency variability of the shallow water model is also found. The wind-driven ocean circulation exhibits subannual and interannual oscillations, which are comparable to those of models in which the conventional numerical methods are used.     

圆形断面管道非满流水动力特性

为探究边壁双曲率对其水动力特性的影响，应用RNG k-ε紊流模型和部分面积体积障碍模拟技术（FAVOR），针对圆管非满流建立三维数值模型，通过物理模型试验结果验证该数值模型的可行性。针对底坡和充满度不同组合条件进行数值试验，分析断面流速、壁面切应力和雷诺切应力分布规律等。研究结果表明：圆管非满流断面上不同垂线的流速分布具有较好的相似性，据此提出了垂线流速分布的抛物线公式；针对不同充满度条件提出了壁面切应力沿湿周分布的统一表达式；雷诺切应力沿垂线服从线性分布规律，充满度越大且距中垂线越远，沿垂线的变化梯度越小；当充满度大于0.5时，由于受二次流影响断面中垂线上雷诺切应力出现负值。     

挡潮闸矩形孔口式鱼道紊流结构及放鱼试验研究

对挡潮闸枢纽中矩形中孔、底孔鱼道中紊流结构进行了较为系统的试验研究,并做了放鱼试验。选择了一种鱼类偏爱流速所对应的流量作为典型流量,考虑了不同的孔口位置(中孔和底孔),用声学多普勒测速仪(ADV)量测了测点的三维瞬时流速及流向,分析了矩形孔口鱼道的三维时均流速分布特征、断面最大流速沿程变化规律、流速矢量场、紊动强度分布及雷诺应力分布。此外,还通过放鱼试验,利用在鱼体植入T形标签和高速摄影机观察了过鱼对象对中孔、底孔的反应情况,分析了过鱼对象与矩形孔口鱼道紊流结构的关系。试验结果表明:水流经中孔形成三维紊动自由射流,经底孔形成三维壁面射流,中孔纵向流速呈高斯分布,而底孔纵向流速则近似为高斯分布,流速由孔口向两侧逐渐减小;中孔和底孔横向流速在位于孔口范围内的纵剖面上沿程减小,孔口之外则变化较小;中孔和底孔垂向流速分布特征表现为在铅垂方向上均存在旋涡;在中孔水平面和纵剖面上,纵向最大流速均沿程衰减;中孔和底孔情形孔口处紊动强度和雷诺应力比非孔口处大得多,而非孔口处不同水深平面上紊动强度和雷诺应力变化趋于平缓;过鱼对象喜爱在紊动强度分布的峰值区和雷诺应力较大变幅区溯游。     

苦草对水流结构的影响研究

为研究沉水植被对水流结构的影响,以苦草为对象,利用实验室水槽实验研究了含淹没水生植被水流时均流速、雷诺应力及紊动能的垂向分布特征。受苦草冠层的影响,时均流速在冠顶以上呈对数分布,且随着流量的增加,冠层倾伏高度降低,对数剖面愈加明显;冠层内部,由于冠层阻流面积在垂向分布上的差异,冠层内时均流速出现逆梯度分布,且在床面附近出现局部流速最大值。雷诺应力在冠顶附近达到最大值,并向水面与床底方向逐渐减小;受逆流速梯度的影响,冠层内部雷诺应力出现负值以及局部最大值。雷诺应力产生的剪切紊动使得紊动能在冠顶处最大,并向水面与床底进行垂向传输;受紊动传输距离的限制,冠层底部以叶片后产生的尾流紊动为主,紊动能较小。     

Turbulent and Macro-turbulent Structures Developed in the Benthic Boundary Layer Downstream of Topographic Features

The characteristics and effects of large-scale flow structures developed in the benthic boundary layer downstream from large topographic features were analysed throughout a tidal cycle. The observed signature of the macro-turbulent features consisted of streamwise modules of low horizontal velocity and high suspended sediment concentration (SSC), alternating with modules of high horizontal velocity and low SSC. These modules extended 10 to 20 m streamwise and exceeded 1 m vertically, and are believed to be related to flow separation effects over large bedforms upstream of the deployment site. The macroscale flow modules intensified the ‘ burst-like ’ turbulent events and favoured sediment transport. ‘ Ejection-like ’ events were magnified during modules of decreasing horizontal velocity and increasing turbidity, whereas ‘ sweep-like ’ events were magnified during modules of increasing horizontal velocity and decreasing SSC. The enhanced turbidity of the macroscale modules may be the result of enhanced upward diffusion of sediment by ejection events, whereas the low-turbidity modules may be induced by increased downward transport of suspended sediment by sweep events. These hypotheses were supported by cross-spectral analysis performed on velocity and suspended sediment concentration time-series recorded at the site. An enhanced (negative) contribution of outward and inward interaction events to the Reynolds stress, compared to those reported in uniform BBLs, resulted in ‘ abnormally ’ low stress values.     

Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations

Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching , or even higher.In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES.This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3×10³Re3×10⁴, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re=10⁵ and 10⁶, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67% when Re is increased from Re=10³ to Re=10⁶, within the computation period, with the most significant increase taking place from Re=3×10³ to Re=10⁵.     

Hydrodynamic shape optimization by high fidelity CFD solver and Gaussian process based response surface method

A computational framework for hydrodynamic shape optimization of complex ship hull form is proposed and applied to improve the calm water performance of the KRISO Container Ship (KCS). The framework relies on three key features: a novel shape morphing method based on a combination of subdivision surfaces and free form deformations, a robust three dimensional viscous computational fluid dynamic solver based on the openFOAM open-source libraries and a Gaussian process-response surface method (GP-RSM) based on ordinary Kriging model which has been created to speed-up the evaluation of the quantity of interest (QoI) of the design process.The accuracy of the hydrodynamic solver is proven by comparing the obtained results against available experimental measurements. A preliminary sensitivity analysis on the mesh size has been carried out aiming at reducing the computational burden required by the CFD predictions. Three GP-RSMs have been trained relying on increasing number of hull designs. Each surrogate model has been cross-validated by both leave-one-out and k-fold techniques. The behaviours of these multi-dimensional surfaces have been analyzed in details by sampling the investigated design space with 10⁷ points according to a Full-Factorial algorithm, highlighting the regions of maximum deviation with respect to the resistance of the reference hull. The three optimum designs provided by the corresponding GP-RSM models have been verified by using high-fidelity CFD simulations with a refined mesh configuration. Calm water resistance, wave patterns and pressure distributions over the selected hull surfaces have been discussed in the light of the generated shape variations.