Abstract—In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, compressible and rotational and irrotational flow. The algorithm is capable of calculating the duct wall geometries from prescribed wall velocity distributions. The equations modeling the flow are expressed using the stream function (x, y) and the function (x, y) as independent variables, where for irrotational flow (x, y) can be recognized as the velocity potential function, for rotational flow (x, y) ceases being the velocity potential function but does remain orthogonal to the stream lines, the x and y are the usual axial and radial coordinates in cylindrical coordinates respectively. The technique described is capable of tackling the so-called inverse problem where the velocity wall distributions are prescribed from which the duct geometry is calculated, as well as the direct problem where the velocity distribution on the pressure and suction surfaces are calculated from prescribed geometries. The two different cases outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively with results for the Neumann boundary condition only included. The axial velocity and the swirl velocity are prescribed such that vorticity is transported through the duct. The governing second order partial differential is coupled with a set of quasi-linear hyperbolic first order partial differential equations with characteristics parallel to the and axes, the numerical solution is thus obtained iteratively using finite differences to approximate derivatives.

Index Terms—Rotational, compressible flow, Hagen-Poiseuille flow, swirl velocity.

V. Pavlika is with University of London, UK (e-mail: vp4@soas.ac.uk).

 

Cite: V. Pavlika, "The Design of Axisymmetric Ducts for Compressible Flow with Vorticity," International Journal of Applied Physics and Mathematics vol. 4, no. 2, pp. 81-88, 2014.