A modified discrete-vortex method algorithm with shedding criterion for aerodynamic coefficients prediction at high angle of attack

Low-order methods require less computing power than classical computational fluid dynamics and can be implemented on a laptop computer, which is needed for engineering tasks. Discrete vortex methods are such low order methods that can describe the unsteady separated flow around an airfoil. After a presentation of the leading edge suction parameter discrete vortex method, a modified algorithm is proposed, in order to reduce the computing cost, and compared with the previous one. Several reference unsteady airfoil motions are discussed in terms of gain in the computation time with comparisons between the previous scheme and the present one. The accuracy of the new method is demonstrated through aerodynamic coefficients. The application of the present discrete vortex method to a transient pitching motion of an airfoil is also presented, in order to understand the leading edge vortex formation, and its implication in terms of lift and drag coefficients. The method is not limited to unsteady or transient motions but can also simulate the flow around a constant angle of attack airfoil. In that case, an original method of fast summation of the vortices located far away from the airfoil, allows a linear dependence of the computation time versus the number of vortices shed, which is a great improvement over the quadratic dependence observed in the classical discrete vortex methods. The development of the aerodynamic coefficients with angle of attack, from values ranging between −10° and 90°, is obtained for a purely two-dimensional flow. In particular, the shape of the lift coefficient of the airfoil in the fully detached flow region is established. Comparisons with relevant experimental or computational fluid dynamics data are discussed in order to grasp the influence of upstream turbulence level and three-dimensional effects in the measured data in the fully detached flow region.     

Turbulent Flow with Embedded Vortical Structures Induced by Vortex Generators in a Cascade

Presented work is the next step after several experimental examinations of vortex generator influence on a flow separation occurring on a model of the NACA 63A421 airfoil with deflected simple flap. In this stage of research the vortices produced by vortex generators (VGs) were studied using Particle Image Velocimetry technique (PIV) and numerical simulations. Vane type VGs with two spacings among VGs pairs in straight channel with turbulent flow were tested. The average velocity flow field, peak of vorticity and circulation decay downstream of VGs were evaluated. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stationary Configurations of Point Vortices on a Cylinder

In this paper we study the problem of constructing and classifying stationary equilibria of point vortices on a cylindrical surface. Introducing polynomials with roots at vortex positions, we derive an ordinary differential equation satisfied by the polynomials. We prove that this equation can be used to find any stationary configuration. The multivortex systems containing point vortices with circulation Γ₁ and Γ₂ (Γ₂ = ?μΓ₁) are considered in detail. All stationary configurations with the number of point vortices less than five are constructed. Several theorems on existence of polynomial solutions of the ordinary differential equation under consideration are proved. The values of the parameters of the mathematical model for which there exists an infinite number of nonequivalent vortex configurations on a cylindrical surface are found. New point vortex configurations are obtained.     

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise charac-     

Relative equilibrium and collapse configurations of four point vortices

Relative equilibrium configurations of point vortices in the plane can be related to a system of polynomial equations in the vortex positions and circulations. For systems of four vortices the solution set to this system is proved to be finite, so long as a number of polynomial expressions in the vortex circulations are nonzero, and the number of relative equilibrium configurations is thereby shown to have an upper bound of 56. A sharper upper bound is found for the special case of vanishing total circulation. The polynomial system is simple enough to allow the complete set of relative equilibrium configurations to be found numerically when the circulations are chosen appropriately. Collapse configurations of four vortices are also considered; while finiteness is not proved, the approach provides an effective computational method that yields all configurations with a given ratio of velocity to position.      

On the problem of falling motion of a circular cylinder and a vortex pair in a perfect fluid

We consider a system consisting of a heavy circular cylinder in the field of gravity interacting dynamically with a vortex pair in a perfect fluid. The circulation about the cylinder is assumed to be zero. It is shown that, unlike the famous Föppl configuration, the vortices cannot be in a relative equilibrium. An asymptotic system and a suitable regularization are considered.     

Transient cavitating flow structure and acoustic analysis of a hydrofoil with whalelike wavy leading edge

