超疏水小球低速入水空泡研究

物体入水问题是一类复杂的流固耦合问题,具有广泛的工程应用背景.物体在跨越自由液面入水的过程中,在一定的条件下,会向水中卷入空气形成空泡,空泡的运动还可能形成指向物体的射流,从而对物体的受力及其运动过程产生影响.超疏水表面能够在物体入水过程中形成多尺度流固耦合作用,进而影响物体的运动和宏观流动现象.而对于小尺度的小球低速入水问题,表面和界面力往往起主导作用.为了在更广的参数空间获得超疏水小球入水空泡类型和小球的运动特性,采用高速摄影实验方法,研究了半径0.175$\sim$10mm的超疏水小球低速入水及空泡动力学行为,获得了小球漂浮振荡、准静态空泡、浅闭合空泡、深闭合空泡和表面闭合空泡5种类型的动力学行为,探讨了这些运动行为与韦伯数We}和邦德数Bo之间的关系,并推导了小球漂浮振荡与下沉现象的无量纲关系.研究结果表明:超疏水小球的入水及空泡动力学行为主要与韦伯数We和邦德数Bo有关.在邦德数Bo $<$ $O$ (10$^{-1})$范围内,表面张力对流动的影响显著,随着韦伯数We}的增大,小球入水及空泡动力学行为依次经历漂浮振荡、准静态闭合、浅闭合、深闭合和表面闭合;在邦德数$O$ (10$^{-1})$<$ Bo} $<$O(1)$范围内,漂浮振荡现象不再发生;当邦德数$Bo>O(1)$后,浅闭合现象也不再发生;小球漂浮振荡与下沉现象的临界关系可以用相似律关系描述.      

旋转圆球入水空泡特性与流场结构的大涡模拟研究

圆球旋转入水过程对于基于先导物投放的新型入水降载方式具有重要研究价值. 采用大涡模拟方法结合均质多相流模型和VOF界面捕捉算法, 对低弗劳德数条件下疏水圆球高速旋转入水的自由运动过程进行了数值模拟, 研究了转速对入水空泡演化、流场结构和水动力特性的影响. 采用动网格与滑移网格技术实现圆球的自由运动, 并基于试验结果对比验证了数值模拟的可靠性与正确性. 旋转运动的升力效应导致圆球入水弹道发生偏转并从水面携带横向楔形射流侵入空泡内部. 采用入水砰击速度与转速进行归一化分析, 结果表明入水转速的增加显著改变了圆球的动力特性: 水平方向的速度和位移以及升力峰值都随入水转速的增加而变大, 但升力峰值受到入水速度的限制; 而垂直方向的速度和加速度以及空泡断裂深度几乎不受转速增加的影响, 并且空泡深闭合发生前圆球转速变化不大. 入水转速的增加也使液面飞溅环和空泡断裂的非对称性增强, 在较低转速时发生空泡表面闭合, 而在较高转速时则发生空泡深闭合. 对于空泡深闭合模式, 入水转速的增加带来更强的横向楔形射流, 并且抑制了空泡断裂产生的高压以及相应涡结构的生成, 致使圆球在入水砰击阶段承受更低的侧向压力.      

A fast mesh deformation method for marine propeller flow

ABSTRACT

A fast mesh deformation method for propeller flow is developed based on the elastic solid method. The flow field of a propeller is assumed to be fulfilled with a kind of pseudo elastic solid which does not influence the flow. The vibration equation for the propeller blade-pseudo elastic solid system is derived. During fluid-structure coupling, the nodal displacement for the blade and the flow mesh is computed by modal superposition of the first several modes. Fluid-structure coupling is performed for a highly skewed propeller. The computing time for the dynamic mesh by the present method is about 0.017% of the computing time by the existing elastic solid method. The computing time for the fluid-structure coupling using the present method is 52% less than the computing time by the existing elastic solid method.     

Force reconstruction based on the union of singular spectral analysis and regularization method

Abstract

