透平叶栅跨音速流动计算中的新型有限元法

将Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法和多级有限元的思想结合起来，构成了在收敛速度和稳定性两方面均较好的新型有限元算法：多级广义有限元。利用这一方法，分别基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和Ｅｕｌｅｒ方程，研究了透平跨音速叶栅无粘流动和粘性流动，并将计算结果与实验结果作了比较。计算结果表明，本方法是透平机械内部跨音速流动计算的强有力的手段。     

Slosh dynamics of inviscid fluids in two‐dimensional tanks of various geometry using finite element method

This paper brings into focus some of the interesting effects arising from the motion of the liquid free surface due to sloshing in partially filled containers of several geometrical shapes in two dimensions. The slosh characteristics that include frequencies, free surface profiles and the hydrodynamic pressure over the container walls have been reported in this study. The equations of motion of the fluid, considered inviscid, are expressed in terms of the pressure variable alone. It is assumed that the frequency of the exciting oscillation is not in the immediate neighborhood of the natural slosh frequency, so that the slope of the free surface is small. Simple harmonic oscillation and earthquake excitations are used as the prescribed boundary conditions. A finite difference‐based iterative time‐stepping technique is employed to advance the solution in the time domain. The paper presents numerical solutions for rectangular, vertically mounted annular cylindrical, trapezoidal and horizontal circular cylindrical containers. Numerical results obtained are compared with the available existing solutions to validate the code developed. The parametric study of the slosh dynamic systems shows the importance of the nature of excitation, fluid height and the geometry of the container. Copyright © 2007 John Wiley & Sons, Ltd.     

纵向参数激励下平动刚-液耦合系统稳定性

纵向参数激励下液体晃动稳定性是航天器动力学中一个广受关注的问题,然而在以往的研究中没有考虑液体晃动与航天器运动之间的耦合作用对系统参数振动稳定性的影响. 建立了用液体晃动等效单摆模型描述的纵向激励下平动刚-液耦合系统的Mathieu方程,采用摄动法确定了耦合系统1/2亚谐波振动和谐波振动的激励幅-频稳定性边界. 研究发现,液体晃动与主刚体横向运动的耦合作用扩大了参数振动不稳定区, 并使其向高频移动,影响的程度随等效晃动质量的减小而减小;液体晃动模态阻尼对1/2亚谐波振动不稳定区的缩小作用远弱于对谐波振动不稳定区的缩小作用. 对耦合系统第1阶液体晃动模态1/2亚谐波振动响应的研究表明:当纵向激励参数在不稳定区内时,可能引起主刚体的纵横耦合振动现象.      

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.

     

电场驱动射流的弯曲失稳分析

本文对电场驱动射流弯曲问题的动力模型作了理论推导和数值分析.在假定射流具有牛顿粘性并作一维流动的基础上,建立了射流运动的控制方程,采用此方程来研究当带电流体穿过电场时受小扰动而产生的失稳现象.将方程无量纲化后,利用简正模态法对控制方程进行了稳定性分析,而后采用伽辽金法求解方程.结果表明,外电场增强、表面电荷量进一步增多(电荷量较大时)、射流轴向流速增大(其值较大时)、射流半径减小以及电荷的分布不均匀性增大,都可以显著增强射流的短波不稳定性;射流偶数阶频率导致了不稳定性,而且阶数越高,不稳定性增长得越快;频率具有虚部,意味着模态中相差的存在;并且射流在短时问内会失去稳定性.     

Simulation of two-dimensional planar flow of viscoelastic fluid

A numerical scheme based on the Finite Element Method has been developed which uses a relaxation factor in the momentum equation with the stresses being evaluated via a streamwise integration procedure. A constitutive equation introduced by Leonov has been used to represent the rheological behavior of the fluid. The convergence of the scheme has been tested on a 2 : 1 abrupt contraction problem by successive mesh refinement for non-dimensional characteristic shear rates, of 5 and 50 for polyisobutylene Vistanex at 27 °C. The recirculation region is shown to increase in size with non-dimensional characteristic shear rate.Theoretical predictions have been compared with the experimental data which include birefringence and pressure loss measurements. In general, the comparisons have been reasonably good and demonstrates the usefulness of the present numerical scheme and the Leonov constitutive equation to describe real polymer flows.     

