Comparison of variants of the bi-conjugate gradient method for compressible Navier-Stokes solver with second-moment closure

Variants of the bi-conjugate gradient (Bi-CG) method are used to resolve the problem of slow convergence in CFD when it is applied to complex flow field simulation using higher-order turbulence models. In this study the Navier-Stokes and Reynolds stress transport equations are discretized with an implicit, total variation diminishing (TVD), finite volume formulation. The preconditioning technique of incomplete lower-upper (ILU) factorization is incorporated into the conjugate gradient square (CGS), bi-conjugate gradient stable (Bi-CGSTAB) and transpose-free quasi-minimal residual (TFQMR) algorithms to accelerate convergence of the overall itertive methods. Computations have been carried out for separated flow fields over transonic bumps, supersonic bases and supersonic compression corners. By comparisons of the convergence rate with each other and with the conventional approximate factorization (AF) method it is shown that the Bi-CGSTAB method gives the most efficient convergence rate among these methods and can speed up the CPU time by a factor of 2·4–6·5 as compared with the AF method. Moreover, the AF method may yield somewhat different results from variants of the Bi-CG method owing to the factorization error which introduces a higher level of convergence criterion.     

二维定常湍流计算中的GMRES算法

在以前工作的基础上将广义极小残差（Generalized Minimum RESidual (GMRES)算法发展到用于求解二维可压Favier平均Navier-Stokes方程组。控制方程经Newton线化处理后构成近似的线性系统,然后采用分别耦合了LUSGS和ILU两种预处理矩阵的GMRES算法求解。Spalart-Allmaras湍流模型被用来封闭流体控制方程组,采用与流体控制方程非耦合的方式,使用LUSGS方法求解。对GMRES算法中矩阵－向量的乘积采用了有限插分方法,从而避免了精确的左端系数矩阵的计算和存储。对预处理矩阵的两种使用方法（左预处理和右预处理）进行了分析和讨论。用两个算例对LUSGS和ILU两各预处理矩阵进行了比较,同时探讨了左预处理和右预处理各自的优缺点。通过对Sajben扩压器和NACA0012有攻角流动的计算,表明带有预处理的GMRES算法在二维定常跨音黏性流动计算中相比于得到广泛应用的DDADI方法具有很大优势,左预处理要优于右预处理。     

A NON-LINEAR ADAPTIVE FULL TRI-TREE MULTIGRID METHOD FOR THE MIXED FINITE ELEMENT FORMULATION OF THE NAVIER–STOKES EQUATIONS

The full adaptive multigrid method is based on the tri-tree grid generator. The solution of the Navier–Stokes equations is first found for a low Reynolds number. The velocity boundary conditions are then increased and the grid is adapted to the scaled solution. The scaled solution is then used as a start vector for the multigrid iterations. During the multigrid iterations the grid is first recoarsed a specified number of grid levels. The solution of the Navier–Stokes equations with the multigrid residual as right-hand side is smoothed in a fixed number of Newton iterations. The linear equation system in the Newton algorithm is solved iteratively by CGSTAB preconditioned by ILU factorization with coupled node fill-in. The full adaptive multigrid algorithm is demonstr ated for cavity flow. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 1037ndash;1047, 1997.     

Inexact Newton method via Lanczos decomposed technique for solving box-constrained nonlinear systems

This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.     

Parallelization of the ILU(0) preconditioner for CFD problems on shared‐memory computers

The use of ILU(0) factorization as a preconditioner is quite frequent when solving linear systems of CFD computations. This is because of its efficiency and moderate memory requirements. For a small number of processors, this preconditioner, parallelized through coloring methods, shows little savings when compared with a sequential one using adequate reordering of the unknowns. Level scheduling techniques are applied to obtain the same preconditioning efficiency as in a sequential case, while taking advantage of parallelism through block algorithms. Numerical results obtained from the parallel solution of the compressible Navier–Stokes equations show that this technique gives interesting savings in computational times on a small number of processors of shared‐memory computers. In addition, it does this while keeping all the benefits of an ILU(0) factorization with an adequate reordering of the unknowns, and without the loss of efficiency of factorization associated with a more scalable coloring strategy. Copyright © 1999 John Wiley & Sons, Ltd.     

