Abstract
In this paper, we develop a stabilized difference finite element (SDFE) method for the 3D steady incompressible Navier-Stokes equations and apply Oseen iterative method to deal with the nonlinear term. Firstly, the finite difference discretization based on finite element pair \(P_1\times P_0\) in the z-direction is used to obtain the finite difference solution \((u^{n}_\tau =\sum ^{l_3}_{k=0}u^{nk}(x,y)\phi _k(z), p^{n}_\tau =\sum ^{l_3}_{k=1}p^{nk}(x,y)\psi _k(z))\), where \((u^{nk}, p^{nk})\) is the solution of 2D linearized Navier-Stokes equations, and then the stabilized finite element discretization based on finite element pair \((P_1,P_1,P_1)\times P_1\) in the (x, y)-plane is used to approximate \((u^{nk}, p^{nk})\), so as to obtain the SDFE solution \((u^{n}_h=\sum ^{l_3}_{k=0}u_{h}^{nk}(x,y)\phi _k(z), p^{n}_h=\sum ^{l_3}_{k=1}p_{h}^{nk}(x,y)\psi _k(z))\) of the 3D linearized Navier-Stokes equations. Our method has the following features. First, difference finite element method overcomes the difficulty of the 3D space discretization. Second, the stabilized method does not require specification of mesh-dependent parameters and retain the symmetry of the original equations. The rigorous stability analysis and error estimate are developed, showing that SDFE method is stable and has optimal convergence. Several numerical tests are presented, confirming the theoretical predictions and verifying the accuracy of the considered method.







