zbMATH — the first resource for mathematics

Automatic differentiation using operator overloading (ADOO) for implicit resolution of hyperbolic single phase and two-phase flow models. (English) Zbl 1453.65049
Summary: Implicit time integration schemes are widely used in computational fluid dynamics to speed-up computations. Indeed, implicit schemes usually allow for less stringent time-step stability constraints than their explicit counterpart. The derivation of an implicit scheme is however a challenging and time-consuming task, increasing substantially with the model equations complexity since this method usually requires fairly accurate evaluation of the spatial scheme’s matrix Jacobian. This article presents a flexible method to overcome the difficulties associated to the computation of the derivatives, based on the forward mode of automatic differentiation using operator overloading (ADOO). Flexibility and simplicity of the method are illustrated through implicit resolution of various flow models of increasing complexity such as the compressible Euler equations, a two-phase flow model in full equilibrium [S. Le Martelot et al., “Towards the direct numerical simulation of nucleate boiling flows”, Int. J. Multiphase Flow 66, 62–78 (2014; doi:10.1016/j.ijmultiphaseflow.2014.06.010)] and a symmetric variant [the second author et al., J. Fluid Mech. 495, 283–321 (2003; Zbl 1080.76062)] of the two-phase flow model of M. R. Baer and J. W. Nunziato [Int. J. Multiphase Flow 12, 861–889 (1986; Zbl 0609.76114)] dealing with mixtures in total disequilibrium.
65D25 Numerical differentiation
65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs
76M12 Finite volume methods applied to problems in fluid mechanics
65-04 Software, source code, etc. for problems pertaining to numerical analysis
76T10 Liquid-gas two-phase flows, bubbly flows
Full Text: DOI
[1] Abgrall, R.; Saurel, R., Discrete equations for physical and numerical compressible multiphase mixtures, J. Comput. Phys., 361-396 (2003) · Zbl 1072.76594
[2] Baer, M. R.; Nunziato, J. W., A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials, Int. J. Multiph. Flow, 861-889 (1986) · Zbl 0609.76114
[3] Balay, S., PETSc Users Manual (2018), Argonne National Laboratory
[4] Barth, T.; Jespersen, D., The Design and Application of Upwind Schemes on Unstructured Meshes (1989), AIAA
[5] Barth, T. J.; Linton, S. W., An Unstructured Mesh Newton Solver for Compressible Fluid Flow and Its Parallel Implementation (1995), AIAA Paper 95-0221
[6] Bilanceri, M., Linearized implicit time advancing and defect correction applied to sediment transport simulations, Comput. Fluids, 63, 82-104 (2012) · Zbl 1365.76138
[7] Bischof, C.; Khademi, P.; Mauer, A.; Carle, A., Adifor 2.0: automatic differentiation of Fortran 77 programs, IEEE Comput. Sci. Eng., 18-32 (1996)
[8] Bramkamp, F. D.; Bücker, H. M.; Rasch, A., Using exact Jacobians in an implicit Newton-Krylov method, Comput. Fluids, 35, 1063-1073 (2006) · Zbl 1177.76213
[9] Briley, W. R.; McDonald, H., Solution of the multi-dimensional compressible Navier-Stokes equations by a generalized implicit method, J. Comput. Phys., 372-397 (1977) · Zbl 0363.76018
[10] Brown, P. N.; Saad, Y., Hybrid Krylov methods for nonlinear systems of equations, SIAM J. Sci. Stat. Comput., 450-481 (1990) · Zbl 0708.65049
[11] Chiapolino, A.; Boivin, P.; Saurel, R., A simple phase transition relaxation solver for liquid-vapor flows, Int. J. Numer. Methods Fluids, 83, 7, 583-605 (2016)
[12] Colonia, S.; Steijl, R.; Barakos, G. N., Implicit implementation of the AUSMP and AUSMP-up schemes, Int. J. Numer. Methods Fluids, 687-712 (2014)
[13] Davis, S. F., Simplified second-order Godunov-type methods, SIAM J. Sci. Stat. Comput., 9, 3, 445-473 (1988) · Zbl 0645.65050
[14] Dennis, J. E.; Schnabel, R. B., Numerical Methods for Unconstrained Optimization and Nonlinear Equations (1983), Prentice-Hall: Prentice-Hall Englewood Cliffs · Zbl 0579.65058
[15] Dubuc, L., Solution of the unsteady Euler equations using an implicit dual time method, AIAA J., 36, 1417-1424 (1998)
[16] Espath, L. F.R.; Linn, R. V.; Awruch, A. M., Shape optimization of shell structures based on NURBS, Int. J. Numer. Methods Eng., 88, 613-636 (2011) · Zbl 1242.74067
[17] Furfaro, D.; Saurel, R., A simple HLLC-type Riemann solver for compressible non-equilibrium two-phase flows, Comput. Fluids, 159-178 (2015) · Zbl 1410.76288
[18] Godunov, S., Finite Difference Method for Numerical Computation of Discontinuous Solutions of the Equations of Fluid Dynamics, Matematicheskii Sbornik, vol. 47, 271-306 (1959), Steklov Mathematical Institute of Russian Academy of Sciences · Zbl 0171.46204
[19] Gottlieb, S.; Shu, C.-W.; Tadmor, E., Strong stability-preserving high-order time discretization methods, SIAM Rev., 89-112 (2001) · Zbl 0967.65098
[20] Griewank, A.; Walther, A., Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation (2000), SIAM · Zbl 0958.65028
[21] Harten, A.; Lax, P. D.; van Leer, B., On upstreaming differencing and Godunov-type, SIAM Rev., 25, 1, 35-61 (1983) · Zbl 0565.65051
[22] Hascoet, L.; Pascual, V., The Tapenade automatic differentiation tool: principles, model, and specification, ACM Trans. Math. Softw. (2013) · Zbl 1295.65026
[23] Hovland, P. D.; McInnes, L. C., Parallel simulation of compressible flow using automatic differentiation and PETSc, Parallel Comput., 27, 503-519 (2001) · Zbl 0972.68165
[24] Jameson, A., Time Dependent Calculations Using Multigrid, with Applications to Unsteady Flows Past Airfoils and Wings (1991), AIAA Paper 91-1596
[25] Jameson, A.; Turkel, E., Implicit scheme and LU-decompositions, Math. Comput., 385-397 (1981) · Zbl 0533.65060
[26] Karypis, G.; Kumar, V., A fast and high quality multilevel scheme for partitioning irregular graphs, SIAM J. Sci. Comput., 20, 1, 359-392 (1998) · Zbl 0915.68129
[27] Kennedy, C. A.; Carpenter, M. H., Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review (2016), NASA Langley Research Center: NASA Langley Research Center Hampton, VA, United States, NASA/TM-2016-219173, L-20470, L-20597, NF1676L-19716
[28] Le Martelot, S.; Saurel, R.; Nkonga, B., Towards the direct numerical simulation of nucleate boiling flows, Int. J. Multiph. Flow, 62-78 (2014)
[29] Liou, M. S., A sequel to AUSM: AUSMP, J. Comput. Phys., 129, 2, 364-382 (1996) · Zbl 0870.76049
[30] Liska, R.; Wendroff, B., Comparison of several difference schemes on 1D and 2D test problems for the Euler equations, SIAM J. Sci. Comput., 995-1017 (2004) · Zbl 1096.65089
[31] Liu, X.; Xia, Y.; Luo, H.; Xuan, L., A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations, Commun. Comput. Phys., 1016-1044 (2016) · Zbl 1373.76096
[32] Martin, M. P.; Candler, G. V., A parallel implicit method for the direct numerical simulation of wall-bounded compressible turbulence, J. Comput. Phys., 153-171 (2006) · Zbl 1088.76021
[33] Mohammadi, B.; Pironneau, O., Shape optimization in fluid mechanics, Annu. Rev. Fluid Mech., 36, 255-279 (2004) · Zbl 1076.76020
[34] Mulder, W. A.; Van Leer, B., Implicit Upwind Methods for the Euler Equations (1983), AIAA Paper 83-1930
[35] Newman, J. C., Overview of sensitivity analysis and shape optimization for complex aerodynamic configurations, J. Aircr., 36, 1, 87-96 (1999)
[36] Pandolfi, M.; Larocca, F., Transonic flow about a circular cylinder, Comput. Fluids, 17, 1, 205-220 (1989)
[37] Park, J. S.; Yoon, S.-H.; Kim, C., Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids, J. Comput. Phys., 788-812 (2010) · Zbl 1185.65150
[38] Pulliam, T. H., Time Accuracy and the Use of Implicit Methods (1993), AIAA Paper 93-3360
[39] Reynolds, R.; Samtaney, R.; Tiedeman, H. C., A fully implicit Newton-Krylov-Schwarz method for tokamak magnetohydrodynamics: Jacobian construction and preconditioner formulation, Comput. Sci. Discov., 5, 1 (2012)
[40] Rinaldi, E.; Pecnik, R.; Colonna, P., Exact Jacobians for implicit Navier-Stokes simulations of equilibrium real gas flows, J. Comput. Phys., 459-477 (2014) · Zbl 1349.76381
[41] Rusanov, V. V., Calculation of interaction of non-steady shock waves with obstacles, USSR Comput. Math. Math. Phys., 1, 2, 304-320 (1962)
[42] Saad, Y.; Schultz, M. H., GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 856-869 (1986) · Zbl 0599.65018
[43] Salas, M. D., Recent developments in transonic Euler flow over a circular cylinder, Math. Comput. Simul., 25, 3, 232-236 (1983)
[44] Saurel, R.; Gavrilyuk, S.; Renaud, F., A multiphase model with internal degrees of freedom: application to shock-bubble interaction, J. Fluid Mech., 283-321 (2003) · Zbl 1080.76062
[45] Saurel, R.; Boivin, P.; Le Métayer, O., A general formulation for cavitating, boiling and evaporating flows, Comput. Fluids, 128, 53-64 (2006) · Zbl 1390.76879
[46] Saurel, R.; Fraysse, F.; Furfaro, D.; Lapébie, E., Multiscale multiphase modeling of detonations in condensed energetic materials, Comput. Fluids, 95-111 (2017) · Zbl 1390.76283
[47] Saurel, R.; Le Métayer, O., A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation, J. Fluid Mech., 431, 239-271 (2001) · Zbl 1039.76069
[48] Saurel, R.; Petitpas, F.; Berry, R. A., Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures, J. Comput. Phys., 1678-1712 (2009) · Zbl 1409.76105
[49] Sod, A. S., A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, J. Comput. Phys., 27, 1, 1-31 (1978) · Zbl 0387.76063
[50] Squire, W.; Trapp, G., Using complex variables to estimate derivatives of real functions, SIAM Rev., 40, 1, 110-112 (1998) · Zbl 0913.65014
[51] Toro, E. F.; Spruce, M.; Speares, W., Restoration of the contact surface in the HLL-Riemann solver, Shock Waves, 4, 1, 25-34 (1994) · Zbl 0811.76053
[52] Utke, J.; Naumann, U.; Fagan, M.; Tallent, N.; Strout, M.; Heimbach, P.; Wunsch, C., OpenAD/F: A Modular Open-Source Tool for Automatic Differentiation of Fortran Codes, ACM Trans. Math. Softw., 34, 1-36 (2008) · Zbl 1291.65140
[53] Venkatakrishnan, V.; Barth, T., Application of Direct Solvers to Unstructured Meshes for the Euler and Navier-Stokes Equations Using Upwind Schemes (1989), AIAA Paper 89-0364
[54] Walther, A.; Griewank, A., Getting started with ADOL-C, (Naumann, U.; Schenk, O., Combinatorial Scientific Computing (2012), Chapman-Hall CRC Computational Science), 181-202
[55] Wengert, R. E., A simple automatic derivative evaluation program, Commun. ACM, 463-464 (1964) · Zbl 0131.34602
[56] Wood, A., A Textbook of Sound (1930), G. Bell and Sons Ltd.: G. Bell and Sons Ltd. London
[57] Xia, Y.; Luo, H.; Frisbey, M.; Nourgaliev, R., A set of parallel, implicit methods for a reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids, Comput. Fluids, 98, 134-151 (2014) · Zbl 1391.76370
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.