No Arabic abstract
We report the results of a complete modal and nonmodal linear stability analysis of the electrohydrodynamic flow (EHD) for the problem of electroconvection in the strong injection region. Convective cells are formed by Coulomb force in an insulating liquid residing between two plane electrodes subject to unipolar injection. Besides pure electroconvection, we also consider the case where a cross-flow is present, generated by a streamwise pressure gradient, in the form of a laminar Poiseuille flow. The effect of charge diffusion, often neglected in previous linear stability analyses, is included in the present study and a transient growth analysis, rarely considered in EHD, is carried out. In the case without cross-flow, a non-zero charge diffusion leads to a lower linear stability threshold and thus to a more unstable low. The transient growth, though enhanced by increasing charge diffusion, remains small and hence cannot fully account for the discrepancy of the linear stability threshold between theoretical and experimental results. When a cross-flow is present, increasing the strength of the electric field in the high-$Re$ Poiseuille flow yields a more unstable flow in both modal and nonmodal stability analyses. Even though the energy analysis and the input-output analysis both indicate that the energy growth directly related to the electric field is small, the electric effect enhances the lift-up mechanism. The symmetry of channel flow with respect to the centerline is broken due to the additional electric field acting in the wall-normal direction. As a result, the centers of the streamwise rolls are shifted towards the injector electrode, and the optimal spanwise wavenumber achieving maximum transient energy growth increases with the strength of the electric field.
Wakes of aircraft and automobiles with relatively flat slanted aftbodies are often characterized by a streamwise-oriented vortex pair, whose strength affects drag and other crucial performance parameters. We examine the stability characteristics of the vortex pair emerging over an abstraction comprised of a streamwise-aligned cylinder terminated with an upswept plane. The Reynolds number is fixed at 5000 and the upsweep angle is increased from 20deg to 32deg. At 20deg, the LES yields a steady streamwise-oriented vortex pair, and the global modes are also stable. At 32deg, the LES displays unsteady flow behavior. Linear analysis of the mean flow reveals different unstable modes. The lowest oscillation frequency is an antisymmetric mode, which is attached to the entire slanted base. At the highest frequency, the mode is symmetric and has the same rotational orientation as the mean vortex pair. Its support is prominent in the rear part of the slanted base and spreads relatively rapidly downstream with prominent helical structures. A receptivity analysis of low- and high-frequency modes suggests the latter holds promise to affect the vortical flow, providing a potential starting point for a control strategy to modify the vortex pair.
The importance of fluid-elastic forces in tube bundle vibrations can hardly be over-emphasized, in view of their damaging potential. In the last decades, advanced models for representing fluid-elastic coupling have therefore been developed by the community of the domain. Those models are nowadays embedded in the methodologies that are used on a regular basis by both steam generators providers and operators, in order to prevent the risk of a tube failure with adequate safety margins. From an R&D point of view however, the need still remains for more advanced models of fluid-elastic coupling, in order to fully decipher the physics underlying the observed phenomena. As a consequence, new experimental flow-coupling coefficients are also required to specifically feed and validate those more sophisticated models. Recent experiments performed at CEA-Saclay suggest that the fluid stiffness and damping coefficients depend on further dimensionless parameters beyond the reduced velocity. In this work, the problem of data reduction is first revisited, in the light of dimensional analysis. For single-phase flows, it is underlined that the flow-coupling coefficients depend at least on two dimensionless parameters, namely the Reynolds number $Re$ and the Stokes number $Sk$. Therefore, reducing the experimental data in terms of the compound dimensionless quantity $V_r=Re/Sk$ necessarily leads to impoverish results, hence the data dispersion. In a second step, experimental data are presented using the dimensionless numbers $Re$ and $Sk$. We report experiments, for a 3x5 square tube bundle subjected to water transverse flow. The bundle is rigid, except for the central tube which is mounted on a flexible suspension allowing for translation motions in the lift direction.
Since their development in 2001, regularised stokeslets have become a popular numerical tool for low-Reynolds number flows since the replacement of a point force by a smoothed blob overcomes many computational difficulties associated with flow singularities (Cortez, 2001, textit{SIAM J. Sci. Comput.} textbf{23}, 1204). The physical changes to the flow resulting from this process are, however, unclear. In this paper, we analyse the flow induced by general regularised stokeslets. An explicit formula for the flow from any regularised stokeslet is first derived, which is shown to simplify for spherically symmetric blobs. Far from the centre of any regularised stokeslet we show that the flow can be written in terms of an infinite number of singularity solutions provided the blob decays sufficiently rapidly. This infinite number of singularities reduces to a point force and source dipole for spherically symmetric blobs. Slowly-decaying blobs induce additional flow resulting from the non-zero body forces acting on the fluid. We also show that near the centre of spherically symmetric regularised stokeslets the flow becomes isotropic, which contrasts with the flow anisotropy fundamental to viscous systems. The concepts developed are used to { identify blobs that reduce regularisation errors. These blobs contain regions of negative force in order to counter the flows produced in the regularisation process, but still retain a form convenient for computations.
We explore the effect of forcing on the linear shear flow or plane Couette flow, which is also the background flow in the very small region of the Keplerian accretion disk. We show that depending on the strength of forcing and boundary conditions suitable for the systems under consideration, the background plane shear flow and, hence, accretion disk velocity profile modifies to parabolic flow, which is plane Poiseuille flow or Couette-Poiseuille flow, depending on the frame of reference. In the presence of rotation, plane Poiseuille flow becomes unstable at a smaller Reynolds number under pure vertical as well as threedimensional perturbations. Hence, while rotation stabilizes plane Couette flow, the same destabilizes plane Poiseuille flow faster and forced local accretion disk. Depending on the various factors, when local linear shear flow becomes Poiseuille flow in the shearing box due to the presence of extra force, the flow becomes unstable even for the Keplerian rotation and hence turbulence will pop in there. This helps in resolving a long standing problem of sub-critical transition to turbulence in hydrodynamic accretion disks and laboratory plane Couette flow.
We perform a three-dimensional, short-wavelength stability analysis on the numerically simulated two-dimensional flow past a circular cylinder for Reynolds numbers in the range $50le Rele300$; here, $Re = U_{infty}D/ u$ with $U_infty$, $D$ and $ u$ being the free-stream velocity, the diameter of the cylinder and the kinematic viscosity of the fluid, respectively. For a given $Re$, inviscid local stability equations from the geometric optics approach are solved on three distinct closed fluid particle trajectories (denoted as orbits 1, 2 & 3) for purely transverse perturbations. The inviscid instability on orbits 1 & 2, which are symmetric counterparts of one another, is shown to undergo bifurcations at $Reapprox50$ and $Reapprox250$. Upon incorporating finite-wavenumber, finite-Reynolds number effects to compute corrected local instability growth rates, the inviscid instability on orbits 1 & 2 is shown to be suppressed for $Relesssim262$. Orbits 1 & 2 are thus shown to exhibit a synchronous instability for $Regtrsim262$, which is remarkably close to the critical Reynolds number for the mode-B secondary instability. Further evidence for the connection between the local instability on orbits 1 & 2, and the mode-B secondary instability, is provided via a comparison of the growth rate variation with span-wise wavenumber between the local and global stability approaches. In summary, our results strongly suggest that the three-dimensional short-wavelength instability on orbits 1 & 2 is a possible mechanism for the emergence of the mode B secondary instability.