scholarly journals Modal and non-modal stability analysis of electrohydrodynamic flow with and without cross-flow

2015 ◽  
Vol 770 ◽  
pp. 319-349 ◽  
Author(s):  
Mengqi Zhang ◽  
Fulvio Martinelli ◽  
Jian Wu ◽  
Peter J. Schmid ◽  
Maurizio Quadrio

We report the results of a complete modal and non-modal linear stability analysis of the electrohydrodynamic flow for the problem of electroconvection in the strong-injection region. Convective cells are formed by the 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 electrohydrodynamics, 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 flow. 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-$\mathit{Re}$Poiseuille flow yields a more unstable flow in both modal and non-modal 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 centreline is broken due to the additional electric field acting in the wall-normal direction. As a result, the centres 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.

1993 ◽  
Vol 252 ◽  
pp. 209-238 ◽  
Author(s):  
Satish C. Reddy ◽  
Dan S. Henningson

In recent work it has been shown that there can be substantial transient growth in the energy of small perturbations to plane Poiseuille and Couette flows if the Reynolds number is below the critical value predicted by linear stability analysis. This growth, which may be as large as O(1000), occurs in the absence of nonlinear effects and can be explained by the non-normality of the governing linear operator - that is, the non-orthogonality of the associated eigenfunctions. In this paper we study various aspects of this energy growth for two- and three-dimensional Poiseuille and Couette flows using energy methods, linear stability analysis, and a direct numerical procedure for computing the transient growth. We examine conditions for no energy growth, the dependence of the growth on the streamwise and spanwise wavenumbers, the time dependence of the growth, and the effects of degenerate eigenvalues. We show that the maximum transient growth behaves like O(R2), where R is the Reynolds number. We derive conditions for no energy growth by applying the Hille–Yosida theorem to the governing linear operator and show that these conditions yield the same results as those derived by energy methods, which can be applied to perturbations of arbitrary amplitude. These results emphasize the fact that subcritical transition can occur for Poiseuille and Couette flows because the governing linear operator is non-normal.


2014 ◽  
Vol 742 ◽  
pp. 254-290 ◽  
Author(s):  
Simon Maretzke ◽  
Björn Hof ◽  
Marc Avila

AbstractNon-normal transient growth of disturbances is considered as an essential prerequisite for subcritical transition in shear flows, i.e. transition to turbulence despite linear stability of the laminar flow. In this work we present numerical and analytical computations of linear transient growth covering all linearly stable regimes of Taylor–Couette flow. Our numerical experiments reveal comparable energy amplifications in the different regimes. For high shear Reynolds numbers$\mathit{Re}$, the optimal transient energy growth always follows a$\mathit{Re}^{2/3}$scaling, which allows for large amplifications even in regimes where the presence of turbulence remains debated. In co-rotating Rayleigh-stable flows, the optimal perturbations become increasingly columnar in their structure, as the optimal axial wavenumber goes to zero. In this limit of axially invariant perturbations, we show that linear stability and transient growth are independent of the cylinder rotation ratio and we derive a universal$\mathit{Re}^{2/3}$scaling of optimal energy growth using Wentzel–Kramers–Brillouin theory. Based on this, a semi-empirical formula for the estimation of linear transient growth valid in all regimes is obtained.


Author(s):  
Anirban Guha ◽  
Ian A. Frigaard

We have investigated the linear stability of plane Couette-Poiseuille flow in the presence of a cross-flow. The base flow is characterised by the cross flow Reynolds number, Ri and the dimensionless wall velocity, k. Corresponding to each k ∈ [0,1], we have observed two ranges of Ri for which the flow is unconditionally linearly stable. In the lower range, we have a stabilisation of long wavelengths leading to a cut-off Ri. In this range, cross-flow stabilisation and Couette stabilisation appear to act via very similar mechanisms in this range, leading to the potential for robust compensatory design of flow stabilisation using either mechanism. As Ri is increased, we see first destabilisation and then stabilisation at very large Ri. The instability is again a long wavelength mechanism. A linear energy analysis reveals that in this range the Reynolds stress becomes amplified, the critical layer is irrelevant and viscous dissipation is completely dominated by the energy production/negation, which approximately balances at criticality.


Author(s):  
Yuxin Liu ◽  
Chaojie Mo ◽  
Lujia Liu ◽  
Qingfei Fu ◽  
Lijun Yang

This article presents the linear stability analysis of an electrified liquid sheet injected into a compressible ambient gas in the presence of a transverse electric field. The disturbance wave growth rates of sinuous and varicose modes were determined by solving the dispersion relation of the electrified liquid sheet. It was determined that by increasing the Mach number of the ambient gas from subsonic to transonic, the maximum growth rate and the dominant wave number of the disturbances were increased, and the increase was greater in the presence of the electric field. The electrified liquid sheet was more unstable than the non-electrified sheet. The increase of both the gas-to-liquid density ratio and the electrical Euler number accelerated the breakup of the liquid sheet for both modes; while the ratio of distance between the horizontal electrode and the liquid-sheet-to-sheet thickness had the opposite effect. High Reynolds and Weber numbers accelerated the breakup of the electrified liquid sheet.


2018 ◽  
Vol 840 ◽  
pp. 5-24 ◽  
Author(s):  
Junho Park ◽  
Paul Billant ◽  
Jong-Jin Baik ◽  
Jaemyeong Mango Seo

