Stability of Thermosolutal Natural Convection in Superposed Fluid and Porous Layers

This article deals with the onset of thermosolutal natural convection in horizontal superposed fluid and porous layers. A linear stability analysis is performed using the one-domain approach. As in the thermal convection case, the results show a bimodal nature of the marginal stability curves where each mode corresponds to a different convective instability. At small wave numbers, the convective flow occurs in the whole cavity (“porous mode”) while perturbations of large wave numbers lead to a convective flow mainly confined in the fluid layer (“fluid mode”). Furthermore, it is shown that the onset of thermosolutal natural convection is characterized by a multi-cellular flow in the fluid region for negative thermal Rayleigh numbers. For positive thermal Rayleigh numbers, the convective flow takes place both in the fluid and porous regions. The influence of the depth ratio and thermal diffusivity ratio is also investigated for a wide range of the thermal Rayleigh numbers.     

Absolute and convective instability of a liquid sheet with transverse temperature gradient

The spatial–temporal instability behavior of a viscous liquid sheet with temperature difference between the two surfaces was investigated theoretically. The practical situation motivating this investigation is liquid sheet heated by ambient gas, usually encountered in industrial heat transfer and liquid propellant rocket engines. The existing dispersion relation was used, to explore the spatial–temporal instability of viscous liquid sheets with a nonuniform temperature profile, by setting both the wave number and frequency complex. A parametric study was performed in both sinuous and varicose modes to test the influence of dimensionless numbers on the transition between absolute and convective instability of the flow. For a small value of liquid Weber number, or a great value of gas-to-liquid density ratio, the flow was found to be absolutely unstable. The absolute instability was enhanced by increasing the liquid viscosity. It was found that variation of the Marangoni number hardly influenced the absolute instability of the sinuous mode of oscillations; however it slightly affected the absolute instability in the varicose mode.     

On the absolute instability in a boundary-layer flow with compliant coatings

The question of absolute instabilities occuring in a boundary-layer flow with compliant coatings is reassessed. Compliant coatings of the Kramer's type are considered. Performing a local, linear absolute/convective stability analysis, a family of spring-backed elastic plates with damping is shown to be absolutely unstable for sufficiently thin plates. The absolute instability arises from the coalescence between an upstream propagating evanescent mode and the Tollmien–Schlichting wave. To reinforce the local, linear stability results the global stability behaviour of the system is investigated, integrating numerically the full nonparallel and nonlinear two-dimensional Navier–Stokes system coupled to the dynamical model. Injecting Gaussian-type, spatially localized flow disturbances as initial conditions, the spatio-temporal evolution of wave packets is computed. The absolute stability behaviour is retrieved in the global system, for a compliant panel of finite length. It is demonstrated numerically that the global stability behaviour of the wall, triggered by finite-end-effects, may be independent of the disturbance propagation in the flow.     

The Absolute Instability of Thin Wakes in an Incompressible/Compressible Fluid

RID="ID=" Communicated by P. HallAbstract:The absolute/convective instability of two-dimensional wakes forming behind a flat plate and near the trailing-edge of a thin wedge-shaped aerofoil in an incompressible/compressible fluid is investigated. The mean velocity profiles are obtained by solving numerically the classical compressible boundary-layer equations with a negative pressure gradient for the flat plate case, and the incompressible triple-deck equations for a thin wedge-shaped trailing-edge. In addition for a Joukowski aerofoil the incompressible mean boundary-layer flow in the vicinity of the trailing-edge is also calculated by solving the interactive boundary-layer equations. A linear stability analysis of the boundary-layer profiles shows that a pocket of absolute instability occurs downstream of the trailing-edge with the extent of the instability region increasing with more adverse pressure gradients. The region of absolute instability persists along the near-wake axis, while the majority of the wake is convectively unstable. For a thin wedge-shaped trailing-edge in an incompressible fluid, a similar stability analysis of the velocity profiles obtained via a composite expansion, also shows the occurrence of absolute instability behind the trailing-edge for a wedge angle greater than a critical value. For increasing values of the wedge angle and for thicker aerofoils, separation takes place near the trailing-edge and the extent of absolute instability increases. Calculations also show that for insulated plates compressibility has a stabilizing effect but cooling the wall destabilizes the flow unlike wall heating.} Received 11 May 1998 and accepted 25 February 1999     

Convective and absolute instabilities in non-Boussinesq mixed convection

