Search Results

Now showing 1 - 7 of 7
  • Item
    Stable computing with an enhanced physics based scheme for the 3d Navier-Stokes equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Case, Michael; Ervin, V.J.; Linke, A.; Rebholz, L.G.; Wilson, N.E.
    We study extensions of an earlier developed energy and helicity preserving scheme for the 3D Navier-Stokes equations and apply them to a more general class of problems. The scheme is studied together with stabilizations of grad-div type in order to mitigate the effect of the Bernoulli pressure error on the velocity error. We prove stability, convergence, discuss conservation properties, and present numerical experiments that demonstrate the advantages of the scheme
  • Item
    An analogue of grad-div stabilization in nonconforming methods for incompressible flows
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Akbas, Mine; Linke, Alexander; Rebholz, Leo G.; Schroeder, Philipp W.
    Grad-div stabilization is a classical remedy in conforming mixed finite element methods for incompressible flow problems, for mitigating velocity errors that are sometimes called poor mass conservation. Such errors arise due to the relaxation of the divergence constraint in classical mixed methods, and are excited whenever the spacial discretization has to deal with comparably large and complicated pressures. In this contribution, an analogue of grad-div stabilization is presented for nonconforming flow discretizations of Discontinuous Galerkin or nonconforming finite element type. Here the key is the penalization of the jumps of the normal velocities over facets of the triangulation, which controls the measure-valued part of the distributional divergence of the discrete velocity solution. Furthermore, we characterize the limit for arbitrarily large penalization parameters, which shows that the proposed nonconforming Discontinuous Galerkin methods remain robust and accurate in this limit. Several numerical examples illustrate the theory and show their relevance for the simulation of practical, nontrivial flows.
  • Item
    On the parameter choice in grad-div stabilization for incompressible flow problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Jenkins, Eleanor W.; John, Volker; Linke, Alexander; Rebholz, Leo G.
    Grad-div stabilization has been proved to be a very useful tool in discretizations of incompressible flow problems. Standard error analysis for inf-sup stable conforming pairs of finite element spaces predicts that the stabilization parameter should be optimally chosen to be O(1). This paper revisits this choice for the Stokes equations on the basis of minimizing the H1( ) error of the velocity and the L2( ) error of the pressure. It turns out, by applying a refined error analysis, that the optimal parameter choice is more subtle than known so far in the literature. It depends on the used norm, the solution, the family of finite element spaces, and the type of mesh. Depending on the situation, the optimal stabilization parameter might range from being very small to very large. The analytic results are supported by numerical examples.
  • Item
    Improving mass conservation in FE approximations of the Navier Stokes equations using continuous velocity fields: a connection between grad-div stabilization and Scott-Vogelius elements
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Case, Michael A.; Ervin, V.J.; Linke, A.; Rebholz, L.G.
    This article studies two methods for obtaining excellent mass conservation in finite element computations of the Navier-Stokes equations using continuous velocity fields. Under mild restrictions, the Scott-Vogelius element pair has recently been shown to be inf-sup stable and have optimal approximation properties, while also providing pointwise mass conservation. We present herein the first numerical tests of this element pair for the time dependent Navier-Stokes equations. We also prove that, again under these mild restrictions, the limit of the grad-div stabilized Taylor-Hood solutions to the Navier-Stokes problem converges to the Scott-Vogelius solution as the stabilization parameter tends to infinity. That is, in this setting, we provide theoretical justification that choosing the parameter large does not destroy the solution. A limiting result is also proven for the general case. Numerical tests are provided which verify the theory, and show how both Scott-Vogelius and grad-div stabilized Taylor-Hood (with large stabilization parameter) elements can provide accurate results with excellent mass conservation for Navier-Stokes approximations.
  • Item
    Improving efficiency of coupled schemes for Navier-Stokes equations by a connection to grad-div stabilitzed projection methods
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Linke, Alexander; Neilan, Michael; Rebholz, Leo G.; Wilson, Nicholas E.
    We prove that in finite element settings where the divergence-free subspace of the velocity space has optimal approximation properties, the solution of Chorin/Temam projection methods for Navier-Stokes equations equipped with grad-div stabilization with parameter , converge to the associated coupled method solution with rate gamma as gamma -> ?. We prove this first for backward Euler schemes, and then extend the results to BDF2 schemes, and finally to schemes with outflow boundary conditions. Several numerical experiments are given which verify the convergence rate, and show how using projection methods in this setting with large grad-div stabilization parameters can dramatically improve accuracy.
  • Item
    On the grad-div stabilization for the steady Oseen and Navier-Stokes equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Ahmed, Naveed
    This paper studies the parameter choice in the grad-div stabilization applied to the generalized problems of Oseen type. Stabilization parameters based on minimizing the H1 (Omega) error of the velocity are derived which do not depend on the viscosity parameter. For the proposed parameter choices, the H1 (Omega) error of the velocity is derived that shows a direct dependence on the viscosity parameter. Differences and common features to the situation for the Stokes equations are discussed. Numerical studies are presented which confirm the theoretical results. Moreover, for the Navier-Stokes equations, numerical simulations were performed on a two-dimensional flow past a circular cylinder. It turns out, for the MINI element, that the best results can be obtained without grad-div stabilization.
  • Item
    Grad-div stabilization for the evolutionary Oseen problem with inf-sup stable finite elements
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Frutos, Javier de; García-Archilla, Bosco; John, Volker; Novo, Julia
    The approximation of the time-dependent Oseen problem using inf-sup stable mixed finite elements in a Galerkin method with grad-div stabilization is studied. The main goal is to prove that adding a grad-div stabilization term to the Galerkin approximation has a stabilizing effect for small viscosity. Both the continuous-in-time and the fully discrete case (backward Euler method, the two-step BDF, and CrankNicolson schemes) are analyzed. In fact, error bounds are obtained that do not depend on the inverse of the viscosity in the case where the solution is sufficiently smooth. The bounds for the divergence of the velocity as well as for the pressure are optimal. The analysis is based on the use of a specific Stokes projection. Numerical studies support the analytical results