scholarly journals Linear stability analysis and direct numerical simulation of two-layer channel flow

2016 ◽  
Vol 798 ◽  
pp. 889-909 ◽  
Author(s):  
Kirti Chandra Sahu ◽  
Rama Govindarajan

We study the stability of two-fluid flow through a plane channel at Reynolds numbers of 100–1000 in the linear and nonlinear regimes. The two fluids have the same density but different viscosities. The fluids, when miscible, are separated from each other by a mixed layer of small but finite thickness, across which the viscosity changes from that of one fluid to that of the other. When immiscible, the interface is sharp. Our study spans a range of Schmidt numbers, viscosity ratios and locations and thicknesses of the mixed layer. A region of instability distinct from that of the Tollmien–Schlichting mode is obtained at moderate Reynolds numbers. We show that the overlap of the layer of viscosity-stratification with the critical layer of the dominant disturbance provides a mechanism for this instability. At very low values of diffusivity, the miscible flow behaves exactly like the immiscible one in terms of stability characteristics. High levels of miscibility make the flow more stable. At intermediate levels of diffusivity however, in both linear and nonlinear regimes, miscible flow can be more unstable than the corresponding immiscible flow without surface tension. This difference is greater when the thickness of the mixed layer is decreased, since the thinner the layer of viscosity stratification, the more unstable the miscible flow. In direct numerical simulations, disturbance growth occurs at much earlier times in the miscible flow, and also the miscible flow breaks spanwise symmetry more readily to go into three-dimensionality. The following observations hold for both miscible and immiscible flows without surface tension. The stability of the flow is moderately sensitive to the location of the interface between the two fluids. The response is non-monotonic, with the least stable location of the layer being mid-way between the wall and the centreline. As expected, flow at higher Reynolds numbers is more unstable.

1978 ◽  
Vol 19 (2) ◽  
pp. 309-318 ◽  
Author(s):  
K. Aderogba ◽  
J.R. Blake

In the paper by Aderogba and Blake [1], the first order approximation to the shape of the interface between two immiscible liquids, at which surface tension acted, was obtained in the case of a point force directed parallel to the interface. In the case of the normal force no solution was obtained due to the logarithmic singularity associated with problems of this type. However, the problem is not physically well-posed in the case of the normal force. Surface tension is not suitable because it cannot alone balance the induced stress on the interface due to a point force. We need an additional force to balance the action of the point force on the interface. The obvious solution is to include a density difference Δρ* between the two fluidswhere and are the densities of the lower and upper fluid respectively.


2011 ◽  
Vol 683 ◽  
pp. 112-148 ◽  
Author(s):  
G. Pfingstag ◽  
B. Audoly ◽  
A. Boudaoud

AbstractWe study the stability of a thin, Newtonian viscous sheet floating on a bath of denser fluid. We first derive a general set of equations governing the evolution of a nearly flat sheet, accounting for geometrical nonlinearities associated with moderate rotations. We extend two classical models by considering arbitrary external body and surface forces; these two models follow from different scaling assumptions, and are derived in a unified way. The equations capture two modes of deformation, namely viscous bending and stretching, and describe the evolution of thickness, mid-surface and in-plane velocity as functions of two-dimensional coordinates. These general equations are applied to a floating viscous sheet, considering gravity, buoyancy and surface tension. We investigate the stability of the flat configuration when subjected to arbitrary in-plane strain. Two unstable modes can be found in the presence of compression. The first one combines undulations of the centre-surface and modulations of the thickness, with a wavevector perpendicular to the direction of maximum applied compression. The second one is a buckling mode; it is purely undulatory and has a wavevector along the direction of maximum compression. A nonlinear analysis yields the long-time evolution of the undulatory mode.


2005 ◽  
Vol 46 (4) ◽  
pp. 507-543 ◽  
Author(s):  
C. E. Watson ◽  
S. R. Otto

AbstractWe consider the stability of high Reynolds number flow past a heated, curved wall. The influence of both buoyancy and curvature, with the appropriate sense, can render a flow unstable to longitudinal vortices. However, conversely each mechanism can make a flow more stable; as with a stable stratification or a convex curvature. This is partially due to their influence on the basic flow and also due to additional terms in the stability equations. In fact the presence of buoyancy in combination with an appropriate local wall gradient can actually increase the wall shear and these effects can lead to supervelocities and the promotion of a wall jet. This leads to the interesting discovery that the flow can be unstable for both concave and convex curvatures. Furthermore, it is possible to observe sustained vortex growth in stably stratified boundary layers over convexly curved walls. The evolution of the modes is considered in both the linear and nonlinear régimes.


2021 ◽  
Vol 5 (3) ◽  
pp. 37
Author(s):  
Hernán Martinelli ◽  
Claudia Domínguez ◽  
Marcos Fernández Leyes ◽  
Sergio Moya ◽  
Hernán Ritacco

In the search for responsive complexes with potential applications in the formulation of smart dispersed systems such as foams, we hypothesized that a pH-responsive system could be formulated with polyacrylic acid (PAA) mixed with a cationic surfactant, Gemini 12-2-12 (G12). We studied PAA-G12 complexes at liquid–air interfaces by equilibrium and dynamic surface tension, surface rheology, and X-ray reflectometry (XRR). We found that complexes adsorb at the interfaces synergistically, lowering the equilibrium surface tension at surfactant concentrations well below the critical micelle concentration (cmc) of the surfactant. We studied the stability of foams formulated with the complexes as a function of pH. The foams respond reversibly to pH changes: at pH 3.5, they are very stable; at pH > 6, the complexes do not form foams at all. The data presented here demonstrate that foam formation and its pH responsiveness are due to interfacial dynamics.