A Newton–Krylov finite volume algorithm for the power-law non-Newtonian fluid flow using pseudo-compressibility technique

An implicit Newton–Krylov finite volume algorithm has been developed for efficient steady-state computation of the power-law non-Newtonian fluid flows. The pseudo-compressibility technique is used for the coupling of continuity and momentum equations. The spatial discretization is central (second-order) for both convective and diffusive terms and the accuracy of the solution is verified. The nine block diagonal Jacobian matrix (needed for implicit formulation) is computed directly through the flux differentiation. Five-diagonal and three-diagonal block matrices (the simplified versions of the main Jacobian matrix) are used with the ILU(0 & 1) and the Thomas linear solvers for preconditioning, respectively. The performance of the Newton-GMRES solver is examined in detail for different preconditioning strategies. The effects of the power-law behavior index and Re number on the convergence rate are also studied. The performance of the Newton-BiCGSTAB and the Newton-GMRES solvers are compared with each other. The results show, the ILU(1)/Newton-GMRES is the most efficient combination that is robust even in high Reynolds number shear-thinning fluid flow cases.     

Reordering and incomplete preconditioning in serial and parallel adaptive mesh refinement and coarsening flow solutions

The effects of reordering the unknowns on the convergence of incomplete factorization preconditioned Krylov subspace methods are investigated. Of particular interest is the resulting preconditioned iterative solver behavior when adaptive mesh refinement and coarsening (AMR/C) are utilized for serial or distributed parallel simulations. As representative schemes, we consider the familiar reverse Cuthill–McKee and quotient minimum degree algorithms applied with incomplete factorization preconditioners to CG and GMRES solvers. In the parallel distributed case, reordering is applied to local subdomains for block ILU preconditioning, and subdomains are repartitioned dynamically as mesh adaptation proceeds. Numerical studies for representative applications are conducted using the object‐oriented AMR/C software system libMesh linked to the PETSc solver library. Serial tests demonstrate that global unknown reordering and incomplete factorization preconditioning can reduce the number of iterations and improve serial CPU time in AMR/C computations. Parallel experiments indicate that local reordering for subdomain block preconditioning associated with dynamic repartitioning because of AMR/C leads to an overall reduction in processing time. Copyright © 2011 John Wiley & Sons, Ltd.     

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, non-linear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these non-linear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modelled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.     

A numerical eigenvalue study of preconditioned non‐equilibrium transport equations

The finite element integration of non‐equilibrium contaminant transport in porous media yields sparse, unsymmetric, real or complex equations, which may be solved by iterative projection methods, such as Bi‐CGSTAB and TFQMR, on condition that they are effectively preconditioned. To ensure a fast convergence, the eigenspectrum of the preconditioned equations has to be very compact around unity. Compactness is generally measured by the spectral condition number. In difficult advection‐dominated problems, however, the condition number may be large and nevertheless, convergence may be good. A numerical study of the preconditioned eigenspectrum of a representative test case is performed using the incomplete triangular factorization. The results show that preconditioning eliminates most of the original complex eigenvalues, and that compactness is not necessarily jeopardized by a large condition number. Quite surprisingly, it is shown that the preconditioned complex problem may have a more compact real eigenspectrum than the equivalent real problem. Copyright © 1999 John Wiley & Sons, Ltd.     

Parallel ILU preconditioning,a priori pivoting and segregation of variables for iterative solution of the mixed finite element formulation of the Navier–Stokes equations

A parallel ILU preconditioning algorithm for the incompressible Navier–Stokes equations has been designed, implemented and tested. The computational mesh is divided into N subdomains which are processed in parallel in different processors. During ILU factorization, matrices and vectors associated with the nodes on the interface between the subdomains are communicated to the equation matrices to the adjacent subdomain. The bases for the parallel algorithm are an appropriate node ordering scheme and a segregation of velocity and pressure degrees of freedom. The inner nodes of the subdomain are numbered first and then the nodes on the interface between the subdomains. To avoid division by zero during the ILU factorization, the equations corresponding to the velocity degrees of freedom are assembled first in the global equation matrix, followed by the equations corresponding to the pressure degrees of freedom. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical approximation of optimal control of unsteady flows using SQP and time decomposition

