Evaluation of Subgrid-scale Models in Large-eddy Simulations of Turbulent Flow in a Centrifugal Pump Impeller

The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.

Applying resolved-scale linearly forced isotropic turbulence in rational subgrid-scale modeling

In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is stationary. Forced isotropic turbulence is often used as a framework for establishing and validating such SGS models based on stationary restrictions, for it generates statistical stationary samples. However, traditional forcing method at low wavenumbers cannot provide an analytic form of forcing term for a complete KEF in physical space, which has been illustrated to be essential in the modeling of such SGS models. Thus, an alternative forcing method giving an analytic forcing term in physical space is needed for rational SGS modeling. Giving an analytic linear driving term in physical space, linearly forced isotropic turbulence should be considered an ideal theoretical framework for rational SGS modeling. In this paper, we demonstrate the feasibility of establishing a rational SGS model with stationary restriction based on linearly forced isotropic turbulence. The performance of this rational SGS model is validated. We, therefore, propose the use of linearly forced isotropic turbulence as a complement to free-decaying isotropic turbulence and low-wavenumber forced isotropic turbulence for SGS model validations.

Mathematical Constraints in Multiscale Subgrid-Scale Modeling of Nonlinear Systems

To shed light on the subgrid-seale （SGS） modeling methodology of nonlinear systems such as the Navier-Stokes turbulence, we define the concepts of assumption and restriction in the modeling procedure, which are shown by generalized derivation of three general mathematical constraints for different combinations of restrictions. These constraints are verified numerically in a one-dimensional nonlinear advection equation. This study is expected to inspire future research on the SGS modeling methodology of nonlinear systems.

A SECOND-ORDER DYNAMIC SUBGRID-SCALE STRESS MODEL

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.

Stochastic modeling of subgrid-scale effects on particle motion in forced isotropic turbulence

The subgrid-scale effects on particle motion were investigated in forced isotropic turbulence by DNS and priorLES methods.In the DNS field,the importance of Kolmogorov scaling to preferential accumulation was validated by comparing the radial distribution functions under various particle Stokes numbers.The prior-LES fields were generated by filtering the DNS data.The subgrid-scale Stokes number(StSGS)is a useful tool for determining the effects of subgrid-scale eddies on particle motion.The subgrid-scale eddies tend to accumulate particles with StSGSb 1 and disperse particles with 1 b StSGSb 10.For particles with StSGS?1,the effects of subgrid-scale eddies on particle motion can be neglected.In order to restore the subgrid-scale effects,the Langevin-type stochastic model with optimized parameters was adopted in this study.This model is effective for the particles with StSGS N 1 while has an adverse impact on the particles with StSGSb 1.The results show that the Langevin-type stochastic model tends to smooth the particle distribution in the isotropic turbulence.

A new dynamic subgrid-scale model using artificial neural network for compressible flow

The subgrid-scale(SGS)kinetic energy has been used to predict the SGS stress in compressible flow and it was resolved through the SGS kinetic energy transport equation in past studies.In this paper,a new SGS eddy-viscosity model is proposed using artificial neural network to obtain the SGS kinetic energy precisely,instead of using the SGS kinetic energy equation.Using the infinite series expansion and reserving the first term of the expanded term,we obtain an approximated SGS kinetic energy,which has a high correlation with the real SGS kinetic energy.Then,the coefficient of the modelled SGS kinetic energy is resolved by the artificial neural network and the modelled SGS kinetic energy is more accurate through this method compared to the SGS kinetic energy obtained from the SGS kinetic energy equation.The coefficients of the SGS stress and SGS heat flux terms are determined by the dynamic procedure.The new model is tested in the compressible turbulent channel flow.From the a posterior tests,we know that the new model can precisely predict the mean velocity,the Reynolds stress,the mean temperature and turbulence intensities,etc.

Rational subgrid-scale modelling: a short survey

We review the previous attempts of rational subgrid-scale （SGS） modelling by employing theKolmogorov equation of filtered quantities. Aiming at explaining and solving the underlyingproblems in these models, we also introduce the recent methodological investigations for therational SGS modelling technique by defining the terms of assumption and restriction. Thesemethodological works are expected to provide instructive criterions for not only the rational SGSmodelling, but also other types of SGS modelling practices.

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.

基于Smagorinsky和Vreman模型的竖通道内旋转热流场大涡模拟

