The Velocity of Viscous Vortex Rings

The motion of vortex rings of small cross section is considered. A formula is derived for the velocity of a ring in an ideal fluid with an arbitrary distribution of vorticity in the core. The effects of viscosity are then examined. A definition of the velocity of an unsteady diffusing ring is given and it is shown that the method used to calculate the ring speed in an ideal fluid can be extended to give the velocity of a vortex ring subject to viscous diffusion.     

Numerical simulation of the interaction between vortex ring and bubble plume

In present paper a three-dimensional Vortex-In-Cell method with two-way coupling effect was developed to study the bubble plume entrainment by a vortex ring. In this method the continuous flow was calculated by the three-dimensional Vortex-In-Cell method and the bubbles are tracked through bubble motion equation. Two-way coupling effect between continuous flow and dispersed bubbles is considered by introducing a vorticity source term, which is induced by the change of void fraction gradient in each computational cell. After validated by the comparison between experimental measurements and simulation results for the motion of vortex rings and the rising velocity of bubble plume, present method is implemented to simulate the interaction between an evolving vortex ring and a rising bubble plume. It was found that there is little effect of the bubble entrainment to the total circulation of vortex ring while the effect of bubble entrainment to the vortex ring structure is quite obvious. The bubble entrainment by the vortex ring not only changed the vorticity distribution in the vortex structure, but also displaced the positions of the vortex cores. The vorticity in the lower vortex core of the vortex ring decreases more than that in the upper vortex core of the vortex ring while the vortex core in the upper part of the vortex ring is displaced to the center of vortex ring by the entrained bubbles. Smaller bubbles are easier to be entrained by the large scale vortex structure and the transportation distance is in inverse proportion to bubble diameter.     

流体动力学方程的三维旋涡解的可叠加性

给出了流体动力学方程的轴对称流动的一类精确解,讨论了一些例子并综述了目前已知的基本涡元精确解。发现三维空间内的某些旋涡解可以叠加成仍然满足非线性方程的新的精确解。由此可以用来分析旋涡的产生、演化和相互作用。此外还讨论了旋涡解的对称性。     

Analysis of Sound Radiated by Axisymmetric Vortex Ring

A theoretical and numerical study using the unsteady, 3D Navier-Stokes equations to generate axisymmetric vortex rings is reported. Increasing the vorticity, the vortex ring transition to a turbulent state are analyzed. After transition to a turbulent stage, the self-similarity of the temporal evolution of the flow is observed. Then we can compare sound radiated by the vortex ring to jet noise, at similar Reynolds number and low Mach number. The agreement between the simulation results and the simplified model is good. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

LES and analyses on the vortex structure behind supersonic MVG with turbulent inflow

In this paper, an implicitly implemented high order large eddy simulation by using the fifth order bandwidth-optimized WENO scheme is applied to make comprehensive studies on ramp flows with and without control at Mach 2.5 and Re_θ = 5760. Flow control in the form of microramp vortex generators (MVG) is applied. The mechanism of vortex ring generation behind MVG has been studied in detail and shear layer instability has been studied and found as the mechanism of K–H vortex ring generation. A series of new observations on the flow around supersonic MVG have been made including inflection points (surface in 3-D), vorticity conservation, interaction of the primary vortex and new generated K–H vortex rings, and the K–H vortex ring structure. The numerical observations have been confirmed by the experimental work.     

Existence and numerical modelling of vortex rings with elliptic boundaries

We revisit in this paper the theory of axisymmetric vortex rings in an ideal fluid. The boundary separating the vortex ring from the external (potential) flow is assumed of elliptic shape. For a given distribution of vorticity in the vortex core, we theoretically put into evidence the critical parameter for the existence of non-trivial solutions, thus confirming the numerical observation of Durst et al. [ZAMP 32 (1981) 156]. A sharp estimation of the critical threshold is analytically derived. Theoretical predictions are confirmed by numerical simulations using finite elements. A new numerical algorithm is presented and shown to display better performances compared to previous published algorithms using finite differences. The convergence of the iterative algorithm is proved using the theory of elliptic partial differential equations with discontinuous nonlinearities.     

