Sensitivity of the thresholds to unsteadiness upon confinement shape in deformed driven cavity
A. Redondo, G. Kasperski and G. Labrosse
Contents |
Introduction
The lid driven cavity problem in non orthogonal geometries is considered with two purposes: one is to develop numerical tools for obtaining spectrally accurate flows in such geometries and then to point out the sensitivity upon the boundary deformations of the transition thresholds to unsteadiness. Accurate solutions of the Navier Stokes equations are indeed obtained thanks to a spectral Chebyshev discretization [2] preconditioned with Finite Volumes [3]. Then, the influence of the deformation upon flow topology and stability properties is shown.
Numerical Method
The Navier-Stokes equations
,
are solved for the incompressible flow of a Newtonian fluid (of kinematic viscosity and density ) confined in deformed driven cavities of height . Taking and as length and velocity scales, the Reynolds number is defined as , being the velocity amplitude of the moving lid. In order to suppress the singularities at the upper corners, the velocity of the lid is regularized according to where is the length of the upper lid and , the coordinate along the upper lid. The lateral deformations are of parabolic type, described by where is the amplitude (positive or negative) of the deformation.
Temporal discretization is based on second-order explicit Adams Bashforth scheme for the non linear terms and an implicit backward Euler for linear terms.
A Projection-Diffusion algorithm is used to uncouple the pressure and velocity fields [1,4], which leads to solve:
a Poisson equation for the pressure,
with Neumann boundary conditions:
where stands for the nonlinear term assumed to be explicitely evaluated in time.
An Helmholtz problem is then solved for each velocity component :
with Dirichlet boundary conditions.
The remaining Poisson and Helmholtz problems are treated using spectral Chebyshev spatial approximation.
Spatial derivatives are evaluated using a Gauss-Lobatto collocation method (), and the time step is fixed to . For each fixed deformation, the Reynolds number is recursively incremented. The criterion for steadiness, based on the vorticity field , is where .
Determination of transition thresholds
These computations led to estimates of the transition thresholds from steady to oscillatory flows for boundary deformations . Results are reported on Fig.1 : transition thresholds lie in between a lower limit corresponding to a steady state and an upper limit for which the flow is oscillatory.
These first results show the sensitivity of the steady flow stability to boundary deformations and a strongly non monotonic behaviour of the thresholds.
The graphic can be read considering three different zones, approximatively defined by , and .
The first zone corresponds to the pinched configuration wherein the critical Reynolds number value is significantly lowered, from 10000 for a straight cavity to 5000. In the second zone, corresponding to straight or slightly deformed cavities, the thresholds are only weakly modified. Refining their determination is necessary to characterize the local behaviour about . The inflated deformation also significantly reduces the threshold values, from 10000 to 7000. In both first and third zones, a non monotonic sensitivity of the thresholds upon is observed, quite clearly for the inflated cavities.
Flow topology
A correlation exists between what has just been described and the flow topologies. The streamfunction is presented for the last observed steady flows, before transition, in Fig.2
The straight cavity is taken as a reference. The flow is made of a big roll driven by the upper
lid, two counter-rotating cells in the left-down (LD) and
right-down (RD) corners, a counter rotating cell close to the
left-up (LU) corner and recirculating Moffatt corner cells.
The pinching deformations tend to stretch the RD cell until it merges the LD cell. The value associated with this junction is located between points C and D of Fig.1 where there is an abrupt evolution of the thresholds.
In a same way, inflating the cavity tends to stretch the LD cell until it reaches the LU cell. The junction is located between points F and G of Fig.1 and seems to be linked to the second abrupt behaviour of the thresholds.
Thus, each zone of Fig.1 corresponds to a different topology of the flow, mainly composed of horizontally organized cells for , and vertically organized ones for . A non-monotonic behaviour in each configuration seems linked to the appearance of new cells (e.g. in cases and ).
References
[1] A. Batoul, H. Khallouf, and G. Labrosse. Une méthode de Résolution Directe (Pseudo-Spectrale) du problème de Stokes 2D/3D Instationnaire. Application à la Cavité Entraînée Carrée. C.R. Acad. Sci. Paris, 319(I):1455-1461, 1994.
[2] C. Canuto, M.Y. Hussaini, A. Quarteroni, and T.A. Zang. Spectral Methods un Fluid Dynamics. Springer Series in Computational Physics. Springer-Verlag, New York, 1988.
[3] V. Delgado, G. Kasperski, and G. Labrosse. Spectrally accurate numerical incompressible flows with deformed boundaries: a finite volume preconditioning method. J. Comput. and Applied Math., 168:145-154,2004.
[4] E. Leriche and G. Labrosse. High-order direct Stokes solvers with or without temporal splitting : numerical investigations of their comparative properties. SIAM J. Scient. Comput., 22(4):1386-1410,2000.