为了确定Vreman亚格子模型对于有侧开缝的竖通道内旋转火焰热流场模拟的适用性,基于一竖通道内旋转火焰实验场景,分别使用基于Smagorinsky和Vreman亚格子模型的大涡模拟技术进行了模拟。通过对所得结果比较确定基于Smagorinsky亚格子模型的大涡模拟技术可以得到与实验吻合的火焰状态,而使用Vreman亚格子模型会在一定程度上低估旋转火焰周围的粘性耗散,从而过高地估计旋转火焰周围的切向速度,使火焰不能保持相对稳定的旋转状态。

基于Vreman亚格子模型的有限开口空间内热驱动流数值模拟

利用基于Vreman亚格子模型的大涡模拟技术对有开口的单室和双室房间内热驱动流进行了数值模拟,利用函数分析法定量分析了模拟结果的准确性,并与Smagorinsky亚格子模型的模拟结果进行了比较。结果表明,Vreman和Smagorinsky亚格子模型的计算结果均能够满足工程的需求,但Vreman亚格子模型在开口附近区域的温度和U速度计算结果在整体上比Smagorinsky亚格子模型更接近实验值;Vreman亚格子模型未像Sma-gorinsky亚格子模型那样过高地估算壁面附近高温区域的粘性耗散;对于单室房间内热烟气层高度的预测,采用Vreman模型得到的计算结果准确性比Smagorinsky亚格子模型提高近50%。

基于全局动态Vreman模型的复杂非定常湍流的大涡模拟

在湍流数值模拟方法中,大涡模拟方法可以提供丰富的大涡旋信息,已逐渐成为复杂湍流问题数值研究的重要方法。而大涡模拟中,最重要的一环是尽量准确地构建能反映流场物理本质特征的亚格子应力模型。基于该思想,将一种新型的大涡模拟亚格子应力模型-Vreman亚格子应力模型用于高雷诺数三维后台阶流动的求解,计算结果与实验结果进行对比分析结果较吻合,验证了该模型的可靠性。这是对该模型用于无任何均匀流动方向的高雷诺数复杂湍流非定常流动的首次检验,计算结果优于基于传统的Smagorinsky涡粘性的动态亚格子模型。

Assessment of subgrid-scale models in wall-modeled large-eddy simulations of turbulent channel flows

Considering the demanding of grid requirements for high-Reynolds-number wall-bounded flow,the wall-modeled large-eddy simulation(WMLES)is an attractive method to deal with near wall turbulence.However,the effect of subgrid-scale(SGS)models for wall-bounded turbulent flow in combination with wall stress models is still unclear.In this paper,turbulent channel flow at Reτ=1000 are numerically simulated by WMLES in conjunction with four different SGS models,i.e.,the wall-adapting local eddy-viscosity model,the dynamic Smagorinsky model,the dynamic SGS kinetic energy model and the dynamic Lagrangian model.The mean velocity profiles are compared with the law of the wall,and the velocity fluctuations are compared with direct numerical simulation data.The energy spectrum of velocity and wall pressure fluctuations are presented and the role of SGS models on predicting turbulent channel flow with WMLES is discussed.

Assessment of the modulated gradient model in decaying isotropic turbulence

A recently introduced nonlinear model undergoes evaluations based on two isotropic turbulent cases:a University of Wiscosion-Madison case at a moderate Reynolds number and a Johns Hopkins University case at a high Reynolds number.The model uses an estimation of the subgrid-scale(SGS) kinetic energy to model the magnitude of the SGS stress tensor,and uses the normalized velocity gradient tensor to model the structure of the SGS stress tensor.Testing is performed for the first case through a comparison between direct numerical simulation(DNS) results and large eddy simulation(LES) results regarding resolved kinetic energy and energy spectrum.In the second case,we examine the resolved kinetic energy,the energy spectrum,as well as other key statistics including the probability density functions of velocities and velocity gradients,the skewness factors,and the flatness factors.Simulations using the model are numerically stable,and results are satisfactorily compared with DNS results and consistent with statistical theories of turbulence.

Subgrid-scale model based on the vorticity gradient tensor for rotating turbulent flows

A new subgrid-scale(SGS)stress model is proposed for rotating turbulent flows,and the new model is based on the traceless symmetric part of the square of the velocity gradient tensor and the symmetric part of the vorticity gradient tensor(or the so-called vorticity strain rate tensor).The new subgrid-scale stress model is taken into account the effect of the vortex motions in turbulence,which is reflected on the anti-symmetric part of the velocity gradient tensor.In addition,the eddy viscosity of the new model reproduces the proper scaling as O(y^3)near the wall.Then,the new SGS model is applied in large-eddy simulation of the spanwise rotating turbulent channel flow.Different simulating cases are selected to test the new model.The results demonstrate that the present model can well predict the mean velocity profiles,the turbulence intensities,and the rotating turbulence structures.In addition,it needs no a second filter,and is convenient to be used in the engineering rotational flows.