In this paper, the open source software OpenFOAM is used to perform a numerical investigation of the cavitating flow around a modified NACA63₄-021 hydrofoil with bioinspired, wavy leading edge, with particular emphasis on study of the interactions between the cavitation and the streamwise vortices and the far-field radiation noise. A modified k-ω shear-stress transport (SST) model coupled with the Schnerr-Sauer cavitation model and the Ffowcs Williams-Hawkings (FW-H) acoustic analogy approach are introduced to the simulation. The transient cavitation structure and the streamwise vortices are captured well and the results show significant interactions between the cavitation and the streamwise vortices. Cavitation can promote fragmentation of the streamwise vortices, while the streamwise vortices cause cavitation inception to occur earlier and bind the cavitation within the trough region by changing the pressure distribution on the hydrofoil. The transformation of the FW-H equation's solution indicates that the monopole noise is directly related to the cavitation volume acceleration and the dipole noise is related to the mechanical force of the hydrofoil on fluids and the rate at which this force changes. The collapse of cavitation cloud and the collision of the re-entrant jet and main flow will cause violent fluctuation of the mechanical force thus produce instantaneous extreme dipole noise values, while the monopole noise is relatively strong in the cavitation collapse stage due to significant cavitation volume acceleration. The time domain hydroacoustic characteristics are similar for the modified hydrofoil and the baseline hydrofoil.     

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode’“dual vortex pair mode” downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f _m in addition to the vortex shedding frequency f _s is identified. Frequency coupling between f _m and f _s occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.     

The critical cross-section of a vortex

Summary This paper discusses the problem of critical-flow cross-sections in vortex flows. It is shown that there are two different types of vortex flows, A-type and B-type vortices (say). An A-type vortex approaches its critical flow state as its cross-sectional area increases and departs from the critical state as the cross-sectional area is decreased. This property is associated with the particular dependence of total pressure and circulation on the stream function, and it holds for both subcritical and supercritical A-type vortices. On the other hand, both subcritical and supercritical B-type vortices approach their critical flow states as their cross-sectional areas are decreased and depart from their critical states for increasing cross-sectional area. As was shown by Benjamin, setting the first variation of the flow force with respect to the stream function equal to zero leads to Euler's equation of motion. The second variation also vanishes if the corresponding flow state is critical. In this case the sign of the third variation decides whether the flow is an A-type or a B-type vortex. Within the framework of inviscid-fluid flow theory the type of a vortex is preserved unless vortex breakdown occurs. Making use of the knowledge that vortex flows are controlled by two different types of critical-flow cross-sections a variety of vortex flow phenomena are investigated, including the two types of inlet vortices that are observed upstream of jet engines, the behavior of vortex valves, the flow characteristics of liquid-fuel atomizers and the bath tub vortex.     

翼涡干扰前缘开孔被动控制数值研究

开孔方法是一种简单的流动被动控制方法.为找到一种有效降低桨涡干扰效应的被动控制方法,以NACA 0012翼型作为研究对象,建立了4种前缘开孔的模型.在不同来流速度、涡的强度和干扰距离条件下,对4种前缘开孔模型和无孔的基准翼型进行了二维平行桨涡干扰（翼涡干扰）数值模拟,对比了升力系数的变化.结果表明:前缘开孔可以降低翼涡干扰效应,但对翼型升力系数有一定的影响;宽度为2.5%弦长的直孔能在翼型升力系数损失较小的情况下有效地降低翼涡干扰效应,且适用范围较广.     

Regularization by Noise for the Point Vortex Model of mSQG Equations

We consider the point vortex model associated to the modified Surface Quasi-Geostrophic(mSQG) equations on the two dimensional torus. It is known that this model is well posed for almost every initial conditions. We show that, when the system is perturbed by a certain space-dependent noise, it admits a unique global solution for any initial configuration. We also present an explicit example for the deterministic system on the plane where three different point vortices collapse.     

Dynamics and self-propulsion of a spherical body shedding coaxial vortex rings in an ideal fluid