The problem of non-Boussinesq mixed convection in a vertical channel formed by two differentially heated infinite plates is investigated and the complete convective/absolute instability boundary is computed for a wide range of physical parameters. A physical insight into the mechanisms causing instabilities is given. In particular, it is shown that the appearance of absolute instability is always dictated by a flow reversal within a channel; however, existence of the flow reversal does not exclude the possibility of convective instability. It is also shown that fluid’s non-linear transport property variations have a dramatic effect on the structure and complexity of spatio-temporal instabilities of the co-existing buoyancy and shear modes as the temperature difference across the channel increases. The validity of the stability results obtained using the procedure described in Suslov (J Comp Phys 212, 188–217, 2006) is assessed using the method of steepest descent. This work was partially supported by a computing grant from the Australian Partnership for Advanced Computing, 2000–2003.     

粘弹性地基上粘弹性输流管道的稳定性分析

从Winkler假设和单轴线性粘弹性本构方程出发,推导了Kelvin-Voigt粘弹性地基上三参量固体模型输流管道的运动微分方程,采用改进的有限差分法,分析了管道和地基的粘弹性参数对输流管道无量纲复频率和无量纲流速之间的变化关系的影响。     

Convective Instability of Oldroyd-B Fluid Saturated Porous Layer Heated from Below using a Thermal Non-equilibrium Model

The linear stability of a viscoelastic fluid saturated densely packed horizontal porous layer heated from below and cooled from above is investigated by considering the Oldroyd-B type fluid. A generalized Darcy model, which takes into account the viscoelastic properties, is employed as momentum equation and a two-field model is used for energy equation each representing solid and fluid phases separately. Linear stability analysis suggests that, there is a competition between the processes of viscous relaxation and thermal diffusion that causes the first convective instability to be oscillatory rather than stationary. Analytical expression for the occurrence of oscillatory onset is obtained, and it is found that the necessary condition for the existence of the same is Λ < 1. Besides, the effect of viscoelastic parameters and the thermal non-equilibrium on the stability of the system is analyzed.     

Conditional and unconditional nonlinear stability for convection induced by absorption of radiation in a porous medium

Convection induced by the selective absorption of radiation is investigated, for the case of an internal heat source that is modelled quadratically with respect to concentration. The growth rate for the linearised system is shown to be real, and a linear instability analysis is performed. To establish conditional and unconditional nonlinear stability results, both the Darcy and Forchheimer models are employed to describe fluid flow. Due to the presence of significant regions of potential subcritical instabilities, the results indicate that linear theory may only be accurate enough to predict the onset of convective motion when the model for the internal heat source is predominantly linear.Received: 6 May 2003, Accepted: 9 August 2003, Published online: 12 December 2003     

Stability of triple diffusive convection in a viscoelastic fluid-saturated porous layer

The triple diffusive convection in an Oldroyd-B fluid-saturated porous layer is investigated by performing linear and weakly nonlinear stability analyses. The condition for the onset of stationary and oscillatory is derived analytically. Contrary to the observed phenomenon in Newtonian fluids, the presence of viscoelasticity of the fluid is to degenerate the quasiperiodic bifurcation from the steady quiescent state. Under certain conditions, it is found that disconnected closed convex oscillatory neutral curves occur, indicating the requirement of three critical values of the thermal Darcy-Rayleigh number to identify the linear instability criteria instead of the usual single value, which is a novel result enunciated from the present study for an Oldroyd-B fluid saturating a porous medium. The similarities and differences of linear instability characteristics of Oldroyd-B, Maxwell, and Newtonian fluids are also highlighted. The stability of oscillatory finite amplitude convection is discussed by deriving a cubic Landau equation, and the convective heat and mass transfer are analyzed for different values of physical parameters.     

Global stability behaviour for the BEK family of rotating boundary layers

Numerical simulations were conducted to investigate the linear global stability behaviour of the Bödewadt, Ekman, von Kármán (BEK) family of flows, for cases where a disc rotates beneath an incompressible fluid that is also rotating. This extends the work reported in recent studies that only considered the rotating-disc boundary layer with a von Kármán configuration, where the fluid that lies above the boundary layer remains stationary. When a homogeneous flow approximation is made, neglecting the radial variation of the basic state, it can be shown that linearised disturbances are susceptible to absolute instability. We shall demonstrate that, despite this prediction of absolute instability, the disturbance development exhibits globally stable behaviour in the BEK boundary layers with a genuine radial inhomogeneity. For configurations where the disc rotation rate is greater than that of the overlying fluid, disturbances propagate radially outwards and there is only a convective form of instability. This replicates the behaviour that had previously been documented when the fluid did not rotate beyond the boundary layer. However, if the fluid rotation rate is taken to exceed that of the disc, then the propagation direction reverses and disturbances grow while convecting radially inwards. Eventually, as they approach regions of smaller radii, where stability is predicted according to the homogeneous flow approximation, the growth rates reduce until decay takes over. Given sufficient time, such disturbances can begin to diminish at every radial location, even those which are positioned outwards from the radius associated with the onset of absolute instability. This leads to the confinement of the disturbance development within a finitely bounded region of the spatial–temporal plane.     

