Multiresolution Schemes for the Reactive Euler Equations

• Published in: Journal of computational physics Vol. 154; no. 1; pp. 197 - 230
• Main Authors: Bihari, Barna L, Schwendeman, Donald
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 01.09.1999
• Subjects: stiff source terms; multiresolution scheme; conservation laws; essentially nonoscillatory interpolation; reaction problems
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1006/jcph.1999.6312

We present multiresolution (MR) schemes for the efficient numerical solution of the one-dimensional system of the reactive Euler equations, which has possibly stiff source terms. The original version of the method was developed by A. Harten (1995, Comm. Pure Appl. Math.48(12), 1305) for homogeneous hyperbolic conservation laws. By computing the cell average MR-representation of the solution, we obtain much information about the solution's regularity. This description of smoothness is then used to reduce the number of direct flux computations as well as the expensive high-order ENO (essentially nonoscillatory) reconstruction both of which are now performed only near discontinuities. Thereby, the numerical solution procedure becomes considerably more efficient. In the present case of the reactive Euler equations, the average efficiency factor measured by counting the number of actual flux computations ranges from about 5 to 12. This is on the same order of, and in some cases comes reasonably close to, actual speed-up factors obtained by code timings, which were between 3 to 5. The MR overhead rate was about 10% for the ENO and 36% for TVD schemes, respectively. The quality of the solution is shown to be the same as that of the finest grid. Detailed numerical and performance results are shown for up to fourth-order accuracy, for source terms ranging from moderate to extremely stiff.