UFR 4-03 Evaluation

From QNET
Jump to navigation Jump to search

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References




Pipe flow - rotating

Underlying Flow Regime 4-03               © copyright ERCOFTAC 2004


Evaluation

Comparison of CFD calculations with Experiments

N

H(pres)

H(MK)

U4-03d32 files image133.gif (pres)

U4-03d32 files image135.gif (RB)

U4-03d32 files image137.gif (pres)

U4-03d32 files image139.gif (RB)

U4-03d32 files image140.gif (KM)

0

1.637

1.625

1.306

1.270

1.

1.

1.

0.5

1.747

1.687

1.440

1.350

0.8386

0.9230

0.8523

1

1.771

1.812

1.503

1.457

0.8290

0.7692

0.7035

2

1.913

1.969

1.660

1.689

0.8247

0.6670

0.5898

Table 1 values of the global quantities in the present simulation indicated by (pres Re=5000),in the experiment by Reich & Beer (1989) (RB Re=5000), in that by Murakami & Kikuyama (1987) (MK Re=10000) and in that by Nikuyama et al (1983) (KM Re=7500)

In Tab.1 the results of Reich & Beer (1989) at Re=5000, those by Nishibori et al. (1987) at Re=10000 and those by Kikuyama et al. (1983) at Re=7500 and at U4-03d32 files image142.gif are compared with the DNS by Orlandi & Fatica (1997). The shape factor H is defined as U4-03d32 files image144.gif

In Tab.1 H and U4-03d32 files image146.gifobtained by the DNS agree satisfactorily well with the experiments. On the contrary the wall friction U4-03d32 files image148.gif reduction in the simulations is lower than that in the experiments, showing a negligible dependence on N for N>0.5.

The lower Re number of the simulations could be a cause of this difference since the experiments by Kikuyama et al. (1983) show a reduced dependence of the friction coefficients when the Reynolds number decreases. However, the Reich & Beer (1989) experiment was at the same Reynolds number thus a further reason for the difference could be the influence of the inlet conditions on the measurements. This is a more plausible reason because Murakami & Kikuyama (1980) showed that the loss coefficient depends on the axial location of the measurements. Perhaps U4-03d32 files image150.gif, the distance where Reich & Beer (1989) took the measurements, is still not sufficient to achieve the condition of fully developed rotating pipe. In addition it should be considered that the direct simulations are performed at constant mass flow rate U4-03d32 files image152.gif; in the experiment, on the other hand, usually the head of the pump is kept constant thus the flow conditions are different and this could produce effects on the friction coefficients.

U4-03d32 files image154.gif
Figure 3 Axial mean velocity normalized by the bulk velocity U4-03d32 files image156.gifas function of the wall distance, (From Feiz et al. 2003 copyright International Journal of Heat and Fluid).

In Figure 3 is plotted the axial mean velocity profile versus the wall distance for various rotation numbers (N=0,1,2). When the pipe is rotating the streamwise velocity increases near the centre and decreases near the wall. The computational mean velocity profile gradually approaches Poiseuille profile due to the stabilizing effect of the centrifugal force. There is a good agreement between the experiment by Reich and Beer (1989), the DNS by Orlandi & Fatica (1997) and the LES by Feiz et al. (2003). The axial velocity profiles, scaled with the bulk velocity, assume at U4-03d32 files image158.gif in Orlandi & Fatica (1997) is about U4-03d32 files image160.gif. A similar behaviour was observed in the experiments by Reich & Beer (1989) and in the simulation by EBN. The value of U4-03d32 files image162.gif found by Orlandi & Fatica (1997) agrees with the value found by Kikuyama et al. (1983) (U4-03d32 files image164.gif) and differs from U4-03d32 files image166.gif and U4-03d32 files image168.gif found respectively by Reich & Beer (1989) and by EBN. Figure 3 shows that the variation with N of the numerical simulations are in qualitative agreement with those measured by Reich & Beer (1989). A better agreement is obtained by scaling the experimental profiles at N>0 with the ratio between the centreline velocity of the experiment and of the simulation at N=0. Since the Reynolds number is the same the reason for the difference between the profiles, at N=0, of the numerical simulation and of Reich & Beer (1989) should be the effect of the entrance conditions in the experiment.

U4-03d32 files image170.gifU4-03d32 files image172.gif

Figure 4 Left: Axial mean velocity normalized by the bulk velocity U4-03d32 files image173.gifas function of wall distance. Right, normalized in wall units, (From Feiz et al. 2003 copyright International Journal of Heat and Fluid)

