Advances in Civil Engineering

Volume 2012 (2012), Article ID 523839, 8 pages

http://dx.doi.org/10.1155/2012/523839

## Unsteady RANS and LES Simulation of an Ideal Rankine Vortex Decay

^{1}Civil Engineering Department, Khulna University of Engineering and Technology, Khulna 9203, Bangladesh^{2}Department of Urban Management, Kyoto University, Kyoto 615-8540, Japan^{3}Faculty and Graduate School of Engineering, Hokkaido University, Sapporo 060-8623, Japan

Received 27 July 2011; Revised 7 December 2011; Accepted 3 February 2012

Academic Editor: Ghassan Chehab

Copyright © 2012 Md. Shahjahan Ali et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The 3D numerical simulation was carried out for an idealized Rankine vortex using nonlinear *k*-*ε* model (one kind of RANS model) and large eddy simulation (LES) techniques. In this 3D simulation, the vortex flow field was given to rotate with the vertical axis in a free surface rectangular domain. In order to investigate the predictability of standard (linear) and non-linear *k*-*ε* models, the decay of a trailing vortex was simulated and compared with previous DNS data. The governing equations for mean velocities and turbulent flows were discretized with the finite volume method based on a staggered grid system. It was observed that in the growth phase as well as in stabilized phase of turbulence, the decay rate of tangential velocity by RANS model was well comparable with LES simulation as well as previous DNS data. However, in the decay phase of turbulence, RANS model showed slightly faster decay of tangential velocity due to its slower decay of turbulence compared to LES or DNS. The patterns as well as magnitudes of secondary currents predicted by RANS and LES models were well comparable to each other.

#### 1. Introduction

The basic types of plane vortices can be classified into two categories: one with slower velocity at center and maximum at sides and the other with maximum at center and minimum velocity at edges. The rotary fluid motion of the first one is called the solid body rotation, since it is similar to the fluid motion filled in a rotating hollow box. On the other hand, if a long circular rod rotates in a fluid with constant velocity around its axis, the fluid velocity is found highest and equal to the velocity of rod at the rod’s surface (due to adhesion), and with increasing distance from the rod, the velocity is diminished in inverse proportional to the distance. Such a fluid motion is called a potential vortex. Fluid motion composed of a potential vortex and solid body rotation is called Rankine vortex (Figure 1) after the fluid dynamicist Rankine. The radial distance from the center to the maximum tangential velocity is called the radius of vortex core. For a steady circular motion without a velocity component normal to the plane of rotation, the Rankine vortex is the only possible vortex whose velocity is zero at the center as well as far away from it. In addition to these basic vortices, there are other time-dependent rotary motions that have azimuthal velocity component as well as radial and axial components.

The existence of such vortices is not only limited to natural phenomena but also in many engineering applications. The Rankine vortex has been used extensively in various studies: for instance, to model the natural phenomena such as hurricanes and tornados, to study the vortex wake hazard caused by flow obstructions, to predict the decay of wing-tip vortices, and so forth. Although some experimental and theoretical studies ([1–3], etc.) as well as DNS simulations [4] were successfully performed to study the structure and decay of such vortex, all of them were performed considering 2D vortex field.

As long as the flow is laminar, the only problem is to solve a system of discretized equations. The accuracy of results depends on both the grid and discretization scheme. The moment that the flow becomes turbulent, a chaotic, random motion is observable. In this case, the instantaneous flow values at a fixed location in the flow field show a fluctuating behavior. Therefore, in order to study such flows numerically, different approaches have been proposed. They have different accuracies, and sometimes their predictability depends on the type of flows. The major approaches to study turbulent flows numerically are Reynolds-averaged Navier-Stokes (RANS), large eddy simulation (LES), and direct numerical simulation (DNS). In this study, the 3D numerical simulations were carried out using non-linear model (one kind of RANS model) and large eddy simulation (LES) techniques. The idea of LES comes from the fact that the small scales show an isotropic behavior regardless of the type of flow and boundary conditions. Hence, if the large scales are resolved, it is accurate enough to model the small scales which are called the subgrid scale or SGS. In order to get reasonable accuracy, the grid should be fine to some extent. Important difference between RANS and LES governing equations is the method of filtering. While in RANS the filtering is performed in time, in LES the filtering is performed in space. Another difference compared to RANS is that, as turbulence is a time-dependent three-dimensional phenomenon, the governing equations in LES also need to solve in unsteady three-dimensional format. Another fact regarding LES is that the grid is supposed to be fine enough to capture some part of scales in the inertial subrange. These differences make the LES much more expensive than RANS models.

