Analysis of flow in the plate‐spiral of a reaction turbine using a streamline upwind Petrov–Galerkin method

The prediction of the flow field in a novel spiral casing has been accomplished. Hydraulic turbine manufacturers are considering the potential of using a special type of spiral casing because of the easier manufacturing process involved in its fabrication. These special spiral casings are known as plate‐spirals. Numerical simulation of complex three‐dimensional flow through such spiral casings has been accomplished using a finite element method (FEM). An explicit Eulerian velocity correction scheme has been deployed to solve the Reynolds‐average Navier–Stokes equations. The simulation has been performed to describe the flow in high Reynolds number (10⁶) regimes. For spatial discretization, a streamline upwind Petrov–Galerkin (SUPG) technique has been used. The velocity field and the pressure distribution inside the spiral casing reveal meaningful results. Copyright © 2000 John Wiley & Sons, Ltd.     

Large eddy simulation of incompressible free round jet with discontinuous Galerkin method

In the paper, discontinuous Galerkin method is applied to simulation of incompressible free round turbulent jet using large eddy simulation with eddy viscosity approach. The solution algorithm is based on the classical projection method, but instead of the solution of the Poisson equation, a parabolic equation is advanced in pseudo‐time, which provides the pressure field ensuring the proper pressure–velocity coupling. For time and pseudo‐time integration, explicit Runge–Kutta method is employed. The computational meshes consist of hexahedral elements with flat faces. Within a given finite element, all flow variables are expressed with modal expansions of the same order (including velocity and pressure). Discretisation of the viscous terms in the Navier–Stokes equations and Laplacian in the Poisson equation is stabilised with mixed finite element approach. The correctness of the solution algorithm is verified in a commonly used test case of laminar flow in 3D lid‐driven cavity. The results of computations of the free jet are compared with experimental and numerical reference data, the latter obtained from the high‐order pseudospectral code. The statistics of centerline flow velocity – mean velocity and its fluctuations – show satisfactory agreement with the reference data. Copyright © 2015 John Wiley & Sons, Ltd.     

A preconditioned semi‐staggered dilation‐free finite volume method for the incompressible Navier–Stokes equations on all‐hexahedral elements

A new semi‐staggered finite volume method is presented for the solution of the incompressible Navier–Stokes equations on all‐quadrilateral (2D)/hexahedral (3D) meshes. The velocity components are defined at element node points while the pressure term is defined at element centroids. The continuity equation is satisfied exactly within each elements. The checkerboard pressure oscillations are prevented using a special filtering matrix as a preconditioner for the saddle‐point problem resulting from second‐order discretization of the incompressible Navier–Stokes equations. The preconditioned saddle‐point problem is solved using block preconditioners with GMRES solver. In order to achieve higher performance FORTRAN source code is based on highly efficient PETSc and HYPRE libraries. As test cases the 2D/3D lid‐driven cavity flow problem and the 3D flow past array of circular cylinders are solved in order to verify the accuracy of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

A study of penalty elements for incompressible laminar flows

A finite element model is developed based on the penalty formulation to study incompressible laminar flows. The study includes a number of new quadrilateral and triangular elements for 2-dimensional flows and a number of new hexahedral and tetrahedral elements for 3-dimensional flows. All elements employ continuous velocity approximations and discontinuous pressure approximations respecting the LBB condition of numerical instability. An incremental Newton–Raphson method coupled with the Broyden method is used to solve the non-linear equations. Several numerical examples (colliding flow, cavity flow, etc.) are presented to assess the efficiency of elements.     

A stable unstructured finite volume method for parallel large-scale viscoelastic fluid flow calculations

A new stable unstructured finite volume method is presented for parallel large-scale simulation of viscoelastic fluid flows. The numerical method is based on the side-centered finite volume method where the velocity vector components are defined at the mid-point of each cell face, while the pressure term and the extra stress tensor are defined at element centroids. The present arrangement of the primitive variables leads to a stable numerical scheme and it does not require any ad-hoc modifications in order to enhance the pressure–velocity–stress coupling. The log-conformation representation proposed in [R. Fattal, R. Kupferman, Constitutive laws for the matrix–logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech. 123 (2004) 281–285] has been implemented in order improve the limiting Weissenberg numbers in the proposed finite volume method. The time stepping algorithm used decouples the calculation of the polymeric stress by solution of a hyperbolic constitutive equation from the evolution of the velocity and pressure fields by solution of a generalized Stokes problem. The resulting algebraic linear systems are solved using the FGMRES(m) Krylov iterative method with the restricted additive Schwarz preconditioner for the extra stress tensor and the geometric non-nested multilevel preconditioner for the Stokes system. The implementation of the preconditioned iterative solvers is based on the PETSc library for improving the efficiency of the parallel code. The present numerical algorithm is validated for the Kovasznay flow, the flow of an Oldroyd-B fluid past a confined circular cylinder in a channel and the three-dimensional flow of an Oldroyd-B fluid around a rigid sphere falling in a cylindrical tube. Parallel large-scale calculations are presented up to 523,094 quadrilateral elements in two-dimension and 1,190,376 hexahedral elements in three-dimension.     

