Higher-order compact finite difference scheme with multigrid algorithm is applied in this paper for solving one-dimensional and two-dimensional inhomogeneous Helmholtz equations. In two-dimensional case, the suggested...
详细信息
Higher-order compact finite difference scheme with multigrid algorithm is applied in this paper for solving one-dimensional and two-dimensional inhomogeneous Helmholtz equations. In two-dimensional case, the suggested scheme has the stencil of twenty one points. An efficient solver multigrid method yields eighth-order accurate approximation on both fine and coarse grids. For the Neumann boundary condition, an eighth-order accurate representation is also developed. The accuracy and efficiency of eighth-order compact difference scheme are exhibited through graphical illustrations and computed results are drafted in tabular form to validate the numerical experiments.
High-order compact difference schemes can achieve higher-order accuracy on uniform grids. However, in some cases these may not achieve the desired accuracy. Therefore, we propose a multigrid method based on high-order...
详细信息
High-order compact difference schemes can achieve higher-order accuracy on uniform grids. However, in some cases these may not achieve the desired accuracy. Therefore, we propose a multigrid method based on high-order compact difference scheme on nonuniform grids. We will use interpolation and restriction operators developed by Ge and Cao (J. Comput. Phys. 230:4051-4070, 2011). The suggested scheme has up to fourth-order accuracy. Lastly, some numerical experiments are given to show the accuracy and performance of the proposed scheme.
For difference elliptic equations, an algorithm based on Fedorenko's multigrid method is constructed. The algorithm is intended for solving three-dimensional boundary value problems for equations with anisotropic ...
详细信息
For difference elliptic equations, an algorithm based on Fedorenko's multigrid method is constructed. The algorithm is intended for solving three-dimensional boundary value problems for equations with anisotropic discontinuous coefficients on parallel computers. Numerical results confirming the performance and parallel efficiency of the multigrid algorithm are presented. These qualities are ensured by using, as a multigrid triad, the standard Chebyshev iteration for coarsest grid equations, Chebyshev-type smoothing explicit iterative procedures, and intergrid transfer operators in problem-dependent form.
A fourth-order compact difference discretization scheme with unequal meshsizes in different coordinate directions is employed to solve a three-dimensional (3D) Poisson equation on a cubic domain. Two multgrid methods ...
详细信息
A fourth-order compact difference discretization scheme with unequal meshsizes in different coordinate directions is employed to solve a three-dimensional (3D) Poisson equation on a cubic domain. Two multgrid methods are developed to solve the resulting sparse linear systems. One is to use the full-coarsening multigrid method with plane Gauss-Seidel relaxation, which uses line Gauss-Seidel relaxation to compute each planewise solution. The other is to construct a partial semi-coarsening multigrid method with the traditional point or plane Gauss-Seidel relaxations. Numerical experiments are conducted to test the computed accuracy of the fourth-order compact difference scheme and the computational efficiency of the multigrid methods with the fourth-order compact difference scheme. (C) 2010 Elsevier Inc. All rights reserved.
In this paper,we will investigate a multigrid algorithm for poroelasticity problem by a new finite element method with homogeneous boundary conditions in two dimensional *** choose N´ed´elec edge element for...
详细信息
In this paper,we will investigate a multigrid algorithm for poroelasticity problem by a new finite element method with homogeneous boundary conditions in two dimensional *** choose N´ed´elec edge element for the displacement variable and piecewise continuous polynomials for the pressure variable in the model *** constructing multigrid algorithm,a distributive Gauss-Seidel iteration method is *** experiments shows that the finite element method achieves optimal convergence order and the multigrid algorithm is almost uniformly convergent to mesh size h and parameter dt on regular meshes.
We propose an efficient multigrid algorithm for solving anisotropic elliptic difference equations. The algorithm is based on using Chebyshev’s explicit iterations at smoothing stages and in solving coarse-grid equati...
详细信息
This paper deals with the task of pricing European basket options in the p-resence of transaction costs. We develop a model that incorporates the illiquidity of the market into the classical two-assets Black-Scholes f...
详细信息
This paper deals with the task of pricing European basket options in the p-resence of transaction costs. We develop a model that incorporates the illiquidity of the market into the classical two-assets Black-Scholes framework. We perform a nu-merical simulation using finite difference method. We consider a nonlinear multigrid method in order to reduce computational costs. The objective of this paper is to inves-tigate a deterministic extension for the Barles' and Soner's model and to demonstrate the effectiveness of multigrid approach to solving a fully nonlinear two dimensional Black-Scholes problem.
Predicting rolling bearing fatigue life requires knowledge of the three-dimensional(3D)stress fields in the roller and raceway near the lubricated *** to the increasingly severe operating conditions,the effect of loca...
详细信息
Predicting rolling bearing fatigue life requires knowledge of the three-dimensional(3D)stress fields in the roller and raceway near the lubricated *** to the increasingly severe operating conditions,the effect of localized features such as surface roughness,subsurface inclusions,and even the crystallographic structure of the material becomes *** such detail requires(locally)extremely dense gridding in simulations,which in 3D is a major *** techniques have been demonstrated to be capable of solving such *** this study,multigrid techniques are shown to further increase the efficiency of the solution by exploiting local grid refinement while maintaining the simplicity of a uniform *** is achieved by employing increasingly finer grids only locally,where the highest resolution is *** are presented for dry contact and elastohydrodynamically lubricated contact cases,circular as well as elliptic,with varying crystallographic structure,and with surface *** results show that the developed algorithm is very well suited for detailed analysis,with also excellent prospects for computational diagnostics involving actual material crystallographic structure from electron backscatter diffraction measurements.
This paper presents a new type of local and parallel multigrid method to solve semilinear elliptic equations. The proposed method does not directly solve the semilinear elliptic equations on each layer of the multigri...
详细信息
This paper presents a new type of local and parallel multigrid method to solve semilinear elliptic equations. The proposed method does not directly solve the semilinear elliptic equations on each layer of the multigrid mesh sequence, but transforms the semilinear elliptic equations into several linear elliptic equations on the multigrid mesh sequence and some low-dimensional semilinear elliptic equations on the coarsest mesh. Furthermore, the local and parallel strategy is used to solve the involved linear elliptic equations. Since solving large-scale semilinear elliptic equations in fine space, which can be fairly time-consuming, is avoided, the proposed local and parallel multigrid scheme will significantly improve the solving efficiency for the semilinear elliptic equations. Besides, compared with the existing multigrid methods which need the bounded second order derivatives of the nonlinear term, the proposed method only requires the Lipschitz continuation property of the nonlinear term. We make a rigorous theoretical analysis of the presented local and parallel multigrid scheme, and propose some numerical experiments to support the theory. (C) 2020 IMACS. Published by Elsevier B.V. All rights reserved.
In this paper, we study a V-cycle multigrid method for linear systems arising from time dependent two-dimensional space-fractional diffusion equations. The coefficient matrices of the linear systems are structured suc...
详细信息
In this paper, we study a V-cycle multigrid method for linear systems arising from time dependent two-dimensional space-fractional diffusion equations. The coefficient matrices of the linear systems are structured such that their matrix-vector multiplications can be computed efficiently. The main advantage using the multigrid method is to handle the space-fractional diffusion equations on non-rectangular domains, and to solve the linear systems with non-constant coefficients more effectively. The main idea of the proposed multigrid method is to employ two banded splitting iteration schemes as pre-smoother and post-smoother. The pre-smoother and the post-smoother take banded splitting of the coefficient matrix under the x-dominant ordering and the y-dominant ordering, respectively. We prove the convergence of the proposed two banded splitting iteration schemes in the sense of infinity norm. Results of numerical experiments for time dependent two-dimensional space-fractional diffusion equations on rectangular, L-shape and U-shape domains are reported to demonstrate that both computational time and iteration number required by the proposed method are significantly smaller than those of the other tested methods. (C) 2017 Elsevier Inc. All rights reserved.
暂无评论