液固两相圆柱绕流尾迹内颗粒扩散分布的离散涡数值研究

基于离散涡方法求得的非定常水流场和颗粒的Lagrange运动方程,数值模拟了稀疏液固两相圆柱绕流尾迹内颗粒的扩散分布．获得了流动的涡谱与3种不同St数颗粒（St=0.25,1.0,40）在流场中的分布．通过引入扩散函数来定量表示颗粒在流场中的纵向扩散强度,并计算得到了不同St数颗粒的扩散函数随时间的变化．数值结果揭示出了液固两相圆柱绕流尾迹中的颗粒扩散分布与颗粒的St数和尾涡结构密切相关:1) 中小St数(St=0.25～4.0)颗粒在运动过程中不能进入涡核区,而在旋涡结构的外沿聚集,且颗粒的St数愈大,其越远离涡核区域;2） 在圆柱绕流尾迹区域内,中小St数(St=0.25～4.0)颗粒的纵向扩散强度随其St数的增大而减小．     

三角翼大迎角流场结构和气动力特性的计算分析研究

采用数值计算方法对亚音速三角翼纵向及带有小侧滑和横侧小扰动情况下的流场结构进行了计算,利用数值计算所得到的大迎角流动流场数据,结合相关的实验研究结果,建立了对大迎角旋涡流场结构进行定量分析的方法．给出了三角翼大迎角情况下相应的气动力、力矩系数,以及机翼前缘分离涡轴线位置和旋涡破裂位置随迎角的变化规律,并对带有小侧滑和横侧小扰动情况下对横侧力矩的影响进行了计算与分析．计算结果表明,在前缘分离涡破裂前的上游旋涡区内,前缘分离涡轴线基本保持为直线,且随着迎角增加,前缘分离涡轴线位置愈靠近翼根,并远离翼面;在前缘分离涡破裂的初始阶段,于旋涡轴线处,压力系数会迅速增加,沿涡轴线方向速度迅速减小,在垂直于流向的截面内,愈靠近涡轴线处,沿涡轴线方向速度愈小,甚至出现负值,说明沿涡轴线方向出现回流．当绕机翼上表面前缘分离涡破裂后,将会导致横侧运动不稳定,如果受到小扰动,将产生横侧力矩发散．     

非定常二维Navier-Stokes方程流函数-涡度形式的特征-混合有限元数值模拟

本文针对非定常的二维Navier-Stokes方程提出了一种基于流函数一涡度形式的数值模拟方法——特征一混合有限元方法，得到了流函数、涡度函数和流场速度的最优阶的L^2误差估计．     

基于高阶空间精度格式和嵌套网格对旋翼尾迹捕捉的改进

发展了一种基于高阶迎风格式和嵌套网格捕捉直升机悬停旋翼涡尾迹的方法．无粘通量采用Roe Reimann求解器,使用改进的5阶加权基本无振荡(WENO)格式对交界面左右状态进行高阶插值,并与MUSCL插值进行比较．为便于捕捉尾迹和实施周期性边界条件,计算采用结构嵌套网格,其中高质量的旋翼网格完全嵌套于背景网格中．当解达到近似收敛后在桨尖涡分布区域对背景网格进行加密,如此经过3次得到优化的背景网格．考虑到WENO格式插值的特点,提出了搜索3层洞边界和人工外边界的方法以便插值的直接进行．用该方法对一跨音速和一亚音速悬停旋翼粘性流场进行了数值计算．数值结果表明:所发展方法对涡尾迹具有很高的捕捉能力;与MUSCL格式相比,WENO格式具有较低的数值耗散．     

湍流边界层内准流向发卡涡生成的数值模拟

采用根据共振三波理论模型建立的初始准二维流场和直接数值模拟方法,对湍流边界层近壁区流向涡的生成进行了研究．计算得到了准流向发卡涡和次生准流向涡结构的生成过程及其主要特征,讨论了它们的产生机理．作为相干结构的主要特征,利用共振三波理论模型对流向涡结构生成与演化过程的研究为湍流边界层内相干结构的研究与流动控制提供了新的途径和思路．     

二维无界自由衰减流的数值研究

