Stability analysis of pipes conveying fluid with fractional viscoelastic model

Divergence and flutter instabilities of pipes conveying fluid with fractional viscoelastic model has been investigated in the present work. Attention is concentrated on the boundaries of the stability. Based on the Euler–Bernoulli beam theory for structural dynamics, viscoelastic fractional model for damping and, plug flow model for fluid flow, equation of motion has been derived. The effects of gravity, and distributed follower forces are also considered. By transferring the equation of motion to the Laplace domain and using the Galerkin method, the characteristic equations are obtained. By solving the eigenvalue problem, frequencies and dampings of the system have been obtained versus flow velocity. Some numerical test cases have been studied with viscoelastic fractional model and the effect of the fractional derivative order and the retardation time is investigated for various boundary conditions.

     

Computation of three-dimensional flow using the Euler equations and a multiple-grid scheme

The unsteady Euler equations are numerically solved using the finite volume one-step scheme recently developed by Ron-Ho Ni. The multiple-grid procedure of Ni is also implemented. The flows are assumed to be homo-enthalpic; the energy equation is eliminated and the static pressure is determined by the steady Bernoulli equation; a local time-step technique is used. Inflow and outflow boundaries are treated with the compatibility relations method of ONERA. The efficiency of the multiple-grid scheme is demonstrated by a two-dimensional calculation (transonic flow past the NACA 12 aerofoil) and also by a three-dimensional one (transonic lifting flow past the M6 wing). The third application presented shows the ability of the method to compute the vortical flow around a delta wing with leading-edge separation. No condition is applied at the leading-edge; the vortex sheets are captured in the same sense as shock waves. Results indicate that the Euler equations method is well suited for the prediction of flows with shock waves and contact discontinuities, the multiple-grid procedure allowing a substantial reduction of the computational time.     

多场作用下输流单层碳纳米管的动力学特性

以非局部弹性理论为基础,采用欧拉-伯努利梁模型,考虑管型区域内滑移边界条件以及碳纳米管的小尺度效应,应用哈密顿原理获得了温度场与轴向磁场共同作用下的输流单层固支碳纳米管（SWCNT）的振动控制方程以及边界条件,依靠微分变换法（DTM法）对此高阶偏微分方程进行求解,通过数值计算研究了多场中单层固支输流碳纳米管的振动与失稳问题。结果表明：温度场、轴向磁场强度、Knudsen数及小尺度参数都会对系统振动频率以及失稳临界流速产生影响。     

输液管的非线性振动、分叉与混沌——现状与展望

围绕输液管的非线性振动、分叉与混沌问题，对近几年来的主要研究工作加以综述并提出预测，其中包括运动方程中非线性项的归纳与讨论、非线性动力分析的一些现代计算方法、定常流速下输液管的分叉与混沌行为、振荡流速下输液管的参数共振以及今后值得进一步研究的某些问题．     

破碎岩体非等温渗流的非线性动力学研究

分别从固体及流体导热的能量方程出发,导出破碎岩体非等温渗流的能量本构方程, 结合渗流的连续性方程、运动方程、状态方程等建立了破碎岩体非等温渗流的一维非线性动力学方程组;结合Mathcad软件计算得到了系统的无量纲化平衡态, 利用逐次亚松弛迭代法分析了对应于不同参数时平衡态的稳定性;指出非等温渗流系统存在鞍结分岔及折叠突变, 与等温渗流相比, 考虑温度场的破碎岩体渗流动力系统更容易发生渗流突变.      

An investigation on internal thermal flow of non-Newtonian fluid

In the present paper an unsteady thermal flow of non-Newtonian fluid is investigated which is of the fiow into axisymmetric mould cavity. In the second part an unsteady thermal flow of upper-convected Maxwell fluid is studied, For the flow into mould cavity the constitutive equation of power-law fluid is used as a rheological model of polymer fluid. The apparent viscosity is considered as a function of shear rate and temperature. A characteristic viscosity is introduced in order to avoid the nonlinearity due to the temperature dependence of the apparent viscosity. As the viscosity of the fluid is relatively high the flow of the thermal fluid can be considered as a flow of fully developed velocity field. However, the temperature field of the fluid fiow is considered as an unsteady one. The governing equations are constitutive equation, momentum equation of steady flow and energy conservation equation of non-steady form. The present system of equations has been solved numerically by the splitting differen     

Generalized potential flow theory and direct calculation of velocities applied to the numerical solution of the Navier-Stokes and the Boussinesq equations

