挤出胀大流动的有限元方法研究

本文对Luo-Tanner提出的流线有限元作了重要的改进,提出了沿通过单元高斯点的流线积分本构方程的方法,迴避了速度梯度在单元边界上间断和出口处应力奇点的困难,同时减少了计算量,对比计算表明,采用压力不连续单元来加强不可压缩性限制能使计算质量和收敛性都得到显著的提高,对Maxwell流体的轴对称挤出胀大流动在Weissenberg数1.2下获得了合理的收敛解。     

前向多翼叶轮流动分析及出口滑移的计算

本文分析了离心风机前向多翼叶轮内的流动特性,计算了跨叶片面内的理想流场和湍流边界层.用三次样条平滑和旋转坐标变换方法解决了用流线曲率法数值解大曲率流场时的收敛问题;用考虑旋转和曲率对流动结构影响的Richardson数修正理想流场.提出了一种予测叶片边界层分离及流动滑移的方法.     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

各向加权异性张力有限元法及其应用

本文利用各向异性张力有限元的特点设计了灵活而有效的各向加权异性有限元算法，采用这种方法，计算了离心叶轮内部准三元流场。     

翼型跨音速流动计算和外形设计的变域变分有限元方法

在气体动力学问题研究中经常会碰到诸如激波、翼型设计等未知界面问题。未知界面的存在为该类问题的理论分析和数值求解带来了很大困难。刘高联针对未知界面问题发展了一种变域变分有限元方法，该方法将未知界面看作是一个变化区域的边界，采用变域变分将未知界面结合在变分泛函中，使其与求解流场的控制方程结合起来，从而将未知界面的求解和流场的求解完全耦合进行，因而是一种处理未知界面的独特工具，极适合于气动外形的设计求解。本文运用变域变分有限元方法对翼型跨音速流动正、反命题进行了数值研究。由于在跨音速翼型绕流中存在激波，所以为了得到压缩激波解，采用了“人工密度”办法。几个算例均得到了满意的计算结果和设计结果，证明了本文方法的有效性和优越性。     

非正交网格控制容积法分析复杂形状池内的流动

介绍一种基于非正交网格控制容积法的数学模型，及其在圆形沉沙池流动研究中的应用．该模型求解轴对称流动的连续方程及时均Ｎ－Ｓ方程，并采用标准ｋ－ε紊流模型，模拟圆形池内的流动．由于采用非正交网格，此计算模型可精确模拟几何形状较复杂的沉沙池内的流动，利用上述模型对某实际沉沙池进行了流场计算，计算所得流场与模型试验实测值符合良好．     

非正交网格控制容积法分析复杂形状池内的流动

介绍一种基于非正交网格控制容积法的数学模型，及其在圆形沉沙池流动研究中的应用．该模型求解轴对称流动的连续方程及时均Ｎ－Ｓ方程，并采用标准ｋ－ε紊流模型，模拟圆形池内的流动．由于采用非正交网格，此计算模型可精确模拟几何形状较复杂的沉沙池内的流动，利用上述模型对某实际沉沙池进行了流场计算，计算所得流场与模型试验实测值符合良好．     

应用变域变分有限元法解平面叶栅反命题

反命题作为一种可变(未知)边界问题近年来得到了广泛的研究。本文给出了亚声速平面叶栅反命题计算的势函数变域变分有限元解法。变域变分通过把可变边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而使可变边界求解和流场分析可以完全耦合进行。本文针对亚声速平面叶栅的反命题，根据泛函的驻值必要条件，介绍了变域变分有限元方法的求解过程，最后给出了两个数值算例。这两个算例均采用四节点矩形单元的插值基函数，第一个算例用于检验程序的可靠性，第二个算例设计了一个给定叶面马赫数分布的叶型，并与试验结果进行对照。     

离心风机子午通道内湍流场数值模拟

由进风口－叶轮－无叶扩压器－蜗壳等部件组成的离心风机通道内流分析是非常复杂的，目前还只能是分别计算各部件内的流场，但必须考虑部件间的相互影响。本文采用轴对称Ｎ－Ｓ方程，根据三维叶轮通道计算给出的叶片力分布，求解了考虑叶片力的进风口－叶轮－无叶扩压器组成的子午通道问题，所得结果可用来给出三维叶轮通道计算的进口条件，并可用于优化设计进风口及叶轮前、后盘形状。该方法已得到实践检验。     

应用有限元法计算离心式叶轮内部流场

运用张量分析,采用半测地坐标系,导入了具有不变形式的流函数,推志了离心式叶轮内部准三元流动任意流面上函数所应满足的偏微分方程和定解条件,并运用伽辽金方法建立了相应的有限元方程。对离心式叶轮回转面上流动进行了计算,数值计算结果正确地揭示了离心式叶轮内部的流动规律。     

离心压缩机叶轮三维有限元强度分析系统

根据离心压缩机叶轮形状的特点,通过有效的网格自动划分和边界条件的处理,采用三维有限元方法进行了强度计算与分析,从而建立了离心叶轮三维有限元强度分析系统。对于叶轮与主轴之间的过盈配合,求出配合处由于变形而产生的弹性支力,以载荷形式施加于配合之处,由于计算时不加入主轴,从而减少了一,通过等厚度圆尖力的计算验证了有限元程序的可靠性,最二,对一个实际的三元流叶轮进行了计算,通过编制的后处理软件,对影响叶轮     

Flow analysis for a double suction centrifugal machine in the pump and turbine operation modes

