In this paper we present an error estimator for unilateral contact problems solved by a Neumann-Neumann domain decomposition algorithm. This error estimator takes into account both the spatial error due to the finite ...
详细信息
In this paper we present an error estimator for unilateral contact problems solved by a Neumann-Neumann domain decomposition algorithm. This error estimator takes into account both the spatial error due to the finite element discretization and the algebraic error due to the domain decomposition algorithm. To differentiate specifically the contribution of these two error sources to the global error, two quantities are introduced: a discretization error indicator and an algebraic error indicator. The effectivity indices and the convergence of both the global error estimator and the error indicators are shown on several examples.
A finite difference domain decomposition algorithm on a non-overlapping non-matching grid for the parabolic equation is discussed. The basic procedure is to define the explicit scheme at the interface points with a la...
详细信息
A finite difference domain decomposition algorithm on a non-overlapping non-matching grid for the parabolic equation is discussed. The basic procedure is to define the explicit scheme at the interface points with a larger mesh spacing H, then the implicit schemes with different mesh spacings are applied on the non-matching subdomains, respectively. The stability bound is released both for the one-dimensional and two-dimensional parabolic problem. Finally, numerical experiments are also presented.
An improved algorithm is presented to overcome the shortage of the domain decomposition algorithm based on equivalence principle for multiple bodies of revolutions (EPA-MBoR). In a dense distribution of non-coaxal bod...
详细信息
An improved algorithm is presented to overcome the shortage of the domain decomposition algorithm based on equivalence principle for multiple bodies of revolutions (EPA-MBoR). In a dense distribution of non-coaxal bodies of revolution (BoRs) containing at least two very close BoRs, no matter how to establish sphere equivalence surfaces (ESs), they are intersecting. It violates the equivalence principle and causes error results. The improvement is introduced to overcome the defect. In the improved algorithm, a sphere ES cluster is built on each BoR instead of a single sphere ES. We can freely adjust the sphere ES clusters on each BoR to avoid the intersection among them. The numerical results verify the above improvement successfully overcomes the deficiency and show the improved algorithm inherits the high efficiency of the previous algorithm.
Schwarz methods are an important type of domaindecomposition methods. Using the Fourier transform, we derive error propagation matrices and their spectral radii of the classical Schwarz alternating method and the add...
详细信息
Schwarz methods are an important type of domaindecomposition methods. Using the Fourier transform, we derive error propagation matrices and their spectral radii of the classical Schwarz alternating method and the additive Schwarz method for the biharmonic equation in this paper. We prove the convergence of the Schwarz methods from a new point of view, and provide detailed information about the convergence speeds and their dependence on the overlapping size of subdomains. The obtained results are independent of any unknown constant and discretization method, showing that the Schwarz alternating method converges twice as quickly as the additive Schwarz method.
In this paper, the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential equation. The explicit group (EG) met...
详细信息
ISBN:
(纸本)9780769550602
In this paper, the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential equation. The explicit group (EG) method is derived from the standard centered seven-point finite difference discretisation formula. We utilize the domaindecomposition technique on this group scheme to divide the tasks involved in solving the equation. The aim of this study is to describe the development of the parallel group iterative scheme under OpenMP programming environment as a way to reduce the computational costs of the solution processes using multiple-core technologies. Numerical experiments are conducted together with their detailed performance analysis. The results will be reported and discussed.
In recent years several implementations of molecular dynamics (MD) codes have been reported on multiple instruction multiple data (MIMD) machines. However, very few implementations of MD codes on single instruction mu...
详细信息
In recent years several implementations of molecular dynamics (MD) codes have been reported on multiple instruction multiple data (MIMD) machines. However, very few implementations of MD codes on single instruction multiple data (SIMD) machines have been reported. The difficulty in using pair lists of nonbonded interactions is the major problem with MD codes for SIMD machines, such that, generally, the full connectivity computation has been used. We present an algorithm, the global cut-off algorithm (GCA), which permits the use of pair lists on SIMD machines. GCA is based on a probabilistic approach and requires the cut-off condition to be simultaneously verified on all nodes of the machine. The MD code used was taken from the GROMOS package;only the routines involved in the pair lists and in the computation of nonbonded interactions were rewritten for a parallel architecture. The remaining calculations were performed on the host computer. The algorithm has been tested on Quadrics computers for configurations of 32, 128, and 512 processors and for systems of 4000, 8000, 15,000, and 30,000 particles. Quadrics was developed by Istituto Nazionale di Fisica Nucleare (INFN) and marketed by Alenia Spazio. (C) 1998 John Wiley & Sons, Inc.
In this paper, the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential equation. The explicit group(EG) meth...
详细信息
In this paper, the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential equation. The explicit group(EG) method is derived from the standard centered seven-point finite difference discretisation formula. We utilize the domaindecomposition technique on this group scheme to divide the tasks involved in solving the equation. The aim of this study is to describe the development of the parallel group iterative scheme under OpenMP programming environment as a way to reduce the computational costs of the solution processes using multiplecore technologies. Numerical experiments are conducted together with their detailed performance analysis. The results will be reported and discussed.
In this paper,the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential *** explicit group(EG)method is derive...
详细信息
In this paper,the parallel implementation of a new explicit group iterative scheme is proposed for the solution of a three dimensional second order telegraph partial differential *** explicit group(EG)method is derived from the standard centered seven-point finite difference discretisation *** utilize the domaindecomposition technique on this group scheme to divide the tasks involved in solving the *** aim of this study is to describe the development of the parallel group iterative scheme under OpenMP programming environment as a way to reduce the computational costs of the solution processes using multiplecore *** experiments are conducted together with their detailed performance *** results will be reported and discussed.
The aim of this paper is to show procedure of mathematical modelling of weight-bearing human joints and numerical simulation of mechanical processes taking place during static loadening. It is important to know stress...
详细信息
ISBN:
(纸本)9788096956241
The aim of this paper is to show procedure of mathematical modelling of weight-bearing human joints and numerical simulation of mechanical processes taking place during static loadening. It is important to know stress and deformation distributions in lower limb for example for application of the total knee or the total hip replacements. Afinite element approximation and domaindecomposition method are used.
暂无评论