Vortex Interaction in a Ducted Flow

The objective of this study is to describe the structure of pipe flow by considering it as a superposition of many axisymmetric vortex rings. In knowing the unsteady gross feature of pipe flow, the investigation on vortex interactions is important. As a first step to the goal, we investigate the nonlinear interaction among vortex rings whose number is three at most. The interaction among vortex rings of equal circulation is here investigated. Momentum and energy conservation of the present vortex ring system are also discussed to know a better understanding of the perturbed pipe flow. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

激波在异种气体中传播及诱导的剪切混合研究

利用二阶迎风TVD格式求解多组分,层流全N－S方程,针对直通道和突扩直通道,研究了马赫数为2和4的激波在H2和空气界面上的传播及诱导的燃料剪切混合,计算结果表明：（1）直通道中,剪切层中的激波阵面要发生畸变,存在对混合起主要作用的卷吸涡,激波马赫数不同,卷吸涡结构和横向混合的尺寸也不同,激波马赫数低,剪切混合效果好,（2）在突扩直通道中,马赫数为2和4的激波在H2中产生不同强度激波,在剪切层中都能产生顺时针,尺度较大的卷吸涡,后台阶增强了剪切层的混合。     

Hydrodynamics of a particle impact on a wall

The problem of a particle impacting on a wall, a common phenomenon in particle-laden flows in the minerals and process industries, is investigated computationally using a spectral-element method with the grid adjusting to the movement of the particle towards the wall. Remeshing is required at regular intervals to avoid problems associated with mesh distortion and the constantly reducing maximum time-step associated with integration of the non-linear convective terms of the Navier–Stokes equations. Accurate interpolation between meshes is achieved using the same high-order interpolation employed by the spectral-element flow solver. This approach allows the full flow evolution to be followed from the initial approach, through impact and afterwards as the flow relaxes. The method is applied to the generic two-dimensional and three-dimensional bluff body geometries, the circular cylinder and the sphere. The principal case reported here is that of a particle colliding normally with a wall and sticking. For the circular cylinder, non-normal collisions are also considered. The impacts are studied for moderate Reynolds numbers up to approximately 1200. A cylindrical body impacting on a wall produces two vortices from its wake that convect away from the cylinder along the wall before stalling while lifting induced wall vorticity into the main flow. The situation for a sphere impact is similar, except in this case a vortex ring is formed from the wake vorticity. Again, secondary vorticity from the wall and particle plays a role. At higher Reynolds number, the secondary vorticity tends to form a semi-annular structure encircling the primary vortex core. At even higher Reynolds numbers, the secondary annular structure fragments into semi-discrete structures, which again encircle and orbit the primary core. Vorticity fields and passive tracer particles are used to characterize the interaction of the vortical structures. The evolution of the pressure and viscous drag coefficients during a collision are provided for a typical sphere impact. For a Reynolds number greater than approximately 1000 for a sphere and 400 for a cylinder, the primary vortex core produced by the impacting body undergoes a short-wavelength instability in the azimuthal/spanwise direction. Experimental visualisation using dye carried out in water is presented to validate the predictions.     

Three-dimensional sand ripples as the product of vortex instability

Three-dimensional sand ripples can be observed under steady liquid flows in both nature and industry. Some examples are the ripples observed on the bed of rivers and in petroleum pipelines conveying sand. Although of importance, the formation of these patterns is not completely understood. There are theoretical and experimental evidence that aquatic ripples grow from two-dimensional bed instabilities, so that a straight vortex is formed just downstream of their crests. The proposition of Raudkivi (2006) [18], that three-dimensionality has its origin in a vortex instability, is employed here. This paper presents a linear stability analysis of the downstream vortex in order to obtain the transverse scales of three-dimensional ripples. The obtained wavelength is compared with experimentally observed ripples.     

