No Arabic abstract
This paper numerically investigates the shear flow between double concentric spherical boundaries rotating differentially, so-called spherical Couette flow, under unstable thermal stratification, focusing on the boundary of the axisymmetric/non-axisymmetric transition in wide gap cases where the inner radius is comparable to the clearance width. While the transition of SCF has been confirmed experimentally in cases without thermal factor, insufficient knowledge on SCF subject to thermal instability, related to geophysical problems especially in wide gap cases, has been accumulated mainly based on numerical analysis; our motivation is to bridge the knowledge gap by a parameter extension. We reconfirm that the transition under no thermal effect is initiated by a disturbance visualised as a spiral pattern with n arms extending from the equatorial zone to the pole in each hemisphere, at the critical Reynolds number, Recr, as previously reported. With increasing thermal factor, the buoyancy effect assists the system rotation to trigger a transition towards non-axisymmetric states, resulting in a relative decrease of Recr. This is in contrast with the result that the system rotation apparently suppresses via Coriolis effect the transition to the thermally convective states at low Reynolds numbers. The present study elucidates that the existence of the axisymmetric state is restricted within a closed area in the extended parameter space, along the boundary of which the spiral patterns observed experimentally in SCF continually connect to the classical spherical Benard convective states.
Equilibrium, traveling wave, and periodic orbit solutions of pipe, channel, and plane Couette flows can now be computed precisely at Reynolds numbers above the onset of turbulence. These invariant solutions capture the complex dynamics of wall-bounded rolls and streaks and provide a framework for understanding low-Reynolds turbulent shear flows as dynamical systems. We present fluid dynamics videos of plane Couette flow illustrating periodic orbits, a close pass of turbulent flow to a periodic orbit, and heteroclinic connections between unstable equilibria.
Motivated by recent experimental and numerical studies of coherent structures in wall-bounded shear flows, we initiate a systematic exploration of the hierarchy of unstable invariant solutions of the Navier-Stokes equations. We construct a dynamical, 10^5-dimensional state-space representation of plane Couette flow at Re = 400 in a small, periodic cell and offer a new method of visualizing invariant manifolds embedded in such high dimensions. We compute a new equilibrium solution of plane Couette flow and the leading eigenvalues and eigenfunctions of known equilibria at this Reynolds number and cell size. What emerges from global continuations of their unstable manifolds is a surprisingly elegant dynamical-systems visualization of moderate-Reynolds turbulence. The invariant manifolds tessellate the region of state space explored by transiently turbulent dynamics with a rigid web of continuous and discrete symmetry-induced heteroclinic connections.
While linear non-normality underlies the mechanism of energy transfer from the externally driven flow to the perturbation field that sustains turbulence, nonlinearity is also known to play an essential role. The goal of this study is to better understand the role of nonlinearity in sustaining turbulence. The method used in this study is implementation in Couette flow of a statistical state dynamics (SSD) closure at second order in a cumulant expansion of the Navier-Stokes equations in which the averaging operator is the streamwise mean. The perturbations are the deviations from the streamwise mean and two mechanisms potentially contributing to maintaining these perturbations are identified. These are parametric perturbation growth arising from interaction of the perturbations with the fluctuating mean flow and transient growth of perturbations arising from nonlinear interaction between components of the perturbation field. By the method of comparing the turbulence maintained in the SSD and in the associated direct numerical simulation (DNS) in which these mechanisms have been selectively included and excluded, parametric growth is found to maintain the perturbation field of the turbulence while the more commonly invoked mechanism of transient growth of perturbations arising from scattering by nonlinear interaction is found to suppress perturbation growth. In addition to verifying that the parametric mechanism maintains the perturbations in DNS it is also verified that the Lyapunov vectors are the structures that dominate the perturbation energy and energetics in DNS. It is further verified that these vectors are responsible for maintaining the roll circulation that underlies the self-sustaining process (SSP) and in particular the maintenance of the fluctuating streak that supports the parametric perturbation growth.
In order to explore the magnetostrophic regime expected for planetary cores, experiments have been conducted in a rotating sphere filled with liquid sodium, with an imposed dipolar magnetic field (the DTS setup). The field is produced by a permanent magnet enclosed in an inner sphere, which can rotate at a separate rate, producing a spherical Couette flow. The flow properties are investigated by measuring electric potentials on the outer sphere, the induced magnetic field in the laboratory frame, and velocity profiles inside the liquid sodium using ultrasonic Doppler velocimetry. The present article focuses on the time-averaged axisymmetric part of the flow. The Doppler profiles show that the angular velocity of the fluid is relatively uniform in most of the fluid shell, but rises near the inner sphere, revealing the presence of a magnetic wind, and gently drops towards the outer sphere. The transition from a magnetostrophic flow near the inner sphere to a geostrophic flow near the outer sphere is controlled by the local Elsasser number. For Rossby numbers up to order 1, the observed velocity profiles all show a similar shape. Numerical simulations in the linear regime are computed, and synthetic velocity profiles are compared with the measured ones. In the geostrophic region, a torque-balance model provides very good predictions. We find that the induced magnetic field varies in a consistent fashion, and displays a peculiar peak in the counter-rotating regime. This happens when the fluid rotation rate is almost equal and opposite to the outer sphere rotation rate. The fluid is then almost at rest in the laboratory frame, and the Proudman-Taylor constraint vanishes, enabling a strong meridional flow. We suggest that dynamo action might be favored in such a situation.
Numerical simulation of Electroconvective vortices behavior in the presence of Couette flow between two infinitely long electrodes is investigated. The two-relaxation-time Lattice Boltzmann Method with fast Poisson solver solves for the spatiotemporal distribution of flow field, electric field, and charge density. Couette cross-flow is applied to the solutions after the electroconvective vortices are established. Increasing cross-flow velocity deforms the vortices and eventually suppresses them when threshold values of shear stress are reached.