Iterative Methods and High-Order Difference Schemes for 2D Elliptic Problems with Mixed Derivative
|
|
- Moris Gordon
- 6 years ago
- Views:
Transcription
1 Iterative Methods and High-Order Difference Schemes for 2D Elliptic Problems with Mixed Derivative Michel Fournié and Samir Karaa Laboratoire MIP, CNRS UMR 5640, Université Paul Sabatier, 118 route de Narbonne, Toulouse Cedex 4, France Department of Mathematics and Statistics, Sultan Qaboos University, P. O. Box 36, Al-Khod 3, Muscat, Sultanate of Oman 1. Introduction We consider the two-dimensional elliptic equation 2 u x u y 2 + β 2 u x y + c u x + d u = fx, y, x, y Ω, 1 y with Dirichlet boundary conditions on Ω, where Ω is a smooth convex domain in IR 2 consisting of a union of rectangular shapes. The coefficients c, d, and β are assumed constant, and β satisfies the ellipticity condition β 2 < 4. The forcing function fx, y as well as the solution ux, y of the problem are assumed sufficiently smooth in Ω. Efficient numerical solution of Equation 1 plays an important role in many areas especially in computational fluid dynamics. Discretizations of 1 using the traditional 5-point difference schemes: either the central difference scheme CDS or the upwind difference scheme UDS, yield unsatisfactory results. The CDS scheme has a truncation error of order Oh 2 but may produce numerical solutions with nonphysical oscillations for large cell Reynolds numbers. The UDS scheme suppresses these oscillations with large artificial viscosity, and stationary iterative methods for solving the resulting linear system are stable converge for large values of Reynolds numbers. However, it is only of first-order accuracy and requires fine discretizations for satisfactory solution resolution. Fine discretization entails increased computational cost, especially for higher dimensional problems. To obtain satisfactory numerical results with reasonable computational cost, several authors proposed improved finite difference discretization schemes that combine the advantages of the second-order central difference scheme high accuracy and the first-order upwind scheme convergence of stationary iterative methods. Among these is the fourth-order compact scheme FCS due to Gupta, Manohar and Stephenson [2] for convection-diffusion problems β = 0 in two 1
2 2 Michel Fournié and Samir Karaa dimensions. The scheme has a 9-point computational stencil using the eight nearest neighboring points of the reference grid point α 6 α 2 α 5 α 3 α 0 α 1, α 7 α 4 α 8 and has a truncation error of order Oh 4. Numerical experiments in [2] showed that this compact scheme has a good numerical stability for large Reynolds numbers, and SOR iterative methods have been found numerically to converge regardless the magnetude of the Reynolds number. This property is very important for implementing multigtrid method which requires similar discretizations on very course grids [3]. Analytic proofs confirming the convergence of some stationary iterative methods with large Reynolds numbers are given in [4] for convection-diffusion problems with constant coefficients. For problems with variable coefficients, several conditions are formulated in [5, 6], under which the iterative methods converge. Comprehensive studies and useful applications of the fourth-order compact scheme in computational fluid dynamics can be found in at least five recent PhD theses [7]-[11]. The fourth-order compact schemes have also been used in the numerical simulation of incompressible Navier-Stokes equations with good results [1,, 13, 14]. Extension of high-order compact schemes to time-dependent problems are reported in [15, 16, 17]. In this paper, we generalize the 9-point compact scheme to the two-dimensional elliptic problem 1 involving a mixed derivative β 0. After derivation of the high-order compact scheme, we conduct experimental study on the numerical solution of the problem discretized by the new scheme and the traditional secondorder central difference scheme. We study the computed accuracy achieved by each scheme and the performance of some iterative methods for solving linear systems arising from the difference schemes. 2. Derivation of High-Order Scheme The idea behind the derivation of the high-order compact scheme is to operate on the differential equations as an auxiliary relation to obtain finite difference approximations for high-order derivatives in the truncation error. Inclusion of these expressions in a central difference method for Equation 1 increase the order of accuracy, typically to Oh 4 while retaining a compact stencil defined by nodes surrounding a grid point. Introducing a uniform grid with mesh spacing h in both the x- and y-directions, the standard central difference approximation to Equation 1 at grid point i, j is simply δ 2 xu ij δ 2 yu ij + βδ x δ y u ij + cδ x u ij + dδ y u ij τ ij = f ij, 2
3 Iterative Methods and High-Order Difference Schemes 3 where δ x and δx 2 resp. δ y and δy 2 denote the first and second order central difference approximations with respect to x resp. with respect to y. The associated truncation error is given by τ ij = h2 [ 2 c 3 u x 3 + u d 3 y 3 4 u x u y 4 ] ij + βh2 6 [ 4 ] u x 3 y + 4 u x y 3 ij +Oh 4. 3 We now seek second-order approximations to the derivatives in 3. Differentiating Equation 1 with respect to x and y respectively yields 3 u x 3 3 u y 3 = 3 u x y 2 + c 2 u x 2 + d 2 u x y + β 3 u x 2 y f x, 4 = 3 u x 2 y + c 2 u x y + d 2 u y 2 + β 3 u x y 2 f y. 5 Differentiating the equation with respect to x and y we obtain 4 u x 3 y + 4 u x y 3 = β 4 u x 2 y 2 + c 3 u x 2 y + d 3 u x y 2 2 f x y. 6 Notice that all the terms in the right hand sides of 4-6 have compact Oh 2 approximations at node i, j. We have for example 3 [ ] u x 2 y = δxδ 2 y u ij h2 2 5 u ij x 2 y u x 4. y Once again by differentiating twice Equation 1 with respect to x and y respectively, we obtain 4 u x 4 = 4 u x 2 y 2 + u c 3 x 3 + d 3 u x 2 y + β 4 u x 3 y 2 f x 2, 7 4 u y 4 = 4 u x 2 y 2 + c 3 u x y 2 + u d 3 y 3 + β 4 u x y 3 2 f y 2. 8 By adding the two equations, we have 4 u x u y 4 = 2 4 u 3 x 2 y 2 + c u x u x y 2 4 u +β x 3 y + 4 u x y 3 3 u + d x 2 y + 3 u y 3 2 f x 2 2 f y 2, 9 which can be approximated to Oh 2 within the 9-point compact stencil. Substituting Equations 4-9 into Equation 3 and simplifying yields an error term of the form
4 4 Michel Fournié and Samir Karaa [ { τ ij = h2 c 3 u x y 2 + u c 2 x 2 + d 2 u x y + β 3 u x 2 y { } +d 3 u x 2 y + c 2 u x y + u d 2 y 2 + β 3 u x y 2 ] +2 4 u x 2 y 2 c 3 u x y 2 d 3 u x 2 y + βh2 + h2 [ ] β 2 u x 2 y 2 + c 3 u x 2 y + d 3 u x y 2 ij [ c f x d f y + 2 f x f y 2 β 2 f x y ij } ] ij + Oh Now, substituting the central Oh 2 approximations to the derivatives in 10 and inserting in 2 yields the following Oh 4 approximation to c 2 h2 h β2 h 2 h 2 = f ij + h2 δ 2 xu ij δ 2 xδ 2 yu ij βh2 βd 2c d 2 h2 h 2 d δ 2 yu ij + βh2 cβ 2d + c β cdh2 6 δxδ 2 y u ij δ x δ y u ij δ x δyu 2 ij + cδ x u ij + dδ y u ij ]. 11 [ c f x d f y + 2 f x f y 2 β 2 f x y ij The fourth-order compact finite difference formula for the mesh point i, j involves the nearest eight neighboring mesh points: 8 8 α l u l = γ l f l, l=0 l=0
5 Iterative Methods and High-Order Difference Schemes 5 where the coefficients α l, l = 0, 1,..., 8, are obtained from 11 α 0 = 20 + d 2 h 2 2 β 2 + c 2 h 2, α 1 = 1 2 α 2 = 1 2 α 3 = 1 2 α 4 = 1 2 α 5 = ch c 2 h 2 2 β 2 2 hβ d, 8 2 β 2 4 dh + d 2 h 2 2 hcβ, 8 + c 2 h 2 2 β ch + 2 hβ d, 8 + d 2 h hcβ + 4 dh 2 β 2, 2 hβ d 2 ch β 2 6 β + cdh hcβ 2 dh, α 6 = hβ d 2 ch 4 2 β 2 6 β + cdh 2 2 hcβ + 2 dh, α 7 = hβ d + 2 ch β 2 6 β + cdh 2 2 hcβ + 2 dh, α 8 = hβ d + 2 ch 4 2 β 2 6 β + cdh hcβ 2 dh. The coefficents γ l, l = 0,, 8 are γ 0 = 4 h 2, γ 1 = 1 4 h2 2 + ch, γ 2 = 1 4 h2 dh 2, γ 3 = 1 4 h2 ch + 2, γ 4 = 1 4 h2 2 + dh, γ 5 = 1 8 h2 β, γ 6 = 1 8 h2 β, γ 7 = 1 8 h2 β, γ 8 = frac18 h 2 β. 3. Numerical Experiments We conduct numerical experiments to test and compare the computed accuracy achieved by the present fourth-order compact FCS scheme and the central scheme. We also examine the numerical performance of a few iterative methods for solving the resulting linear systems. A similar study has been carried out
6 6 Michel Fournié and Samir Karaa in [18] for convection-diffusion problems. We consider the following two test problems c = Re, Test Problem 1. d = Re, ux, y = sinπx + sinπy + sinπx sinπy. Test Problem 2. c = Re, d = Re, ux, y = [e Re1 x + e Rey 2]/[e Re 1]. The problems are solved on a unit square Ω = 0, 1 0, 1 using a uniform meshsize h in the x- and y-directions. The right-hand side fx, y and the Dirichlet boundary conditions are prescribed to satisfy the given exact solution ux, y. Re is referred to as the Reynolds number. In contrast to Problem 1, the exact solution of Problem 2 is related to the Reynolds number. Even for moderate Re Re 100 the solution is nearly identically zero in Ω except at the boundary layers of thickness O1/Re near x = 0 and y = 1. This example is used to test how well the discretization schemes resolve the thin boundary layers. In our experiments, Re is allowed to vary between 1 and 10 6 and β between -2 and Computed accuracy We first test the computed accuracy that was achieved by each discretization scheme. We select GMRES preconditioned by ILUTτ,s to solve the resulting linear systems in all test runs. The parameters τ and s are selected for each test so that full convergence is reached. The computed error is the maximum absolute error over the discretized grid points, and the accuracy order is obtained by comparing the errors after refining the meshsize, i.e., ln 2 e h=1/64 /e h=1/8. All computer programs were coded in Fortran programming language and were run on a SunBlade100 machine. Test results for Problems 1 and 2, with differents values of h, Re, and β are listed in Tables 1 and 2, respectively. The numerical results show that the magnitude of Re affects the computed accuracy of the discretization schemes inversely, and that β has a little affect. For both test problems with Re < 10 3, the fourth-order compact scheme yields much better solution than the central difference scheme, and both schemes seem to maintain their typical orders of accuracy. For Problem 1 with large Re and without boundary layers, the central difference scheme produces accurate solution and succeeds to maintain its order of accuracy. The fourth-order scheme produces also accurate solution, but fails to preserve its high accuracy especially when Re 10 5, and becomes typically
7 Iterative Methods and High-Order Difference Schemes 7 Table 1. Maximum errors and order of accuracy of different discretization schemes Problem 1. Central Difference Scheme Fourth-Order Scheme Re h 1/64 h 1/8 order h 1/64 h 1/8 order β = β = β = of order 2. In [19], Berger et all. showed that in the 1D case, the approximation error of the fourth-order compact scheme is OReh 4 if Reh < 1 and Oh 2 if Reh > 1. Our numerical results show a similarity. We observe that the approximation error in the 2D case is Oh 2 when Reh 2 > 1. For Problem 2 with thin boundary layers, the central difference scheme behaved badly due to the large numerical oscillations. The fourth-order scheme yields satisfactory solution only for Re For Re 10 3, the error produced by the FCS scheme is small but not satisfactory. A careful examination of the computed solution showed that, in contrast to the CDS case, numerical oscillations are not presented in the solution. For large Re, no scheme yielded acceptable results. This boundary layer problem is common to all uniform discretization schemes. The accuracy degradation is mainly caused by the insufficient number of grid points inside the boundary layers. Special treatments, such
8 8 Michel Fournié and Samir Karaa Table 2. Maximum errors and order of accuracy of different discretization schemes Problem 2. Central Difference Scheme Fourth-Order Scheme Re h 1/64 h 1/8 order h 1/64 h 1/8 order β = β = β = as local adaptive mesh refinement techniques or graded mesh techniques, are then needed to resolve the thin boundary layers. The use of nonuniform grids for convection-diffusion problems with boundary layers have been investigated by Gupta, Manohar and Stephenson [2], and Spotz and Carey [20]. It has been found that the fourth-order accuracy can be recovered by employing the graded mesh techniques. A numerical comparison with the central and upwind schemes in such situations is given in [21]. 3.2 Performance of iterative methods We now examine the numerical performance of a few iterative methods for solving the sparse linear systems arising from the discretization schemes. The
9 Iterative Methods and High-Order Difference Schemes 9 first iterative method to be tested is the line Gauss-Seidel method along the x- axis. The grid points with the same y-index are solved simultaneously. One line Gauss-Seidel iteration involves the solution of N tridiagonal linear subsystems of size N each, where N is the number of interior grid points in each direction. Although stationary iterative methods are seldom used as stand alone solvers, they are usually employed as basic components in building modern iterative methods such as the multigrid methods [22] and the Krylov subspace methods [23]. In such cases, basic iterative methods may be used as smoothers in multigrid methods or as preconditioners in Krylov subspace methods. Understanding the behaviors of these iterative methods for solving sparse linear systems arising from the fourth-order compact discretization is practically interesting. Table 3. Number of line Gauss-Seidel iterations for solving the linear systems discretized by the central and compact difference schemes with h = 1/64 Problem 1. β = 0.2 β = 1 β = 1.8 β = 1 β = 1.8 Re CDS FCS CDS FCS CDS FCS CDS FCS CDS FCS indicates divergence. In Table 3 we list the number of Gauss-Seidel iterations for solving the linear systems arising from Problem 1 discretized by the central and compact difference schemes with h = 1/64 and with different values of Re and β. The experimental data show that Gauss-Seidel with the central scheme diverges for Re 10 3 and for all the values of β listed in the table. Several numerical tests reveal that divergence occurs for all values of β in the interval 2, 2, and similar results has been obtained in the case of Problem 2. Gauss-Seidel method with the fourth-order compact scheme is found to converge in all cases, regardless of the magnetude of Re, but its performance deteriorates as Re increases, as well as for small Re. The Table also reveals that the number of Gauss-Seidel iterations for both schemes are close, which means that the amounts of work in both cases do not differ much. The number of Gauss-Seidel iterations as a function of Re with different values of β is shown in Figure 1 for the two test problems. In summary, the numerical experiments can be explained as that the Gauss-Seidel method is more stable with the compact scheme, slightly affected by the values of β when Re 10 2, but heavily affected
10 10 Michel Fournié and Samir Karaa by the magnitude of Re. Its performance is usually good when Re is between 10 2 and We notice that in the non-elliptic case for example when β = 3, the Gauss-Seidel method does not converge for all Reynolds numbers β =1.8 β =1 β = 1 β = 1.8 Problem β =1.8 β =1 β = 1 β = 1.8 Problem Number of GS iterations Number of GS iterations Re Re Figure 1. Number of line Gauss-Seidel iterations as a function of Re for different values of β, for solving the test problems with the compact scheme, h = 1/64. We next solve the linear systems arising from the discretization schemes using a standard multigrid method. The multigrid method has been shown to be very effective for solving discretized elliptic problems [22]. In the present computations, we use a standard multigrid V-cycle algorithm, with the point Gauss-Seidel method as a smoother. We perform one relaxation on each level before projection and after interpolation. Standard full-weighting and bilinear interpolation are employed as the intergrid transfer operators. Standard coarsening technique is used the meshsize of the coarse grid doubles that of the fine grid and all possible levels are exployed, that is the coarsest level has a meshsize h = 1/2. For a description of the multigrid method with different cycling algorithms, we refer to the books [22] and [24]. We choose the finest grid with h = 1/64, the initial guess is ux, y = 0. The iterations are terminated when the 2-norm of the residual is reduced by 10 10, i.e., r n 2 / r 1 2 < Various studies of the multigrid method for solving convection-diffusion equations discretized by the fourth-order compact scheme can be found in [10, 3, 25]. Convergence results for the multigrid method with different Re and β are presented in Tables 4 and 5 for Problems 1 and 2, respectively. The data show that the multigrid method is more stable with the fourth-order compact scheme.
11 Iterative Methods and High-Order Difference Schemes 11 It always converges when Re 10 3 for all values of β with a few exceptions at Re = 10 2, while it diverges with the central scheme when Re 10 2 in all cases. The use of the red-black Gauss-Seidel as a smoother did not improve the results. We believe that the smoother is not the main factor for divergence with the compact scheme, and that suitable intergrid transfer operators with appropriate scaling of the residual, see for instance [26] and [10], can make the multigrid method converge in all cases with the compact scheme. Table 4. Number of Multigrid iterations for solving Problem 1, h = 1/64. β = 0.2 β = 1 β = 1.8 β = 1 β = 1.8 Re CDS FCS CDS FCS CDS FCS CDS FCS CDS FCS indicates divergence. Table 5. Number of Multigrid iterations for solving Problem 2, h = 1/64. β = 0.2 β = 1 β = 1.8 β = 1 β = 1.8 Re CDS FCS CDS FCS CDS FCS CDS FCS CDS FCS indicates divergence. We finally solve the linear systems arising from the discretization schemes using GMRES [27] preconditioned by the incomplete LU preconditioner ILUTτ,s. In full GMRES implementation the storage requirements grow quadratically with the number of iterations. Hence, in practice it is often necessary to use a restarted version, GMRESm, where m indicates the selected dimension of the Krylov subspace. In our case GMRES is restarted after every 20 iterations. The ILUT preconditioner uses a dual truncation dropping strategy developed in [28]. The amount of fill-ins is controlled by two parameters, a threshold drop tolerance τ and a fill-in number s. All entries with magnetude less than s multiplied by the average value of absolute values of the current row are dropped. In addition,
12 Michel Fournié and Samir Karaa Table 6. Number of preconditioned GMRES20 iterations and the amount of fill-ins in ILUT for the different schemes with h = 1/64 and β = 1. Problem 1 Problem 2 CDS FCS CDS FCS Re s Iter. s Iter. s Iter. s Iter Table 7. Number of preconditioned GMRES20 iterations and the amount of fill-ins in ILUT for the different schemes with h = 1/8 and β = 1. Problem 1 Problem 2 CDS FCS CDS FCS Re s Iter. s Iter. s Iter. s Iter indicates divergence. only the largest s entries in each row of the L and U factors are retained. When testing with this preconditioner, it is convenient to set τ to a fixed value, and to vary s to produce a preconditioner with desired amount of fill-ins [29]. In this way, the amount of storage can be predicted in advance and the trade-off between memory cost and computing time can be balanced more easily. In our computations, the dropping tolerenace τ = 10 4 is fixed, and the number of fill-ins s is initially chosen equal to 5. If convergence is not reached within 200 iterations, then s is increased by 5 each time untill we have convergence. The number of preconditioned GMRES iterations and the corresponding parameter s of fill-ins in ILUT for solving the linear systems with h = 1/64 and h = 1/8 are listed in Tables 6 and 7, respectively. We notice that for Problem
13 Iterative Methods and High-Order Difference Schemes 13 2 with Re = 10 5 and h = 1/8, we decreased the dropping tolerance τ to 10 5 in order to achieve convergence with the central difference scheme. The tables indicate that with the fourth-order scheme a small storage space is needed; s = 5 was always enough except for Problem 2 when Re However, with the central scheme, s exceeds 50 in all cases with Re We also remark that when we have the same value for s, GMRES has the least iteration counts with the fourth-order scheme. Since the setup of the ILUT preconditioner requires a certain amount of preprocessing costs depending in part on τ and s, we compare in Figure 2 the total CPU time in seconds, including the setup of the ILUT preconditioner, for solving Problem 2 with h = 1/8. The figure shows a significant growth of the CPU time in the case of the CDS scheme with large Re, while the CPU time remained less than 9 seconds through the entire simulation with the FCS scheme. We can say that preconditioned GMRES is more robust with the FCS scheme, and more effective in terms of computational cost and storage requirement PGMRES with CDS PGMRES with FCS Total CPU time in seconds Figure 2. Comparison of total CPU time in seconds delievered by preconditioned GMRES for solving Problem 2 with the central and the fourth-order schemes h = 1/8. Re 4. Concluding remarks
14 14 Michel Fournié and Samir Karaa We derived a fourth-order compact finite difference scheme for a 2D elliptic problem with a mixed derivative and constant coefficients. Numerical experiments are conducted to test its high accuracy and compare it with the standard central difference scheme. The numerical performances of a few iterative methods for solving the sparse linear systems arising from the discretization schemes are examined. The Gauss-Seidel and multigrid methods showed good stablility with the compact scheme and preconditioned GMRES was found more effective in terms of computational cost and storage requirement. We finally point out that, in the case where the coefficients in the governing equation 1 are not constant, we can derive a compact fourth-order scheme using techniques different from those presented in Section 2. Acknowledgments : S. Karaa s research was supported by Sultan Qaboos University under Grant IG/SCI/DOMAS/05/10 References 1. M. M. Gupta, R. P. Manohar and J. W. Stephenson, A fourth-order, cost effective and stable finite difference scheme for the convection-diffusion equation, In Numerical Properties and Methodologies in Heat Transfer, Washington, DC, 1983, Hemisphere Publishing Corp, M. M. Gupta, R. P. Manohar and J. W. Stephenson, A single cell high-order scheme for the convection-diffusion equation with variable coefficients, Int. J. Numer. Methods Fluids, , M. M. Gupta, J. Kouatchou and J. Zhang, A compact multigrid solver for convectiondiffusion equations, J. Comput. Phys., , S. Karaa and J. Zhang, Analysis of stationary iterative methods for the discrete convectiondiffusion equation with a 9-point compact scheme, J. Comput. Appl. Math., , No. 2, S. Karaa and J. Zhang, A note on convergence of line iterative methods for a nine-point matrix, Appl. Math. Lett., , No. 4, S. Karaa and J. Zhang, Convergence and performance of iterative methods for solving variable coefficient convection-diffusion equation with a fourth-order compact difference scheme, Comput. Math. Appl., , No. 3-4, J. Kouatchou, High-Order Multigrid Techniques for Partial Differential Equations, Ph.D. thesis, The George Washington University, Washington, DC, M. Li, Numerical Solution for the Incompressible Navier-Stokes Equations, Ph.D. thesis, Simon Frasher University, British Columbia, Canada, W. F. Spotz, High-Order Compact Finite Difference Schemes for Computational Mechanics, Ph.D. thesis, University of Texas at Austin, Austin, TX, J. Zhang, Multigrid Acceleration Techniques and Applications to the Numerical Solution of Partial Differential Equations, Ph.D. thesis, The George Washington University, Washington, DC, M. Fournié, Construction et analyse de schémas compacts d ordre élevé pour des problémes fortement convectifs. Applications à la simulation de semi-conducteurs, Ph.D. thesis, Paul Sabatier University, Toulouse, France, 1999.
15 Iterative Methods and High-Order Difference Schemes 15. S. C. R. Dennis and J. D. Hudson, Compact h 4 finite difference approximations to the operators of Navier-Stokes type, J. Comput. Phys., , M. Li, T. Tang and B. Fornberg, A compact fourth-order finite difference scheme for the steady incompressible Navier-Stokes equations, Int. J. Numer. Methods Fluids, , W. F. Spotz and C. F. Carey, High-order compact scheme for the steady stream-function vorticity equations, Int. J. Numer. Methods Eng., , W. F. Spotz and C. F. Carey, Extension of high-order compact schemes to time-dependent problems, Numer. Methods Partial Differential Equations, , No. 6, J. C. Kalita, D. C. Dalal and A. K. Dass, A class of higher order compact schemes for the unsteady two-dimensional convection-diffusion equation with variable convection coefficients, Internat. J. Numer. Methods Fluids, , No., S. Karaa and J. Zhang, High order ADI method for solving unsteady convection-diffusion problems, J. Comput. Phys., , No. 1, J. Zhang, Preconditioned iterative methods and finite difference schemes for convectiondiffusion, Appl. Math. Comput., , A. E. Berger, J. M. Solomon, M. Ciment and B. C. Weinberg, Generalized OCI schemes for boundary layer problems, Math. Comput., , W. F. Spotz and C. F. Carey, Formulation and experiments with high-order compact schemes for nonuniform grids, Int. J. Numer. Methods for Heat and Fluid Flow, , L. Ge and J. Zhang, Accuracy, robustness, and efficiency comparison in iterative computation of convection diffusion equation with boundary layers, Numer. Methods Partial Differential Equations, , P. Wesseling, An Introduction to Multigrid Methods, Wiley, Chichester, England, Y. Saad, Iterative Methods for Sparse Linear Systems, PWS Publishing, New York, NY, W. L. Briggs. A Multigrid Tutorial, SIAM, Philadelphia, PA, A. L. Pardhanani and W. F. Spotz and G. F. Carey, A stable multigrid strategy for convection-diffusion using high order compact discretization, Electron. Trans. Numer. Anal., , P. M. de Zeeuw, Matrix-dependent prolongations and restrictions in a blackbox multigrid solver, J. Comput. Appl. Math., , No. 1, Y. Saad and M. H. Schultz, GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Statist. Comput., , Y. Saad, ILUT: a dual threshold incomplete LU preconditioner, Numer. Linear Algebra Appli., , No. 4, A. Chapman, Y. Saad and L. Wigton, High-order ILU preconditioners for CFD problems, Internat. J. Numer. Methods Fluids, , No. 6,
High-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
More informationAMS Mathematics Subject Classification : 65F10,65F50. Key words and phrases: ILUS factorization, preconditioning, Schur complement, 1.
J. Appl. Math. & Computing Vol. 15(2004), No. 1, pp. 299-312 BILUS: A BLOCK VERSION OF ILUS FACTORIZATION DAVOD KHOJASTEH SALKUYEH AND FAEZEH TOUTOUNIAN Abstract. ILUS factorization has many desirable
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 39, pp. 3-45, 01. Copyright 01,. ISSN 1068-9613. ETNA A COMBINED FOURTH-ORDER COMPACT SCHEME WITH AN ACCELERATED MULTIGRID METHOD FOR THE ENERGY EQUATION
More informationMultigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme Accelerated by Krylov Subspace
201, TextRoad Publication ISSN: 2090-27 Journal of Applied Environmental and Biological Sciences www.textroad.com Multigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationMULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW
MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW VAN EMDEN HENSON CENTER FOR APPLIED SCIENTIFIC COMPUTING LAWRENCE LIVERMORE NATIONAL LABORATORY Abstract Since their early application to elliptic
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccs Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary c 2008,2010
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary
More informationANALYSIS AND COMPARISON OF GEOMETRIC AND ALGEBRAIC MULTIGRID FOR CONVECTION-DIFFUSION EQUATIONS
ANALYSIS AND COMPARISON OF GEOMETRIC AND ALGEBRAIC MULTIGRID FOR CONVECTION-DIFFUSION EQUATIONS CHIN-TIEN WU AND HOWARD C. ELMAN Abstract. The discrete convection-diffusion equations obtained from streamline
More informationIncomplete LU Preconditioning and Error Compensation Strategies for Sparse Matrices
Incomplete LU Preconditioning and Error Compensation Strategies for Sparse Matrices Eun-Joo Lee Department of Computer Science, East Stroudsburg University of Pennsylvania, 327 Science and Technology Center,
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationChapter 5. Methods for Solving Elliptic Equations
Chapter 5. Methods for Solving Elliptic Equations References: Tannehill et al Section 4.3. Fulton et al (1986 MWR). Recommended reading: Chapter 7, Numerical Methods for Engineering Application. J. H.
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationA New Multilevel Smoothing Method for Wavelet-Based Algebraic Multigrid Poisson Problem Solver
Journal of Microwaves, Optoelectronics and Electromagnetic Applications, Vol. 10, No.2, December 2011 379 A New Multilevel Smoothing Method for Wavelet-Based Algebraic Multigrid Poisson Problem Solver
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationTHE EFFECT OF MULTIGRID PARAMETERS IN A 3D HEAT DIFFUSION EQUATION
Int. J. of Applied Mechanics and Engineering, 2018, vol.23, No.1, pp.213-221 DOI: 10.1515/ijame-2018-0012 Brief note THE EFFECT OF MULTIGRID PARAMETERS IN A 3D HEAT DIFFUSION EQUATION F. DE OLIVEIRA *
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationMathematics and Computer Science
Technical Report TR-2007-002 Block preconditioning for saddle point systems with indefinite (1,1) block by Michele Benzi, Jia Liu Mathematics and Computer Science EMORY UNIVERSITY International Journal
More informationThe flexible incomplete LU preconditioner for large nonsymmetric linear systems. Takatoshi Nakamura Takashi Nodera
Research Report KSTS/RR-15/006 The flexible incomplete LU preconditioner for large nonsymmetric linear systems by Takatoshi Nakamura Takashi Nodera Takatoshi Nakamura School of Fundamental Science and
More informationA greedy strategy for coarse-grid selection
A greedy strategy for coarse-grid selection S. MacLachlan Yousef Saad August 3, 2006 Abstract Efficient solution of the very large linear systems that arise in numerical modelling of real-world applications
More information2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More informationAggregation-based algebraic multigrid
Aggregation-based algebraic multigrid from theory to fast solvers Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire CEMRACS, Marseille, July 18, 2012 Supported by the Belgian FNRS
More informationPreconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners
Preconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, Georgia, USA
More informationK.S. Kang. The multigrid method for an elliptic problem on a rectangular domain with an internal conductiong structure and an inner empty space
K.S. Kang The multigrid method for an elliptic problem on a rectangular domain with an internal conductiong structure and an inner empty space IPP 5/128 September, 2011 The multigrid method for an elliptic
More informationITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS
ITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS YOUSEF SAAD University of Minnesota PWS PUBLISHING COMPANY I(T)P An International Thomson Publishing Company BOSTON ALBANY BONN CINCINNATI DETROIT LONDON MADRID
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationSplitting Iteration Methods for Positive Definite Linear Systems
Splitting Iteration Methods for Positive Definite Linear Systems Zhong-Zhi Bai a State Key Lab. of Sci./Engrg. Computing Inst. of Comput. Math. & Sci./Engrg. Computing Academy of Mathematics and System
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationNewton-Krylov-Schwarz Method for a Spherical Shallow Water Model
Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Chao Yang 1 and Xiao-Chuan Cai 2 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, yang@mail.rdcps.ac.cn
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationMATHEMATICAL RELATIONSHIP BETWEEN GRID AND LOW PECLET NUMBERS FOR THE SOLUTION OF CONVECTION-DIFFUSION EQUATION
2006-2018 Asian Research Publishing Network (ARPN) All rights reserved wwwarpnjournalscom MATHEMATICAL RELATIONSHIP BETWEEN GRID AND LOW PECLET NUMBERS FOR THE SOLUTION OF CONVECTION-DIFFUSION EQUATION
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationAlgebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov
More informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationA new 9-point sixth-order accurate compact finite difference method for the Helmholtz equation
A new 9-point sixth-order accurate compact finite difference method for the Helmholtz equation Majid Nabavi, M. H. Kamran Siddiqui, Javad Dargahi Department of Mechanical and Industrial Engineering, Concordia
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationThe Effect of the Schedule on the CPU Time for 2D Poisson Equation
Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014. The Effect of the Schedule on the CPU Time for 2D Poisson Equation Fabiane de Oliveira, State University of Ponta Grossa, Department of Mathematics and
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationImplementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows
Short Training Program Report Implementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows Julian Koellermeier RWTH Aachen University Supervisor: Advisor: Prof. Thierry Magin von
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners for
More informationAn efficient multigrid solver based on aggregation
An efficient multigrid solver based on aggregation Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Graz, July 4, 2012 Co-worker: Artem Napov Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationFinding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems
Finding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems AMSC 663-664 Final Report Minghao Wu AMSC Program mwu@math.umd.edu Dr. Howard Elman Department of Computer
More informationA Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows
Journal of Scientific Computing, Vol. 16, No. 1, March 2001 ( 2001) A Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows Ming Li 1 and Tao Tang 2 Received January 23,
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More informationA fast method for the solution of the Helmholtz equation
A fast method for the solution of the Helmholtz equation Eldad Haber and Scott MacLachlan September 15, 2010 Abstract In this paper, we consider the numerical solution of the Helmholtz equation, arising
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
More informationBLOCK ILU PRECONDITIONED ITERATIVE METHODS FOR REDUCED LINEAR SYSTEMS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 12, Number 3, Fall 2004 BLOCK ILU PRECONDITIONED ITERATIVE METHODS FOR REDUCED LINEAR SYSTEMS N GUESSOUS AND O SOUHAR ABSTRACT This paper deals with the iterative
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationGeneralized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,
Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More information1. Introduction. In this work we consider the solution of finite-dimensional constrained optimization problems of the form
MULTILEVEL ALGORITHMS FOR LARGE-SCALE INTERIOR POINT METHODS MICHELE BENZI, ELDAD HABER, AND LAUREN TARALLI Abstract. We develop and compare multilevel algorithms for solving constrained nonlinear variational
More informationA Numerical Study of Some Parallel Algebraic Preconditioners
A Numerical Study of Some Parallel Algebraic Preconditioners Xing Cai Simula Research Laboratory & University of Oslo PO Box 1080, Blindern, 0316 Oslo, Norway xingca@simulano Masha Sosonkina University
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationAvailable online: 19 Oct To link to this article:
This article was downloaded by: [Academy of Mathematics and System Sciences] On: 11 April 01, At: 00:11 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 107954
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More information9.1 Preconditioned Krylov Subspace Methods
Chapter 9 PRECONDITIONING 9.1 Preconditioned Krylov Subspace Methods 9.2 Preconditioned Conjugate Gradient 9.3 Preconditioned Generalized Minimal Residual 9.4 Relaxation Method Preconditioners 9.5 Incomplete
More informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More information