The influence of turbulence on a columnar vortex with axial flow

The interaction between a columnar vortex and external turbulence is investigated numerically. A q -vortex is immersed in an initially isotropic homogeneous turbulence field, which itself is produced numerically by a direct numerical simulation of decaying turbulence. The formation of turbulent eddies around the columnar vortex and the vortex-core deformations are studied in detail by visualizing the flow field. In the less-stable case with q = –1.5, small thin spiral structures are formed inside the vortex core. In the unstable case with q = –0.45, the linear unstable modes grow until the columnar vortex make one turn. Its growth rate agrees with that of the linear analysis of Mayer and Powell[1]. After two turns of the vortex, the secondary instability is excited, which causes collapse of the columnar q -vortex and the sudden appearance of many fine scale vortices. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lyapunov orbits in the n-vortex problem

In the reduced phase space by rotation, we prove the existence of periodic orbits of the n-vortex problem emanating from a relative equilibrium formed by n unit vortices at the vertices of a regular polygon, both in the plane and at a fixed latitude when the ideal fluid moves on the surface of a sphere. In the case of a plane we also prove the existence of such periodic orbits in the (n + 1)-vortex problem, where an additional central vortex of intensity κ is added to the ring of the polygonal configuration.     

Numerical simulation of vortex induced vibrations and its control by suction and blowing

The vortex formation and shedding behind bluff structures is influenced by fluid flow parameters such as, Reynolds number, surface roughness, turbulence level, etc. and structural parameters such as, mass ratio, frequency ratio, damping ratio, etc. When a structure is flexibly mounted, the Kármán vortex street formed behind the structure gives rise to vortex induced oscillations. The control of these flow induced vibrations is of paramount practical importance for a wide range of designs. An analysis of flow patterns behind these structures would enable better understanding of wake properties and their control. In the present study, flow past a smooth circular cylinder is numerically simulated by coupling the mass, momentum conservation equations along with a dynamical evolution equation for the structure. An active flow control strategy based on zero net mass injection is designed and implemented to assess its efficacy. A three actuator system in the form of suction and blowing slots are positioned on the cylinder surface. A single blowing slot is located on the leeward side of the cylinder, while two suction slots are positioned at an angle α = 100°. This system is found to effectively annihilate the vortex induced oscillations, when the quantum of actuations is about three times the free stream velocity. The dynamic adaptability of the proposed control strategy and its ability to suppress vortex induced oscillations is verified. The exact quantum of actuation involved in wake control is achieved by integrating a control equation to decide the actuator response in the form of a closed loop feed back system. Simulations are extended to high Reynolds number flows by employing eddy viscosity based turbulence models. The three actuator system is found to effectively suppress vortex induced oscillations.     

Numerical study on the dynamics of a jet-vortex interaction

The interaction of a vortex and a jet plays an important role for many industrial processes such as carbon black formation or combustion in diesel engines. The knowledge of physics of these phenomena is crucial for engineers, but also for scientists who wish to study this complex flow pattern. In this research, we numerically investigate cases where a cylindrical reactor is charged with gas injected through five ports. The first one is located along the main axis of the cylinder and this leads to the formation of the main jet. The other four ports are situated along the side walls such that the gas entering the cylinder tangentially through them causes a vortex to be formed. The objective of this paper is to show the fundamental physical phenomena and also how the initial and boundary conditions influence the results. It is shown how the process is influenced not only by modifying the velocity of injected gas, but also by the orientation and position of the vortex inlet nozzles.     

The Zero Divisor Graph of 2 × 2 Matrices Over a Field

A zero divisor graph, Γ(R), is formed from a ring R by having each element of Z(R) \ {0} to be a vertex in the graph and having two vertices u and v adjacent if the corresponding elements from the ring are nonequal and have product equal to zero. In this paper, the structure of the zero-divisor graph of 2 × 2 matrices over a field, Γ(M2(F)), are completely determined.     

