Published online by Cambridge University Press: 12 April 2006
A semi-analytical numerical study was performed to simulate the development of a vortex ring in a stratified and/or shearing environment. Practical applications of results of this study can be found in turbulence modelling and in studies of plumes and wakes. The objective is to follow exactly the evolution of a vortex ring so that the three-dimensional vortex-stretching mechanisms due to stratification and the shear effects, respectively, can be understood.
The basic formulation consists of the solution of the vorticity equation in a stratified medium. The approach adopted is unique in that discrete vortex elements are used and arbitrary nonlinear interactions are allowed (therefore three-dimensional effects) among various vorticity generators. One of the two fundamental assumptions in this approach is that the vorticity is allowed to be generated only along the density discontinuity. The second assumption is that, while the vorticity carried by the vortex ring is modelled by vortex elements tangential to the vortex loop (which was a vortex ring initially), the vorticity generated by stratification effects is modelled by long vortex lines parallel to the axis of the vortex ring. This limits the validity of the present calculation to high Froude number flow.
Numerical stability is guaranteed by the finite core radius for each discrete vortex element and uniform spacing between them; the former is determined by consideration of the momentum integral over the vortex-ring plane. The latter is determined by a cubic spline interpolation method which conserves the circulation and centroids of the vorticity. The velocity of each vortex element is determined by the discretized Biot-Savart law, and motion of the vortex loop is calculated by a predicator-corrector time integration method.
Calculations were carried out for both momentum-carrying and momentumless vortex rings. A particular two-dimensional case gives good agreement with Kármán's theory. The evolution of the vortex loop reveals a process in which only the vorticity normal to the stratification is conserved; the remaining vorticity is dissipated through a simulated viscous dissipation. Evolution of a vortex loop on a shear layer reveals a vortex-loop rotation rate equal to the velocity shear, and a twisting motion due to the Magnus force which can lead to the turbulence energy cascade phenomenon. Numerical results demonstrate effects of each individual vorticity source and observed phenomena can be explained.