We describe a model for the dynamic interaction of a sphere with uniform density and a system of coaxial circular vortex rings in an ideal fluid of equal density. At regular intervals in time, a constraint is imposed that requires the velocity of the fluid relative to the sphere to have no component transverse to a particular circular contour on the sphere. In order to enforce this constraint, new vortex rings are introduced in a manner that conserves the total momentum in the system. This models the shedding of rings from a sharp physical ridge on the sphere coincident with the circular contour. If the position of the contour is fixed on the sphere, vortex shedding is a source of drag. If the position of the contour varies periodically, propulsive rings may be shed in a manner that mimics the locomotion of certain jellyfish. We present simulations representing both cases.     

Relative Equilibria of the (1+N)-Vortex Problem

We examine existence and stability of relative equilibria of the n-vortex problem specialized to the case where N vortices have small and equal circulation and one vortex has large circulation. As the small circulation tends to zero, the weak vortices tend to a circle centered on the strong vortex. A special potential function of this limiting problem can be used to characterize orbits and stability. Whenever a critical point of this function is nondegenerate, we prove that the orbit can be continued via the Implicit Function Theorem, and its linear stability is determined by the eigenvalues of the Hessian matrix of the potential. For N≥3 there are at least three distinct families of critical points associated to the limiting problem. Assuming nondegeneracy, one of these families continues to a linearly stable class of relative equilibria with small and large circulation of the same sign. This class becomes unstable as the small circulation passes through zero and changes sign. Another family of critical points which is always nondegenerate continues to a configuration with small vortices arranged in an N-gon about the strong central vortex. This class of relative equilibria is linearly unstable regardless of the sign of the small circulation when N≥4. Numerical results suggest that the third family of critical points of the limiting problem also continues to a linearly unstable class of solutions of the full problem independent of the sign of the small circulation. Thus there is evidence that linearly stable relative equilibria exist when the large and small circulation strengths are of the same sign, but that no such solutions exist when they have opposite signs. The results of this paper are in contrast to those of the analogous celestial mechanics problem, for which the N-gon is the only relative equilibrium for N sufficiently large, and is linearly stable if and only if N≥7.     

The Motion of a Vortex Pair in a Stratified Atmosphere

The movement of a horizontal vortex pair through an inhomogeneous fluid is considered. The problem is formulated first for the case when the ambient fluid is uniform, the fluid moving with the vortex pair has a different density, and the motion is supposed laminar and inviscid. An approximate solution is obtained, which predicts that the distance between the vortices stays constant and the vortices accelerate at a constant rate. This solution is then applied to motion in a stratified atmosphere and it is found that the vortices oscillate vertically with a frequency and amplitude depending on the initial conditions and the stratification. Finally, approximate equations are constructed to describe the effects of turbulent entrainment into the fluid moving with the vortex pair, and an estimate of the damping is obtained.     

Finite-time Collapse of Three Point Vortices in the Plane

We investigate the finite-time collapse of three point vortices in the plane utilizing the geometric formulation of three-vortexmotion from Krishnamurthy, Aref and Stremler (2018) Phys. Rev. Fluids 3, 024702. In this approach, the vortex system is described in terms of the interior angles of the triangle joining the vortices, the circle that circumscribes that triangle, and the orientation of the triangle. Symmetries in the governing geometric equations of motion for the general three-vortex problem allow us to consider a reduced parameter space in the relative vortex strengths. The well-known conditions for three-vortex collapse are reproduced in this formulation, and we show that these conditions are necessary and sufficient for the vortex motion to consist of collapsing or expanding self-similar motion. The geometric formulation enables a new perspective on the details of this motion. Relationships are determined between the interior angles of the triangle, the vortex strength ratios, the (finite) system energy, the time of collapse, and the distance traveled by the configuration prior to collapse. Several illustrative examples of both collapsing and expanding motion are given.     

Wakes and vortex streets behind a localized force: Numerical simulations

Two-dimensional wakes behind a body force acting inside a small circular area are investigated using direct numerical simulations. The flows induced by a single force are asymptotically related to the far-field wakes of a bluff body but belong to a wider class of flows because the problem contains an extra control parameter. Stable (almost parallel) wakes as well as regular vortex streets similar to those observed in the wakes of bluff bodies were obtained in our simulations. The behavior of the frequency of vortex shedding in the unstable wakes is described in detail for different values of the main control parameters of the flow, namely the amplitude of the forcing, the velocity of the stream, the size of the forcing area and the kinematic viscosity of fluid. Two different regimes of vortex shedding were observed in the space of these control parameters. Transition between the regimes is characterized by a rapid drop in frequency. The relation between the flows generated by a force and those past a circular cylinder is identified.     

