Influence of swirl on the stability of a rod in annular leakage flow

The influence of swirl (flow rotation) on the stability of a rod in annular leakage flow is investigated. Under the assumption of laminar flow and plane vibrations (no whirling), it is shown that the swirl acts, in effect, as an elastic foundation with negative foundation stiffness, the magnitude being proportional to the mean circumferential flow rate squared. Consequently, swirl always lowers the critical axial flow speed in case of divergence instability of a rod of finite length. Numerical analysis is needed to predict the effect of swirl in case of flutter instability of a finite rod; this is not performed here. However, for the flutter-like instability of travelling waves in an infinite rod-channel system, it is shown analytically that swirl again always lowers the critical axial flow speed. Finally, it is found that by circumferential flow alone, the travelling waves are extinguished at a certain flow rate, followed by a divergence-like instability.     

流体黏性对输流管道运动方程及临界流速的影响

考虑实际流体黏性引起的管内流速非均匀分布,针对层流和两种不同的湍流流态,对理想流体情况下输流管道运动方程中的离心力项进行了修正,得到的修正系数分别为1.333(圆管层流)、1.020(光滑管壁圆管湍流)和1.037～1.055(粗糙管壁圆管湍流).根据修正后的运动方程得到的上述3种情况下的发散失稳临界流速比理想流体流动情况下依次分别低13.4%,1.0%和1.8%～2.6%.流体黏性对输流管道运动方程及临界流速的影响只与流态有关,雷诺数决定流态,而黏性系数通过雷诺数间接起作用.     

Interaction between a cantilevered-free flexible plate and ideal flow

We develop a new computational model of the linear fluid–structure interaction of a cantilevered flexible plate with an ideal flow in a channel. The system equation is solved via numerical simulations that capture transients and allow the spatial variation of the flow–structure interaction on the plate to be studied in detail. Alternatively, but neglecting wake effects, we are able to extract directly the system eigenvalues to make global predictions of the system behaviour in the infinite-time limit. We use these complementary approaches to conduct a detailed study of the fluid–structure system. When the channel walls are effectively absent, predictions of the critical velocity show good agreement with those of other published work. We elucidate the single-mode flutter mechanism that dominates the response of short plates and show that the principal region of irreversible energy transfer from fluid to structure occurs over the middle portion of the plate. A different mechanism, modal-coalescence flutter, is shown to cause the destabilisation of long plates with its energy transfer occurring closer to the trailing edge of the plate. This mechanism is shown to allow a continuous change to higher-order modes of instability as the plate length is increased. We then show how the system response is modified by the inclusion of channel walls placed symmetrically above and below the flexible plate, the effect of unsteady vorticity shed at the trailing edge of the plate, and the effect of a rigid surface placed upstream of the flexible plate. Finally, we apply the modelling techniques in a brief study of upper-airway dynamics wherein soft-palate flutter is considered to be the source of snoring noises. In doing so, we show how a time-varying mean flow influences the type of instability observed as flow speed is increased and demonstrate how localised stiffening can be used to control instability of the flexible plate.     

Stability of wave forms of motion of a viscous fluid layer under tangential stress

We study the stability of wave flow of a viscous incompressible fluid layer subjected to tangential stress and an inclined gravity force with respect to long-wave disturbances.An asymptotic solution is constructed for the equations of the disturbed motion and the problem is reduced to the study of a second-order ordinary differential equation. It is shown that after loss of stability by a Poiseuille flow the laminar nature of the flow is not destroyed, but the form of the free surface acquires a wave-like profile. The Poiseuille regime is stable for low Reynolds numbers. The critical Reynolds number for wave flow is found, and the stability and instability regions are determined.     

Analysis of coupled-mode flutter of pipes conveying fluid on the elastic foundation

ＩｎｔｒｏｄｕｃｔｉｏｎＦｌｕｉｄｉｎｄｕｃｅｄｖｉｂｒａｔｉｏｎｅｘｉｓｔｓｉｎｍａｎｙｅｎｇｉｎｅｅｒｉｎｇｆｉｅｌｄｓ.Ｔｈｅｖｉｂｒａｔｉｏｎａｎｄｓｔａｂｉｌｉｔｙｏｆｐｉｐｅｃｏｎｖｅｙｉｎｇｆｌｕｉｄｉｓａｔｙｐｉｃａｌｅｘａｍｐｌｅ.Ｍａｎｙｓｃｈｏｌａｒｓａｔｈｏｍｅａｎｄａｂｒｏａｄｈａｖｅａｌｗａｙｓｂｅｅｎｉｎｔｅｒｅｓｔｅｄｉｎｔｈｉｓｓｕｂｊｅｃｔａｎｄｍａｄｅａｌｏｔｏｆｓｔｕｄｉｅｓｏｆｉｔ.Ｐａｒｔｉｃｕｌａｒｌｙｄｕｒｉｎｇｒｅｃｅｎｔｄｅｃａｄｅｓ,ｓｏｍｅｒｅ…     

