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

        ?

        Some notes on numerical simulation and error analyses of the attached turbulent cavitating flow by LES *

        2018-05-14 01:43:32XinpingLong龍新平YunLong龍云WentingWang王文婷HuaiyuCheng程懷玉BinJi季斌
        水動力學研究與進展 B輯 2018年2期

        Xin-ping Long (龍新平), Yun Long (龍云), Wen-ting Wang (王文婷), Huai-yu Cheng (程懷玉),Bin Ji (季斌)

        1. State Key Laboratory of Water Resources and Hydro power Engineering Science, School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China

        2. Key Laboratory of Jet Theory and New Technology of Hubei Province, Wuhan University, Wuhan 430072,China

        It is a great challenge for both experiment and numerical simulation to investigate the complicated transient cavitating flow. It is noted that numerical simulation of cavitating flow has achieved a remarkable progress in the last two decades. Li et al.[1,2]propose an improved all-speed Roe scheme for LES and obtain satisfactory results. The re-entrant jet[3]and shock wave[4]has attracted much attention. Peng et al.[5]observe the U-type fl ow structures around the hydrofoils in the cavitation tunnel. Ji et al.[6], Wang et al.[7], Long et al.[8]and Cheng et al.[9]study transient cavitating flow structure by high performance numerical methods and three-dimensional Lagrangian technology.

        All these researches enhance the understanding of the transient cavitating flow. However, less attention are paid to assessing the accuracy of these numerical results with systematic methodology and procedures. Verification and validation (V&V) is a basic procedure to evaluate the accuracy and reliability of numerical simulation, but the practical use of LES V&V is a big challenge[10]. Related studies of LES V&V are ongoing[11]. Freitag and Klein[12]have proposed a systematic grid and model variation method to estimate the numerical and modeling error.Xing[11]proposes a general framework for V&V of LES to get the numerical, modeling and their coupling errors. The investigations by Dutta and Xing[13]have achieved promising achievements.

        Inspired by the aforementioned studies, this letter carries out the large eddy simulation and error analysis of the attached turbulent cavitating flow around the Clark-Y hydrofoil. Comparisons between available experiments[14]and numerical results have been made. Quantitative analyses of LES errors are presented and cavitation influence on LES errors is also discussed.

        The chord lengthCof hydrofoil is 0.07 m and angle of attack is 8°. The length, width and height of the computational domain are 10.0C, 0.3Cand 2.7C, respectively. The boundary inlet velocity is set asU∞=10m/s , and outlet boundary is fixed at the static pressure pout, deriving from the cavitation number σ=0.8 (σ=(pout-pv)/(0.5ρU∞2)), wherepvis the saturated vapor pressure andρis the liquid density). Detailed information about the setup can see Ref. [15]. Case 1 is refined with the same mesh topology from Case 3, and Case 2 is the same as Case 1 except for the change of Smagorinsky constant.Detail information about mesh, model constant and calculations setup has been displayed in Table 1.

        Table 1 Mesh and model constant and calculation setup

        Herein, the LES errors estimation method, proposed by Freitag and Klein[12], has been employed.Some revised thoughts by Xing[11]have been adopted in this paper considering that the time step and grid spacing changing simultaneously. The numerical error and modeling error are estimated with power series expansions as a function ofh*and Δ:

        whereSCis the numerical benchmark,S1,S2,S3denotes the standards LES solution, the LES solution with a modified model contribution and solution on a coarser mesh, respectively.cNandcMare the unknown coefficient. The local spatial and temporal resolutionh*[11]is defined as

        where Δtis the time step andhis the local mesh size. Herehis equal to the filter width Δ. The model variation factor α=1.3. The time step and grid spacing share the same refinement ratior=1.6.pNandpMhas a large effect on the LES errors estimations.pN=1.7 andpM=1.5 were recommended by Dutta and Xing[16], whilepN=2 andpM=2/3 by Freitag and Klein[12]have been investigated and applied in some researches. Therefore,pN=2 andpM=2/3 are adopted in this paper.

        The numerical error and modeling error for case1 are corresponding to the first term and second term on the right hand sides of Eq. (1). The normalized averaged stream wise velocityu/U∞is chosen as the target variable to gain the numerical and modeling error. The velocity and cavitation data are obtained from the medial surface along the span wise location.

        Figure 1 shows the comparison of experimental and numerical snapshots of cavity shedding within one typical cycle. The experimental results are from Ref. [14]. Numerical results are used to represent the unsteady cavity evolution colored with vertical velocity. From Figs. 1(a)-1(d), it can be observed that the numerical results can capture the unsteady cavity evolution features, including the cavity inception,development and shedding along the hydrofoil suction face. The periodic cavity shedding behavior can lead to complex and turbulent flow structures. This can cause huge difficulties to accurately simulate the unsteady cavitating flow even with well performance numerical methods. Owing to the extensive use of LES in academic research in recent years, the reliability or risk assessment has become an essential step for LES results in cavitation. At present, it has the possibility to carry out quantitative LES errors analysis in cavitation simulation. Verification and validation of LES will be a necessary procedure in the future.

        Fig.1 (Color online)Timeevolution of cavity shedding in one typical cycle for Case1 (Left: Experimental results[14].Right: Cavitation patterns by iso-surface of αv=0.1)

        Figure 2 presents the average values versus different flow time about the normalized stream wise velocityu/U∞at monitoring point A (x/c=0.2,y/c=0.09) and point B (x/c=0.8,y/c=0.16). It can be seen that the variations ofu/U∞at the ending of statistical flow time are clearly smaller compared to the beginning of flow time. This ensures that the selection of final averaged velocity in this paper has almost no effect on obtaining the LES errors.

        Fig.2 Average values versus different flow time about the normalized stream wise velocities at monitoring points

        The magnitudes of numerical error and modeling error for Case 1 are shown in Fig.3 for the unsteady cavitating flow. From the results, it can be seen that the modeling error and numerical error are almost on the same order of magnitude in most situations, while the modeling error is often a bit larger if compared to numerical error. For each chord wise location in Fig.3,the numerical error and modeling error sometimes partially have opposite signs and they can offset each other. In Fig.3(a), the errors are large and squeezed in very little space at vertical direction, but sharply decrease to almost zero. The cavitation occurs generally around this narrow location, and the accuracy of the inner cavitating flow structure at this narrow location is a great difficulty for numerical simulation,also for experiments. This might be in charge of the error distribution in Fig.3(a). From Figs. 3(b)-3(e),the numerical and modeling errors start to oscillate along the vertical direction. The errors magnitudes and fluctuation levels gradually become larger fromx/c=0.4 tox/c=1.0. It has the similar tendency that the errors decrease to almost zero when the distance is far away from foil surface. The cavitation area becomes thicker from Figs. 3(b), 3(e). The velocity fluctuate sviolently due to the periodic cavity shedding,which can cause much large scale vortexstructure. It has a huge impact on the cavitation simulation. This might result in the errors changing in a large range. From above analyses, it can be con-cluded that the LES errors are significantly influenced by the unsteady cavitation.

        Fig.3 The numerical errors and modeling errors of averaged u/ U∞ at five chord wise locations for Case 1

        Large eddy simulation and error analyses of the attached turbulent cavitating flow around the Clark-Y hydrofoil has been made in this letter. Quantitative LES errors analyses of cavitation simulation are carried out and cavitation influence on LES errors is also discussed. The main conclusions can be summarized as follows:

        (1) The numerical results can capture the periodic cavity shedding behavior and they are in well agreement with the experiments. The shedding of cavity can lead to complex and turbulent flow structures,which results in great difficulties to model the details of the turbulent cavitating flow even with well performance numerical methods.

        (2) In this simulation of the transient cavitating flow, the modeling error and numerical error are on the same order of magnitude in most situations, while the modeling error is often a little bit larger than numerical error. The numerical error and modeling error sometimes have the opposite sign and they can partially offset each other.

        (3) The numerical error and modeling error are significantly influenced by the cavitation. The unsteady cavitation can extend the magnitudes and oscillation levels of LES numerical error and modeling error.

        It still has a long way to go for the quantitative LES reliability estimation with V&V. Owing to the extensive use of LES in academic research in recent years, the LES reliability estimation should become an essential part to LES results in cavitation. Due to that Smagorinsky constant cannot change arbitrarily, the LES error estimation method adopted in this paper needs more investigations. The three-equation method proposed by Dutta and Xing[16]may give more promising results. The LES error estimation applied in the unsteady cavitation simulation with the guideline of Ref. [16] will be the focus of our future investigations.

        [1] Li X. S., Li X. L. All-speed Roe scheme for the large eddy simulation of homogeneous decaying turbulence [J].International Journal of Computational Fluid Dynamics,2016, 30(1): 69-78.

        [2] Li X. S., Xu J. Z., Gu C. W. Preconditioning method and engineering application of large eddy simulation [J].Science in ChinaSeries G-Physics, Mechanics and Astronomy, 2008, 51(6): 667-677.

        [3] Wu Q., Huang B., Wang G. et al. The transient characteristics of cloud cavitating flow over a flexible hydrofoil[J].International Journal of Multiphase Flow, 2018, 99:162-173.

        [4] Cui P., Zhang A. M., Wang S. et al. Ice breaking by a collapsing bubble [J].Journal of Fluid Mechanics, 2018,841: 287-309.

        [5] Peng X. X., Ji B., Cao Y. et al. Combined experimental observation and numerical simulation of the cloud cavitation with U-type flow structures on hydrofoils [J].International Journal of Multiphase Flow, 2016, 79:10-22.

        [6] Ji B., Luo X., Arndt R. E. A. et al. Large eddy simulation and theoretical investigations of the transient cavitating vortical flow structure around a NACA66 hydrofoil [J].International Journal of Multiphase Flow, 2015, 68:121-134.

        [7] Wang Y., Xu C., Wu X. et al. Ventilated cloud cavitating flow around a blunt body close to the free surface [J].Physical Review Fluids, 2017, 2(8): 084303.

        [8] Long X., Cheng H., Ji B. et al. Large eddy simulation and Euler-Lagrangian coupling investigation of the transient cavitating turbulent flow around a twisted hydrofoil [J].International Journal of Multiphase Flow, 2018, 100:41-56.

        [9] Cheng H. Y., Long X. P., Ji B. et al. 3-D Lagrangian based investigations of the time-dependent cloud cavitating flows around a Clark-Y hydrofoil with special emphasis on shedding process analysis [J].Journal of Hydrodynamics,2018, 30(1): 122-130.

        [10] Long Y., Long X., Ji B. et al. Verification and validation of URANS simulations of the turbulent cavitating flow around the hydrofoil [J].Journal of Hydrodynamics, 2017,29(4): 610-620.

        [11] Xing T. A general framework for verification and validation of large eddy simulations [J].Journal of Hydrodynamics, 2015, 27(2): 163-175.

        [12] Freitag M., Klein M. An improved method to assess the quality of large eddy simulations in the context of implicit filtering [J].Journal of Turbulence, 2006, 7(40): 1-11.

        [13] Dutta R., Xing T. Quantitative solution verification of large eddy simulation of channel flow [C].Proceedings of the 2nd Thermal and Fluid Engineering Conference and 4th International Workshop on Heat Transfer, Las Vegas,USA, 2017.

        [14] Huang B., Young Y. L., Wang G. et al. Combined experimental and computational investigation of unsteady structure of sheet/cloud cavitation [J].Journal of Fluids Engineering, 2013, 135(7): 071301.

        [15] Ji B., Long Y., Long X. P. et al. Large eddy simulation of turbulent attached cavitating flow with special emphasis on large scale structures of the hydrofoil wake and turbulence-cavitation interactions [J].Journalof Hydrodynamics,2017, 29(1): 27-39.

        [16] Dutta R., Xing T. Five-equation and robust three-equation methods for solution verification of large eddy simulation[J].Journal of Hydrodynamics,2018, 30(1): 23-33.

        91热爆在线精品| 国产激情内射在线影院| 四虎影院在线观看| 国产精品不卡无码AV在线播放| 久久精品国产亚洲av日韩一| www夜插内射视频网站| 97人妻碰碰视频免费上线| 欧美人与动zozo| 亚洲一区二区女优视频| 在线观看一级黄片天堂| 精品成人av一区二区三区| 亚洲最大在线精品| 亚洲一区二区三区麻豆| 亚洲成av人片一区二区密柚| 国产成年无码v片在线| 四虎成人精品国产一区a| 国产偷闻女邻居av在线观看| 噜噜综合亚洲av中文无码| 久久久久国产一区二区三区| 日本岛国精品中文字幕| 日本加勒比精品一区二区视频 | 亚洲精品女人天堂av麻| 免费国产线观看免费观看| 亚洲av无码一区二区二三区 | 国产女合集小岁9三部| 夫妻一起自拍内射小视频| 91免费国产| 国产一区二区av在线免费观看| 亚洲精品成人片在线观看精品字幕| 日日摸夜夜添夜夜添无码免费视频 | 91国视频| 日本红怡院东京热加勒比| 亚洲一二三区免费视频| 正在播放东北夫妻内射| 欧美中文字幕在线看| 国产黄色一区二区三区av| 少妇愉情理伦片| 欧美日韩中文制服有码| 日本在线观看一区二区视频| 一边捏奶头一边高潮视频| aaaaa级少妇高潮大片免费看|