轴对称体垂直出水的非线性数值解

根据格林定理,本文采用严格的非线性边界条件,对任意形状的轴对称体的垂直出水,包括水下运动、接近液面、穿出液面的全过程进行了数值计算。计算结果与实验结果比较表明,非线性计算结果是可信的。     

自由液压流体流动问题的数值分析研究

对国内外具有自由液面的流体流动问题研究的数值求解方法进行了较为全面的评述。文中首先讨论了自由液面流体流动问题的数值求解方法，包括差分法、有限元法和边界元法等，然后说明了自由液面的跟踪和网格更新等问题。     

微重力作用下旋转轴对称容器内液体静表面形状分析

)引言随着液体火箭和充液卫星的不断发展,在失重或微重环境下,容器内液体晃动问题引起了普遍的重视。理正从理论与实验两方面进行研究。一般地说,在静平衡液面上提出边界条件比较方便;对于微重情况,自由液面已不 ...     

带弹性隔板的矩形贮箱中流弹耦合振动问题

研究了带弹性隔板的部分充液的矩形贮箱中液体晃动与弹性隔板振动的耦合问题．对防晃板的上下两部分中的液体分别采用不同的速度势函数,确定了双势函数满足的复杂边界条件和防晃板振动满足的边界条件,采用能量法导出了耦合频率的方程;通过数值计算得到了基频随隔板位置及刚度变化规律．     

基于SPH法的二维矩形液舱晃荡研究

液体晃荡是一种复杂的流体运动现象,自由液面的存在使得该现象具有很强的非线性和随机性。针对二维矩形液舱在不同振幅水平激励下的纵荡问题,应用SPH法对其进行了数值研究。首先计算了小振幅激励下的纵荡,计算结果分别与线性理论解、文献VOF法结果及文献SPH法结果作了比较分析,以验证所建数值模型的合理性;然后计算了液舱在大振幅水平激励下的纵荡,着重分析了不同振幅下液体晃荡的速度向量图、液面波动时程、压强波动时程、动量波动时程以及波动的频谱图,并将计算所得液面波动结果与小振幅激励下的液面波动结果作了比较。分析结果表明,在大振幅水平激励下,液面波动的波峰值较小振幅下的结果有较为明显的增大,而波谷值则无过大的变化,总体波动幅值比小振幅下的结果大;随着激励幅值的增大,液面波动幅值呈现明显增大的趋势,压强的整体波动幅值也呈增大趋势,动量波动的均值亦有明显增大;波动能量随着激励幅值的增大而增大并向第一阶频率区域集中。SPH法对处理液体大幅晃荡这种具有自由表面大变形的问题有十分优越的特性。     

贮腔类三维自由液面动力学问题数值研究

讨论了贮腔类三维自由液面动力学问题的数值研究,将任意的拉格朗日-欧拉运动学描述关系引入到系统的控制方程中,采用任意的拉格朗日-欧拉描述跟踪自由液面,推导了自由面上结点的法向矢量计算公式。采用Galerkin余量法推导了Navier-Stokes方程的空间离散有限元方程,采用三维自由液面上微分几何理论推导了表面张力计算公式。数值研究中考虑了接触角效应,最后进行了三维数值算例分析。     

圆筒形贮腔中微重力液体非线性晃动的数值模拟

本文讨论低重力液体在圆筒形贮腔中的非线性晃动问题。将ALE（任意的拉格朗日－欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkinlk加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式。推导了自由液面上法向矢量的计算公式。模拟了圆筒形贮腔中低重力液体的非线性晃动,并得到了自由液面、波高变化、压力响应等非线性动力特性。揭示了微重力液体非线性晃动的重要特征并将所得结论与现有的实验结果进行了比较。从而证实了本文方法的有效性与正确性。     

Bingham流体与固壁、液面斜撞击的数值模拟

在交错网格中用MAC方法求解二维不可压N S方程 ,对圆截面液柱与固壁、液面斜撞击问题进行了数值模拟 ,得到了自由面随时间演化的图像。主要考察了 :(1)固壁和液面 ;(2 )牛顿、非牛顿流体 ;(3)碰撞入射角 ;(4 )Bingham流体近似本构式中参数q0 、K对计算结果的影响。结果表明 :液柱与液面碰撞形成的自由面更复杂。碰撞初期 ,Bingham流体和水的自由面相似 ;但碰撞后期 ,Bingham流体的自由面相对简单。对液柱与液面碰撞自由面的影响较大。本文条件下 ,当K 2 0时 ,K对自由面的影响不大 ;当q0 增大时 ,自由面变得相对简单。     

平面炸药对刚性飞片的正向与反向驱动问题

