Skip to main content Accessibility help
×
Home
Hostname: page-component-568f69f84b-jtg5s Total loading time: 0.26 Render date: 2021-09-20T15:44:50.366Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "metricsAbstractViews": false, "figures": true, "newCiteModal": false, "newCitedByModal": true, "newEcommerce": true, "newUsageEvents": true }

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

Published online by Cambridge University Press:  03 June 2015

E. Ferrer*
Affiliation:
ETSIA – School of Aeronautics, Universidad Politécnica de Madrid, Pza. Cardenal Cisneros 3, Madrid, E-28040, Spain
D. Moxey*
Affiliation:
Department of Aeronautics, Imperial College London, London, SW7 2AZ, UK
R. H. J. Willden*
Affiliation:
Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK
S. J. Sherwin*
Affiliation:
Department of Aeronautics, Imperial College London, London, SW7 2AZ, UK
Get access

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.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[3]Babuška, I.The finite element method with Lagrangian multipliers. Numerische Mathematik, 20(3):179192,1973.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[8]Chorin, A.J.Numerical solution of the Navier-Stokes equations. Mathematics of Computation, 22:742762,1968.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[10]Codina, R.Pressure stability in fractional step finite element methods for incompressible flows. Journal of Computational Physics, 140(1):112140,2001.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.3.0.CO;2-U>CrossRefGoogle Scholar
[20]Guermond, J.L. and Shen, J.Velocity-correction projection methods for incompressible flows. SIAM Journal of Numerical Analysis, 41:112134, January 2003.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[23]Johnson, C.Numerical solutions of partial differential equations by the finite element method. Cambridge University Press, 1987.Google Scholar
[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.CrossRefGoogle Scholar
[25]Karniadakis, G.E. and Sherwin, S.J.Spectral h/p element methods for computational fluid dynamics. Oxford Science Publications, 2005.CrossRefGoogle Scholar
[26]Ladyzhenskaya, O.A.The mathematical theory of viscous incompressible flow. Gordon and Breach, New York, 2nd edition, 1969.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[29]Orszag, S.A., Israeli, M., and Deville, M.Boundary conditions for incompressible flows. Journal of Scientific Computing, 1(1):75111,1986.CrossRefGoogle Scholar
[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.Google Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
33
Cited by

Send article to Kindle

To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

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

Send article to Dropbox

To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

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

Send article to Google Drive

To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

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

Reply to: Submit a response

Please enter your response.

Your details

Please enter a valid email address.

Conflicting interests

Do you have any conflicting interests? *