A formulation based on three scalar functions or potentials is applied to analyse the Navier-Stokes and Boussinesq equations in three dimensions. In this formulation an explicit expression for the pressure exists, the so-called generalized Bernoulli equation. Therefore the scalar functions formulation may be considered as a generalization of the well-known potential flow and Bernoulli theory for irrotational fluid motion. The many advantages of this formulation applied to three-dimensional Navier-Stokes and Boussinesq flow will be discussed, and a numerical example is given as an illustration.     

Dynamics of a pipe conveying fluid flexibly restrained at the ends

The subject of this paper is the study of dynamics and stability of a pipe flexibly supported at its ends and conveying fluid. First, the equation of motion of the system is derived via the extended form of Hamilton׳s principle for open systems. In the derivation, the effect of flexible supports, modelled as linear translational and rotational springs, is appropriately considered in the equation of motion rather than in the boundary conditions. The resulting equation of motion is then discretized via the Galerkin method in which the eigenfunctions of a free-free Euler–Bernoulli beam are utilized. Thus, a general set of second-order ordinary differential equations emerges, in which, by setting the stiffness of the end-springs suitably, one can readily investigate the dynamics of various systems with either classical or non-classical boundary conditions. Several numerical analyses are initially performed, in which the eigenvalues of a simplified system (a beam) with flexible end-supports are obtained and then compared with numerical results, so as to verify the equation of motion, in its simplified form. Then, the dynamics of a pinned-free pipe conveying fluid is systematically investigated, in which it is found that a pinned-free pipe conveying fluid is generally neutrally stable until it becomes unstable via a Hopf bifurcation leading to flutter. The next part of the paper is devoted to studying the dynamics of a pinned-free pipe additionally constrained at the pinned end by a rotational spring. A wide range of dynamical behaviour is seen as the mass ratio of the system (mass of the fluid to the fluid+pipe mass) varies. It is surprising to see that for a range of rotational spring stiffness, an increase in the stiffness makes the pipe less stable. Finally, a pipe conveying fluid supported only by a translational and a rotational spring at the upstream end is considered. For this system, the critical flow velocity is determined for various values of spring constants, and several Argand diagrams along with modal shapes of the unstable modes are presented. The dynamics of this system is found to be very complex and often unpredictable (unexpected).     

黏弹流体流动的数值模拟研究进展

综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。     

Subsonic-Sonic Limit of Approximate Solutions to Multidimensional Steady Euler Equations

A compactness framework is established for approximate solutions to subsonic-sonic flows governed by the steady full Euler equations for compressible fluids in arbitrary dimension. The existing compactness frameworks for the two-dimensional irrotational case do not directly apply for the steady full Euler equations in higher dimensions. The new compactness framework we develop applies for both non-homentropic and rotational flows. One of our main observations is that the compactness can be achieved by using only natural weak estimates for the mass balance and the vorticity, along with the Bernoulli law and the entropy relation, through a more delicate analysis on the phase space. As direct applications, we establish two existence theorems for multidimensional subsonic-sonic full Euler flows through infinitely long nozzles.     

On the case when steady converging/diverging flow of a non-Newtonian fluid in a round cone permits an exact solution

This communication considers the steady converging/diverging flow of a non-Newtonian viscous power-law fluid in a round cone. The motion is driven by a sink/source of mass at the origin. It is shown that the problem permits exact similarity solution for a particular value (n=4/3) of the fluid index. In this case a complete set of governing equations can be reduced to an ordinary differential equation, which is solved numerically for different values of the main non-dimensional parameters (the cone angle and the dimensionless sink/source intensity).     

Magnetohydrodynamic transient flows of a non-Newtonian fluid

Some exact solutions of the time-dependent partial differential equations are discussed for flows of an Oldroyd-B fluid. The fluid is electrically conducting and incompressible. The flows are generated by the impulsive motion of a boundary or by application of a constant pressure gradient. The method of Laplace transform is applied to obtain exact solutions. It is observed from the analysis that the governing differential equation for steady flow in an Oldroyd-B fluid is identical to that of the viscous fluid. Several results of interest are obtained as special cases of the presented analysis.     

广义Maxwell黏弹性流体在两平板间的非定常流动

将分数阶微积分运算引入Maxwell黏弹性流体的本构方程，研究了黏弹性流体在两平板问的非定常流动．对于广义Maxwell黏弹性流体的分数阶导数模型，导出了对时间具有分数阶导数的特殊运动方程，利用分数阶微积分的Laplace变换理论，得到了流动的解析解．     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.     