本文采用特征线法和有限差分的数值计算方法,对平面炸药在不同端面起爆和不同边界条件下(正向与反向)驱动刚性飞片的问题进行了研究。所得结论是:为要提高炸药能量的利用效率,惰性炮塞的存在是个关键。在一定条件下,反向驱动比正向驱动可使飞片更快加速,以致在爆炸初期,在飞片运动的一定距离内,反向驱动可使飞片速度提高约15—20％。飞片的极限速度与Gurney理论进行了比较,二者符合一致。     

可压缩材料平面应变滑移线场的数值解和边界条件

简要介绍了可压缩材料平面应变滑移线理论。导出了滑移线场应力和速度数值解。处理了应力边界条件，阐明了应力和速度间断规律。举例说明了边界条件处理和数值解步骤。     

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier–Stokes–Boussinesq equations governing the transport of momentum, mass and heat in a non-isothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a non-orthogonal co-ordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.     

Basic flows of a micropolar liquid

Steady plane shear flow, rotation flow, plane wave disturbances, oscillatory boundary layers, and free surface waves in a micropolar liquid are investigated.     

The stability of the motion of a liquid,due to thermocapillary forces

A study is made of the stability of the plane-parallel flow of a viscous liquid in a layer with a free boundary, under weightless conditions. The motion of the liquid is due to the dependence of the surface tension on the temperature. An exact solution for an unperturbed boundary is obtained by the same method used in [1], but with a more general boundary condition for the temperature. A study of the stability was carried out by the method of small vibrations, taking account of the perturbation of the free boundary. The article discusses the asymptotic behavior of long waves at small Reynolds numbers, and the conditions for instability are found.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 94–98, November–December, 1971.     

半浮区液桥振荡对流的微重力实验

利用日本微重力中心８００ｍ落井装置，完成了半浮区液桥振荡对流的微重力实验，对振荡对流的典型物理量诸如内部温度、流场、自由面边缘变化及表面波进行了综合测量。实验结果给出了振荡对流由地球重力环境向微重力环境的过渡，以及不同几何参数半浮区液桥的振荡特征，并首次获得了微重力环境下热毛细对流的表面波位形及边缘振荡特征．     

Downward Flow of a Nonisothermal thin Liquid Film with Variable Viscosity

Equations of a thinfilm flow with linear dependences of viscosity and surface tension of the liquid on temperature are derived. The impact of various factors on the shape of the free boundary of the film are numerically analyzed.     

Numerical simulation of unsteady inviscid flows with free surfaces

The boundary element method is used to calculate numerically the unteady flow of a capillary liquid associated with the interaction of an expanding gas cavity and the free surface of the liquid.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–7, March–April, 1990.     

Circulation flow around airfoils by a steady plane-parallel flow of a heavy liquid of finite depth with a free surface

Steady problems of a circulation flow around bodies by a flow of a heavy liquid bounded by a free surface and a straight bottom are solved. The method of complex boundary elements is used, which is based on the integral Cauchy formula written for a complex-conjugate velocity. Results of numerical calculations of the flow around a circular contour and the Joukowski airfoil are presented. Shapes of the free surface and the most important hydrodynamic characteristics of the process (velocity circulation over the airfoil and the lifting force and its moment relative to the sharp edge of the airfoil) are given. Kemerovo State University, Kemerovo 650043. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 101–110, May–June, 2000.     

Free surface asymptotics for thermocapillary flow at high marangoni numbers

Formal asymptotic expansions of the solution of the steady-state problem of incompressible flow in an unbounded region under the influence of a given temperature gradient along the free boundary are constructed for high Marangoni numbers. In the boundary layer near the free surface the flow satisfies a system of nonlinear equations for which in the neighborhood of the critical point self-similar solutions are found. Outside the boundary layer the slow flow approximately satisfies the equations of an inviscid fluid. A free surface equation, which when the temperature gradient vanishes determines the equilibrium free surface of the capillary fluid, is obtained. The surface of a gas bubble contiguous with a rigid wall and the shape of the capillary meniscus in the presence of nonuniform heating of the free boundary are calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 61–67, May–June, 1989.     

不定边界问题的边界元和有限元数学模型

对于不可压缩有势流动,有两类典型的不定边界或可动边界问题,即定常型不定边界如常水位闸门出流、过水坝或过水闸、水堰绕流及射流等,非定常型的不定边界如装液容器内的流体晃动问题.它们的共同特点是自由面位置与形状事先是的,需在计算过程中调整网格作自由面拟合.且由于自由面条件呈强烈非线性,给不确定数值计算带来困难.本文综述了两类不定边界问题的有限元和边界元模式,简述了笔者的一些计算经验.