×

Domain decomposition based exponential time differencing method for fluid dynamics problems with smooth solutions. (English) Zbl 1458.76070

Summary: The exponential time differencing method (ETD) has been recently introduced for gas dynamics problems, through which significant speedup against the popular TVD Runge-Kutta scheme can be gained. In this paper, a local ETD time advancing method is introduced in solving fluid dynamics problems with smooth solutions, which has the natural characteristic of parallelism. The application of four ETD schemes in conjunction with weighted essentially non-oscillation spatial discrete scheme is discussed, in which case the nearly analytic Jacobian matrix is obtained using the chain rule, and the matrix exponential computation is approximated through Krylov subspace projection method. For large scale simulations, the local version of four ETD schemes are introduced on the basis of overlapping domain decomposition method. Addictive Schwarz iteration is applied to decouple the multi-domain problem, and stationary problems are solved on the sub-domains in every time step. Numerical tests against several one- and two-dimensional problems demonstrate that the present method is accurate and efficient. Schwarz iteration does not increase the computation load much, and the numerical error of local ETD method will converge when overlap size is large enough.

MSC:

76M20 Finite difference methods applied to problems in fluid mechanics
65M55 Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs

Software:

Expokit
Full Text: DOI

References:

[1] David, A. P., An exponential method of numerical integration of ordinary differential equations, Commun ACM, 6, 491-493 (1963) · Zbl 0117.11204
[2] Beylkin, G.; Keiser, J. M.; Vozovoi, L., A new class of time discretization schemes for the solution of nonlinear PDEs, J Comput Phys, 147, 362-387 (1998) · Zbl 0924.65089
[3] Hochbruck, M.; Ostermann, A., Explicit exponential runge-kutta methods for semilinear parabolic problems, SIAM J Numer Anal, 43, 1069-1090 (2005) · Zbl 1093.65052
[4] Cox, S. M.; Matthnews, P. C., Exponential time differencing for stiff systems, J Comput Phys, 176, 430-455 (2002) · Zbl 1005.65069
[5] Zhang, X.; Yang, D. H.; Song, G. J., A nearly analytical exponential time difference method for solving 2d seismic wave equation, Earthq Sci, 27, 57-77 (2014)
[6] Li, S. J.; Wang, Z. J.; Ju, L. L.; Luo, L. S., Explicit large time stepping with a second-order exponential time integrator scheme for unsteady and steady flows, 55th AIAA ASM (2018)
[7] Courant, R.; Friedrichs, K.; Lewy, H., On the partial difference equations of mathematical physics, IBM J Res Dev, 11, 215-234 (1928) · Zbl 0145.40402
[9] Cai, X. C.; Sarkis, M., A restricted addictive schwarz preconditioner for general sparse linear systems, SIAM J Sci Comput, 21, 792-797 (1999) · Zbl 0944.65031
[10] Hoang, T. T.; Ju, L. L.; Wang, Z., Overlapping localized exponential time differencing methods for diffusion problems, Comm Math Sci (2017)
[11] Jiang, G. S.; Shu, C. W., Efficient implementation of weighted ENO schemes, J Comput Phys, 126, 202-228 (1996) · Zbl 0877.65065
[12] Saad, Y., Analysis of some krylov subspace approximations to the matrix exponential operator, SIAM J Numer Anal, 29, 209-228 (1992) · Zbl 0749.65030
[13] Roger, B. S., Expokit: a software package for computing matrix exponentials, ACM Trans Math Softw, 24, 130-156 (1998) · Zbl 0917.65063
[14] Wang, Z. J.; Fidkowski, K. J.; Abgrall, R., High-order CFD methods: current status and perspective, Inter J Numer Meths Fluids, 72, 811-845 (2013) · Zbl 1455.76007
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.