Flow modelling in a straight hard-walled duct with two rectangular axisymmetric narrowings

A method for modelling the flow in a rigid-walled duct with two narrowings has been developed. It has the second order of accuracy in the spatial and the first order of accuracy in the temporal coordinates, provides high stability of the solution, and compared to the similar methods requires much less computational time to obtain a result. According to the method, the stream function and the vorticity are introduced initially, and consequently the transition from the governing equations, as well as the initial and boundary conditions to the proper relationships for the introduced variables is performed. The obtained relationships are rewritten in a non-dimensional form. After that a computational domain and a uniform computational mesh are chosen, and the corresponding discretization of the non-dimensional relationships is performed. Finally, the linear algebraic equations obtained as a result of the discretization are solved.