Modal part analysis of slosh of the liquid in a cylindrical container in the case of small bond number

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｔｉｏｎｏｆｔｈｅｉｎｖｉｓｃｉｄ,ｉｒｒｏｔａｔｉｏｎａｌａｎｄｉｎｃｏｍｐｒｅｓｓｉｂｌｅｆｌｕｉｄｉｎａｃｙｌｉｎｄｒｉｃａｌｃｏｎｔａｉｎｅｒｗｉｔｈｒａｄｉｕｓａｉｓｇｏｖｅｒｎｅｄｂｙＦｉｇ．１Ｓｌｏｓｈｏｆａｌｉｑｕｉｄｉｎａｃｙｌｉｎｄｒｉｃａｌｃｏｎｔａｉｎｅｒ２＝０　　（ｉｎＶ）,（１）ｒ＝０　　（ｏｎｒ＝ａ）,（２）ｚ＝０　　（ｏｎｚ＝－ｈ）,（３）ηｔ＋·η＝ｚ　　（ｏｎｚ＝η）,（４）ｔ＋１２｜…     

Three‐dimensional sloshing: A consistent finite element approach

This paper presents a new computational methodology based on Legendre's polynomials to predict the slosh and acoustic motion in nearly incompressible fluids in both rigid and flexible structures with free surface. Here, we have used a finite element formulation based on Lagrangian frame of reference to model the fluid motion derived using Hamiltonian equation of the fluid system. We formulated three hexahedral finite elements based on strain fields expressed in terms of extended Legendre's polynomials. Sloshing and acoustic motion of liquid is investigated using these newly formulated elements and inf–sup test is performed on these new elements to check the performance of these elements in modeling sloshing under two severe constraints, namely incompressibility and irrotationality. Comparisons of slosh and acoustic frequencies, and mode shapes with exact solutions are given. Dynamic analysis with earthquake and harmonic kind of forcing function is carried out to validate the formulated hexahedral elements to analyze the sloshing response. Numerical results obtained with these new finite elements, and with the present finite element formulation of the mathematical model agree well with the exact solution and as well as with published experimental literature. Copyright © 2010 John Wiley & Sons, Ltd.     

填隙二阶流体下圆球平行于壁平移的粘性阻力

离散元法是分析散体力学行为的数值方法。存在填隙流体时，颗粒之间或颗粒与壁之间产生的法向挤压力和切向阻力、阻力矩，是湿颗粒离散元法的理论基础。二阶流体是以微小偏离牛顿流体本构而考虑时间影响的一种流体。它具有常粘度，并且第一和第二法向应力差正比于剪切率的平方。根据Reynolds润滑理论，采用小参数法，导出了存在填隙二阶流体时，圆球沿平行于平壁缓慢移动时流体的速度场和压力方程，进而求出切向阻力和阻力矩的解析解。有趣的是在推导时所得的速度场和压力方程形式比牛顿流体要复杂得多，但最终结果表明圆球沿平行于平壁移动时因填隙二阶流体引起的切向阻力和阻力矩与牛顿流体时的结果相同。     

Dewar瓶内液体晃动的近似计算方法

本文研究Dewar瓶内液体的晃动特征问题，将液体晃动的微分方程边值问题转换为具有积分形式的泛函极值问题，在此基础上建立了旋转对称容器液体晃动特征问题的有限元数值计算方法和通用程序，并计算了Dewar球，球腔，圆柱腔液体的晃动固有频率，结果显示了Dewar瓶综合了球腔和圆柱腔液体晃动的特征，它将液体晃动频率控制在一窄带范围之内。利用该程序，本文又计算了Dewar瓶液体晃动的单摆模型，为了便于工程应用，本文提出了一种对Dewar瓶建立等效圆柱容器液体晃动模型的方法，给出了模型的参数变换公式，从而可以通过对圆柱窝器液体晃动的解析计算得到Dewar瓶的液体晃动结果。利用数值计算结果，本文验算了近似计算方法的有效性和近似程序。     

