This paper introduces a kind of parallel multigrid method for solving Steklov eigenvalue problem based on the multilevel correction method. Instead of the common costly way of directly solving the Steklov eigenvalue p...
详细信息
This paper introduces a kind of parallel multigrid method for solving Steklov eigenvalue problem based on the multilevel correction method. Instead of the common costly way of directly solving the Steklov eigenvalue problem on some fine space, the new method contains some boundary value problems on a series of multilevel finite element spaces and some steps of solving Steklov eigenvalue problems on a very low dimensional space. The linear boundary value problems are solved by some multigrid iteration steps. We will prove that the computational work of this new scheme is truly optimal, the same as solving the corresponding linear boundary value problem. Besides, this multigrid scheme has a good scalability by using parallel computing technique. Some numerical experiments are presented to validate our theoretical analysis.
Standard multigrid methods are not so effective for equations with highly oscillatory coefficients. New coarse grid operators based on homogenized operators are introduced to restore the fast convergence rate of multi...
详细信息
Standard multigrid methods are not so effective for equations with highly oscillatory coefficients. New coarse grid operators based on homogenized operators are introduced to restore the fast convergence rate of multigrid methods. Finite difference approximations are used for the discretization of the equations. Convergence analysis is based on the homogenization theory. Proofs are given for a two-level multigrid method with the homogenized coarse grid operator for two classes of two-dimensional elliptic equations with Dirichlet boundary conditions.
We present a self-consistent ab initio simulation method to calculate coherent quantum transport through a molecule connected to metal electrodes in the linear-response regime. Density-functional theory (DFT) is appli...
详细信息
We present a self-consistent ab initio simulation method to calculate coherent quantum transport through a molecule connected to metal electrodes in the linear-response regime. Density-functional theory (DFT) is applied to the metal-molecule-metal system. The molecule and the metal electrodes are treated on the same footing as one extended molecule. The Full Approximation Scheme (FAS) nonlinear multigrid technique is used to accelerate convergence in a nonorthogonal localized orbital basis. The Landauer formula is employed to calculate the current with the transmission function obtained from a Green's function calculation. The current-voltage characteristics of a benzene-1,4-dithiolate (BDT) extended molecule are studied as an example, and our results are compared to other theoretical calculations. We also show that a recently formulated constrained-current formalism is invariant to a reversal in the imposed current. Hence, the predicted voltage drop must be zero. This suggests the theory must be modified to properly treat possible nonlinearities in the nonzero current case.
Finite-element multigrid method is applied to two-dimensional elastic wave equation in frequency domain. To avoid nonsensical increase in computational cost caused by global grid refinement, an adaptive algorithm is d...
详细信息
Finite-element multigrid method is applied to two-dimensional elastic wave equation in frequency domain. To avoid nonsensical increase in computational cost caused by global grid refinement, an adaptive algorithm is designed. Finally, numerical tests are carried out. (C) 2002 Elsevier Science Inc. All rights reserved.
Recently, nonlocal models attract the wide interest of scientists. They mainly come from two applied scientific fields: peridynamics and anomalous diffusion. Even though the matrices of the algebraic equation correspo...
详细信息
Recently, nonlocal models attract the wide interest of scientists. They mainly come from two applied scientific fields: peridynamics and anomalous diffusion. Even though the matrices of the algebraic equation corresponding to the nonlocal models are usually Toeplitz (denote a(0) as the principal diagonal element, a(1) as the trailing diagonal element, etc). There are still some differences for the models in these two fields. For the model of anomalous diffusion, a(0)/a(1) is uniformly bounded;most of the time, a(0)/a(1) of the model for peridynamics is unbounded as the step size h tends to zero. Based on the uniform boundedness of a(0)/a(1), the convergence of the two-grid method is well established [R. H. Chan, Q.-S. Chang, and H.-W. Sun, SIAM J. Sci. Comput., 19 (1998), pp. 516-529;H. Pang and H. Sun, J. Comput. Phys., 231 (2012), pp. 693-703;M. H. Chen, Y. T. Wang, X. Cheng, and W. H. Deng, BIT, 54 (2014), pp. 623-647]. This paper provides the detailed proof of the convergence of the two-grid method for the nonlocal model of peridynamics. Some special cases of the full multigrid and the V-cycle multigrid methods are also discussed. The numerical experiments are performed to verify the convergence.
Purpose As an important part of the rotor system, the damping coefficient of ball bearing has a great influence on the dynamic characteristics of the system. This study aims to propose a theoretical calculation method...
详细信息
Purpose As an important part of the rotor system, the damping coefficient of ball bearing has a great influence on the dynamic characteristics of the system. This study aims to propose a theoretical calculation method and an experimental test method to obtain the damping coefficient of ball bearing. Design/methodology/approach Based on Hertzian contact theory and elastohydrodynamic lubrication theory, the point contact oil film damping analysis model of ball bearing is established. The comprehensive damping calculation method considering external radial load, centrifugal force, ball spin, rotational speed and lubricating oil film is derived. The multigrid method is used to obtain the oil film pressure and thickness distribution in the contact zone. The variation trend of comprehensive damping with bearing radial load, rotational speed, oil film thickness and viscosity is analyzed. The test platform is designed and the influencing factors of damping are tested. Findings The validity of the model and reliability of the test device are verified by comparing the good consistency obtained in the work. The results show that the comprehensive damping of ball bearing increases with the increase of radial load and decreases with the increase of rotational speed. Originality/value At present, the existing bearing damping model can achieve approximate calculation of damping, but the factors considered in these models are not comprehensive enough. Besides, few studies exist regarding test platform of bearing damping, and a perfect test plan has not yet been formed. In this paper, the comprehensive damping calculation model of ball bearing is improved, and a complete experimental scheme is proposed to provide reference for the comprehensive damping theory and experimental research of bearing. Peer review The peer review history for this article is available at:
In this paper, we introduce an intergrid transfer operator which is based on the covolume of nodes in a P-1 nonconforming multigrid method and prove the convergence of the multigrid method with this intergrid transfer...
详细信息
In this paper, we introduce an intergrid transfer operator which is based on the covolume of nodes in a P-1 nonconforming multigrid method and prove the convergence of the multigrid method with this intergrid transfer operator. This intergrid transfer operator needs fewer computations and neighborhood node values than previous operators, which is a good property for parallelization. The P-1 nonconforming multigrid method with this intergrid transfer operator is suitable for solving problems with Robin boundary conditions and nonlinear problems with bound constraints on solutions. (C) 2004 IMACS. Published by Elsevier B.V. All rights reserved.
The gradient vector flow (GVF) deformable model was introduced by Xu and Prince as an effective approach to overcome the limited capture range problem of classical deformable models and their inability to progress int...
详细信息
The gradient vector flow (GVF) deformable model was introduced by Xu and Prince as an effective approach to overcome the limited capture range problem of classical deformable models and their inability to progress into boundary concavities. It has found many important applications in the area of medical image processing. The simple iterative method proposed in the original work on GVF, however, is slow to converge. A new multigrid method is proposed for GVF computation on 2D and 3D images. Experimental results show that the new implementation significantly improves the computational speed by at least an order of magnitude, which facilitates the application of GVF deformable models in processing large medical images.
作者:
Pan, KejiaHe, DongdongHu, HonglingCent S Univ
Sch Math & Stat Changsha 410083 Peoples R China Tongji Univ
Sch Aerosp Engn & Appl Mech Shanghai 200092 Peoples R China Hunan Normal Univ
Coll Math & Comp Sci Key Lab High Performance Comp & Stochast Informat Minist Educ China Changsha 410081 Hunan Peoples R China
Extrapolation cascadic multigrid (EXCMG) method is an efficient multigrid method which has mainly been used for solving the two-dimensional elliptic boundary value problems with linear finite element discretization in...
详细信息
Extrapolation cascadic multigrid (EXCMG) method is an efficient multigrid method which has mainly been used for solving the two-dimensional elliptic boundary value problems with linear finite element discretization in the existing literature. In this paper, we develop an EXCMG method to solve the three-dimensional Poisson equation on rectangular domains by using the compact finite difference (FD) method with unequal meshsizes in different coordinate directions. The resulting linear system from compact FD discretization is solved by the conjugate gradient (CG) method with a relative residual stopping criterion. By combining the Richardson extrapolation and tri-quartic Lagrange interpolation for the numerical solutions from two-level of grids (current and previous grids), we are able to produce an extremely accurate approximation of the actual numerical solution on the next finer grid, which can greatly reduce the number of relaxation sweeps needed. Additionally, a simple method based on the midpoint extrapolation formula is used for the fourth-order FD solutions on two-level of grids to achieve sixth-order accuracy on the entire fine grid cheaply and directly. The gradient of the numerical solution can also be easily obtained through solving a series of tridiagonal linear systems resulting from the fourth-order compact FD discretizations. Numerical results show that our EXCMG method is much more efficient than the classical V-cycle and W-cycle multigrid methods. Moreover, only few CG iterations are required on the finest grid to achieve full fourth-order accuracy in both the -norm and -norm for the solution and its gradient when the exact solution belongs to . Finally, numerical result shows that our EXCMG method is still effective when the exact solution has a lower regularity, which widens the scope of applicability of our EXCMG method.
In this paper, we present an efficient multigrid (MG) algorithm for solving the three-dimensional variable coefficient diffusion equation in cylindrical coordinates. The multigrid V-cycle combines a semi-coarsening in...
详细信息
In this paper, we present an efficient multigrid (MG) algorithm for solving the three-dimensional variable coefficient diffusion equation in cylindrical coordinates. The multigrid V-cycle combines a semi-coarsening in azimuthal direction with the red-black Gauss-Seidel plane (radial-axial plane) relaxation. On each plane relaxation, we further semi-coarsen the axial direction with red-black line relaxation in the radial direction. We also prove the convergence of two-level MG with plane Jacobi relaxation. Numerical results show that the present multigrid method indeed is scalable with the mesh size. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.
暂无评论