VIBRATIONS OF CIRCULAR CYLINDRICAL SHELLS WITH NONUNIFORM CONSTRAINTS,ELASTIC BED AND ADDED MASS. PART II: SHELLS CONTAINING OR IMMERSED IN AXIAL FLOW

The model introduced in Part I of the present study is extended to take into account a flowing fluid, a mean radial pressure and initial pre-stress in circular cylindrical shells. The axial flow can be external, internal or annular and is described by the potential theory for inviscid and incompressible fluid. The computer program DIVA has been developed. It takes into account all the following complicating effects on the vibrations of circular cylindrical shells: (i) nonuniform boundary conditions around the shell edges including elastic boundary conditions; (ii) fluid–structure interaction including both flowing and quiescent fluids; (iii) internal, external and annular fluids; (iv) effect of a mean radial pressure and initial pre-stress; (v) elastic bed of partial extension in circumferential and longitudinal directions; (vi) intermediate constraints; (vii) added masses. It can be considered the most complete computer program specifically dedicated to dynamics of circular cylindrical shells. The Flügge theory of shells is used to describe the shell deformations. The system has been proved to be conservative for any combination of boundary conditions with restrained displacement at the shell ends. Numerical results show that shells clamped at the upstream end and simply supported at the downstream end have a larger critical velocity than simply supported shells, solving the paradox of Horáček and Zolotarev.     

Stability analysis of loaded coaxial cylindrical shells with internal fluid flow

We present numerical results for dynamical stability of loaded coaxial shells of revolution interacting with the internal fluid flow. The motion of the incompressible fluid is described in the framework of the theory of frictionless potential flow, whereas the static load acting on the shells is caused by the steady forces of viscous drag arising in the viscous turbulent flow in a closed channel. For shells with different boundary conditions, we study how the stability boundary is affected by the value of the gap between the shells for different versions of the outer shell rigidity and fluid flow. We show that, as in the case of unloaded coaxial shells, there is a significant deviation from the previous numerical and analytical results.     

VIBRATIONS OF CIRCULAR CYLINDRICAL SHELLS WITH NONUNIFORM CONSTRAINTS,ELASTIC BED AND ADDED MASS; PART I: EMPTY AND FLUID-FILLED SHELLS

Circular cylindrical shells with non uniform edge constraints (with zero radial and circumferential displacement) are investigated, including riveted shells. The linear modes of simply supported shells vibrating in vacuo are used as admissible functions, and the solution is obtained with the artificial spring method. The Flügge theory of shells is used and in-plane inertia is retained. Any shell constraint other than simple supports can be studied with the proposed method. Complicating effects due to the contained inviscid fluid, elastic bed of partial axial and angular dimensions, intermediate constraints and added mass are considered. The convergence of the method is numerically investigated and the effect of the number of rivets (clamped arcs) on shell modes is studied.     

Forces on oscillating bodies in viscous fluid

This paper describes a method for determining the fluid forces on oscillating bodies in viscous fluid when the corresponding flow problem has been solved using the finite element method. These forces are characterized by the concept of added mass, added damping and added force. Numerical results are obtained for several example body shapes. Comparison is made with exact analytical results and other finite element results for the limiting cases of Stoke's flow and inviscid flow, and good agreement is obtained. The results for finite values of the body amplitude parameter β show the appearance of added force from the steady streaming component of the flow for asymmetric bodies. Results are also obtained for the associated flow where the fluid remote from a fixed body is oscillating.     

Effect of geometric imperfections on non-linear stability of circular cylindrical shells conveying fluid

Circular cylindrical shells conveying incompressible flow are addressed in this study; they lose stability by divergence when the flow velocity reaches a critical value. The divergence is strongly subcritical, becoming supercritical for larger amplitudes. Therefore the shell, if perturbed from the initial configuration, has severe deformations causing failure much before the critical velocity predicted by the linear threshold. Both Donnell's non-linear theory retaining in-plane displacements and the Sanders-Koiter non-linear theory are used for the shell. The fluid is modelled by potential flow theory but the effect of steady viscous forces is taken into account. Geometric imperfections are introduced and fully studied. Non-classical boundary conditions are used to simulate the conditions of experimental tests in a water tunnel. Comparison of numerical and experimental results is performed.     

