圆形贫化电炉熔池内电场和温度场的数值计算

从麦克斯韦方程组出发推导了圆形三电极贫化电炉熔池内电传导方程，提出了炉渣湍流导热系数概念和估算方法．运用这一方法可大大简化温度场的计算过程．作者还研制出圆柱坐标系下电炉熔池内电场、温度场解析程序，可迅速解出炉内电场、温度场的分布     

离心泵叶轮内部清水湍流的动态大涡模拟

提出旋转坐标下的单相液体湍流的二阶双系数动态亚格子应力大涡模拟模型并加以验证。考虑亚格子应力的对称性和量纲一致性 ,在基于应变率张量和旋转率张量的不可约量及 Smagorinsky模型和亥姆霍兹定理的基础上 ,提出亚格子应力应表示为可解的应变率张量和旋转率张量的函数。在贴体坐标系下 ,利用交错网格系统和有限体积法离散湍流控制方程 ,采用 SIMPL EC方法求解方程。水泵叶轮中湍流的流速分布规律和压力分布等计算结果都与试验结果基本吻合 ,从而证明了此模型的可行性 ,以及计算方法和程序的可靠性     

对流边界层顶部夹卷速度参数化的分析研究

应用 1979年 7月美国CIRCE试验边界层观测资料对两个夹卷速度参数化方案进行了对比验证 ,结果表明 :基于简单夹卷层结构模型的参数化方案we/w =A·Ri- 1 只适用于自由对流边界层 ,当对流边界层中存在较强机械湍流时 ,该方案变得不再适用 ;而基于“侵吞模型”的参数化方案we/w =B·θ /(Γzi)则不仅适用于自由对流边界层 ,对于含有较强机械湍流的对流边界层也能较好地符合 ,因而具有较强的适用性 .当对流边界层受机械湍流和浮力湍流共同驱动时 ,仅以浮力对流速度w 表征混合层垂直湍流速度尺度是不恰当的 ,还应该包含机械湍流的贡献 .对比验证表明 ,这一改进非常有效 ,不仅使物理意义更加合理 ,也使夹卷速度参数化方案与实际大气状况符合得更好     

Application of a line implicit method to fully coupled system of equations for turbulent flow problems

For unstructured finite volume methods, we present a line implicit Runge–Kutta method applied as smoother in an agglomerated multigrid algorithm to significantly improve the reliability and convergence rate to approximate steady-state solutions of the Reynolds-averaged Navier–Stokes equations. To describe turbulence, we consider a one-equation Spalart–Allmaras turbulence model. The line implicit Runge–Kutta method extends a basic explicit Runge–Kutta method by a preconditioner given by an approximate derivative of the residual function. The approximate derivative is only constructed along predetermined lines which resolve anisotropies in the given grid. Therefore, the method is a canonical generalisation of point implicit methods. Numerical examples demonstrate the improvements of the line implicit Runge–Kutta when compared with explicit Runge–Kutta methods accelerated with local time stepping.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods

In this paper, we present a finite element method with a residual‐based artificial viscosity for simulation of turbulent compressible flow, with adaptive mesh refinement based on a posteriori error estimation with sensitivity information from an associated dual problem. The artificial viscosity acts as a numerical stabilization, as shock capturing, and as turbulence capturing for large eddy simulation of turbulent flow. The adaptive method resolves parts of the flow indicated by the a posteriori error estimates but leaves shocks and turbulence under‐resolved in a large eddy simulation. The method is tested for examples in 2D and 3D and is validated against experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling effects of subgrid-scale mixture fraction variance in LES of a piloted diffusion flame

A posteriori analysis of the statistics of two large-eddy simulation (LES) solutions describing a piloted methane–air (Sandia D) flame is performed on a series of grids with progressively increased resolution reaching about 10.5 million cells. Chemical compositions, density and temperature fields are modelled with a steady flamelet approach and parametrised by the mixture fraction. The difference between the LES solutions arises from a different numerical treatment of the subgrid scale (SGS) mixture fraction variance – an important quantity of interest in non-premixed combustion modelling. In the first case (model I), the variance transport equation is solved directly, while in the second (model II), an equation for the square of the mixture fraction is solved, and the variance is computed from its definition. The comparison of the LES solutions is based on the convergence properties of their statistics with respect to the turbulence resolution length scale. The dependence of the LES statistics is analysed for velocity and the mixture fraction fields, and tested for convergence. For the most part, the statistics converge for the finest grids, but the variance of the mixture fraction shows some residual grid dependence in the high-gradient regions of the jet near field. The SGS variance given by model I exhibits realisability everywhere, whereas in regions of the flame model II is non-realisable, predicting negative variances. Furthermore, the LES statistics of model I exhibit superior convergence behaviour.