修正的SST k-ω模型在云状空化流动计算中的应用研究

采用基于Rayliegh-Plesset方程的空化模型及实验结果,评价了修正的SST k-ω湍流模型在云状空化流动计算中的应用。采用不同的修正系数,分别计算了绕Clar k-y型水翼的云状空化流动,获得了随时间变化的空化形态和升、阻力等流场及动力特性。通过与实验结果的对比表明,修正后的模型可以更准确地捕捉云状空化区域的空穴形态和空泡脱落的非定常细节;密度函数中指数n的选取对计算所得的空穴长度、升阻力和主要频谱分布均有影响,并给出了n的合理取值范围,为后续水工模型的空化数值计算提供了参考依据。

Numerical Simulations of Cavitation Flows around Clark-Y Hydrofoil

Cavitation is a complex flow phenomenon including unsteady characteristics, turbulence, gas-liquid two-phase flow. This paper provides a numerical investigation on comparing the simulation performance of three different models in OpenFOAM-Merkle model, Kunz model and Schnerr-Sauer model, which is helpful for understanding the cavitation flow. Considering the influence of vapor-liquid mixing density on turbulent viscous coefficient, the modified SST k-ω model is adopted in this paper to increase the computing reliability. The InterPhaseChangeFoam solver is utilized to simulate the two-dimensional cavitation flow of the Clark-Y hydrofoil with three cavitation models. The hydrodynamic performance including lift coefficient, drag coefficient and cavitation flow shape of the hydrofoil is analyzed. Through the comparison of the numerical results and experimental data, it is found that the Schnerr-Sauer model can get the most accurate results among the three models. And from the simulation point of water and water vapor mixing, the Merkle model has the best water and water vapor mixing simulation.

TWO MODIFICATORY K-ε TURBULENCE MODELS FOR TURBULENT SWIRLING FLOWS

Since the standard K-ε model used to predict the strongly swirling flowleads to a large deviation from experimental results, it is necessary to introduce modification tothe standard K-ε model. Based on the algebraic Reynolds stress model and Bradshaw's turbulentlength scale modification conception, we present two modified K-ε models. To investigate thebehaviour of the modified turbulence models, they are used to predict two representative turbulentswirling flows. The computational results, after compared with the experimental data, show that themodified K-ε models substantially improve the prediction of the standard K-ε model for theturbulent swirling flows.

基于黏性修正SST k-ω模型的水翼空化流数值模拟计算

该研究利用开源平台Open FOAM中的interPhaseChangeFOAM求解器,应用SchnerrSauer空化模型,并结合黏性修正后的SSTk-ω模型,对二维水翼NACA0015及三维扭翼TWIST-11N进行了空化流数值模拟,从升力系数和空化脱落频率等方面对计算结果进行比较,验证了修正SST k-ω模型的可靠性,并对空化脱落机理进行了分析。首先对二维水翼NACA0015,应用了基于LES方法的Smagorinsky模型、基于RANS方法的标准SST k-ω模型和黏性修正SST k-ω模型三种不同的湍流模型,结果表明修正SSTk-ω模型能准确预报二维非定常空化的脱落频率,并且发现水翼尾部的回射流是空化脱落的主要原因。之后将修正SST k-ω模型应用于三维水翼TWIST-11N,发现应用该模型依旧能准确模拟出三维水翼的非定常空化,结果和Ji等人的结果十分接近。通过进一步对三维空化流的流场进行分析,发现三维空化流会出现主脱落及次脱落共两次空化脱落,分别由尾部回射流和侧向流导致。

CFD Prediction of the Airflow at a Large-Scale Wind Farm above a Steep, Three-Dimensional Escarpment

The Duogu Wind Farm, China Huadian Group Corporation’s first wind project in Yunnan, China, has been approved by the Provincial Development and Reform Commission. The acquired site is in Mengzi, in the south-east of Yunnan Province. The developer has deployed thirty-three 1.5 MW turbines in this wind farm (49.5 MW), and the total cost of construction has been estimated to be CNY449.7 million ($69.61 million). The present study compared the prediction accuracy of two CFD software packages for simulating flow over an escarpment with a steep slope. The two software packages were: 1) Open FOAM (Turbulence model: SST k-ω RANS), which is a free, open source CFD software package developed by Open CFD Ltd at the ESI Group and distributed by the Open FOAM Foundation and 2) RIAM-COMPACT (Turbulence model: Standard Smagorinsky LES), which has been developed by the lead author of the present paper. Generally good agreement was obtained between the results from the simulations with Open FOAM and RIAM-COMPACT.

Experimental and Numerical Investigations of the Hydrodynamic Characteristics, Twine Deformation, and Flow Field Around the Netting Structure Composed of Two Types of Twine Materials for Midwater Trawls