THE INFLUENCE OF THE NEW TECHNOLOGY OF FIXED SURFACE PERFORATION ON CASING STRENGTH1)

用于非常规储层的定面射孔工艺可有效地降低水力裂缝的起裂压力,提高近井筒的泄油面积;但该工艺对套管强度的影响研究较少。本文基于流固耦合理论,建立了地层$\!$-$\!$-$\!$水泥环$\!$-$\!$-$\!$射孔套管的应力耦合模型,退化为裸套管与API标准进行比对,效验了模型,并发现地层约束条件下套管承受内压的能力提高了1.48$\sim$2.16倍,在此基础上,给出了套管尺寸、射孔直径及射孔相位角对套管承压能力的变化规律。计算表明：螺旋射孔方位角为45$^\circ$时,定面射孔间夹角为30$^\circ$时,套管内径越小,套管壁厚越厚,则套管的承压能力越高;相同参数下,定面射孔工艺比螺旋射孔工艺引起的套管承压能力略低。     

Consistent vs. reduced integration penalty methods for incompressible media using several old and new elements

The frequently used reduced integration method for solving incompressible flow problems ‘a la penalty’ is critically examined vis-a-vis the consistent penalty method. For the limited number of quadrilateral and hexahedral elements studied, it is shown that the former method is only equivalent to the latter in certain special cases. In the general case, the consistent penalty method is shown to be more accurate. Finally, we demonstrate significant advantages of a new element, employing biquadratic (2-D) or triquadratic (3-D) velocity and linear pressure over that using the same velocity but employing bilinear (2-D) or trilinear (3-D) pressure approximation.     

圆筒内横向流体与弹性管耦合的分域求解算法

海洋丛式井组隔水管和换热器管束等在流体作用下,会诱导管束振动及碰撞,导致管束断裂失效。将弹性管离散成梁单元,采用非线性结构动力学方程描述;将圆筒流体域离散成实体单元,采用计算流体动力学方程描述。在流固耦合界面处,推导了界面位移、速度和载荷计算公式及收敛判断准则,建立了圆筒内横向流体与弹性管耦合的分域求解算法。算例表明,其分域与全域求解计算结果吻合较好,本文算法为复杂流体域内多根管束的振动及碰撞问题求解提供行之有效的计算方法。     

A family of multi‐point flux approximation schemes for general element types in two and three dimensions with convergence performance

A family of flux‐continuous, locally conservative, control‐volume‐distributed multi‐point flux approximation (CVD‐MPFA) schemes has been developed for solving the general geometry‐permeability tensor pressure equation on structured and unstructured grids. These schemes are applicable to the full‐tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full‐tensor flow approximation. The family of flux‐continuous schemes is characterized by a quadrature parameterization. Improved numerical convergence for the family of CVD‐MPFA schemes using the quadrature parameterization has been observed for structured and unstructured grids in two dimensions. The CVD‐MPFA family cell‐vertex formulation is extended to classical general element types in 3‐D including prisms, pyramids, hexahedra and tetrahedra. A numerical convergence study of the CVD‐MPFA schemes on general unstructured grids comprising of triangular elements in 2‐D and prismatic, pyramidal, hexahedral and tetrahedral shape elements in 3‐D is presented. Copyright © 2011 John Wiley & Sons, Ltd.     

On the stability of bubble functions and a stabilized mixed finite element formulation for the Stokes problem

In this paper we investigate the relationship between stabilized and enriched finite element formulations for the Stokes problem. We also present a new stabilized mixed formulation for which the stability parameter is derived purely by the method of weighted residuals. This new formulation allows equal‐order interpolation for the velocity and pressure fields. Finally, we show by counterexample that a direct equivalence between subgrid‐based stabilized finite element methods and Galerkin methods enriched by bubble functions cannot be constructed for quadrilateral and hexahedral elements using standard bubble functions. Copyright © 2008 John Wiley & Sons, Ltd.     

密闭腔体声-结构耦合系统的动力灵敏度分析

以密闭空腔为对象,开展了声-结构耦合系统的动力分析和灵敏度计算,为系统性态优化设计提供理论和算法基础。分别把结构和声场进行离散化,推导了声-结构耦合系统的有限元方程,求解了耦合系统的频率和声压级响应。在此基础上,以结构尺寸为设计变量,计算了耦合系统的固有频率和声压级响应的灵敏度,解决了声-结构耦合系统动力灵敏度的数值算法问题。     

螺旋槽干气密封槽形参数的协调优化

应用PH线性化方法、迭代法,近似求解了螺旋槽内稳态流动场的非线性雷诺方程,求得了气体动压和速度分布的解析解.利用多目标优化方法构建了气膜刚度与泄漏量之比的协调函数,并对该目标函数进行了近似求解,获得了最佳的螺旋槽几何参数值.     

An approximately divergence-free 9-node velocity element (with variations) for incompressible flows

Explicit basis functions are constructed for 9-node biquadratic velocity fields which guarantee that a weak form of the continuity equation is satisfied. The corresponding pressure approximations are either piecewise constant, piecewise linear or piecewise bilinear. These results are extended to give bases for bilinear velocity/piecewise constant pressure elements and also to some three-dimensional brick elements.     

