The study considers the solution of a number of problems of supersonic and hypersonic gas dynamics using a model that takes into account the dissociation and ionization of air. The results of verification and validation of the developed numerical method using various difference schemes (the Roe scheme, Rusanov scheme, AUSM scheme) for discretizing convective flows are presented. The formulation of the mathematical model for high-temperature air uses the presence of equilibrium chemical reactions of dissociation and ionization. For this purpose, at high incoming flow velocities, the Kraiko model is applied, which includes equilibrium chemical reactions in air at high temperatures. To discretize the basic equations, the finite volume method on an unstructured grid is applied. One of the features of the constructed mathematical model is the implementation of the transition between physical and conservative variables. Relationships are given, with the help of which the transition from conservative variables to physical ones and vice versa is carried out when using the high-temperature air model. To ensure the stability of numerical calculations, an entropy correction is introduced. The decrease in entropy in the solution of hyperbolic equations is excluded by introducing an artificial viscosity according to Neumann, as well as by using the Godunov method with an exact solution of the Riemann problem and methods based on the approximate solution of the problem of the decay of an arbitrary discontinuity. A number of problems of supersonic gas dynamics (supersonic flow in a channel with a straight step and supersonic flow around a sphere) are numerically solved taking into account high-temperature effects. The criteria for the accuracy of numerical calculations related to the location of shock-wave structures are discussed. The calculated shock-wave structure of the flow is compared with the data available in the literature, as well as with calculations using the perfect gas model. Some results of numerical calculations are compared with the available experimental data. The shock-wave flow patterns obtained in the framework of the inviscid model, which takes into account the effect of viscosity and its dependence on temperature, and the turbulent flow model are compared. On the basis of numerical simulation data, the influence of viscous effects on the flow characteristics in a channel with a straight step and hypersonic flow around a sphere is considered. The influence of various numerical factors on the shape of the bow shock and the presence of fluctuations in the solution behind the shock is emphasized. As part of the work, a computational module was prepared for the commercial package Ansys Fluent, implemented with the help of user programming tools. The prepared module expands the standard capabilities of commercial software focused on solving computational gas dynamics problems, and is available to Ansys Fluent users for solving hypersonic aerodynamics problems. The developed means of numerical simulation can be useful in the design and optimization of hypersonic aircraft.