In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.1...In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.134:108344,2022).The convergence analysis of the proposed method is established,which can improve the existing results.Numerical examples show that the proposed method is efficient with the two-step technique.展开更多
In this paper,we present a modulus-based multisplitting iteration method based on multisplitting of the system matrix for a class of weakly nonlinear complementarity problem.And we prove the convergence of the method ...In this paper,we present a modulus-based multisplitting iteration method based on multisplitting of the system matrix for a class of weakly nonlinear complementarity problem.And we prove the convergence of the method when the system matrix is an H_(+)-matrix.Finally,we give two numerical examples.展开更多
In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent im...In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent implicit fixed-point equation, then introduces a smoothing function to obtain its approximation solutions. The convergence analysis of the algorithm was given, and the efficiency of the algorithms was verified by numerical experiments.展开更多
To develop an efficient and robust aerodynamic analysis method for numerical optimization designs of wing and complex configuration, a combination of matrix preconditioning and multigrid method is presented and invest...To develop an efficient and robust aerodynamic analysis method for numerical optimization designs of wing and complex configuration, a combination of matrix preconditioning and multigrid method is presented and investigated. The time derivatives of three-dimensional Navier-Stokes equations are preconditioned by Choi-Merkle preconditioning matrix that is originally designed for two-dimensional low Mach number viscous flows. An extension to three-dimensional viscous flow is implemented, and a method improving the convergence for transonic flow is proposed. The space discretizaition is performed by employing a finite-volume cell-centered scheme and using a central difference. The time marching is based on an explicit Rtmge-Kutta scheme proposed by Jameson. An efficient FAS multigrid method is used to accelerate the convergence to steady-state solutions. Viscous flows over ONERA M6 wing and M100 wing are numerically simulated with Mach numbers ranging from 0.010 to 0.839. The inviscid flow over the DLR-F4 wing-body configuration is also calculated to preliminarily examine the performance of the presented method for complex configuration. The computed results are compared with the experimental data and good agreement is achieved. It is shown that the presented method is efficient and robust for both compressible and incompressible flows and is very attractive for aerodynamic optimization designs of wing and complex configuration.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
In this paper we study the theoretical properties of multigrid algorithm for discretization of the Poisson equation in 2D using a mortar element method under the assumption that the triangulations on every subdomain a...In this paper we study the theoretical properties of multigrid algorithm for discretization of the Poisson equation in 2D using a mortar element method under the assumption that the triangulations on every subdomain are uniform.We prove the convergence of the W-cycle with a sufficiently large number of smoothing steps.The variable V-cycle multigrid preconditioner are also available.展开更多
Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomp...Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomposition (ILU) method and a combined multigrid ILU method are used to solve the linear system.Systematic computations using the ILU method have shown that the CPU time can be reduced to 30% to 40% of that using an incomplete Gaussian elimination method. In the proposed multigrid ILU method an averaged restriction and a piecewise constant prolongation are used.The construction of the coefficient matrix at coarse levels is based on geometrical considerations.It turns out that the condition of the relative consistency is fulfilled.Comparison computations have shown that nearly the same results were obtained.However,due to additional CPU time needed for the execution of the matrix vector products in the restriction and the prolongation proceses of the multigrid method,a further reduction of the total CPU time could not be reailized.展开更多
In this paper, an optimal V-cycle multigrid method for some conforming and nonconforming plate elements are constructed. A new method dealing with nonnested multigrid methods is presented.
A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-bal...A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-balanced finite volume discretization of the steady-state SWEs is solved by the Newton method as the outer iteration and a geometric multigrid method with the block symmetric Gauss-Seidel smoother as the inner iteration. The proposed Newton multigrid method makes use of the local residual to regularize the Jacobian matrix of the Newton iteration, and can handle the steady- state problem with wet/dry transition. Several numerical experiments are conducted to demonstrate the efficiency, robustness, and well-balanced property of the proposed method. The relation between the convergence behavior of the Newton multigrid method and the distribution of the eigenvalues of the iteration matrix is detailedly discussed.展开更多
In this paper, we obtained the numerical solutions of the modified regularized long-wave (MRLW) equation, by using the multigrid method and finite difference method. The solitary wave motion, interaction of two and th...In this paper, we obtained the numerical solutions of the modified regularized long-wave (MRLW) equation, by using the multigrid method and finite difference method. The solitary wave motion, interaction of two and three solitary waves, and development of the Maxwellian initial condition into solitary waves are studied using the proposed method. The numerical solutions are compared with the known analytical solutions. Usingerror norms and conservative properties of mass, momentum and energy, accuracy and efficiency of the mentioned method will be established through comparison with other techniques.展开更多
This paper presents a numerical method for PDE-constrained optimization problems. These problems arise in many fields of science and engineering including those dealing with real applications. The physical problem is ...This paper presents a numerical method for PDE-constrained optimization problems. These problems arise in many fields of science and engineering including those dealing with real applications. The physical problem is modeled by partial differential equations (PDEs) and involve optimization of some quantity. The PDEs are in most cases nonlinear and solved using numerical methods. Since such numerical solutions are being used routinely, the recent trend has been to develop numerical methods and algorithms so that the optimization problems can be solved numerically as well using the same PDE-solver. We present here one such numerical method which is based on simultaneous pseudo-time stepping. The efficiency of the method is increased with the help of a multigrid strategy. Application example is included for an aerodynamic shape optimization problem.展开更多
The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consis...The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.展开更多
In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NAC...In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.展开更多
In this paper, a V-cycle multigrid method is presented for a Hermite rectangular element. By defining proper mesh-dependent inner product and transfer operator, we obtain its convergence property and the uniform conve...In this paper, a V-cycle multigrid method is presented for a Hermite rectangular element. By defining proper mesh-dependent inner product and transfer operator, we obtain its convergence property and the uniform convergence rate independent of mesh size and level are established.展开更多
This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same ...This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same as that of the SOR method and the CPU time is much less than that of the latter one.展开更多
In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids....In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.展开更多
Numerical solutions of the modified equal width wave equation are obtained by using the multigrid method and finite difference method. The motion of a single solitary wave, interaction of two solitary waves and develo...Numerical solutions of the modified equal width wave equation are obtained by using the multigrid method and finite difference method. The motion of a single solitary wave, interaction of two solitary waves and development of the Maxwellian initial condition into solitary waves are studied using the proposed method. The numerical solutions are compared with the known analytical solutions. Using error norms and conservative properties of mass, momentum and energy, accuracy and efficiency of the mentioned method will be established through comparison with other methods.展开更多
This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were di...This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were discretized by Finite Volume Method(FVM)with second-order approximation scheme and deferred correction.Moreover,the algebraic method and the differential method were used to generate the non-orthogonal structured curvilinear grids.Furthermore,the influence of some parameters of geometric multigrid method,as well as lexicographical Gauss–Seidel(Lex-GS),η-line Gauss–Seidel(η-line-GS),Modified Strongly Implicit(MSI)and modified incomplete LU decomposition(MILU)solvers on the Central Processing Unit(CPU)time was investigated.Therefore,several parameters of multigrid method and solvers were tested for the problem,with the use of nonorthogonal structured curvilinear grids and multigrid method,resulting in an algorithm with the combination that achieved the best results and CPU time.The geometric multigrid method with Full Approximation Scheme(FAS),V-cycle and standard coarsening ratio for this problem were utilized.This article shows how to calculate the coordinates transformation metrics in the coarser grids.Results show that the MSI and MILU solvers are the most efficient.Moreover,theMSI solver is faster thanMILU for both grids generators;and the solutions are more accurate for the Burgers problem with grids generated using elliptic equations.展开更多
Cascadic multigrid technique for mortar Wilson finite element method ofhomogeneous boundary value planar linear elasticity is described and analyzed. Firstthe mortar Wilson finite element method for planar linear elas...Cascadic multigrid technique for mortar Wilson finite element method ofhomogeneous boundary value planar linear elasticity is described and analyzed. Firstthe mortar Wilson finite element method for planar linear elasticity will be analyzed,and the error estimate under L2 and H1 norm is optimal. Then a cascadic multigridmethod for the mortar finite element discrete problem is described. Suitable grid trans-fer operator and smoother are developed which lead to an optimal cascadic multigridmethod. Finally, the computational results are presented.展开更多
In this paper a cascadic multigrid algorithm for the mortar finite element approximation of the semilinear elliptic problem is proposed,and corresponding theorems aregiven,which display the error estimate and the comp...In this paper a cascadic multigrid algorithm for the mortar finite element approximation of the semilinear elliptic problem is proposed,and corresponding theorems aregiven,which display the error estimate and the computational complexity of the method.展开更多
基金supported by the Scientific Computing Research Innovation Team of Guangdong Province(no.2021KCXTD052)the Science and Technology Development Fund,Macao SAR(no.0096/2022/A,0151/2022/A)+3 种基金University of Macao(no.MYRG2020-00035-FST,MYRG2022-00076-FST)the Guangdong Key Construction Discipline Research Capacity Enhancement Project(no.2022ZDJS049)Technology Planning Project of Shaoguan(no.210716094530390)the ScienceFoundation of Shaoguan University(no.SZ2020KJ01).
文摘In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.134:108344,2022).The convergence analysis of the proposed method is established,which can improve the existing results.Numerical examples show that the proposed method is efficient with the two-step technique.
基金This work was supported by the National Natural Science Foundation of China(Grant No.11771275)the Science and Technology Program of Shandong Universities(No.J16LI04).
文摘In this paper,we present a modulus-based multisplitting iteration method based on multisplitting of the system matrix for a class of weakly nonlinear complementarity problem.And we prove the convergence of the method when the system matrix is an H_(+)-matrix.Finally,we give two numerical examples.
文摘In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent implicit fixed-point equation, then introduces a smoothing function to obtain its approximation solutions. The convergence analysis of the algorithm was given, and the efficiency of the algorithms was verified by numerical experiments.
文摘To develop an efficient and robust aerodynamic analysis method for numerical optimization designs of wing and complex configuration, a combination of matrix preconditioning and multigrid method is presented and investigated. The time derivatives of three-dimensional Navier-Stokes equations are preconditioned by Choi-Merkle preconditioning matrix that is originally designed for two-dimensional low Mach number viscous flows. An extension to three-dimensional viscous flow is implemented, and a method improving the convergence for transonic flow is proposed. The space discretizaition is performed by employing a finite-volume cell-centered scheme and using a central difference. The time marching is based on an explicit Rtmge-Kutta scheme proposed by Jameson. An efficient FAS multigrid method is used to accelerate the convergence to steady-state solutions. Viscous flows over ONERA M6 wing and M100 wing are numerically simulated with Mach numbers ranging from 0.010 to 0.839. The inviscid flow over the DLR-F4 wing-body configuration is also calculated to preliminarily examine the performance of the presented method for complex configuration. The computed results are compared with the experimental data and good agreement is achieved. It is shown that the presented method is efficient and robust for both compressible and incompressible flows and is very attractive for aerodynamic optimization designs of wing and complex configuration.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金This research was supported by the National Natural Science Foundation of China under grant 10071015
文摘In this paper we study the theoretical properties of multigrid algorithm for discretization of the Poisson equation in 2D using a mortar element method under the assumption that the triangulations on every subdomain are uniform.We prove the convergence of the W-cycle with a sufficiently large number of smoothing steps.The variable V-cycle multigrid preconditioner are also available.
文摘Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomposition (ILU) method and a combined multigrid ILU method are used to solve the linear system.Systematic computations using the ILU method have shown that the CPU time can be reduced to 30% to 40% of that using an incomplete Gaussian elimination method. In the proposed multigrid ILU method an averaged restriction and a piecewise constant prolongation are used.The construction of the coefficient matrix at coarse levels is based on geometrical considerations.It turns out that the condition of the relative consistency is fulfilled.Comparison computations have shown that nearly the same results were obtained.However,due to additional CPU time needed for the execution of the matrix vector products in the restriction and the prolongation proceses of the multigrid method,a further reduction of the total CPU time could not be reailized.
基金The rescarch was supported by the Doctoral Point Foundation of chinese Universities and NSF
文摘In this paper, an optimal V-cycle multigrid method for some conforming and nonconforming plate elements are constructed. A new method dealing with nonnested multigrid methods is presented.
基金Project supported by the National Natural Science Foundation of China(Nos.91330205and 11421101)the National Key Research and Development Program of China(No.2016YFB0200603)
文摘A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-balanced finite volume discretization of the steady-state SWEs is solved by the Newton method as the outer iteration and a geometric multigrid method with the block symmetric Gauss-Seidel smoother as the inner iteration. The proposed Newton multigrid method makes use of the local residual to regularize the Jacobian matrix of the Newton iteration, and can handle the steady- state problem with wet/dry transition. Several numerical experiments are conducted to demonstrate the efficiency, robustness, and well-balanced property of the proposed method. The relation between the convergence behavior of the Newton multigrid method and the distribution of the eigenvalues of the iteration matrix is detailedly discussed.
文摘In this paper, we obtained the numerical solutions of the modified regularized long-wave (MRLW) equation, by using the multigrid method and finite difference method. The solitary wave motion, interaction of two and three solitary waves, and development of the Maxwellian initial condition into solitary waves are studied using the proposed method. The numerical solutions are compared with the known analytical solutions. Usingerror norms and conservative properties of mass, momentum and energy, accuracy and efficiency of the mentioned method will be established through comparison with other techniques.
文摘This paper presents a numerical method for PDE-constrained optimization problems. These problems arise in many fields of science and engineering including those dealing with real applications. The physical problem is modeled by partial differential equations (PDEs) and involve optimization of some quantity. The PDEs are in most cases nonlinear and solved using numerical methods. Since such numerical solutions are being used routinely, the recent trend has been to develop numerical methods and algorithms so that the optimization problems can be solved numerically as well using the same PDE-solver. We present here one such numerical method which is based on simultaneous pseudo-time stepping. The efficiency of the method is increased with the help of a multigrid strategy. Application example is included for an aerodynamic shape optimization problem.
基金supported by the Natural Science Foundation of Hubei Province(CN)(Grant No.2019CFB693)the Research Foundation of the Education Department of Hubei Province(CN)(Grant No.B2019003)the open Foundation of the Key Laboratory of Metallurgical Equipment and Control of Education Ministry(CN)(Grant No.2015B14).
文摘The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.
基金Supported by National 863 Plan Project of Ministry of Science and Technology of China under Grant No. 2006AA09Z354National Natural Science Foundation of China under Grant No. 10672101.
文摘In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.
基金Supported by NSF of China(10971203)Supported by the NSF of the education Department of Henan Province (2009A110017)
文摘In this paper, a V-cycle multigrid method is presented for a Hermite rectangular element. By defining proper mesh-dependent inner product and transfer operator, we obtain its convergence property and the uniform convergence rate independent of mesh size and level are established.
文摘This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same as that of the SOR method and the CPU time is much less than that of the latter one.
文摘In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.
文摘Numerical solutions of the modified equal width wave equation are obtained by using the multigrid method and finite difference method. The motion of a single solitary wave, interaction of two solitary waves and development of the Maxwellian initial condition into solitary waves are studied using the proposed method. The numerical solutions are compared with the known analytical solutions. Using error norms and conservative properties of mass, momentum and energy, accuracy and efficiency of the mentioned method will be established through comparison with other methods.
文摘This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were discretized by Finite Volume Method(FVM)with second-order approximation scheme and deferred correction.Moreover,the algebraic method and the differential method were used to generate the non-orthogonal structured curvilinear grids.Furthermore,the influence of some parameters of geometric multigrid method,as well as lexicographical Gauss–Seidel(Lex-GS),η-line Gauss–Seidel(η-line-GS),Modified Strongly Implicit(MSI)and modified incomplete LU decomposition(MILU)solvers on the Central Processing Unit(CPU)time was investigated.Therefore,several parameters of multigrid method and solvers were tested for the problem,with the use of nonorthogonal structured curvilinear grids and multigrid method,resulting in an algorithm with the combination that achieved the best results and CPU time.The geometric multigrid method with Full Approximation Scheme(FAS),V-cycle and standard coarsening ratio for this problem were utilized.This article shows how to calculate the coordinates transformation metrics in the coarser grids.Results show that the MSI and MILU solvers are the most efficient.Moreover,theMSI solver is faster thanMILU for both grids generators;and the solutions are more accurate for the Burgers problem with grids generated using elliptic equations.
基金The project is supported by the Special Funds for Major State Basic Research Projects G19990328 and the National Natural Science Foundation of China(No.10071015)
文摘Cascadic multigrid technique for mortar Wilson finite element method ofhomogeneous boundary value planar linear elasticity is described and analyzed. Firstthe mortar Wilson finite element method for planar linear elasticity will be analyzed,and the error estimate under L2 and H1 norm is optimal. Then a cascadic multigridmethod for the mortar finite element discrete problem is described. Suitable grid trans-fer operator and smoother are developed which lead to an optimal cascadic multigridmethod. Finally, the computational results are presented.
基金supported by Educational Commission of Guangdong Province,China(No.2012LYM-0066)the National Social Science Foundation of China(No.14CJL016)
文摘In this paper a cascadic multigrid algorithm for the mortar finite element approximation of the semilinear elliptic problem is proposed,and corresponding theorems aregiven,which display the error estimate and the computational complexity of the method.
