Abstract
In this paper, a two-level Newton iterative method is proposed for the stationary micropolar fluid equations. Firstly, the original equations are solved on a coarse grid based on Newton-type linearization. Then, the simplified linearized equations are solved on a fine grid. The stability and error estimates of the method are given in the theoretical part. The results of the theoretical analysis show that when the coarse mesh size \(\varvec{H}\) and fine mesh size \(\varvec{h}\) satisfy the relation \(\varvec{h=O(H^{2}})\), the two-level Newton iterative method can achieve an optimal convergence rate. Finally, the effectiveness and applicability of the method are verified by some numerical experiments.













Similar content being viewed by others
Availability of data
All data is available.
Code availability
Not applicable.
References
Eringen, A.: Simple micro fluids. Int. J. Eng. Sci. 2, 205–217 (1964)
Eringen, A.: Theory of Micropolar Fluids. J. Math. Mech. 16, 1–18 (1966)
Łukaszewicz, G.: Micropolar fluids: theory and applications. Birkhäuser, Boston, MA, USA (1999)
Ariman, T.: On the analysis of blood flow. J. Biomech. 4, 185–192 (1971)
Galdi, G., Rionero, S.: A note on the existence and uniqueness of solutions of the micropolar fluid equations. Int. J. Eng. Sci. 2, 105–108 (1977)
Łukaszewicz, G.: On non-stationary flows of incompressible asymmetric fluids. Math. Meth. Appl. Sci. 13, 219–232 (1990)
Dong, B., Li, J., Wu, J.: Global well-posedness and large-time decay for the 2D micropolar equations. J. Differ. Equ. 262, 3488–3523 (2017)
Ortega-Torres, E., Rojas-Medar, M.: Optimal error estimate of the penalty finite element method for the micropolar fluid equations. Numer. Funct. Anal. Optim. 29, 612–637 (2008)
Nochetto, R., Salgado, A., Tomas, I.: The micropolar Navier-Stokes equations: a priori error analysis. Math. Mod. Meth. Appl. Sci. 24, 1237–1264 (2014)
Jiang, Y., Yang, Y.: Analysis of some projection methods for the incompressible fluids with microstructure. J. Korean Math. Soc. 55, 471–506 (2018)
Maimaiti, H., Liu, D.: Pressure-correction projection methods for the time-dependent micropolar fluids. Int. J. Numer. Meth. Fluids 94, 377–393 (2021)
Ren, Y., Liu, D.: A pressure correction projection finite element method for the 2D/3D time-dependent thermomicropolar fluid problem. Comput. Math. Appl. 136, 136–150 (2023)
Liu, J., Liu, D.: Numerical analysis and comparison of four stabilized finite element methods for the steady micropolar equations. Entropy 24, 454 (2022)
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. Eng. 198, 1351–1359 (2009)
Dong, X., He, Y., Zhang, Y.: Convergence analysis of three finite element iterative methods for the 2D/3D stationary incompressible magnetohydrodynamics. Comput. Methods Appl. Mech. Eng. 276, 287–311 (2014)
Xing, X., Liu, D.: Numerical analysis and comparison of three iterative methods based on finite element for the 2D/3D stationary micropolar fluid equations. Entropy 24, 628 (2022)
Xu, J.: A novel two two-grid method for semilinear elliptic equations, SIAM. J. Sci. Comput. 15, 231–237 (1994)
Xu, J.: Two-grid discretization techniques for linear and nonlinear PDEs. SIAM. J. Numer. Anal. 33, 1759–1777 (1996)
Layton, W.: A two-level discretization method for the Navier-Stokes equations. Comput. Math. Appl. 26, 33–38 (1993)
Layton, W., Lenferink, W.: Two-level Picard and modified Picard methods for the Navier-Stokes equations. Appl. Math. Comput. 69, 263–274 (1995)
Layton, W., Tobiska, L.: A two-level method with backtraking for the Navier-Stokes equations. SIAM. J. Numer. Anal. 35, 2035–2054 (1998)
Xu, H., He, Y.: Some iterative finite element methods for steady Navier-Stokes equations with different viscosities. J. Comput. Phys. 232, 136–152 (2013)
He, Y., Zhang, Y., Shang, Y., Xu, H.: Two-level Newton iterative method for the 2D/3D steady Navier-Stokes equations. Numer. Methods Partial. Differ. Equ. 28, 1620–1642 (2012)
He, Y., Li, J., Yang, X.: Two-level penalized finite element methods for the stationary Navier-Stoke equations. Int. J. Syst. Sci. 2, 131–143 (2006)
Girault, V., Lions, J.: Two-grid finite element scheme for the transient NavierCStokes problem. Math. Model. Numer. Anal. 35, 945–80 (2001)
He, Y.: Two-level method based on finite element and Crank-Nicolson extrapolation for the time-dependent Navier-Stokes equations. SIAM. J. Numer. Anal. 41, 1263–1285 (2003)
Layton, W., Lenferink, H., Peterson, J.: A two-level Newton finite element algorithm for approximating electrically conducting incompressible fluid flows. Comput. Math. Appl. 28, 21–31 (1994)
Dong, X., He, Y.: Two-level Newton iterative method for the 2D/3D stationary incompressible magneto-hydrodynamics. J. Sci. Comput. 63, 426–451 (2015)
Su, H., Feng, X., Zhao, J.: Two-level penalty Newton iterative method for the 2D/3D stationary incompressible magnetohydrodynamics equations. J. Sci. Comput. 70, 1144–1179 (2017)
Zhang, G., Zhang, Y., He, Y.: Two-level coupled and decoupled parallel correction methods for stationary incompressible magnetohydrodynamics. J. Sci. Comput. 65, 920–939 (2015)
Girault, V., Raviart, P.: Finite element methods for Navier-Stokes equations: Theory and Algorithms. Springer, Berlin (1986)
Prakash, J., Sinha, P.: Lubrication theory for micropolar fluids and its application to a journal bearing. Int. J. Eng. Sci. 13, 217–232 (1975)
Hecht, F., Pironneau, O., Le Hyaric, A., Ohtsuka, K.: Freefem++ Manual. Available online at: https://freefem.org/ (2012). Accessed 1 Oct 2023
Adams, R.: Sobolev Space. Academic Press, New York (1975)
Heywood, J., Rannacher, R.: Finite element approximation of the nonstationary Navier-Stokes problem. Part I. Regularity of solutions and second-order error estimates for spatial discretization. SIAM. J. Numer. Anal. 19, 275–311 (1982)
Cai, M., Mu, M., Xu, J.: Numerical solution to a mixed Navier-Stokes/Darcy model by the two-grid approach. SIAM. J. Numer. Anal. 47, 3325–3338 (2009)
Cai, M., Huang, P., Mu, M.: Some multilevel decoupled algorithms for a mixed Navier-Stokes/Darcy model. Adv. Comput. Math. 44, 115–145 (2018)
Huang, P., Cai, M., Wang, F.: A Newton type linearization based two grid method for coupling fluid flow with porous media flow. Appl. Numer. Math. 106, 182–198 (2016)
Eichel, H., Tobiska, L., Xie, H.: Supercloseness and superconvergence of stabilized low-order finite element discretizations of the Stokes Problem. Math. Comp. 80, 697–722 (2011)
John, V.: Finite element methods for incompressible flow problems. Springer International Publishing (2016)
Funding
National Natural Science Foundation of China (No. 12061075); Key Laboratory of Xinjiang Province (No. 2022D04014). National Key Research and Development Program of China (No. 2023YFB3001604)
Author information
Authors and Affiliations
Contributions
Formal analysis: XX and DL. Methodology: XX and DL.
Corresponding author
Ethics declarations
Ethics approval
Not applicable.
Consent to participate
All authors participated.
Consent for publication
All authors agreed to publish.
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
1.1 The Proof of Theorem 4
Proof
The first inequality (27) can be obtained by using Lemma 2.1.
Then, let us prove (28). Subtracting (25) from (3), we have
Taking \((\textbf{v},\varvec{\uppsi })=(R_{\mu }\textbf{u}-\textbf{u}_{\mu },Q_{\mu }\varvec{\upomega }-\varvec{\upomega }_{\mu }), q=\rho _{\mu }p-p_{\mu }\) in (53), we can derive that
For the left-hand side of (54), by using the (11)–(13), we have
Moreover, from (12), (13) and (18), we have
According to (54), we can get the equation of pressure
Combining the three equations above, we have
Then, by using (11)–(13) and discrete inf-sup condition (21), we can obtain
The \(L^2(\Omega )^d\) norm error estimate can be obtained by the dual technique. Firstly, we construct the dual problem. The error equation (53) is rewritten as
If we denote \((\textbf{u}-\textbf{u}_{\mu },\varvec{\upomega }-\varvec{\upomega }_{\mu })=(\textbf{v},\varvec{\uppsi })\), \(p-p_{\mu }=q\), \((\textbf{v}_{\mu },\varvec{\uppsi }_{\mu })=(\textbf{w},\varvec{\upphi })\), and \(q_{\mu }=s\), the dual problem can be given: find \(((\textbf{w},\varvec{\upphi }),s) \in \textbf{W} \times M\) such that
where \(((\textbf{v},\varvec{\uppsi }),\textbf{G})=B((\textbf{v},\varvec{\uppsi }),(\textbf{v},\varvec{\uppsi }),(\textbf{w},\varvec{\upphi }))\), \(\textbf{G}:= (G_{1},G_{2}) \in L^{2}(\Omega )^{d}\times L^{2}(\Omega )^{d}\).
According to Theorem 1, there exists a unique solution \(((\textbf{w},\varvec{\upphi }),s) \in \textbf{W} \times M\) of the dual problem (61), and it satisfies
Further, the projection operator \(((R_{\mu }\textbf{w}, Q_{\mu }\varvec{\upphi }),\rho _{\mu }s)\) satisfies
We set \(\textbf{G}=(\textbf{u}-\textbf{u}_{\mu },\varvec{\upomega }-\varvec{\upomega }_{\mu }), \ (\textbf{v},\varvec{\uppsi })=(\textbf{u}-\textbf{u}_{\mu },\varvec{\upomega }-\varvec{\upomega }_{\mu })\), and \(q=p-p_{\mu }\) in (61), it yields that
Let \((\textbf{v},\varvec{\uppsi })=(R_{\mu }\textbf{w}, Q_{\mu }\varvec{\upphi })\), \(q=\rho _{\mu }s\) in (3) and (25), we have
Making a difference between the above two formulas, we get
By using (11), (13), (57), and (64) to simplify the right-hand side of (67), we have
Finally, we can obtain
The proof is completed. \(\square \)
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) 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
Xing, X., Liu, D. A two-level iterative method with Newton-type linearization for the stationary micropolar fluid equations. Numer Algor 97, 475–501 (2024). https://doi.org/10.1007/s11075-023-01711-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-023-01711-w
Keywords
- Two-level method
- Newton iterative method
- Micropolar fluid equations
- Finite element method
- Error estimates