Active decoupling of the axisymmetric body wake response to a pitching motion

Controlled interactions between fluidic actuators and the cross flow over the aft end of a wire-mounted axisymmetric wind tunnel bluff body model (Re_D=2.3·10⁵) are exploited for modification of the near wake dynamics, and the consequent global aerodynamic loads. Actuation is effected using an array of four aft-facing synthetic jet modules through narrow, azimuthally-segmented slots that are equally distributed around the perimeter of the tail end. The model is supported by eight wires, each including a miniature inline force transducer for measurements of the time-resolved tension. The model’s position is varied in a prescribed trajectory by synchronous activation of shape memory alloy (SMA) segments in each of the mounting wires, and the aerodynamic forces and moments are manipulated over a range of pitch attitude. The effectiveness of the flow control approach is demonstrated by decoupling of the wake response from the body’s pitch motion at a low pitch frequency (k=0.013). It is shown that, under the active control, the wake symmetry can be restored or its asymmetry can be amplified.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical prediction of a turbulent curved wake and comparison with experimental data

Numerical studies of the curved wake of a NACA 0012 airfoil of chord length 0.150 m are presented. The airfoil is placed in air at 10 m/s in the straight section of a duct of 0.457 m × 0.457 m cross‐section followed by a 90° bend with a mean radius‐to‐height ratio of 1.17. The trailing edge is located at one chord length upstream of the bend entry plane. The authors' own measurements are used to define the boundary conditions and for comparison with the predicted results. The numerical models are based on the time‐averaged, three‐dimensional conservation equations of fluid flow, incorporating the k–ε, RNG k–ε, realizable k–ε and the Reynolds stress turbulence models. The results show that the models are capable of predicting the effects of curvature on the wake development. However, quantitative differences between prediction and experiment exist. The results obtained using the Reynolds stress model show better agreement with the experimental data, compared with the k–ε based models, but not consistently for all parameters. There are also better predictions by the RNG k–ε and realizable k–ε models compared with the standard k–ε model. The predicted results using the RNG k–ε are closer to experimental data than the realizable k–ε. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical simulation of a low-Reynolds-number turbulent wake behind a flat plate

A numerical simulation of a plane turbulent wake at a very low Reynolds number has been performed using finite volume methods. The wake was produced by allowing two turbulent boundary layers, simulated separately in advance, to interact downstream of the trailing edge of a thin flat plate. A number of innovative numerical techniques were required in the simulation, such as the provision of fully turbulent time-dependent inflow data from a separate simulation, advective outflow boundary conditions and the approximate representation of an internal solid surface by a method which is computationally efficient. The resulting simulation successfully reproduced many of the statistical properties of the turbulent near-wake flow at low Reynolds number.     

Structure of wake of a sharp-edged bluff body in a shallow channel flow

The flow field downstream of a bluff body in a typical open channel flow was explored by two-dimensional particle image velocimetry. Measurements are obtained in horizontal planes at the near-bed, mid-depth and near-surface locations downstream of the body up to a streamwise distance of 10D, where D is the width of the body. The dimensionless streamwise defect velocity profile of the wake flow matches well with the data of a previous investigation and does not reflect any dependency on the distance from the bed. However, the nature of development of the recirculation region is found to be different at the three vertical locations. The time-averaged streamline pattern indicates the existence of a unique nodal pattern close to the bed. The variation of the half-width is also found to be affected by the presence of the bed and the free surface. The bed friction arrests the transverse growth of the shear layer, and the free-surface helps to redistribute the turbulent kinetic energy in the streamwise and transverse directions. Swirling strength analysis is carried out to compare the behavior and statistics of the vortex population in the vertical direction. The prevailing magnitude of the swirling strength is found to be different at the three vertical locations. Bed friction assists to dissipate vorticity rapidly, and therefore reduces the probability of appearance of strong vortices close to the bed.     

Direct numerical simulation of turbulent wake behind a surface-mounted square cylinder

In this research, direct numerical simulation has been performed to study the turbulent wake behind a wall-mounted square cylinder with aspect ratio 4 and Reynolds number 12 000 (based on the free-stream velocity and obstacle side length) in a developing boundary layer. Owing to the relatively high Reynolds number and high aspect ratio of the cylinder tested, the wake is wide spread behind the cylinder and exhibits complex and energetic vortex motions. The lateral and tip vortex shedding patterns at different frequencies, coherent structures downstream of the obstacle, the production rate and distribution of turbulent kinetic energy, and the instantaneous pressure distribution in the wake region have been thoroughly investigated. In order to validate the numerical results, the first- and second-order flow statistics obtained from the simulations have been carefully compared against available wind-tunnel measurement data.     

