A numerical study of flow past a cylinder with cross flow and inline oscillation

Two-dimensional fluid flow around an oscillating circular cylinder is studied numerically at different values of oscillation frequency and amplitude. A novel finite element method which uses discretization along the characteristic line is used for simulation. The solver is coupled to a mesh movement scheme using the Arbitrary Lagrangian-Eulerian (ALE) formulation to account for body motion in the flow field. Two cases of cylinder motion have been studied, cross flow and inline oscillation. In both cases, occurrence of lock on is investigated and the bounds of the lock on region are determined. A comparison of the numerical results with the experimental data indicates that 2D simulation is valid up to Re = 300. Beyond that, 3D effects appear. By using flow visualization, effect of a cylinder oscillation on the flow field and wake pattern has been studied. Also, variation of the mean drag coefficient against the oscillation parameters is discussed. The numerical results are in good agreement with the experimental data available in the literature.     

A numerical investigation of fluid flow over a rotating cylinder with cross flow oscillation

This paper studies a two-dimensional incompressible viscous flow past a rotating cylinder with cross flow oscillation using a finite element method based on the characteristic based split (CBS) algorithm to solve governing equations including full Navier–Stokes and continuity equations. Dynamic unstructured triangular grid is used employing lineal and torsional spring analogy which is coupled with the solver by an Arbitrary Lagrangian–Eulerian (ALE) formulation. After verifying the accuracy of the numerical code, simulations are conducted for the flow past a rotating cylinder with cross flow oscillation at moderate Reynolds numbers of 50, 100, and 200 considering different non-dimensional rotational speeds based on the free-stream velocity in the range 0–2.5, and various oscillating amplitudes and frequencies. Effects of the oscillation and rotation of the cylinder on the vortex shedding both in lock-on and non-lock-on regions, the mean drag and lift coefficients, and the Strouhal number are investigated in detail. It is found that similar to the fixed cylinder beyond a critical non-dimensional rotational speed the vortex shedding is highly suppressed. In addition, by increasing the rotational speed of the cylinder, the lift coefficient increases while decreasing the drag coefficient. However, in the vortex lock-on region both the lift and the drag coefficients increase significantly.     

Large eddy simulation of pulsating flow over a circular cylinder at subcritical Reynolds number

The pulsating cross-flow over a single circular cylinder at the subcritical Reynolds number Re_D = 2580 is studied with the large eddy simulation (LES) technique using the standard Smagorinsky model as well as a dynamic model in which the test filtered quantities are evaluated through a truncated Taylor series expansion. The filtered equations are discretised using the finite volume method in an unstructured, collocated grid arrangement with a second-order accurate method, in both space and time. The predictions are compared against very detailed experiments for mean velocities and Reynolds stresses that were performed in a duct of cross-section 72 mm × 72 mm using the PIV technique. The effects of mesh refinement close to the cylinder as well as of subgrid scale model are also examined. The numerical predictions are in very good agreement with the measurements in terms of mean as well as turbulence quantities. The instantaneous flow patterns of the flow field are examined and the effect of the external flow pulsation on the wake characteristics such as vortex formation length, vortex strength, Strouhal number as well as the lift and drag coefficients is quantified. The vortex formation length is decreased while the mean drag, as well as the rms values of the drag and lift coefficients increase significantly under pulsating flow conditions. The performance of the LES technique is analysed in the light of the wake characteristics.     

Computation of shape and drag of a deforming elastic body under fluid dynamic force

This paper describes numerical simulations for the shape and its drag of an elastic body deforming under the fluid dynamic force. The simulations were carried out by coupling the Navier-Stokes equations and the equations of motion of the elastic body. The equations of motion are formulated for an elastic shell model which is composed of material particles connected with elastic springs and dampers. The relation between deforming elastic body shape in response to the fluid dynamic force and its drag force was investigated under the constraint of constant volume and fixed center of gravity of the elastic body for incompressible and compressible supersonic flows. In these simulations, an initial shape of the elastic body is a circular cylinder and starts deformation under the fluid dynamic force.     