A double suction centrifugal machine has been studied, both experimentally and numerically, operating as a pump and as a turbine. Experimentally, the static performance of the machine working as a pump was obtained. These measurements were compared with equivalent numerical results from a URANS calculation. As a second step, the numerical results have been exploited to get detailed information about the flow in both operating modes (pump and turbine). The main goals of the study are, first, the validation of the numerical procedure proposed and second, the possible turbine operation of the impeller, which could point out a wider working range for the machine. The first aspect is handled by detailed analysis in the pump mode, according to previous experience of the research group. The second objective is obtained by using the numerical model to explore the flow fields obtained, when working in an inverse mode. Therefore, the presented results join the use of a numerical methodology and the turbine mode of operation for a centrifugal impeller, providing insight into the flow characteristics. When working as a pump, the flow at the suction side is characterized by the existence of an inlet tongue, which tends to enforce a uniform flow for the nominal conditions. For the turbine mode, flow patterns in the impeller, volute and suction regions are carefully investigated. The influence of the specific geometrical arrangement is also considered for this operation mode. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite Element Computations of Flow around a Wall-mounted Cube

In this paper, we present the application of a finite element scheme to full three-dimensional incompressible flow around a cube mounted on the wall in a channel. This scheme is based on the Petrov-Galerkin weak formulation using exponential weighting functions. The incompressible Navier-Stokes equations are numerically integrated in time by using a fractional step strategy with a second-order accurate Adams-Bashforth scheme. The workability and validity of the present approach are demonstrated through the results of streamlines and pressure coefficients in the flow field up to high Reynolds number regimes.     

Streamline computations for porous media flow including gravity

In this paper we consider porous media flow without capillary effects. We present a streamline method which includes gravity effects by operator splitting. The flow equations are treated by an IMPES method, where the pressure equation is solved by a (standard) finite element method. The saturation equation is solved by utilizing a front tracking method along streamlines of the pressure field. The effects of gravity are accounted for in a separate correction step. This is the first time streamlines are combined with gravity for three-dimensional (3D) simulations, and the method proves favourable compared to standard splitting methods based on fractional steps. By our splitting we can take advantage of very accurate and efficient 1D methods. The ideas have been implemented and tested in a full field simulator. In that context, both accuracy and CPU efficiency have tested favourably.     

Unsteady flow calculation in a radial flow centrifugal pump with spiral casing

The prediction of the two-dimensional unsteady flow established in a radial flow centrifugal pump is considered. Assuming the fluid incompressible and inviscid, the velocity field is represented by means of source and vorticity surface distributions as well as a set of point vortices. Using this representation, a grid-free (Lagrangian) numerical method is derived based on the coupling of the boundary element and vortex particle methods. In this context the source and vorticity surface distributions are determined through the non-entry boundary condition together with the unsteady Kutta condition. In order to satisfy Kelvin's theorem, vorticity is shed at the trailing edges of the impeller blades. Then the vortex particle method is used to approximate the convection of the free vorticity distribution. Results are given for a pump configuration experimentally tested by Centre Technique des Industries Mécaniques (CETIM). Comparisons between predictions and experimental data show the capability of the proposed method to reproduce the main features of the flow considered.     

PIV measurement of internal flow characteristics of very low specific speed semi-open impeller

Detailed particle-image velocimetry (PIV) measurements of flow fields inside semi-open impellers have been performed to understand better the internal flow patterns that are responsible for the unique performance of these centrifugal pumps operated in the range of very low specific speed. Two impellers, one equipped with six radial blades (impeller A) and the other with four conventional backward-swept blades (impeller B), are tested in a centrifugal pump designed to be operated at a non-dimensional specific speed of n_s=0.24. Complex flow patterns captured by PIV are discussed in conjunction with the overall pump performance measured separately. It is revealed that impeller A achieves higher effective head than impeller B even though the flow patterns in impeller A are more complex, exhibiting secondary flows and reverse flows in the impeller passage. It is shown that both the localized strong outward flow at the pressure side of each blade outlet and the strong outward through-flow along the suction side of each blade are responsible for the better head performance of impeller A.     

A novel experimental facility for measuring internal flow of Solid-liquid two-phase flow in a centrifugal pump by PIV

PIV is one of the most effective flow measurement methods, but for solid-liquid two-phase flow in a centrifugal pump with complex internal structure and high speed rotation, many problems exist in the experimental facility design, such as solid particle releasement, synchronization issues among camera, laser and pump impeller, and methods of avoiding introducing bubbles into flow passage, etc. To solve these critical problems, a new experimental facility without agitation was designed for internal solid-liquid two-phase flow test in a centrifugal pump by PIV. The facility was tested by measuring two-phase flow of water-glass beans and water-rape seeds in a centrifugal slurry pump. The new test facility has many obvious advantages, such as the simpler structure without the agitating device, less noise in the test process, less power consumption, less number of particles to get a certain particle volume concentration, and the decrease of the experiment cost. The result shows that the new facility can realize a stable solid-liquid two-phase flow measurement with high reliability in a centrifugal pump.     

离心泵叶轮内密相液固两相湍流的数值模拟

利用作者建立的描述密相液固两相湍流的 KET模型和推导的基本控制方程组 ,在处理壁面边界条件时考虑了颗粒和叶片的相互碰撞作用 ,对离心泵叶轮内密相液固两相流动进行了数值模拟 ,得到了泵叶轮内两相流动的一些规律 ,为液固两相流泵的设计提供了一定的理论依据。     

Measurement of rotating flows using PIV and image derotation

This paper describes the use of a rotating all-mirror image derotator system, high-speed video and particle image velocimetry (PIV) to visualise and quantitatively examine the flow patterns between the blades of a centrifugal impeller. The flow field relative to the moving centrifugal impeller is presented. Published online: 13 December 2002