高浓度悬浮固粒对射流界面的粘性稳定性影响

该文首次利用双流体模型和扰动速度势理论,推得含高浓度悬浮固粒的射流界面粘性稳定性方程和对应的固气扰动速度比值方程．通过数值计算,得到了不同雷诺数及固粒属性的射流界面粘性稳定性曲线和对应的固气扰动速度比值曲线．在分析和比较所得的粘性稳定性曲线的基础上,得到了流场雷诺数及固粒特性对射流界面粘性稳定性影响的结论．同时,通过分析所得的固气扰动速度比值曲线,得到了流场雷诺数及固粒等效斯托克斯数对固粒跟随气流的扰动 性能的影响的结论．这些结论是首次在计入气流的粘性的条件下得到的,不同于文献［８］和文献［１０］相关的囿于无粘情形的研究,对于两相射流发展的认识和工程实际中实施对两相射流场的人工控制有重要意义．     

含高浓度悬浮固粒运动射流稳定性研究

首先用连续介质耦合模型推导出含高浓度悬浮固粒运动射流的空间稳定性方程,然后借助渐进分析法和欧拉守恒差分格式,用有限差分数值解法得到不同流向位置、流场雷诺数、固粒属性和喷射装置运动速度时流场的稳定性特征曲线,说明喷射装置的反向运动使不稳定扰动频率范围扩大,正向运动则相反。固粒抑制流场的不稳定性,随着固粒等效斯托克斯数的减小,这种效应增强。这些结论,对于两相运动射流发展的认识有重要意义。     

悬浮固粒对二维混合层流动失稳特性的影响*

本文在不可压缩二维混合层流动方程的基础之上,通过添加固粒的作用项,推导得到了修正的瑞利方程;然后用数值计算方法解其特征方程,得到了悬浮固粒的质量密度、固粒和气流的速度比值以及Stokes数不同时二维混合层流动中扰动频率与空间增长率的关系曲线,给出了关于悬浮固粒对流场失稳特性影响的几个重要结论。     

气固两相固射流场涡结构影响固粒扩散的研究

用三维离散涡丝方法模拟轴对称圆射流场涡结构的发展,所得结果在某个周向位置上与用二维点涡方法计算的结果符合较好。然后采用单向淹合模型模拟固粒在圆射流中的运动,说明当固粒Ｓｔ数远小于１时,固粒受流场的作用较明显,当Ｓｔ数为１时,固粒主要分布在涡结构的周围,分布较均匀;当Ｓｔ数远大于１时,固粒受流场的影响较弱,当对涡环沿周向施以五个波长的扰动时,固粒扩散的范围较宽,固粒的扩散与扰动的振幅成正比,中所得     

气固两相圆射流场涡结构影响固粒扩散的研究

用三维离散涡丝方法模拟轴对称圆射流场涡结构的发展·所得结果在某个周向位置上与用二维点涡方法计算的结果符合较好·然后采用单向耦合模型模拟固粒在圆射流中的运动,说明当固粒Ｓｔ数远小于１时,固粒受流场的作用较明显,当Ｓｔ数为１时,固粒主要分布在涡结构的周围,分布较均匀;当Ｓｔ数远大于１时,固粒受流场的影响较弱·当对涡环沿周向施以五个波长的扰动时,固粒扩散的范围较宽·固粒的扩散与扰动的振幅成正比·文中所得结论与一些实验结果相符,对实际应用有指导意义     

两相流中柱状固粒对流体湍动特性影响的研究

对含柱状固粒的两相流场,建立了包含柱状固粒对流场影响的流体脉动速度方程,在求解脉动速度方程的基础上,经平均得到流体的湍流强度和雷诺应力.将该方法用于槽流湍流场的求解,并与单相流实验结果进行了比较.计算中变化柱状固粒的参数,给出了固粒的体积分数、长径比、松驰时间对流场湍动特性的影响,说明粒子对流场的湍动特性起着抑制作用,其抑制的程度与粒子的体积分数、长径比成正比,与粒子的松弛时间成反比.     

二维气固两相混合层中固粒对流场影响的研究

采用双向耦合模型对含有固粒的二维气固两相混合层流场进行了研究。流场采用拟谱方法直接数值模拟,固粒采用颗粒-轨道模型,在考虑流场对固粒作用的同时,考虑固粒对流场的反作用。结果发现固粒的浓度和Stokes数对流场影响明显。固粒的作用使涡量扩散加快,并阻碍流场的变化,减弱了流场中拟序结构的强度,缩短涡的生存期;固粒在流场中的分布规律与单相耦合所得结果相似。     

经过修正的层次流流动的流动稳定性问题(Ⅱ)——经过修正的平行剪切流的线性稳定性研究

本文根据文[1]给出的经过修正的层流流动的流动稳定性理论及平行剪切流中平均速度的一类修正剖面,研究了平行剪切流的线性稳定性性质,对于平面Couette流动和圆管Poiseuill流动,首次得到了二维扰动和轴对称扰动也能造成失稳的结果,并给出了这两种流动在某种定义下的中性曲线.     

柱状固粒两相边界层流场稳定性研究

