Consistency of iterative operator‐splitting methods: Theory and applications

In this article, we describe a different operator‐splitting method for decoupling complex equations with multidimensional and multiphysical processes for applications for porous media and phase‐transitions. We introduce different operator‐splitting methods with respect to their usability and applicability in computer codes. The error‐analysis for the iterative operator‐splitting methods is discussed. Numerical examples are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

On the convergence of operator splitting for the Rosenau–Burgers equation

We present convergence analysis of operator splitting methods applied to the nonlinear Rosenau–Burgers equation. The equation is first splitted into an unbounded linear part and a bounded nonlinear part and then operator splitting methods of Lie‐Trotter and Strang type are applied to the equation. The local error bounds are obtained by using an approach based on the differential theory of operators in Banach space and error terms of one and two‐dimensional numerical quadratures via Lie commutator bounds. The global error estimates are obtained via a Lady Windermere's fan argument. Lastly, a numerical example is studied to confirm the expected convergence order.     

Efficient numerical methods for pricing American options under stochastic volatility

Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two‐dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M‐matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one‐dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid‐based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower. ©John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

An efficient method for fractional nonlinear differential equations by quasi‐Newton's method and simplified reproducing kernel method

An efficient method for nonlinear fractional differential equations is proposed in this paper. This method consists of 2 steps. First, we linearize the nonlinear operator equation by quasi‐Newton's method, which is based on Fréchet derivative. Then we solve the linear fractional differential equations by the simplified reproducing kernel method. The convergence of the quasi‐Newton's method is discussed for the general nonlinear case as well. Finally, some numerical examples are presented to illustrate accuracy, efficiency, and simplicity of the method.     

The modulus‐based matrix splitting algorithms for a class of weakly nonlinear complementarity problems

In this paper, we study a class of weakly nonlinear complementarity problems arising from the discretization of free boundary problems. By reformulating the complementarity problems as implicit fixed‐point equations based on splitting of the system matrices, we propose a class of modulus‐based matrix splitting algorithms. We show their convergence by assuming that the system matrix is positive definite. Moreover, we give several kinds of typical practical choices of the modulus‐based matrix splitting iteration methods based on the different splitting of the system matrix. Numerical experiments on two model problems are presented to illustrate the theoretical results and examine the numerical effectiveness of our modulus‐based matrix splitting algorithms. Copyright © 2016 John Wiley & Sons, Ltd.     

Generalized skew‐Hermitian triangular splitting iteration methods for saddle‐point linear systems

A generalized skew‐Hermitian triangular splitting iteration method is presented for solving non‐Hermitian linear systems with strong skew‐Hermitian parts. We study the convergence of the generalized skew‐Hermitian triangular splitting iteration methods for non‐Hermitian positive definite linear systems, as well as spectrum distribution of the preconditioned matrix with respect to the preconditioner induced from the generalized skew‐Hermitian triangular splitting. Then the generalized skew‐Hermitian triangular splitting iteration method is applied to non‐Hermitian positive semidefinite saddle‐point linear systems, and we prove its convergence under suitable restrictions on the iteration parameters. By specially choosing the values of the iteration parameters, we obtain a few of the existing iteration methods in the literature. Numerical results show that the generalized skew‐Hermitian triangular splitting iteration methods are effective for solving non‐Hermitian saddle‐point linear systems with strong skew‐Hermitian parts. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical approximation of the conservative Allen–Cahn equation by operator splitting method

In this paper, a second‐order fast explicit operator splitting method is proposed to solve the mass‐conserving Allen–Cahn equation with a space–time‐dependent Lagrange multiplier. The space–time‐dependent Lagrange multiplier can preserve the volume of the system and keep small features. Moreover, we analyze the discrete maximum principle and the convergence rate of the fast explicit operator splitting method. The proposed numerical scheme is of spectral accuracy in space and of second‐order accuracy in time, which greatly improves the computational efficiency. Numerical experiments are presented to confirm the accuracy, efficiency, mass conservation, and stability of the proposed method. Copyright © 2017 John Wiley & Sons, Ltd.     

Splitting method for an inverse source problem in parabolic differential equations: Error analysis and applications

In this work, we present a numerical method based on a splitting algorithm to find the solution of an inverse source problem with the integral condition. The source term is reconstructed by using the specified data and by employing the Lie splitting method, we decompose the equation into linear and nonlinear parts. Each subproblem is solved by the Fourier transform and then by combining the solutions of subproblems, the solution of the original problem is computed. Moreover, the framework of strongly continuous semigroup (or C₀-semigroup) is employed in error analysis of operator splitting method for the inverse problem. The convergence of the proposed method is also investigated and proved. Finally, some numerical examples in one, two, and three-dimensional spaces are provided to confirm the efficiency and capability of our work compared with some other well-known methods.     

A domain decomposition method based on the iterative operator splitting method

