scholarly journals Numerical stability analysis of a vortex ring with swirl

2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès
Keyword(s):  
Numerical Simulation ◽  
Growth Rate ◽  
Direct Numerical Simulation ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Base Flow ◽  
Critical Layer ◽  
Navier Stokes ◽  
Instability Mode ◽  
Navier Stokes Equations

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.

scholarly journals Controlling the onset of turbulence by streamwise travelling waves. Part 2. Direct numerical simulation

2010 ◽  
Vol 663 ◽  
pp. 100-119 ◽  
Author(s):  
BINH K. LIEU ◽  
RASHAD MOARREF ◽  
MIHAILO R. JOVANOVIĆ
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Stokes Equations ◽  
Travelling Waves ◽  
Base Flow ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Flow Equations ◽  
Onset Of Turbulence

This study builds on and confirms the theoretical findings of Part 1 of this paper (Moarref & Jovanović, J. Fluid Mech., 2010, doi:10.1017/S0022112010003393). We use direct numerical simulation of the Navier–Stokes equations to assess the efficacy of blowing and suction in the form of streamwise travelling waves for controlling the onset of turbulence in a channel flow. We highlight the effects of the modified base flow on the dynamics of velocity fluctuations and net power balance. Our simulations verify the theoretical predictions of Part 1 that the upstream travelling waves promote turbulence even when the uncontrolled flow stays laminar. On the other hand, the downstream travelling waves with parameters selected in Part 1 are capable of reducing the fluctuations' kinetic energy, thereby maintaining the laminar flow. In flows driven by a fixed pressure gradient, a positive net efficiency as large as 25 % relative to the uncontrolled turbulent flow can be achieved with downstream waves. Furthermore, we show that these waves can also relaminarize fully developed turbulent flows at low Reynolds numbers. We conclude that the theory developed in Part 1 for the linearized flow equations with uncertainty has considerable ability to predict full-scale phenomena.

On Direct Numerical Simulation of Turbulent Flows

2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
High Performance Computing ◽  
Turbulent Flows ◽  
High Performance ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Turbulent Shear ◽  
Navier Stokes Equations ◽  
Performance Computing

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.

Direct Numerical Simulation of Effects of Small Angle of Incidence on Honji Instability

2011 ◽  
Author(s):  
Kun Yang ◽  
Liang Cheng ◽  
Hongwei An ◽  
Ming Zhao
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Oscillatory Flow ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Angle Of Incidence ◽  
Navier Stokes Equations ◽  
Small Incidence ◽  
Inclination Angles

This paper concerns Honji instability generated around a circular cylinder in an oscillatory flow with a small oblique angle. In this study, direct numerical simulation has been conducted for an oscillatory flow past a stationary cylinder with small incidence angles (α) of 5° and 10° at KC number of 2 and β number of 200. The three-dimensional Navier-Stokes equations are solved using the Petrov-Galerkin finite element method. Flow structures around the cylinder are visualized through using streamlines, velocity vectors and vorticity contours. Honji instability has been captured at both chosen inclination angles. However Honji vortex pairs are asymmetric at α = 5° and 10° due to the inclination of the oscillation direction and can only be observed during the flow reversal. It is also found that the flow inclination appears to suppress the three-dimensional instability.

Direct numerical simulation of controlled transition in a flat-plate boundary layer

1995 ◽  
Vol 298 ◽  
pp. 211-248 ◽  
Author(s):  
U. Rist ◽  
H. Fasel
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Flat Plate ◽  
Stokes Equations ◽  
Plate Boundary ◽  
Numerical Data ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Flat Plate Boundary Layer

The three-dimensional development of controlled transition in a flat-plate boundary layer is investigated by direct numerical simulation (DNS) using the complete Navier-Stokes equations. The numerical investigations are based on the so-called spatial model, thus allowing realistic simulations of spatially developing transition phenomena as observed in laboratory experiments. For solving the Navier-Stokes equations, an efficient and accurate numerical method was developed employing fourth-order finite differences in the downstream and wall-normal directions and treating the spanwise direction pseudo-spectrally. The present paper focuses on direct simulations of the wind-tunnel experiments by Kachanov et al. (1984, 1985) of fundamental breakdown in controlled transition. The numerical results agreed very well with the experimental measurements up to the second spike stage, in spite of relatively coarse spanwise resolution. Detailed analysis of the numerical data allowed identification of the essential breakdown mechanisms. In particular, from our numerical data, we could identify the dominant shear layers and vortical structures that are associated with this breakdown process.