The stably stratified Taylor–Couette flow is investigated experimentally and numerically through linear stability analysis. In the experiments, the stability threshold and flow regimes have been mapped over the ranges of outer and inner Reynolds numbers: $-2000<Re_{o}<2000$ and $0<Re_{i}<3000$, for the radius ratio $r_{i}/r_{o}=0.9$ and the Brunt–Väisälä frequency $N\approx 3.2~\text{rad}~\text{s}^{-1}$. The corresponding Froude numbers $F_{o}$ and $F_{i}$ are always much smaller than unity. Depending on $Re_{o}$ (or equivalently on the angular velocity ratio $\unicode[STIX]{x1D707}=\unicode[STIX]{x1D6FA}_{o}/\unicode[STIX]{x1D6FA}_{i}$), three different regimes have been identified above instability onset: a weakly non-axisymmetric mode with low azimuthal wavenumber $m=O(1)$ is observed for $Re_{o}<0$ ($\unicode[STIX]{x1D707}<0$), a highly non-axisymmetric mode with $m\sim 12$ occurs for $Re_{o}>840$ ($\unicode[STIX]{x1D707}>0.57$) while both modes are present simultaneously in the lower and upper parts of the flow for $0\leqslant Re_{o}\leqslant 840$ ($0\leqslant \unicode[STIX]{x1D707}\leqslant 0.57$). The destabilization of these primary modes and the transition to turbulence as $Re_{i}$ increases have been also studied. The linear stability analysis proves that the weakly non-axisymmetric mode is due to the centrifugal instability while the highly non-axisymmetric mode comes from the strato-rotational instability. These two instabilities can be clearly distinguished because of their distinct dominant azimuthal wavenumber and frequency, in agreement with the recent results of Park et al. (J. Fluid Mech., vol. 822, 2017, pp. 80–108). The stability threshold and the characteristics of the primary modes observed in the experiments are in very good agreement with the numerical predictions. Moreover, we show that the centrifugal and strato-rotational instabilities are observed simultaneously for $0\leqslant Re_{o}\leqslant 840$ in the lower and upper parts of the flow, respectively, because of the variations of the local Reynolds numbers along the vertical due to the salinity gradient.


2007 ◽  
Vol 577 ◽  
pp. 211-239 ◽  
Author(s):  
C. NOUAR ◽  
N. KABOUYA ◽  
J. DUSEK ◽  
M. MAMOU

The receptivity problem of plane Bingham–Poiseuille flow with respect to weak perturbations is addressed. The relevance of this study is highlighted by the linear stability analysis results (spectra and pseudospectra). The first part of the present paper thus deals with the classical normal-mode approach in which the resulting eigenvalue problem is solved using the Chebychev collocation method. Within the range of parameters considered, the Poiseuille flow of Bingham fluid is found to be linearly stable. The second part investigates the most amplified perturbations using the non-modal approach. At a very low Bingham number (B ≪ 1), the optimal disturbance consists of almost streamwise vortices, whereas at moderate or large B the optimal disturbance becomes oblique. The evolution of the obliqueness as function of B is determined. The linear analysis presented also indicates, as a first stage of a theoretical investigation, the principal challenges of a more complete nonlinear study.


2017 ◽  
Vol 828 ◽  
pp. 812-836 ◽  
Author(s):  
Marc A. Regan ◽  
Krishnan Mahesh

The stability of low-speed jets in cross-flow (JICF) is studied using tri-global linear stability analysis (GLSA). Simulations are performed at a Reynolds number of 2000, based on the jet exit diameter and the average velocity. A time stepper method is used in conjunction with the implicitly restarted Arnoldi iteration method. GLSA results are shown to capture the complex upstream shear-layer instabilities. The Strouhal numbers from GLSA match upstream shear-layer vertical velocity spectra and dynamic mode decomposition from simulation (Iyer & Mahesh, J. Fluid Mech., vol. 790, 2016, pp. 275–307) and experiment (Megerian et al., J. Fluid Mech., vol. 593, 2007, pp. 93–129). Additionally, the GLSA results are shown to be consistent with the transition from absolute to convective instability that the upstream shear layer of JICFs undergoes between $R=2$ to $R=4$ observed by Megerian et al. (J. Fluid Mech., vol. 593, 2007, pp. 93–129), where $R=\overline{v}_{jet}/u_{\infty }$ is the jet to cross-flow velocity ratio. The upstream shear-layer instability is shown to dominate when $R=2$, whereas downstream shear-layer instabilities are shown to dominate when $R=4$.


2019 ◽  
Author(s):  
Edwin Khoo ◽  
Hongbo Zhao ◽  
Martin Z. Bazant

We study the linear stability analysis of time-dependent electrodeposition in a charged random porous medium, whose pore surface charges can generally be of any sign, that is flanked by a pair of planar metal electrodes. Discretization of the linear stability problem results in a generalized eigenvalue problem for the dispersion relation that is solved numerically, which agrees well with the analytical approximation obtained from a boundary layer analysis valid at high wavenumbers. Under galvanostatic conditions in which an overlimiting current is applied, in the classical case of zero pore surface charges, the voltage and electric field at the cathode diverge when the bulk electrolyte concentration there vanishes at Sand's time. The same phenomenon happens for positive surface charges but at a time earlier than Sand's time. In contrast, negative surface charges allow the electrochemical system to sustain an overlimiting current via surface conduction past Sand's time, keeping the voltage and electric field bounded. Therefore, at Sand's time, negative surface charges greatly reduce the electrode surface instabilities while zero and positive surface charges magnify them. We compare theoretical predictions for overall electrode surface stabilization from the linear stability analysis with published experimental data for copper electrodeposition in cellulose nitrate membranes and demonstrate good agreement between theory and experiment. We also use the linear stability analysis as a tool to analyze how the crystal grain size changes with duty cycle during pulse electroplating.


Sign in / Sign up

Export Citation Format

Share Document