Symmetry Reductions of Equations for Solute Transport in Soil

Solute transport in saturated soil is represented by anonlinear system consisting of a Fokker–Planck equation coupled toLaplace's equation. Symmetries, reductions and exact solutions are foundfor two dimensional transient solute transport through some nontrivialwedge and spiral steady water flow fields. In particular, the mostgeneral complex velocity potential is determined, such that the soluteequation admits a stretching group of transformations that wouldnormally be possessed by a point source solution.     

A hybrid vortex method for the simulation of three‐dimensional flows

This paper presents an integral vorticity method for solving three‐dimensional Navier–Stokes equations. A finite volume scheme is implemented to solve the vorticity transport equation, which is discretized on a structured hexahedral mesh. A vortex sheet algorithm is used to enforce the no‐slip boundary condition through a vorticity flux at the boundary. The Biot–Savart integral is evaluated to compute the velocity field, in conjunction with a fast algorithm based on multipole expansion. This method is applied to the simulation of uniform flow past a sphere. Copyright © 2007 John Wiley & Sons, Ltd.     

螺旋槽干气密封微尺度流动场的近似计算及其参数优化

应用PH线性化方法、迭代法,近似求解了螺旋槽内稳态微尺度流动场的非线性雷诺方程,求得了气体动压和速度分布的解析解。继而利用多目标优化方法构建了气膜刚度与泄漏量之比的协调函数,并对该目标函数进行了近似求解,获得了最佳的螺旋槽几何参数值。     

液压作用下套管柱临界屈曲载荷计算的能量法

推导了内外液压作用下套管柱微弯时的能量平衡公式,并计算了其屈曲临界载荷计算公式。得到的结果低于目前已发表的结果。内外液压对套管稳定性的影响相当于在套管底端作用一个附加的轴向力和一个沿轴线均布的线载荷。对含有内压的两端封堵的细长管柱,内压对临界屈曲载荷没有影响。结果可供油井设计及作业参考。     

A particle swarm optimization approach for hexahedral mesh smoothing

There are some approaches for all‐hexahedral mesh quality improvement by means of node‐movement while preserving the connectivity. Among these methods, the most easily implemented and well known one is the Laplacian smoothing method; however, for this method mesh quality improvement is not guaranteed in all cases, and this approach might cause inverted elements especially in concave regions. In this work, a method for the improvement of hexahedral mesh shape‐quality without causing inverted elements is proposed; which is based on optimization of an objective function calculated by means of the individual qualities of hexahedral elements in the mesh. The shape‐quality for each hexahedral element is defined via the condition number of the relevant element. The numerical optimization scheme is the particle swarm optimization method, which originated from observations related to the social behaviors of bird, insect, or fish colonies. The purpose of this paper is to discuss the applicability of this approach to mesh smoothing. Some examples are given in order to demonstrate the applicability. Copyright © 2008 John Wiley & Sons, Ltd.     

Experimental analysis of the swirling flow in a Francis turbine draft tube: Focus on radial velocity component determination

The draft tube of a hydraulic turbine is the component where the flow exiting the runner is decelerated, thereby converting the excess of kinetic energy into static pressure. In the case of machine refurbishment of an existing power plant, most of the time only the runner and the guide vanes are currently modified. For financial and safety reasons, the spiral casing and the draft tube are seldom redesigned, even if these components present some undesirable behaviour. In some cases, the installation of an upgraded runner leads to a peculiar and undesirable efficiency drop as the discharge is increased above the best efficiency point value. It is found to be related to a corresponding sudden variation in the draft tube pressure recovery coefficient at the same discharge.The swirling flow exiting the runner is complex and highly turbulent. The radial velocity is rarely measured because a quite complicated measurement setup is needed. However, this velocity component is greatly needed in order to properly initialize the numerical simulations, and its influence is important in spite of its small magnitude. Velocity measurements downstream of the runner include radial component made at CREMHyG (Grenoble) by LDV, and PIV techniques are presented. An analytical formulation for this velocity component based on the formulation for the conical diffuser and on the three vortices structure is proposed and compared with measurements.     

A stabilized finite element method for the Stokes problem based on polynomial pressure projections

A new stabilized finite element method for the Stokes problem is presented. The method is obtained by modification of the mixed variational equation by using local L² polynomial pressure projections. Our stabilization approach is motivated by the inherent inconsistency of equal‐order approximations for the Stokes equations, which leads to an unstable mixed finite element method. Application of pressure projections in conjunction with minimization of the pressure–velocity mismatch eliminates this inconsistency and leads to a stable variational formulation. Unlike other stabilization methods, the present approach does not require specification of a stabilization parameter or calculation of higher‐order derivatives, and always leads to a symmetric linear system. The new method can be implemented at the element level and for affine families of finite elements on simplicial grids it reduces to a simple modification of the weak continuity equation. Numerical results are presented for a variety of equal‐order continuous velocity and pressure elements in two and three dimensions. Copyright © 2004 John Wiley & Sons, Ltd.