Search Results

Now showing 1 - 3 of 3
  • Item
    On a reduced sparsity stabilization of grad-div type for incompressible flow problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Linke, Alexander; Rebholz, Leo
    We introduce a new operator for stabilizing error that arises from the weak enforcement of mass conservation in finite element simulations of incompressible flow problems. We show this new operator has a similar positive effect on velocity error as the well-known and very successful grad-div stabilization operator, but the new operator is more attractive from an implementation standpoint because it yields a sparser block structure matrix. That is, while grad-div produces fully coupled block matrices (i.e. block-full), the matrices arising from the new operator are block-upper triangular in two dimensions, and in three dimensions the 2,1 and 3,1 blocks are empty. Moreover, the diagonal blocks of the new operators matrices are identical to those of grad-div. We provide error estimates and numerical examples for finite element simulations with the new operator, which reveals the significant improvement in accuracy it can provide. Solutions found using the new operator are also compared to those using usual grad-div stabilization, and in all cases, solutions are found to be very similar.
  • Item
    Efficient linear solvers for incompressible flow simulations using Scott-Vogelius finite elements
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Cousins, Benjamin; Le Borne, Sabine; Linke, Alexander; Rebholz, Leo G.; Wang, Zhen
    Recent research has shown that in some practically relevant situations like multi-physics flows [11] divergence-free mixed finite elements may have a significantly smaller discretization error than standard nondivergence-free mixed finite elements. In order to judge the overall performance of divergence-free mixed finite elements, we investigate linear solvers for the saddle point linear systems arising in ((Pk)d; Pdisc k-1 )) Scott-Vogelius finite element implementations of the incompressible Navier-Stokes equations. We investigate both direct and iterative solver methods. Due to discontinuous pressure elements in the case of Scott-Vogelius elements, considerably more solver strategies seem to deliver promising results than in the case of standard mixed finite elements like Taylor-Hood elements. For direct methods, we extend recent preliminary work using sparse banded solvers on the penalty method formulation to finer meshes, and discuss extensions. For iterative methods, we test augmented Lagrangian and H
  • Item
    Simulation of composite materials by a Network FEM with error control
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Eigel, Martin; Peterseim, Daniel
    A novel Finite Element Method (FEM) for the computational simulation in particle reinforced composite materials with many inclusions is presented. It is based on a specially designed mesh consisting of triangles and channel-like connections between inclusions which form a network structure. The total number of elements and, hence, the number of degrees of freedom are proportional to the number of inclusions. The error of the method is independent of the possibly tiny distances of neighbouring inclusions. We present algorithmic details for the generation of the problem adapted mesh and derive an efficient residual a posteriori error estimator which enables to compute reliable upper and lower error bounds. Several numerical examples illustrate the performance of the method and the error estimator. In particular, it is demonstrated that the (common) assumption of a lattice structure of inclusions can easily lead to incorrect predictions about material properties.