Direct Numerical Simulation of particulate flow with heat transfer

Numerical simulations of heat transfer in non-isothermal particulate flows are important to better understand the flow pattern. The complexity of numerical algorithms coupling the heat and mass transfer and the considerable computational resources required limit the number of such direct simulations that can be reasonably performed. We suggest a Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method to compute the temperature distribution and the heat exchange between the fluid and solid phases. The Boussinesq approximation is considered for the flow/temperature fields coupling. We employ a Finite Element Method (FEM) to solve the fluid flow conservation equations for mass, momentum and energy. The motion of particles is computed by a Discrete Element Method (DEM). On each particle, heat transfer is solved using a FEM. For each class of particles, we generate a single FEM grid and translate/rotate it at each time step to match the physical configuration of each particle. Distributed Lagrange multipliers for both the velocity and temperature fields are introduced to treat the fluid/solid interaction. This work is an extension of the method we proposed in Yu et al. (2006). Two two-dimensional (2D) test cases are proposed to validate the implementation by comparing our computational results with those reported in the literature. Finally, the sedimentation of a single sphere in a semi-infinite channel is presented and the results are discussed.     

A regular boundary element method for fluid flow

The Boundary Element Method is now well established as a valid numerical technique for the solution of field problems, equal to the Finite Element Method in generality and surpassing it in computational efficiency in some cases.¹ In this paper is presented a 'Regular Boundary Element Method' as applied to inviscid laminar fluid flow problems. It involves the formation of a system of regular integral equations obtained by moving the singularity outside the domain of the given problem. It is also shown that non-conforming elements may be used whereby freedoms are not defined at the geometric nodes under the boundary element discretization. A linear element is developed here; higher order variants could easily be defined. Satisfactory numerical results have been obtained using the proposed regular method with both conventional (continuous across the boundary) and non-conforming boundary elements for two-dimensional inviscid laminar fluid flow problems having regular and singular solutions.     

Effect of baffles on sloshing modulated fluid force and torque fluctuations on the dewar driven by gravity gradient acceleration in microgravity

The dynamical behavior of fluids affected by the asymmetric gravity gradient acceleration is studied. In particular the effect of surface tension on partially liquid filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank with and without a baffle have been investigated. Results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient acceleration indicate that the gravity gradient acceleration is equivalent to the combined effect of a twisting force and torsional moment acting on the spacecraft. As viscous force (between liquid and solid interface), and surface tension force (between liquid-vapor-solid interface) greatly contribute to the damping effect of slosh wave excitation, a rotating dewar with baffle provides more areas of liquid-solid and liquid-vapor-solid interfaces than that of a rotating dewar without baffle. Results show that the damping effect provided by the baffle reduces the amplitude of slosh wave excitation, lowers the fluid force, torque, and the moment arm of fluid torque fluctuations than that without baffle, and also lowers the degree of asymmetry in the liquid-vapor distribution.     

流体与结构相互作用问题的FE-SPH耦合模拟

应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题.流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理.为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体.采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟.模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的.     

叶轮机械迴转面流场的有限元近似解方法的研究和数值计算

本文提出一种新的概念:即在应用有限元处理不同类型的流体力学、传热学问题时,应该根据各自问题的物理本质,采用不同型式的基函数(shape function),并在文[1]的基础上,应用有限元近似解方法(FEAM)对直角弯管流动进行了数值计算和精度分析,还对带长短叶片的前向离心通风机的迴转面流场进行了数值计算和分析,在迴转面流场计算中采用了一种新的简化方法,解决了在非矩形单元上求取近似解的困难,计算结果表明在稀疏网格下能保持较高精度,从而明显地提高了计算速度。     

A Simplified Quadrature Element Method to compute the natural frequencies of multispan beams and frame structures

In this paper a new version of the Modified Quadrature Element Method (MQEM) is proposed. Like MQEM, the proposed method overcomes the drawback of the distance δ of the Quadrature Element Method (QEM) without introducing further degrees of freedom at the ends of the element as in the Differential Quadrature Element Method (DQEM), but it makes the computational cost of the stiffness matrix (and the mass matrix) lighter and uses a general procedure to generate the sampling points distribution. The method here presented has been applied to compute the fundamental frequencies of some structures.     