Steady and unsteady flow past a rotating circular cylinder at low Reynolds numbers

Results of calculations of the steady and unsteady flows past a circular cylinder which is rotating with constant angular velocity and translating with constant linear velocity are presented. The motion is assumed to be two-dimensional and to be governed by the Navier-Stokes equations for incompressible fluids. For the unsteady flow, the cylinder is started impulsively from rest and it is found that for low Reynolds numbers the flow approaches a steady state after a large enough time. Detailed results are given for the development of the flow with time for Reynolds numbers 5 and 20 based on the diameter of the cylinder. For comparison purposes the corresponding steady flow problem has been solved. The calculated values of the steady-state lift, drag and moment coefficients from the two methods are found to be in good agreement. Notable, however, are the discrepancies between these results and other recent numerical solutions to the steady-state Navier-Stokes equations. Some unsteady results are also given for the higher Reynolds numbers of 60, 100 and 200. In these cases the flow does not tend to be a steady state but develops a periodic pattern of vortex shedding.     

Time-dependent solutions of the viscous incompressible flow past a circular cylinder by the method of series truncation

Semi-analytic solutions of the Navier-Stokes equations are calculated for two-dimensional, symmetrical, viscous incompressible flow past a circular cylinder. The stream and vorticity functions are expanded in the finite Fourier series and then substituted in the Navier-Stokes equations. This led to a system of coupled parabolic partial differential equations which are solved numerically. More terms of the series are required as Reynolds number increases and the present calculations were terminated at Reynolds number 600 with 60 terms of Fourier series. The results are compared with similar calculations and experimental data for Reynolds numbers 60, 100, 200, 500, 550 and 600. At the termination of the calculations for Reynolds numbers 60 and 100, the separation angle, the wake length, the drag coefficient, and the vorticity distributions around the surface were very close to their steady-state values. A secondary vortex appeared on the surface of the cylinder in the case of Reynolds numbers 500, 550 and 600. The wake length, the drag coefficient and the separation angle differ slightly at a given instant in the case of Reynolds numbers 500, 550 and 600.     

A Coupled Legendre-Laguerre Spectral-Element Method for the Navier-Stokes Equations in Unbounded Domains

A coupled Legendre-Laguerre spectral element method is proposed for the Stokes and Navier-Stokes equations in unbounded domains. The method combines advantages of the high accuracy of the Laguerre-spectral method for unbounded domains and the geometric flexibility of the spectral-element method. Rigorous stability and error analysis for the Stokes problem is carried out. Numerical results indicate that the proposed method is very effective for some realistic flow problems in unbounded domains, such as flows passing a circular cylinder.     

Evolutionary techniques versus swarm intelligences: application in reservoir release optimization

In this paper, a nonlinear reservoir release optimization problem has been solved by using four optimization tools with various combinations of input parameters that are generally used in this research field. A comparison has been made between evolutionary methods [genetic algorithm (GA)] and swarm intelligences [particle swarm optimization (PSO) and artificial bee colony (ABC) optimization] in searching the optimum reservoir release policy. From the historical recorded data, the monthly inflow was categorized into three states: high, medium and low. As a guideline for the decision maker, an optimum release curve was generated for each month showing the release options with a variety of different storage conditions. GA (real and binary), ABC optimization and PSO algorithm have been used as optimization tools with the same formulation and objective function for all the methods. For verification of the models, a simulation is done by using 264 monthly historical inflow data. Different indices such as reliability, vulnerability and resiliency were calculated in order to check the performance and risk analysis purposes. The results show that the most recently developed ABC optimization technique provides the best results in meeting demands, avoiding wastage of water and in handling critical period of low flows.     

Numerical simulation of an oscillating cylinder in a cross-flow at low Reynolds number: Forced and free oscillations