Author(s):  
P. G. Drazin

ABSTRACTSome aspects of generation of water waves by wind and of turbulence in a heterogeneous fluid may be described by the theory of hydrodynamic stability. The technical difficulties of these problems of instability have led to obscurities in the literature, some of which are elucidated in this paper. The stability equation for a basic steady parallel horizontal flow under the influence of gravity is derived carefully, the undisturbed fluid having vertical variations of density and viscosity. Methods of solution of the equation for large Reynolds numbers and for long-wave disturbances are described. These methods are applied to simple models of wind blowing over water and of fresh water flowing over salt water.


2002 ◽  
Vol 452 ◽  
pp. 163-187 ◽  
Author(s):  
C. L. BURCHAM ◽  
D. A. SAVILLE

A liquid bridge is a column of liquid, pinned at each end. Here we analyse the stability of a bridge pinned between planar electrodes held at different potentials and surrounded by a non-conducting, dielectric gas. In the absence of electric fields, surface tension destabilizes bridges with aspect ratios (length/diameter) greater than π. Here we describe how electrical forces counteract surface tension, using a linearized model. When the liquid is treated as an Ohmic conductor, the specific conductivity level is irrelevant and only the dielectric properties of the bridge and the surrounding gas are involved. Fourier series and a biharmonic, biorthogonal set of Papkovich–Fadle functions are used to formulate an eigenvalue problem. Numerical solutions disclose that the most unstable axisymmetric deformation is antisymmetric with respect to the bridge’s midplane. It is shown that whilst a bridge whose length exceeds its circumference may be unstable, a sufficiently strong axial field provides stability if the dielectric constant of the bridge exceeds that of the surrounding fluid. Conversely, a field destabilizes a bridge whose dielectric constant is lower than that of its surroundings, even when its aspect ratio is less than π. Bridge behaviour is sensitive to the presence of conduction along the surface and much higher fields are required for stability when surface transport is present. The theoretical results are compared with experimental work (Burcham & Saville 2000) that demonstrated how a field stabilizes an otherwise unstable configuration. According to the experiments, the bridge undergoes two asymmetric transitions (cylinder-to-amphora and pinch-off) as the field is reduced. Agreement between theory and experiment for the field strength at the pinch-off transition is excellent, but less so for the change from cylinder to amphora. Using surface conductivity as an adjustable parameter brings theory and experiment into agreement.


2005 ◽  
Vol 19 (28n29) ◽  
pp. 1547-1550
Author(s):  
YOULIANG CHENG ◽  
XIN LI ◽  
ZHONGYAO FAN ◽  
BOFEN YING

Representing surface tension by nonlinear relationship on temperature, the boundary value problem of linear stability differential equation on small perturbation is derived. Under the condition of the isothermal wall the effects of nonlinear surface tension on stability of heat transfer in saturated liquid film of different liquid low boiling point gases are investigated as wall temperature is varied.


2006 ◽  
Vol 17 (01) ◽  
pp. 65-73 ◽  
Author(s):  
SHIRO SAWADA

The optimal velocity model which depends not only on the headway but also on the relative velocity is analyzed in detail. We investigate the effect of considering the relative velocity based on the linear and nonlinear analysis of the model. The linear stability analysis shows that the improvement in the stability of the traffic flow is obtained by taking into account the relative velocity. From the nonlinear analysis, the relative velocity dependence of the propagating kink solution for traffic jam is obtained. The relation between the headway and the velocity and the fundamental diagram are examined by numerical simulation. We find that the results by the linear and nonlinear analysis of the model are in good agreement with the numerical results.


1998 ◽  
Vol 358 ◽  
pp. 357-378 ◽  
Author(s):  
M. NAGATA

The stability of nonlinear tertiary solutions in rotating plane Couette flow is examined numerically. It is found that the tertiary flows, which bifurcate from two-dimensional streamwise vortex flows, are stable within a certain range of the rotation rate when the Reynolds number is relatively small. The stability boundary is determined by perturbations which are subharmonic in the streamwise direction. As the Reynolds number is increased, the rotation range for the stable tertiary motions is destroyed gradually by oscillatory instabilities. We expect that the tertiary flow is overtaken by time-dependent motions for large Reynolds numbers. The results are compared with the recent experimental observation by Tillmark & Alfredsson (1996).


1994 ◽  
Vol 258 ◽  
pp. 131-165 ◽  
Author(s):  
Peter W. Duck ◽  
Gordon Erlebacher ◽  
M. Yousuff Hussaini

The linear stability of compressible plane Couette flow is investigated. The appropriate basic velocity and temperature distributions are perturbed by a small-amplitude normal-mode disturbance. The full small-amplitude disturbance equations are solved numerically at finite Reynolds numbers, and the inviscid limit of these equations is then investigated in some detail. It is found that instabilities can occur, although the corresponding growth rates are often quite small; the stability characteristics of the flow are quite different from unbounded flows. The effects of viscosity are also calculated, asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional regimes to the problem occur when the wave speed of the disturbances approaches the velocity of either of the walls, and these regimes are also analysed in some detail. Finally, the effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place of impermeability) is investigated, and shown to yield results common to both bounded and unbounded flows.


Sign in / Sign up

Export Citation Format

Share Document