The work presented in this paper can be classified into twofold. In the first part, unsteady numerical simulation with 3D velocity field was carried out for a turbulent trailing vortex using non-linear model. To assess the predictability of non-linear model, general features of vortex decay are compared with previous DNS data. The simulation was also performed using standard model and comparison is shown with non-linear model. In the second part, the 3D numerical simulations were carried out for an idealized Rankine vortex using non-linear model and large eddy simulation (LES). 3D staggered grid was used for both RANS and LES simulations. Comparison is shown for the general flow features such as temporal change of vortex decay, radial distribution of tangential velocity, and water surface profile.

#### 2. Numerical Model and Simulation Details

##### 2.1. Nonlinear Model

Following 3D flow equations [5–7] were used in unsteady RANS calculation for the chosen flow field.

Continuity equation:

Momentum equation:

-equation:

-equation:

Turbulence model (nonlinear):

Here,

In (5),
The values of the model constants are = 0.09, = 0.005, = 0.0068, = 0.008, = 0.004, = −0.003, = 0.00005, = 0.00005, = 0.00025, = 0.40, = 0.0, = −0.13, = 0.01, = 0.003, = 1.0, = 1.3, = 1.44, and * *= 1.92. In the standard model, = 0.09 and = 0.0.

##### 2.2. Large Eddy Simulation

The filtered equations were developed from the Navier-Stokes equations. Consider that is a velocity in direction that contains resolved scale part and subgrid-scale part ; then ; similarly for pressure. Substituting these two decompositions into Navier-Stokes equations, the filtered equations of LES can be obtained as follows:

Continuity equation.

Momentum equation: Here, is the subgrid scale Reynolds stress defined by The most commonly used subgrid-scale turbulence model is Smagorisky model proposed by Smagorinsky [8]. In this model, the residual stress takes the Boussinesq eddy viscosity form: Here, the resolved strain rate is .

The eddy viscosity is determined by the Smagorisky model: Δ is the filter width related to grid size calculated as the geometric average of the grid spacing in three directions, . is the Smagorisky model constant, whose value varies from 0.065 to 0.20 [9]; 0.15 was used in the present simulation.

##### 2.3. Computational Schemes and Flow Domain

The basic equations used in LES were three-dimensional, time-dependent, filtered Navier stokes equations. In nonlinear model, the turbulent stress terms were solved using two additional equations for turbulent kinetic energy () and its dissipation rate (). The governing equations for mean velocities and turbulent flows were discretized with the finite volume method based on a staggered grid system. For the momentum equation, convective and diffusive fluxes were approximated with quick and central difference schemes, respectively. The hybrid central upwind scheme was used for the and equations. Time advancement was achieved by Adam-Bashforth scheme of second-order accuracy, in each equation. The basic equations were discretized as fully explicit forms and solved successively with the time increment in step by step. The pressure field was solved using iterative procedure at each time step. Three cases have been studied. The initial conditions and simulation details are explained as follows, and the simulation parameters are shown in Table 1.

###### 2.3.1. Case 1: Simulation of Trailing Vortex by Unsteady RANS

For the -vortex, the direct numerical simulation results were reported by Qin [4]. The simulation for case 1 was performed with the same initial condition of Qin. Since the axial velocity of this vortex is not a function of axial distance but of radial distance, 2D numerical grid was used to simulate 3D flow field. The details of the initial flow field are given as follows.

Tangential velocity:

Radial velocity:

Axial velocity: Here, is the scaling velocity, related to the initial mean tangential flow, defined as Here, = 0.009 and = 1.1209 were used in Qin’s [4] DNS simulation. As reported by Lessen et al. [10], the stability of -vortex can be related to the value of swirl number . In DNS simulation, the initial swirl number was chosen to be unity. The time is nondimensionalized by at . Here, is the radial distance where the tangential velocity contains the peak value in the initial flow field (at ).

Variable grid spacing was used with dense grid at centre and coarser towards the boundary. The variable mesh size with radial distance for case 1 is shown in Figure 2. The computational domain consists of 100 grids in each ( and ) direction. The size of domain was taken sufficiently large () to overcome the interference of boundary in the vortex decay process. Cartesian grid was used for computation, and the results are presented in cylindrical coordinate using the geometric conversion. The value of was taken as 0.01 sec. The truncation error was allowed as 0.000001 and the maximum number of iterations per time step was set to be 200.

###### 2.3.2. Case 2 and Case 3: Simulation of an Ideal Rankine Vortex by Unsteady RANS and LES

These two cases were performed for the same vortex flow field using different simulation methods: case 2 by RANS and case 3 by LES. The initial condition of the vortex is given as follows.

Tangential velocity:

Radial velocity:

Axial velocity: Here, is the scaling velocity, related to the initial mean tangential flow as defined in (16). Like case 1, in these two cases, the initial swirl number was also chosen to be unity.

