UFR 4-03 Test Case: Difference between revisions

From KBwiki
Jump to navigation Jump to search
Line 61: Line 61:
In cylindrical coordinates, as mentioned before, the equations for the velocity components are singular at r=0. If staggered "fluxes" are used (see Figure 2), one can set [[Image:U4-03d32_files_image072.gif]]at the centreline, and the [[Image:U4-03d32_files_image074.gif]]  equation can be discretized at the points ([[Image:U4-03d32_files_image076.gif]]) without any approximation. The viscous terms in the radial direction at the first point close to the wall, for the [[Image:U4-03d32_files_image078.gif]] and [[Image:U4-03d32_files_image079.gif]] equations have been approximated as done by Rai and Moin (1991) for the plane channel. The discretization of the [[Image:U4-03d32_files_image080.gif]] and [[Image:U4-03d32_files_image082.gif]] equations at the first grid position near the axis does not require any approximation for the radial derivatives. Periodicity has been imposed in the stream-wise direction in a domain of different length [[Image:U4-03d32_files_image084.gif]] depending on the rotation number N.
In cylindrical coordinates, as mentioned before, the equations for the velocity components are singular at r=0. If staggered "fluxes" are used (see Figure 2), one can set [[Image:U4-03d32_files_image072.gif]]at the centreline, and the [[Image:U4-03d32_files_image074.gif]]  equation can be discretized at the points ([[Image:U4-03d32_files_image076.gif]]) without any approximation. The viscous terms in the radial direction at the first point close to the wall, for the [[Image:U4-03d32_files_image078.gif]] and [[Image:U4-03d32_files_image079.gif]] equations have been approximated as done by Rai and Moin (1991) for the plane channel. The discretization of the [[Image:U4-03d32_files_image080.gif]] and [[Image:U4-03d32_files_image082.gif]] equations at the first grid position near the axis does not require any approximation for the radial derivatives. Periodicity has been imposed in the stream-wise direction in a domain of different length [[Image:U4-03d32_files_image084.gif]] depending on the rotation number N.


Finite differences have the advantages to use arbitrary coordinates transformation to cluster grid points in regions of high gradients. An analytical transformation by the hyperbolic tangent allowed to have the first grid point at a distance [[Image:U4-03d32_files_image086.gif]]from the wall. This clustering is necessary to fully resolve the thin vortical layers near the wall responsible for the wall friction and for the turbulence production. Although the Reynolds numbers is low the grid spacing in the central region of the pipe is not small enough to consider the flow fully resolved. However the grid spacing was close to those used in the full channel by KMM. Moreover it is clear that the resolution for the circular pipe is better than in the plane channel, since the polar coordinates reduce the size [[Image:U4-03d32_files_image088.gif]]of the physical grid moving towards the axis. In the three-dimensional case, in the limit of [[Image:U4-03d32_files_image090.gif]]energy is conserved by the discretized equations. The check of the accurate treatment of the singularity was done by letting a dipole to move across the axis where large variations of radial velocity occur. The simulation showed that the dipolar structure maintains its shape and that in the inviscid case the energy was conserved.
Finite differences have the advantages to use arbitrary coordinates transformation to cluster grid points in regions of high gradients. An analytical transformation by the hyperbolic tangent allowed to have the first grid point at a distance [[Image:U4-03d32_files_image086.gif]]from the wall. This clustering is necessary to fully resolve the thin vortical layers near the wall responsible for the wall friction and for the turbulence production. Although the Reynolds numbers is low the grid spacing in the central region of the pipe is not small enough to consider the flow fully resolved. However the grid spacing was close to those used in the full channel by KMM. Moreover it is clear that the resolution for the circular pipe is better than in the plane channel, since the polar coordinates reduce the size [[Image:U4-03d32_files_image088.gif]]of the physical grid moving towards the axis. In the three-dimensional case, in the limit of [[Image:U4-03d32_files_image090.gif]] energy is conserved by the discretized equations. The check of the accurate treatment of the singularity was done by letting a dipole to move across the axis where large variations of radial velocity occur. The simulation showed that the dipolar structure maintains its shape and that in the inviscid case the energy was conserved.


Feiz, Oul-Rouis & Lauriat (2003) performed the Large Eddy Simulation with two subgrid models, a dynamic model and the usual Smagorinsky model. In the Smagorinsky model, the subgrid scale eddy viscosity is related to the deformation of the resolved velocity field as:
Feiz, Oul-Rouis & Lauriat (2003) performed the Large Eddy Simulation with two subgrid models, a dynamic model and the usual Smagorinsky model. In the Smagorinsky model, the subgrid scale eddy viscosity is related to the deformation of the resolved velocity field as:


[[Image:U4-03d32_files_image092.gif]] [[Image:U4-03d32_files_image094.gif]]. This subgrid model, largely used in LES of isotropic turbulence produced good results. When it was applied to inhomogeneous and in particular to wall bounded flows the constant was modified. The dynamic model provides a methodology for determining an appropriate local value of the Smagorinsky coefficient. The model was proposed by Germano ''et al.'' (1991) and later improved by Lilly (1992). The turbulent viscosity is expressed using an eddy viscosity assumption as [[Image:U4-03d32_files_image096.gif]]with the value of the constant [[Image:U4-03d32_files_image098.gif]]computed during the simulation from the flow variables.
[[Image:U4-03d32_files_image092.gif]] where [[Image:U4-03d32_files_image094.gif]]. This subgrid model, largely used in LES of isotropic turbulence produced good results. When it was applied to inhomogeneous and in particular to wall bounded flows the constant was modified. The dynamic model provides a methodology for determining an appropriate local value of the Smagorinsky coefficient. The model was proposed by Germano ''et al.'' (1991) and later improved by Lilly (1992). The turbulent viscosity is expressed using an eddy viscosity assumption as [[Image:U4-03d32_files_image096.gif]]with the value of the constant [[Image:U4-03d32_files_image098.gif]]computed during the simulation from the flow variables.


[[Image:U4-03d32_files_image100.gif]] [[Image:U4-03d32_files_image102.gif]] [[Image:U4-03d32_files_image104.gif]]. The constant could be positive or negative. The positive values are linked to energy flowing from large to small scales, and the negative to energy going from small to large scales (backward energy transfer). A grid (65, 39, 65), non uniform in the radial direction and uniform in the other two directions was used. The pipe length was [[Image:U4-03d32_files_image106.gif]] [[Image:U4-03d32_files_image108.gif]]and two different rotation rates [[Image:U4-03d32_files_image110.gif]] [[Image:U4-03d32_files_image112.gif]]  
[[Image:U4-03d32_files_image100.gif]] [[Image:U4-03d32_files_image102.gif]] [[Image:U4-03d32_files_image104.gif]]. The constant could be positive or negative. The positive values are linked to energy flowing from large to small scales, and the negative to energy going from small to large scales (backward energy transfer). A grid (65, 39, 65), non uniform in the radial direction and uniform in the other two directions was used. The pipe length was [[Image:U4-03d32_files_image106.gif]] [[Image:U4-03d32_files_image108.gif]]and two different rotation rates [[Image:U4-03d32_files_image110.gif]] [[Image:U4-03d32_files_image112.gif]]  

Revision as of 13:55, 8 March 2009


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References




Pipe flow - rotating

Underlying Flow Regime 4-03               © copyright ERCOFTAC 2004


Test Case

Brief description of the study test case

In the three-dimensional case, expressing the Navier-Stokes equations in primitive variables reduces the number of operations necessary to integrate the equations numerically. In cylindrical coordinates, when the variables U4-03d32 files image020.gif are introduced, the continuity equations becomes:

U4-03d32 files image022.gif

It is clear that the quantities U4-03d32 files image024.gifdo not all have the same physical dimensions; this could be a matter of confusion. We wish to point out that these variables have been introduced for numerical reasons only. As discussed later, the introduction of U4-03d32 files image026.gif simplifies the treatment of the singularity at r=0 and U4-03d32 files image028.gif permits better accuracy near the axis. The Navier-Stokes equations in terms of the new variables, in a reference frame rotating with the pipe wall, are:

U4-03d32 files image030.gif

 

U4-03d32 files image032.gif

 

U4-03d32 files image034.gif

The equations nondimensionalized by the centreline stream-wise velocity U4-03d32 files image036.gif of the laminar Poiseuille flow and by the pipe radius R give the dimensionless numbers U4-03d32 files image038.gif and U4-03d32 files image040.gif. In this paper rather than the Rossby number U4-03d32 files image042.gif usually used in geophysical problems, the rotation number U4-03d32 files image044.gif has been adopted. This choice has been dictated by the fact that N was previously used in studies of turbulent flows in rotating pipes. The mean pressure gradient in the streamwise direction maintains the bulk velocity U4-03d32 files image046.gifconstant; it is evaluated by integrating the discretized U4-03d32 files image048.gif equation in the whole domain. U4-03d32 files image050.gif is related to the mean wall friction velocity U4-03d32 files image052.gif by U4-03d32 files image054.gif .

