An energy criterion for the linear stability of conservative flows

2000 ◽  
Vol 402 ◽  
pp. 329-348 ◽  
Author(s):  
P. A. DAVIDSON

We investigate the linear stability of inviscid flows which are subject to a conservative body force. This includes a broad range of familiar conservative systems, such as ideal MHD, natural convection, flows driven by electrostatic forces and axisymmetric, swirling, recirculating flow. We provide a simple, unified, linear stability criterion valid for any conservative system. In particular, we establish a principle of maximum action of the formformula herewhere η is the Lagrangian displacement,e is a measure of the disturbance energy, T and V are the kinetic and potential energies, and L is the Lagrangian. Here d represents a variation of the type normally associated with Hamilton's principle, in which the particle trajectories are perturbed in such a way that the time of flight for each particle remains the same. (In practice this may be achieved by advecting the streamlines of the base flow in a frozen-in manner.) A simple test for stability is that e is positive definite and this is achieved if L(η) is a maximum at equilibrium. This captures many familiar criteria, such as Rayleigh's circulation criterion, the Rayleigh–Taylor criterion for stratified fluids, Bernstein's principle for magnetostatics, Frieman & Rotenberg's stability test for ideal MHD equilibria, and Arnold's variational principle applied to Euler flows and to ideal MHD. There are three advantages to our test: (i) d2T(η) has a particularly simple quadratic form so the test is easy to apply; (ii) the test is universal and applies to any conservative system; and (iii) unlike other energy principles, such as the energy-Casimir method or the Kelvin–Arnold variational principle, there is no need to identify all of the integral invariants of the flow as a precursor to performing the stability analysis. We end by looking at the particular case of MHD equilibria. Here we note that when u and B are co-linear there exists a broad range of stable steady flows. Moreover, their stability may be assessed by examining the stability of an equivalent magnetostatic equilibrium. When u and B are non-parallel, however, the flow invariably violates the energy criterion and so could, but need not, be unstable. In such cases we identify one mode in which the Lagrangian displacement grows linearly in time. This is reminiscent of the short-wavelength instability of non-Beltrami Euler flows.

2019 ◽  
Vol 865 ◽  
Author(s):  
Cheng-Nian Xiao ◽  
Inanc Senocak

We investigate the stability of the Prandtl model for katabatic slope flows using both linear stability theory and direct numerical simulations. Starting from Prandtl’s analytical solution for uniformly cooled laminar slope flows, we use linear stability theory to identify the onset of instability and features of the most unstable modes. Our results show that the Prandtl model for parallel katabatic slope flows is prone to transverse and longitudinal modes of instability. The transverse mode of instability manifests itself as stationary vortical flow structures aligned in the along-slope direction, whereas the longitudinal mode of instability emerges as waves propagating in the base-flow direction. Beyond the stability limits, these two modes of instability coexist and form a complex flow structure crisscrossing the plane of flow. The emergence of a particular form of these instabilities depends strongly on three dimensionless parameters, which are the slope angle, the Prandtl number and a newly introduced stratification perturbation parameter, which is proportional to the relative importance of the disturbance to the background stratification due to the imposed surface buoyancy flux. We demonstrate that when this parameter is sufficiently large, then the stabilising effect of the background stratification can be overcome. For shallow slopes, the transverse mode of instability emerges despite meeting the Miles–Howard stability criterion of $Ri>0.25$. At steep slope angles, slope flow can remain linearly stable despite attaining Richardson numbers as low as $3\times 10^{-3}$.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


2010 ◽  
Vol 656 ◽  
pp. 417-447 ◽  
Author(s):  
ANIRBAN GUHA ◽  
IAN A. FRIGAARD

We present a detailed study of the linear stability of the plane Couette–Poiseuille flow in the presence of a crossflow. The base flow is characterized by the crossflow Reynolds number Rinj and the dimensionless wall velocity k. Squire's transformation may be applied to the linear stability equations and we therefore consider two-dimensional (spanwise-independent) perturbations. Corresponding to each dimensionless wall velocity, k ∈ [0, 1], two ranges of Rinj exist where unconditional stability is observed. In the lower range of Rinj, for modest k we have a stabilization of long wavelengths leading to a cutoff Rinj. This lower cutoff results from skewing of the velocity profile away from a Poiseuille profile, shifting of the critical layers and the gradual decrease of energy production. Crossflow stabilization and Couette stabilization appear to act via very similar mechanisms in this range, leading to the potential for a robust compensatory design of flow stabilization using either mechanism. As Rinj is increased, we see first destabilization and then stabilization at very large Rinj. The instability is again a long-wavelength mechanism. An analysis of the eigenspectrum suggests the cause of instability is due to resonant interactions of Tollmien–Schlichting waves. 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. The stabilization at very large Rinj appears to be due to decay in energy production, which diminishes like Rinj−1. Our study is limited to two-dimensional, spanwise-independent perturbations.


Author(s):  
Alan Elcrat ◽  
Bartosz Protas

In this investigation, we revisit the question of linear stability analysis of two-dimensional steady Euler flows characterized by the presence of compact regions with constant vorticity embedded in a potential flow. We give a complete derivation of the linearized perturbation equation which, recognizing that the underlying equilibrium problem is of a free-boundary type, is carried out systematically using methods of shape-differential calculus. Particular attention is given to the proper linearization of contour integrals describing vortex induction. The thus obtained perturbation equation is validated by analytically deducing from it stability analyses of the circular vortex, originally due to Kelvin, and of the elliptic vortex, originally due to Love, as special cases. We also propose and validate a spectrally accurate numerical approach to the solution of the stability problem for vortices of general shape in which all singular integrals are evaluated analytically.


