Accelerating Vector Iteration Methods

1986 ◽  
Vol 53 (2) ◽  
pp. 291-297 ◽  
Author(s):  
M. Papadrakakis

This paper describes a technique for accelerating the convergence properties of iterative methods for the solution of large sparse symmetric linear systems that arise from the application of finite element method. The technique is called partial preconditioning process (PPR) and can be combined with pure vector iteration methods, such as the conjugate gradient, the dynamic relaxation, and the Chebyshev semi-iterative methods. The proposed triangular splitting preconditioner combines Evans’ SSOR preconditioner with a drop-off tolerance criterion. The (PPR) is attractive in a FE framework because it is simple and can be implemented at the element level as opposed to incomplete Cholesky preconditioners, which require a sparse assembly. The method, despite its simplicity, is shown to be more efficient on a set of test problems for certain values of the drop-off tolerance parameter than the partial elimination method.

Author(s):  
Henry Piehl ◽  
Aleksandar-Saša Milaković ◽  
Sören Ehlers

Shipping in ice-covered regions has gained high attention within recent years. Analogous to weather routing, the occurrence of ice in a seaway affects the selection of the optimal route with respect to the travel time or fuel consumption. The shorter, direct path between two points—which may lead through an ice-covered area—may require a reduction of speed and an increase in fuel consumption. A longer, indirect route, could be more efficient by avoiding the ice-covered region. Certain regions may have to be avoided completely, if the ice thickness exceeds the ice-capability of the ship. The objective of this study is to develop a computational method that combines coastline maps, route cost information (e.g., ice thickness), transport task, and ship properties to find the optimal route between port of departure, A, and port of destination, B. The development approach for this tool is to formulate the transport task in the form of a potential problem, solve this equation with a finite element method (FEM), and apply line integration and optimization to determine the best route. The functionality of the method is first evaluated with simple test problems and then applied to realistic transport scenarios.


1969 ◽  
Vol 91 (3) ◽  
pp. 524-533 ◽  
Author(s):  
M. M. Reddi

A minimum principle for the transient incompressible Reynold’s equation, with the natural boundary conditions of prescribed pressure, as well as flow, is presented. The finite element method is introduced as the numerical counterpart of the Rayleigh-Ritz procedure. Flow computation is shown to be a natural corollary of the integral principle. Solutions of several test problems are presented.


2020 ◽  
Author(s):  
T Katsaounis ◽  
Dimitrios Mitsotakis ◽  
G Sadaka

© 2020 Elsevier Inc. In this paper we consider the numerical solution of Boussinesq-Peregrine type systems by the application of the Galerkin finite element method. The structure of the Boussinesq systems is explained and certain alternative nonlinear and dispersive terms are compared. A detailed study of the convergence properties of the standard Galerkin method, using various finite element spaces on unstructured triangular grids, is presented. Along with the study of the Peregrine system, a new Boussinesq system of BBM-BBM type is derived. The new system has the same structure in its momentum equation but differs slightly in the mass conservation equation compared to the Peregrine system. Further, the finite element method applied to the new system has better convergence properties, when used for its numerical approximation. Due to the lack of analytical formulas for solitary wave solutions for the systems under consideration, a Galerkin finite element method combined with the Petviashvili iteration is proposed for the numerical generation of accurate approximations of line solitary waves. Various numerical experiments related to the propagation of solitary and periodic waves over variable bottom topography and their interaction with the boundaries of the domains are presented. We conclude that both systems have similar accuracy when approximate long waves of small amplitude while the Galerkin finite element method is more efficient when applied to BBM-BBM type systems.


2004 ◽  
Vol 261-263 ◽  
pp. 537-542
Author(s):  
Shui Lin Wang ◽  
Xia Ting Feng ◽  
Xiu Run Ge

In finite element method, the order of complete polynomial of the interpolation function is related to the number of nodes in the element. This paper presents a four-node quadrilateral element with quadratic function on it. The presented displacement functions maintain C0-continuity. Meanwhile, the element stiffness matrix is derived from the displacement functions. Test problems show that high accuracy can be achieved by the use of the new displacement function on the element.


Sign in / Sign up

Export Citation Format

Share Document