U4-03d32 files image056.jpg

Figure 1 Flow Configuration.

The flow configuration is shown in Figure 1.

Test Case Experiments

Durst Jovanovic & Sender (1995) carried out LDA measurements in the near wall region of a turbulent pipe flow. This test experiment has been chosen to validate the CFD at N=0 (no-rotation). The experimental facility consist of a closed--loop pipe flow installation driven by a screw conveyor pump operating in suction mode. This screw conveyor pump was chosen because of its very low flow-rate pulsations. Electronic control of the pump was available which allowed flow--rate setting with a tolerance of 1%. Reynolds number ranged from 7442 to 30.000. In order to permit laser-Doppler measurements, the test section was made from a Duran-50 glass pipe. To avoid interference of the pipe with the light beams, the test fluid was a two mixed DIESEL oil in order to be refractive-index-matched to the Duran glass. Heating and cooling system controlled the flow temperature. The laser--Doppler is composed of a 30 mW helium neon laser and double Brag cell transmission optics.

For low rotation rates the experiment by Kikuyama et al. (1983) was considered. The velocity profile at the inlet section was made uniform and with low turbulence levels (.3%) with a honeycomb and 6 screens and a nozzle with aspect ratio 12.3:1. A step ring of 1mm height was installed at the exit section of the nozzle to ensure to flow be turbulent in the boundary layer. Hot wires probes were used to measure velocities at x/D=2.7,5.7,9.7,15.5,22.5 (x is the distance from the inlet section) for U4-03d32 files image058.gif. The probes were inserted through a small hole drilled in a narrow stationary ring 10mm in length placed between two rotating sections.

For higher rotation rates CFD results have been compared with the experiments carried out by Reich & Beer (1989). In the experiment by Reich & Beer (1989), in order to obtain fully developed flow conditions, the air was passed through an air filter into the inlet section, which is a non--rotating pipe with a diameter of D=50mm and a dimensionless length of L/D=50. The length of the ensuing rotating test section could be varied in the range 0<L/D<120. At the end of the rotating pipe a narrow non rotating ring of 3mm length was inserted for measurements.

A cylindrical three-hole aerodynamic probe was inserted into the rotating ring and transversed in the radial direction, in order to obtain axial and tangential velocity profiles. Behind the non--rotating ring a second rotating pipe section (L/D=10), a mixing chamber for measuring of temperature, a downstream non-rotating pipe, a Venturi nozzle for measuring the flow rate and a centrifugal blower were arranged. The rotating tube wall temperature was measured with chromel-alumel thermocouples. To prevent vibration due to rotation, the pipe was supported by a series of ball bearings located at intervals of L/D=20. The flow-rate Reynolds number was varied in the range U4-03d32 files image060.gif and the rotational Reynolds number in the range U4-03d32 files image062.gif

CFD Methods

Orlandi & Fatica (1997) have performed DNS of the turbulent rotating pipe flow. The numerical scheme described in detail in the paper by Verzicco and Orlandi (1996) is here recalled briefly. The system of equations was solved by a fractional step method (Rai and Moin 1991). In a first step a non-solenoidal velocity field U4-03d32 files image064.gif. A scalar quantity U4-03d32 files image066.gif is introduced to project the non-solenoidal field onto a solenoidal one. The large band matrix associated to the elliptic equation for U4-03d32 files image066.gif is reduced to a tridiagonal matrix by periodic FFTs in the azimuthal and in the axial directions. The updated pressure is computed from the scalar U4-03d32 files image066.gif. A third order Runge-Kutta scheme, described by Rai and Moin (1991), was used to advance in time, through three sub-steps. The viscous terms are treated implicitly by the Crank-Nicolson scheme.

U4-03d32 files image070.gif
Figure 2 Typical grid cell

