Simulations of axisymmetric, inviscid swirling flows in circular pipes with various geometries

The numerical simulations of the dynamics of high Reynolds number ( Re>100,000 ) swirling flows in pipes with varying geometries of engineering applications continues to be a challenging computational problem, specifically when vortex-breakdown zones or wall-separation regions naturally evolve in the flows. To tackle this challenge, the present paper describes a simulation scheme of the evolution of inviscid, axisymmetric and incompressible swirling flows in expanding or contracting pipes. The integration in time of the circulation together with azimuthal vorticity uses an explicit, first-order accurate finite-difference scheme with a second-order accurate upwind difference formulation in the axial and radial directions. The Poisson solver for advancing in time the spatial distribution of the stream function as a function of azimuthal vorticity uses a second-order accurate over-relaxation difference scheme. No additional numerical steps are needed for computing the natural evolution of flows including the dynamics to states with slow-speed recirculation zones along the pipe centerline or attached to pipe wall. This numerical method shows convergence of the computed results with mesh refinement for various swirl levels and pipe geometry variations. The computed results of time-asymptotic states also present agreement with available theoretical predictions of steady vortex flows in diverging or contracting pipes. In addition, comparison with available experimental data demonstrates that the present algorithm accurately predicts instability processes and long-term mean-flow dynamics of vortex flows in pipes at high Re . The inviscid-flow simulation results support the theoretical predictions and clarify the nature of high- Re flow evolution.