Sloshing dynamics-modulated fluid angular momentum and moment fluctuations driven by orbital gravity gradient and jitter accelerations in microgravity

The dynamical behavior of spacecraft propellant affected by the asymmetric combined gravity gradient and jitter acceleration are studied. In particular the effect of surface tension on partially-filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank has been investigated. Three different cases of orbital accelerations: (a) gravity gradient-dominated, (b) equally weighted between gravity gradients and jitter, and (c) gravity jitter-dominated accelerations are considered. Fluctuations of angular momentum, fluid moment and bubble mass center caused by slosh wave excitations driven by gravity gradient and jitter accelerations are also investigated.     

Nonlinear vibrations of axially moving multi-supported strings having non-ideal support conditions

In this study, nonlinear vibrations of an axially moving multi-supported string have been investigated. The main difference of this study from the others is in that there are non-ideal supports allowing minimal deflections between ideal supports at both ends of the string. Nonlinear equations of the motion and boundary conditions have been obtained using Hamilton’s Principle. Dependence of the equations of motion and boundary conditions on geometry and material of the string have been eliminated by non-dimensionalizing. Method of multiple scales, a perturbation technique, has been employed for solving the equations of motion. Axial velocity has been assumed a harmonically varying function about a constant value. Axially moving string has been investigated in three regions. Vibrations have been examined for three different cases of the velocity variation frequency. Stability has been analyzed and stability boundaries have been established for the principal parametric resonance case. Effects of the non-ideal support conditions on stability boundaries and vibration amplitudes have been investigated.

     

A time splitting fictitious domain algorithm for fluid–structure interaction problems (A fictitious domain algorithm for FSI)

A Finite Element Method in mixed Eulerian and Lagrangian formulation is developed to allow direct numerical simulations of dynamical interaction between an incompressible fluid and a hyper-elastic incompressible solid. A Fictitious Domain Method is applied so that the fluid is extended inside the deformable solid volume and the velocity field in the entire computational domain is resolved in an Eulerian framework. Solid motion, which is tracked in a Lagrangian framework, is imposed through the body force acting on the fluid within the solid boundaries. Solid stress smoothing on the Lagrangian mesh is performed with the Zienkiewicz–Zhu patch recovery method. High-order Gaussian integration quadratures over cut elements are used in order to avoid sub-meshing within elements in the Eulerian mesh that are intersected by the Lagrangian grid. The algorithm is implemented and verified in two spatial dimensions by comparing with the well validated simulations of solid deformation in a lid driven cavity and periodic elastic wall deformation driven by a time-dependent flow. It shows good agreement with the numerical results reported in the literature. In 3-D the method is validated against previously reported numerical simulations of 3-D rhythmically contracting alveolated ducts.     

气固两相流介尺度LBM-DEM模型

LBM-DEM耦合方法通常是指一种颗粒流体系统直接数值模拟算法,即是一种不引入经验曳力模型的计算方法,颗粒尺寸通常比计算网格的长度大一个量级,颗粒的受力通过表面的粘性力与压力积分获得,其优点是能描述每个颗粒周围的详细流场,产生详细的颗粒-流体相互作用的动力学信息,可以探索颗粒流体界面的流动、传递和反应的详细信息及两相相互作用的本构关系,但其缺点是计算量巨大,无法应用于真实流化床过程模拟。本文针对气固流化床中的流体以及固体颗粒间的多相流体力学行为,建立了一种稠密气固两相流的介尺度LBMDEM模型,即LBM-DEM耦合的离散颗粒模型,实现在颗粒尺度上流化床的快速离散模拟。该耦合模型采用格子玻尔兹曼方法(LBM)描述气相的流动和传递行为,离散单元法(DEM)用于描述颗粒相的运动,并利用能量最小多尺度(EMMS)曳力解决气固耦合不成熟问题,以提高其模拟精度。通过经典快速流态化的模拟,验证了介尺度LBM-DEM耦合模型的有效性。模拟结果表明介尺度LBM-DEM模型是一种探索实验室规模气固系统的有力手段。