A numerical simulation of the flow past a circular cylinder which is able to oscillate transversely to the incident stream is presented in this paper for a fixed Reynolds number equal to 100. The 2D Navier-Stokes equations are solved by a finite volume method with an industrial CFD code in which a coupling procedure has been implemented in order to obtain the cylinder displacement. A preliminary work is first conducted for a fixed cylinder to check the wake characteristics for Reynolds numbers smaller than 150 in the laminar regime. The Strouhal frequency f_S and the aerodynamic coefficients are thus controlled among other parameters. Simulations are then performed with forced oscillations characterized by the frequency ratio F = f₀/f_S, where f₀ is the forced oscillation frequency, and by the adimensional amplitude A. The wake characteristics are analyzed using the time series of the fluctuating aerodynamic coefficients and their power spectral densities (PSD). The frequency content is then linked to the shape of the phase portraits and to the vortex shedding mode. By choosing interesting couples (A, F), different vortex shedding modes have been observed, which are similar to those of the Williamson-Roshko map. A second batch of simulations involving free vibrations (so-called vortex-induced vibrations or VIV) is finally carried out. Oscillations of the cylinder are now directly induced by the vortex shedding process in the wake and therefore, the time integration of the motion is realized by an explicit staggered algorithm which provides the cylinder displacement according to the aerodynamic charges exerted on the cylinder wall. Amplitude and frequency response of the cylinder are thus investigated over a wide range of reduced velocities to observe the different phenomena at stake. In particular, the vortex shedding modes have also been related to the frequency response observed and our results at Re = 100 show a very good agreement with other studies using different numerical approaches.     

Space-time discretizing optimal DRP schemes for flow and wave propagation problems

The present paper reports constrained optimization of explicit Runge–Kutta (RK) schemes, coupled with optimal upwind compact scheme to achieve dispersion relation preservation (DRP) property for high performance computing. Essential ideas of optimization employed in arriving at the proposed time integration scheme are extension of the earlier work reported in Rajpoot et al. (J Comput Phys 2010;229:3623–51). This is in turn an application of the correct error evolution equation in Sengupta et al. (J Comput Phys 2007;226:1211–8). Resultant DRP scheme demonstrated the idea for explicit spatial central difference schemes. Present work is similar, extending it for near-spectral accuracy compact schemes. Practical utility of the developed method is demonstrated by solution of model problems and for flow problems by solving Navier–Stokes equation, some of which cannot be solved by conventional schemes, as the problem of rotary oscillation of cylinder.Developed method is calibrated with: (i) flow past a circular cylinder performing rotary oscillation at Re = 150 and (ii) flow inside a 2D lid-driven cavity (LDC) at Reynolds numbers of Re = 1000 and Re = 10,000. Quantitative and qualitative comparisons show excellent match for rotary oscillation cylinder cases with the experimental results of Thiria et al. (J Fluid Mech 2006;560:123–47). Results for LDC for Re = 1000 are compared with that in Botella & Peyret (Comp Fluids 1998;27:421–33) and results for Re = 10,000 are compared with recent published ones showing triangular vortex in the core.     

Effect of tube spacing on the vortex shedding characteristics of laminar flow past an inline tube array: A numerical study

The effect of tube spacing on the vortex shedding characteristics and fluctuating forces in an inline cylinder array is studied numerically. The examined Reynolds number is 100 and the flow is laminar. The numerical methodology and the code employed to solve the Navier-Stokes and continuity equations in an unstructured finite volume grid are validated for the case of flow past two tandem cylinders at four spacings. Computations are then performed for a six-row inline tube bank for eight pitch-to-diameter ratios, s, ranging from 2.1 to 4. At the smallest spacing examined (s = 2.1) there are five stagnant and symmetric recirculation zones and weak vortex shedding activity occurs only behind the last cylinder. As s increases, the symmetry of the recirculation zones breaks leading to vortex shedding and this process progressively moves upstream, so that for s = 4 there is clear shedding from every row. For any given spacing, the shedding frequency behind each cylinder is the same. A critical spacing range between 3.0 and 3.6 is identified at which the mean drag as well as the rms lift and drag coefficients for the last three cylinders attain maximum values. Further increase to s = 4 leads to significant decrease in the force statistics and increase in the Strouhal number. It was found that at the critical spacing there is 180° phase difference in the shedding cycle between successive cylinders and the vortices travel a distance twice the tube spacing within one period of shedding.     

