Thermodynamics and dynamics of surface phases

The paper presents and discusses the set of field-equations for an interface between two immiscible, uncharged, unpolarized, non-elastic and non micro-polar fluids modelled as a geometrical surface S in arbitrary motion and with thermodynamic and dynamic properties. The crucial assumptions are those implied by a local formulation of the equilibrium thermodynamics of S. The set of equations derived comprises: integral and local forms of the surface balance-equation for a scalar or vectorial extensive property; continuity, normal and tangential momentum and energy equations; surface entropy production; linear phenomenological relations for irreversible surface diffusions and non-equilibrium exchanges of extensive properties between surface and volume phases. Momentum and energy conservation equations are formulated also in terms of a suitably defined surface total enthalpy. Available previous formulations are either not as complete or hold only in more particular cases. Whenever applicable, their results are compared with the present ones and commented upon.     

基于有限元法的大气层外弹道式目标表面温度场分析

获取目标表面温度场是进行红外特征分析的重要前提。为确定大气层外弹道式目标的表面温度场分布,建立了有限元模型。根据目标温度场的轴对称分布特点,在柱坐标系内建立了二维瞬态热传递模型,从而降低了计算量。与当前大部分的有限元软件和文献不同,严格分析了温度非均匀性有限单元的辐射热损失,并推导出轴对称三结点三角形单元的辐射热损失公式。为了能够使用Galerkin法求解时间微分方程组,将辐射热损失视为与时间相关的热载荷项,并用Newton\|Raphson法迭代求解处理后得到的非线性方程组。一个简单的数值实验表明了所提出方法的有效性,最后应用该法求解了目标飞行全程的表面温度场分布。
     

Problems related to non-linear variations of surface tension

Generally, the tension at a liquid/gas interface is decreasing with temperature. For aqueous solution of alcohol, alkylammonium chloride in water, some binary metallic alloys and nematic liquid crystals the surface tension presents an extremum as a function of temperature. The influence of this property on the hydrodynamic stability and on the surface tension induced flows is studied as well theoretically as experimentally.     

An approach to the autoignition of a turbulent mixture

This paper considers the turbulent homogeneous mixing of two reactants undergoing a one step, second order, irreversible, exothermic chemical reaction with a rate constant of the Arrhenius type. A statistically stationary turbulent velocity field is assumed given and unaffected by mass or heat production due to the chemical reaction. Relative density fluctuations are neglected. A Hopf-like functional formalism is presented, with application to both statistically inhomogeneous and statistically homogeneous flows. Single and double point probability density function differential equations are derived from those functional equations. The limit of very large activation energies is considered; a low degree of statistical correlation between temperature and concentration fields during the ignition period is hypothesized. After making use of the homogeneity assumption a closure problem is still present due to the nonlocalness of the molecular diffusion term. The problem is rendered closed by assuming a Gaussian conditional expected value for the temperature at a point given the temperature at a neighboring point. The closure is seen to preserve very important mathematical and physical properties. A linear first order hyperbolic differential equation with variable coefficients for the probability density function of the temperature field is obtained. A second Damköhler number based on Taylor's microscale turns out to be an important controlling parameter. A numerical integration for different values of the second Damköhler number and the initial stochastic parameters is carried out. The mixture is seen to evolve towards an eventual thermal runaway, the detailed behavior however being different for different systems. Some peculiarities during the ignition period evolution are uncovered.     

Stability of large flexible damped spacecraft modeled as elastic continua

Vibrational stability of a large flexible, structurally damped spacecraft subject to large rigid body rotations is analysed modelling the system as an elastic continuum. Using solution of rigid body attitude motion under torque free conditions and modal analysis, the vibrational equations are reduced to ordinary differential equations with time-varying coefficients. Stability analysis is carried out using Floquet theory and Sonin-Polya theorem. The cases of spinning and non-spinning spacecraft idealized as a flexible beam plate undergoing simple structural vibration are analysed in detail. The critical damping required for stabilization is shown to be a function of the spacecraft's inertia ratio and the level of disturbance.     

自由浮动空间机械臂的最优控制

自由浮动空间机械臂所满足的动量矩守恒方程是一个非完整约束条件。首先把铰关节角速度作为输入，建立起系统的状态方程。然后把输入表示成待定系数的Ｆｏｕｒｉｅｒ 基组合，从而将最优控制转化为对代数问题的求解。最后用非线性最小二乘法迭代求出同时使载体和铰关节达到预定位形，而且使输入能量最小的最优控制。仿真算例说明，本文的方法是有效的。     

“希望一号”卫星热平衡试验的误差分析

为研究卫星热平衡试验误差,建立“希望一号”卫星热试验的联合热数学模型,对热平衡试验进行模拟,并与理想状态进行比较,发现仪器试验温度值在低温工况中偏低13～20 ℃,高温工况中偏低11～17 ℃;高温工况中-Y舱板表面到达热流不均匀度高达27.3%。此外,分析了热试验中热沉的温度、表面发射率、与卫星的表面积比,支架的导热、辐射漏热,电缆的导热、辐射漏热,以及热流计温度测量误差对卫星温度场的影响。     