Nettings are complex flexible structures used in various fisheries.Understanding the hydrodynamic characteristics,de-formation,and the flow field around nettings is important to design successful fishing gear.This study investigated the hydrodynamic characteristics and deformation of five nettings made of polyethylene and nylon materials in different attack angles through numeri-cal simulation and physical model experiment.The numerical model was based on the one-way coupling between computational fluid dynamics(CFD)and large deflection nonlinear structural models.Navier-Stokes equations were solved using the finite volume ap-proach,the flow was described using the k-ωshear stress turbulent model,and the large deflection structural dynamic equation was derived using a finite element approach to understand the netting deformation and nodal displacement.The porous media model was chosen to model the nettings in the CFD solver.Numerical data were compared with the experimental results of the physical model to validate the numerical models.Results showed that the numerical data were compatible with the experimental data with an average relative error of 2.34%,3.40%,6.50%,and 5.80%in the normal drag coefficients,parallel drag coefficients,inclined drag coefficients,and inclined lift coefficients,respectively.The hydrodynamic forces of the polyethylene and nylon nettings decreased by approxi-mately 52.56%and 66.66%,respectively,with decreasing net solidity.The drag and lift coefficients of the nylon netting were appro-ximately 17.15%and 6.72%lower than those of the polyethylene netting.A spatial development of turbulent flow occurred around the netting because of the netting wake.However,the flow velocity reduction downstream from the netting in the wake region in-creased with increasing attack angle and net solidity.In addition,the deformation,stress,and strain on each netting increased with in-creasing solidity ratio.

3D Simulation of Storm Surge Disaster Based on Scenario Analysis

The occurrence of storm surge disaster is often accompanied with floodplain, overflow, dike breach and other complex phenomena, while current studies on storm surge flooding are more concentrated on the 1D/2D numerical simulation of single disaster scenario(floodplain, overflow or dike breach), ignoring the composite effects of various phenomena. Therefore, considering the uncertainty in the disaster process of storm surge, scenario analysis was firstly proposed to identify the composite disaster scenario including multiple phenomena by analyzing key driving forces, building scenario matrix and deducing situation logic. Secondly, by combining the advantages of k-ω and k-ε models in the wall treatment, a shear stress transmission k-ω model coupled with VOF was proposed to simulate the 3D flood routing for storm surge disaster. Thirdly, risk degree was introduced to make the risk analysis of storm surge disaster. Finally, based on the scenario analysis, four scenarios with different storm surge intensity(100-year and 200-year frequency) were identified in Tianjin Binhai New Area. Then, 3D numerical simulation and risk map were made for the case.

An optimization method on runner blades in bulb turbine based on CFD analysis

In this paper an optimization method of the runner blades in a bulb turbine based on CFD analysis is proposed.In the method the main scales of the turbine including guide vane,runner and draft are maintained.Only the runner blades are modified based on the present method.In the optimization method the runner blade is expressed by spline surface with a gather of coordinate points.The B-spline curve is used to keep the modified blades smooth.In order to make the blade optimization simple and ef- ficient,one of the coordinates is fixed and only the angles of the points are changed according to different modification purposes.Three main optimization principles based on flow diagnosis are presented here.These three principles are all based on the CFD analysis of the internal flow in bulb turbine.For the purpose of method verification,the optimization method is used in a model bulb turbine.A three dimensional steady turbulent computation is carried out through the whole passage including the bulb body,guide vanes,runner and draft tube of the bulb turbine under seven different work conditions.An SST k-ωturbulence model is used during the CFD analysis and the performance of the turbine can be achieved.The runner blade is optimized according to the three optimization principles based on flow diagnosis.The CFD analysis is conducted again on the optimized turbine and another modification is needed if the new turbine can’t satisfy the required performance.Comparison of the computational results between the original turbine and an optimized one indicates that the optimization method is practical and does improve the performance of the bulb turbine.

Numerical simulations of turbulent flows in aeroramp injector/gas-pilot flame scramjet

To uncover the internal flow characteristics in an ethylene-fueled aeroramp injector/gaspilot（ARI/G-P）flame scramjet,a Reynolds-averaged Navier-Stokes（RANS）solver is constructed under a hybrid polyhedral cell finite volume frame.The shear stress transport（SST）k-x model is used to predict the turbulence,while the Overmann’s compressibility corrected laminar flamelet model is adopted to simulate the turbulent combustion.Nonreactive computations for Case 1（G-P jet on）,Case 2（ARI jets on）,and Case 3（both ARI and G-P jets on）were conducted to analyze the mixing mechanism,while reactive Cases 4–7 at equivalent ratios of 0.380,0.278,0.199 and0.167 respectively were calculated to investigate the flame structure and combustion modes.The numerical results are compared well to those of the experiments.It is shown that the G-P jet plays significant role in both the fuel/air mixing and flame holding processes;the combustion for the four reactive cases takes place intensively in the regions downstream of the ARI/G-P unit;Cases 4 and 5are under subsonic combustion mode,whereas Cases 6 and 7 are mode transition critical and supersonic combustion cases,respectively;the mode transition equivalent ratio is approximately 0.20.

Investigation of the fluid flow in an isolated rotor-stator system with a peripheral opening

This paper deals with an experimental, theoretical and numerical study of a turbulent flow with separated boundary layers between a rotor and a stator. The system is not subjected to any superimposed radial flow. The periphery of the cavity is opened to the atmosphere so that the solid body rotation for infinite discs is not always observed. Emphasis was placed on develop- ment of an asymptotic approach and a step-by-step method to compute the radial distribution of the core swirl ratio and the static pressure on the stator side. The theory also includes the radial and axial velocities in the core region. The numerical simulation has been conducted with the commercial CFD code Fluent 6.1. The k- SST turbulence model is used, with the assumption of 2D-axisymmetric and steady flow. CFD validations have been performed by comparison of the numerical results with the corresponding theoretical results. Numerical and experimental results are in good agreement with analytical solutions.