In this paper, we present numerical approximations of optimal control of unsteady flow problems using sequential quadratic programming method (SQP) and time domain decomposition. The SQP method is considered superior due to its fast convergence and its ability to take advantage of existing numerical techniques for fluid flow problems. It iteratively solves a sequence of linear quadratic optimal control problems converging to the solution of the non‐linear optimal control problem. The solution to the linear quadratic problem is characterized by the Karush–Kuhn–Tucker (KKT) optimality system which in the present context is a formidable system to solve. As a remedy various time domain decompositions, inexact SQP implementations and block iterative methods to solve the KKT systems are examined. Numerical results are presented showing the efficiency and feasibility of the algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

A preconditioned alternating inner-outer iterative solution method for the mixed finite element formulation of the Navier-Stokes equations

In the present work a new iterative method for solving the Navier-Stokes equations is designed. In a previous paper a coupled node fill-in preconditioner for iterative solution of the Navier-Stokes equations proved to increase the convergence rate considerably compared with traditional preconditioners. The further development of the present iterative method is based on the same storage scheme for the equation matrix as for the coupled node fill-in preconditioner. This storage scheme separates the velocity, the pressure and the coupling of pressure and velocity coefficients in the equation matrix. The separation storage scheme allows for an ILU factorization of both the velocity and pressure unknowns. With the inner-outer solution scheme the velocity unknowns are eliminated before the resulting equation system for the pressures is solved iteratively. After the pressure unknown has been found, the pressures are substituted into the original equation system and the velocities are also found iteratively. The behaviour of the inner-outer iterative solution algorithm is investigated in order to find optimal convergence criteria for the inner iterations and compared with the solution algorithm for the original equation system. The results show that the coupled node fill-in preconditioner of the original equation system is more efficient than the coupled node fill-in preconditioner of the reduced equation system. However, the solution technique of the reduced equation system revals properties which may be advantageous in future solution algorithms.     

Multilayer open flow model: A simple pressure correction method for wave problems

This paper describes a pressure correction method for single‐ and multilayer open flow models. The method does not require any complex procedures to solve the discretization of the Poisson equation and is distinguished by a high computational efficiency. The algorithm can easily be adapted to irregular meshes and parallelized. Parabolic interpolation of the pressure profile is used for the free surface. The discretization of the Poisson equation is written in a matrix form, allowing its usage also in the case of basic function expansion of the depth pressure profile. The paper presents the results of algorithm verification where experimental data sensitive to the numerical dissipation of the calculation model was used. Iteration convergence is high including problems with dry‐bed flooding. The complete described technique of pressure correction is implemented in OpenCL on the GPU. Computation time for a test problem solved using CPU and GPU is compared.     

用拟压缩性方法求解不可压Euler方程

用拟压缩性方法和Ｊａｍｅｓｏｎ的有限体积算法求解了二维和三维定常可可压Ｅｕｌｅｒ方程。分别采用显、隐式时间离散推进求解;分析了人工粘性的阶数对定常解收敛性的影响,应用该方法计算了单个翼型和翼身组合体的低速绕流,结果与实验吻合较好。     

A parallel finite element sliding mesh technique for the simulation of viscous flows in agitated tanks

A parallel sliding mesh algorithm for the finite element simulation of viscous fluid flows in agitated tanks is presented. Lagrange multipliers are used at the sliding interfaces to enforce the continuity between the fixed and moving subdomains. The novelty of the method consists of the coupled solution of the resulting velocity–pressure‐Lagrange multipliers system of equations by an ILU(0)‐QMR solver. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid pivoting problems in the ILU(0) algorithm. To handle the convective term, both the Newton–Raphson scheme and the semi‐implicit linearization are tested. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid the failure of the ILU(0) algorithm due to the lack of pivoting. Furthermore, this approach is versatile enough so that it allows partitioning of sliding and fixed subdomains if parallelization is required. Although the sliding mesh technique is fairly common in CFD, the main advantage of the proposed approach is its low computational cost due to the inexpensive and parallelizable calculations that involve preconditioned sparse iterative solvers. The method is validated for Couette and coaxial stirred tanks. Copyright © 2011 John Wiley & Sons, Ltd.     

