O triangulacích bez tupých úhlů
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian.
On a parallel implementation of the mortar element method
On a Parallel Implementation of the Mortar Element Method
We discuss a parallel implementation of the domain decomposition method based on the macro-hybrid formulation of a second order elliptic equation and on an approximation by the mortar element method. The discretization leads to an algebraic saddle- point problem. An iterative method with a block- diagonal preconditioner is used for solving the saddle- point problem. A parallel implementation of the method is emphasized. Finally the results of numerical experiments are presented.
On adaptive BDDC for the flow in heterogeneous porous media
We study a method based on Balancing Domain Decomposition by Constraints (BDDC) for numerical solution of a single-phase flow in heterogeneous porous media. The method solves for both flux and pressure variables. The fluxes are resolved in three steps: the coarse solve is followed by subdomain solves and last we look for a divergence-free flux correction and pressures using conjugate gradients with the BDDC preconditioner. Our main contribution is an application of the adaptive algorithm for selection...
On algebraic multilevel methods for non-symmetric systems - convergence results.
On monotone and Schwarz alternating methods for nonlinear elliptic PDEs
The Schwarz alternating method can be used to solve elliptic boundary value problems on domains which consist of two or more overlapping subdomains. The solution is approximated by an infinite sequence of functions which results from solving a sequence of elliptic boundary value problems in each of the subdomains. In this paper, proofs of convergence of some Schwarz alternating methods for nonlinear elliptic problems which are known to have solutions by the monotone method (also known as the method...
On Monotone and Schwarz Alternating Methods for Nonlinear Elliptic PDEs
The Schwarz alternating method can be used to solve elliptic boundary value problems on domains which consist of two or more overlapping subdomains. The solution is approximated by an infinite sequence of functions which results from solving a sequence of elliptic boundary value problems in each of the subdomains. In this paper, proofs of convergence of some Schwarz alternating methods for nonlinear elliptic problems which are known to have solutions by the monotone method (also known as the method...
On nonoverlapping domain decomposition methods for the incompressible Navier-Stokes equations
In this paper, a Dirichlet-Neumann substructuring domain decomposition method is presented for a finite element approximation to the nonlinear Navier-Stokes equations. It is shown that the Dirichlet-Neumann domain decomposition sequence converges geometrically to the true solution provided the Reynolds number is sufficiently small. In this method, subdomain problems are linear. Other version where the subdomain problems are linear Stokes problems is also presented.
On nonoverlapping domain decomposition methods for the incompressible Navier-Stokes equations
In this paper, a Dirichlet-Neumann substructuring domain decomposition method is presented for a finite element approximation to the nonlinear Navier-Stokes equations. It is shown that the Dirichlet-Neumann domain decomposition sequence converges geometrically to the true solution provided the Reynolds number is sufficiently small. In this method, subdomain problems are linear. Other version where the subdomain problems are linear Stokes problems is also presented.
On selection of interface weights in domain decomposition methods
Different choices of the averaging operator within the BDDC method are compared on a series of 2D experiments. Subdomains with irregular interface and with jumps in material coefficients are included into the study. Two new approaches are studied along three standard choices. No approach is shown to be universally superior to others, and the resulting recommendation is that an actual method should be chosen based on properties of the problem.
On the accuracy of multigrid truncation error estimates.
On the convergence of multigrid methods for flow problems.
On the equivalence of primal and dual substructuring preconditioners.
On the multi-grid iteration for the eigenvalue problem and the degree of interpolation wich it requires (I).
On the representation of functions and finite difference operators on adaptive dyadic grids.
On the Schwarz algorithms for the elliptic exterior boundary value problems
Tuning the alternating Schwarz method to the exterior problems is the subject of this paper. We present the original algorithm and we propose a modification of it, so that the solution of the subproblem involving the condition at infinity has an explicit integral representation formulas while the solution of the other subproblem, set in a bounded domain, is approximated by classical variational methods. We investigate many of the advantages of the new Schwarz approach: a geometrical convergence...
On the Schwarz algorithms for the Elliptic Exterior Boundary Value Problems
Tuning the alternating Schwarz method to the exterior problems is the subject of this paper. We present the original algorithm and we propose a modification of it, so that the solution of the subproblem involving the condition at infinity has an explicit integral representation formulas while the solution of the other subproblem, set in a bounded domain, is approximated by classical variational methods. We investigate many of the advantages of the new Schwarz approach: a geometrical convergence...