To overcome the problem of some ill-posed inverse problem of force reconstruction, which is caused by the noise in the measured responses and small singular values of the structure, a technology of force reconstruction based on a hybrid method of singular spectral analysis (SSA) and the Landweber regularization method is proposed in this study for the first time. The SSA is used to filter the structural response before using Landweber regularization. A new choice method of phase space reconstruction dimension is theoretically proposed, and the minimum embedding dimension is determined by the concept of optimizing difference spectrum theory. The feasibility of this method was demonstrated through three kinds of force reconstructions. The numerical simulation results and an acoustic vibration experiment demonstrated that the proposed method is more effective than the traditional method.

Communicated by Wei-Chau Xie.     

Asymptotic modeling of the impact of a spherical indenter on an elastic half-space

The impact of a rigid sphere on a homogeneous, isotropic elastic half-space in the absence of friction and adhesion is considered. The influence of the superseismic stage immediately following the moment of first contact upon the impact process is investigated in the frame of the Hertzian impact theory. The first order asymptotic approximation for the contact force in a three-dimensional dynamic contact problem with the slowly moving contact zone boundary is obtained and the corresponding asymptotic model of impact is developed. The motion of the indenter as it indents and rebounds from the elastic medium is analytically described. Explicit formulas are derived for the peak indentation depth, contact time, and rebound velocity as functions of the initial impact velocity, indenter mass, and characteristics of the elastic half-space.     

Implicit Calculations of an Aeroelasticity Problem

The purpose of this work is to show that a linearized implicit scheme for the flow resolution can be an efficient and accurate method for solving fluid-structure interaction. The fluid is modeled by the Euler equations in two dimensions and the structure by a one (free piston) or a two (NACA0012 airfoil) degrees of freedom system. The schemes are developed using a finite volume/finite element formulation and, stating the moving boundary problem in the space-time domain, the Riemann solver is generalized in a suitable manner. Assuming a modal decomposition for the structure's response, an analytical solution to the equation of motion is obtained.

The effects of the linearized implicit scheme on the aeroelastic response are demonstrated on the free piston and the NACA 0012 airfoil problems. In the latter case, we focus on the capability of the linearized implicit scheme to accurately predict the stability limit of the coupled response (wing flutter analysis). Although the above analysis is performed using a rigid transformation, a robust moving mesh strategy is presented for more general 2-D and 3-D deformations.     

Steady Flow of a Navier–Stokes Liquid Past an Elastic Body

We perform a mathematical analysis of the steady flow of a viscous liquid, L{\mathcal{L}} , past a three-dimensional elastic body, B{\mathcal{B}} . We assume that L{\mathcal{L}} fills the whole space exterior to B{\mathcal{B}} , and that its motion is governed by the Navier–Stokes equations corresponding to non-zero velocity at infinity, v _∞. As for B{\mathcal{B}} , we suppose that it is a St. Venant–Kirchhoff material, held in equilibrium either by keeping an interior portion of it attached to a rigid body or by means of appropriate control body force and surface traction. We treat the problem as a coupled steady state fluid-structure problem with the surface of B{\mathcal{B}} as a free boundary. Our main goal is to show existence and uniqueness for the coupled system liquid-body, for sufficiently small |v _∞|. This goal is reached by a fixed point approach based upon a suitable reformulation of the Navier–Stokes equation in the reference configuration, along with appropriate a priori estimates of solutions to the corresponding Oseen linearization and to the elasticity equations.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis

Based on linearized 2-D Navier-Stokes equation, a Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis is proposed. Under assumption of incompressibility for the fluid, the corresponding equivalent boundary integral equation in terms of the potential function and stream function is first established by Lamb's transform in the Laplace transform domain. It enables us to solve impact water problems in frequency domain by the boundary element method, in which the effect of viscous flow on the dynamic response can be taken into account. Then a complete solution of the problem under consideration in time domain is obtained by means of Durbin's formulas for the numerical inversion of the Laplace transform. Finally, a practical example is given to confirm the validity of the present method. Project supported by the National Defence Foundation of Science & Technology of China (No. J14. 8. 1. JW0515).     

ANALYSIS OF DYNAMIC RESPONSE OF AN IMPACTED ELASTIC PLATE