结构性软土弹塑性模型的隐式算法实现

对于考虑软土结构性的高度非线性弹塑性本构模型,在采用Newton-CPPM隐式算法对模型进行数值实现的过程中容易出现Jacobian矩阵奇异和不收敛问题。为此,本文提出了两种改进隐式算法。考虑到Newton-CPPM隐式算法是局部收敛性算法,因此引入大范围收敛的同伦延拓算法对Newton-CPPM算法的迭代初值进行改进,形成了同伦-Newton-CPPM算法。考虑到Newton-CPPM隐式算法单个迭代步的计算量过大,因此借鉴显式算法的思想提出一种两阶段迭代算法,第一阶段先求出一致性参数,第二阶段采用类似于显示算法的方法进行回代得出状态变量的值。然后,以考虑软土结构性的SANICLAY模型为例,从弹塑性本构模型的组成和算法的特点两个角度分析了引起Jacobian矩阵奇异和不收敛问题的原因,并且在单单元计算的基础上,对全显式算法、传统隐式算法和两种改进隐式算法在计算收敛性、计算精度和计算效率方面进行了对比。最后,将同伦-Newton-CPPM算法和传统隐式算法用于地基承载力多单元计算中,结果表明该算法能够有效地解决Jacobian矩阵奇异和不收敛问题。      

Two‐level consistent splitting methods based on three corrections for the time‐dependent Navier–Stokes equations

Three kinds of two‐level consistent splitting algorithms for the time‐dependent Navier–Stokes equations are discussed. The basic technique of two‐level type methods for solving the nonlinear problem is first to solve a nonlinear problem in a coarse‐level subspace, then to solve a linear equation in a fine‐level subspace. Hence, the two‐level methods can save a lot of work compared with the one‐level methods. The approaches to linearization are considered based on Stokes, Newton, and Oseen corrections. The stability and convergence demonstrate that the two‐level methods can acquire the optimal accuracy with the proper choice of the coarse and fine mesh scales. Numerical examples show that Stokes correction is the simplest, Newton correction has the best accuracy, while Oseen correction is preferable for the large Reynolds number problems and the long‐time simulations among the three methods. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

A comparative study of pressure correction and block-implicit finite volume algorithms on parallel computers

A comparison of a new parallel block-implicit method and the parallel pressure correction procedure for the solution of the incompressible Navier–Stokes equations is presented. The block-implicit algorithm is based on a pressure equation. The system of non-linear equation s is solved by Newton's method. For the solution of the linear algebraic systems the Bi-CGSTAB algorithm with incomplete lower–upper (ILU) decomposition of the matrix is applied. Domain decomposition serves as a strategy for the parallelization of the algorithms. Different algorithms for the parallel solution of the linear system of algebraic equations in conjunction with the pressure correction procedure are proposed. Three different flows are predicted with the parallel algorithms. Results and efficiency data of the block-implicit method are compared with the parallel version of the pressure correction algorithm. The block-implicit method is characterized by stable convergence behaviour, high numerical efficiency, insensitivity to relaxation parameters and high spatial accuracy. © 1997 John Wiley & Sons, Ltd.     

A segregated implicit solution algorithm for 2D and 3D laminar incompressible flows

A segregated algorithm for the solution of laminar incompressible, two- and three-dimensional flow problems is presented. This algorithm employs the successive solution of the momentum and continuity equations by means of a decoupled implicit solution method. The inversion of the coefficient matrix which is common for all momentum equations is carried out through an approximate factorization in upper and lower triangular matrices. The divergence-free velocity constraint is satisfied by formulating and solving a pressure correction equation. For the latter a combined application of a preconditioning technique and a Krylov subspace method is employed and proved more effecient than the approximate factorization method. The method exhibits a monotonic convergence, it is not costly in CPU time per iteration and provides accurate solutions which are independent of the underrelaxation parameter used in the momentum equations. Results are presented in two- and three-dimensional flow problems.