DAMPING OF VERTICALLY EXCITED SURFACE WAVE IN WEAKLY VISCOUS FLUID

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.     

圆平板挤压膜实验及其参数识别

完成了简谐激励下水中平行圆板的挤压膜振动实验．用最小二乘法识别出非线性粘性挤压膜力模型中的４个系数．不同挤压膜厚、不同频率和不同振幅情况的数值模拟结果与实验值吻合较好．研究结果表明，所用模型可较好地描述粘性挤压膜运动，识别出的系数在一定范围内给出良好精度的数值模拟．     

Fully coupled flow-induced vibration of structures under small deformation with GMRES method

Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual （GMRES） method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.     

激波管实验和果冻实验界面不稳定性数值计算

运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学粘性以及热流部分对流场的影响,将可压缩多介质流体动力学高精度欧拉计算方法MFPPM发展到适用于NS(Nayier Stokes)方程的可压缩多介质粘性流体计算方法MVPPM.通过与界面不稳定性实验结果的比较,来检验计算方法的正确性,并验证计算程序的有效性.主要包括一个激波管实验和两个果冻实验,即英国AWE(Atomic Weapons Establishment)激波管实验和LLNL实验室的果冻环实验以及中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室进行的爆轰驱动下单层果冻界面不稳定性实验.研究结果表明:数值模拟结果与实验测量结果以及对应时刻的实验图象均吻合较好.     

Vibration characteristics of FGM circular cylindrical shells filled with fluid using wave propagation approach

The vibration characteristics of a functionally graded material circular cylindrical shell filled with fluid are examined with a wave propagation approach. The shell is filled with an incompressible non-viscous fluid. Axial modal dependence is approximated by exponential functions. A theoretical study of shell vibration frequencies is analyzed for simply supported-simply supported, clamped-simply supported, and clamped-clamped boundary conditions with the fluid effect. The validity and the accuracy of the present method are confirmed by comparing the present results with those available in the literature. Good agreement is observed between the two sets of results.     

非结构网格生成Bowyer—Watson方法的改进

完善了应用Ｄｅｌａｎｎａｙ原理生成非结构网格的Ｂｏｗｙｅｒ－Ｗａｔｓｏｎ方法。通过求解简化Ｐｏｉｓｓｉｏｎ方程，灵活地控制网格的疏密，通过网格层次的区分，提高了网络的贴体性能；并通过分区，实现了适用于粘流计算的网格生成。     

A HYBRID METHOD FOR ANALYZING THE DYNAMIC RESPONSES OF CAVITIES OR SHELLS BURIED IN AN ELASTIC HALF-PLANE

In this paper, based on a variational formalism which originally proposed by Mei [1] for infinite elastic medium and extended by Yeh, et al. [2,3] for elastic half-plane, a hybrid method which combines the finite element and series expansion method is implemented to solve the diffraction of plane waves by a cavity buried in an elastic half-plane. The finite domain which encloses all inhomogeneities including the cavity can be easily formulated by finite element methods. The unknown boundary data obtained by subtracting the known free fields from the total fields which include the boundary nodal displacements and tractions at the interface between the finite domain and the surrounding elastic half-plane are not independent of each other and can be correlated through aseries repre sentation. Due to the continuity condition at the interface, the same series representation is still valid for the exterior elastic half-plane to represents the scattered wave. The unknown coefficients of this series are treated as generalized coordinates and can be easily formulated by the same variational principle. The expansion function of the series is composed of basis function. Each basis function is constructed from the basis function for an infinite plane by superimposing an additional homogeneous reflective term to satisfy both traction free conditions at ground surface and radiation conditions at infinity. The numerical results are made against those obtained by boundary element methods, and good agreements are found.     

在粘性流体中船行波的速度场

从Navier-Stokes方程出发，导得了水上运动物体船行波的波高和水平速度二重积分的表达式。其次采用Lighthill两步法将上述的积分表达式化粉代数表达式。与海水作为理想流体所得的结果相比，消除或减弱了原有的奇性，提高了船行波的数字资料的精度。     

弱粘性流体中垂直激励表面波的阻尼效应

在竖直振动的圆柱形容器中，将Navier-Stokes方程线性化，利用两时间尺度奇异摄动展开法研究了弱粘性流体的单一自由面驻波运动．整个流场被分为外部势流区和内部边界层区两部分，对两部分区域分别求解，得到包含阻尼项和外驱动影响的线性振幅方程．利用稳定性分析，得到形成稳定表面波的条件，给出了临界曲线．此外，还获得了阻尼系数的解析表达式．最后，将线性阻尼加到理想流体条件下所得到的色散关系中对其进行修正，理论结果证明修正后的驱动频率更加接近实验的结果．通过计算发现，当驱动的频率较低时，流体的粘性对表面波模式选择有重要影响，而表面张力的影响不明显；但当驱动频率较高时，流体的表面张力起主要作用，而流体的粘性影响甚小．     

Characteristics of Flow around an Impulsively Rotating Square Cylinder via LB-DF/FD Method

Flow around an impulsively rotating square cylinder in a viscous fluid in range of 1 ≤ Re ≤ 300 is numerically investigated by the previously developed LB-DF/FD method, which combines the lattilce Boltzmann method （LBM） and direct-forcing fictitious domain （DF/FD） scheme. Results show that in total three kinds of transient characteristics depending on Re are observed： 1 ≤ Re ≤ 20, four vortices arising from the corners of the square cylinder separate from the surfaces and gradually become stable without vortex integration or shedding; 20 ≤ Re 〈 100, vortices integration is observed when they grow long enough, then separated from each other; Re ≥ 100, vortex shedding takes place in this regime. The shedding vortex joins the downstream vortex to form a new one. It is also found that vortex shedding happens more than one time when Re ≥ 160. Furthermore, each vortex shedding induces a fluctuation in the torque exerted on the cylinder.