Nonlinear development of two-dimensional hydroelastic instability in a turbulent boundary layer on an elastic coating

Nonlinear evolution of hydroelastic instability arising in the flow past a coating of a rubber-type material by a turbulent boundary layer of an incompressible fluid is studied. A nonlinear dispersion equation for two-dimensional, quasi-monochromatic, low-amplitude waves is derived. The Prandtl equations for the mean (over the waviness period) boundary-layer flow are solved in the approximation of local similarity and by direct numerical integration. Evolution of unstable waves in time is studied on the basis of the Landau equation, which is derived separately for the instability of fast waves (flutter) and the quasi-static instability (divergence). The calculation results are compared with available experimental data. Institute of Applied Physics, Russian Academy of Sciences, Nizhnii Novgorod 603600. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 4, pp. 69–80, July–August, 2000.     

On the stability of the shear flow of a viscoelastic fluid with slip along the fixed wall

In this paper we solve the time-dependent shear flow of an Oldroyd-B fluid with slip along the fixed wall. We use a non-linear slip model relating the shear stress to the velocity at the wall and exhibiting a maximum and a minimum. We assume that the material parameters in the slip equation are such that multiple steady-state solutions do not exist. The stability of the steady-state solutions is investigated by means of a one-dimensional linear stability analysis and by numerical calculations. The instability regimes are always within or coincide with the negative-slope regime of the slip equation. As expected, the numerical results show that the instability regimes are much broader than those predicted by the linear stability analysis. Under our assumptions for the slip equation, the Newtonian solutions are stable everywhere. The interval of instability grows as one moves from the Newtonian to the upper-convected Maxwell model. Perturbing an unstable steady-state solution leads to periodic solutions. The amplitude and the period of the oscillations increase with elasticity.     

The pendulum-slosh problem: Simulation using a time-dependent conformal mapping

Suspending a rectangular vessel which is partially filled with fluid from a single rigid pivoting pole produces an interesting theoretical model with which to investigate the dynamic coupling between fluid motion and vessel rotation. The exact equations for this coupled system are derived with the fluid motion governed by the Euler equations relative to the moving frame of the vessel, and the vessel motion governed by a modified forced pendulum equation. The nonlinear equations of motion for the fluid are solved numerically via a time-dependent conformal mapping, which maps the physical domain to a rectangle in the computational domain with a time dependent conformal modulus. The numerical scheme expresses the implicit free-surface boundary conditions as two explicit partial differential equations which are then solved via a pseudo-spectral method in space. The coupled system is integrated in time with a fourth-order Runge–Kutta method. The starting point for the simulations is the linear neutral stability contour discovered by Turner et al. (2015, Journal of Fluid & Structures 52, 166–180). Near the contour the nonlinear results confirm the instability boundary, and far from the neutral curve (parameterized by longer pole lengths) nonlinearity is found to significantly alter the vessel response. Results are also presented for an initial condition given by a superposition of two sloshing modes with approximately the same frequency from the linear characteristic equation. In this case the fluid initial conditions generate large nonlinear vessel motions, which may have implications for systems designed to oscillate in a confined space or on the slosh-induced-rolling of a ship.     

Dynamics of pipes conveying fluid with non-uniform turbulent and laminar velocity profiles

Previous analytical work on stability of fluid-conveying pipes assumed a uniform velocity profile for the conveyed fluid. In real fluid flows, the presence of viscosity leads to a sheared region near the wall. Earlier studies correctly note that viscous forces do not affect the dynamics of the system since these forces are balanced by pressure drop in the conveyed fluid. Although viscous shear has not been ignored in these studies, a uniform velocity profile assumes that the sheared region is infinitely thin. Prior analysis was extended to account for a fully developed non-uniform profile such as would be encountered in real fluid flows. A modified, highly tractable equation of motion was derived, which includes a single additional parameter to account for the true momentum of the fluid. This empirical parameter was determined by numerical analysis over the Reynolds number range of interest. The stability of cantilever pipes conveying fluid with two types of non-uniform velocity profile was assessed. In the first case, the profile was a function of Reynolds number and transition to turbulence occurred before the onset of flutter instability. This case had stability properties similar to the uniform velocity case except in specific narrow regions of the parameter space. The second case required that the Reynolds number be such that the flow was always laminar. For this case, lower fluid velocity was required to achieve instability, and the oscillation frequency at instability was considerably lower over much of the parameter space, compared to the uniform case.     