ＡＮＡＬＹＳＩＳＯＦＤＹＮＡＭＩＣＲＥＳＰＯＮＳＥＯＦＡＮＩＭＰＡＣＴＥＤＥＬＡＳＴＩＣＰＬＡＴＥ￥(尹邦信)ＹｉｎＢａｎｇｘｉｎ（ＤｅｐａｒｔｍｅｎｔｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇ，ＳｏｕｔｈｗｅｓｔｉｎｓｔｉｔｕｔｅｏｆＴｅｃｈｎｏｌｏｇｙ，...     

A simplified method for time domain simulation of cross-flow vortex-induced vibrations

A new method for time domain simulation of cross-flow vortex-induced vibrations of slender circular cylindrical structures is developed. A model for the synchronization between the lift force and structure motion is derived from already established data for the cross-flow excitation coefficient. The proposed model is tested by numerical simulations, and the results are compared to experimental observations. When a sinusoidal cross-flow motion is given as input to the algorithm, the generated force time series are generally in good agreement with experimental measurements of cross-flow force in phase with cylinder velocity and acceleration. The model is also utilized in combination with time integration of the equation of motion to simulate the cross-flow vibration of a rigid cylinder. The resulting amplitude and frequency of motion as functions of reduced velocity are compared to published experimental results. In combination with the finite element method, the model is used to simulate cross-flow vibrations of a flexible cylinder in shear flow. Comparison with experiments shows that the model is capable of reproducing important quantities such as frequency, mode and amplitude, although some discrepancies are seen. This must be expected due to the complexity of the problem and the simple form of the present method.     

刚性球体三维高速垂直自由入水载荷

入水结构体在从空中弹道转入水下弹道的入水阶段，其周围的流体将呈现出强非线性性质，本文针对传统基于Wagner理论的结构体入水载荷计算模型不能很好描述流体三维流动的情况，基于无黏不可压流体流动模型，考虑流体弹性，采用微元边界运动等效方法对运动边界进行分段分析，计及入水过程中系统的动能损失，根据能量守恒，对刚性球体高速垂直自由入水过程中流体的三维流动进行了理论分析，建立了基于无黏不可压弹性流体的刚性球体垂直高速入水载荷计算模型，并与基于多介质任意拉格朗日欧拉方法的有限元模型进行了对比分析，验证了该方法的可行性。基于此模型，本文进一步分析了入水载荷的影响因素。该方法提供了一种计算结构体垂直高速入水载荷的思路，具有一定的理论意义和工程应用价值。     

Interaction of a shock with a sphere suspended in a vertical shock tube

Shock wave interaction with a sphere is one of the benchmark tests in shock dynamics. However, unlike wind tunnel experiments, unsteady drag force on a sphere installed in a shock tube have not been measured quantitatively. This paper presents an experimental and numerical study of the unsteady drag force acting on a 80 mm diameter sphere which was vertically suspended in a 300 mm x 300 mm vertical shock tube and loaded with a planar shock wave of M _s = 1.22 in air. The drag force history on the sphere was measured by an accelerometer installed in it. Accelerometer output signals were subjected to deconvolution data processing, producing a drag history comparable to that obtained by solving numerically the Navier-Stokes equations. A good agreement was obtained between the measured and computed drag force histories. In order to interpret the interaction of shock wave over the sphere, high speed video recordings and double exposure holographic interferometric observations were also conducted. It was found that the maximum drag force appeared not at the time instant when the shock arrived at the equator of the sphere, but at some earlier time before the transition of the reflected shock wave from regular to Mach reflection took place. A negative value of the drag force was observed, even though for a very short duration of time, when the Mach stem of the transmitted shock wave relfected and focused at the rear stagnation point of the sphere.Received: 31 March 2003, Accepted: 7 July 2003, Published online: 2 September 2003     

Fundamental solutions for axi-symmetric translational motion of a microstretch fluid

