Skip to main content Accessibility help
×
Home

Stability of Projection Methods for Incompressible Flows Using High Order Pressure-Velocity Pairs of Same Degree: Continuous and Discontinuous Galerkin Formulations

  • E. Ferrer (a1), D. Moxey (a2), R. H. J. Willden (a3) and S. J. Sherwin (a2)

Abstract

This paper presents limits for stability of projection type schemes when using high order pressure-velocity pairs of same degree. Two high order h/p variational methods encompassing continuous and discontinuous Galerkin formulations are used to explain previously observed lower limits on the time step for projection type schemes to be stable [18], when h- or p-refinement strategies are considered. In addition, the analysis included in this work shows that these stability limits do not depend only on the time step but on the product of the latter and the kinematic viscosity, which is of particular importance in the study of high Reynolds number flows. We show that high order methods prove advantageous in stabilising the simulations when small time steps and low kinematic viscosities are used.

Drawing upon this analysis, we demonstrate how the effects of this instability can be reduced in the discontinuous scheme by introducing a stabilisation term into the global system. Finally, we show that these lower limits are compatible with Courant-Friedrichs-Lewy (CFL) type restrictions, given that a sufficiently high polynomial order or a mall enough mesh spacing is selected.

Copyright

Corresponding author

References

Hide All
[1]Arnold, D.N., Brezzi, F., Cockburn, B., and Marini, L.D.Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM Journal of Numerical Analysis, 39(5):17491779, 2001.
[2]Auteri, F., Guermond, J.L., and Parolini, N.Role of the LBB condition in weak spectral projection methods: 405. spournal of Computational Physics, 174(1):405420,2001.
[3]Babuška, I.The finite element method with Lagrangian multipliers. Numerische Mathematik, 20(3):179192,1973.
[4]Badia, S. and Codina, R.Algebraic pressure segregation methods for the incompressible Navier-Stokes equations. Archives of Computational Methods in Engineering, 15:152, 2007.
[5]Brezzi, F.An the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. ESAIM: Mathematical Modelling and Numerical Analysis, 8:129151,1974.
[6]Brezzi, F. and Bathe, K.J.A discourse on the stability conditions for mixed finite element formulations. Computational Methods in Applied Mechanics and Engineering, 82:2757, 1990.
[7]Brezzi, F., Cockburn, B., Marini, L.D., and Süli, E.Stabilization mechanisms in discontinuous Galerkin finite element methods. Computer Methods in Applied Mechanics and Engineering, 195(2528):32933310, 2006.
[8]Chorin, A.J.Numerical solution of the Navier-Stokes equations. Mathematics of Computation, 22:742762,1968.
[9]Cockburn, B., Kanschat, G., and Schötzau, D.An equal-order DG method for the incompressible Navier-Stokes equations. Journal of Scientific Computing, 40:188210,2009.
[10]Codina, R.Pressure stability in fractional step finite element methods for incompressible flows. Journal of Computational Physics, 140(1):112140,2001.
[11]Deville, M.O., Fischer, P.F., and Mund, E.H.High-order methods for incompressible fluid flow. Cambridge monographs on applied and computational mathematics. Cambridge University Press, 2002.
[12]Duffy, M.G.Quadrature over a pyramid or cube of integrands with a singularity at a vertex. SIAM Journal on Numerical Analysis, 19(6):12601262,1982.
[13]Elman, H.C., Silvester, D.J., and Wathen, A.J.Finite elements and fast iterative solvers: with applications in incompressible fluid dynamics. Oxford University Press, 2005.
[14]Epshteyn, Y. and Rivière, B.Estimation of penalty parameters for symmetric interior penalty Galerkin methods. Journal of Computational and Applied Mathematics, 206:843872, September 2007.
[15]Ferrer, E.A high order Discontinuous Galerkin - Fourier incompressible 3D Navier-Stokes solver with rotating sliding meshes for simulating cross-flow turbines. PhD thesis, University of Oxford, 2012.
[16]Ferrer, E. and R. Willden, H.J.A high order discontinuous Galerkin - Fourier incompressible 3D Navier-Stokes solver with rotating sliding meshes. Journal of Computational Physics, 231(21):70377056,2012.
[17]Ferrer, E. and Willden, R.H.J.A high order discontinuous Galerkin finite element solver for the incompressible Navier-Stokes equations. Computers & Fluids, 46(1):224230,2011.
[18]Guermond, J.L., Minev, P., and Shen, J.An overview of projection methods for incompressible flows. Computer Methods in Applied Mechanics and Engineering, 195(4447):60116045, 2006.
[19]Guermond, J.L. and Quartapelle, L.An stability and convergence of projection methods based on pressure Poisson equation. International Journal of Numerical Methods in Fluids, 26:10391053,1998.
[20]Guermond, J.L. and Shen, J.Velocity-correction projection methods for incompressible flows. SIAM Journal of Numerical Analysis, 41:112134, January 2003.
[21]Hartmann, R.Numerical analysis of higher order discontinuous Galerkin finite element methods. In Deconinck, H., editor, VKI LS 2008-08: CFD - ADIGMA course on very high order discretization methods, Oct. 1317, 2008. yon Karman Institute for Fluid Dynamics, Rhode Saint Genese, Belgium, 2008.
[22]Hu, Ning, Guo, Xian-Zhong, and Kat, I. NormanBounds, z.for eigenvalues and condition numbers in the p-version of the finite element method. Mathematics of Computation, 67(224):14231450,1998.
[23]Johnson, C.Numerical solutions of partial differential equations by the finite element method. Cambridge University Press, 1987.
[24]Karniadakis, G.E., Israeli, M., and Orszag, S.A.High-order splitting methods for incompressible Navier-Stokes equations. Journal of Computational Physics, 97:414443,1991.
[25]Karniadakis, G.E. and Sherwin, S.J.Spectral h/p element methods for computational fluid dynamics. Oxford Science Publications, 2005.
[26]Ladyzhenskaya, O.A.The mathematical theory of viscous incompressible flow. Gordon and Breach, New York, 2nd edition, 1969.
[27]Maday, Y., Meiron, D., Patera, A., and Ronquist, E.Analysis of iterative methods for the steady and unsteady Stokes problem: Application to spectral element discretizations. SIAM Journal on Scientific Computing, 14(2):310337,1993.
[28]Maday, Y., Patera, A., and Rønquist, E.M.An operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow. Journal of Scientific Computing, 5(4):263292,1990.
[29]Orszag, S.A., Israeli, M., and Deville, M.Boundary conditions for incompressible flows. Journal of Scientific Computing, 1(1):75111,1986.
[30]Rannacher, R.An Chorin’s projection method for the incompressible Navier-Stokes equations. volume 1530 of Lecture Notes in Mathematics, pages 167183. sppringer Berlin Heidelberg, 1992.
[31]Rivière, B.Discontinuous Galerkin methods for solving elliptic and parabolic equations: theory and implementation. Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 2008.
[32]Shahbazi, K., Fischer, P.F., and Ethier, C.R.A high-order discontinuous Galerkin method for the unsteady incompressible Navier-Stokes equations. Journal of Computational Physics, 222(1):391407, 2007.
[33]Shen, J.Pseudo-compressibility methods for the unsteady incompressible Navier-Stokes equations. In Guo, B, editor, Proceedings of the 1994 Beijing Symposium on Nonlinear Evolution Equations and Infinite Dynamical Systems, pages 6878. sphongShan University Press, 1997.
[34]Steinmoeller, D.T., Stastna, M., and Lamb, K.G.A short note on the discontinuous galerkin discretization of the pressure projection operator in incompressible flow. Journal of Computational Physics, 251(0):480486, 2013.
[35]Temam, A.R.Une methodé d’approximation de la solution des équations de Navier-Stokes. Bulletin de la Société Mathematique de France, 96:115152,1968.

Keywords

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed