In this article, we study theoretically and numerically the interaction of a vortex induced by a rotating cylinder with a perpendicular plane. We show the existence of weak solutions to the swirling vortex models by using the Hopf extension method, and by an elegant contradiction argument, respectively. We demonstrate numerically that the model could produce phenomena of swirling vortex including boundary layer pumping and two-celled vortex that are observed in potential line vortex interacting with a plane and in a tornado.

Vortex and vorticity are two correlated but fundamentally different concepts which have been the central issues in fluid mechanics research. Vorticity has rigorous mathematical definition (curl of velocity), but no clear physical meaning. Vortex has clear physical meaning (rotation) but no rigorous mathematical definition. For a long time, many people treat them as a same thing. However, based on our high-order direct numerical simulation (DNS), we found that first, “vortex” is not “vorticity tube” or “vortex tube” which is widely defined as a bundle of vorticity lines without any vorticity line leak. Actually, vortex is an open area for vorticity line penetration. Second, vortex is not necessarily congregation of vorticity lines, but dispersion in many 3-dimensional cases. Some textbooks say that vortex cannot end inside the flow field but must end on the solid wall (and/or boundaries). Our DNS observation and many other numerical results show almost all vortices are ended inside the flow field. Finally, a more theoretical study shows that neither vortex nor vorticity line can attach to the solid wall and they must be detached from the wall.

Computational fluid dynamics (CFD) has been used by numerous researchers for the simulation of flows around wind turbines. Since the 2000s, the experiments of NREL phase VI blades for blind comparison have been a de-facto standard for numerical software on the prediction of full scale horizontal axis wind turbines (HAWT) performance. However, the characteristics of vortex structures in the wake, whether for modeling the wake or for understanding the aerodynamic mechanisms inside, are still not thoroughly investigated. In the present study, the flow around N-REL phase VI blades was numerically simulated, and the results of the wake field were compared with the experimental ones of a one-to-eight scaled model in a low-speed wind tunnel. A good agreement between simulation and experimental results was achieved for the evaluation of overall performances. The simulation captured the complete formation procedure of tip vortex structure from the blade. Quantitative analysis showed the streamwise translation movement of vortex cores. Both the initial formation and the damping of vorticity in near wake field were predicted. These numerical results showed good agreements with the measurements. Moreover, wind tunnel wall effects were also investigated on these vortex structures, and it revealed further radial expansion of the helical vortical structures in comparison with the free-stream case.

Dynamics and instability of a vortex ring impinging on a wall were investigated by means of large eddy simulation for two vortex core thicknesses corresponding to thin and thick vortex rings. Various fundamental mechanisms dictating the flow behaviors, such as evolution of vortical structures, formation of vortices wrapping around vortex rings, instability and breakdown of vortex rings, and transition from laminar to turbulent state, have been studied systematically. The evolution of vortical structures is elucidated and the formation of the loop-like and hair-pin vortices wrapping around the vortex rings (called briefly wrapping vortices) is clarified. Analysis of the enstrophy of wrapping vortices and turbulent kinetic energy (TKE) in flow field indicates that the formation and evolution of wrapping vortices are closely associated with the flow transition to turbulent state. It is found that the temporal development of wrapping vortices and the growth rate of axial flow generated around the circumference of the core region for the thin ring are faster than those for the thick ring. The azimuthal instabilities of primary and secondary vortex rings are analyzed and the development of modal energies is investigated to reveal the flow transition to turbulent state. The modal energy decay follows a characteristic –5/3 power law, indicating that the vortical flow has become turbulence. Moreover, it is identified that the TKE with a major contribution of the azimuthal component is mainly distributed in the core region of vortex rings. The results obtained in this study provide physical insight of the mechanisms relevant to the vortical flow evolution from laminar to turbulent state.

