The loss of stability in ground based experiments in liquid bridges

In this study our primary goal is to investigate the loss of stability of the steady convection in deformable liquid bridges. Definitely, the deformation plays an important role in the transition process from the steady axi-symmetric 2-D basic state to the 3-D periodical one. As it was shown experimentally by many researchers, the critical Marangoni number is very sensitive to the volume of liquid.

It seems to be proved^1–4, that the critical wave number for high Prandtl fluids (for example, silicone oil of different viscosities. Pr > 50) is m = 1 for aspect ratios Γ = height/radius ≥ 1.0 regardless of the free surface shape. But the data from our last experiments show, that by changing surrounding conditions around the liquid bridge we can change the critical wave number. Particularly, by placing the liquid bridge of diameter 2r = 6mm into another cylindrical volume of the diameter 2R = 12mm kept at constant temperature, the critical mode is switched from m = 1 to m = 2.     

Effect of liquid bridge form on oscillatory thermocapillary convection under 1 g and μg conditions

This study analyses the effect of temperature difference between hot and cool disk (ΔT), and non-dimensional liquid bridge volume (V/V_o) on the transition process from steady thermocapillary convection to periodic or chaotic thermocapillary convection in a liquid bridge modeled after the floating zone method under normal gravity and microgravity conditions. From normal gravity and drop shaft experiments, the difference of the regime of the steady state and the oscillatory state was clarified on the ΔT−V/V_o plane under 1 g and μg conditions. A gap or stability region was observed in the specific V/V_o range under 1 g conditions. In the gap or stable region, after the gravity changed from 1 g to μg conditions, the temperature signals showed oscillation. From these results, the critical temperature difference under the μg conditions appeared to be smaller than that under the 1 g conditions. Temperature signals were defined as 6 different types of states. The various temperature oscillatory state regimes were obtained on a ΔT−V/V_o plane under 1 g and μg conditions. Under μg conditions, in these experimental conditions, all temperature oscillatory states exhibited only the Periodic state.     

Oscillatory Thermocapillary Flows in Simulated Floating-zones with Time-dependent Temperature Boundary Conditions

This study deals with numerical simulations of the Maxus sounding rocket experiment on oscillatory Marangoni convection in liquid bridges. The problem is investigated through direct numerical solution of the non-linear, time-dependent, three-dimensional Navier–Stokes equations. In particular, a liquid bridge of silicon oil 2[cs] with a length L=20 [mm] and a diameter D=20 (mm) is considered. A temperature difference ΔT=30 [K] is imposed between the supporting disks, by heating the top disk and cooling the bottom one with different rates of ramping. The results show that the oscillatory flow starts as an ‘axially running wave', but after a transient time the instability is described by the dynamic model of a ‘standing wave', with an azimuthal spatial distribution corresponding to m=1 (where m is the critical wave number). After the transition, the disturbances become larger and the azimuthal velocity plays a more important role and the oscillatory field is characterized by a travelling wave. The characteristic times for the onset of the different flow regimes are computed for different rates of ramping.     

Oscillatory thermocapillary flows in simulated floating-zones with time-dependent temperature boundary conditions

This study deals with numerical simulations of the Maxus sounding rocket experiment on oscillatory Marangoni convection in liquid bridges. The problem is investigated through direct numerical solution of the non-linear, time-dependent, three-dimensional Navier-Stokes equations. In particular, a liquid bridge of silicon oil 2[cs] with a length L = 20 [mm] and a diameter D = 20 (mm) is considered. A temperature difference ΔT = 30 [K] is imposed between the supporting disks, by heating the top disk and cooling the bottom one with different rates of ramping. The results show that the oscillatory flow starts as an ‘axially running wave’, but after a transient time the instability is described by the dynamic model of a ‘standing wave’, with an azimuthal spatial distribution corresponding to m = 1 (where m is the critical wave number). After the transition, the disturbances become larger and the azimuthal velocity plays a more important role and the oscillatory field is characterized by a travelling wave. The characteristic times for the onset of the different flow regimes are computed for different rates of ramping.     

固体发动机人工脱粘的二维实验与数值模拟

针对固体火箭发动机前封头人工脱粘缝隙,设计了二维实验装置,用来模拟发动机点火过程中前封头人工脱粘位置受到点火冲击后的应力-应变情况.同时,采用流固耦合的数值方法,对比计算了实验工况,计算得到人工脱粘根部的应力-应变与实验吻合很好.通过数值研究发现,点火初期燃气对人工脱粘部位的冲击,会引起缝隙内部的压强振荡及装药表面的振...     

The effect of vibrations on heterogeneous fluids: Some studies in weightlessness

Under weightlessness, the effect of harmonic vibrations can easily induce average motions in fluids with density inhomogeneities. We will consider more particularly pure fluids near their gas–liquid critical point, where the temperature variation of all important parameters (e.g. interfacial tension, density difference) follows universal, scaling laws. We will thus study, below the critical point, the ordering of the gas–liquid interfaces (at equilibrium and during the phase separation) and, above the critical point, some thermovibrational instabilities. These studies have been performed in the Russian MIR station and in the MiniTexus and Maxus rockets of ESA and under magnetic compensation of gravity on earth.     

