The mixed finite element multigrid method for stokes equations.

The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q2-Q1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results.

On a cubically convergent iterative method for matrix sign.

We propose an iterative method for finding matrix sign function. It is shown that the scheme has global behavior with cubical rate of convergence. Examples are included to show the applicability and efficiency of the proposed scheme and its reciprocal.

On a derivative-free variant of king's family with memory.

The aim of this paper is to construct a method with memory according to King's family of methods without memory for nonlinear equations. It is proved that the proposed method possesses higher R-order of convergence using the same number of functional evaluations as King's family. Numerical experiments are given to illustrate the performance of the constructed scheme.

On a new iterative scheme without memory with optimal eighth order.

The purpose of this paper is to derive and discuss a three-step iterative expression for solving nonlinear equations. In fact, we derive a derivative-free form for one of the existing optimal eighth-order methods and preserve its convergence order. Theoretical results will be upheld by numerical experiments.

On a new three-step class of methods and its acceleration for nonlinear equations.

A class of derivative-free methods without memory for approximating a simple zero of a nonlinear equation is presented. The proposed class uses four function evaluations per iteration with convergence order eight. Therefore, it is an optimal three-step scheme without memory based on Kung-Traub conjecture. Moreover, the proposed class has an accelerator parameter with the property that it can increase the convergence rate from eight to twelve without any new functional evaluations. Thus, we construct a with memory method that increases considerably efficiency index from 8(1/4) ≈ 1.681 to 12(1/4) ≈ 1.861. Illustrations are also included to support the underlying theory.

Numerical Investigation of the Effect of Unsteadiness on Three-Dimensional Flow of an Oldroyb-B Fluid.

A spectral relaxation method used with bivariate Lagrange interpolation is used to find numerical solutions for the unsteady three-dimensional flow problem of an Oldroyd-B fluid with variable thermal conductivity and heat generation. The problem is governed by a set of three highly coupled nonlinear partial differential equations. The method, originally used for solutions of systems of ordinary differential equations is extended to solutions of systems of nonlinear partial differential equations. The modified approach involves seeking solutions that are expressed as bivariate Lagrange interpolating polynomials and applying pseudo-spectral collocation in both independent variables of the governing PDEs. Numerical simulations were carried out to generate results for some of the important flow properties such as the local skin friction and the heat transfer rate. Numerical analysis of the error and convergence properties of the method are also discussed. One of the benefits of the proposed method is that it is computationally fast and gives very accurate results after only a few iterations using very few grid points in the numerical discretization process.