In the present work, a new type of coupled compact difference scheme has been proposed for the solution of computational acoustics and flow problems. The proposed scheme evaluates the first, the second and the fourth derivative terms simultaneously. Derived compact difference scheme has a significant spectral resolution and a physical dispersion relation preserving (DRP) ability over a considerable wavenumber range when a fourth order four stage Runge-Kutta scheme is used for the time integration. Central stencil has been used for the present numerical scheme to evaluate spatial derivative terms. Derived scheme has the capability of adding numerical diffusion adaptively to attenuate spurious high wavenumber oscillations responsible for numerical instabilities. The DRP nature of the proposed scheme across a wider wavenumber range provides accurate results for the model wave equations as well as computational acoustic problems. In addition to the attractive feature of adaptive diffusion, present scheme also helps to control spurious reflections from the domain boundaries and is projected as an alternative to the perfectly matched layer (PML) technique.

Interaction of vortex rings with solid is an important research topic of hydrodynamic. In this study, a multiple-relaxation time (MRT) lattice Boltzmann method (LBM) is used to investigate the flow of a vortex ring impacting spheroidal particles. The MRT-LBM is validated through the cases of vortex ring impacting a flat wall. The vortex evolution due to particle size, the aspect ratio of a prolate particle, as well as Reynolds (Re) number are discussed in detail. When the vortex ring impacting a stationary sphere, the primary and secondary vortex rings wrap around each other, which is different from the situation of the vortex ring impacting a plate. For the vortex ring impacting with a prolate spheroid, the secondary vortex ring stretches mainly along the long axis of the ellipsoid particle. However, it is found that after the vortex wrapping stage, the primary vortex recovers along the short axis of the particle faster than that in the long axis, i.e., the primary vortex ring stretches mainly along the short axis of the particle. That has never been address in the literature.

A vortex ring impacting a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re = 4 × 104 based on the initial diameter and translational speed of the vortex ring. The effects of bump height and vortex core thickness for thin and thick vortex rings on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The boundary vorticity flux is analyzed to reveal the mechanism of the vorticity generation on the bump surface. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Further, the analysis of turbulent kinetic energy reveals the transition from laminar to turbulent state. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the flow evolution and the flow transition to turbulent state.

The swimming of a 3D fish-like body with finlets is numerically investigated at Re = 1000 (the Reynolds number is based on the uniform upstream flow and the length of the fish-like body). The finlets are simply modeled as thin rigid rectangular plates that undulate with the body. The wake structures and the flow around the caudal peduncle are studied. The finlets redirect the local flow across the caudal peduncle but the vortical structures in the wake are almost not affected by the finlets. Improvement of hydrodynamic performance has not been found in the simulation based on this simple model. The present numerical result is in agreement with that of the work of Nauen and Lauder [J. Exp. Biol., 204 (2001), pp. 2251-2263] and partially supports the hypothesis of Webb [Bull. Fish. Res. Bd. Can., 190 (1975), pp. 1-159].

Cooling methods are needed for turbine blade tips to ensure a long durability and safe operation. A common way to cool a tip is to use serpentine passages with 180-deg turn under the blade tip-cap taking advantage of the three-dimensional turning effect and impingement like flow. Improved internal convective cooling is therefore required to increase the blade tip lifetime. In the present study, augmented heat transfer of an internal blade tip with pin-fin arrays has been investigated numerically using a conjugate heat transfer method. The computational domain includes the fluid region and the solid pins as well as the tip regions. Turbulent convective heat transfer between the fluid and pins, and heat conduction within pins and tip are simultaneously computed. The main objective of the present study is to observe the effect of the pin material on heat transfer enhancement of the pin-finned tips. It is found that due to the combination of turning, impingement and pin-fin crossflow, the heat transfer coefficient of a pin-finned tip is a factor of 2.9 higher than that of a smooth tip at the cost of an increased pressure drop by less than 10%. The usage of metal pins can reduce the tip temperature effectively and thereby remove the heat load from the tip. Also, it is found that the tip heat transfer is enhanced even by using insulating pins having low thermal conductivity at low Reynolds numbers. The comparisons of overall performances are also included.