Surface effects on the free vibration behavior of postbuckled circular higher-order shear deformable nanoplates including geometrical nonlinearity

This investigation deals with the free vibration characteristics of circular higher-order shear deformable nanoplates around the postbuckling configuration incorporating surface effects. Using the Gurtin–Murdoch elasticity theory, a size-dependent higher-order shear deformable plate model is developed which takes account all surface effects including surface elasticity, surface stress and surface density. Geometrical nonlinearity is considered based on the von Karman type nonlinear strain–displacement relationships. Also, in order to satisfy the balance conditions between bulk and surfaces of nanoplate, it is assumed that the normal stress is distributed cubically through the thickness of nanoplate. Hamilton?s principle is utilized to derive non-classical governing differential equations of motion and related boundary conditions. Afterwards, an efficient numerical methodology based on a generalized differential quadrature (GDQ) method is employed to solve numerically the problem so as to discretize the governing partial differential equations along various edge supports using Chebyshev–Gauss–Lobatto grid points and pseudo arc-length continuation technique. A comparison between the results of present non-classical model and those of the classical plate theory is conducted. It is demonstrated that in contrast to the prebuckling domain, for a specified value of axial load in the postbuckling domain, increasing the plate thickness leads to higher frequencies.     

热屏绝热型热流计可行性分析

热平衡试验是航天器热控设计的重要环节,外热流模拟和测量是热平衡试验的两个关键因素。两者的不准确将给热平衡试验带来不可预估的影响。因此热流测量技术是热平衡试验的关键技术之一。目前使用的热流计测量低热流密度时存在较大误差,为此研制一种高精度热流计——热屏绝热型热流计来测量热流。该热流计主要包括敏感面和热屏两部分。敏感面用来测量热流,热屏则为敏感面提供一个绝热环境。其关键技术是在热屏和敏感面充分隔热的基础上,采用跟踪控温的方法使热屏和敏感面的温度相同,从而使两者之间的热交换可以忽略。文章利用热分析计算分析了该热流计其敏感面和热屏的等温性能满足热流测量要求;并利用试验证明了跟踪控温技术可行,热屏和敏感面温度一致,从而说明了该热流计测量热流的可行性。     

On an analytical model of the aerothermodynamics of fluid–structure interaction

An analytical model of the aerothermodynamics of fluid–structure interaction through an ablative wall is developed, involving: (i) the solution to the Navier–Stokes equations in the fluid; (ii) the solution to the equations of thermoelasticity in the solid wall; (iii) the matching of the two solutions across the fluid–wall interface by continuity of normal and tangential stresses and jump of the heat flux equal to the heat of ablation. These three matching conditions determine the three constants of integration, besides the boundary conditions: (i) in the free stream far above the wall, the fluid velocity and temperature are given; (ii) in the backing material below the wall the temperature is given and displacement is assumed to be zero. In order to arrive at a simple analytical solution the simplest geometry is chosen: two-dimensional, with a flat fluid–wall interface, and all quantities in the fluid and structure depending only on distance from the wall. Extensions to unsteady and three-dimensional cases would be more complex. The present theory specifies: (i) the velocity, pressure and density in the fluid; (ii) the displacement vector, and strain and stress tensors in the structure; (iii) the temperature and heat flux in the fluid and structure, including discontinuities across the fluid–wall interface. The model is applied to examples of (a) hypersonic flight and (b) atmospheric re-entry.     

空天飞行器六自由度数学建模研究

研究了空天飞行器超声速和高超声速飞行条件下六自由度仿真模型,该模型包含了完整的六自由度动力学方程和运动方程。气动力和力矩系数是迎角、马赫数及控制舵面偏角的函数;发动机模型为吸气发动机和变推力火箭发动机的组合推进装置;飞行器的质心;惯性矩是飞行器质量的时变函数。所得结果可以用于未来高超声速飞行器或新一代单级入轨运载器轨迹优化、姿态控制等问题的概念设计和仿真研究。     

高超声速流动壁面催化复合气动加热特性

针对高超声速流动壁面催化特性,计算了不同壁面催化复合系数条件下的球锥驻点热环境。引入了经验证的数值求解Navier Stokes方程的方法,在不同壁温500K～2500K的条件下分别分析了O 2和N 2气体在壁面处的催化复合气动加热特性,得到如下结论：（1） 原子复合放热将提高近壁面温度梯度,改变近壁面组分分布;原子复合放热一部分加热飞行器形成组分扩散热流,一部分加热近壁气体提高近壁温度梯度。（2） 在壁面催化复合系数较小时,原子复合放热主要转化为组分扩散加热,对于不同壁面温度,壁面催化复合系数α<0.1时,单一气体反应组分扩散热流小于总热流的20%。     