Similar content being viewed by others
Data Availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].
References
Temam, R.: Navier-Stokes Equations, Theory and Numerical Analysis, 3rd edn. North-Holland, Amsterdam, New York and Oxford (1984)
Girault, V., Raviart, P.A.: Finite element methods for Navier-Stokes equations. Springer-Verlag, Berlin, Theory and algorithms (1986)
Heywood, J.G., Rannacher, R.: Finite element approximation of the nonstationary Navier-Stokes problem I. Regularity of solutions and second-order spatial discretization,. SIAM J. Numer. Anal. 19, 275–311 (1982)
Heywood, J.G., Rannacher, R.: Finite element approximation of the nonstationary Navier-Stokes problem. II. Stability of solutions and error estimates uniform in time,. SIAM J. Numer. Anal. 23, 750–777 (1986)
Heywood, J.G., Rannacher, R.: Finite element approximation of the nonstationary Navier-Stokes problem. III. Smoothing property and higher order estimates for spatial discretization,. SIAM J. Numer. Anal. 25(3), 489–512 (1988)
Heywood, J.G., Rannacher, R.: Finite element approximation of the nonstationary Navier-Stokes problem, IV: Error analysis for second order time discretizafion. SIAM J. Numel. Anal. 27(2), 353–384 (1990)
Quarteroni, A., Valli, A.: Numerical Approximation of Partial Differential Equations, vol. 23. Springer, (2008)
Glowinski, R.: Finite Element Methods for Incompressible Viscous Flow. Elsevier (2003)
Elman, H.C., Silvester, D.J., Wathen, A.J.: Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics. Oxford University Press (2014)
Layton, W.: A two level discretization method for the Navier-Stokes equations. Comput. Math. Appl. 26(2), 33–38 (1993)
He, Y., Wang, A., Chen, Z., Li, K.: An optimal nonlinear Galerkin method with mixed finite elements for the steady Navier-Stokes equations. Numer. Methods Part. Differ. Equ. 19(6), 762–775 (2003)
He, Y., Wang, A., Mei, L.: Stabilized finite element methods for the stationary Navier-Stokes equations. J. Engrg. Math. 51(4), 367–380 (2005)
He, Y., Li, K.: Two-level stabilized finite element methods for the steady Navier-Stokes equations. Computing 74(4), 337–351 (2005)
He, Y., Li, J.: Convergence of three iterative methods based on the finite element discretization for the stationary Navier-Stokes equations. Comput. Methods Appl. Mech. Engrg. 198, 1351–1359 (2009)
Hu, X., Mu, L., Ye, X.: A weak Galerkin finite element method for the Navier-Stokes equations on polytopal meshes. J. of Comput. Appl. Math. 362, 614–625 (2019)
Mu, L.: A pressure-robust weak Galerkin finite element method for Navier-Stokes equations, arXiv preprint arXiv:2011.11526 (2020)
Irisarri, D., Hauke, G.: Stabilized virtual element methods for the unsteady incompressible Navier-Stokes equations. Calcolo 56(4), 1–21 (2019)
Gatica, G.N., Munar, M., Sequeira, F.A.: A mixed virtual element method for the Navier-Stokes equations. Math. Models Methods Appl. Sci. 28(14), 2719–2762 (2018)
Liu, X., Chen, Z.: The nonconforming virtual element method for the Navier-Stokes equations. Adv. Comput. Math. 45(2), 51–74 (2019)
Chen, H., Li, K., Wang, S.: A dimension split method for the incompressible Navier-Stokes equations in three dimensions. Internat. J. Numer. Methods Fluids. 73(5), 409–435 (2013)
Chen, H., Li, K., Chu, Y., Chen, Z., Yang, Y.: A dimension splitting and characteristic projection method for three-dimensional incompressible flow. Discrete Contin. Dyn. Syst. Ser. B. 24(1), 127–147 (2019)
He, R., Feng, X., Chen, Z.: \(H^1\)-Superconvergence of a difference finite element method based on the \(P_{1}\)-\(P_{1}\)-conforming element on non-uniform meshes for the 3D Poisson equation. Math. Comp. 87(312), 1659–1688 (2018)
Feng, X., He, R., Chen, Z.: Superconvergence in \(H^{1}\)-norm of a difference finite element method for the heat equation in a 3D spatial domain with almost-uniform mesh. Numer. Algorithms. 86(1), 357–395 (2021)
Feng, X., Lu, X., He, Y.: Difference finite element method for the 3D steady Stokes equations. Appl. Numer. Math. 173, 418–433 (2022)
Feng, X., Lu, X., He, Y.: Difference finite element method for the 3D steady Navier-Stokes equations, SIAM J. Numel. Anal. (revised)
Medjo, T.T., Temam, R.: The two-grid finite difference method for the primitive equations of the ocean. Nonlinear Anal. 69(3), 1034–1056 (2008)
He, Y.: First order decoupled method of the 3D primitive equations of the ocean I: time discretization. J. Math. Anal. Appl. 412(2), 895–921 (2014)
He, Y.: Second order decoupled implicit/explicit method of the 3D primitive equations of the ocean I: time discretization. Int. J. Numer. Anal. Model. 12, 1–30 (2015)
He, Y., Zhang, Y., Xu, H., Chen, Z.: First-order decoupled finite element method of the 3D primitive equations of the ocean. SIAM J. Sci. Comput. 38, A273–A301 (2016)
He, Y., Xu, H., Chen, Z.: A second-order decoupled implicit/explicit method of the 3D primitive equations of ocean II: finite element spatial discretization. Int. J. Numer. Methods Eng. 108(7), 750–789 (2016)
Pei, Y.: Continuous data assimilation for the 3D primitive equations of the ocean, arXiv preprint arXiv:1805.06007. (2018)
He, Y.: On the solutions of the 3D steady and unsteady primitive equations of the ocean. J. Math. Anal. Appl. 491(1), 124243 (2020)
Babuska, I.: The finite element method with Lagrangian multipliers. Numer. Math. 20, 179–192 (1973)
Brezzi, F.: On the Existence, Uniqueness and Approximation of Saddle-point Problems arising from Lagrangian Multipliers,. Rev. Franc. Autom. Inf. Rech. Oper. Ser. Rouge 8(R–2), 129–151 (1974)
Brezzi, F., Pitkäranta, J.: On the stabilisation of finite element approximations of the Stokes problems. In: Hackbusch, W.: (eds) Efficient Solutions of Elliptic Systems, pp 11–19. Springer, (1984)
Brezzi, F., Douglas, J., Jr.: Stabilized mixed methods for the Stokes problem. Numer. Math. 53, 225–235 (1988)
Bochev, P.B., Dohrmann, C.R., Gunzburger, M.D.: Stabilization of low-order mixed finite elements for the Stokes equations. SIAM J. Numer. Anal. 44, 82–101 (2006)
Lu, X., Huang, P., Feng, X., He, Y.: A stabilized difference finite element method for the 3D steady Stokes equations, Appl. Math. Comput. 430, 127270 (2022)
Li, J., He, Y.: A stabilized finite element method based on two local Gauss integrations for the Stokes equations. J. Comp. Appl. Math. 214, 58–65 (2008)
Li, J., Chen, Z.: A new stabilized finite element method for the transient Navier-Stokes equations. Comput. Methods Appl. Mech. Engrg. 197, 22–35 (2007)
Hecht, F.: New development in freefem++. J. Numer. Math. 20(3-4), 251–265 (2012)
Acknowledgements
The authors would like to thank the editor and reviewers for their valuable comments and suggestions that greatly contributed to improving the quality of the present manuscript.
Funding
The work of Pengzhan Huang was supported by Natural Science Foundation (NSF) of China (No.11861067). The work of Xinlong Feng was supported by the Research Fund from Key Laboratory of Xinjiang Province, China (No. 2020D04002), the Scientific Research Plan of Universities in Xinjiang, China (No. XJEDU2020Y001, No.XJEDU2020I001) and the NSF of China (No.U19A2079, No.12071406). The work of Yinnian He was supported by the NSF of China (No. 11771348 and 12026257).
Author information
Authors and Affiliations
Contributions
All authors contributed to the study conception and design. Numerical analysis was performed by Xiaoli Lu and Yinnian He. The numerical simulation part was performed by Xiaoli Lu, Pengzhan Huang and Xinlong Feng. The first draft of the manuscript was written by Xiaoli Lu and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of Interest
The authors have no relevant financial or non-financial interests to disclose.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Lu, X., Huang, P., Feng, X. et al. A Stabilized Difference Finite Element Method for the 3D Steady Incompressible Navier-Stokes Equations. J Sci Comput 92, 104 (2022). https://doi.org/10.1007/s10915-022-01928-2
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-022-01928-2
Keywords
- 3D steady Navier-Stokes equations
- Oseen iterative method
- Stabilized difference finite element method
- Finite element pair \(((P_{1}, P_{1}, P_{1})\times P_{1})\times (P_{1}\times P_{0})P_{1}\)
- Optimal convergence