scholarly journals Embedded direct numerical simulation for aeronautical CFD

2001 ◽  
Vol 105 (1046) ◽  
pp. 193-198 ◽  
Author(s):  
N. D. Sandham ◽  
M. Alam ◽  
S. Morin
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Euler Equations ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Modified Equations ◽  
Separation Bubbles

Abstract A method is proposed by which a direct numerical simulation of the compressible Navier-Stokes equations may be embedded within a more general aeronautical CFD code. The method may be applied to any code which solves the Euler equations or the Favre-averaged Navier-Stokes equations. A formal decomposition of the flowfield is used to derive modified equations for use with direct numerical simulation solvers. Some preliminary applications for model flows with transitional separation bubbles are given.

scholarly journals Simulation and visualization of formation of vortex ring, its propagation and transportation of passive scalar

2021 ◽  
pp. 182-199
Author(s):  
К.Н. Волков ◽  
В.Н. Емельянов ◽  
И.Е. Капранов
Keyword(s):  
Numerical Simulation ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Passive Scalar ◽  
Navier Stokes ◽  
Time Interval ◽  
Navier Stokes Equations ◽  
Volume Method ◽  
Traditional Approaches ◽  
Similar Theory

Рассматривается численное моделирование газодинамических процессов, сопровождающих формирование и распространение вихревых колец, получаемых при помощи поршневого генератора. Обсуждается влияние характеристик вихревого кольца на перенос пассивной примеси. Для численных расчетов применяются нестационарные уравнения Навье--Стокса, для дисукретизации которых применяется метод конечных объемов. Результаты численного моделирования позволяют получить геометрические и динамические характеристики вихревого кольца, которые соответствуют автомодельному теории вихревого кольца и экспериментальным данным. Помимо традиционных подходов к визуализации вихревых течений, основанных на построении линий уровня различных характеристик потока, для визуализации вихревых структур применяются инварианты тензора градиента скорости и метод показателей Ляпунова на конечном промежутке времени. Numerical simulation of gas-dynamic processes accompanying the formation and propagation of vortex rings obtained using a piston generator is considered. The influence of the characteristics of the vortex ring on the transfer of the passive particles is discussed. Unsteady Navier--Stokes equations are used for numerical calculations, and finite volume method is applied to their discretization. The results of numerical simulation make it possible to obtain the geometric and dynamic characteristics of the vortex ring, which correspond to the self-similar theory of the vortex ring and experimental data. In addition to traditional approaches to visualization of vortex flows based on the construction of contours of various flow quantities, invariants of the velocity gradient tensor and the method of Lyapunov exponents over a finite time interval are used to visualize vortex structures.

SMOOTHED PROFILE METHOD TO SIMULATE COLLOIDAL PARTICLES IN COMPLEX FLUIDS

2009 ◽  
Vol 20 (09) ◽  
pp. 1457-1465 ◽  
Author(s):  
RYOICHI YAMAMOTO ◽  
YASUYA NAKAYAMA ◽  
KANG KIM
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Colloidal Particles ◽  
Stokes Equations ◽  
Electrolyte Solutions ◽  
Rigid Body Dynamics ◽  
Spherical Particles ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Charged Colloids

A new direct numerical simulation scheme, called "Smoothed Profile (SP) method," is presented. The SP method, as a direct numerical simulation of particulate flow, provides a way to couple continuum fluid dynamics with rigid-body dynamics through smoothed profile of colloidal particle. Our formulation includes extensions to colloids in multicomponent solvents such as charged colloids in electrolyte solutions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier–Stokes equations so that the electro-hydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters.

Sign in / Sign up

Export Citation Format

Share Document