A unified theory for turbulent wake flows described by eddy viscosity

In this paper, we consider the conservation laws for the far downstream wake equations described by eddy viscosity. A basis of conserved vectors is constructed. The well-known conserved quantities for the turbulent classical wake and the turbulent wake of a self-propelled body are obtained by integrating the corresponding conservation law across the wake and imposing the boundary conditions. For the wake of a self-propelled body the additional condition that the drag on the body is zero and is required to obtain the conserved quantity. A third conservation law, which possibly belongs to another type of wake, is discovered. The Lie point symmetry associated with the conserved vector is used to obtain the invariant solution and a typical velocity profile for this wake is provided. This wake appears to have common properties with the other two well-known wakes. We then analyse the invariant solutions to all three wake problems and prove that a simple mathematical relationship exists between them thus unifying the theory for turbulent wake flows.     

Mechanisms of wake deflection angle change behind a heaving airfoil

An immersed-boundary numerical method is applied to simulate the wake downstream of a two-dimensional heaving airfoil. A switch of vortex pattern is found to be the major reason that a deflected asymmetric wake reverses its deflection angle. Parameters of the heaving airfoil and flow that influence the onset and location of the vortex switching are discussed. While the previous literature deliberately discussed the wake deflection in the near wake region, this study shows that the deflection angle can change from the near wake to far wake regions. A cross-flow effective phase velocity is introduced to analyze the already-formed asymmetric wake behind the airfoil. A vortex dipole model and the related vortex dynamics analysis are developed to show that the change of the distance between the vortices is the primary factor that leads to the vortex pattern switching in the far wake.     

Three-dimensional evolution of vortices and early features of coherent structure in the turbulent wake behind a 2-D circular cylinder

3-D evolution of Kármán vortex filaments and vortex filaments in braid regions in the turbulent wake of a 2-D circular cylinder is investigated numerically based on inviscid vortex dynamics by analyzing the response of the initially 2-D spanwise vortex filaments to periodic spanwise disturbance of varying magnitude, wavelength and initial phase angles. Our results reveal a kind of 3-D vortex system in the wake which consists of large scale horseshoe-shaped vortices and small scale γ-shaped vortex filaments as well as vortex loops. The mechanism and the dynamic process about the generation of streamwise vortical structure and the 3-D coherent structure are reported. currently published in the Chinese Edition of Acta Mechanica Sinica, Vol.25, No.3, 1993 The project supported by National Natural Science Foundation of China and the National Basic Research Project “Nonlinear Science”     

Numerical modeling of a turbulent flow behind a heated grid

A numerical model of the dynamics of turbulence and temperature fluctuations behind a heated grid located in a wind tunnel is constructed on the basis of closed Kármán-Howarth and Corrsin equations. Results calculated by this model are in reasonable agreement with available experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 118–126, May–June, 2009.     

An efficient parallel high-order compact scheme for the 3D incompressible Navier–Stokes equations

This article provides a strategy for solving incompressible turbulent flows, which combines compact finite difference schemes and parallel computing. The numerical features of this solver are the semi-implicit time advancement, the staggered arrangement of the variables and the fourth-order compact scheme discretisation. This is the usual way for solving accurately turbulent incompressible flows. We propose a new strategy for solving the Helmholtz/Poisson equations based on a parallel 2d-pencil decomposition of the diagonalisation method. The compact scheme derivatives are computed with the parallel diagonal dominant (PDD) algorithm, which achieves good parallel performances by introducing a bounded numerical error. We provide a new analysis of its effect on the numerical accuracy and conservation features. Several numerical experiments, including two simulations of turbulent flows, demonstrate that the PDD algorithm maintains the accuracy and conservation features, while conserving a good parallel performance, up to 4096 cores.     

A numerical study of the wake behind a blunt‐trailing‐edged body. Part 2: the topology of the flow

