A mathematical model and a numerical model for hyperbolic mass transport in compressible flows

Abstract

A number of contributions have been made during the last decades to model pure-diffusive transport problems by using the so-called hyperbolic diffusion equations. These equations are used for both mass and heat transport. The hyperbolic diffusion equations are obtained by substituting the classic constitutive equation (Fick’s and Fourier’s law, respectively), by a more general differential equation, due to Cattaneo (C R Acad Sci Ser I Math 247:431–433, 1958). In some applications the use of a parabolic model for diffusive processes is assumed to be accurate enough in spite of predicting an infinite speed of propagation (Cattaneo, C R Acad Sci Ser I Math 247:431–433, 1958). However, the use of a wave-like equation that predicts a finite velocity of propagation is necessary in many other calculations. The studies of heat or mass transport with finite velocity of propagation have been traditionally limited to pure-diffusive situations. However, the authors have recently proposed a generalization of Cattaneo’s law that can also be used in convective-diffusive problems (Gomez, Technical Report (in Spanish), University of A Coruna, 2003; Gomez et al., in An alternative formulation for the advective-diffusive transport problem. 7th Congress on computational methods in engineering. Lisbon, Portugal, 2004a; Gomez et al., in On the intrinsic instability of the advection–diffusion equation. Proc. of the 4th European congress on computational methods in applied sciences and engineering (CDROM). Jyvaskyla, Finland, 2004b) (see also Christov and Jordan, Phys Rev Lett 94:4301–4304, 2005). This constitutive equation has been applied to engineering problems in the context of mass transport within an incompressible fluid (Gomez et al., Comput Methods Appl Mech Eng, doi:10.1016/j.cma.2006.09.016, 2006). In this paper we extend the model to compressible flow problems. A discontinuous Galerkin method is also proposed to numerically solve the equations. Finally, we present some examples to test out the performance of the numerical and the mathematical model.