Coriolis induced compressibility effects in
rotating shear layers
Bernard J. Geurts1,2, Darryl D. Holm3,4, and Arkadiusz K. Kuczaj1
1
Multiscale Modeling and Simulation, NACM, J.M. Burgers Center, Faculty EEMCS, University of Twente, P.O. Box 217, 7500 AE Enschede, the Netherlands: b.j.geurts@utwente.nl
2
Anisotropic Turbulence, Fluid Dynamics Laboratory, Department of Applied Physics, Eindhoven University of Technology, P.O. Box 513, 5300 MB Eindhoven, the Netherlands
3
Mathematics Department, Imperial College London, SW7 2AZ, London, UK
4
CCS2, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
Rotation about a fixed axis introduces a competition between two- and three-dimensional tendencies in a turbulent flow [1, 2, 3]. At strong rotation rates, this competition expresses itself, e.g., by suppressing fluid motion along the axis of rotation. This yields a reduced mixing-efficiency in a rotating frame of reference that is particularly relevant in the context of atmospheric and oceanic flows. The exchange of gases between atmosphere and oceans, the transport of heat and the spreading of pollutants or large-scale plankton-populations, are all significantly affected by rotation in large-scale environ-mental flows.
The canonical problem of flow in a horizontal temporal mixing layer, sub-jected to rotation about a vertical axis is investigated with the use of direct numerical simulation. The compressible mixing layer is considered, composed of two counter-flowing horizontal slabs of fluid [4]. Periodic boundary condi-tions are adopted in the horizontal direccondi-tions and free-slip condicondi-tions apply in the vertical far-field. High order finite volume discretization and explicit, compact storage, four-stage Runge-Kutta time-stepping are used.
Qualitatively, in the non-rotating case the flow shows vigorous mixing in the vertical direction leading to a complex three-dimensional flow. A snap-shot of the vertical vorticity component in the turbulent regime is shown in Fig. 1(a). Adding rotation to this flow may completely alter the transitional and developed flow. As a result of rotation, the flow breaks up into vertically oriented cells that collectively undergo an oscillatory motion (cf. Fig. 1(b)). The frequency of this oscillation increases with rotation rate.
The effect of rotation was investigated at a Reynolds number Re = 50, based on the initial momentum-thickness of the shear layer, and a low convec-tive Mach number of M = 0.2. The strength of the rotation is characterized
2 Bernard J. Geurts, Darryl D. Holm, and Arkadiusz K. Kuczaj
(a) (b)
Fig. 1.Snapshot of the vertical component of vorticity in the developed regime at t = 90 for (a) the non-rotating case Ro = ∞ and (b) rotating flow at Ro = 10. Positive and negative values are shown in different colors.
by the Rossby number Ro = ur/(2Ωrℓr) where ur, Ωr and ℓr are reference
velocity, rotation-rate and length-scale used in the non-dimensionalization. At high Rossby numbers Ro, i.e., low rotation rates, the flow develops very similarly to incompressible flow. However, at sufficiently low Rossby numbers Ro < 1, rotation was found to induce significant compressibility effects. These give rise to rapidly varying small scale flow-features. By comparing results ob-tained at a variety of time-steps, we verified that the explicit Runge-Kutta time-stepping properly captures these rapid variations. The spatial resolution was found to be adequate at 2563
grid-cells for the purpose of studying the decay of the kinetic energy. Results at 2563
grid-cells were compared with lower resolutions at 1923
and 1283
, confirming the accuracy.
Rotation has an indirect influence on the decay-rate of kinetic energy E = hu · ui/2. Here h·i denotes averaging over the flow-domain. The evolu-tion of the kinetic energy may be written as
dE
dt = D − W ; D = hp∇ · ui ; W = hS : S/(2Re)i (1) where p denotes pressure, u the velocity field, S the rate of strain tensor and Re the Reynolds number. The viscous dissipation is represented in W while the compressible pressure/velocity-divergence is given by D. In (1) there is no direct effect of rotation. Rather, rotation effects may be distinguished in the dynamics of D and W . By investigating the various contributions to dD/dt and dW/dt we may identify the basic mechanisms that govern the kinetic energy dynamics, particularly the role of the Coriolis forces. The dominant mechanisms at various Ro in the dynamics of the viscous dissipation and the pressure/velocity-divergence can be extracted from the DNS-data.
Coriolis induced compressibility effects 3 0 10 20 30 40 50 60 70 80 90 100 7.5 8 8.5 9 9.5 10x 10 4 (a) 00 10 20 30 40 50 60 70 80 90 100 1 2 3 4 5 6 7 8 9 10x 10 4 (b) E E t t
Fig. 2.Decay of kinetic energy at low (a) and high (b) rotation rates. In (a) Ro = ∞ (solid), Ro = 10 (dashed), Ro = 5 (dash-dotted); in (b) Ro = 1 (solid), Ro = 0.5 (dashed) and Ro = 0.2 (dash-dotted).
The marked differences in the decay of kinetic energy in slowly and in rapidly rotating shear layers are shown in Fig. 2. At high Ro an increase in the rotation-rate yields a decrease in the decay rate (cf. Fig. 2(a)). This reduction in decay-rate and the strict monotonicity of dE/dt are also found in incompressible flow. For sufficiently rapid rotation (in this flow Ro . 1; Fig. 2(b)) explicit compressibility effects become dominant. The behavior of the decay-rate is completely reversed and an increase in rotation-rate yields a strong increase in the ‘average decay-rate’. Moreover, the kinetic energy is no longer a monotonously decaying function of time. Next to the strictly dis-sipative viscous contributions, the pressure/velocity-divergence correlations are found to become dynamically significant at low Ro. The observed rapid fluctuations are a purely compressible effect. A complete analysis is topic of current investigations.
A.K.K. gratefully acknowledges financial support from the Dutch Foundation for Fundamental Research of Matter (FOM) and from the Turbulence Working Group at Los Alamos National Laboratory - USA.
References
1. Mansour, NN, Cambon, C, Speziale, CG (1992) in Studies in Turbulence (ed. Gatski, TB, Sarkar, S, Speziale, CG). Springer.
2. Hossain, M (1994) Phys. Fluids 6:1077.
3. Smith, LM, Waleffe, F (1999) Phys. Fluids. 6:1608. 4. Geurts, BJ, Holm, DD (2006) J. of Turbulence. 7 1 - 33.