Iterative solvers for quadratic discretizations of the generalized Stokes problem

We present in this paper various iterative methods for the solution of large linear and non‐linear systems resulting from the discretization of the generalized Stokes problem. A second‐order (O(h²)) P₂‐P₁ mixed finite element is used for the approximation of the velocity and the pressure. Solution strategies based on conjugate gradient‐like methods, the Uzawa's and Arrow–Hurwicz's methods are presented. Schur complement methods are also explored in the context of a hierarchical decomposition of the velocity field. The ever present preconditioning problem is also addressed. The performance of these iterative methods will be discussed on complex flows of industrial interest. Copyright © 2004 John Wiley & Sons, Ltd.     

Segregated finite element algorithms for the numerical solution of large-scale incompressible flow problems

This paper presents results of an ongoing research program directed towards developing fast and efficient finite element solution algorithms for the simulation of large-scale flow problems. Two main steps were taken towards achieving this goal. The first step was to employ segregated solution schemes as opposed to the fully coupled solution approach traditionally used in many finite element solution algorithms. The second step was to replace the direct Gaussian elimination linear equation solvers used in the first step with iterative solvers of the conjugate gradient and conjugate residual type. The three segregated solution algorithms developed in step one are first presented and their integrity and relative performance demonstrated by way of a few examples. Next, the four types of iterative solvers (i.e. two options for solving the symmetric pressure type equations and two options for solving the non-symmetric advection–diffusion type equations resulting from the segregated algorithms) together with the two preconditioning strategies employed in our study are presented. Finally, using examples of practical relevance the paper documents the large gains which result in computational efficiency, over fully coupled solution algorithms, as each of the above two main steps are introduced. It is shown that these gains become increasingly more dramatic as the complexity and size of the problem is increased.     

Melnikov's method for chaos of a two-dimensional thin panel in subsonic flow with external excitation

In this brief communication, Melnikov's method is adopted to study the chaotic behaviors of a two-dimensional thin panel subjected to subsonic flow and external excitation. The nonlinear governing equations of the subsonic panel system are reduced to a series of ordinary differential equations by using Galerkin method. The critical parameters for chaos are obtained. It is found that the critical parameters obtained by the theoretical analysis are in agreement with the numerical simulations. The method suggested in this paper can also be extended for other fluid-structure dynamic systems, such as the fluid-conveying system.     

自由面势流问题边界元法的数值色散误差

以连续及离散Fourier分析研究自由面势流问题边界元法的数值色散误差,并从理论上探讨有关计算中数值色散误差的改善问题.研究表明:对于该问题的数值色散误差而言,重要的在于以问题相应的离散算子考察计及各种数值手段后的总体色散误差,而非仅考虑该数值手段自身的数值色散误差大小.高阶面元、自由面域外奇点或适当的耦合方法是降低有关问题算子总体色散误差的较好选择.     

A higher-order boundary element method for three-dimensional potential problems

A method of eliminating the singularities involved in boundary element methods for three-dimensional potential problems is presented and the non-singular expressions of integrals on an element on which the singular point is situated are given for linear and quadratic interpolation functions. Numerical examples are compared with analytical solutions to show that the higher-order interpolations have better precision.     

Local pressure preconditioning method for steady incompressible flows

The convergence and accuracy characteristics of the preconditioned incompressible Euler and Navier–Stokes equations are studied. An object-oriented C++ numerical code has been developed for solving the inviscid and viscous, steady, incompressible flows problems. The code is based on the cell-centred finite volume method. In this scheme, two-dimensional incompressible Euler and Navier–Stokes equations are modified by a robust artificial compressibility (AC) and a local preconditioning matrix of pressure-sensor type. The preconditioned equations are solved with the Jameson's numerical approach, i.e. artificial dissipation and artificial viscosity terms under the form of a fourth- and second-order derivative, respectively. An explicit four-stage Runge–Kutta integration algorithm is applied to obtain the steady-state condition. The computed results include the steady-state solution of flow past the NACA-hydrofoils and a circular cylinder in free stream, for which the numerical results are compared with numerical works of other researchers. Good agreement is observed. The effects of AC parameter, artificial viscosity and dissipation factor, and local preconditioning coefficient on convergence rate and solution accuracy are tested by computing flow over the NACA0012 hydrofoil. In addition, some important design criteria of a preconditioner, such as stiffness reduction, hyperbolicity, symmetrisability, accuracy preservation for M → 0, and M-property have been examined analytically.     

The simple layer potential method of fundamental solutions for certain biharmonic problems

A novel formulation of the method of fundamental solutions for the numerical solution of plane biharmonic problems, based on the simple layer potential representation of Fichera, is presented. The applicability and accuracy of the method are demonstrated by examining its performance on a set of practical problems arising in Stokes fluid flow.     

无穷域势流问题的有限元／差分线法混合求解

提出了一种将有限元和差分线法相结合求解无穷域势流问题的算法。用两同心圆将求解域划分为存在重叠的有限和无限两个区域,在有限和无限域上分别用有限元和差分线法求解Laplace方程边值问题。用差分线法推导出的关系式修正有限元方程,求解该方程组从而得到原问题的解。本算法将求解无穷域问题转化为代数特征值问题和有限域内线性方程组的...     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

THE SOLUTION FOR TRANSIENT TWO-PHASE FLOW BY SPLIT FLUX VECTOR METHOD

In this paper the transient two-phase flow equations and their eigenvalues are first intreduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difJerence operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.     

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。     

The lift computation for an oscillating flat plate in incompressible potential flow

