离心泵启动过程内部瞬态流动的二维数值模拟

为分析离心泵启动过程外部瞬态特性的内流机理,建立求解叶轮启动和加速过程泵内部非定常流动的数值模拟方法,以二维离心泵为模型进行启动过程内部流动的数值模拟研究.以绝对坐标系描述伞场流动,采用动网格方法实现离心泵启动过程中叶轮加速旋转引起的流场变形,选用人涡模拟(LES)描述湍流.分析了叶轮启动过程巾内部非定常流动结构和演化过程及其与外部瞬态特性的关系,并与准稳态假设下的流动模拟结果进行了对比分析,定性地分析了离心泵在启动过程中的瞬态效应,验证了采用动网格方法求解启动过程瞬态流动的可行性.     

非定常叶顶间隙泄漏流动和换热的数值研究

通过数值方法研究了带叶顶间隙的某一级半透平中的非定常流动和换热问题.数值模拟采用标准k-ω两方程湍流模型,求解非定常雷诺平均N-S方程.动叶顶部间隙取为0.4 mm.分析了动静干涉对动叶顶部间隙内泄漏流动与换热的影响.结果表明,周期性通过的上游静叶尾迹和通道涡足动叶通道中非定常现象的主要来源.流场的波动主要存在于叶顶吸力面侧中间弦长附近.叶顶换热系数波动主要存在于两个位置,一是叶顶吸力面侧,一是叶顶主泄漏通道.叶顶表面面积平均传热系数非定常计算的时均结果与定常计算获得的结果偏差小于2%.     

低压高负荷涡轮叶栅的大涡模拟

了解预测转捩与层流分离泡,对于低压高负荷叶片设计特别重要.为此,采用大涡模拟方法,对典型高负荷叶栅T106进行了模拟.计算描述了层流分离泡,给出了壁面压力、壁面摩擦力系数的分布,从而准确预测了T106压力面上的层流分离、转捩、湍流再附现象.计算同时刻画了层流分离泡的非定常性质,给出了瞬时压力,进行了频谱分析,确定主频为2500 Hz,并观察到了瞬时涡等有意义的现象.计算同时表明二维大涡模拟计算结果可以复现出湍流的许多重要现象并且精度可以满足工程需要.     

基于时间的驱动方腔流的高精度多重网格方法数值模拟

提出了一种求解二维非定常不可压缩Navier-Stokes方程组的全隐二阶时间推进和空间四阶差分紧致格式,在每一个时间步上,采用多重网格的全近似格式(FAS)加速其迭代收敛过程。利用该方法对驱动方腔流动问题进行了直接数值模拟,结果显示对于Re≤5000,本文在粗网格上(64×64等分和128×128等分)即可得到较为准确的定常层流解;对于Re=7500和10000,由于更多二次涡的出现,本文在256×256等分网格上同样可得到与前人的结果相吻合的非定常周期性解。     

基于Omega涡识别方法的离心泵内非定常流动数值评价

离心泵内部非定常流动与流场中生成的旋涡密切相关.本文采用Omega涡识别方法研究离心泵流场涡结构及非定常流动特性,阈值选择为0.51.分析得出LES湍流模型相对比其他两种湍流模型涡识别结果更为细致,对叶片流道中的小涡结构可以很好地捕捉.设计工况下,LES模型识别出的涡面积数值是DDES模型的1.15倍,是SSTk-ω模型的1.55倍.压力脉动受流场空间位置与工况影响较大,其主要激励为叶频,并且压力脉动强度与涡面积密切相关.叶轮出口处的涡结构周期性脱落,在隔舌处会撞击隔舌而破裂.离心泵动静干涉的产生与叶轮出口处的涡结构脱落有关.     

湍流热对流近底板流态与温度边界层特性

本文采用湍流热对流的并行直接数值模拟(PDM-DNS),计算了系列Ra数的二维方腔和三维扁方腔的Rayleigh-Bénard热对流.针对平均场计算结果,选取Ra=10~9,10~(10),5×10~(10),讨论了二维方腔和展向平均三维扁方腔热对流流动特性.发现二维方腔和三维扁方腔流动中都存在大尺度环流和角涡,而且随着Ra数的增加,大尺度环流的形状变圆,角涡尺寸变小.在二维方腔流动中,大尺度环流呈椭圆形,有四个角涡,而展向平均三维扁方腔流动中,大尺度环流呈梭形,只有两个角涡.由于角涡特性的不同,二维方腔流动中羽流向上运动的范围比展向平均三维扁方腔流动更广,造成二维流动局部区域温度分布高温层厚度变大.温度边界层厚度λ_θ与Ra数之间存在标度关系,二维方腔和三维扁方腔热对流温度边界层变化的标度指数基本一致,标度关系的系数稍有不同.     

基于大涡模拟下风力机翼型非定常转捩流动的动态模态分解分析

