Fast solver for systems of axisymmetric ring vortices

Abstract

A method which is capable of efficient calculation of the axisymmetric flow field produced by a large system of ring vortices is presented in this report. The system of ring vortices can, in turn, be used to model body surfaces and wakes in incompressible unsteady axisymmetric flow fields. This method takes advantage of source point and field point series expansions which enables one to make calculations for interactions between groups of vortices which are in well separated spatial domains rather than having to consider interactions between every pair of vortices. In this work, series expansions for the stream function of the ring vortex system are obtained. Such expansions explicitly contain the radial and axial velocity components. A Fortran computer code RSOLV has been written to execute the fast solution technique to calculate the stream function and the axial and radial velocity components at points in the flow field. Test cases have been run to optimize the code and to benchmark the truncation errors and CPU time savings associated with the method. Non-dimensional truncation errors for the stream function and total velocity field are on the order of 5 {times} 10{sup {minus}5} and 3 {times} 10{sup {minus}3} respectively. Single precision accuracymore » produces errors in these quantities up to about 1 {times} 10{sup {minus}5}. For 100 vortices in the field, there is virtually no CPU time savings with the fast solver. For 10,000 vortices in the flow, the fast solver obtains solutions in about 1% to 3% of the time required for the direct solution technique. Simulations of vortices with square and circular cores were run in order to obtain expressions for the self-induced velocities of such vortices. 8 refs., 26 figs.« less