The fundamental solution for the axi-symmetrictranslational motion of a microstretch fluid due to a concentrated point body force is obtained.A general formula for thedrag force exerted by the fluid on an axi-symmetric rigid particle translating in it is then deduced.As an application to theobtained drag formula,this paper has discussed the problemof creeping translational motion of a rigid sphere in a microstretch fluid.The slip boundary condition on the surfaceof the spherical particle is applied.The drag force and theother physical quantities are obtained and represented graphically for various values of the micropolarity and slip parameters.     

Fluid-structure impact analysis with a mixed method of arbitrary Lagrangian-Eulerian BE and FE

A mixed method of arbitrary Lagrangian-Eulerian boundary element and finite element method (ALE-BE-FE method) is proposed for solving fluid-structure impact problems, in which the effect of structural deformation due to hydrodynamic pressure is taken into account. In addition, this method also enables us to analyze the influence of nonlinear free surface conditions on the impact response. Two numerical examples of an impacting cylinder and an impacting wedge into an initially calm water treated as 2-D problems are presented. It shows that the proposed method is effective to obtain a fluid-structure impact solution.This project is financially supported by the National Education Foundation of China.     

流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

研究柔性体撞击问题的子结构离散方法

本文用子结构方法研究刚性小球和均质柔性杆的纵向撞击以及和均质柔性梁的横向撞击问题,导出了用模态坐标表示的动力学方程,通过对刚性小球和柔性杆的纵向撞击的仿真教育处发现,在总单元数相同的情况下,取8个子结构比较合理,在此基础上对刚性小球和均质柔性梁的横向撞击的特性进行研究,发现撞击力在变化过程中会产生上下波动,当梁的弹性模量增加时,撞击力增大,撞击时间缩短。     

Accelerating motion of a sphere in a maxwell fluid

The translatory accelerating motion of a sphere due to an arbitrarily applied force in an unlimited Maxwell fluid is considered. The exact solutions for the velocity of the sphere for three particular types of accelerating motion are presented. The first is for a falling sphere; the second is for the decelerating motion of a sphere after the force which maintains the sphere with a constant velocity is removed; the third is for the motion of the sphere subjected to an impulsive force. The exact solutions are expressed in terms of real, regular, definite integrals which can be evaluated by numerical technique. Also presented are the asymptotic solutions for the velocity of the sphere in all three cases which are valid for small values of time.     

一种强耦合预估-校正浸入边界法

为克服传统浸入边界法的质量不守恒缺陷,提出了一种用于可压缩流固耦合问题的强耦合预估-校正浸入边界法。通过阐述一般流固耦合系统的矩阵表示,推导了流固耦合系统的强耦合Gauss-Seidel迭代格式,进一步导出预估-校正格式,提出了预估-校正浸入边界法。该方法使用无耦合边界模型对流体进行预估,将流固耦合边界视为自由面,固体原本占据的空间初始化为零质量的单元,允许流体自由穿过耦合边界。对于流体的计算,使用带有minmod限制器的二阶MUSCL有限体积格式和基于Zha-Bilgen分裂的AUSM+-up方法,配合三阶Runge-Kutta格式推进时间步。在校正步骤中,通过一组质量守恒的输运规则来实现输运过程。输运算法可概括为将边界内侧的流体进行标记,根据标记顺序以均匀方式分割和移动流体,产生一个指向边界外侧的流动,最后在边界附近施加速度校正保证无滑移条件。标记和输运算法避免了繁琐的对截断单元的几何处理,确保了算法易于实现。对于固体的计算,分别采用一阶差分格式和隐式动力学有限元格式求解刚体和线弹性体,并利用高斯积分获得固体表面的耦合力。使用预估-校正浸入边界法计算了一维问题和二维问题。在一维活塞问题中,获得了压力分布、相对质量历史和误差曲线,并与其他方法进行了对比。在二维的激波冲击平板问题中,获得了数值模拟纹影和平板结构的挠度历史,并与实验结果进行了对比。研究表明,该方法区别于传统的虚拟网格方法和截断单元方法,能够精确地维持流场的质量守恒并易于实现,且具有一阶收敛精度,能够较准确地预测激波绕射后的流场以及平板在激波作用下的挠度,为开发流固耦合算法提供了一种新的思路。