Abstract: With ideas of vector splitting method and technique of coordinate rotation, a D1Q4 compressible lattice Boltzmann model is used to construct convection flux solvers of 3D Navier-Stokes equations with chemical reactions. Three-dimensional chemical non-equilibrium flow is solved by finite volume method of Navier-Stokes equations. A time-splitting method is applied to resolve stiff problem in computation of chemical reaction flows. Three classic examples of supersonic chemical non-equilibrium flows are simulated. It shows that:At high Mach numbers, numerical simulation using three-dimensional convection flux solvers based on D1Q4 model has no non-physical solution, and can distinguish physical phenomena such as shock, combustion wave in chemical non-equilibrium flow field. Accuracy and resolution are high,which verifys reliability of the three-dimensional convection flux solver. It broadens application of D1Q4 model, and offers a new method for calculation of hypersonic chemical non-equilibrium flows.

Key words: Lattice Bohzmann method, D1Q4 model, 3D Navier-Stokes equations, chemical non-equilibrium flow