Simulations were performed using 3D Cartesian grid. The size of flow domain for both the cases was (in two lateral, and directions) (in axial/depth-wise, direction). For LES grids were used in , , and directions, respectively. For RANS simulation, the grid numbers were reduced to . 3D staggered variable grid system was used for both RANS and LES. The other parameters for the simulation are shown in Table 1.

#### 3. Results and Discussions

At first, the decay of a trailing vortex simulated by standard and non-linear model is compared with previous DNS result. Then the decay of an ideal Rankine vortex predicted by nonlinear model and large eddy simulation (LES) is compared. The comparative discussion of simulated results is presented here in after.

##### 3.1. Comparison of RANS Prediction with Qin’s DNS Data

###### 3.1.1. Comparison of Temporal Decay

To investigate the predictability of standard (linear) and non-linear model, the decay of trailing vortex as defined in case 1 was simulated with the same initial conditions of Qin [4] using both models. In this section, the simulated results are compared with Qin’s DNS data. Comparison of temporal decay of maximum tangential velocity is shown in Figure 3. It is clearly seen that the standard model predicts much faster decay than that of non-linear model. Although the non-linear model shows better comparison with DNS calculation, the decay rate is still faster than that of DNS prediction, especially at about > 4.0.

The radial distribution of circumferential velocity, circulation, and axial velocity is compared with the DNS calculation at time = 3.72, as shown in Figures 4, 5, and 6, respectively. At this age the vortex has already overcome the effect of initial conditions and the turbulent flow field became saturated and stabilized gaining a peak value of the turbulent kinetic energy after passing its growth periods. In these figures, the non-linear model shows much better agreement with DNS data compared to standard model. Even for non-linear model, the magnitude and position of maximum tangential velocity show slightly faster decay than DNS data. Such discrepancy is also seen in the comparison of circulation profile. The cause of deviation is explained by comparing the temporal decay of vortex with the growth/decay of turbulence in the next section.

###### 3.1.2. Temporal Growth/Decay of Turbulence

Production of turbulence in the vortex field is due to the tangential and axial velocity gradients. In the vortex lifetime, five different growth/decay rates of turbulent kinetic energy as reported in DNS calculation [4] were also reflected in the simulation of non-linear model. Figure 7 presents the growth/decay of turbulence kinetic energy with approximated five different time periods. Initially (zone (i)), the vortex changes very slowly, as it requires adjustment of any unphysical nature of initial conditions. In time zone (ii), an exponential growth of turbulent kinetic energy is seen. The growth slows down in the next time period of zone (iii). It is found that the growth of axial velocity is significantly higher in time zone (ii) and (iii). The turbulent kinetic energy reaches its peak value and remains almost constant throughout the short period of zone (iv). It reveals that in this stabilization period the flow field becomes saturated and cannot support additional turbulence. Finally, the decay of turbulence is started as predicted by most of the previous researches (such as [2, 11], etc.). The logarithmic plot depicts that the decay rate in zone (v) is much slower than the growth rate in zone (ii). The decay of velocity field slows down as the turbulence decay period starts.

Comparing Figure 3 with Figure 7, we can conclude that in the growth phase as well as in the stabilized phase of turbulence, the decay rate of tangential velocity by nonlinear model is well comparable with DNS result. However, in the decay phase of turbulence, RANS model shows faster decay of tangential velocity due to its slower decay of turbulence compared to DNS.

Figure 8 depicts the comparison of the radial distribution of turbulent kinetic energy using three models. It is seen from the figure that the model overpredicts the kinetic energy at the peak compared to DNS. This is the cause for the faster decay of tangential velocity in the model (as observed in Figure 4).

##### 3.2. Decay of an Ideal Rankine Vortex by LES and RANS Simulations

In this 3D simulation, the vortex flow field was given to rotate with the vertical axis in a free surface rectangular domain. The temporal change of velocity at a point is shown in Figure 9 for both LES and RANS simulations. In addition to temporal decay of velocity, the typical velocity fluctuations are observed in LES results; on the other hand, RANS simulation shows a mean decay profile.

###### 3.2.1. Comparison of Temporal Decay

Figure 10 shows the radial distribution of tangential velocity predicted by LES and RANS models at = 3.06 and = 4.48. Although the velocity profiles at = 3.06 show good comparison, RANS model shows faster decay than LES at = 4.48.

The variation of water surface for different times predicted by RANS and LES is shown in Figure 11. The water surface is found to be depressed near the central area with minimum water depth at the centre of vortex. Since the strength of vortex decreases due to the decay of tangential velocity, the amount of depression in water surface is also decreased with time. Similar to Figure 10, although the depression of water surface at = 3.06 shows good comparison, RANS model predicted less depression in water surface than that of LES at = 4.48.

