Open AccessImplementation of kinetic mechanisms of methane combustion by example of expanding functionality of physicochemical libraries in conjunction with reactingPimpleCentralFoam solver
Author(s) -
D. S. Koov,
V. Yu. Gidaspov,
Sergei Strijhak
Publication year - 2021
Publication title -
trudy instituta sistemnogo programmirovaniâ ran/trudy instituta sistemnogo programmirovaniâ
Language(s) - English
Resource type - Journals
eISSN - 2220-6426
pISSN - 2079-8156
DOI - 10.15514/ispras-2021-33(5)-17
Subject(s) - solver , combustion , methane , shock tube , mach number , mechanics , thermodynamics , shock (circulatory) , shock wave , flow (mathematics) , computational fluid dynamics , viscosity , chemistry , ignition system , mathematics , physics , mathematical optimization , medicine , organic chemistry
The possibility of using reduced combustion mechanisms for hydrocarbon fuels in solvers developed and used at ISP RAS is investigated. These mechanisms contain smaller number of stages and substances appearing in them, but they allow obtaining results in a good agreement with experimental data in a much shorter calculation time. The comparison is made with the results obtained with using the Moscow Aviation Institute solvers. A modified mechanism of methane combustion is considered. It can be extended to describe the chemical reaction processes in other hydrocarbon-oxygen mixtures. The choice of methane is due to the prospects of this fuel at the present time. As the first test problem, a standard chemFoam is used. This solver was designed to demonstrate the occurrence of chemical reactions in a computational domain consists of one cell only. The ignition delay and the parameter values of the thermodynamic equilibrium state reached are taken as the comparison criteria. The second test problem is the flow modeling in a shock tube after the shock wave has been reflected from the wall. This problem is considered in a three-dimensional domain using the ISP RAS reactingPimpleCentralFoam solver. Results were compared with ones obtained by grid-characteristic and Godunov's methods in one-dimensional nonstationary calculations. The effects of viscosity, thermal conductivity and diffusion are not taken into account. The distributions of the flow parameters behind the reflected shock wave are obtained. Results are analyzed depending on the value of the falling shock wave Mach number. Estimates of the possible application of this reduced mechanism are given.