A K-εTWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.

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.

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.

Theoretical and numerical study of hydraulic characteristics of orifice energy dissipator

Different factors affecting the efficiency of the orifice energy dissipator were investigated based on a series of theoretical analyses and numerical simulations. The main factors investigated by dimension analysis were identified, including the Reynolds number （Re）, the ratio of the orifice diameter to the inner diameter of the pipe （ did ）, and the ratio of distances between orifices to the inner diameter of the pipe （ LID ）. Then, numerical simulations were conducted with a k-ε two-equation turbulence model. The calculation results show the following： Hydraulic characteristics change dramatically as flow passes through the orifice, with abruptly increasing velocity and turbulent energy, and decreasing pressure. The turbulent energy appears to be low in the middle and high near the pipe wall. For the energy dissipation setup with only one orifice, when Re is smaller than 105, the orifice energy dissipation coefficient K increases rapidly with the increase of Re. When Re is larger than l05, K gradually stabilizes. As diD increases, K and the length of the recirculation region L1 show similar variation patterns, which inversely vary with diD. The function curves can be approximated as straight lines. For the energy dissipation model with two orifices, because of different incoming flows at different orifices, the energy dissipation coefficient of the second orifice （K2） is smaller than that of the first. If LID is less than 5, the K value of the LID model, depending on the variation of/（2, increases with the spacing between two orifices L, and an orifice cannot fulfill its energy dissipation function. If LID is greater than 5, K2 tends to be steady; thus, the K value of the LID model gradually stabilizes. Then, the flow fully develops, and L has almost no impact on the value of K.

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.

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.