The results of the Large Eddy Simulation in Figure 4 have been obtained with the dynamic model. This, indeed, has been shown to reproduce the results of DNS more accurately (Figure 4). The axial velocities using LES wit Smagorinsky model and the simulation without subgrid-scale model are in less agreement with the DNS. The reason of this discrepancy may be the choice of the value of the constant U4-03d32 files image175.gif (here 0.15). One of the problems with the Smagorinsky model is that the appropriate value of the coefficient U4-03d32 files image176.gif is different in different flow regimes. In particular, it is zero in laminar flows, and it is attenuated near the walls compared with its value (0.15) in high Reynolds number free turbulent flows. Indeed the model works well for the isotropic turbulence. For inhomogeneous flows the Smagorinsky model is too dissipative (it transfers too much energy to the residual motion Pope 2000).

U4-03d32 files image178.jpg
Figure 5 Predictions of various two-equation models and second-order closures (LRR is the Launder et al. 1975 model) for the axial mean velocity in axially rotating pipe flow. Comparisons are made with the experimental data of Imao et al. (1996), (from Speziale et al. 2000 copyright J. Fluid Mech).

In Figure 5, the predictions of several two-equation models for the axial mean velocity in an axially rotating pipe are displayed for a non-dimensional rotation rate of N= 0:5 and compared with the experimental data of Imao et al. (1996) as well as with the predictions of several second-order closures. It is clear from these calculations that, consistent with the results derived in this paper, the U4-03d32 files image180.gif model does not respond to the rotation, incorrectly yielding a mean velocity that is independent of the rotation rate of the pipe (the nonlinear U4-03d32 files image182.gif model yields similar results which are not shown for simplicity). On the other hand, the two-dimensional explicit algebraic stress model based on the SSG model (Gatski & Speziale 1993) responds to the rotation and predicts a reasonably good axial mean velocity in agreement with experiments as shown in Figure 5. However, all of these two-equation models predict that the mean swirl velocity U4-03d32 files image184.gif is zero relative to the rotating pipe. Although the two-dimensional explicit algebraic stress model does not predict the presence of a mean swirl velocity relative to the rotating pipe, it is not that serious from an engineering standpoint. This mean swirl velocity only constitutes approximately a 15% effect, relative to the axial mean velocity.

U4-03d32 files image186.jpg
Figure 6 Predictions of the cubic algebraic models of Wallin & Johansson (1997) for the mean velocity in an axially rotating pipe compared with the experimental data ● of Imao et al. (1996). (a) Axial mean velocity and (b) mean swirl velocity relative to the inertial frame. Taken from Wallin & Johansson (1997), (from Speziale et al. 2000 copyright J. Fluid Mech).

In Figure 6, it is shown that with a three-dimensional explicit algebraic stress model it is possible to predict the presence of a mean swirl velocity along with a rotationally dependent axial mean velocity. For this purpose, calculations are taken from Wallin & Johansson (1997) where comparisons were made with the experimental data of Imao et al. (1996).

U4-03d32 files image188.jpg
Figure 7 Predictions of various second-order closure models for the axial mean velocity and mean swirl velocity relative to the rotating pipe in axially rotating pipe flow. Comparisons are made with the experimental data of Imao et al. (1996). (a) Axial mean velocity and (b) mean swirl velocity relative to the rotating frame, (from Speziale et al. 2000 copyright J. Fluid Mech).

The results obtained for the axial mean velocity and mean swirl velocity, relative to the rotating pipe, obtained from second-order closures are displayed in Figure 7 side by side for the axially rotating pipe (N = 0:5). It is clear from the results for the mean velocity that the Launder et al. model is able to predict both mean effects reasonably well in an axially rotating pipe while the SSG model yields results that are even better. It should be remembered that second-order closures are equivalent to an algebraic model with a quartic nonlinearity when the algebraic stress model approximation is made in three dimensions (Gatski & Speziale 1993). This explains why both three-dimensional explicit algebraic stress models and second-order closure models do well.

