We propose an efficient matrixdecomposition Method of Fundamental Solutions algorithm for the solution of certain two-dimensional linear elasticity problems. In particular, we consider the solution of the Cauchy-Navi...
详细信息
We propose an efficient matrixdecomposition Method of Fundamental Solutions algorithm for the solution of certain two-dimensional linear elasticity problems. In particular, we consider the solution of the Cauchy-Navier equations in circular domains subject to Dirichlet boundary conditions, that is when the displacements are prescribed on the boundary. The proposed algorithm is extended to the case of annular domains. Numerical experiments for both types of problems are presented.
The method of fundamental solutions is a boundary-type meshless method for the solution of certain elliptic boundary value problems. By exploiting the structure of the matrices appearing when this method is applied to...
详细信息
The method of fundamental solutions is a boundary-type meshless method for the solution of certain elliptic boundary value problems. By exploiting the structure of the matrices appearing when this method is applied to certain three-dimensional potential problems, we develop an efficient matrix decomposition algorithm for their solution. (C) 2003 Elsevier Ltd. All rights reserved.
matrix decomposition algorithms (MDAs) employing fast Fourier transforms are developed for the solution of the systems of linear algebraic equations arising when the finite element Galerkin method with piecewise Hermi...
详细信息
matrix decomposition algorithms (MDAs) employing fast Fourier transforms are developed for the solution of the systems of linear algebraic equations arising when the finite element Galerkin method with piecewise Hermite bicubics is used to solve Poisson's equation on the unit square. Like their orthogonal spline collocation counterparts, these MDAs, which require O(N (2)logN) operations on an NxN uniform partition, are based on knowledge of the solution of a generalized eigenvalue problem associated with the corresponding discretization of a two-point boundary value problem. The eigenvalues and eigenfunctions are determined for various choices of boundary conditions, and numerical results are presented to demonstrate the efficacy of the MDAs.
We consider a multi-level method of fundamental solutions for solving polyharmonic problems governed by 4N N u = 0 , N is an element of N\{1} in both two and three dimensions. Instead of approximating the solution wit...
详细信息
We consider a multi-level method of fundamental solutions for solving polyharmonic problems governed by 4N N u = 0 , N is an element of N\{1} in both two and three dimensions. Instead of approximating the solution with linear combinations of N fundamental solutions, we show that, with appropriate deployments of the source points, it is possible to employ an approximation involving only the fundamental solution of the operator 4N N . To determine the optimal position of the source points, we apply the recently developed effective condition number method. In addition, we show that when the proposed technique is applied to boundary value problems in circular or axisymmetric domains, with appropriate distributions of boundary and source points, it lends itself to the application of matrix decomposition algorithms. The results of several numerical tests are presented and analysed.
The framework has been made for a Movie Rating System to enable us to enjoy our time fruitfully. Film recommendation systems enable users to sort and select the best film fans to see without going through thousands of...
详细信息
A fast numerical algorithm is proposed for calculating the difference field radar cross-section of an object above a rough surface. The electric-field integral equation of the difference of the induced current on the ...
详细信息
A fast numerical algorithm is proposed for calculating the difference field radar cross-section of an object above a rough surface. The electric-field integral equation of the difference of the induced current on the rough surface and the induced current on the object is derived, which is solved by an iterative solver. The characteristic basis functions are used to calculate the induced current on the rough surface, which is part of the right-hand side of the system. It is observed that the coupling matrices (nonself-block interactions of the rough surface and the interactions between the object and the rough surface) are rank deficient. The matrix decomposition algorithm is used and extended to compress the coupling matrices. Numerical examples show that the characteristic basis functions method combined with the matrix decomposition algorithm can greatly reduce the CPU time and memory cost. The accuracy and efficiency are discussed by numerical examples.
The plane waves method is employed for the solution of Dirichlet and Neumann boundary value problems for the homogeneous Helmholtz equation in two- and three-dimensional domains possessing radial symmetry. The appropr...
详细信息
The plane waves method is employed for the solution of Dirichlet and Neumann boundary value problems for the homogeneous Helmholtz equation in two- and three-dimensional domains possessing radial symmetry. The appropriate selection of collocation points and unitary direction vectors in the method leads to circulant and block circulant coefficient matrices in two and three dimensions, respectively. We propose efficient matrix decomposition algorithms which make use of fast Fourier transforms for the solution of the systems resulting from such a discretization. In conjunction with the method of particular solutions, the method is extended to the solution of inhomogeneous axisymmetric Helmholtz problems. Several numerical examples are presented. (C) 2016 Elsevier Ltd. All rights reserved.
We study the computation of the orthogonal spline collocation solution of a linear Dirichlet boundary value problem with a nonselfadjoint or an indefinite operator of the form Lu = Sigmaa(ij) (x) u(xixj) + Sigma b(i)(...
详细信息
We study the computation of the orthogonal spline collocation solution of a linear Dirichlet boundary value problem with a nonselfadjoint or an indefinite operator of the form Lu = Sigmaa(ij) (x) u(xixj) + Sigma b(i)(x) u(xi) + c(x) u. We apply a preconditioned conjugate gradient method to the normal system of collocation equations with a preconditioner associated with a separable operator, and prove that the resulting algorithm has a convergence rate independent of the partition step size. We solve a problem with the preconditioner using an efficient direct matrix decomposition algorithm. On a uniform N x N partition, the cost of the algorithm for computing the collocation solution within a tolerance epsilon is O(N-2 lnN| ln epsilon|).
We solve the variable coefficient Cauchy-Navier equations of elasticity in the unit square, for Dirichlet and Dirichlet-Neumann boundary conditions, using second order finite difference schemes. The resulting linear s...
详细信息
We solve the variable coefficient Cauchy-Navier equations of elasticity in the unit square, for Dirichlet and Dirichlet-Neumann boundary conditions, using second order finite difference schemes. The resulting linear systems are solved by the preconditioned conjugate gradient (PCG) method with preconditioners corresponding to to the Laplace operator. The multiplication of a vector by the matrices of the resulting systems and the solution of systems with the preconditioners are performed at optimal and nearly optimal costs, respectively. For the case of Dirichlet boundary conditions, we prove the second order accuracy of the scheme in the discrete norm, symmetry of the resulting matrix and its spectral equivalence to the preconditioner. For the case of Dirichlet-Neumann boundary conditions, we prove symmetry of the resulting matrix. Numerical tests demonstrating the convergence properties of the schemes and PCG are presented.
Orthogonal spline collocation is implemented for the numerical solution of two-dimensional Helmholtz problems with discontinuous coefficients in the unit square. A matrix decomposition algorithm is used to solve the c...
详细信息
Orthogonal spline collocation is implemented for the numerical solution of two-dimensional Helmholtz problems with discontinuous coefficients in the unit square. A matrix decomposition algorithm is used to solve the collocation matrix system at a cost ofO(N-2 log N)on anN x Npartition of the unit square. The results of numerical experiments demonstrate the efficacy of this approach, exhibiting optimal global estimates in various norms and superconvergence phenomena for a broad spectrum of wave numbers.
暂无评论