Stability of convective motion in a vertical layer in the presence of longitudinal vibrations

The influence of vibrations of a cavity containing a fluid on the convective stability of the equilibrium has been investigated on a number of occasions [1]. The stability of convective flows in a modulated gravity field has not hitherto been studied systematically. There is only the paper of Baxi, Arpaci, and Vest [2], which contains fragmentary data corresponding to various values of the determining parameters of the problem. The present paper investigates the linear stability of convective flow in a vertical plane layer with walls at different temperatures in the presence of longitudinal harmonic vibrations of the cavity containing the fluid. It is assumed that the frequency of the vibrations is fairly high; the motion is described by the equations of the averaged convective motion. The stability boundaries of the flow with respect to monotonic perturbations in the region of Prandtl numbers 0 ? P ? 10 are determined. It is found that high-frequency vibrations have a destabilizing influence on the convective motion. At sufficiently large values of the vibration parameter, the flow becomes unstable at arbitrarily small values of the Grashof number, this being due to the mechanism of vibrational convection, which leads to instability even under conditions of weightlessness, when the main flow is absent [3, 4].     

The flow-induced instability of long hanging pipes

The effect of increasing length on the stability of a hanging fluid-conveying pipe is investigated. Experiments show that there exists a critical length above which the flow velocity necessary to cause flutter becomes independent of the pipe length. The fluid-structure interaction is thus modelled by following the work of Bourrières and of Paı&#x0308;doussis. Computations using a standard Galerkin method confirm this evolution. A short pipe model is then considered, where gravity plays a negligible role. Transition between this short length model and the asymptotic situation is found to occur where a local stability criterion is satisfied at the upstream end of the pipe. For longer pipes, a model is proposed where the zone of stable waves is totally disregarded. Comparison of these models with experiments and computations show a good agreement over all ranges of mass ratios between the flowing fluid and the pipe.     

Magnetic effect on instability and nonlinear stability in a reacting fluid

We study the problem of convective movement of a reacting solute in a viscous incompressible fluid occupying a plane layer and subjected to a vertical magnetic field. The thresholds for linear instability are found and compared to those derived by a global nonlinear energy stability analysis.     

Convective instability in a two-layer system of reacting fluids with concentration-dependent diffusion

The development of convective instability in a two-layer system of miscible fluids placed in a narrow vertical gap has been studied theoretically and experimentally. The upper and lower layers are formed with aqueous solutions of acid and base, respectively. When the layers are brought into contact, the frontal neutralization reaction begins. We have found experimentally a new type of convective instability, which is characterized by the spatial localization and the periodicity of the structure observed for the first time in the miscible systems. We have tested a number of different acid–base systems and have found a similar patterning there. In our opinion, it may indicate that the discovered effect is of a general nature and should be taken into account in reaction–diffusion–convection problems as another tool with which the reaction can govern the movement of the reacting fluids. We have shown that, at least in one case (aqueous solutions of nitric acid and sodium hydroxide), a new type of instability called as the concentration-dependent diffusion convection is responsible for the onset of the fluid flow. It arises when the diffusion coefficients of species are different and depend on their concentrations. This type of instability can be attributed to a variety of double-diffusion convection. A mathematical model of the new phenomenon has been developed using the system of reaction–diffusion–convection equations written in the Hele–Shaw approximation. It is shown that the instability can be reproduced in the numerical experiment if only one takes into account the concentration dependence of the diffusion coefficients of the reagents. The dynamics of the base state, its linear stability and nonlinear development of the instability are presented. It is also shown that by varying the concentration of acid in the upper layer one can achieve the occurrence of chemo-convective solitary cell in the bulk of an almost immobile fluid. Good agreement between the experimental data and the results of numerical simulations is observed.     

A minimal model for flow control on an aerofoil using a poro-elastic coating

A minimal model is obtained for vortex-shedding from an aerofoil with a porous coating of flow-compliant feather-like actuators, in order to better understand this passive way to achieve flow control. This model is realized by linearly coupling a minimal-order model for vortex-shedding from the same aerofoil without any such coating with an equation for the poro-elastic coating, here modelled as a linear damped oscillator. The various coefficients in this model, derived using perturbation techniques, aid in our understanding of the physics of this fluid–structure interaction problem. The minimal model for a coated aerofoil indicates the presence of distinct regimes that are dependent on the flow and coating characteristics. The models and the parametric studies performed provide insight into the selection of optimal coating parameters, to enable flow control at low Reynolds numbers.