The Immersed Interface Method for Navier-Stokes Equations with Interfaces in Cylindrical Coordinates

Abstract

Many three-dimensional physical applications can be better analyzed and solved using the cylindrical coordinates. In this paper, the immersed interface method (IIM) tailored for Navier-Stokes equations involving interfaces under the cylindrical coordinates has been developed. Note that, while the IIM has been developed for Stokes equations in the cylindrical coordinates assuming the axis-symmetry in the literature, there is a gap in dealing with Navier-Stokes equations, where the non-linear term includes an additional component involving the coordinate

\phi

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varphi $$\end{document}

, even if the geometry and force term are axis-symmetric. Solving the Navier-Stokes equations in cylindrical coordinates becomes challenging when dealing with interfaces that feature a discontinuous pressure and a non-smooth velocity, in addition to the pole singularity at

r=0

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r=0$$\end{document}

. In the newly developed algorithm, we have derived the jump conditions under the cylindrical coordinates. The numerical algorithm is based on a finite difference discretization on a uniform and staggered grid in the cylindrical coordinates. The finite difference scheme is standard away from the interface but is modified at grid points near and on the interface. As expected, the method is shown to be second-order accurate for the velocity. The developed new IIM is applied to the solution of some related fluid dynamic problems with interfaces.

Keywords

Immersed interface method (IIM) / Navier-Stokes equations / Axis-symmetric interface problem / Staggered grid / Pole singularity / Finite difference method

References

[1] Bell, J.B., Colella, P., Glaz, H.M.: A second-order projection method for the incompressible Navier-Stokes equations. J. Comput. Phys. 85(2), 257–283 (1989)

[2] Brown, D.L., Cortez, R., Minion, M.L.: Accurate projection methods for the incompressible Navier-Stokes equations. J. Comput. Phys. 168(2), 464–499 (2001)

[3] Ito, K., Li, Z., Wan, X.: Pressure jump conditions for Stokes equations with discontinuous viscosity in 2D and 3D. Methods Appl. Anal. 13, 199–214 (2006)

[4] Kim, J., Moin, P.: Application of a fractional-step method to incompressible Navier-Stokes equations. J. Comput. Phys. 59(2), 308–323 (1985)

[5] Lai, M.C., Li, Z.: A remark on jump conditions for the three-dimensional Navier-Stokes equations involving an immersed moving membrane. Appl. Math. Lett. 14(2), 149–154 (2001)

[6] Landau L.D., Lifshitz E.M.: Fluid Mechanics, Course of Theoretical Physics, vol 6, 2nd ed. Pergamon Press, Oxford (1987)

[7] LeVequeRJ, LiZ. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. Anal., 1994, 3141019-1044.

[8] LeVequeRJ, LiZ. Immersed interface methods for Stokes flow with elastic boundaries or surface tension. SIAM J. Sci. Comput., 1997, 183709-735.

[9] LiY, WilliamsSA, LaytonAT. A hybrid immersed interface method for driven Stokes flow in an elastic tube. Numer. Math. Theory Methods Appl., 2013, 64600-616.

[10] Li, Z., Ito, K.: The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics). Society for Industrial and Applied Mathematics, Philadelphia, PA, USA (2006)

[11] Li, Z., Lai, M.C.: The immersed interface method for the Navier-Stokes equations with singular forces. J. Comput. Phys. 171(2), 822–842 (2001)

[12] MittalR, IaccarinoG. Immersed boundary methods. Annu. Rev. Fluid Mech., 2005, 37: 239-261.

[13] Neustupa, J., Pokorný, M.: Axisymmetric flow of Navier-Stokes fluid in the whole space with non-zero angular velocity component. In: Proceedings of Partial Differential Equations and Applications (Olomouc, 1999), pp. 469–481 (2001)

[14] NganguiaH, YoungYN, LaytonAT, HuW-F, LaiM-C. An immersed interface method for axisymmetric electrohydrodynamic simulations in Stokes flow. Commun. Comput. Phys., 2015, 182429-449.

[15] PeskinCS. The immersed boundary method. Acta Numer., 2002, 11: 479-517.

[16] RuizJ, ChenJ, LiZ. The IIM in polar coordinates and its application to electro capacitance tomography problems. Numer. Algorithms, 2011, 573405-423.

[17] Ruiz, J., Li, Z.: The immersed interface method for axis-symmetric problems and application to the Hele-Shaw flow. Appl. Math. Comput. 264, 179–197 (2015)

[18] SchumannU, SweetR. A direct method for the solution of Poisson’s equation with Neumann boundary conditions on a staggered grid of arbitrary size. J. Comput. Phys., 1976, 20: 171-182.

[19] SwarztrauberP. A direct method for the discrete solution of separable elliptic equations. SIAM J. Numer. Anal., 1974, 11: 1136-1150.

[20] Swarztrauber, P., Sweet, R.: Efficient FORTRAN subprograms for the solution of elliptic equations. Technical Report TN/IA-109, National Center for Atmospheric Research (1975)

[21] SweetR. A cyclic reduction algorithm for solving block tridiagonal systems of arbitrary dimensions. SIAM J. Numer. Anal., 1977, 14: 706-720.

[22] Van KanJ. A second-order accurate pressure-correction scheme for viscous incompressible flow. SIAM J. Sci. Stat. Comput., 1986, 73870-891.

[23] ZeeuwPD. Matrix-dependent prolongations and restrictions in a blackbox multigrid solver. J. Comput. Appl. Math., 1990, 3311-27.

Funding

Simonsen Foundation(633724)

Innovative Research Group Project of the National Natural Science Foundation of China(12261070)

New England Biolabs Foundation Este trabajo es resultado de la estancia financiada por la Fundación Séneca-Agencia de Ciencia y Tecnología de la Región de Murcia con cargo al Programa Regional de Movilidad, Colaboración Internacional e Intercambio de Cono(21728/EE/22)

RIGHTS & PERMISSIONS

Shanghai University