In cylindrical coordinates, as mentioned before, the equations for the velocity components are singular at r=0. If staggered "fluxes" are used (see Figure 2), one can set U4-03d32 files image072.gifat the centreline, and the U4-03d32 files image074.gif equation can be discretized at the points (U4-03d32 files image076.gif) without any approximation. The viscous terms in the radial direction at the first point close to the wall, for the U4-03d32 files image078.gif and U4-03d32 files image079.gif equations have been approximated as done by Rai and Moin (1991) for the plane channel. The discretization of the U4-03d32 files image080.gif and U4-03d32 files image082.gif equations at the first grid position near the axis does not require any approximation for the radial derivatives. Periodicity has been imposed in the stream-wise direction in a domain of different length U4-03d32 files image084.gif depending on the rotation number N.

Finite differences have the advantages to use arbitrary coordinates transformation to cluster grid points in regions of high gradients. An analytical transformation by the hyperbolic tangent allowed to have the first grid point at a distance U4-03d32 files image086.giffrom the wall. This clustering is necessary to fully resolve the thin vortical layers near the wall responsible for the wall friction and for the turbulence production. Although the Reynolds numbers is low the grid spacing in the central region of the pipe is not small enough to consider the flow fully resolved. However the grid spacing was close to those used in the full channel by KMM. Moreover it is clear that the resolution for the circular pipe is better than in the plane channel, since the polar coordinates reduce the size U4-03d32 files image088.gifof the physical grid moving towards the axis. In the three-dimensional case, in the limit of U4-03d32 files image090.gif energy is conserved by the discretized equations. The check of the accurate treatment of the singularity was done by letting a dipole to move across the axis where large variations of radial velocity occur. The simulation showed that the dipolar structure maintains its shape and that in the inviscid case the energy was conserved.

Feiz, Oul-Rouis & Lauriat (2003) performed the Large Eddy Simulation with two subgrid models, a dynamic model and the usual Smagorinsky model. In the Smagorinsky model, the subgrid scale eddy viscosity is related to the deformation of the resolved velocity field as:

U4-03d32 files image092.gif where U4-03d32 files image094.gif. This subgrid model, largely used in LES of isotropic turbulence produced good results. When it was applied to inhomogeneous and in particular to wall bounded flows the constant was modified. The dynamic model provides a methodology for determining an appropriate local value of the Smagorinsky coefficient. The model was proposed by Germano et al. (1991) and later improved by Lilly (1992). The turbulent viscosity is expressed using an eddy viscosity assumption as U4-03d32 files image096.gifwith the value of the constant U4-03d32 files image098.gifcomputed during the simulation from the flow variables.

U4-03d32 files image100.gif U4-03d32 files image102.gif U4-03d32 files image104.gif. The constant could be positive or negative. The positive values are linked to energy flowing from large to small scales, and the negative to energy going from small to large scales (backward energy transfer). A grid (65, 39, 65), non uniform in the radial direction and uniform in the other two directions was used. The pipe length was U4-03d32 files image106.gif U4-03d32 files image108.gifand two different rotation rates U4-03d32 files image110.gif U4-03d32 files image112.gif

The Reynolds averaged form of the Navier-Stokes equations, under fully developed conditions (U4-03d32 files image114.gifU4-03d32 files image116.gif

U4-03d32 files image118.gif

A non-zero mean swirl velocity is generated by the shear component U4-03d32 files image120.gif of the Reynolds stress tensor. In order to guarantee that this component is non-zero and that the flow is not unidirectional, as indicated by experiments, non-zero U4-03d32 files image122.gif Reynolds shear stress must be generated by the axial mean velocity U4-03d32 files image124.gif et al. (2000) discussed the capability of several closure models to predict a swirl velocity and an axial velocity depending on U4-03d32 files image126.gif. Here, the main conclusions are summarised, details are in the paper.

The standard and nonlinear U4-03d32 files image128.gifmodels, like all traditional two equation models with conventional near-wall treatments, are incapable of predicting the main physical features of turbulent flow in an axially rotating pipe.

Since the two-dimensional explicit algebraic stress model does not generate a non-zero Reynolds shear stress U4-03d32 files image129.gif,there will be no mean swirl velocity. However, unlike the standard and nonlinear U4-03d32 files image131.gifmodels, it does predict a rotationally dependent axial mean velocity, which is the more significant effect.

The three-dimensional form of the explicit algebraic stress models is much more complicated. A mean swirl velocity arises from a cubic nonlinearity in conventional explicit algebraic stress models.

A non-zero mean swirl velocity will be predicted in axially rotating pipe flow by virtually all existing second order closures by the presence of Coriolis terms.

© copyright ERCOFTAC 2004



Contributors: Stefano Leonardi - Universita di Roma 'La Sapienza'


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References