A GENERAL FORMULA FOR CALCULATING FORCES ON A 2-D ARBITRARY BODY IN INCOMPRESSIBLE FLOW

In the present paper, a general integral equation is presented to calculate the forces exerted on a two-dimensional (2-D) body of arbitrary shape immersed in unsteady, incompressible flows. By finding the general solutions of a set of Laplace equations with particular boundary conditions, the equation can be simplified to produce a simplified formula for calculating the forces. The simplified formula consists of three parts, representing contributions from different physical phenomena: added mass force and/or inertial force in inviscid flow, the force caused by the deformation of fluid and viscosity and the force caused by the convection of fluid with nonzero circulation. It can be applied to any 2-D arbitrary body in viscous or inviscid, steady or unsteady incompressible flow. As the formula excludes either temporal derivatives of velocity or spatial derivatives of vorticity in the flow field, the numerical errors contained in the numerical solution of velocity and vorticity fields will not be magnified, and therefore the resulting force calculated is more accurate. Most importantly, the formula presents an alternative method for obtaining the added mass of a 2-D body of arbitrary shape accelerating in a fluid. For bodies of simple shape, such as a circle, ellipse and plate, the added masses predicted using the present method are in agreement with that obtained by conventional methods. For bodies of complex shape, the present method only requires the calculation of the first two coefficients of the conformal transformation and cross-sectional area.     

The Opposing Effect of Viscous Dissipation Allows for a Parallel Free Convection Boundary-Layer Flow Along a Cold Vertical Flat Plate

External free convection boundary-layer flows are usually treated by neglecting the effect of viscous dissipation. This assumption always results in a non-parallel flow, besides a strong parallel component also a weak transversal component of the (steady) velocity field occurs. The present paper shows, however, that the weak opposing effect of the buoyancy forces due to heat release by viscous dissipation, can give rise along a cold vertical plate adjacent to a fluid saturated porous medium to a strictly parallel steady free convection flow. This boundary-layer flow is described by an algebraically decaying exact analytical solution of the basic balance equations.     

圆柱薄壳稳定性的一个修正理论

著名的唐乃尔(Donnell)——穆什塔利的简化壳体理论只能较精确地适用于较短圆柱壳稳定性计算.其近似性误差随长度与半径之比的增加而增大.本文考虑了横向切力的影响,对非完善型圆柱壳体推导了几何非线性理论的基本方程,建立了对各种长度半径比的圆柱壳体稳定性均适用的修正理论.     

GENERALIZATION OF KIRCHHOFF KINETIC ANALOGY TO THIN ELASTIC SHELLS 1)

将弹性细杆的"Kirchhoff动力学比拟"方法推广到弹性薄壳,使弹性薄壳的变形在物理概念上和刚体的运动对应, 在数学表述上等同,从而可以用刚体动力学的理论和方法研究弹性薄壳的变形,为连续的弹性薄壳提供新的离散化方法. 在直法线假设下,在弹性中面上构筑空间正交轴系, 此轴系沿坐标线"运动"的角速度构成两自变量的弯扭度. 沿两个坐标线的弯扭度表达了弹性薄壳的变形和位形,证明了弯扭度之间以及弯扭度与中面切矢间的相容关系. 用Euler角和Lam$\acute{e}$系数表达了非完整约束和中面位形的微分方程,用弯扭度和Lam$\acute{e}$系数表达了应变和应力以及内力及其本构方程.导出了用分布内力集度表达的弹性薄壳在变形后位形上的平衡偏微分方程组,方程的形式与刚体动力学的Euler方程和弹性细杆的Kirchhoff方程具有相似性,实现了Kirchhoff动力学比拟对弹性薄壳的推广.总结了弹性薄壳静力学和刚体动力学以及弹性细杆静力学在概念上的比拟关系.最后给出了一个算例. 为研究弹性薄壳的变形和运动提供新的建模方法和研究思路.也可进一步推广到弹性薄壳动力学.     

Kirchhoff动力学比拟对弹性薄壳的推广

将弹性细杆的"Kirchhoff动力学比拟"方法推广到弹性薄壳,使弹性薄壳的变形在物理概念上和刚体的运动对应,在数学表述上等同,从而可以用刚体动力学的理论和方法研究弹性薄壳的变形,为连续的弹性薄壳提供新的离散化方法.在直法线假设下,在弹性中面上构筑空间正交轴系,此轴系沿坐标线"运动"的角速度构成两自变量的弯扭度.沿两个...     

