Method solving hydrodynamics equations using velocity variables

An effective numerical method for solving Navier-Stokes equations and those for porous body models in velocity variables is suggested. The nonexplicit difference analog of the solvable system of equations is obtained by its integration on the net with separated velocity and pressure components and using an approximation down the flow for convective members. To solve it numerically, a trial run method is used that is characterized by an exact solution of a continuity equation at each calculation step. The suggested method is absolutely stable, has a high convergence rate as compared to the methods in velocity-pressure variables and may be used for calculating turbulent flows and flows in the regions with heterogeneous properties. The above method is realized as DUPT program unit for BESM-6 and ES-1060 computers, and is aimed to calculate liquid flow in multi-coupled regions with varied properties, for example, in nuclear reactor coolant circuits.