水平板自然对流换热的非线性特性

采用SIMPLE算法,QUICK差分方案,对封闭方腔内水平板自然对流换热进行了数值模拟.数值结果显示,低Ra数时流动和换热处于稳态,当Rayleigh数超过某一临界值时,流动和换热就会发生非稳态振荡,此时流动和换热表现出非对称性.对不同Rayleigh数,流动和换热通过单周期分岔从稳态过渡到非稳态,并通过倍周期分岔过渡到混沌.在混沌区,仍然会出现周期性窗口,并且数值结果与初始条件有关.     

水平空气层自然对流换热的分岔和振荡

本文用SIMPLE算法对底部加热的水平空气层的自然对流换热进行了数值计算,研究了这种空气层的流动与换热数值解的振荡和分岔问题。结果表明,对流与换热存在分岔情况。分岔存在一个临界Ra。分岔的临界值与Pr相关,随着Pr的增大,其相应的临界Ra也增大。但当Ra取到5×10~6,这种空气层的对流和换热没有发生振荡。     

锌锅中低Pr流体混合对流的自维持振荡

本文对以带钢连续热镀锌为背景抽象出的锌锅中低Pr流体的流动和传热进行了数值模拟.数值结果显示,对Re=0的纯自然对流,Ra在104和105间时数值解由稳态解分岔为振荡解;对于Ra=0的纯强制对流,Re在4×103和5x1003之间时,数值解从稳态解分岔为振荡解;当Ra和Re均小于自然对流和强制对流单一机制作用时发生振荡的临界值时,混合对流的流动和传热为稳态;当Ra和Re中的一个参数大于单一机制作用发生振荡的临界值时,另一个参数由小于增加到大于单一机制作用发生振荡的临界值时,数值解由振荡解变为稳态解.     

方腔内竖直板自然对流分岔现象的数值研究

采用非稳态数学模型通过数值模拟研究了具有对称结构的封闭方腔内竖直板自然对流换热中的非线性现象。数值结果表明,Ra在4×104-7×104之间,会出现定常解、周期性振荡解和非周期性振荡解,即解存在分岔;初始条件与解的分岔相关;网格Pe数的考核表明,本文所考虑的解的振荡不是由于数值解的不稳定而引起的振荡,而是客观存在的物理振荡。     

锌锅中低Pr流体混合对流的数值模拟

以带钢连续热镀锌生产工艺为背景,对抽象出的低Pr数流体混合对流流动和换热模型进行了数值模拟,给出了在不同Re、Ra及Ri时的流场和温度场.数值结果表明,当Re、Ra都不等于0时,在所考虑的参数范围内,流动和换热受自然对流和强制对流两种机理控制.Re不变,增大Ri,自然对流作用加强,并且当Ri增加到一定值时,流动和换热发生振荡.所给出的速度相图显示,对应不同的Re、Ra及Ri,流动和换热会出现稳态解、周期性振荡解和混沌.     

紧致方法对流动换热及静态分岔的模拟

发展了基于投影法的紧致方法求解流动换热问题,对顶盖驱动流和侧壁加热的方腔内自然对流换热问题进行了数值模拟。与其它传统方法相比,紧致方法能在较少的网格结点下获得精度较高的计算结果。进一步,采用所发展的紧致方法对不同工况下的Rayleigh-Benard对流及其静态分岔现象进行了数值模拟。数值计算结果表明当长宽比变大时,底部努塞尔数会有小幅度增加。当长宽比为8时,用所发展的紧致方法不同的初场可以得出三种不同的流场和温度场。与基于QUICK格式的SIMPLE算法相比,所发展的紧致方法可以多预测一种静态分岔现象。     

周期性边界条件下多孔介质方腔内自然对流换热数值研究

对充满多孔介质的倾斜方腔,在其上下壁面绝热、右侧壁面维持恒温To及左侧壁面温度基于To随时间按正弦规律变化的情况下,采用Brinkman扩展达西模型和SIMPLER算法对方腔内的自然对流与换热进行了非稳态数值模拟,方腔倾角α的范围为0°～90°,Pr数为1,Ra数为106.计算研究不同的振荡频率和方腔倾角对方腔内对流换热的影响.数值模拟结果表明:振荡频率f为60π,方腔倾角为43°时,方腔内的换热最强.     