This project focuses on the wake behind a two‐dimensional blunt‐trailing‐edged body. The body has an elliptical nose followed by a straight section that ends in a blunt base. For this geometry, the boundary layer separation points are defined and fixed. This permits a better assessment of the vital influence of the boundary layers upon the wake. It is hoped that the study of this special geometry will contribute to a better understanding of the wake flow. The main aim of this article is to study the topology of the wake by means of direct numerical simulation (DNS). The idea is to mimic some experimental work techniques. In an experiment, by means of proper synchronization, one is able to arrive at the mean‐in‐the‐phase fields. In the present work, similar data are obtained through the use of numerical tools. The topology of the wake flow in terms of the movement of the coherent structures, fields of Reynolds stresses, position of saddles and centers, is obtained for certain values of the body's aspect ratio. A comparison with the circular cylinder experimental data is also provided. An important result is the finding of some specific stress distributions around the main structures. The stress level lines are a hint of where the energy transfer of the mean flow to the turbulence is more intense. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical prediction of turbulent flow over a two-dimensional ridge

Predictions are presented of the two-dimensional turbulent flow over a triangular ridge. The time-averaged Reynolds equations are written in an orthogonal curvilinear co-ordinate system and transformed to finite difference form after discretization in physical space. Turbulence is simulated by the two-equation κ-ε model of turbulence. In the first part of the paper the basics of the numerical method are presented and in the second part comparisons are made between predictions and available laboratory data. Therefore the validity and reliability of the method as well as its flexibility in treating complex recirculating flows are assessed. Results of engineering significance are presented of the effect of the ridge slope on the length of the recirculation region and on the overspeed factor on top of the ridge.     

An improved approximation technique to obtain numerical solution of a class of two-point boundary value similarity problems in fluid mechanics

A simple and efficient approximate numerical technique is presented to obtain solutions to a wide class of two-point boundary value similarity problems in fluid mechanics. This technique is based on the common finite difference method with central differencing, a tridiagonal matrix manipulation and an iterative procedure. The technique described in this paper has been successfully applied to three different representative similarity problems of fluid mechanics. Each one of these problems is described by a coupled, non-linear system of three ordinary differential equations and has already been solved elsewhere using a different numerical method. So, the obtained numerical results, by our efficient numerical technique, permit a comparative study and show the accuracy and the effectiveness of this technique.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical modeling of the turbulent viscous-gas flow past a cylinder with a circular vortex cell in the presence of suction from the central body surface

The compressibility effect on the cylinder drag reduction due to air suction through the surface of a central body in a circular vortex cell is estimated on the basis of the solution of the steady Reynolds equations closed by the shear stress transfer model, together with the continuity, energy, and state equations.     

Quasi-Euler-Poinsot motion of a rigid body

The phase-plane method of nonlinear oscillation is used to discuss the influence of the small dissipation upon the Euler-Poinsot motion of a rigid body about a fixed point. The equations of phase coordinates are applied instead of Eulerian equations, and the global characteristics of the motion of rigid body are analysed according to the distribution and the type of the singular points. A Chaplygin's sphere on a rough plane, a rigid body in viscous medium and one with a cavity filled with viscous fluid are discussed as examples. It is shown that the motions of rigid bodies dissipated by various physical factors have a common qualitative character. The rigid body tends to make a permanent rotation about the principal axis of the largest moment of inertia. The transitive process can change from oscillatory to aperiodic with the decrease in dissipation.     

Stress concentration of an ellipsoidal inclusion of revolution in a semi-infinite body under biaxial tension

Summary This paper deals with the stress concentration problem of an ellipsoidal inclusion of revolution in a semi-infinite body under biaxial tension. The problem is formulated as a system of singular integral equations with Cauchy-type or logarithmic-type singularities, where unknowns are densities of body forces distributed in the r- and z-directions in semi-infinite bodies having the same elastic constants as the ones of the matrix and inclusion. In order to satisfy the boundary conditions along the ellipsoidal boundary, four fundamental density functions proposed in [24, 25] are used. The body-force densities are approximated by a linear combination of fundamental density functions and polynomials. The present method is found to yield rapidly converging numerical results for stress distribution along the boundaries even when the inclusion is very close to the free boundary. The effect of the free surface on the stress concentration factor is discussed with varying the distance from the surface, the shape ratio and the elastic modulus ratio. The present results are compared with the ones of an ellipsoidal cavity in a semi-infinite body.accepted for publication 11 November 2003     

Verification of a mixed high‐order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high‐order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier–Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity–velocity formulation for a two‐dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high‐order compact and non‐compact finite‐differences from fourth‐order to sixth‐order of accuracy. The other scheme used spectral methods instead of finite‐difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed‐order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright © 2009 John Wiley & Sons, Ltd.