Transient state analysis of separated flow around a sphere

Transient state solutions of the Navier-Stokes equations were obtained for incompressible flow around a sphere accelerating from zero initial velocity to its terminal free falling velocity. By assuming rotational symmetry about the axis in the direction of motion, the Navier-Stokes equations and the continuity equation were simplified in terms of vorticity and stream function. The instantaneous acceleration of the falling sphere was calculated by considering the difference between the gravitational force and the drag force in a transient state. A set of implicit finite difference equations was developed. In order to obtain accurate information around the body, an exponential transformation along the radial direction was used to provide finer meshes in the vicinity of the surface of the sphere. The vorticity equation was solved by an alternating direction implicit (ADI) method while the stream function equation was solved by a successive over-relaxation (SOR) method. Simultaneous solutions were obtained. Transient state solutions were compared with steady state solutions for Reynolds numbers up to 300. Separations first occurred at a Reynolds number 20 for steady state flows and at Reynolds numbers 22·46 and 28·24 for transient state flows with terminal Reynolds numbers of 100 and 300, respectively. Separation angles, sizes of separation regions, and drag coefficents were calculated for both steady and unsteady states. Good agreement was obtained with existing experimental data in the steady state.     

A finite volume method to solve the 3D Navier–Stokes equations on unstructured collocated meshes

A new method to solve the Navier–Stokes equations for incompressible viscous flows and the transport of a scalar quantity is proposed. This method is based upon a fractional time step scheme and the finite volume method on unstructured meshes. The governing equations are discretized using a collocated, cell-centered arrangement of velocity and pressure. The solution variables are stored at the cell-circumcenters. Theoretical results and numerical properties of the scheme are provided. Predictions of lid-driven cavity flow, flows past a cylinder and heat transport in a cylinder are performed to validate the method.     

Flow about a spherically blunted body in a supersonic wake

Stationary flows about the nose of a spherically blunted body in a supersonic wake-type, nonuniform oncoming stream are investigated. The flow calculation is carried out using the Navier-Stokes equations which are solved by means of the implicit finite-difference, unidirectional scheme. Oncoming stream nonuniformity, rarefaction and wall cooling effects on the shock layer structure and the distributions of the drag and heat transfer parameters along the body surface are investigated. The solutions of the full Navier-Stokes equations are compared with the results obtained using the reduced Navier-Stokes and Euler equations.     

Calculation of internal viscous flows in axisymmetric ducts at moderate to high reynolds numbers

A formal procedure is presented for parabolizing the Navier-Stokes equations for two dimensional internal flows in artitrary ducts with streamline curvature and streamline divergence. This procedure consists of constructing an orthogonal coordinate system from the potential flow solution and parabolizing the corresponding Navier-Stokes equations using the thin channel approximation. Theoretical arguments based on the fundamental existence theorem and numerical examples are used to demonstrate the validity of this procedure. A new method using conformal mapping based on the Schwartz-Christoffel transformation is presented which obtains the solution of the inverse potential flow problem in a direct manner. A numerical solution algorithm based on the two point box scheme is presented for solving the viscous flow equations. This method is shown to be accurate and stable for flows at moderate to high Reynolds numbers over a wide range of conditions.     

Universal instability modes in internal and external flows