An improved dynamic subgrid-scale model and its application to large eddy simulation of rotating channel flows

A new dynamic subgrid-scale (SGS) model, which is proved to satisfy the principle of asymptotic material frame indifference (AMFI) for rotating turbulence, is proposed based on physical and mathematical analysis. Comparison with direct numerical simulation (DNS) results verifies that the new SGS model is effective for large eddy simulation (LES) on rotating turbulent flow. The SGS model is then applied to the LES of the spanwise rotating turbulent channel flow to investigate the rotation effect on turbulence characteristics, budget terms in the transport equations of resolved Reynolds stresses, and flow structures near the wall regions of the rotating channel.

A dynamic subgrid-scale model for the large eddy simulation of stratified flow

A new dynamic subgrid-scale (SGS) model, including subgrid turbulent stress and heat flux models for stratified shear flow is proposed by using Yoshizawa’ s eddy viscosity model as a base model. Based on our calculated results, the dynamic subgrid-scale model developed here is effective for the large eddy simulation (LES) of stratified turbulent channel flows. The new SGS model is then applied to the large eddy simulation of stratified turbulent channel flow under gravity to investigate the coupled shear and buoyancy effects on the near-wall turbulent statistics and the turbulent heat transfer at different Richardson numbers. The critical Richardson number predicted by the present calculation is in good agreement with the value of theoretical analysis.

An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.

混流式水轮机尾水管近壁湍流特性和流场结构研究

采用滑移网格技术以及基于Vreman亚格子模型的全局动态大涡模拟方法,得到了HLA551-LJ-43型水轮机小开度工况下全流道内速度、压力以及涡量的分布,同时捕捉到尾水区域高强度的近壁湍流特性以及独特的分离流动现象。计算结果表明,当高转速、小开度工况时,流道内流体的圆周速度占有绝对优势,流体进入尾水管后'被甩出'呈旋转下泄趋势,在锥管段和弯管段内产生了明显的分流,同时在尾水管中心区域形成一个较大的空腔涡带。该涡带的流动形式与常规工况差异较大,反映出小开度工况偏离最优工况最远,流道中流体的漩涡、脱流、分离间断和回流等各种水力不稳定现象更为剧烈和复杂,具有独特的流动形式,容易导致机组异常振动的发生。

混流式水轮机跨尺度流道内复杂湍流的数值模拟

由于水轮机主流道和密封间隙流道空间尺度相差较大,为同时获得水轮机流道和密封通道内的流动特性以及水封间隙泄漏特性,充分考虑叶道流和间隙流的相互干扰,采用非一致网格界面插值技术以及基于Vreman亚格子模型的全局动态大涡模拟方法,不仅得到主流道内速度、压力以及涡量的分布,同时捕捉到了微通道间隙内的流动特征以及不同间隙宽度对密封渗漏流动特性及顶盖压力分布的影响。计算结果表明,主叶道内的强旋拧涡带,密封间隙内旋涡结构以及尾水涡带都直接影响转轮的动态平衡,流体从间隙进入梳齿空腔后产生很强的旋涡而耗散部分动能,而且不同尺度密封间隙流道内产生的旋涡结构不同,耗散能量差异较大,反映出密封渗漏流量的不同。

顶部侧开缝的倾斜通道热流场及开口中性面研究

为了掌握通道在上壁面不同宽度的开缝条件下,不同倾斜角度的通道内热流场特性及两端开口处中性面高度的变化规律,对一200 cm长的通道内一靠近低端开口一侧火源产生的热流场进行了实验与基于Vreman亚格子模型的大涡模拟研究.通过实验值与数值模拟结果的对比,验证了构建的数学模型对于上壁面有侧开缝的倾斜通道内热流场模拟的有效性;确定了通道上壁面的开缝对于通道内热流场的温度分层特性、通道两端开口处U速度以及中性面高度随着通道倾角变化的规律影响不大;对于一定的通道倾斜角度,随着通道上壁面侧开缝宽度增大,火源两端温度分布的对称性增强,通道上壁面附近的温度降低,通道内的高温区域向下壁面的偏移幅度减弱,通道两端开口处靠近上壁面附近的空气流速逐渐降低,两端开口处中性面高度升高;上壁面有开缝的通道两端开口处中性面的高度曲线随着通道倾斜角度的增大有一个交叉,交叉发生在通道倾斜角度小于5°时.