有限长圆柱壳中轴对称弹性瞬态波

有限长的计及剪切变形和转动惯性的弹性圆柱壳的轴对称运动方程经过Laplace变换后转化为一组相空间中的方程。对该方程组作了一些适当的处理后,应用广义射线法,得到了相空间中位移和内力的射线法表达式。采用快速Fourier变换作Laplacl逆变换,即可得到圆柱壳受轴对称冲击载荷时的弹性瞬态波解。     

圆柱壳在侧向非对称冲击载荷下的塑性动力屈曲

采用塑性动力屈曲能量准则，对圆柱壳非均匀局部径向冲击下的塑性动力屈曲进行了研究。文中分别采用弹性的角弹簧和位移弹簧模拟非屈曲部分对屈曲部分的影响，推导了有关模态和临界速度的冲击公式，并与有关的实验结果进行了比较，分析了两种边界对塑性动力屈曲的影响。     

Von Kàrmàn problem for a rotating permeable disk

The von Kàrmàn problem of steady viscous incompressible flow in a half-space beneath a permeable, porous disk of infinite radius, uniformly rotating it its own plane, is considered. The Bevers-Joseph condition which takes account for viscous fluid slip relative to the plane disk is taken for the boundary condition for the tangential velocity component. This makes it possible to consider the flow only outside the porous disk. The disk permeability effect on the boundary layer is studied, in particular, its influence on the dependence of the moment of viscous forces distributed over the disk about the axis of rotation on the permeability and the pressure difference.     

Non-Darcy free convection in a thermally stratified porous medium along a vertical plate with variable heat flux

Free convection along an impermeable vertical plate embedded in a thermally stratified, fluid-saturated porous medium is analyzed. The wall heat flux is varied in a power-law form. The non-Darcian effects, such as solid-boundary viscous resistances, high-flow-rate inertia forces, near wall nonuniform porosity distribution and thermal dispersion, have been considered in the present study. Due to the variation of porosity in the near wall region, the stagnant thermal conductivity also varies accordingly. The nonsimilar system of transformed equations is solved with Keller's Box method. It is shown that the thermal stratification effect and the higher value of the exponent m can increase the local Nusselt number. Also the non-Darcian and thermal dispersion effects significantly influence the velocity and temperature profiles and local Nusselt number.     

NONLINEAR VIBRATION RESPONSE AND BIFURCATION OF CIRCULAR CYLINDRICAL SHELLS UNDER TRAVELING CONCENTRATED HARMONIC EXCITATION

The nonlinear vibration of a cantilever cylindrical shell under a concentrated harmonic excitation moving in a concentric circular path is proposed. Nonlinearities due to large-amplitude shell motion are considered, with account taken of the effect of viscous structure damping. The system is discretized by Galerkin's method. The method of averaging is developed to study the nonlinear traveling wave responses of the multi-degrees-of-freedom system. The bifurcation phenomenon of the model is investigated by means of the averaged system in detail. The results reveal the change process and nonlinear dynamic characteristics of the periodic solutions of averaged equations.     

Study of modeling unsteady blade row interaction in a transonic compressor stage part 1: code development and deterministic correlation analysis

The average-passage equation system (APES) provides a rigorous mathematical framework for accounting for the unsteady blade row interaction through multistage compressors in steady state environment by introducing deterministic correlations (DC) that need to be modeled to close the equation system.The primary purpose of this study is to provide insight into the DC characteristics and the influence of DC on the time-averaged flow field of the APES.In Part 1 of this two-part paper,firstly a 3D viscous unsteady and time-averaging flow CFD solver is developed to investigate the APES technique.Then steady and unsteady simulations are conducted in a transonic compressor stage.The results from both simulations are compared to highlight the significance of the unsteady interactions.Furthermore,the distribution characteristics of DC are studied and the DC at the rotor/stator interface are compared with their spatial correlations (SC).Lastly,steady and time-averaging (employing APES with DC) simulations for the downstream stator alone are conducted employing DC derived from the unsteady results.The results from steady and time-averaging simulations are compared with the time-averaged unsteady results.The comparisons demonstrate that the simulation employing APES with DC can reproduce the time-averaged field and the 3D viscous time-averaging flow solver is validated.     