Figure 12 shows the calculated temporal decay of maximum tangential velocity for both LES and non-linear model. It is found that although the prediction of RANS model shows good agreement with LES in initial times, some deviation is observed at about > 4.0. Such discrepancy in the decay phase of turbulence is also observed in the comparison of non-linear model predictions with previous DNS results of trailing vortex decay (case 1). Therefore, similar to previous case, it can be concluded that in the growth phase as well as in stabilized phase of turbulence, the decay rate of tangential velocity by RANS model is well comparable with LES. However, in the decay phase of turbulence, RANS model shows slightly faster decay of tangential velocity due to its slower decay of turbulence compared to LES.

###### 3.2.2. Secondary Currents

The secondary currents calculated by both the models are shown in Figure 13. The flow vectors of secondary current show three circulation cells in each side of the centerline. Secondary current inside the vortex core is found relatively weaker than the outside of core. It is explained that (in Figure 1) the Rankine vortex is a combination of nearly solid body rotation and potential vortex. For both type of vortices, the individual patterns of secondary currents are reported by Laugt [12]. The pattern of secondary currents predicted by LES as well as RANS is found well agreed to Laugt. In the vortex core, the secondary current shows the pattern similar to that of solid body rotation; outside the core, the pattern is similar to that of potential vortex. The patterns as well as magnitudes of secondary currents predicted by two models are well comparable to each other.

#### 4. Conclusions

3D numerical simulations were carried out for an ideal Rankine vortex using both RANS (non-linear ) and LES models. In the growth phase as well as in stabilized phase of turbulence, the decay rate of tangential velocity by RANS model was found to be well comparable with LES. However, in the decay phase of turbulence, RANS model showed faster decay of tangential velocity compared to LES due to its slower decay of turbulence. The patterns as well as magnitudes of secondary currents predicted by two models are in good agreement with each other. The water surface was found to be depressed near the central area of the vortex with a minimum water depth observed at the centre of vortex.

In order to investigate the predictability of standard (linear) and non-linear models, the decay of a trailing vortex was simulated and compared with previous DNS data. The comparison of non-linear model prediction with previous DNS result showed similar phenomenon as observed in the comparison with LES as described previously. Due to the inherent deficiency of standard model, it predicted much faster decay than that of non-linear model.

#### References

- W. R. C. Phillips, “The turbulent trailing vortex during roll-up,”
*Journal of Fluid Mechanics*, vol. 105, pp. 451–467, 1981. View at Google Scholar · View at Scopus - M. S. Uberoi, “Mechanisms of decay of laminar and turbulent vortices,”
*The Journal of Fluid Mechanics*, vol. 90, no. 2, pp. 241–255, 1979. View at Google Scholar · View at Scopus - P. G. Saffman, “The structure and decay of turbulent trailing vortices,”
*Archives of Mechanics*, vol. 26, no. 3, pp. 423–439, 1974. View at Google Scholar · View at Scopus - J. H. Qin,
*Numerical simulation of a turbulent axial vortex*, Ph.D. thesis, Purdue University, West Lafayette, Ind, USA, 1998. - I. Kimura, W. S. J. Uijttewaal, T. Hosoda, and M. S. Ali, “URANS computations of shallow grid turbulence,”
*Journal of Hydraulic Engineering*, vol. 135, no. 2, pp. 118–131, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. S. Ali, T. Hosoda, and I. Kimura, “Unsteady simulation of turbulent axial vortex using non-linear $k-\epsilon $ model,”
*Journal of Applied Mechanics*, vol. 11, pp. 821–832, 2008. View at Google Scholar - I. Kimura and T. Hosoda, “A non-linear k-
*ε*model with realizability for prediction of flows around bluff bodies,”*International Journal for Numerical Methods in Fluids*, vol. 42, no. 8, pp. 813–837, 2003. View at Publisher · View at Google Scholar · View at Scopus - J. Smagorisky, “General circulation experiments with primitive equations-I: the basic experiment,”
*Monthly Weather Review*, vol. 91, pp. 99–165, 1963. View at Google Scholar - J. H. Ferziger and M. Peric,
*Computational Methods for Fluid Dynamics*, Springer publications, New York, NY, USA, 2nd edition, 1999. - M. Lessen, P. J. Singh, and F. Paillet, “The stability of a trailing line vortex. Part I: inviscid theory,”
*Journal of Fluid Mechanics*, vol. 63, no. 4, pp. 753–763, 1974. View at Publisher · View at Google Scholar - G. K. Bachelor, “Axial flow in trailing line vortices,”
*The Journal of Fluid Mechanics*, vol. 20, no. 4, pp. 645–658, 1964. View at Google Scholar · View at Scopus - H. J. Laugt,
*Vortex Flow in Nature and Technology*, John Wiley & Sons, New York, NY, USA, 1983.