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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
文摘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.
基金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.
文摘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.
文摘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.
文摘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.
基金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.
