6 results
Search Results
Now showing 1 - 6 of 6
- ItemStabilizing 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.
- ItemStable 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
- ItemCollision in a cross-shaped domain(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Linke, AlexanderIn 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.
- ItemImproving 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.
- ItemOn 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, AlexanderAccording 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.
- ItemOptimal and pressure-independent L2 velocity error estimates for a modified Crouzeix-Raviart Stokes element with BDM reconstructions(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Brennecke, ChristianNearly all inf-sup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the velocity error become pressure-dependent, while divergence-free mixed finite elements deliver pressure-independent estimates. A recently introduced new variational crime using lowest-order Raviart-Thomas velocity reconstructions delivers a much more robust modified Crouzeix-Raviart element, obeying an optimal pressure-independent discrete H1 velocity estimate. Refining this approach, a more sophisticated variational crime employing the lowest-order BDM element is proposed, which also allows proving an optimal pressure-independent L2 velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the Navier-Stokes case.