Simulation of Aerodynamic Performances of Flexible Flapping Wing Airfoils北大核心CSCD

与固定翼相比,在低速、小Reynolds数条件下,扑翼飞行具有显著的气动性能优势,受到越来越多的重视。然而,目前对扑翼翼型的研究以刚性翼型为主,对柔性翼型气动性能认识还不清楚。该文建立了柔性椭圆翼型的流固耦合仿真模型,分析了不同风速、迎角下柔性椭圆翼型的周围流场、变形以及气动性能。仿真结果表明,较刚性翼型,柔性翼型延缓了尾涡脱落时间,有效降低升力扰动振荡频率;柔性翼型显著抑制了尾流流场的扰动,降低升力扰动振荡幅值,合适的弹性模量翼型使得扰动振荡完全消除。研究结果可为软飞行器气动设计提供参考。     

Invariant tori in perturbed three vortex motion

The motion of three point vortices in an ideal fluid in a plane comprises a Hamiltonian dynamical system – one that is completely integrable, so it exhibits numerous periodic orbits, and quasiperiodic orbits on invariant tori. Certain perturbations of three vortex dynamics, such as three vortex motion in a half-plane, are also Hamiltonian, but not completely integrable. Yet these perturbed systems may still have periodic trajectories and invariant tori close to those for the unperturbed dynamics. Extending recent work by the authors [4], invariant 2-tori approximating those for the unperturbed system are located and analyzed using a combination of classical analysis, asymptotics, and Hamiltonian methods. It is shown that the results and approximation methods used are applicable to several perturbations of three vortex dynamics such as three vortices in a half-plane, the restricted four vortex problem in the plane, and three coaxial vortex rings in 3-space. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hydrodynamic Vortex on Surfaces

The equations of motion for a system of point vortices on an oriented Riemannian surface of finite topological type are presented. The equations are obtained from a Green’s function on the surface. The uniqueness of the Green’s function is established under hydrodynamic conditions at the surface’s boundaries and ends. The hydrodynamic force on a point vortex is computed using a new weak formulation of Euler’s equation adapted to the point vortex context. An analogy between the hydrodynamic force on a massive point vortex and the electromagnetic force on a massive electric charge is presented as well as the equations of motion for massive vortices. Any noncompact Riemann surface admits a unique Riemannian metric such that a single vortex in the surface does not move (“Steady Vortex Metric”). Some examples of surfaces with steady vortex metric isometrically embedded in \(\mathbb {R}^3\) are presented.     

Homostrophic vortex interaction under external strain,in a coupled QG-SQG model

The interaction between two co-rotating vortices, embedded in a steady external strain field, is studied in a coupled Quasi-Geostrophic — Surface Quasi-Geostrophic (hereafter referred to as QG-SQG) model. One vortex is an anomaly of surface density, and the other is an anomaly of internal potential vorticity. The equilibria of singular point vortices and their stability are presented first. The number and form of the equilibria are determined as a function of two parameters: the external strain rate and the vertical separation between the vortices. A curve is determined analytically which separates the domain of existence of one saddle-point, and that of one neutral point and two saddle-points. Then, a Contour-Advective Semi-Lagrangian (hereafter referred to as CASL) numerical model of the coupled QG-SQG equations is used to simulate the time-evolution of a sphere of uniform potential vorticity, with radius R at depth −2H interacting with a disk of uniform density anomaly, with radius R, at the surface. In the absence of external strain, distant vortices co-rotate, while closer vortices align vertically, either completely or partially (depending on their initial distance). With strain, a fourth regime appears in which vortices are strongly elongated and drift away from their common center, irreversibly. An analysis of the vertical tilt and of the horizontal deformation of the internal vortex in the regimes of partial or complete alignment is used to quantify the three-dimensional deformation of the internal vortex in time. A similar analysis is performed to understand the deformation of the surface vortex.