Numerical Computation of the Stability of Plane Poiseuille Flow

1978 ◽  
Vol 45 (1) ◽  
pp. 13-18 ◽  
Author(s):  
L. Wolf ◽  
Z. Lavan ◽  
H. J. Nielsen

The hydrodynamic stability of plane Poiseuille flow to infinitesimal and finite amplitude disturbances is investigated using a direct numerical technique. The governing equations are cast in terms of vorticity and stream function using second-order central differences in space. The vorticity equation is used to advance the vorticity values in time and successive over-relaxation is used to solve the stream function equation. Two programs were prepared, one for the linearized and the other for the complete disturbance equations. Results obtained by solving the linearized equations agree well with existing solutions for small disturbances. The nonlinear calculations reveal that the behavior of a disturbance depends on the amplitude and on the wave number. The behavior at wave numbers below and above the linear critical wave number is drastically different.

1972 ◽  
Vol 51 (2) ◽  
pp. 403-416 ◽  
Author(s):  
C. H. Green ◽  
C. H. Ellen

This paper examines the linear stability of antisymmetric disturbances in incompressible plane Poiseuille flow between identical flexible walls which undergo transverse displacements. Using a variational approach, an approximate solution of the problem is formulated in a form suitable for computational evaluation of the (complex) wave speeds of the system. A feature of this formulation is that the varying boundary conditions (and the Orr-Sommerfeld equation) are satisfied only in the mean; this reduces the labour involved in determining the approximate solution for a variety of wall conditions without increasing the difficulty of obtaining solutions to a given accuracy. In this paper the symmetric stream function distribution across the channel is represented by a series of cosines whose coefficients are determined by the variational solution. Comparisons with previous work, both for the flexible-wall and rigid-wall problems, show that the method gives results as accurate as those obtained previously by other methods while new results, for flexible walls, indicate the presence of a higher wave-number stability boundary which joins the distorted Tollmien-Schlichting stability boundary at lower wave-numbers. In some cases this upper unstable region, which is characterized by large amplification rates, may determine the critical Reynolds number of the system.


The stability of plane Poiseuille flow in a channel forced by a wavelike motion on one of the channel walls is investigated. The amplitude Є of this forcing is taken to be small. The most dangerous modes of forcing are identified and it is found in general the critical Reynolds number is changed by O (Є) 2 . However, we identify two particular modes of forcing which give rise to decrements of order Є 2/3 and Є in the critical Reynolds number. Some types of forcing are found to generate sub critical stable finite amplitude perturbations to plane Poiseuille flow. This contrasts with the unforced case where the only stable solution is the zero amplitude solution. The forcing also deforms the unstable subcritical limit cycle solution from its usual circular shape into a more complicated shape. This has an effect on the threshold amplitude ideas suggested by, for example, Meksyn & Stuart (1951). It is found that the phase of disturbances must also be considered when finding the amplitude dependent critical Reynolds numbers.


1967 ◽  
Vol 29 (1) ◽  
pp. 31-38 ◽  
Author(s):  
Chaim L. Pekeris ◽  
Boris Shkoller

Stuart (1960) has developed a theory of the stability of plane Poiseuille flow to periodic disturbances of finite amplitude which, in the neighbourhood of the neutral curve, leads to an equation of the Landau (1944) type for the amplitude A of the disturbance: \[ d|A|^2/dt = k_1|A|^2 - k_2|A|^4. \] If k2 is positive in the supercritical region (R > RC) where k1 is positive, then, according to Stuart, there is a possibility of the existence of periodic solutions of finite amplitude which asymptotically approach a constant value of (k1/k2)½. We have evaluated the coefficient k2 and found that there indeed exists a zone in the (α, R)-plane where it is positive. This is the zone inside the dashed curve shown in figure 1, with the region of instability predicted by the linear theory included inside the ‘neutral curve’. Stuart's theory and Eckhaus's generalization thereof could apply in the overlapping zone just above the lower branch of the neutral curve.


1999 ◽  
Vol 67 (4) ◽  
pp. 834-837 ◽  
Author(s):  
R. E. Khayat ◽  
N. Ashrafi

The purely elastic stability and bifurcation of the one-dimensional plane Poiseuille flow is determined for a large class of Oldroyd fluids with added viscosity, which typically represent polymer solutions composed of a Newtonian solvent and a polymeric solute. The problem is reduced to a nonlinear dynamical system using the Galerkin projection method. It is shown that elastic normal stress effects can be solely responsible for the destabilization of the base (Poiseuille) flow. It is found that the stability and bifurcation picture is dramatically influenced by the solvent-to-solute viscosity ratio, ε. As the flow deviates from the Newtonian limit and ε decreases below a critical value, the base flow loses its stability. Two static bifurcations emerge at two critical Weissenberg numbers, forming a closed diagram that widens as the level of elasticity increases. [S0021-8936(00)00703-0]


1983 ◽  
Vol 133 ◽  
pp. 179-206 ◽  
Author(s):  
P. K. Sen ◽  
D. Venkateswarlu

In this work a study has been made of the Stuart (1960)–Watson (1960) formalism as applied to plane Poiseuille flow. In particular, the higher-order Landau coefficients have been calculated for the Reynolds & Potter (1967) method and for the Watson (1960) method. The results have been used to study the convergence of the Stuart–Landau series. A convergence curve in the (α, R)-plane has been obtained by using suitable Domb–Sykes plots. In the region of poor convergence of the series, and also in a part of the divergent region of the series, it has been found that the Shanks (1955) method, using the em1 transformation, serves as a very effective way of finding the proper sum of the series, or of finding the proper antilimit of the series. The results for the velocity calculations at R = 5000 are in very good agreement with Herbert's (1977) Fourier-truncation method using N = 4. The Watson method and the Reynolds & Potter method have also been compared inthe subcritical and supercritical regions. It is found in the supercritical region that there is not much difference in the results by the ‘true problem’ of Watson and the ‘false problem’ of Reynolds & Potter when the respective series in both methods are summed by the Shanks method. This fact could possibly be capitalized upon in the subcritical region, where the Watson method is difficult to apply.


Sign in / Sign up

Export Citation Format

Share Document