无界区域上的流体运动是流体力学中的热点和难点问题.采用传统的扩大计算区域算法和新发展的基于无界区域的Hermite基函数算法对二维无界区域的自由衰减流动进行研究. 结果发现,对于只存在相同符号涡的初始流场而言,两种方法都可以得出正确的结果;而对于正负涡都存在的初始流场,传统方法即便利用非常大的计算区域也无法进行正确的长时间模拟,但是新方法却能高效求解.对算例的Hermite算法数值模拟验证了理论解Oseen涡的存在.     

螺旋型旋风分离器两相流场的数值模拟

对螺旋型旋风分离器进行了两相流场的三维数值模拟．气体流场通过求解三维N-S方程得到,湍流模型采用了雷诺应力模型．计算结果表明,旋风分离器内部的流场分为两部分:螺旋通道内比较稳定的流场和筒体中心区域的复合涡结构流场．对颗粒运动轨迹的计算表明,颗粒在入口处的初始位置对颗粒分离有比较显著的影响．同时得到了不同入口速度下颗粒的分级效率曲线,并给出了气体流量对旋风分离器性能的影响,结果显示:气体流量的增加会提高分离效率,但同时导致压力损失的急剧增加．     

浅水横流中异重冲击射流的大尺度涡结构

采用RNG湍流模型对浅水横流中异重冲击射流的大尺度涡结构进行了详细的数值研究．分析了冲击区滞止点上游壁面涡结构和近区Scarf涡结构的尺度、形成机理和演化特征．计算得到了上游壁面涡的特征尺度,结果表明上游壁面涡具有高度的三维性,其特征尺度依赖于流速比和环境水深．近区Sarf涡结构对横流冲击射流的横向浓度分布具有重要的影响．当流速比相对较小时,在底层壁射流与环境横流的横向边界附近出现明显的高浓度聚集现象,计算结果表明Scarf涡结构对这一高浓度聚集区的形成起主导作用．     

不可压缩Navier-Stokes方程数值模拟

建立不可压缩Navier-Stokes方程的Crank-Nicolson有限差分方法,数值模拟了在初始正弦波下的二维水槽内流体受到倾斜激励时流场和涡的演变机制.数值结果表明,初始波为正弦波时,流场内出现一个单涡,单涡下沉变成了两个小涡,两个小涡消失后流场内部出现三条规则的流场带,最后这三条流场带演变成一个尖涡,尖涡在周围流体的作用下演变成一个单涡,最后单涡在自由面消失,当耗散系数和Reynolds数增大时,流场和涡演化的周期变小.     

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

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

Heat transfer enhancement for boundary layer flow over a wedge by the use of electric fields

The effect of an electric field is investigated for heat transfer properties in a laminar, incompressible, non-isothermal boundary layer gas flow over a wedge. The governing boundary layer equations are reduced to an ordinary differential equation system using similarity transformations. The reduced equations are solved numerically for different values of electric and flow field parameters characterizing the ratio of electric force to fluid inertia force, Joule heating and ion kinetic work. For specific electric field function forms, leading to similarity solutions, velocity boundary layers are observed to become thinner and heat transfer properties are shown to be enhanced near the wall. The level of enhancement is controlled by the electric body force with additional effects of Joule heating and ion kinetic work on the bulk flow. The effects of low and high Prandtl numbers are also demonstrated. Heat transfer enhancement is observed to increase with increasing Prandtl number.     

细长体截面绕流中的一种临界状态

应用拓扑分析的方法研究了细长体截面绕流拓扑结构的演变过程．指出随着细长体背涡的发展,导致截面流场的拓扑结构发生变化,会出现一种临界流动状态．在这种临界流态下,流场中会出现一种高阶奇点．这种高阶奇点的指数为-3/2．这种高阶奇点是结构不稳定的,稍有扰动就会产生分叉,使流场的拓扑结构发生变化．     

Navier–Stokes equations with degenerate viscosity,vacuum and gravitational force