The effect of vibrations on heterogeneous fluids: Some studies in weightlessness

Under weightlessness, the effect of harmonic vibrations can easily induce average motions in fluids with density inhomogeneities. We will consider more particularly pure fluids near their gas–liquid critical point, where the temperature variation of all important parameters (e.g. interfacial tension, density difference) follows universal, scaling laws. We will thus study, below the critical point, the ordering of the gas–liquid interfaces (at equilibrium and during the phase separation) and, above the critical point, some thermovibrational instabilities. These studies have been performed in the Russian MIR station and in the MiniTexus and Maxus rockets of ESA and under magnetic compensation of gravity on earth.     

Vibration and pulsation processes in feed systems of liquid rocket engines

Pulsation and vibration process in high speed centrifugal pumps and feed lines were investigated. Linear and quadratic relations between pulsation parameters and main factors were obtained. The investigation of turbopump low frequency oscillations showed that the intensity of pulsations is proportional to the pump specific linear size. Linear mathematic models with lumped parameters fairly well demonstrate the vibration influence on low frequency flow pressure vibrations in pumps and feed lines. It is shown that, in some conditions, the pressure oscillation amplitude distribution by the vibrating feed lines length may be characterized by the presence of a node in an intermediate section. A one-dimensional mathematic model based on wave equations proved to be suitable for calculations of vibration and pulsation parameters within a wide frequency range.     

液氧贮箱增压过程研究

采用数值模拟方法对液氧贮箱增压过程进行研究.贮箱内流场采用流体体积函数（VOF）多相流模型考虑,选择标准双方程k-ε湍流模型分析湍流效应,气液两项之间的热量、质量转移通过自定义程序（UDF）求解.获得了贮箱压力、排液流量、气垫温度、液氧温度对贮箱内流场温度分布的影响.计算结果表明,在稳定增压过程中,贮箱液面无扰动,贮箱内温度分层分布;各参数变化时,对贮箱内温度分布的影响主要是温度梯度的变化,并且各工况下液面附近和扩散器附近温度梯度基本相同.     

液体火箭发动机涡轮泵转子弯扭耦合振动研究

为获得转子振动特性,针对液体火箭发动机涡轮泵转子系统建立了其在密封流体激振作用下的弯扭耦合动力学模型。通过数值仿真和试验研究了涡轮泵转子系统弯扭耦合振动的动力学特性,结果显示在密封流体激励作用下弯扭耦合振动的非线性特性显著。还研究了偏心距对涡轮泵转子系统弯扭耦合振动的影响。本研究可为液体火箭发动机涡轮泵转子的结构设计、诊断与维护提供可靠信息。     

The Impact of Variable Accelerations on Crystal Growth onboard Spacecraft by the Floating Zone Method

The numerical investigation of the impact of time-dependent accelerations (vibrations) on the flow and heat and mass transfer in the melt is carried out for the case of modeling the crystal growth by the floating zone method under conditions of microgravity that exist onboard spacecraft. The effects of the Archimedean buoyancy force and of vibrations of the free surface of fluid are considered separately. When solving the problem of the effect of the vibrations of the free surface of fluid, the previously obtained data were used. It is shown that vibrations of the free surface have a much stronger effect on the processes under consideration than the buoyancy. Some problems that are related to the newly discovered effects are discussed. The use of vibroprotected systems and a rotating magnetic field can help solve these problems. We plan to continue our investigations in future spacecraft experiments, in particular, at the International Space Station, which is under construction at the moment.     

航天器再入全过程轴对称烧蚀热防护数值仿真研究

对航天器再入全过程轴对称烧蚀热防护进行了全过程数值仿真研究.采用修正Lees驻点热流密度方法和参考焓方法计算再入热流密度.采用JANAF模型计算烧蚀率.利用有限元法计算钝锥体再入航天器烧蚀层在移动边界条件下的轴对称温度场.采用碳化层-热解面-原始材料的轴对称碳化烧蚀模型;推导了热解气体流量计算方法.针对再入飞行大热流密度条件下,用有限元方法求解瞬态温度场时会产生的时间和空间上解的振荡问题.通过分析温度振荡现象产生的原因,采用集中热容矩阵向后差分方法解决振荡问题.计算结果表明,在时间步长选择合适的情况下,求解集中热容矩阵能够很好地解决数值振荡问题,同时烧蚀率和温度场计算比较准确.     

液体大幅晃动情形的航天器刚液耦合动响应仿真分析

基于任意拉格朗日-欧拉有限元(ALE)方法处理液体大幅晃动问题,并系统地建立了航天器刚-液耦合动力学的数值模型.通过液体受惯性力、刚体受流体作用力与力矩来实现系统的耦合,采用交错算法对液体晃动模块与刚体动力学模块进行迭代求解.将刚-液耦合系统的数值计算结果与物理试验及液体刚化计算结果进行了比较,验证了本文方法的适用性和...     

小卫星瞬态外热流下动态传热特性分析