Monopropellant droplet decomposition for large activation energies

The quasi-steady decomposition of a monopropellant droplet in a quiescent atmosphere is analyzed, in the limit of large activation energy, for a gas phase decomposition reaction of the Arrhenius type.An analytical relation is obtained for the pre-exponential rate constant, or the Damköhler number, as a function of the droplet vaporization rate. The curve giving the vaporization rate in terms of the Damköhler number has, for low values of the temperature at infinity, an S shaped form which exhibits ignition-extinction characteristics.With increasing values of the droplet vaporization rate above the pure vaporization value we pass from a nearly frozen regime, under which ignition conditions occur due to the large temperature sensitivity of the reaction rate, to a complete decomposition regime, in which practically all the fuel is decomposed in a thin reactive-diffusion zone, separated from the droplet surface by a transport region of lower temperature. For sufficiently large Damköhler numbers, the thin reaction zone is located close to the droplet surface, so that the flame structure becomes planar: a linear relation between the droplet radius and time is obtained in this case. The non planar effects are responsible for the extinction of the flame at low Damköhler numbers and ambient temperatures below the adiabatic flame temperature. For sufficiently low values of the ambient temperature, an intermediate regime exists in which the enhancement of the vaporization rate is due only to the fraction of the vaporized fuel decomposed at a thin reaction zone, not far from the droplet, where the temperature reaches its maximum value; the remaining fuel is decomposed very far from the droplet.     

大型空间电站在轨展开与组装动力学与控制

针对大型空间电站在轨展开与组装过程中复杂的动力学环境，将桁架类大型空间电站视为一个刚柔多体系统，采用自然坐标法对刚性构件建模，基于绝对节点坐标方法描述柔性桁架结构，编写了一套动力学仿真软件，实现了大型空间电站在轨展开与组装动力学的精确仿真。仿真结果表明，采用摆线运动插值函数作为控制方程有利于提高多级展开过程的稳定性；在结构组装过程中，减小组装速度，增大组装机构阻尼系数与劲度系数，有利于提高结构稳定性。最终集成了一套针对大型空间电站的动力学仿真平台搭建方法，为大型空间电站在轨展开与组装动力学预示提供了理论依据和方法指导。     

The heliospheric magnetic field and its relation to the temperature,density, and velocity of solar plasma: Experimental evidence

Using daily and hourly data on solar plasma parameters at the Ulysses spacecraft orbit and at 1 AU it is demonstrated that there is a simple relationship between plasma temperature and density with the heliospheric magnetic field (HMF). A mathematical expression connecting HMF with plasma temperature and density is suggested. Correlation coefficients and regression equations for measured and calculated magnetic fields are presented for the 1990–2009 period according to Ulysses spacecraft data and for 2003–2010 at 1 AU (OMNI database). The roles played by density, temperature, and high-speed solar wind streams in forming the magnetic-field peaks are demonstrated using hourly data of OMNI2 and Ulysses.     

航天器热平衡温度预测的粒子群算法

为解决航天器热平衡试验温度外推预测方法中存在的问题,提出了热平衡温度预测的粒子群优化算法。根据热平衡温度模型,建立待优化的目标函数以及热平衡温度粒子群算法的计算流程。依据某型号卫星热平衡试验的前期温度数据,对热平衡极限温度以及瞬时温度进行预测,并与实测数据进行了比对。计算结果与实测数据相比具有较好的一致性,所提出的方法可有效地应用于航天器热平衡试验温度预测。     

Numerical modelling of the interface conditions in coupled marangoni flows

The implementation in numerical codes of the coupling conditions between bulk phases adjacent to an interface exhibiting Marangoni effect is examined. The conditions are detailed in the case of a vorticity, stream-function, temperature formulation of the field equations. Numerical examples are reported for particular cases.     

利用热平衡试验稳态数据修正对接机构热模型

根据热平衡试验数据与计算结果的比较,分析两者存在差异的原因,用热平衡试验稳态工况数据修正对接机构热模型中的单元划分和热网络方程。研究了对热网络方程修正的方法,基于参数化和节点群概念改进参数化节点群整合方法,以明显减少需修正系数的数量,并与分级修正和局部修正综合,提高热网络修正方法的实用性。算例表明：通过热模型修正使温度计算值与试验值的偏差达到允许的要求,表明该热模型修正方法有效。     

Effects of variable transport properties on thermal Marangoni flows

The effects of temperature dependence of the diffusion coefficients on the solution fields of Thermal Marangoni Flows in micro-gravity is analyzed by using an unsteady—variable mesh—two dimensional—finite difference numerical code to solve the full Navier-Stokes equations.A typical cylindrical floating zone geometry is analyzed and the estimation of the deviations caused by considering a temperature-dependent variable viscosity is verified with respect to the most important field parameters (streamlines, isotherms, isobars, Nusselt heating curves, velocity and temperature profiles).Comparison with TEXUS experiment is finally performed.