Correction of local elasticity for nonlocal residuals: application to Euler–Bernoulli beams

Complications exist when solving the field equation in the nonlocal field. This has been attributed to the complexity of deriving explicit forms of the nonlocal boundary conditions. Thus, the paradoxes in the existing solutions of the nonlocal field equation have been revealed in recent studies. In the present study, a new methodology is proposed to easily determine the elastic nonlocal fields from their local counterparts without solving the field equation. This methodology depends on the iterative-nonlocal residual approach in which the sum of the nonlocal fields is treaded as a residual field. Thus, in this study the corrections of the local linear and nonlinear elastic fields for the nonlocal residuals in materials are presented. These corrections are formed based on the general nonlocal theory. In the context of the general nonlocal theory, two distinct nonlocal parameters are introduced to form the constitutive equations of isotropic elastic continua. In this study, it is demonstrated that the general nonlocal theory outperforms Eringen’s nonlocal theory in accounting for the impacts of the material’s Poisson’s ratio on its mechanics. To demonstrate the effectiveness of the proposed approach, the corrections of the local static bending, vibration, and buckling characteristics of Euler–Bernoulli beams are derived. Via these corrections, bending, vibration, and buckling behaviors of simple-supported nonlocal Euler–Bernoulli beams are determined without solving the beam’s equation of motion.     

Stability of fluid conveying nanobeam considering nonlocal elasticity

In this study, the nonlocal Euler–Bernoulli beam theory is employed in the vibration and stability analysis of a nanobeam conveying fluid. The nanobeam is assumed to be traveling with a constant mean velocity along with a small harmonic fluctuation. In the considered analysis, the effects of the small-scale of the nanobeam are incorporated into the equations. By utilizing Hamilton’s principle, the nonlinear equations of motion including stretching of the neutral axis are derived. Damping effect is considered in the analysis. The closed form approximate solution of nonlinear equations is solved by using the multiple scale method, a perturbation technique. The effects of the different value of the nonlocal parameters, mean speed value and ratios of fluid mass to the total mass as well as effects of the simple–simple and clamped–clamped boundary conditions on the linear and nonlinear frequencies, stability, frequency–response curves and bifurcation point are presented numerically and graphically. The solvability conditions are obtained for the three distinct cases of velocity fluctuation frequency. For all cases, the stability areas of system are constructed analytically.     

油藏数值模拟发展现状

油藏数值模拟是以多相流体在多孔介质中流动的理论为基础,目前已广泛应用于油藏工程.本文就流体流动方程、数值计算方法、软件的开发和应用等方面,阐述了所谓黑油模型发展现状,同时介绍了组分模型和热采模型等其他类型模型.还就国内如何开展油藏模拟工作提出了建议.      

Fluid structure interaction of cantilever micro and nanotubes conveying magnetic fluid with small size effects under a transverse magnetic field

Micro and nanotubes have found major application in fluidic systems as channels for conveying fluid. In some micro and nanofluidic applications such as drug delivery, a transverse magnetic field can be used to guide the fluid flow by generating an axial force in the flow direction. An important issue in the design of micro and nanofluidic systems is the structural vibration caused by the fluid flow. In the current study, we investigate the effect of transverse magnetic field on the vibration of cantilever micro and nanotubes conveying fluid by considering the small size effects. We couple the nonlocal Euler–Bernoulli beam model with Navier–Stokes theory to determine a fluid structure interaction (FSI) model for the vibration analysis of the system. We modify the FSI governing equation by driving a velocity correction factor to consider the effect of transverse magnetic field on the fluid flow’s pattern through the tube. Then, we use the Galerkin’s method to obtain the frequency diagrams for the instability analysis of the system. We show that the transverse magnetic field can have a substantial effect on the dynamics of tube conveying fluid by increasing the system’s natural frequencies and critical flow velocity which contributes to the flutter instability. We also discover that although the transverse magnetic field plays a crucial role on dynamics of microstructures, its effect on the dynamics of nanotubes is not significant and can be ignored.     

THE TEACHING OF FLUID MECHANICS WITH REPRESENTATION OF FLUID FLOW

为了在课堂教学中清晰阐述流体的运动,在作者多年的流体力学课程教学经验和体会的基础上,应用连续介质假设,对流体力学课程中关于流体运动描述的Lagrange 方法和Euler 方法及其相互转换的教学方法进行了探讨。将抽象的概念,形象地从时空角度,讲述流体作为连续介质的运动规律。分析结果能够对学生学习这些内容提供帮助,同时为其他教师对流体力学课程的讲授提供参考。