Free oscillations and response of a viscous liquid in a rigid circular cylindrical tank

The natural damped frequencies and the response to translational excitation of a viscous liquid in a circular cylindrical container are obtained. The response of the liquid surface elevation above its equilibrium position, as well as the viscous liquid force in x-direction due to translational excitation in this direction have been numerically evaluated. In this analysis the side wall adhesive boundary conditions have been satisfied, while only the normal bottom condition has been observed. This makes the results of the analysis applicable to liquid height ratio h/a>1.     

Frictionless liquid sloshing in circular cylindrical container configurations

If the free liquid surface in a circular cylindrical container is covered with solid structural parts, the natural frequencies and the response of the liquid in the container is drastically changed. This is also expressed by a much lower flow behavior inside the container. The benefit of a partly covered free surface plane lies in the shifting of the natural propellant frequencies above and away from the control frequency of a space vehicle, and in the decrease of the sloshing masses participating in the dynamic motion of the system. The fundamental natural frequency of a non-viscous and incompressible liquid has been determined for increasing width of an annular ring baffle attached to the tank wall. With increasing structural coverage of the free surface the natural frequencies increase considerably. The response to translational excitation has also been numerically evaluated and shows the shifting of the resonances to higher values. The procedure presented here may without any difficulties also be applied to various arbitrary coverages of the free liquid surface. It may also be used as a test case for purely numerical methods, such as finite element methods and others.     

圆柱容器液体晃动问题的数值计算

本文将二维VOF方法推广到三维情形,研究了圆柱容器遭受水平强迫振荡时的液体晃动问题。在数值模型中,采用了改进的施主与受主单元体方法计算流体在单元体之间的交换。同时对处理任意容器边界的部分单元体方法进行了改进,使之满足边界面上的速度边界条件。应用改进后的VOF方法得到的结果与试验结果符合较好,能够成功地数模出由于动力不稳定性引起的圆柱容器遭受纵向水平振荡时液体表面的旋转运动响应。     

基于SPH方法的圆柱形贮箱液体晃动研究

在飞行器燃油供给方式的选择和飞行姿态控制调整中,必须充分考虑贮箱内燃油的晃动情况及其对贮箱冲击力的影响。然而液体晃动问题的高度非线性使得其研究极具挑战性。采用光滑质点流体动力学方法对圆柱形贮箱内液体的晃动进行建模仿真,将晃动液体的控制微分方程组离散到一系列流体粒子上,从而将不同物理量函数的连续积分转换成支持域内粒子的离散化求和形式,通过求解线性方程组得到数值解。在固壁边界上布置一组虚粒子,以对内部粒子产生排斥力来模拟固壁边界条件,同时采用狄利克雷边界条件方法模拟自由液面,实现对自由液面的追踪。将仿真的液面变化与试验结果作了对比,吻合较好。得到了贮箱壁面上不同监测点的压强变化以及贮箱与液体之间的相互作用力,并说明了结果的准确性。数值计算方法对液体晃动问题的研究提供一定参考价值。     

失重环境下圆柱形容器内稳态旋转液体的液面形状

1.引言 在航天高技术中,由于星体充液量的增大以及飞行器姿态指向精度的严格要求,从而液体晃动问题受到了国内外科技界的广泛重视。飞行器在外层空间正常飞行时,液体处于失重状态,此时,不能忽视表而张力的作用,静止自由液面呈特殊曲面,其形状直接影响着液体的动力学特性由于卫星在入轨后处于稳定的自旋状态,这时液面由液体表面     

凝固过程液相区内晶体形成的模拟研究

通过NH_4C1水溶液二维凝固过程的模拟实验,研究了凝固过程中液相区内的晶体形成条件。结果表明,随凝固条件和时间的不同,型壁上晶体的游离、表面凝固产生“晶雨”及枝晶熔断均可在液相区形成晶体。     

凹槽孔槽深和动量比对气膜冷却特性的影响

采用瞬态液晶全表面测量技术测量了凹槽孔的气膜冷却特性,研究了动量比和槽深对冷却效率的影响,并和圆柱形孔的测量结果实施了比较.结果表明:冷气在凹槽内横向传播,凹槽抑制了冷气射流分离,凹槽孔比圆柱形孔提供更多的二维气膜.槽深孔径比为0.75是最佳的.随着动量比的增加,冷却效率值逐渐趋于稳定,冷却效率仅有轻微增加.      

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

高超声速弯曲激波压缩侧压式进气道数值研究