翼型绕翼流动对风力机整机性能产生重要影响.本文基于大涡模拟方法,得到风力机翼型非定常转捩流动.通过对压力流场的动态模态分解(DMD)分析,发现翼型层流分离泡的生成和发展为流动主要非定常特征,且该特征具有主要频率.预估得到高增长率T-S扰动波频率与DMD模态频率接近,发现高频变化的分离泡由边界层分离点下游不远处的T-S扰动波诱导K-H不稳定性而主导,以及中低频的DMD模态表现为T-S扰动波所引起的湍流结构.对比不同攻角下的DMD分解结果,发现上表面分离泡会逐渐向前移动,长度变短;而下表面分离泡略微后移,长度变化不大.     

带标量修正的气膜冷却各向异性湍流模型

本文采用作者提出的带标量修正的代数各向异性湍流模型,对平板带复合角气膜冷却在带横槽/不带横槽情况下的流动传热进行了研究.通过比较不同工况下的流动传热预测结果,讨论了湍流动量、标量输运对于气膜冷却预测的影响.各向异性修正的雷诺应力项能减弱各向同性湍流模型所带来的肾形涡系强度过预测的问题。而各向异性的湍流标量扩散项则能有效地将雷诺应力和温度梯度对于湍流标量扩散系数的影响体现出来。结果表明带标量修正的各向异性湍流模型能很好地提高气膜冷却的预测精度.     

基于DDES模拟的叶顶间隙流湍流特性研究

叶顶间隙泄漏流使得叶片顶部流动变得复杂,对间隙泄漏流的流动结构进行精细地捕捉并探讨其湍流特性有利于更深入的了解间隙泄漏流的流动机理,为控制泄漏损失提供依据。本文利用延迟脱体涡模拟方法对叶顶间隙泄漏流动进行非定常数值模拟;然后从叶顶间隙流的流动涡结构演变,湍动能和各向异性等方面研究叶顶间隙流的湍流特性;利用POD模态分解的方法对湍流特性做进一步分析,对分流场的基本规律、演变特点进行了详细研究;最后结合熵生成率,对泄漏流流动结构造成的损失进行了分析。研究发现,泄漏流在下游区存在着强湍流特性,涡尺度受不同的雷诺应力影响,且大尺度流动结构造成主要的耗散损失。     

凹腔非定常特性的数值模拟

采用基于Menter SST两方程湍流模型的DES方法,数值模拟开式凹腔在跨声速条件下的非定常流动特性.计算凹腔底部和后壁面上的点的声压级频谱以及总声压级,证明在第二噪声模态上的声压级最大.     

Computational study of the laminar to turbulent transition over the SD7003 airfoil in ground effect

In this work, we present a numerical study of the laminar-turbulence transition flow around a symmetrical air-foil at a low Reynolds number in free flow and near the ground surface at different angles of attack. Finite volume method is employed to solve the unsteady Reynolds-averaged Navier–Stokes (RANS) equation. In this way, the Transition SST turbulence model is used for modeling the flow turbulence. Flow around the symmetrical airfoil SD7003 is numerically simulated in free stream and near the ground surface. Our numerical method can detect different aspects of flow such as adverse pressure gradient, laminar separation bubble and laminar to turbulent transition onset and the numerical results are in good agreement with the experimental data.     

Combined effects of viscous dissipation and MHD on free convection flow past a semi-infinite vertical plate with variable surface temperature in the presence of heat source

An unsteady MHD laminar viscous dissipative fluid flow past a semi-infinite vertical plate with variable surface temperature in the presence of heat source is considered in the present analysis. The present approach transforms the governing boundary layer equations into nondimensional form using the appropriate nondimensional quantities, which is valid in the free convection region. The resulting governing equations are solved numerically using the Crank–Nicolson method, an efficient implicit finite-difference scheme. Numerical results are obtained and presented in the form of local as well as average shearing stress, local and average heat transfer rate, velocity and temperature during the transient period. The present results are compared with the available results in the literature and are found to be in an excellent agreement.     

On flow unsteadiness induced by structural vibration

In this paper, the effects of structural vibration on flow unsteadiness are investigated numerically. A fully coupled model, that solves the unsteady flow equations as well as the dynamic equations of the structure, is used. Numerical experiments are carried-out for flow over a backward-facing step, where a large number of numerical and experimental data exist for comparisons. The flexible structure is upstream of the step and is excited by a plane acoustic wave from the side opposite to the flow. Three Reynolds number cases are studied: 300 for a laminar flow, 3000 for a transitional flow and 15 000 for a turbulent flow. The results obtained are in good agreement with experimental observations and show the strong coupling between structural vibration and the resulting flow unsteadiness.     

Qualitative comparison between numerical and experimental results of unsteady flow in a radial diffuser pump