封闭圆内开缝圆自然对流换热的振荡特性

本文通过数值计算探讨了封闭圆内开缝圆自然对流换热的振荡特性。数值计算以整个圆为计算区域,采用了非稳态的数学模型和具有QUICK差分格式的SIMPLE算法。在相同条件下计算结果和实验结果符合很好。数值结果显示, 当几何结构一定时,Rayleigh数Ra小于某个临界值时,流动和换热处于稳态,并且关于垂直中心线对称;Ra大于这个临界值时,流动和换热是振荡的,非对称的。数值实验还表明,流动和换热出现振荡时的临界Rayleigh数Rac与开缝圆的开缝度有关,且流动和换热的振荡会出现对称振荡和非对称振荡两种情形。     

自然对流的LB方法模拟及其非线性特性

采用LB方法和QUICK差分方法模拟了方腔内竖直平板自然对流和底部加热方腔内自然对流换热问题.实现了LB方法对具有孤立体的封闭空间内耦合的流动和传热问题的数值模拟,所得数值结果与QUICK差分方法的数值结果及已有的烟可视化实验结果一致;两种方法对底部加热方腔内自然对流问题的预测结果均出现了静态分岔和动态分岔;根据底部加热方腔内自然对流换热数值结果给出的极限环型速度相图和功率谱表明,两种方法得到的数值结果的非线性特性一致.     

倾角对方腔内热对流非线性演化与分岔的影响

倾斜封闭腔内对流换热问题是非线性非平衡系统中研究的热点问题之一.本文采用高精度数值方法对倾斜方腔内流体热对流进行了直接数值模拟,研究了腔体倾角在0°—180°之间变化时,倾角的不同变化过程对流场非线性演化、传热效率以及流动分岔的影响.所考虑的Rayleigh数范围为10~3—10~6.结果表明:表征传热效率的Nusselt数对Rayleigh数、Prandtl数及倾斜角度均具有较强依赖性,在较高Rayleigh数时, Nusselt数会在80°和100°附近产生较大幅度的变化;高Rayleigh数下流场及温度场的演变更为复杂,腔体内存在1—3个对流强度不等的涡卷;低Rayleigh数下腔体倾角接近90°时流动状态为热传导状态.当腔体倾角介于70°—110°之间时,在Rayleigh数Ra∈(4949, 314721)内存在解的两条稳定分支.     

液封液桥内振荡热毛细对流的三维数值模拟

为了了解微重力下液封液桥内热毛细对流的基本特性,利用有限差分法进行了非稳态三维数值模拟,液桥高为(1-3)mm,直径为2mm和3 mm,液封外直径为(4-7)mm。模拟结果表明,当Marangoni数较小时,液封液桥内的热毛细对流为稳定的轴对称运动,当Marangoni数超过某一临界值后,流动将转化为三维振荡流动;为此,确定了发生振荡的临界Marangoni数,分析了各种条件下热毛细对流的振荡特性,计算了相应的振荡频率。     

Modelling of a closed two-dimensional reactor bounded by a wavy wall

In this paper, we propose a model for a two-dimensional closed reactor bounded by a wavy wall. The left, right and top walls of the reactor are assumed to be flat surfaces while the bottom wall is a wavy surface. In order to formulate a model for such a reactor, we introduce a coordinate transformation into the dimensionless equations of a rectangular closed domain. Then the resulting equations illustrate the phenomena for a closed reactor bounded by a wavy wall. We solve these equations using the finite difference method. The astonishing results are that the intensity of streamlines and the maximum temperature within the reactor significantly increase with an increase of the number of waves in the bottom wall, the amplitude of waves and the Frank-Kamenetskii number. Converse characteristics are observed for higher values of the enlargement of a wave. Moreover, larger Rayleigh number induces stronger vortices in the flow field and reduces the maximum temperature. The Nusselt number at the bottom wavy wall is found to increase for higher values of the Frank-Kamenetskii number and the amplitude of a wave. A transition from the steady-state to the oscillatory convection is identified for a certain value of the Frank-Kamenetskii number. However, for a low value of the Rayleigh number, there occurs a transition from the steady-state to an explosion for increasing value of the Frank-Kamenetskii number. Results also demonstrate that the critical value of the Frank-Kamenetskii number, for which a transition from the steady-state to the oscillatory convection occurs, is higher for increasing values of the number of waves, the enlargement of a wave and the amplitude of a wave.     