为揭示小卫星瞬态外热流下的动态传热特性规律,以小卫星双层集总参数模型为研究对象,推导得到动态热平衡方程。类比阻尼振荡系统,采用时频变换和传递函数分析的新思路,对温度与热流波动量间的幅值特性和相位特性变化规律进行理论研究,并利用数值方法进行验证。结果表明:推导获得了小卫星传热系统自然频率和阻尼比的热参数准则式,并证明了小卫星传热系统在热激励下振荡特性为过阻尼。阻尼比和频率比的增大及热流静位移的减小均可降低温度的波动幅度,不同频率比范围下阻尼比对热流的波动量与温度波动量间的相位差的影响呈相反规律,数值结果与解析分析结果一致。可为低热惯性小卫星的热控优化设计提供理论参考。     

超声电喷推力器的驻波特性

为研究超声电喷推力器(UAET)驻波尺寸特性的相关变化规律,建立发射极液面的液体振动数值模型,并开展驻波尺寸测量试验对数值模型进行验证。对比驻波间距的测量值与计算值发现,两者在变化趋势上可认为一致,且计算误差在6.55%以内。在此基础上,利用数值模型对不同波源频率、振幅下的驻波相关参数进行计算。实验结果表明,随波源频率升高,驻波间距、高度以及半径均下降,其影响机制在于频率主要对波纹形成时间产生影响;而振幅的升高对间距几乎不产生影响,仅会导致驻波高度和半径升高,其影响机制在于振幅对波纹形成应力产生影响。     

Theoretical study about the Marangoni convection in a liquid column in zero-gravity

The thermal Marangoni effect on the surface of a liquid bridge induces a convection inside the liquid. For an imposed arbitrary periodic axial circumferential temperature distribution on the liquid surface the velocity distributions in radial-, angular- and axial direction are determined theoretically by solving the linearized Navier-Stokes equations. Of particular interest is the effect of the viscosity parameter $\text{v}\text{a}{}^2$ and axial wave length to diameter ratio $\text{l}\text{a}$. It was found that the increase of viscosity decreases the magnitude of the velocity distributions and that for small axial wave length to diameter ratios the radial- and axial velocities exhibit peak values close to the free surface of the liquid. This is in a less pronounced way also true for the angular velocity, which shows for increasing moderate values $\text{l}\text{a}\text{(0 ≤}\text{l}\text{a}\text{≤ 2}$) a strong increase in magnitude and for larger axial wavelength a decrease again. For increasing axial wavelength the peak value of the radial- and axial velocity shifts towards the center of the liquid bridge, of which for a further increase a decrease of the magnitude appears.     

大功率、高转速、高扬程涡轮泵振动分析与减振研究

比功率和能量密度都很高的火箭发动机涡轮泵,工作时引起振动的基本原因之一是较大的损失功作用在质量较小的产品上,减振必须是以提高组件的效率降低其振动比载荷为首选措施,同时改变其刚度与强度。组件中动、静件叶片之间的能量转换是引起流体压力脉动的主要原因,避免压力脉动的频率与转子转动的倍频耦合,特别是较低的倍频,是降低产品耦合振动的关键。流体密封间隙内的激振力对产品振动的影响很复杂,但激振力的主频与转子的某一固有频率接近时,将会对转子激起很大的振动,改变密封间隙内流体激振力的频率是抑制流—固耦合振动的主要方向。     

超声速进气道喘振的机理研究

应用数值模拟方法对中心锥中心进气混压式进气道的喘振现象进行了研究。在数值计算的基础上,根据进气道出口截面每个网格点的压力、密度、速度等参数计算了进气道喘振过程中流量系数和总压恢复系数随时间的变化情况。同时给出了在喘振过程中激波振荡的振幅、频率、对应的波系图案。并根据进气道头部分离涡的发展情况以及进气道内通道中状态参数的变化情况对喘振产生的机理进行了分析,认为进气道头部分离涡对喘振的产生起关键的作用。     

Lift force acting on the cylinder in viscous liquid under vibration

The averaged lift force acting on a solid in viscous liquid in case of translational or rotary vibrations of a cavity is experimentally investigated. Experiments are performed with a heavy cylinder; different types of vibrations, translational and rotational, are investigated. It is found that the vibrations excite the mean lift force which could provide the suspension of solid even in gravity field. The repulsion lift force acts on the body near the walls of the vibrating cavity. It is caused by the viscous hydrodynamic interaction of oscillating body with the wall and is significant at a distance comparable with the thickness of Stokes layer. The intensification of vibration results in the excitation of tangential lift force, which is caused by the brake of the symmetry of the body's oscillations with respect to the cavity wall. In case of rotary vibrations the lift force of high intensity manifests itself in the bulk of the cavity due to the interaction of body with the oscillating shear flow. The mean dynamics of the solid body in a cavity under the rotary vibrations is determined by the combined action of two averaged vibrational effects—levitation of the body in the oscillating shear flow and hydrodynamic interaction of the body with the wall. In case of translational vibrations the dynamics of the solid is mainly determined by interaction with the walls. The experiments demonstrate that the vibrations have strong mean effect on the bodies in liquid; they could be used for efficient control of solid inclusions in microgravity and must be taken into account in space experiments and technologies.