The topic of this work is to obtain discrete Sobolev inequalities for piecewise constant functions, and to deduce Lp error estimates on the approximate solutions of convection diffusion equations by finite volume schemes.
Our studies are motivated by a desire to model long-time simulations of possible scenarios for a waste disposal. Numerical methods are developed for solving the arising systems of convection-diffusion-dispersion-reaction equations, and the received results of several discretization methods are presented. We concentrate on linear reaction systems, which can be solved analytically. In the numerical methods, we use large time-steps to achieve long simulation times of about 10 000 years. We propose...
The idea of replacing a divergence constraint by a divergence boundary condition is investigated. The connections between the formulations are considered in detail. It is shown that the most common methods of using divergence boundary conditions do not always work properly. Necessary and sufficient conditions for the equivalence of the formulations are given.
By re-examining the arguments and counterexamples in I. Babuška, A. K. Aziz (1976) concerning the well-known maximum angle condition, we study the convergence behavior of the linear finite element method (FEM) on a family of distorted triangulations of the unit square originally introduced by H. Schwarz in 1880. For a Poisson problem with polynomial solution, we demonstrate arbitrarily slow convergence as well as failure of convergence if the distortion of the triangulations grows sufficiently fast....
Three non-overlapping domain decomposition methods are proposed for the numerical approximation of time-harmonic Maxwell equations with damping (i.e., in a conductor). For each method convergence is proved and, for the discrete problem, the rate of convergence of the iterative algorithm is shown to be independent of the number of degrees of freedom.
Three non-overlapping domain decomposition methods are proposed for the numerical approximation of time-harmonic Maxwell equations with damping (i.e., in a conductor). For each method convergence is proved and, for the discrete problem, the rate of convergence of the iterative algorithm is shown to be independent of the number of degrees of freedom.
This paper reviews the basic mathematical ideas and convergence analysis of domain decomposition methods. These are parallel and scalable iterative methods for the efficient numerical solution of partial differential equations. Two examples are then presented showing the application of domain decomposition methods to large-scale numerical simulations in computational mechanics and electrocardiology.
We present a parallel solution algorithm for the transient heat equation in one and two spatial dimensions. The problem is discretized in space by the lowest-order conforming finite element method. Further, a one-step time integration scheme is used for the numerical solution of the arising system of ordinary differential equations. For the latter, the parareal method decomposing the time interval into subintervals is employed. It leads to parallel solution of smaller time-dependent problems. At...
An axisymmetric second order elliptic problem with mixed boundary conditions is considered. The shape of the domain has to be found so as to minimize a cost functional, which is given in terms of the cogradient of the solution. A new dual finite element method is used for approximate solutions. The existence of an optimal domain is proven and a convergence analysis presented.