Efficiency of an auto-propelled flapping airfoil

The present study deals with an investigation of the flow aerodynamic characteristics and the propulsive velocity of a system equipped with a nature inspired propulsion system. In particular, the study is aimed at studying the effect of the flapping frequency on the flow behavior. We consider a NACA0014 airfoil undergoing a vertical sinusoidal flapping motion. In contrast to nearly all previous studies in the literature, the present work does not impose any velocity on the inlet flow. During each iteration the outer flow velocity is computed after having determined the forces exerted on the airfoil. Forward motion may only be produced by flapping motion of the airfoil. This is more consistent with the physical phenomenon. The non-stationary viscous flow around the flapping airfoil is simulated using Ansys-Fluent 12.0.7. The airfoil movement is achieved using the deformable mesh technique and an in-house developed User Define Function (UDF). Our results show the influence of flapping frequency and amplitude on both the airfoil velocity and the propulsive efficiency. The resulting motion is contrasts to the applied forces. In the present study, the frequency ranges from 0.1 to 20 Hz while the airfoil amplitude values considered are: 10%, 17.5%, 25% and 40%.     

Nonlinear vibration of a rotating laminated composite circular cylindrical shell: traveling wave vibration

In this paper, the large-amplitude (geometrically nonlinear) vibrations of rotating, laminated composite circular cylindrical shells subjected to radial harmonic excitation in the neighborhood of the lowest resonances are investigated. Nonlinearities due to large-amplitude shell motion are considered using the Donnell’s nonlinear shallow-shell theory, with account taken of the effect of viscous structure damping. The dynamic Young’s modulus which varies with vibrational frequency of the laminated composite shell is considered. An improved nonlinear model, which needs not to introduce the Airy stress function, is employed to study the nonlinear forced vibrations of the present shells. The system is discretized by Galerkin’s method while a model involving two degrees of freedom, allowing for the traveling wave response of the shell, is adopted. The method of harmonic balance is applied to study the forced vibration responses of the two-degrees-of-freedom system. The stability of analytical steady-state solutions is analyzed. Results obtained with analytical method are compared with numerical simulation. The agreement between them bespeaks the validity of the method developed in this paper. The effects of rotating speed and some other parameters on the nonlinear dynamic response of the system are also investigated.     

Development of a lattice Boltzmann method for two‐layered shallow‐water flow

In this paper the dynamics of a two‐layered liquid, made of two immiscible shallow‐layers of different density, has been investigated within the framework of the lattice Boltzmann method (LBM). The LBM developed in this paper for the two‐layered, shallow‐water flow has been obtained considering two separate sets of LBM equations, one for each layer. The coupling terms between the two sets have been defined as external forces, acted on one layer by the other. Results obtained from the LBM developed in this paper are compared with numerical results obtained solving the two‐layered, shallow‐water equations, with experimental and other numerical results published in literature. The results are interesting. First, the numerical results obtained by the LBM and by the shallow‐water model can be considered as equivalent. Second, the LBM developed in this paper is able to simulate motion conditions on nonflat topography. Third, the agreement between the LBM (and also shallow‐water model) numerical results and the experimental results is good when the evolution of the flow does not depend on the viscosity, that is, during the initial phase of the flow, dominated by gravity and inertia forces. When the viscous forces dominate the evolution of the flow the agreement between numerical and experimental results depends strongly on the viscosity; it is good if the numerical LBM viscosity has the same order of magnitude of the liquid's kinematic viscosity. Copyright © 2012 John Wiley & Sons, Ltd.     

Two-layer flow of magnetic fluids

The plane two-layer flow of viscous incompressible fluids with differentmagnetic properties between two horizontal rigid planes in a nonuniform traveling magnetic field is investigated. An arbitrary nonuniform periodic traveling magnetic field generates wavelike changes in the interface between the media and fluid flows with nonzero flow rates. From the given magnetic field the interface shape and the velocities, pressures and mean flow rates of the fluids are calculated. The cases of a cosine magnetic force and of the magnetic field created by ferromagnetic cylinders moving in a uniform magnetic field are considered. The effect of various parameters on the mean flow rates of the fluids is investigated.