1972 ◽  
Vol 39 (3) ◽  
pp. 717-722 ◽  
Author(s):  
H. H. E. Leipholz

Using Galerkin’s method it is shown that in the domain of divergence, the nonconservative system of the follower-load type is always more stable than the corresponding conservative system. Hence, for nonconservative systems of the divergence type, the critical load of the corresponding conservative system becomes a lower bound for the buckling load, and the energy criterion remains sufficient for predicting stability. Moreover, it is proven that even for more general nonconservative systems, the energy criterion is sufficient under certain restrictions.


Author(s):  
Elia Merzari ◽  
Paul Fischer ◽  
W. David Pointer

Buoyancy-driven systems are subject to several types of flow instabilities. To evaluate the performance of such systems it is becoming increasingly crucial to be able to predict the stability of a given base flow configuration. Traditional Modal Linear stability Analysis requires the solution of very large eigenvalue systems for three-dimensional flows, which make this problem difficult to tackle. An alternative to modal Linear stability Analysis is the use of adjoint solvers [1] in combination with a power iteration [2]. Such methodology allows for the identification of an optimal disturbance or forcing and has been recently used to evaluate the stability of several isothermal flow systems [2]. In this paper we examine the extension of the methodology to non-isothermal flows driven by buoyancy. The contribution of buoyancy in the momentum equation is modeled through the Boussinesq approximation. The method is implemented in the spectral element code Nek5000. The test case is the flow is a two-dimensional cavity with differential heating and conductive walls and the natural circulation flow in a toroidal thermosiphon.


1982 ◽  
Vol 37 (8) ◽  
pp. 816-829
Author(s):  
E. Rebhan

A family of ideal MHD equilibria is considered introducing the concept of a driving parameter λ the increase of which beyond a certain threshold λ0 drives the plasma from a linearly stable to an unstable state. Using reductive perturbation theory, the nonlinear ideal MHD equations of motion are expanded in the neighbourhood of λ0 with respect to a small parameter ε. An appropriate scaling for the expansions is derived from the linear eigenmode problem. Integrability conditions for the reduced nonlinear equations yield nonlinear amplitude equations for the marginal mode. Nonlinearly, the instabilities are either oscillations about bifurcating equilibria, or they are explosive. In the latter case, the stability limit depends on the amplitude of the perturbation and is shifted into the linearly stable regime. Generally bifurcation of dynamically connected equilibria is observed at λ0


2012 ◽  
Vol 711 ◽  
pp. 27-39 ◽  
Author(s):  
Bo-Fu Wang ◽  
Dong-Jun Ma ◽  
Cheng Chen ◽  
De-Jun Sun

AbstractThe instabilities and transitions of flow in a vertical cylindrical cavity with heated bottom, cooled top and insulated sidewall are investigated by linear stability analysis. The stability boundaries for the axisymmetric flow are derived for Prandtl numbers from 0.02 to 1, for aspect ratio $A$ ($A= H/ R= \mathrm{height} / \mathrm{radius} $) equal to 1, 0.9, 0.8, 0.7, respectively. We found that there still exists stable non-trivial axisymmetric flow beyond the second bifurcation in certain ranges of Prandtl number for $A= 1$, $0. 9$ and 0.8, excluding the $A= 0. 7$ case. The finding for $A= 0. 7$ is that very frequent changes of critical mode (azimuthal Fourier mode) of the second bifurcation occur when the Prandtl number is changed, where five kinds of steady modes $m= 1, 2, 8, 9, 10$ and three kinds of oscillatory modes $m= 3, 4, 6$ are presented. These multiple modes indicate different flow structures triggered at the transitions. The instability mechanism of the flow is explained by kinetic energy transfer analysis, which shows that the radial or axial shear of base flow combined with buoyancy mechanism leads to the instability results.


1999 ◽  
Vol 388 ◽  
pp. 171-195 ◽  
Author(s):  
CHRISTIAN PANKIEWITZ ◽  
ECKART MEIBURG

The influence of a non-monotonic viscosity–concentration relationship on miscible displacements in porous media is studied for radial source flows and the quarter five-spot configuration. Based on linear stability results, a parametric study is presented that demonstrates the dependence of the dispersion relations on both the Péclet number and the parameters of the viscosity profile. The stability analysis suggests that any displacement can become unstable provided only that the Péclet number is sufficiently high. In contrast to rectilinear flows, for a given end-point viscosity ratio an increase of the maximum viscosity generally has a destabilizing effect on the flow. The physical mechanisms behind this behaviour are examined by inspecting the eigensolutions to the linear stability problem. Nonlinear simulations of quarter five-spot displacements, which for small times correspond to radial source flows, confirm the linear stability results. Surprisingly, displacements characterized by the largest instability growth rates, and consequently by vigorous viscous fingering, lead to the highest breakthrough recoveries, which can even exceed that of a unit mobility ratio flow. It can be concluded that, for non-monotonic viscosity profiles, the interaction of viscous fingers with the base-flow vorticity can result in improved recovery rates.


1988 ◽  
Vol 43 (12) ◽  
pp. 1009-1016 ◽  
Author(s):  
D. Lortz

The stability of axisymmetric ideal MHD equilibria which are symmetric with respect to the equatorial plane is considered. It is found that for external axisymmetric modes which are antisymmetric with respect to the equatorial plane and for profiles such that the current density vanishes at the free plasma boundary the stability problem reduces to a classical interior-exterior scalar eigenvalue problem. Because of the separation property the resulting stability condition is necessary and sufficient and is thus more stringent than criteria derived by choosing special test functions, e.g. the vertical shift condition.


Sign in / Sign up

Export Citation Format

Share Document