采用壁面马赫数呈线性分布的曲面压缩系统改进参考侧压式进气道的顶板,得到弯曲激波压缩侧压式进气道,并与参考侧压式进气道进行了比较.数值研究结果表明:设计状态无黏时曲面压缩顶板壁面马赫数分布与给定的马赫数分布基本一致,并且有黏时其壁面压力分布也与二维曲面的基本相同;同参考侧压式进气道相比,顶板采用曲面压缩能够一定程度地改善壁面压力分布,使其末端压力梯度变化平缓;并且非设计状态下的性能也得到有效地改善,特别是来流马赫数为4时,其流量系数提高6.0%、达到0.799,喉道截面总压恢复系数提高1.9%;来流马赫数为5时,其流量系数提高5.2%、达到0.909,喉道截面总压恢复系数提高3.2%.随着攻角增大,该进气道流量捕获能力增强、隔离段出口截面流场畸变减小,但喉道截面总压恢复系数下降剧烈.      

Hydrodynamic model of film liquid evaporation on a cylindrical wall

A mathematical model of hydrodynamic processes at film evaporation of a liquid running off a cylindrical vertical wall under gravity is presented.     

斜面兜孔圆柱滚子轴承的高速防打滑特性

为解决航空发动机主轴圆柱滚子轴承的打滑问题,提出一种具有斜面兜孔结构的圆柱滚子轴承,并对其高速防打滑特性展开研究。利用Hertz接触理论及弹流润滑、流体动力润滑的经验公式,建立了滚子-内外套圈滚道、兜孔-滚子、保持架-外圈引导面的接触模型。运用多体运动学及牛顿-欧拉动力学理论,建立了轴承径向平面内的三自由度动力学模型,在此基础上,利用龙格库塔数值积分法进行了轴承的动力学数值仿真,探讨了兜孔前后壁倾角对轴承打滑的影响规律,分析了兜孔前壁倾角为5°、后壁倾角为10°下轴承的高速打滑特性,结果表明：兜孔的前后壁倾角对轴承打滑有显著影响,通过对其进行优化,可有效抑制轴承的高速打滑。     

离心泵双密封环平衡系统泄漏特性预测方法研究

为了准确预测离心泵双密封环平衡系统的泄漏特性,通过分析该平衡系统的几何特征和液体流动特点,对其主要过流部件进行简化,将泵腔液体流动简化为由圆周剪切流和径向压差流叠加而成的轴对称二维粘性层流运动,同时将密封环间隙液体流动简化为两平行平板间的定常不可压缩粘性流体层流运动.然后,利用N-S方程分别建立了泵腔和密封环间隙的液体...     

通道曲面的柱面逼近方法及其在叶轮插铣中的应用

为了使叶轮粗加工后的通道表面与其设计曲面尽可能接近,将整体叶轮的定轴插铣理论模型概括为通道曲面的最大内接柱面优化问题。同时,为简化该问题的计算难度,将最大内接柱面的优化问题转化为柱面直母线方向和柱面准线的确定两个步骤。提出了确定内接柱面直母线方向的最小二乘模型与解算方法,在此基础上,将通道曲面沿柱面的直母线方向进行投影,通过相关曲线间的裁剪操作,获得了最大内接柱面的准线。以最大内接柱面的准线作为插铣加工时的边界曲线,以柱面直母线方向作为插铣加工时平头刀的轴线方向,给出了计算刀心位置和安全高度的迭代算法。最后通过数值算例验证了所提方法的正确性。该方法将有助于叶轮粗加工效率的提高,同时为叶轮的精加工预留较均匀的加工余量。     

定常三维分离准则探讨

本文从奇点理论拓扑规律出发,说明在垂直于分离线及物面的截面上的截面流线有半鞍点存在,进而说明Lighthill等人认为的极限流线要在分离线以前抬起来是不合适的。合理的说法应该是极限流线总是同剪切应力线一致。只有在进入奇点以后,极限流线才抬起来离开物面。本文最后表明,单纯把二维分离准则用到三维分离的情况是不够的,并给出三维分离准则的表达式。     

剪切气流驱动液滴运动的受力分析

实验研究了剪切流驱动的液滴在固体表面上起始运动的受力机理。工作中使用一系列液体和固体表面来获得不同的液滴接触角,并在小型风洞中进行实验。实验中对液滴的启动气流速度进行了测量,并综合各种起始时刻的参数信息,建立了一个关于液滴接触线表面张力和剪切气流拖拽力平衡的数学模型,揭示了液滴脱落时刻的受力情况。所建立的模型更适合液滴1变形情况,但对于其它类似情况的剪切气流驱动液滴运动也能够进行合理的描述。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