Stability Analysis of Maxwell Viscoelastic Pipes Conveying Fluid with Both Ends Simply Supported

ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｓｉｍｐｌｙｓｕｐｐｏｒｔｅｄｐｉｐｅｓｃｏｎｖｅｙｉｎｇｆｌｕｉｄａｒｅｎａｍｅｄａｓｇｙｒｏｓｃｏｐｉｃｃｏｎ ｓｅｒｖａｔｉｖｅｓｙｓｔｅｍｂｅｃａｕｓｅｉｔｓｅｎｅｒｇｙａｔｔｈｅｅｘｉｔｉｓｅｑｕａｌｔｏｔｈａｔａｔｔｈｅｅｎｔｅｒ[1].Ｔｈｉｓｓｙｓｔｅｍｗａｓｓｔｕｄｉｅｄｂｙｓｏｍｅｓｃｈｏｌａｒｓａｔｈｏｍｅａｎｄａｂｒｏａｄ .Ｐａｉｄｏｕｓｓｉｓ[2 ]ｓｔｕｄｉｅｄｔｈｅｐｒｏｂｌｅｍｏｆｄｙｎａｍｉｃｓａｎｄｓｔａｂｉ…     

Flow-induced flutter instability of cantilever carbon nanotubes

Carbon nanotubes are finding significant application to nanofluidic devices. This work studies the influence of internal moving fluid on free vibration and flow-induced flutter instability of cantilever carbon nanotubes based on a continuum elastic model. Since the flow-induced vibration of cantilever pipes is non-conservative in nature, cantilever carbon nanotubes conveying fluid are damped with decaying amplitude for flow velocity below a certain critical value. Beyond this critical flow velocity, flutter instability occurs and vibration becomes amplified with growing amplitude. Our results indicate that internal moving fluid substantially affects vibrational frequencies and the decaying rate of amplitude especially for longer cantilever carbon nanotubes of larger innermost radius at higher flow velocity, and the critical flow velocity for flutter instability in some cases may fall within the practical range. On the other hand, a moderately stiff surrounding elastic medium (such as polymers) can significantly suppress the effect of internal moving fluid on vibrational frequencies and suppress or eliminate flutter instability within the practical range of flow velocity.     

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.     

Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profiles

Considering the non-uniformity of the flow velocity distribution in fluid-conveying pipes caused by the viscosity of real fluids, the centrifugal force term in the equation of motion of the pipe is modified for laminar and turbulent flow profiles. The flow-profile-modification factors are found to be 1.333, 1.015–1.040 and 1.035–1.055 for laminar flow in circular pipes, turbulent flow in smooth-wall circular pipes and turbulent flow in rough-wall circular pipes, respectively. The critical flow velocities for divergence in the above-mentioned three cases are found to be 13.4%, 0.74–1.9% and 1.7–2.6%, respectively, lower than that with plug flow, while those for flutter are even lower, which could reach 36% for the laminar flow profile. By introducing two new concepts of equivalent flow velocity and equivalent mass, fluid-conveying pipe problems with different flow profiles can be solved with the equation of motion for plug flow.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

Stability of plane poiseuille flow of a highly elastic liquid

The stability of fully developed pressure driven plane laminar flow of a Maxwell fluid has been studied using linear hydrodynamic stability theory. Elasticity is destabilizing in the inertial regime, but the flow is found to be stable to infinitesimal disturbances at low Reynolds numbers. This result contradicts previous calculations, which predicted a low Reynolds number flow instability at a critical recoverable shear of order unity. The previous calculations were carried out using less accurate numerical methods; the eigenvalue problem which must be solved is a delicate one, requiring sophisticated umerical techniques in order to avoid the calculation of spurious unstable modes.This work has direct bearing on the question of the mechanism of a low Reynolds number extrusion instability known as “melt fracture”. It is observed that the intensity of melt fracture increases with increasing die length for high density polyethylene, and it is therfore believed by some experimentalists that fully-developed die flow is unstable for this polymer above a critical recoverable shear. The analysis appears to be at variance with this interpretation of the experimental results.     

In-plane dynamics of a fluid-conveying corrugated pipe supported at both ends

