Original articlesInstability and acoustic fields of the Rankine vortex as seen from long-term calculations with the tenth-order multioperators-based scheme
Introduction
The multioperators-based prescribed-order schemes for fluid dynamics proposed in [20], [21] and described in [22], [23], [25], [26] provide a good opportunity to perform long-term calculations with high resolution of small details of flows described by the Euler or the Navier–Stokes equations. In particular, they fit neatly into the requirements of DNS, LES and computational aeroacoustics.
Using multioperators-based schemes, some shear layers instability problems were numerically simulated. The direct numerical simulation of 2D turbulence was reported in [16]. The unsteady behavior of hot subsonic jets generating vortex rings with self-sustained oscillations was investigated in [15]. Very clear pictures of screech waves generated by slightly supersonic underexpanded jets were presented in [27]. The previous experience has shown that the multioperators technique allows one to describe fine flow details using comparatively modest numbers of grid points.
In the present study, we continue to use very high order and high resolution schemes for investigating sound generation due to flow instabilities. We consider now the development of the instability of the Rankine vortex (that is the vortex with constant vorticity inside a circle and the zero one outside it) in the framework of the Euler equations. It is well known that the flow field of the vortex is the exact steady state solutions of both incompressible and compressible Euler equations. However the solution is not stable in the latter case and the Rankine vortex is an oscillatory system emitting sound. There are numerous publications considering its instability, the natural method of attack being the linear theory (see [2], [3], [9]). The Rankine vortex has drawn considerable attention in the context of sound scattering problems [1], [4], [6], [7], [9], [10], [11] in the theory of sound interaction with low Mach number flows.
In the early work [9], the vorticity boundary was initially () disturbed in the form where are the polar coordinates, is the vortex radius, is a small parameter while the integer is the mode number. It was shown by Kelvin that in the incompressible case one has for the stable solution where is a constant vorticity inside the circle. There is no sound radiation in that case. However the solution is unstable if the compressible Euler equations are considered.
In [2], [3], [11], a general form of the compressible Rankine vortex perturbations was used to describe the disturbed flow fields. The resulting eigenvalue problem solved numerically in [2], [3] produced modified frequencies and the increments of the unstable th modes for various Mach numbers.
It was found in particular [9] that the increment of exponential growth of the most unstable mode is where is the Mach number defined by being the ambient sound speed. The eigenvalue problems for the linearized Euler equations were solved numerically in [17] for various vorticity distributions, flow field calculations based on the linear theory being carried out.
Obviously, the linear theory can describe the behavior of unstable modes up to relatively small values of time . So the question arise what happens with the Rankine vortex (in general, with an isolated vortex) described by the Euler equations for large values of . The answer to the question is relevant to real life cases of viscous flows once the characteristic time of the vortex diffusion is considerably greater than .
There are certain difficulties when attempting to numerically simulate the behavior of unstable isolated vortices during large time intervals. First of all, the amount of harmonics resolved for example by a finite difference scheme with mesh size is limited by the relation where is the wave number. Thus refined meshes are required to meet that condition for describing the impact of high-order modes on the instability development. Unfortunately, the condition is not sufficient for getting accurate solutions for large due to accumulation with time of phase and (possibly) amplitude errors introduced by any scheme. Thus only schemes with small errors of that type can successfully accomplish it, the smaller being the errors, the greater being the admissible values.
An attempt to perform direct numerical simulations of the flow generated by the disturbed vortex with the boundary given by Eq. (1) for was reported in [28]. The second-order scheme [8] was applied to the Euler equations, the initial data being based on the well known exact solution for the incompressible elliptic vortex. It was found that the obtained flow field fluctuations at least for 100 rotations periods are characterized by the frequencies and the growth increments which agree well with the value from [2] and the value defined by Eq. (3) respectively. In-depth study of various isolated vortices presented in [3] contains also the results of a numerical simulation of the instability. The calculations based on the Euler equations for the local Mach number equal to 1.5 were aimed at clarifying the nonlinear stage of the flow field development. It was found that the initially disturbed by the second eigenmode vortex splits after some time moments into two corotating vortices.
Both calculations [3], [28] were limited to relatively modest time intervals (less than 1000 dimensionless time units ). They do not reveal the details of the modes excitations. To the best of our knowledge, the complete description of the vortex behavior and sound radiation up to considerably larger values of is lacking in the relevant literature.
The high-order multioperators-based schemes provide a good opportunity of long time direct numerical simulations in the area of computational aeroacoustic and turbulence modelings. It is due to their potential for producing very small phase and amplitude errors not only in the range of long and medium waves harmonics, but in the range of the shortest waves supported by meshes as well.
The main aim of the present study is to investigate the scenario of the unstable behavior of the subsonic Rankine vortex and its sound radiation during large time intervals using the advantages of high resolution dispersion-preserving multioperators-based schemes. The Rankine vortex was chosen because there are in-depth studies concerning its linear instability thus allowing to compare numerical and analytical data.
Both numerical method and the obtained results are presented in the following Sections.
In Section 2, the multioperators-based scheme used in the calculations is briefly outlined, the emphasize being placed on its optimization. The numerical data describing the vortex instability development and the resulting sound radiation are presented in Section 2. Finally, Section 3 the summarizes the findings and main conclusions. Appendix A Details of the multioperators, Appendix B Optimization procedures contain some details of the multioperators and the optimization procedure.
Section snippets
Multioperators-based scheme
To perform the long-time calculations resolving small-scale flow details of the Rankine vortex, the tenth-order multioperator and the ninth-order dissipation multioperator were used to approximate inviscid terms of the Navier–Stokes equations. They are particular realizations of the approximation to a target linear operator looking as where is a th order compact approximation to defined on a uniform mesh with mesh size while
Numerical results
The analysis of the obtained numerical solutions allows to cover the following topics:
-
describing the instability onset and development seen when using the present highly accurate scheme;
-
comparisons with linear theory;
-
estimates of sound radiations for different Mach numbers;
-
estimates of viscosity influence;
-
comparisons with the solutions obtained with low-order schemes.
Before presenting them, it should be emphasized that the linear theory of the Rankine vortex
Summary and discussion
A version of multioperators-based schemes for fluid dynamics with inversions of two-point operators is presented. It differs from the schemes described previously (see for example [25], [26]) in that it does not use upwind–downwind pairs of multioperators. Instead, a central multioperators are constructed for both convection and dissipation terms. In this way, higher orders can be achieved for the same numbers of basis operators. Four values of the parameter and hence four basis operators were
Acknowledgments
The study was performed by a grant from the Russian Science Foundation (project No. 14-11-00775). Funding was provided through Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences. The authors would like to thank Professor Victor F. Kopiev for drawing attention to the problem, providing helpful information as well as for valuable comments during fruitful discussions.
References (28)
On the scattering of sound by a rectilinear vortex
J. Sound Vib.
(1999)- et al.
On long-wave sound scattering by a rankine vortex: non-resonant and resonant cases
J. Sound Vib.
(2010) - et al.
Dispersion-relation-preserving finite difference schemes for computational acoustics
J. Comput. Phys.
(1993) Multioperator high-order compact upwind methods for cfd parallel calculations
On multioperators principle for constructing arbitrary-order difference schemes
Appl. Numer. Math.
(2003)Development of arbitrary-order multioperators-based schemes for parallel calculations. 1. Higher-than-fifth-order approximations to convection terms
J. Comput. Phys.
(2007)Development of arbitrary-order multioperators-based schemes for parallel calculations. 2. Families of compact approximations with two-diagonal inversions and related multioperators
J. Comput. Phys.
(2008)- et al.
Fast calculations of screech using highly accurate multioperators-based schemes
Appl. Acoust.
(2013) - et al.
On the statement of the problem of sound scattering by a cylindrical vortex
Acoust. Phys.
(2008) - et al.
Acoustic destabilization of a vortices
Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci.
(1979)
Instabilities of two-dimensional inviscid compressible vortices
J. Fluid Mech.
The scattering of sound waves by a vortex: numerical simulations and analytical solutions
J. Fluid Mech.
Accuracy evaluation of unsteady CFD numerical schemes by vortex preservation
Comput. & Fluids
Low-frequency acoustic scattering from a trailing vortex
J. Acoust. Soc. Amer.
Cited by (8)
General scenario and fine details of compressible Gaussian vortex unforced instability
2021, European Journal of Mechanics, B/FluidsCitation Excerpt :It is because the wave length corresponding the second mode frequency is greater than the sizes of the computational domains used in the calculations. Judging from the calculations for the Rankine vortex [6] and the linear theory [4], the second mode frequency at the final stage is expected to be not far from a linear function of the maximal swirl velocity. Moreover, it was found in [6] that it is approximately equal to the Mach number corresponding to that velocity which is about 0.06 in the present case.
Parallel Implementation of the 16th-Order Multioperator Scheme: Application to Problems of the Instability of Vortices and Boundary Layers
2023, Mathematical Models and Computer SimulationsExcitation and Development of Instability in a Compressible Boundary Layer as Obtained in High-Order Accurate Numerical Simulation without Introducing Artificial Perturbations
2022, Computational Mathematics and Mathematical PhysicsUsing 16-th Order Multioperators-Based Scheme for Supercomputer Simulation of the Initial Stage of Laminar-Turbulent Transitions
2021, Communications in Computer and Information ScienceOn 32nd-Order Multioperators-Based Schemes for Fluid Dynamics
2021, Smart Innovation, Systems and Technologies