On the stability of discrete tripole,quadrupole, Thomson’ vortex triangle and square in a two-layer/homogeneous rotating fluid

A two-layer quasigeostrophic model is considered in the f-plane approximation. The stability of a discrete axisymmetric vortex structure is analyzed for the case when the structure consists of a central vortex of arbitrary intensity Γ and two/three identical peripheral vortices. The identical vortices, each having a unit intensity, are uniformly distributed over a circle of radius R in a single layer. The central vortex lies either in the same or in another layer. The problem has three parameters (R, Γ, α), where α is the difference between layer thicknesses. A limiting case of a homogeneous fluid is also considered.A limiting case of a homogeneous fluid is also considered.The theory of stability of steady-state motions of dynamic systems with a continuous symmetry group G is applied. The two definitions of stability used in the study are Routh stability and G-stability. The Routh stability is the stability of a one-parameter orbit of a steady-state rotation of a vortex multipole, and the G-stability is the stability of a three-parameter invariant set O _G , formed by the orbits of a continuous family of steady-state rotations of a multipole. The problem of Routh stability is reduced to the problem of stability of a family of equilibria of a Hamiltonian system. The quadratic part of the Hamiltonian and the eigenvalues of the linearization matrix are studied analytically.The cases of zero total intensity of a tripole and a quadrupole are studied separately. Also, the Routh stability of a Thomson vortex triangle and square was proved at all possible values of problem parameters. The results of theoretical analysis are sustained by numerical calculations of vortex trajectories.     

Numerical discovery and experimental confirmation of vortex ring generation by microramp vortex generator

An implicitly implemented large eddy simulation (ILES), by using the modified fifth order WENO scheme, is applied to study the flow around the microramp vortex generator (MVG) at Mach 2.5 and Re_θ = 5760. A series of new discoveries on the flow around supersonic MVG have been made by the UTA LES team including source of the momentum deficit, inflection points (surface in 3-D), Kelvin–Helmholtz instability and vortex ring generation. Most of the new discoveries, which were made by the UTA LES team and presented in 2009, were confirmed by experiment conducted by the UTA experiment team in 2010. A new 5-pair-vortex-tube model near the MVG is given based on the ILES observation.     

Vortex damping of sloshing in tanks with baffles

The problem of damping the sloshing in tanks with sharp-edged baffles (thin inserts which partially span a longitudinal or transverse cross-section) is considered. Separation of the boundary layer and the formation of vortices occur at these sharp edges. It is assumed that the domains where there is significant vortex motion of the fluid are localized in small neighbourhoods of the sharp edges of the baffles. The non-linear vortex damping is determined from the distribution of the velocity intensity factors at these sharp edges in the same way as the linear damping, caused by the dissipation of energy in a boundary layer close to a wall, is determined from the fluid velocity distribution on the walls of a cavity. Both of the above-mentioned distributions are calculated by solving the same boundary-value problem on the oscillations of an ideal fluid. The second of the distributions characterizes the singular properties of the solutions of this problem on particular lines. A method based on the variation of the area of the baffles, which simplifies the calculation of the velocity intensity factors is described. The distinctive features arising when the method of finite elements is used are considered. The results of numerical calculations of the damping of sloshing in a cylindrical tank with a ring baffle are compared with experimental data.     

Uniform Flow Past an Axially Uniform Viscous Vortex

It is believed that the flow past a tornado causes the formation of smaller vortices which produce the “suction spots” observed along the path of destruction. Here we develop a greatly simplified mathematical model to investigate this phenomenon. An axially uniform vortex is developed by visualizing a circular tube with uniform surface suction of fluid possessing circulation at infinity. This vortex is then perturbed by a uniform flow past it. An inner asymptotic expansion of an E^1/3 radial boundary layer is matched to an outer expansion to obtain a solution. The results show that a stagnation point developing into a secondary vortex is formed in a free shear layer at critical flow conditions. However, it is difficult to apply our results quantitatively because of the difficulty in comparing the axially uniform vortex with a real tornado vortex.