Comparison between numerical simulation and experimental results for unsteady flow field in a radial diffuser pump is presented for the design operating point. The numerical result is obtained by solving three-dimensional, unsteady Reynolds-averaged Navier-Stokes equations by the commercial CFD code CFX-10 withk-ω based shear stress transport turbulence model. Two-dimensional PIV measurements are conducted to acquire the experiment result. The phase-averaged velocity and turbulent kinetic energy fields are compared in detail between the results by the two methods in the impeller, diffuser and return channel regions. The qualitative comparison between CFD and PIV results is quite good in the phase-averaged velocity field. Although the turbulence level by PIV is higher than that by CFD generally, the main turbulence features are nearly the same. Furthermore, the blade orientation effect and other associated unsteady phenomena are also examined, in order to enhance the understanding on impeller-diffuser interaction in a radial diffuser pump.     

Effect of couple stresses on transient MHD poiseuille flow

The unsteady laminar flow of an electrically conducting viscous fluid between parallel insulating plates subject to a transverse magnetic field is considered. The plates are fixed and flow is due to a constant pressure gradient. The induced field is taken into account. The fluid is incompressible and of couple stress type. The defining equations are coupled and numerical solutions for different values of couple stress parameter are obtained. The velocity and induced magnetic field profiles are sketched as functions of time, Hartmann number, and magnetic Prandtl number. The velocity decreases with increase in couple stress parameter.     

分数阶Oldroyd-B黏弹性Poiseuille流的Laplace数值反演分析

采用Laplace数值反演的Stehfest算法研究了分数阶Oldroyd-B粘弹性流体在两平板间非定常的Poiseuille流动问题.首先,通过数值解与近似解析解的比较验证了Stehfest算法的有效性.其次,运用Stehfest算法对平板Poiseuille流动进行了研究,揭示了分数阶黏弹性平板流的速度过冲和应力过冲现象,指出这些现象对分数导数的阶数存在明显的依赖性.同时,数值结果表明,整数阶本构方程仅仅是分数阶本构方程的特例,分数阶本构方程较整数阶本构方程具有更广泛的适用性。     

An implicit high-order spectral difference approach for large eddy simulation

The filtered fluid dynamic equations are discretized in space by a high-order spectral difference (SD) method coupled with large eddy simulation (LES) approach. The subgrid-scale stress tensor is modelled by the wall-adapting local eddy-viscosity model (WALE). We solve the unsteady equations by advancing in time using a second-order backward difference formulae (BDF2) scheme. The nonlinear algebraic system arising from the time discretization is solved with the nonlinear lower–upper symmetric Gauss–Seidel (LU-SGS) algorithm. In order to study the sensitivity of the method, first, the implicit solver is used to compute the two-dimensional (2D) laminar flow around a NACA0012 airfoil at Re = 5 × 10⁵ with zero angle of attack. Afterwards, the accuracy and the reliability of the solver are tested by solving the 2D “turbulent” flow around a square cylinder at Re = 10⁴ and Re =  2.2 × 10⁴. The results show a good agreement with the experimental data and the reference solutions.     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

Numerical simulation of transient processes in hydroturbines

A method for calculation of unsteady 3D turbulent flows in hydro turbines of power plants developed for simulation of the transient processes is presented herein. The method is based on joint solution to the Reynolds-averaged Navier—Stokes equations for incompressible fluid flow written for moving mesh, the equation of runner rotation and the system of 1D equations describing propagation of elastic hydraulic shock in the flow domain. The exchange in flow parameters between the penstock and hydro turbine regions is considered in this approach. Results of transient processes simulation are presented for several modes: start-up to the turbine regime, instantaneous load shedding, and output power decrease. Comparison with experimental data is performed.     

Numerical Study of Unsteady MHD Flow and Entropy Generation in a Rotating Permeable Channel with Slip and Hall Effects

This article investigates an unbiased analysis for the unsteady two-dimensional laminar flow of an incompressible, electrically and thermally conducting fluid across the space separated by two infinite rotating permeable walls.The influence of entropy generation, Hall and slip effects are considered within the flow analysis. The problem is modeled based on valid physical arguments and the unsteady system of dimensionless PDEs (partial differential equations) are solved with the help of Finite Difference Scheme. In the presence of pertinent parameters, the precise movement of the flow in terms of velocity, temperature, entropy generation rate, and Bejan numbers are presented graphically, which are parabolic in nature. Streamline profiles are also presented, which exemplify the accurate movement of the flow. The current study is one of the infrequent contributions to the existing literature as previous studies have not attempted to solve the system of high order non-linear PDEs for the unsteady flow with entropy generation and Hall effects in a permeable rotating channel. It is expected that the current analysis would provide a platform for solving the system of nonlinear PDEs of the other unexplored models that are associated to the two-dimensional unsteady flow in a rotating channel.