从运动方程和本构方程出发,推导得到了含柱状粒子两相流场的修正Orr-Sommerfeld方程,然后在边界层流场中,采用数值计算方法,得到了含柱状粒子流场的稳定性中性曲线,给出了流场失稳的临界雷诺数.结果表明在所述情况下,柱状粒子对流场起着抑制失稳的作用,而且抑制的程度随着柱状粒子体积分数和长径比的增加而提高.     

气固两相流动系统基本特征的大规模模拟

运用基于颗粒运动分解思想的硬球模型模拟了气固两相流动系统的基本特征，讨论了操作条件对其的影响，其模拟结果表明：与国际上现有 用型模型相比，基于颗粒运动分解思想的硬球模型具有更为真实的优点，它能够定性而形象地复现气固流化系统的实验特征。     

纤维悬浮槽流空间模式稳定性分析

采用扰动的空间发展模式而非通常的时间发展模式,对含有悬浮纤维的槽流进行了线性稳定性分析。建立了适用于纤维悬浮流的稳定性方程并针对较大范围的流动Re数及扰动波角频率进行了数值求解。计算结果表明,纤维轴向抗拉伸力与流体惯性力之比H可以反映纤维对流动稳定性的影响。H增大使临界Re数升高,对应的扰动波数减小,扰动空间衰减率增加,扰动速度幅值的峰值降低,不稳定扰动区域缩小,长波扰动所受影响相对较大。纤维的存在抑制了流场的失稳。     

流体饱和多孔介质中波传播问题的有限元分析

采用基于混合物理论的多孔介质模型,给出流体饱和两相多孔介质波动问题的有限元分析方法·　采用罚方法导出的有限元动力方程,时间积分可采用显式和隐式积分两种方案·　用编制的有限元程序分析了一维柱体在跃阶载荷和脉冲载荷作用下的波传播问题,得到该两种瞬态载荷作用下固体和流体相位移、速度以及固体相有效应力和孔隙压力随时间的变化关系,并对波的传播现象进行了分析·　所得结果与理论相吻合·     

慢扩张旋转流的稳定性分析（Ⅰ）—理论研究

本文研究了无粘不可压慢扩张旋转流的稳定性问题，采用多重尺度展开法对有慢扩张的旋转流的非对称扰动进行浅化稳定性研究，导出了零阶及一阶扰动模所应满足的微分方程及由于慢扩张引起振幅变化的控制方程，将Ｐｌａｓｃｈｋｏ关于慢扩张喷流的结果推广到具有慢扩张的旋转流情况。     

The limiting properties of piecewise-potential subcritical and critical jet streams of an ideal gas

The limiting properties of subcritical and critical (with Mach numbers M ⩽ 1) plane-parallel jet streams are investigated in the approximation of an ideal (inviscid and non-heat-conducting) gas. Chaplygin's equation is used with the pressure and the angle of inclination of the velocity as the independent variables, which are measured from the limiting values corresponding to the cross-section of the equalizing of the jet with respect to these variables. The stream function in the neighbourhood of the “equalizing cross-section” is represented in the form of an expansion in powers of the “distance to the origin of the coordinates” (in the plane of the independent variables) with coefficients which depend on a previously unknown combination of the independent variables. The limiting property of the flow, i.e. the position of the equalizing cross-section at a finite of infinite distance, is defined by the leading terms of the expansion. A symmetric potential jet and symmetric piecewise-potential jets flowing out into a submerged space in the case of subcritical and critical pressure drops are considered as examples. The critical pressures of the potential parts of the composite jet can be different or identical (including when the thermodynamic properties of the gases are different).     

Electrohydrodynamic instability of two superposed viscous miscible streaming fluids

The linear stability of two dielectric viscous fluids separated by a horizontal interface is investigated. The interface admits heat and mass transfer. The system is stressed by a normal periodic electric field producing surface charges at the interface. The effect of surface tension, small viscosity, velocity streaming and gravity on the critical surface charges density and on corresponding electric field are analyzed. The contribution of viscosity with the existence of surface charges and streaming are discussed. The investigation includes the stability analysis of the presence of the periodic electric field as well as the constant one. It is found that the presence of the surface charges made by the normal electric field play a dual role in the stability criterion, which shows some analogy with the nonlinear theory of stability. Some previous studies are compared using appropriate data. The marginal state of stability is also considered. It is found that the surface charges vanish under certain conditions. This study shows that the mass and heat transfer parameter has a destabilizing effect whether the electric field is static or periodic. Parametric excitation of the electrohydrodynamic (EHD) surface waves is analyzed in the case of Rayleigh–Taylor (R–T) instability. The transition curves are obtained by means of Whittaker's technique. The analytical results are numerically confirmed.     

STUDY OF CRYSTAL GROWTH AND SOLUTE PRECIPITATION THROUGH FRONT TRACKING METHOD

Crystal growth and solute precipitation is a Stefan problem. It is a free boundary problem for a parabolic partial differential equation with a time-dependent phase interface. The velocity of the moving interface between solute and crystal is a local function. The dendritic structure of the crystal interface, which develops dynamically, requires high resolution of the interface geometry. These facts make the Lagrangian front tracking method well suited for the problem. In this paper, we introduce an upgraded version of the front tracking code and its associated algorithms for the numerical study of crystal formation. We compare our results with the smoothed particle hydrodynamics method （SPH） in terms of the crystal fractal dimension with its dependence on the Damkohler number and density ratio.