60m3卧式液氮贮罐预冷过程试验

为进一步了解大型卧罐预冷过程，搭建了60 m3卧式液氮贮罐预冷试验台，对小流量下贮罐预冷过程的降温和罐体应变特性开展了研究。结果表明：预冷过程罐内气体温度首先整体迅速降低，然后缓慢下降，且呈现分层现象；预冷初期液氮在罐底难以积累；靠近底部的罐壁降温过程分为三个阶段，首先与低温氮气进行自然对流传热，温度缓慢下降，液位增长到相应高度后与液氮进行沸腾换热从而温度迅速下降，最后稳定在液氮温度；对于最终液位以上的罐壁，一直维持着平稳的降温速率；贮罐轴向应变全为负值，随着预冷过程进行轴向应变随之增大，与液氮接触的局部罐壁轴向应变会迅速增加。该项试验的成功进行有力地补充了国内大型卧罐预冷试验数据的空白，为低温贮罐可靠性及寿命预测等相关研究提供数据支撑。     

大型复合材料箱体成型模具设计

文 摘 为了制造合格的大型复合材料箱体，本文制定了箱体成型模具的设计和制造方案。某复合材料箱体属于大型长轴类产品，其成型模具挠曲变形及脱模力大。模具主体采用通长轴、周向辐板、轴向筋板和蒙皮组合的结构形式。模具通过无缝钢管、模具骨架及蒙皮形成的三维网状结构有效保证模具刚度，通过带斜度工作平面及变圆角设计大大降低脱模力，通过在顶出端使用更厚的钢板提升其在脱模时的可靠性。仿真结果表明，简支工况下模具最大应力值为42.9MPa，最大变形量为0.51mm；脱模工况下模具最大应力值为188.6MPa，最大变形量为0.5mm。根据模具结构制定了分别组焊、整体装配的加工工艺路线并完成了模具制造。仿真及模具实际使用结果表明，该模具满足复合材料箱体的生产需求。     

Tracing the magnetic topology of coronal mass ejection events by Ulysses/HI-SCALE energetic particle observations in and out of the ecliptic

In January 2000, the Ulysses spacecraft observed an ICME event at 43° S heliographic latitude and ∼ 4.1 AU. We use electron (E _e>38 keV) observations to trace the topology of the IMF embedded within the ICME. The still controversial issue of whether ICMEs have been detached from the solar corona or are still magnetically anchored to it when they arrive at the spacecraft is tackled. An in ecliptic ICME event is also presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

Accurate measurements of local skin friction coefficient using hot-wire anemometry

The practicality and accuracy of many existing methods of local skin friction measurement suffer when the boundary layer flow under consideration is non-canonical. Such shortcomings are exacerbated in three-dimensional flows, by the necessity to map local c_f over a wider area in order to characterise fully the contribution to global skin friction. These problems have led the authors to seek novel experimental methods of c_f measurement. The technique proposed herein utilises velocity measurements made using hot-wire anemometry combined with accurate positioning of the sensor element in respect to the test surface. In essence it is proposed that the local skin friction can be evaluated via a single velocity measurement made at a known wall-normal distance within the linear region of the viscous sublayer. This technique relies on accurate probe positioning, and two methods of achieving this are outlined. A study of the hot-wire characteristics in near-wall proximity has revealed a previously unnoticed feature corresponding to probe–wall contact. It is shown that this anomaly can be used as a positional flag to accurately locate the aerodynamic origin of the hot-wire sensor. A second technique using a laser triangulation displacement sensor is also outlined. Both positional techniques are shown to offer positioning to a sufficient level of accuracy for the proposed c_f measurement technique. Single-point local c_f measurement is tested experimentally, demonstrating the improved repeatability and standard error as predicted by initial error analysis. In this way it is shown that a single $\text{90}\text{s}$ velocity sample coupled with accurate wall positioning can define local c_f to a standard error of σ_cf≈1.0%. Analysis of error contributions reveals that longer sampling periods can realise even greater accuracy. The proposed technique is also used to measure local c_f in a three-dimensional boundary layer where micro-vortex generators have introduced large-scale spanwise distortions. This is an example of an application in a non-canonical boundary layer, and initial results show that the method is capable of providing repeatable comparative spanwise-averaged skin friction results to within approximately ±1%. The technique is immediately applicable to researchers using hot-wire anemometry in low Reynolds number flow over electrically non-conductive test surfaces.