The dynamics and stability of fluid-conveying corrugated pipes are investigated. The flow velocity is assumed to harmonically vary along the pipe rather than with time. The dimensionless equation is discretized with the differential quadrature method(DQM). Subsequently, the effects of the mean flow velocity and two key parameters of the corrugated pipe, i.e., the amplitude of the corrugations and the total number of the corrugations, are studied. The results show that the corrugated pipe will lose stability by flutter even if it has been supported at both ends. When the total number of the corrugations is sufficient, this flutter instability occurs at a micro flow velocity. These phenomena are verified via the Runge-Kutta method. The critical flow velocity of divergence is analyzed in detail. Compared with uniform pipes, the critical velocity will be reduced due to the corrugations, thus accelerating the divergence instability. Specifically,the critical flow velocity decreases if the amplitude of the corrugations increases. However, the critical flow velocity cannot be monotonously reduced with the increase in the total number of the corrugations. An extreme point appears, which can be used to realize the parameter optimization of corrugated pipes in practical applications.     

温度场作用下输流悬臂碳纳米管的颤振失稳分析

以非局部弹性理论为基础,采用欧拉-伯努利梁模型,考虑碳纳米管的小尺度效应,应用哈密顿原理获得了温度场作用下的输流悬臂单层碳纳米管（SWCNT）的振动控制方程以及边界条件,依靠微分变换法（DTM法）对此高阶偏微分方程进行求解,通过数值计算研究了温度场中悬臂单层输流碳纳米管的振动与颤振失稳问题。结果表明：管内流体流速、温度场中温度变化情况与小尺度参数都会对系统振动频率以及颤振失稳临界流速产生影响。其中,小尺度效应将会降低悬臂输流系统的稳定性,使系统更为柔软;而高温场与低温场对系统动态失稳的影响不同,低温场中随温度变化值的增加,系统的稳定性提高;高温场这一作用效果恰好与之相反。     

Effect of transverse stream velocity in an incompressible boundary layer on the stability of the laminar flow regime

The stability of the laminar flow regime in the boundary layer developed on a wall is increased considerably by the relatively slight extraction of fluid from the wall [1–4]. In the theoretical study of this phenomenon, all the investigators known to the present authors have taken into account only the increase in the fullness of the velocity profile in the boundary layer with suction. Computations of the stability characteristics have been made on the assumption that there are no transverse velocities in the laminar boundary layer.We present below an analysis of the stability of the laminar boundary layer in the presence of a constant transverse velocity in the near-wall region (suction). The calculations made show the existence for a given velocity profile in the boundary layer of a relative suction velocity v=v such that with suction velocities greater than v the flow remains stable at all Reynolds numbers, while the method used in the cited references gives a definite finite critical Reynolds number, equal in our notation to the Reynolds number at v=0, for each relative suction velocity.It was found that with suction of fluid from the boundary layer the region of instability has finite dimensions, i.e., there exist lower and upper critical Reynolds numbers. The flow is stable if its Reynolds number is less than the lower, or greater than the upper values of the critical Reynolds number.The instability region diminishes with increase in the relative suction velocity, and at a value of this velocity which is specific for each value of the velocity profile the instability region degenerates into a point-the flow becomes absolutely stable. Thus, with distributed suction it is advisable to increase the relative suction velocity only to a definite magnitude corresponding to disappearance of the instability region. The computational results presented make it possible to estimate this velocity for velocity profiles ranging from a Blasius profile to an asymptotic profile. Specific calculations were made for a family of Wuest profiles, since under actual conditions with suction there always exists a starting segment of the boundary layer [1, 2].     

Aeroelastic analysis of an idealized landing gear door in subsonic potential flow

Landing gear doors on aircraft have experienced flutter during preliminary flight testing. While designs vary widely, landing gear doors are typically plate-like structures with a relatively rigid actuator attached to their inside surface. To better understand the aeroelasticity of landing gear doors, this study investigates the aeroelastic stability of an idealized model. The model consists of a hinged plate with an interior constraint approximating the actuator attachment. The plate is subject to uniform flow, and an unsteady vortex lattice model is coupled to the structural model to predict critical flow velocities. The location and footprint area of the internal constraint, along with plate aspect and mass ratios, are varied to investigate a large parameter space. Results reveal that the critical flow speed and instability mechanism are sensitive to the postulated actuator placement. In general, flutter is the dominant mode of instability when the actuator is postulated in the leading quarter of the plate. In other postulated locations, divergence dominates. However, the exact shape and location of the boundary between flutter and divergence is configuration dependent and found to be especially sensitive to changes in aspect ratio.     

Non-conservative stability of intermediate spring supported columns with an elastically restrained end support

