In this work a hybrid scheme based on the PISO-algorithm and Kurganov-Tadmor's numerical scheme is proposed. This scheme utilizes compressible PISO method for coupling between velocity and pressure and Kurganov-Tadmor scheme for formulation of non-oscillating convective fluxes. Compressible and incompressible regimes of developed model are switched with blending function depending on local Mach and CFL number. A numerical scheme is implemented by means of OpenFOAM ver. 2.3.0 as pisoCentralFoam independent solver. Investigation of the mathematical model was conducted and exemplified by test cases. Proposed scheme can be used for wide range of Mach numbers from 0.01 to 3 or higher The mesh convergence was analyzed. Comparison of the results with the experimental and analytic data was carried out. The solver was tested in a parallel mode on a computer cluster

A large number of available computational fluid dynamics codes include tools for analysis of computational aeroacoustics problems. Such tools are proprietary as well as the codes themselves. Actually, the level of development of such codes as OpenFOAM makes it possible to implement enough opportunities for complication of physical models and increasing the scale of the issues described. In our paper, we develop a dynamic library libAcoustics which may be compiled independently of any modules of the main OpenFOAM package and the type of solvers being used in the model. The implemented Curle's analogy in library makes it possible to obtain acoustic spectra under the conditions of turbulent flows around arbitrary solids in a medium moving at a low velocity of flow. Calculation of the acoustic field was made for 3D test case of Cylinder – NACA 0012 Wing Profile configuration. The analogy allows user to define solvers settings through standard user I/O dictionaries of the OpenFOAM. A complete implementation of the analogy is capable of producing parallel computation. The libAcoustics library is free and is available for download on demand

Results of QGD (quasi-gas dynamics) equations approximation using OpenFOAM Finite Volume Method library are presented. The new developed solver QGDFoam compared to standard OpenFOAM solver rhoCentralFoam. Two approaches for the approximation of QGD terms are considered and compared - (1) least square method and (2) face normal derivative. Qualitative analysis showed that the first approach is more accurate compared with the second one and is preferable for reproducing complex shock interactions during unsteady flow. It was shown that the second approach produces the solution, similar to rhoCentralFoam solver with upwind interpolation. Being more diffusive than rhoCentralFoam solver with TVD schemes, QGDFoam, allows, however, to reproduce flows at wider Mach number range without introducing too much numerical diffusion comparing to upwind interpolation. Parallel scalability of developed solver has been checked. It scales for meshes down to approximately 1000 cells per node.