Search Results

Now showing 1 - 3 of 3
  • Item
    Stabilizing poor mass conservation in incompressible flow problems with large irrotational forcing and application to thermal convection
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2011) Galvin, Keith J.; Linke, Alexander; Rebholz, Leo G.; Wilson, Nicholas E.
    We consider the problem of poor mass conservation in mixed finite element algorithms for flow problems with large rotation-free forcing in the momentum equation. We provide analysis that suggests for such problems, obtaining accurate solutions necessitates either the use of pointwise divergence-free finite elements (such as Scott-Vogelius), or heavy grad-div stabilization of weakly divergence-free elements. The theory is demonstrated in numerical experiments for a benchmark natural convection problem, where large irrotational forcing occurs with high Rayleigh numbers.
  • Item
    Collision in a cross-shaped domain
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Linke, Alexander
    In the numerical simulation of the incompressible Navier-Stokes equations different numerical instabilities can occur. While instability in the discrete velocity due to dominant convection and instability in the discrete pressure due to a vanishing discrete LBB constant are well-known, instability in the discrete velocity due to a poor mass conservation at high Reynolds numbers sometimes seems to be underestimated. At least, when using conforming Galerkin mixed finite element methods like the Taylor-Hood element, the classical grad-div stabilization for enhancing discrete mass conservation is often neglected in practical computations. Though simple academic flow problems showing the importance of mass conservation are well-known, these examples differ from practically relevant ones, since specially designed force vectors are prescribed. Therefore we present a simple steady Navier-Stokes problem in two space dimensions at Reynolds number 1024, a colliding flow in a cross-shaped domain, where the instability of poor mass conservation is studied in detail and where no force vector is prescribed.
  • Item
    On the role of the Helmholtz-decomposition in mixed methods for incompressible flows and a new variational crime
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Linke, Alexander
    According to the Helmholtz decomposition, the irrotational parts of the momentum balance equations of the incompressible Navier-Stokes equations are balanced by the pressure gradient. Unfortunately, nearly all mixed methods for incompressible flows violate this fundamental property, resulting in the well-known numerical instability of poor mass conservation. The origin of this problem is the lack of L2-orthogonality between discretely divergence-free velocities and irrotational vector fields. In order to cure this, a new variational crime using divergence-free velocity reconstructions is proposed. Applying lowest order Raviart-Thomas velocity reconstructions to the nonconforming Crouzeix-Raviart element allows to construct a cheap flow discretization for general 2d and 3d simplex meshes that possesses the same advantageous robustness properties like divergence-free flow solvers. In the Stokes case, optimal a-priori error estimates for the velocity gradients and the pressure are derived. Moreover, the discrete velocity is independent of the continuous pressure. Several detailed linear and nonlinear numerical examples illustrate the theoretical findings.