The initial aim of this work was the estimation of the lift acting on a flat plate performing small oscillations in a plane uniform stream by means of a simplified model based on one or at the most two lumped vortices, and the assessment of its results by comparison to those that were exact. The model was found to work well up to a reduced frequency of about 1 or 2, above which the results diverged from those that were correct. In order to improve the model, its behaviour at very high frequencies was then investigated, discovering: (i) that if the number of lumped vortices is greater than one the possibility to impose all boundary conditions is subject to certain geometrical constraints; (ii) that the asymptotical behaviour is not the right one. A straightforward extension of this conclusion to the exact case of a continuous, vorticity distribution simulating the motion of the plate and to the classical equation describing it leads apparently to an incorrect result. The reason for the discrepancy is found in the singularity displayed by the integral equation which cannot be reproduced by the discrete model. It this therefore concluded that the latter can be trusted at low and middle frequencies but its extension to higher ones is fundamentally uncorrect.

---

Sommario Lo scopo iniziale di questo lavoro era il calcolo approssimato della portanza agente su una lamina piana soggetta a piccole oscillazioni in una corrente bidimensionale uniforme per mezzo di un modello semplificato basato su uno o al più due vortici concentrati, e il confronto dei risultati con quelli esatti. Il modello risulta funzionare bene per frequenze ridotte inferiori ad 1 o 2, sopra le quali, tuttavia, i risultati si allontanano da quelli corretti. Per migliorarlo si è allora studiato il suo comportamento alle frequenze molto alte, scoprendo che: (i) la possibilità di imporre tutte le condizioni al contorno quando il numero dei vortici concentrati è superiore a uno è soggetta a certe limitazioni sulla configurazione geometrica; (ii) che il comportamento asindotico non è quello corretto. Un'estensione automatica di questa conclusione al caso esatto in cui il moto della lamina è simulato da una distribuzione continua di vorticità e alla classica equazione che lo descrive sembra condurre ad un risultato errato. La ragione di questa discrepanza viene individuata nella singolarità contenuta nell'equazione integrale, che non può essere riprodotta dal modello discreto. Se ne conclude perciò che esso è utilizzabile alle basse e medie frequenze ma che una sua estensione alle alte è fondamentalmente errata.

     

Application of the finite point method to high‐Reynolds number compressible flow problems

In this work, the finite point method is applied to the solution of high‐Reynolds compressible viscous flows. The aim is to explore this important field of applications focusing on two main aspects: the easiness and automation of the meshless discretization of viscous layers and the construction of a robust numerical approximation in the highly stretched clouds of points resulting in such domain areas. The flow solution scheme adopts an upwind‐biased scheme to solve the averaged Navier–Stokes equations in conjunction with an algebraic turbulence model. The numerical applications presented involve different attached boundary layer flows and are intended to show the performance of the numerical technique. The results obtained are satisfactory and indicative of the possibilities to extend the present meshless technique to more complex flow problems. Copyright © 2014 John Wiley & Sons, Ltd.     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

On shock layer method for the hypersonic flow problems

Chernyi’s series method[1] is not proper for the case that(γ-l)/(γ+l)<<2/(γ+1)×M²sin²β (γ=c_p/c_v-adiabatic index number, M-Much number, β-shock incidence). In this paper, we only suppose that in the neighbour of the shock, there exists a shock layer in which the density of the gas is very big, but we do not remove the case that (γ-1)/(γ+1)<<2/(γ+1)M²sin²β.     

Comparison of three solvers for viscoelastic fluid flow problems

The computational efficiency of three numerical schemes has been examined for the solution of a linearized system of equations resulting from the finite element discretization of a viscoelastic fluid flow problem. The first scheme is a modified frontal solver, which solves the linear system of equations directly. The other two, one based on a biconjugate gradient stabilized (BiCGStab) method and another based on a generalized minimal residual (GMRES) method, are iterative schemes. The stick-slip problem and the four-to-one contraction problem were analyzed and the viscoelastic fluid was assumed to obey the Oldroyd-B model. The two iterative schemes are superior to the direct scheme in terms of CPU time consumed and the BiCGStab scheme is even faster than the GMRES scheme. The range of convergence for both iterative schemes is compatible with that of the direct scheme.     

A three-dimensional fictitious domain method for incompressible fluid flow problems

A new Galerkin finite element method for the solution of the Navier–Stokes equations in enclosures containing internal parts which may be moving is presented. Dubbed the virtual finite element method, it is based upon optimization techniques and belongs to the class of fictitious domain methods. Only one volumetric mesh representing the enclosure without its internal parts needs to be generated. These are rather discretized using control points on which kinematic constraints are enforced and introduced into the mathematical formulation by means of Lagrange multipliers. Consequently, the meshing of the computational domain is much easier than with classical finite element approaches. First, the methodology will be presented in detail. It will then be validated in the case of the two-dimensional Couette cylinder problem for which an analytical solution is available. Finally, the three-dimensional fluid flow inside a mechanically agitated vessel will be investigated. The accuracy of the numerical results will be assessed through a comparison with experimental data and results obtained with a standard finite element method. © 1997 John Wiley & Sons, Ltd.     

A total linearization method for solving viscous free boundary flow problems by the finite element method

In this paper a total linearization method is derived for solving steady viscous free boundary flow problems (including capillary effects) by the finite element method. It is shown that the influence of the geometrical unknown in the totally linearized weak formulation can be expressed in terms of boundary integrals. This means that the implementation of the method is simple. Numerical experiments show that the iterative method gives accurate results and converges very fast.