亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Numerical simulation of flow through circular array of cylinders using porous media approach with non-constant local inertial resistance coefficient*

        2017-03-09 09:09:58JieminZhan詹杰民WenqingHu胡文清WenhaoCai蔡文豪YejunGong龔也君ChiwaiLi李志偉

        Jie-min Zhan (詹杰民), Wen-qing Hu (胡文清), Wen-hao Cai (蔡文豪), Ye- jun Gong (龔也君), Chi-wai Li (李志偉)

        1.Department of Applied Mechanics and Engineering, School of Engineering, Sun Yat-Sen University, Guangzhou 510275, China, E-mail: cejmzhan@vip.163.com

        2.Department of Civil and Environmental Engineering, Hong Kong Polytechnic University, Hong Kong, China

        (Received November 29, 2016, Revised December 15, 2016)

        Numerical simulation of flow through circular array of cylinders using porous media approach with non-constant local inertial resistance coefficient*

        Jie-min Zhan (詹杰民)1, Wen-qing Hu (胡文清)1, Wen-hao Cai (蔡文豪)1, Ye- jun Gong (龔也君)1, Chi-wai Li (李志偉)2

        1.Department of Applied Mechanics and Engineering, School of Engineering, Sun Yat-Sen University, Guangzhou 510275, China, E-mail: cejmzhan@vip.163.com

        2.Department of Civil and Environmental Engineering, Hong Kong Polytechnic University, Hong Kong, China

        (Received November 29, 2016, Revised December 15, 2016)

        Aquatic vegetation zone is now receiving an increasing attention as an effective way to protect the shorelines and riverbeds. To simulate the flow through the vegetation zone, the vegetation elements are often simplified as equidistant rigid cylinders, and in the whole zone, the porous media approach can be applied. In this study, a non-constant inertial resistance coefficient is introduced to model the unevenly distribution of the drag forces on the cylinders, and an improved porous media approach is applied to one circular array of cylinders positioned in a 2-D flume. The calculated velocity profile is consistent with the experimental data.

        Vegetation zone, cylinder array, porous media, inertial resistance coefficient

        The aquatic vegetation is widely researched due to the significant influence on ecosystem and hydraulics. On the aspect of laboratory experiments, a series of classic experiments were conducted by Ghisalberti and Nepf[1]and Nezu and Onitsuka[2]. In their studies, the vegetation is simulated by rigid bodies, and the whole vegetation zone is regarded as an array of solid cylinders. The tracking dye is used to track the flow trajectory, showing the complex flow pattern, turbulence and diffusion in the flows through the vegetation.

        For the numerical simulation, two main methods[3]are widely employed to simulate the flow through the vegetation. The first method, which resolves the flow around each vegetation unit, is called the multi-body model[4]. Each vegetation unit is generallysimplified as a cylinder and each one of them is taken into consideration in the model. The second method, the porous model[5], treats the entire vegetation zone as a porous zone and the impact of the vegetation on the flow is governed by the relevant parameters of the porous zone. Compared with the porous model, the multi-body model considers the influence of each single vegetation unit, therefore, is more accurate and time consuming. Relatively speaking, as the porous model regards the vegetation as one porous zone, the computational efficiency is better, therefore, it is more widely used in the real engineering applications.

        For the flow of incompressible, homogeneous and viscous fluid, the governing Navier-Stokes equations can be written as

        wherevis the velocity vector,ρis the fluid density,pis the pressure,νis the kinematic viscosity,Fbrepresents the force acting on the unit mass and it is generally related to the gravitational accelerationg.

        A hybrid RANS/LES turbulence model, the scale adaptive simulation (SAS) model, is employed. The SAS model is based on the SST model, which behaves very much like thek-ωmodel away from the wall and serves to control the turbulence level in the near wall region. But the length scale produced by the SST model is too large. To avoid this shortcoming, in the SAS model, the von Karman length scale is introduced into the turbulence scale equation, such that the turbulence length-scale can be resolved correctly. The transport equations for the SAS model are defined as[8]

        The additional source termQSASis defined as

        The porous media are modeled by adding one source term to the momentum equations. The source term is composed of two parts: a viscous loss term and an inertial loss term. For homogeneous porous media, the source term is defined as

        whereαis the permeability andC2is the inertial resistance factor. This momentum sink contributes to the pressure gradient in the porous cell.

        In the convectional porous model, the inertial resistance factor is assumed to be constant. However, it is in conflict with the real situation, especially, when the size of the porous region is large. In this study, two assumptions are made in the porous region of the cylinder array: the local inertial resistance is proportional to the local pressure drop, the pressure head loss along the flow direction can be neglected. With these assumptions, the local inertial resistance can be redefined as

        whereC2iis the inertial resistance coefficient in celli, and2is the averaged coefficient value in the area of the circular patches.Ptotalis the total pressure on the inlet boundary,Piis the local cell static pressure,is the time-averaged drag force of the circular patches,is the velocity magnitude,vxiis the velocity in the flow direction,Siis the cross sectional area of the local cell normal to the flow direction,Apatchis the area of the circular patches,∑Aiis the summation of the areas of the local cells, and Δniis the thickness of the local cell.Ptotalis a constant when the inlet velocity remains the same.

        In order to validate the feasibility and the rationality of the improved porous method with the adoption of the non-constant local inertial resistance coefficient, the calculated velocity field is compared to the measurements of the experiments conducted by Zong et al.[9]. The computational domain is shown in Fig.1, where the circular vegetation zone with diameterD= 0.22 m is centered atx=0. The diameter of each circular dot (vegetation unit) isd=0.006 m. The flume’s width is 1.2 m, and the water depth is 0.133 m. Additionally, two monitoring lines (marked in red in Fig.1) are set at the appropriate positions to obtain the relative velocity profiles measured in the experiments. The time-averaged horizontal velocityuis measured at liney=0, while the time-averaged vertical velocityis measured at liney=D/2, respectively. The dimensionless constantΦ=N(d/D)2(whereNrepresents the number of circular patches) is defined as the ratio of the sum of all circular patch areas to the total vegetation region area, representing the plant density of the vegetation. In the porous model, the value ofΦrepresents the porosity, andΦ=0.03both for the multi-body and the porous models in this study.

        Fig.1 The schematic diagram of computational model

        Fig.2 (Color online) Contour maps of the local inertial resistance coefficients

        Figure 2 is the contour maps of the local inertial resistance coefficients. It is calculated based on the SAS multi-body model. As is described above, the porous model has a smaller size, that requires an interpolation of the C2 distribution inside the porous zone.

        Figure 3 shows the time-averaged horizontal velocity profile along the liney=0. Similarly, Fig.4 shows the time-averaged vertical velocity profile along the liney=D/2.U∞is the inlet velocity,Dis the diameter of the circular patches. Both results of the multi-body model and the porous model fit well with the experimental data outside the circular region of the cylinder array. The multi-body model performs better in the prediction of the vertical velocity than the porous model away from the cylinders. Moreover, the multi-body model predicts the fluctuation of the flow velocity inside the region of the cylinder array, caused by impacting effect of the flow through the dense array of cylinders. Note that the distribution of thecylinders is unknown in the experiment, so the measured velocity in the region of the cylinder array is not shown in Figs.3, 4.

        Fig.3 (Color online) Time-averaged horizontal velocity profiles alongy=0

        Fig.4 (Color online) Time-averaged vertical velocity profiles alongy=D/2

        The vorticity fields simulated by the multi-body and porous models are shown in Figs.5, 6, where Fig.6 shows the details of the vorticity field around the cylinder array. As expected, the porous model can not capture the small scale vortices in the near-field, but still gives some details of the disturbance inside the porous zone, as shown in Fig.6. However, in the farfield, the porous model gives a similar vorticity field as the multi-body model.

        Fig.5 (Color online) Vorticity fields of multi-body (a) and porous (b) models

        Fig.6 (Color online) Contours of transient vorticity magnitude near porous zone

        Compared to the multi-body method, the improved porous model can give a proper result over the flow field while the simulation cost is much less than the multi-body model, such that it is suitable for the large-scale engineering applications.

        [1] Ghisalberti M., Nepf H. M. Mixing layers and coherent structures in vegetated aquatic flows [J].Journal of Geophysical Research: Oceans, 2002, 107(2): 1-11.

        [2] Nezu I., Onitsuka K. Turbulent structures in partly vegetated open-channel flows with LDA and PIV measurements [J].Journal of Hydraulic Research, 2001, 39(6): 629-642.

        [3] Yu L. H., Zhan J. M., Li Y. S. Numerical investigation of drag force on flow through circular array of cylinders [J].Journal of Hydrodynamics, 2013, 25(3): 330-338.

        [4] Zhang H., YANG J. M., XIAO L. F. et al. Large-eddy simulation of the flow past both finite and infinite circular cylinders atRe=3900[J].Journal of Hydrodynamics, 2015, 27(2): 195-203.

        [5] Su X., Li C. W. Large eddy simulation of free surface turbulent flow in partly vegetated open channels [J].International Journal for Numerical Methods in Fluids, 2002, 39(10): 919-937.

        [6] Yu L. H., Zhan J. M., Li Y. S. Numerical simulation of flow through circular array of cylinders using multi-body and porous models [J].Coastal Engineering Journal, 2014, 56(3): 1450014.

        [7] Wang Q., Li M., Xu S. Experimental study on vortex induced vibration (VIV) of a wide-D-section cylinder in a cross flow [J].Theoretical and Applied Mechanics Letters, 2015, 5(1): 39-44.

        [8] Menter F. R., Egorov Y. The scale-adaptive simulation method for unsteady turbulent flow predictions. Part 1: Theory and model description [J].Flow, Turbulence and Combustion, 2010, 85(1): 113-138.

        [9] Zong L., Nepf H. Vortex development behind a finite porous obstruction in a channel [J].Journal of Fluid Mechanics, 2012, 691: 368-391.

        * Project supported by the Special Fund of Marine-Fishery Science-Technology Extension in Guangdong Province (Grant No. A201401B08).

        Biography:Jie-min Zhan (1963-), Male, Ph. D., Professor

        Ye-jun Gong, E-mail:yejungong@126.com

        国产一区二区三区亚洲avv| 黑人巨大无码中文字幕无码| 亚洲日韩乱码中文无码蜜桃臀| 国产精品开放小视频| 国产成人精品三级在线影院| 中文字幕高清一区二区| 在线观看视频日本一区二区| 国产一区二区女内射| 九九在线视频| 国产一区精品二区三区四区| 亚洲国产精品无码久久一区二区| 人人妻人人狠人人爽天天综合网| 日本熟妇人妻xxxxx视频| 亚洲天堂在线播放| 日本一区二区三区看片| 亚洲六月丁香色婷婷综合久久| 免费国产在线精品一区二区三区免| 少妇被粗大的猛进出69影院| а√天堂资源8在线官网在线| 久久亚洲成a人片| 国产三级国产精品国产专播| 色中文字幕在线观看视频| 成人a级视频在线观看| 久久国产精品波多野结衣av| 国产一级淫片免费播放电影| 日本免费一区精品推荐| 日本一二三区在线观看视频| 狠狠色噜噜狠狠狠777米奇小说 | 美女爽好多水快进来视频| 中文字幕中乱码一区无线精品| 久久天堂精品一区二区三区四区| 中文字幕无线码免费人妻| 豆国产95在线 | 亚洲| 亚洲av色香蕉一区二区三区蜜桃 | 人妻少妇无乱码中文字幕| 中文国产乱码在线人妻一区二区| 国产对白国语对白| 比比资源先锋影音网| 91超碰在线观看免费| 日本在线观看一区二区三区视频| 久久99精品久久久大学生|