Two prototypical external and internal flows have been studied which display linear temporal instability followed by nonlinear saturation taking the flows to a new equilibrium state. Direct simulation results are obtained using a specific formulation and numerical methods with very high accuracy. These results are analyzed via proper orthogonal decomposition (POD), which reveal similar modes for flow past a circular cylinder and flow inside a lid-driven cavity, indicating universality of such modes. Unlike many other efforts on reduced order modeling via POD, here the emphasis has been on understanding the physical aspect of the flow instability which requires very high accuracy of the simulation. Then, the obtained POD modes are related to the instability modes (in the classical sense of defining the latter) and new generic types of instability modes are identified in the studied external and internal flows. These new modes have been reported for flow past a circular cylinder [Sengupta TK, Singh N, Suman VK. Dynamical system approach to instability of flow past a circular cylinder. J Fluid Mech 2010;656:82–115] which help one in understanding the instability sequence and the relative importance of these modes in the flow evolution starting from an impulsive start. Present comparative study, furthermore, reveals universality of such instability modes by showing their presence for the flow inside a lid-driven cavity as well. Despite seeming dissimilarities between these two flows, similarities between the instability portrait of these two flows suggest universality of such modes. From the equilibrium amplitude of vorticity time-series, we establish the presence of multiple modes and multiple bifurcation sequences for these flows in parameter space. Existing theory due to Landau and Stuart that considers only a single dominant mode and its self-interaction does not explain all these features. We invoke a multi-modal interaction model in the cited reference above, termed as Landau–Stuart–Eckhaus (LSE) equation in recognition of Eckhaus’ work in modifying the classical Stuart–Landau equation. We also show that the new instability modes do not follow either the classical Stuart–Landau or the newly proposed LSE model equations and for this reason we call these as anomalous modes. Two specific classes of anomalous modes are identified and classified in the present work. Empirical expressions for the evolution of these anomalous modes are presented and their unambiguous role during instabilities is discussed.     

Computations of optimal cylinder flow control in weakly conductive fluids

A nonlinear adjoint-based optimal control approach of cylinder wake by electromagnetic force has been investigated numerically in the paper. A cost functional representing the balance of the regulated quantities with different weights and interaction parameter N (Lorentz force) has been constituted, where the regulated quantities related with flow and force are taken as targets of regulation and the Lorentz force, (as interaction parameter N), is taken as a control input. Based on the cost functional and Navier-Stokes equations, the corresponding adjoint equations have been derived and the sensitivity of the cost functional is found to be a simple function of the adjoint stream function in the adjoint field. For the different regulations, the forms of optimal control rules are similar while the adjoint equations are different. The receding-horizon predictive control setting is employed to discuss the optimal control problems. Under the action of optimal N(t), the flow separation is suppressed fully, so that the oscillations of drag and lift are suppressed and the total drag coefficient decreases dramatically. For the different regulations, the control effects have some differences due to the different values of optimal inputs corresponding to the different adjoint flow fields.     

Numerical study of interference between simple-shape bodies in hypersonic flows

Hypersonic rarefied-gas flows near two side-by-side plates and cylinders, toroidal balloon, plate and cylinder over a plane surface, and plate behind a cylinder in argon, nitrogen, oxygen, and carbon dioxide have been studied numerically using the direct simulation Monte-Carlo technique under the transition flow conditions at Knudsen numbers from 0.004 to 10. Strong influences of the geometrical factor (the ratio of a distance between bodies to a body length) and the Knudsen number on the flow structure about the bodies (shock-wave shapes, the configuration of subsonic flow zones), skin friction, pressure distribution, lift, and drag have been found.     

On restraining the production of small scales of motion in a turbulent channel flow

Since most turbulent flows cannot be computed directly from the incompressible Navier-Stokes equations, a dynamically less complex mathematical formulation is sought. In the quest for such a formulation, we consider nonlinear approximations of the convective term that preserve the symmetry and conservation properties. In particularly, the energy, enstrophy (in 2D) and helicity are conserved. The underlying idea is to restrain the convective production of small scales in an unconditional stable manner, meaning that the approximate solution cannot blow up in the energy-norm (in 2D also: enstrophy-norm). The numerical algorithm used to solve the governing equations preserves the symmetry and conservation properties too. The resulting simulation method is successfully tested for a turbulent channel flow (Re_τ = 180 and 395).     

The numerical solution of the Navier-Stokes equations for laminar incompressible flow past a paraboloid of revolution

A numerical methods is presented for the solution of the Navier-Stokes equations for flow past a paraboloid of revolution. This method is based upon the ideas of van de Vooren and collaborators [1,2]. The flow field has been computed for a large range of Reynolds numbers. Results are presented for the skinfriction and the pressure together with their respective drag coefficients. The total drag has been checked by means of an application of the momentum theorem.