On numerical solution of compressible flow in time-dependent domains

Feistauer, Miloslav ; Horáček, Jaromír ; Kučera, Václav ; Prokopová, Jaroslava

Mathematica Bohemica, Tome 137 (2012) no. 1, pp. 1-16

Abstract (VO)

The paper deals with numerical simulation of a compressible flow in time-dependent 2D domains with a special interest in medical applications to airflow in the human vocal tract. The mathematical model of this process is described by the compressible Navier-Stokes equations. For the treatment of the time-dependent domain, the arbitrary Lagrangian-Eulerian (ALE) method is used. The discontinuous Galerkin finite element method (DGFEM) is used for the space semidiscretization of the governing equations in the ALE formulation. The time discretization is carried out with the aid of a linearized semi-implicit method with good stability properties. We present some computational results for the flow in a channel, representing a model of glottis and a part of the vocal tract, with a prescribed motion of the channel walls at the position of vocal folds.

Détail

DOI : 10.21136/MB.2012.142782

Classification : 65M60, 76M10, 76N15
Keywords: compressible Navier-Stokes equations; arbitrary Lagrangian-Eulerian method; discontinuous Galerkin finite element method; interior and boundary penalty; semi-implicit time discretization; biomechanics of voice

Bibliographie

[1] Arnold, D. N.: An interior penalty finite element method with discontinuous elements. SIAM J. Numer. Anal. 19 (1982), 742-760. | DOI | MR | Zbl

[2] Arnold, D. N., Brezzi, F., Cockburn, B., Marini, D.: Discontinuos Galerkin methods for elliptic problems. Discontinuous Galerkin methods. Theory, Computation and Applications. Lecture Notes in Computational Science and Engineering 11 Cockburn, B., et al. Springer, Berlin (2000), 89-101. | MR

[3] Arnold, D. N., Brezzi, F., Cockburn, B., Marini, D.: Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal. 39 (2001), 1749-1779. | DOI | MR

[4] Bassi, F., Rebay, S.: High-order accurate discontinuous finite element solution of the 2{D} Euler equations. J. Comput. Phys. 138 (1997), 251-285. | DOI | MR | Zbl

[5] Baumann, C. E., Oden, J. T.: A discontinuous $hp$ finite element method for the Euler and Navier-{S}tokes equations. Int. J. Numer. Methods Fluids 31 (1999), 79-95. | DOI | MR | Zbl

[6] Davis, T. A., Duff, I. S.: A combined unifrontal/multifrontal method for unsymmetric sparse matrices. ACM Transactions on Mathematical Software 25 (1999), 1-20. | DOI | MR | Zbl

[7] Dolejší, V.: Semi-implicit interior penalty discontinuous Galerkin methods for viscous compressible flows. Commun. Comput. Phys. 4 (2008), 231-274. | MR

[8] Dolejší, V., Feistauer, M.: A semi-implicit discontinuous Galerkin finite element method for the numerical solution of inviscid compressible flow. J. Comput. Phys. 198 (2004), 727-746. | DOI | MR | Zbl

[9] Dolejší, V., Feistauer, M.: Error estimates of the discontinuous Galerkin method for nonlinear nonstationary convection-diffusion problems. Numer. Func. Anal. Optimiz. 26 (2005), 349-383. | DOI | MR | Zbl

[10] Dolejší, V., Feistauer, M., Hozman, J.: Analysis of semi-implicit DGFEM for nonlinear convection-diffusion problems on nonconforming meshes. Comput. Methods Appl. Mech. Engrg. 196 (2007), 2813-2827. | DOI | MR | Zbl

[11] Dolejší, V., Feistauer, M., Schwab, C.: On some aspects of the discontinuous Galerkin finite element method for conservation laws. Math. Comput. Simul. 61 (2003), 333-346. | DOI | MR | Zbl

[12] Feistauer, M., Dolejší, V., Kučera, V.: On the discontinuous Galerkin method for the simulation of compressible flow with wide range of Mach numbers. Comput. Vis. Sci. 10 (2007), 17-27. | DOI | MR

[13] Feistauer, M., Felcman, J., Straškaraba, I.: Mathematical and Computational Methods for Compressible Flow. Clarendon Press, Oxford (2003). | MR

[14] Feistauer, M., Kučera, V.: On a robust discontinuous Galerkin technique for the solution of compressible flow. J. Comput. Phys. 224 (2007), 208-221. | DOI | MR | Zbl

[15] Krivodonova, L., Berger, M.: High-order accurate implementation of solid wall boundary conditions in curved geometries. J. Comput. Phys. 211 (2006), 492-512. | DOI | MR | Zbl

[16] Neustupa, J.: Existence of a weak solution to the Navier-Stokes equation in a general time-varying domain by the Rothe method. Math. Meth. Appl. Sci. 32 (2009), 653-683. | DOI | MR | Zbl

[17] Nomura, T., Hughes, T. J. R.: An arbitrary Lagrangian-{E}ulerian finite element method for interaction of fluid and a rigid body. {Comput. Methods Appl. Mech. Engrg.} 95 (1992), 115-138. | DOI | Zbl

[18] Prokopová, J.: Numerical solution of compressible flow. Master thesis, Charles University, Praha (2008).

[19] Punčochářová, P., Fürst, J., Kozel, K., Horáček, J.: Numerical solution of compressible flow with low Mach number through oscillating glottis. Proceedings of the 9th International Conference on Flow-Induced Vibration (FIV 2008), Praha, Institute of Thermomechanics AS CR (2008), 135-140. | MR

[20] Sváček, P., Feistauer, M., Horáček, J.: Numerical simulation of flow induced airfoil vibrations with large amplitudes. J. Fluids Structures 23 (2007), 391-411. | DOI

[21] Vegt, J. J. W. van der, Ven, H. van der: Space-time discontinuous Galerkin finite element method with dynamic grid motion for inviscid compressible flow. J. Comput. Phys. 182 (2002), 546-585. | DOI | MR

[22] Vijayasundaram, G.: Transonic flow simulation using upstream centered scheme of Godunov type in finite elements. J. Comput. Phys. 63 (1986), 416-433. | DOI | MR

[23] Zolésio, J. P.: Approximation for the wave equation in a moving domain. Proceedings of the conference Control of Partial Differential Equations. IFIP WG 7.2, Marcel Dekker, Lect. Notes Pure Appl. Math. 165, New York (1994), 271-279. | MR | Zbl