圆内开缝圆自然对流的Ruelle-Takens混沌道路

使用热格子Boltzmann方法针对圆内开缝圆自然对流的流动与换热进行数值模拟，通过相空间、功率谱等进行非线性动力学特性分析，研究其流动与换热的稳定性.结果表明：随着瑞利数Ra的增加，流场的相图从开始稳定的平衡点经历Hopf分岔后转变为极限环，表明流场进入一个倍周期性振荡状态；随着瑞利数进一步增加，稳定的极限环分岔为二维环面，系统相空间结构复杂化；当瑞利数Ra大于某一临界值时，二维环面分岔突变进入混沌状态，系统在相空间中出现非常复杂的轨线结构.总体上，通过系统不同瑞利数所对应的非线性动力学特性的表现形式，表明系统经过Ruelle-Takens道路到达混沌，展现出自然对流从稳定的流动和换热发展到非线性运动特征的混沌历程.     

Transition from periodic to non-periodic oscillation observed in a mathematical model of bioconvection by motile micro-organisms

Bioconvection observed in a culture of motile micro-organisms was analyzed numerically. The governing equations are the Navier-Stokes equations with the Boussinesq approximation and a diffusion equation for the motile micro-organism. A transition from a static condition to periodic oscillation was observed according to the increase of the Rayleigh number. It was found that the system of bioconvection could be led into chaotic conditions via a single-frequency oscillatory behavior to a sequence of period-doubling bifurcations by increasing the Rayleigh number, which is analogous to Bénard convection.     

MAC类方法计算方腔内自然对流

本文总结了用MAC类方法计算二维方腔内自然对流传热的经验.改进了MAC方法的边界处理并提高其精度。用Leapfrog Dufoft-Frankel格式代替MAC方法中的FTCS格式,来放宽稳定性限制。给出高Raylei吵数10⁶的数值管果和比较。讨论了初始温度分布对不定常过程的影响,不定常过程的流动特性符合理论结论。最后还讨论了流出边界处理。     

Instability of single-frequency operation in semilinear photorefractive coherent oscillators

The transition of the single-frequency oscillation of a semilinear photorefractive coherent oscillator for sufficiently large coupling strengths into two-frequency oscillation is predicted and is observed experimentally. The critical value of the coupling strength at which the bifurcation occurs is a function of pump intensity ratio and cavity losses. The supercritical bifurcation in the oscillation spectrum is analogous to the second-order phase transition.     

On the initiation of a corona discharge near the surface of a nonlinearly oscillating liquid layer on the surface of a charged hailstone

An expression is derived for the electric field strength near a wet hailstone in an approximation quadratic in the oscillation amplitude of a charged liquid layer on its surface. It is found that the electric field strength in a small neighborhood of the capillary wave crests grows with the number of a mode governing the initial deformation of the equilibrium (spherical) shape of the liquid layer. Even if the charge is small (when the Rayleigh parameter of the hailstone equals one-hundredth of the value critical for stability against the self-charge), the electric field near the hailstone is high enough for initiating a corona discharge in its vicinity.     

Deformation and nonlinear corrections to critical conditions of charged drop instability in an electrostatic suspension

Nonlinear asymptotic analysis of a charged drop placed in electrostatic and gravitational fields reveals a correction to the oscillation frequency and, accordingly, to the critical Rayleigh parameter. The analysis uses approximations quadratic in oscillation amplitude and linear in dimensionless equilibrium deformation of the drop. The correction is found to be proportional to the product of the oscillation amplitude and deformation. It is natural to name this correction deformational. In computations of the third order of smallness in oscillation dimensionless amplitude, a correction to the frequency and Rayleigh parameter appears, which is due to a nonlinear interaction between oscillation modes. This correction is larger than the deformational one in magnitude. Deformational corrections can be eliminated by experimenting under no-gravity conditions, but corrections due to the nonlinearity of hydrodynamic equations cannot be eliminated in this way. It is these corrections that are responsible for a critical Rayleigh parameter measurement inaccuracy.