Extended local fourier analysis for multigrid optimal smoothing, coarse grid correction, and preconditioning

Multigrid methods are fast iterative solvers for partial di erential equations. Especially for elliptic equations they have been proven to be highly e cient. For problems with nonelliptic and nonsymmetric features--as they often occur in typical real life applications--a rigorous mathematical theory is generally not available. For such situations Fourier smoothing and two-grid analysis can be considered as the main analysis tools to obtain quantitative convergence estimates and to optimize different multigrid components like smoothers or inter-grid transfer operators. In general, it is difficult to choose the correct multigrid components for large classes of problems. A popular alternative to construct a robust solver is the use of multigrid as a preconditioner for a Krylov subspace acceleration method like GMRES. Our contributions to the Fourier analysis for multigrid are two-fold. Firstly we extend the range of situations for which the Fourier analysis can be applied. More precisely, the Fourier analysis is generalized to k-grid cycles and to multigrid as a preconditioner. With a k-grid analysis it is possible to investigate real multigrid effects which cannot be captured by the classical two-grid analysis. Moreover, the k-grid analysis allows for a more detailed investigation of possible coarse grid correction difficulties. Additional valuable insight is obtained by evaluating multigrid as a preconditioner for GMRES. Secondly we extend the range of discretizations and multigrid components for which detailed Fourier analysis results exist. We consider four well-known singularly perturbed model problems to demonstrate the usefulness of the above generalizations: The anisotropic Poisson equation, the rotated anisotropic diffusion equation, the convection diffusion equation with dominant convection, and the driven cavity problem governed by the incompressible Navier Stokes equations. Each of these equations represents a larger class of problems with similar features and complications which are of practical relevance. With the help of the newly developed Fourier analysis methods, a comprehensive study of characteristic difficulties for singular perturbation problems can be performed. Based on the insights from this analysis it is possible to identify remedies resulting in an improved multigrid convergence. The theoretical considerations are validated by numerical test calculations

Fast iterative solvers for convection-diffusion control problems

In this manuscript, we describe effective solvers for the optimal control of stabilized convection-diffusion problems. We employ the local projection stabilization, which we show to give the same matrix system whether the discretize-then-optimize or optimize-then-discretize approach for this problem is used. We then derive two effective preconditioners for this problem, the �first to be used with MINRES and the second to be used with the Bramble-Pasciak Conjugate Gradient method. The key components of both preconditioners are an accurate mass matrix approximation, a good approximation of the Schur complement, and an appropriate multigrid process to enact this latter approximation. We present numerical results to demonstrate that these preconditioners result in convergence in a small number of iterations, which is robust with respect to the mesh size h, and the regularization parameter β, for a range of problems

HIGH ACCURACY MULTISCALE MULTIGRID COMPUTATION FOR PARTIAL DIFFERENTIAL EQUATIONS

Scientific computing and computer simulation play an increasingly important role in scientific investigation and engineering designs, supplementing traditional experiments, such as in automotive crash studies, global climate change, ocean modeling, medical imaging, and nuclear weapons. The numerical simulation is much cheaper than experimentation for these application areas and it can be used as the third way of science discovery beyond the experimental and theoretical analysis. However, the increasing demand of high resolution solutions of the Partial Differential Equations (PDEs) with less computational time has increased the importance for researchers and engineers to come up with efficient and scalable computational techniques that can solve very large-scale problems. In this dissertation, we build an efficient and highly accurate computational framework to solve PDEs using high order discretization schemes and multiscale multigrid method. Since there is no existing explicit sixth order compact finite difference schemes on a single scale grids, we used Gupta and Zhang’s fourth order compact (FOC) schemes on different scale grids combined with Richardson extrapolation schemes to compute the sixth order solutions on coarse grid. Then we developed an operator based interpolation scheme to approximate the sixth order solutions for every find grid point. We tested our method for 1D/2D/3D Poisson and convection-diffusion equations. We developed a multiscale multigrid method to efficiently solve the linear systems arising from FOC discretizations. It is similar to the full multigrid method, but it does not start from the coarsest level. The major advantage of the multiscale multigrid method is that it has an optimal computational cost similar to that of a full multigrid method and can bring us the converged fourth order solutions on two grids with different scales. In order to keep grid independent convergence for the multiscale multigrid method, line relaxation and plane relaxation are used for 2D and 3D convection diffusion equations with high Reynolds number, respectively. In addition, the residual scaling technique is also applied for high Reynolds number problems. To further optimize the multiscale computation procedure, we developed two new methods. The first method is developed to solve the FOC solutions on two grids using standardW-cycle structure. The novelty of this strategy is that we use the coarse level grid that will be generated in the standard geometric multigrid to solve the discretized equations and achieve higher order accuracy solution. It is more efficient and costs less CPU and memory compared with the V-cycle based multiscale multigrid method. The second method is called the multiple coarse grid computation. It is first proposed in superconvergent multigrid method to speed up the convergence. The basic idea of multigrid superconvergent method is to use multiple coarse grids to generate better correction for the fine grid solution than that from the single coarse grid. However, as far as we know, it has never been used to increase the order of solution accuracy for the fine grid. In this dissertation, we use the idea of multiple coarse grid computation to approximate the fourth order solutions on every coarse grid and fine grid. Then we apply the Richardson extrapolation for every fine grid point to get the sixth order solutions. For parallel implementation, we studied the parallelization and vectorization potential of the Gauss-Seidel relaxation by partitioning the grid space with four colors for solving 3D convection-diffusion equations. We used OpenMP to parallelize the loops in relaxation and residual computation. The numerical results show that the parallelized and the sequential implementation have the same convergence rate and the accuracy of the computed solutions

Application of Operator Splitting Methods in Finance

Financial derivatives pricing aims to find the fair value of a financial contract on an underlying asset. Here we consider option pricing in the partial differential equations framework. The contemporary models lead to one-dimensional or multidimensional parabolic problems of the convection-diffusion type and generalizations thereof. An overview of various operator splitting methods is presented for the efficient numerical solution of these problems. Splitting schemes of the Alternating Direction Implicit (ADI) type are discussed for multidimensional problems, e.g. given by stochastic volatility (SV) models. For jump models Implicit-Explicit (IMEX) methods are considered which efficiently treat the nonlocal jump operator. For American options an easy-to-implement operator splitting method is described for the resulting linear complementarity problems. Numerical experiments are presented to illustrate the actual stability and convergence of the splitting schemes. Here European and American put options are considered under four asset price models: the classical Black-Scholes model, the Merton jump-diffusion model, the Heston SV model, and the Bates SV model with jumps