U4-03d32 files image190.gif
Figure 8 Averaged streamwise velocity profiles in wall units: Orlandi & Fatica (1997): U4-03d32 files image191.gif U4-03d32 files image193.gif U4-03d32 files image194.gif N=0.5, U4-03d32 files image196.gif N=1.0, U4-03d32 files image198.gif N=2.0, ▲ N=0, × Eggels et al (1994), (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

In Figure 8 the mean velocity profiles, in wall units, obtained with the DNS by Orlandi & Fatica (1997), show the drag reduction through the upward shifting of the log law. As usual, U4-03d32 files image200.gif indicates the distance from the wall in wall units, in the pipe it is defined as U4-03d32 files image202.gif. While at N<1 the size of the buffer region remains unchanged at N=2 the buffer region almost completely disappears and a first log region starts just after the viscous region. It corresponds to the region where the turbulent energy is constant.

U4-03d32 files image204.jpg
Figure 9 Profiles of rms velocity fluctuations for N=0, ○ U4-03d32 files image206.gif; □ U4-03d32 files image208.gif; ▲ U4-03d32 files image210.gif. Lines present simulations, Full symbols Eggels et al. (1994), open symbols Durst et al. (1995), (from Orlandi 1997, copyright Physics of Fluids).

Orlandi & Fatica (1997) before analysing the rotating cases, tested the numerical method for a non-rotating pipe. The rms agrees closely with those of Eggels et al. (1994) and with those of Durst Jovanovic & Sender (1995) (Figure 9).

U4-03d32 files image212.jpg U4-03d32 files image214.jpg

Figure 10 Root Mean Square profiles of azimuthal, radial and axial velocity for N=1 (left) and of axial velocity for varius N (right) , (From Feiz et al. 2003 copyright International Journal of Heat and Fluid).

Figure 10 illustrates the rms velocity profiles for N=1 in the LES by Feiz et al. (2003) and in the DNS by Orlandi & Fatica (1997). It can be observed that the rotation of the wall has large effects on the rms, these being more pronounced for streamwise rms velocity. Similar observations have been reported by Eggels et al. (1994). The fluctuating velocity components with both subgrid scale models are quite close to the DNS. However, the numerical results, obtained using the dynamic model, are slightly better. The streamwise rms velocity of Feiz et al. (2003) and of Orlandi & Fatica (1997) is shown for various rotation numbers N (Figure 10b). Near the wall, the peak is reduced when N increases. For N=1 and 2, the distributions tend to have the same values. Orlandi & Fatica (1997) claimed that this is a tendency towards isotropization of turbulence when the rotation rate increases. For all the rotational rates the dynamic model gives results closer to DNS.

U4-03d32 files image216.jpg

Figure 11 Radial velocity in the laboratory reference frame U4-03d32 files image191.gifN=0.5, U4-03d32 files image217.gifN=1.0, U4-03d32 files image218.gifN=2.0, U4-03d32 files image220.gif theoretical U4-03d32 files image222.gif○ N=0.5, ▲ N=1, □ N=2 symbols experiment by Kikuyama et al. (1983), (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

When U4-03d32 files image224.gif is scaled with the wall rotational velocity U4-03d32 files image226.gif Reich & Beer (1989) claimed that it is independent of N and of Re. The same outcome was found by Kikuyama et al. (1983) in their large number of experiments. In the DNS by Orlandi & Fatica (1997), performed in the rotating frame, to have the same scaling as in the experiment the linear solid body rotation velocity U4-03d32 files image228.gif must be added to U4-03d32 files image230.gif. Moreover, the dimensional U4-03d32 files image232.gif should be scaled with the centre line velocity. By this operation the profiles of U4-03d32 files image234.gif (Figure 11) show a slight dependence on N. In Figure 11 there are also the values by Kikuyama at Re=10000 reported in the paper by Hirai et al. (1988). In the original paper by Kikuyama et al. (1983) it was difficult to read the data, and those at N=0.5 were given only for few radial locations. A very careful observation of the data at N=0.5, however, shows that the symbols do not exactly coincide with those at N=1. The values of Reich & Beer (1989) were not reported since these coincided exactly with the theoretical curve U4-03d32 files image236.gif

The profiles of U4-03d32 files image238.gif by Orlandi & Fatica (1997) do not collapse on a single curve, as in Reich & Beer (1989), however they are in a good agreement with Kikuyama’s data. Also EBN in the direct simulation, for two values of N, found a not perfect collapse of the profiles. On the contrary, the EBN profiles collapsed in the higher Reynolds number LES simulations at N=0.71. From these considerations we could conjecture that the discrepancies with the experimental results of Reich & Beer (1989) could be due, once more, to the effects of the entrance condition and not to Re, since in the experiments the independence of Re was found.

We checked that the lack of similarity, in Fig.3 did not depend on the insufficient time integration to reach the steady state. Moreover the averages were done with fields at U4-03d32 files image240.gifconstant. The independence of the radial resolution was also investigated, simulations with 49 and 97 points in the radial directions produced the same U4-03d32 files image242.gif profiles. The check that our U4-03d32 files image244.gif at the steady state satisfies the relationship U4-03d32 files image246.gif confirmed the assertion in EBN, that near the wall the U4-03d32 files image247.gif profile must go to zero with a zero slope. Eq.(5) shows that if the theoretical relationship U4-03d32 files image249.gif holds even close to the wall, U4-03d32 files image251.gif is different from zero at the wall.

The Kikuyama measurements at N=0.5 show for y<0.2 a tendency towards the linear profile, even if the measurement point closest to the wall is on the theoretical curve. Our believe is that the measurements near the wall are difficult when the wall rotates, as it is confirmed by some inconsistency of the data which does not appear in the central region. Inside the channel at N=0.5 the data of Kikuyama and the present ones do not coincide with the theoretical curve. By increasing N our results smoothly approach the theoretical curve while the Kikuyama data at N=1 and N=2 lie on the theoretical curve. The smooth increase of U4-03d32 files image253.gifwith N of the numerical simulation can be explained as a viscous correction that is greater as smaller is N.

6.1 ACCURACY CHECK

Finite-differences permit in a very simple way to handle free-slip and no-slip boundary conditions and thus to verify whether the non-uniform spacing maintains the energy conservation in the in-viscid limit. This check was done by performing a coarse simulation (65 39 33) for U4-03d32 files image255.gifand with a free-slip boundary. The calculation was started from a viscous unresolved simulation at Re=4900, corresponding to U4-03d32 files image257.gif, and it was let free to evolve to a statistical steady state. In this case, the mean pressure gradient drops to zero and the mean axial velocity profile tends to be constant across the pipe.

U4-03d32 files image259.jpg

Figure 12 a) Time history of total and turbulent energies; b) profiles of the normal stresses of (U4-03d32 files image260.gif i=z), (U4-03d32 files image191.gif i=r) and (U4-03d32 files image261.gifU4-03d32 files image263.gif) velocity components, (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

Figure 12a shows that the total energy remains constant and that the total turbulent energy decreases due to the reduction of turbulent energy production caused by the diminishing of the mean shear.

Figure 12b, in addition, shows that the radial profiles of each of the normal turbulent stresses reach the condition of almost isotropic turbulence within a large part of the pipe. The U4-03d32 files image265.gif assumption is the cause of the anisotropy near the free-slip wall. This check proves that second order centred finite-differences, for the nonlinear terms, do not introduce any spurious numerical viscosity, which for positive values could lead to unphysical dissipation and for negative values to unphysical production of energy. The influence of aliasing errors for finite difference methods is not as important as for pseudo-spectral methods. Kravchenko & Moin (1996) have analyzed these errors for different schemes for the non-linear terms and they have shown that aliasing errors are negligible when staggered variables are used. We recall that, truncation and aliasing errors are very important for large eddy simulations but not for direct simulations, where all the significant scales must be resolved.

U4-03d32 files image267.jpg

Figure 13 Mean axial velocity profiles in wall coordinates, (U4-03d32 files image268.gif log law), (U4-03d32 files image191.gif 33 × 39 × 17), (U4-03d32 files image269.gif 65 × 39 × 33), (U4-03d32 files image270.gif 65 × 39 × 65), (▲ 129 × 49 × 129), (▼129 × 97 × 257) (× Eggels et al. 1994), (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

A test of the accuracy of the viscous terms discretization and of the entire numerical method was performed by a grid refinement check. In the non-rotating case this check was done, starting from a very coarse grid (32 38 16) and by doubling the number of points in the U4-03d32 files image272.gif and U4-03d32 files image274.gif directions until reaching the grid (128 96 128). A non-uniform grid in r with an enhanced clustering near the wall, located the first point at U4-03d32 files image276.gif. The present results are compared with those by EUW, which were validated by comparison with PIV and LDA measurements and with the direct simulation of KMM in a plane channel.

For all grids Figure 13a-b show that some sort of logarithmic profile is achieved, but the values of the von Karman constant are function of the resolution. In Figure 13a the full profile is shown, while in Figure 13b only the part within the buffer and the log regions is shown to emphasize the differences with the log law U4-03d32 files image278.gifand with the numerical results by EUW.

By <> we denote space averages in the two homogeneous directions z and U4-03d32 files image280.gifand time averages every 10 time steps for a long physical time. The averaging time depends on the case studied, for N=0 the averages were performed for 200 time units, and this time was found to be sufficient. The 128 48 128 simulation, which uses fewer points than the EUW simulation, is in a very good agreement with their results. The good agreement is due to the non-uniform grid, which allows 48 points to be sufficient to locate the first point closer to the wall than the 96 equidistant points in EUW.

From physical arguments the changes in the log law should be expected; the decrease of the number of points causes a worse resolution of the near-wall stream-wise vortices. Since these vortices are responsible for the turbulent wall friction, a decrease in their strength produces a drag reduction.

U4-03d32 files image282.jpg

Figure 14 Second order turbulence statistics,(U4-03d32 files image191.gif 33 × 39 × 17), (U4-03d32 files image283.gif 65 × 39 × 33), (U4-03d32 files image284.gif65 × 39 × 65),(▲ 129 × 49 × 129),(129 × 97 × 257) × Eggels et al. 1994), (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

This effect is reflected in Figure 14 by the rms profiles showing a decrease of the U4-03d32 files image286.gif and U4-03d32 files image288.gif, and an increase of the U4-03d32 files image290.gif. The increase of the fluctuations in the stream-wise direction is due to a reduction on the mean velocity centreline. Figure 14a-d show that the refinement in U4-03d32 files image292.gif has a greater influence than in z. The physical reason is that the better resolution of the azimuthal gradients of U4-03d32 files image294.gif brings to a better description of the sweep and ejection events and hence to greater values of U4-03d32 files image296.gif

The observation that stream-wise vorticity produces the high- and low-speed streaks, and thus the turbulent wall friction was proposed by KMM based on observations made in the plane channel direct simulation, and was also proved by Orlandi & Jimenez (1994) by a two-dimensional model. The results in Figure 14 show that the refinement in the axial direction is less important, in fact the profiles by 128 48 64 and by 128 48 128 points do not largely differ from those by EUW by 128 96 256 points. In Figure 14d the Reynolds stress profiles of U4-03d32 files image298.gif confirm that the coarsest simulation gives a sort of drag reduction. The reduction of U4-03d32 files image300.gif is mainly due to a decrease of the U4-03d32 files image301.gif as in all situations in which a drag reduction is achieved. For the rotating case the check of the numerics was performed only at N=2. In the rotating case the check must be done not only on the resolution but also on the length of the computational domain in the stream-wise direction. In fact, if the length of the pipe was too short it could affect the long helical structures near the wall and at the centre of pipe. The near wall structures are smaller than those at the centre and then have shorter time scales. The statistical steady state is obtained with few fields near the wall, on the other hand it is necessary a larger number of fields to have satisfactory profiles in the central part. From these arguments it should be expected that the grid resolution and the pipe length should effect more the central than the wall region. We performed simulations on a pipe of length U4-03d32 files image303.gif by a grid 128 48 128 and 128 96 128 to check the radial resolution. Since a non-uniform grid was used this check is not so important, as it was shown for N=0. Two further simulations with U4-03d32 files image305.gif by 128 128 128 and 128 128 256 points were performed to investigate the effects of the length and the grid resolution in the stream-wise direction. From these simulations it was decided to use the simulation at N=2 by a grid 128 96 257 in a pipe of length U4-03d32 files image307.gif to analyze the results. With 257 points in z, the resolution in wall units is U4-03d32 files image309.gifsmall enough to consider the present simulation as a fully direct simulation even in the central region.

U4-03d32 files image311.jpg

Figure 15 Second order turbulence statistics,(U4-03d32 files image191.gif 129 × 49 × 129 U4-03d32 files image313.gif), (U4-03d32 files image314.gif 129 × 97 × 129 U4-03d32 files image316.gif), (U4-03d32 files image317.gif 129 × 129 × 129 U4-03d32 files image319.gif), (U4-03d32 files image320.gif 257 × 129 × 129 U4-03d32 files image322.gif), (U4-03d32 files image323.gif 257 × 97 × 129 U4-03d32 files image325.gif), (from Orlandi & Fatica 1997 copyright J. Fluid Mech.).

Figure 15a-d show that when the pipe is rotating the length of the pipe is an important parameter in the direct simulations.

The final comment of the grid refinement checks is that very coarse calculations in 3D give acceptable results to reproduce some of the aspects of wall turbulent flows without introducing any sub-grid scale (s.g.s) model. Although a very coarse simulation without s.g.s. model could be more or less useful, from a physical point of view accurate LES simulations are always more necessaries.

However, since this numerical viscosity produces energy spectra decreasing reasonably well at high wave numbers, it is possible that these truncation errors affect the flow in the same manner as s.g.s. models do.

© copyright ERCOFTAC 2004



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


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References