The non-conservative stability of an intermediate spring supported uniform column clastically restrained at one end and subjected to a follower force at the other unsupported end is studied. It is found that when the intermediate spring support is far from the unsupported end, the instability mechanism is flutter. As the intermediate spring support approaches the unsupported end, the instability mechanism is changed from flutter to divergence with the increase of intermediate spring stiffness. For the hinged-intermediate and guided-intermediatc spring supported columns, the critical buckling load of flutter instability will first decrease, then increase as the intermediate spring stiffness is increased. Nevertheless, when the instability mechanism is divergence, the critical buckling load depends on the location of the intermediate spring support only, whereas for the clamped-intermediate spring supported column the critical buckling load of divergence instability decreases monotonically to a fixed value as the intermediate spring stiffness is increased. Finally, the influence of elastic end restraints on the stability of the column is also investigated.     

Flow instability in a curved porous channel formed by two concentric cylindrical surfaces

Stability of laminar flow in a curved channel formed by two concentric cylindrical surfaces is investigated. The channel is occupied by a fluid saturated porous medium; the flow in the channel is driven by a constant azimuthal pressure gradient. The momentum equation takes into account two drag terms: the Darcy term that describes friction between the fluid and the porous matrix, and the Brinkman term, which allows imposing the no-slip boundary condition at the channel walls. An analytical solution for the basic flow velocity is obtained. Numerical analysis is carried out using the collocation method to investigate the onset of instability leading to the development of a secondary motion in the form of toroidal vortices. The dependence of the critical Dean number on porosity and the channel width is analyzed.     

On the Lower Critical Reynolds Number for Flow in a Circular Pipe

Flow in a circular pipe is investigated experimentally at Reynolds numbers higher than that at which the resistance coefficients calculated from the Blasius formula for laminar flow and from the Prandtl formula for turbulent flow are equal. The corresponding Reynolds number based on the mean-flow velocity and the pipe diameter is about 1000. The experiments were performed at a high level of inlet pulsations produced by feeding gas into the pipe through a hole with a diameter several times smaller than the pipe diameter. In our experiments the critical Reynolds number was determined as the value, independent of the distance from the inlet, at which the ratio of the axial to the mean-flow velocity as a function of the Reynolds number deviated from 2. At the maximum ratio of the pipe cross-sectional area to the area of the hole through which the gas entered the pipe, equal to 26, the critical Reynolds number was about 2300. After a fivefold increase in the hole area the critical Reynolds number increased by approximately 4%.At Reynolds numbers below 2000, after at a high level of the inlet pulsations an almost laminar flow had developed in the pipe, a perturbation was introduced by inserting a diametrically oriented cylindrical rod with a diameter 10–20 times smaller than the pipe diameter. In these experiments, at Reynolds numbers higher than 1000, at a distance from the rod equal to 50 pipe diameters the axial to mean-flow velocity ratio was less than 2, approaching this value again at large distances from the rod. The insertion of the rod led to a decrease in the critical Reynolds number by approximately 12%.     

Numerical simulation of confined laminar jet flows

Two‐dimensional incompressible jet development inside a duct has been studied in the laminar flow regime, for cases with and without entrainment of ambient fluid. Results have been obtained for the flow structure and critical Reynolds number values for steady asymmetric jet development and for the onset of temporal oscillations, at various values of the duct‐to‐jet width ratio (aspect ratio). It is found that at low aspect ratios and Reynolds numbers, jet development inside the duct is symmetric. For larger aspect ratios and Reynolds numbers, the jet flow at steady state becomes asymmetric with respect to the midplane, and for still higher values, it becomes oscillatory with respect to time. When entrainment is present, the instabilities of asymmetric development and temporal oscillations occur at a much higher critical Reynolds number for a given aspect ratio, indicating that the stability of the jet flow is higher with entrainment. Copyright © 2000 John Wiley & Sons, Ltd.     

Local similarity solution for the flow of a “second-grade” viscoelastic fluid above a moving plate

The boundary layer flow of a viscoelastic fluid of the second-grade type over a rigid continuous plate moving through an otherwise quiescent fluid with constant velocity U is studied. Assuming the flow to be laminar and two-dimensional, local similarity solution is found with fluid's elasticity and plate's withdrawal speed as the main variables. Results are presented for velocity profiles, boundary layer thickness, wall skin friction coefficient and fluid entrainment in terms of the local Deborah number. A marked formation of boundary layer is predicted, even at low Reynolds numbers, provided the Deborah number is sufficiently large. The boundary layer thickness and the wall skin friction coefficient are found to scale with fluid's elasticity—both decreasing the higher the fluid's elasticity. The amount of fluid entrained is also predicted to decrease whenever a fluid exhibits elastic behavior.