In this article a new approach is proposed for constructing a domain decomposition method based on the iterative operator splitting method. The convergence properties of such a method are studied. The main feature of the proposed idea is the decoupling of space and time. We present a multi-iterative operator splitting method that combines iteratively the space and time splitting. We confirm with numerical applications the effectiveness of the proposed iterative operator splitting method in comparison with the classical Schwarz waveform relaxation method as a standard method for domain decomposition. We provide improved results and convergence rates.     

On convergence and complexity of the modified forward‐backward method involving new linesearches for convex minimization

In optimization theory, convex minimization problems have been intensively investigated in the current literature due to its wide range in applications. A major and effective tool for solving such problem is the forward‐backward splitting algorithm. However, to guarantee the convergence, it is usually assumed that the gradient of functions is Lipschitz continuous and the stepsize depends on the Lipschitz constant, which is not an easy task in practice. In this work, we propose the modified forward‐backward splitting method using new linesearches for choosing suitable stepsizes and discuss the convergence analysis including its complexity without any Lipschitz continuity assumption on the gradient. Finally, we provide numerical experiments in signal recovery to demonstrate the computational performance of our algorithm in comparison to some well‐known methods. Our reports show that the proposed algorithm has a good convergence behavior and can outperform the compared methods.     

Analysis of the operator splitting scheme for the Cahn‐Hilliard equation with a viscosity term

In this paper, we consider a second‐order fast explicit operator splitting method for the viscous Cahn‐Hilliard equation, which includes a viscosity term αΔu_t (α ∈ (0, 1)) described the influences of internal micro‐forces. The choice α = 0 corresponds to the classical Cahn‐Hilliard equation whilst the choice α = 1 recovers the nonlocal Allen‐Cahn equation. The fundamental idea of our method is to split the original problem into linear and nonlinear parts. The linear subproblem is numerically solved using a pseudo‐spectral method, and thus an ordinary differential equation is obtained. The nonlinear one is solved via TVD‐RK method. The stability and convergence are discussed in L²‐norm. Numerical experiments are performed to validate the accuracy and efficiency of the proposed method. Besides, a detailed comparison is made for the dynamics and the coarsening process of the metastable pattern for various values of α. Moreover, energy degradation and mass conservation are also verified.     

An efficient pseudo‐spectral Legendre–Galerkin method for solving a nonlinear partial integro‐differential equation arising in population dynamics

The pseudo‐spectral Legendre–Galerkin method (PS‐LGM) is applied to solve a nonlinear partial integro‐differential equation arising in population dynamics. This equation is a competition model in which similar individuals are competing for the same resources. It is a kind of reaction–diffusion equation with integral term corresponding to nonlocal consumption of resources. The proposed method is based on the Legendre–Galerkin formulation for the linear terms and interpolation operator at the Chebyshev–Gauss–Lobatto (CGL) points for the nonlinear terms. Also, the integral term, which is a kind of convolution, is directly computed by a fast and accurate method based on CGL interpolation operator, and thus, the use of any quadrature formula in its computation is avoided. The main difference of the PS‐LGM presented in the current paper with the classic LGM is in treating the nonlinear terms and imposing boundary conditions. Indeed, in the PS‐LGM, the nonlinear terms are efficiently handled using the CGL points, and also the boundary conditions are imposed strongly as collocation methods. Combination of the PS‐LGM with a semi‐implicit time integration method such as second‐order backward differentiation formula and Adams‐Bashforth method leads to reducing the complexity of computations and obtaining a linear algebraic system of equations with banded coefficient matrix. The desired equation is considered on one and two‐dimensional spatial domains. Efficiency, accuracy, and convergence of the proposed method are demonstrated numerically in both cases. Copyright © 2012 John Wiley & Sons, Ltd.     

Diagonal and Toeplitz splitting iteration methods for diagonal‐plus‐Toeplitz linear systems from spatial fractional diffusion equations

The finite difference discretization of the spatial fractional diffusion equations gives discretized linear systems whose coefficient matrices have a diagonal‐plus‐Toeplitz structure. For solving these diagonal‐plus‐Toeplitz linear systems, we construct a class of diagonal and Toeplitz splitting iteration methods and establish its unconditional convergence theory. In particular, we derive a sharp upper bound about its asymptotic convergence rate and deduct the optimal value of its iteration parameter. The diagonal and Toeplitz splitting iteration method naturally leads to a diagonal and circulant splitting preconditioner. Analysis shows that the eigenvalues of the corresponding preconditioned matrix are clustered around 1, especially when the discretization step‐size h is small. Numerical results exhibit that the diagonal and circulant splitting preconditioner can significantly improve the convergence properties of GMRES and BiCGSTAB, and these preconditioned Krylov subspace iteration methods outperform the conjugate gradient method preconditioned by the approximate inverse circulant‐plus‐diagonal preconditioner proposed recently by Ng and Pan (M.K. Ng and J.‐Y. Pan, SIAM J. Sci. Comput. 2010;32:1442‐1464). Moreover, unlike this preconditioned conjugate gradient method, the preconditioned GMRES and BiCGSTAB methods show h‐independent convergence behavior even for the spatial fractional diffusion equations of discontinuous or big‐jump coefficients.     