The global existence of weak solutions to the compressible Navier–Stokes equations with vacuum attracts many research interests nowadays. For the isentropic gas, the viscosity coefficient depends on density function from physical point of view. When the density function connects to vacuum continuously, the vacuum degeneracy gives some analytic difficulties in proving global existence. In this paper, we consider this case with gravitational force and fixed boundary condition. By giving a series of a priori estimates on the solution coping with the degeneracy of vacuum, gravitational force and boundary effect, we give global existence and uniqueness results similar to the case without force and boundary. Copyright © 2006 John Wiley & Sons, Ltd.     

Exact and approximate expressions for resistance coefficients of a colloidal sphere approaching a solid surface at intermediate Reynolds numbers

A mathematical model has been developed to describe the force of liquid flow acting on a colloidal spherical particle as it approaches a solid surface at intermediate-Reynolds-number-flow regime. The model has incorporated bispherical coordinates to determine a stream function for the flow disturbed by the sphere. The stream function was then used to derive the flow force on the particle as a function of the inter-surface separation distance. The force equation was related to the modified Stokes equation to obtain an exact analytical expression for the correction factor to the Stokes law. Finally, a rational approximation is presented, which is in good agreement with the exact numerical result, and can be readily applied to more general particle–surface interactions involving short-range hydrodynamics associated with colloidal particles in the near vicinity of a large solid collector surface at intermediate Reynolds number of the supporting flow.     

Stress and heat flux stabilization for viscous compressible medium equations with a nonmonotone state function

We study a system of quasilinear equations describing one-dimensional flow of a viscous compressible heat-conducting medium with a nonmonotone state function and mass force. The large-time behavior of solutions is considered for arbitrarily large initial data. In spite of possible nonuniqueness and discontinuity of the stationary solution, we prove L²-stabilization for the stress and heat flux as t → ∞ along with corresponding global energy estimates for them. The new method of proof utilizes a combination of energy type equalities for the stress and heat flux. Consequently, H¹-stabilization of the velocity and temperature along with global estimates for their derivatives are valid as well.     

细长体截面绕流的稳定性及扰动的影响

应用拓扑结构的稳定性理论,分析了细长旋成体截面绕流的结构稳定性．在分析时取极限流线作为流场的内边界,并证明极限流线的鞍点-鞍点连接是拓扑结构稳定的A·D2通过分析发现,由于旋成体背涡的发展,导致截面流场拓扑结构变化,由稳定对称旋涡流态变成不稳定对称旋涡流态．此时流场中存在空间的鞍点-鞍点连接的不稳定拓扑结构,在小扰动下出现分叉,变成稳定非对称旋涡流态,形成非对称背涡．并应用开折理论分析了扰动对流场结构的影响．     

蛋白质空间结构预测的一种优化模型及算法

用理论方法预测蛋白质结构有两个难点,第一是要有一个合理的势函数,第二是要有一个有效的寻优方法找到势函数的全局极小点,本文采用联合残基力场建立了蛋白质空间结构预测模型,然后用我们给出的一种改进模拟退火算法搜索势函数的全局极小点,对脑啡肽的空间结构进行了预测和分析。     

The critical cross-section of a vortex

Summary This paper discusses the problem of critical-flow cross-sections in vortex flows. It is shown that there are two different types of vortex flows, A-type and B-type vortices (say). An A-type vortex approaches its critical flow state as its cross-sectional area increases and departs from the critical state as the cross-sectional area is decreased. This property is associated with the particular dependence of total pressure and circulation on the stream function, and it holds for both subcritical and supercritical A-type vortices. On the other hand, both subcritical and supercritical B-type vortices approach their critical flow states as their cross-sectional areas are decreased and depart from their critical states for increasing cross-sectional area. As was shown by Benjamin, setting the first variation of the flow force with respect to the stream function equal to zero leads to Euler's equation of motion. The second variation also vanishes if the corresponding flow state is critical. In this case the sign of the third variation decides whether the flow is an A-type or a B-type vortex. Within the framework of inviscid-fluid flow theory the type of a vortex is preserved unless vortex breakdown occurs. Making use of the knowledge that vortex flows are controlled by two different types of critical-flow cross-sections a variety of vortex flow phenomena are investigated, including the two types of inlet vortices that are observed upstream of jet engines, the behavior of vortex valves, the flow characteristics of liquid-fuel atomizers and the bath tub vortex.