The motion of particle clouds formed by dumping dredged material into quiescent waters is experimentally and numerically studied. In the numerical model, the particle phase is modeled by the dispersion model, and turb...The motion of particle clouds formed by dumping dredged material into quiescent waters is experimentally and numerically studied. In the numerical model, the particle phase is modeled by the dispersion model, and turbulence is calculated by the large eddy simulation. The governing equations, including the filtered Navier-Stokes equations and mass transport equation, are solved based on the operator-splitting algorithm and an implicit cubic spline interpolation scheme. The eddy viscosity is evaluated by the modified Smagorinsky model including the buoyancy term. Comparisons of main flow characteristics, including shape, size, average density excess, moving speed and the amount of particles deposited on the bed, between experimental and computational results show that the numerical model well predicts the motion of the cloud from the falling to spreading stage. The effects of silt-fence on the motion of the particle cloud are also investigated.展开更多
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 ...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.展开更多
The most common method to determine the coefficient of Smagorinsky model now is to employ the Germano identity,however it is too complex and expensive in numerical calculation. In this letter we propose a new dynamic ...The most common method to determine the coefficient of Smagorinsky model now is to employ the Germano identity,however it is too complex and expensive in numerical calculation. In this letter we propose a new dynamic formula for determining the coefficient,which is based on the Kolmogorov equation of filtered velocity in physical space.The simplified formula is quite easy to implement in calculation.It is then verified in both homogeneous isotropic turbulence and wall-bounded turbulence by A Priori and A Posteriori tests.展开更多
In this work, a discrete particle model （DPM） was applied to investigate the dynamic characteristics in a gas-solid bubbling fluidized bed of binary solid particles. The solid phase was simulated by the hard- sphere...In this work, a discrete particle model （DPM） was applied to investigate the dynamic characteristics in a gas-solid bubbling fluidized bed of binary solid particles. The solid phase was simulated by the hard- sphere discrete particle model. The large eddy simulation （LES） method was used to simulate the gas phase. To improve the accuracy of the simulation, an improved sub-grid scale （SGS） model in the LES method was also applied. The mutative Smagorinsky constant case was compared with the previously published experimental data. The simulation by the mutative Smagorinsky constant model exhibited better agreement with the experimental data than that by the common invariant Smagorinsky constant model. Various restitution coefficients and different compositions of binary solids were investigated to determine their influences on the rotation characteristics and granular temperatures of the particles. The particle translational and rotational characteristic distributions were related to certain simulation parameters.展开更多
In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds...In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds number fluid flow applications, a subgrid turbulence model for LBM was introduced based on standard Smagorinsky subgrid model and Lattice Bhatnagar-Gross-Krook (LBGK) model. The subgrid LBGK model was subsequently used to simulate the two-dimensional driven cavity flow at high Reynolds numbers. The simulation results including distribution of stream lines, dimensionless velocities distribution, values of stream function, as well as location of vertex center, were compared with benchmark solutions, with satisfactory agreements.展开更多
A Large Eddy Simulation (LES) technique was applied to solve the turbulentchannel flow for Re_τ = 150 . Three types of turbulence models are employed, such as theSmagorinsky model, the Dynamic Sub-Grid Scale(SGS) mod...A Large Eddy Simulation (LES) technique was applied to solve the turbulentchannel flow for Re_τ = 150 . Three types of turbulence models are employed, such as theSmagorinsky model, the Dynamic Sub-Grid Scale(SGS) model and the Generalized Normal Stress (GNS)model. The simulated data in time series for the LES were averaged in both time and space to carryout the statistical analysis. Results of LES were compared with that of a DNS. As an application, aLES technique was used for 2D body in order to check the validation by investigating the turbulentvortical motion around the afterbody with a slant angle.展开更多
A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution di...A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.展开更多
Parallel implementation of algorithm of numerical solution of Navier-Stokes equations for large eddy simulation (LES) of turbulence is presented in this research. The Dynamic Smagorinsky model is applied for sub-grid ...Parallel implementation of algorithm of numerical solution of Navier-Stokes equations for large eddy simulation (LES) of turbulence is presented in this research. The Dynamic Smagorinsky model is applied for sub-grid simulation of turbulence. The numerical algorithm was worked out using a scheme of splitting on physical parameters. At the first stage it is supposed that carrying over movement amount takes place only due to convection and diffusion. Intermediate field of velocity is determined by method of fractional steps by using Thomas algorithm (tridiaginal matrix algorithm). At the second stage found intermediate field of velocity is used for determination of the field of pressure. Three dimensional Poisson equation for the field of pressure is solved using upper relaxation method. Moreover various ways of geometrical decomposition for parallel numerical solution of three dimensional Poisson equations are investigated.展开更多
When the sediment and the dissolved matter laden in the river meet a clear water in reservoirs, the turbid water will plunge and spread into the clear water, forming the turbidity current and influencing the water qua...When the sediment and the dissolved matter laden in the river meet a clear water in reservoirs, the turbid water will plunge and spread into the clear water, forming the turbidity current and influencing the water quality and the life of the reservoir. Due to the unsteady nature of the flood, the turbidity current is unsteady. In the present study, we use the MIKE 3 computational fluid dynamics code to simulate continuous and discontinuous turbidity currents on a flat slope. With the model used by us, the turbulence is divided into two parts: the horizontal turbulence and the vertical turbulence, which are separately modeled by the Smagorinsky model and our model to capture the anisotropic turbulence. In this model, the sediment settling and deposition are considered. The simulation results concerning the flume water surface level, the front velocity and sediment concentration profiles are found consistent with the experimental data, particularly, for the sediment concentration profiles with an absolute mean error of 0.026 kg/m3and the root mean square error of 0.046 kg/m3. This finding suggests that this model can be used to well predict the turbidity current on the flat slope.展开更多
This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow...This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow using two different nonlinear eddy viscosity functional forms to capture anisotropic dissipation mechanism,namely,(i)the Smagorinsky model using the local strain rate field,and(ii)the Leith model using the gradient of the vorticity field.The proposed models are applied to the one-layer and two-layer wind-driven quasigeostrophic ocean circulation problems,which are standard prototypes of more realistic ocean dynamics.Results show that both models capture the quasi-stationary ocean dynamics and provide the physical level of eddy viscosity distribution without using any a priori estimation.However,it is found that slightly less dissipative results can be obtained by using the dynamic Leith model.Two-layer numerical experiments also reveal that the proposed dynamic models automatically parameterize the subgrid-scale stress terms in each active layer.Furthermore,the proposed scale-aware models dynamically provide higher values of the eddy viscosity for smaller resolutions taking into account the local resolved flow information,and addressing the intimate relationship between the eddy viscosity coefficients and the numerical resolution employed by the quasigeostrophic models.展开更多
基金This study was supported by the Grant-in-Aid for Science Research of the Ministry of Education and Culture, Japan, under the Grant No. 08455232.
文摘The motion of particle clouds formed by dumping dredged material into quiescent waters is experimentally and numerically studied. In the numerical model, the particle phase is modeled by the dispersion model, and turbulence is calculated by the large eddy simulation. The governing equations, including the filtered Navier-Stokes equations and mass transport equation, are solved based on the operator-splitting algorithm and an implicit cubic spline interpolation scheme. The eddy viscosity is evaluated by the modified Smagorinsky model including the buoyancy term. Comparisons of main flow characteristics, including shape, size, average density excess, moving speed and the amount of particles deposited on the bed, between experimental and computational results show that the numerical model well predicts the motion of the cloud from the falling to spreading stage. The effects of silt-fence on the motion of the particle cloud are also investigated.
文摘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.
文摘The most common method to determine the coefficient of Smagorinsky model now is to employ the Germano identity,however it is too complex and expensive in numerical calculation. In this letter we propose a new dynamic formula for determining the coefficient,which is based on the Kolmogorov equation of filtered velocity in physical space.The simplified formula is quite easy to implement in calculation.It is then verified in both homogeneous isotropic turbulence and wall-bounded turbulence by A Priori and A Posteriori tests.
基金financially supported by the National Natural Science Foundation of China(Grant No.51322601)the National Natural Science Foundation of China-China National Petroleum Corporation Joint Fund of Petrochemical Engineering(U1162122)the Fundamental Research Funds for the Central Universities (Grant No.HIT.BRETIV.201315)
文摘In this work, a discrete particle model （DPM） was applied to investigate the dynamic characteristics in a gas-solid bubbling fluidized bed of binary solid particles. The solid phase was simulated by the hard- sphere discrete particle model. The large eddy simulation （LES） method was used to simulate the gas phase. To improve the accuracy of the simulation, an improved sub-grid scale （SGS） model in the LES method was also applied. The mutative Smagorinsky constant case was compared with the previously published experimental data. The simulation by the mutative Smagorinsky constant model exhibited better agreement with the experimental data than that by the common invariant Smagorinsky constant model. Various restitution coefficients and different compositions of binary solids were investigated to determine their influences on the rotation characteristics and granular temperatures of the particles. The particle translational and rotational characteristic distributions were related to certain simulation parameters.
文摘In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds number fluid flow applications, a subgrid turbulence model for LBM was introduced based on standard Smagorinsky subgrid model and Lattice Bhatnagar-Gross-Krook (LBGK) model. The subgrid LBGK model was subsequently used to simulate the two-dimensional driven cavity flow at high Reynolds numbers. The simulation results including distribution of stream lines, dimensionless velocities distribution, values of stream function, as well as location of vertex center, were compared with benchmark solutions, with satisfactory agreements.
文摘A Large Eddy Simulation (LES) technique was applied to solve the turbulentchannel flow for Re_τ = 150 . Three types of turbulence models are employed, such as theSmagorinsky model, the Dynamic Sub-Grid Scale(SGS) model and the Generalized Normal Stress (GNS)model. The simulated data in time series for the LES were averaged in both time and space to carryout the statistical analysis. Results of LES were compared with that of a DNS. As an application, aLES technique was used for 2D body in order to check the validation by investigating the turbulentvortical motion around the afterbody with a slant angle.
文摘A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.
文摘Parallel implementation of algorithm of numerical solution of Navier-Stokes equations for large eddy simulation (LES) of turbulence is presented in this research. The Dynamic Smagorinsky model is applied for sub-grid simulation of turbulence. The numerical algorithm was worked out using a scheme of splitting on physical parameters. At the first stage it is supposed that carrying over movement amount takes place only due to convection and diffusion. Intermediate field of velocity is determined by method of fractional steps by using Thomas algorithm (tridiaginal matrix algorithm). At the second stage found intermediate field of velocity is used for determination of the field of pressure. Three dimensional Poisson equation for the field of pressure is solved using upper relaxation method. Moreover various ways of geometrical decomposition for parallel numerical solution of three dimensional Poisson equations are investigated.
基金Project supported by the National Natural Science Foundation of China(Grant No.51579164)the National Key R&D Program of China(Grant No.2016YFC0502207)
文摘When the sediment and the dissolved matter laden in the river meet a clear water in reservoirs, the turbid water will plunge and spread into the clear water, forming the turbidity current and influencing the water quality and the life of the reservoir. Due to the unsteady nature of the flood, the turbidity current is unsteady. In the present study, we use the MIKE 3 computational fluid dynamics code to simulate continuous and discontinuous turbidity currents on a flat slope. With the model used by us, the turbulence is divided into two parts: the horizontal turbulence and the vertical turbulence, which are separately modeled by the Smagorinsky model and our model to capture the anisotropic turbulence. In this model, the sediment settling and deposition are considered. The simulation results concerning the flume water surface level, the front velocity and sediment concentration profiles are found consistent with the experimental data, particularly, for the sediment concentration profiles with an absolute mean error of 0.026 kg/m3and the root mean square error of 0.046 kg/m3. This finding suggests that this model can be used to well predict the turbidity current on the flat slope.
文摘This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow using two different nonlinear eddy viscosity functional forms to capture anisotropic dissipation mechanism,namely,(i)the Smagorinsky model using the local strain rate field,and(ii)the Leith model using the gradient of the vorticity field.The proposed models are applied to the one-layer and two-layer wind-driven quasigeostrophic ocean circulation problems,which are standard prototypes of more realistic ocean dynamics.Results show that both models capture the quasi-stationary ocean dynamics and provide the physical level of eddy viscosity distribution without using any a priori estimation.However,it is found that slightly less dissipative results can be obtained by using the dynamic Leith model.Two-layer numerical experiments also reveal that the proposed dynamic models automatically parameterize the subgrid-scale stress terms in each active layer.Furthermore,the proposed scale-aware models dynamically provide higher values of the eddy viscosity for smaller resolutions taking into account the local resolved flow information,and addressing the intimate relationship between the eddy viscosity coefficients and the numerical resolution employed by the quasigeostrophic models.