High‐order characteristic‐tracking strategy for simulation of a nonlinear advection–diffusion equation

In this study, new high‐order backward semi‐Lagrangian methods are developed to solve nonlinear advection–diffusion type problems, which are realized using high‐order characteristic‐tracking strategies. The proposed characteristic‐tracking strategies are second‐order L‐stable and third‐order L(α)‐stable methods, which are based on a classical implicit multistep method combined with a error‐correction method. We also use backward differentiation formulas and the fourth‐order finite‐difference scheme for diffusion problem discretization in the temporal and spatial domains, respectively. To demonstrate the adaptability and efficiency of these time‐discretization strategies, we apply these methods to nonlinear advection–diffusion type problems such as the viscous Burgers' equation. Through simulations, not only the temporal and spatial accuracies are numerically evaluated but also the proposed methods are shown to be superior to the compared existing characteristic‐tracking methods under the same rates of convergence in terms of accuracy and efficiency. Finally, we have shown that the proposed method well preserves the energy and mass when the viscosity coefficient becomes zero.     

一类求解单调变分不等式的隐式方法

1．引言变分不等式是一个非常有趣。非常困难的数学问题［"］．它具有广泛的应用（例如，数学规划中的许多基本问题都可以归结为一个变分不等式问题），因而得到深入的研究并有了不少算法［1，2，5－8，17－21］．对线性单调变分不等式，我们最近提出了一系列投影收缩算法Ig－13］．本文考虑求解单调变分不等式其中0CW是一闭凸集，F是从正p到自身的一个单调算子，一即有我们用比（·）表示到0上的投影．求解单调变分不等式的一个简单方法是基本投影法［1，6］，它的迭代式为然而，如果F不是仿射函数，只有当F一致强单调且LIPSChitZ连续…     

Splitting methods for Hamilton‐Jacobi equations

We explain how the exploitation of several kinds of operator splitting methods, both local and global in time, lead to simple numerical schemes approximating the solution of nonlinear Hamilton‐Jacobi equations. We review the existing local methods which have been used since the early 80's and we introduce a new method which is global in time. We show some numerical experiments. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Extended Projection Methods for Monotone Variational Inequalities

In this paper, we prove that each monotone variational inequality is equivalent to a two-mapping variational inequality problem. On the basis of this fact, a new class of iterative methods for the solution of nonlinear monotone variational inequality problems is presented. The global convergence of the proposed methods is established under the monotonicity assumption. The conditions concerning the implementability of the algorithms are also discussed. The proposed methods have a close relationship to the Douglas–Rachford operator splitting method for monotone variational inequalities.     

Parameter estimation in the Hermitian and skew‐Hermitian splitting method using gradient iterations

This article presents enhancement strategies for the Hermitian and skew‐Hermitian splitting method based on gradient iterations. The spectral properties are exploited for the parameter estimation, often resulting in a better convergence. In particular, steepest descent with early stopping can generate a rough estimate of the optimal upper bound. This is better than an arbitrary choice since the latter often causes stability problems or slow convergence. In addition, delayed gradient methods are considered as inner solvers for the splitting method. Experiments verify the effectiveness of the proposed estimation strategies and show that delayed gradient methods are competitive with conjugate gradient in low precision.     

An optimal sixth‐order quartic B‐spline collocation method for solving Bratu‐type and Lane‐Emden–type problems

This paper is concerned with the numerical solutions of Bratu‐type and Lane‐Emden–type boundary value problems, which describe various physical phenomena in applied science and technology. We present an optimal collocation method based on quartic B‐spine basis functions to solve such problems. This method is constructed by perturbing the original problem and on a uniform mesh. The method has been tested by four nonlinear examples. In order to show the advantage of the new method, numerical results are compared with those obtained by some of the existing methods, such as normal quartic B‐spline collocation method and the finite difference method (FDM). It has been observed that the order of convergence of the proposed method is six, which is two orders of magnitude larger than the normal quartic B‐spline collocation method. Moreover, our method gives highly accurate results than the FDM.     

Modelling and numerical approach to a class of metal‐forming problems—Quasi‐steady case

A class of quasi‐steady metal‐forming problems, with rigid‐plastic, incompressible, strain and strain‐rate dependent material model and with unilateral frictionless and nonlinear, nonlocal Coulomb's frictional contact conditions is considered. A coupled variational formulation, constituted of a variational inequality, with nonlinear and nondifferentiable terms, and a strain evolution equation, is derived and under a restriction on the material characteristics and using a variable stiffness parameters method with time retardation, existence, uniqueness and convergence results are obtained and presented. An algorithm, combining this method and the finite element method, is proposed and applied for solving an example strip drawing problem. Copyright © 2011 John Wiley & Sons, Ltd.