平行平板流动腔脉动流切应力的计算

高度远小于横向和纵向几何尺寸的矩形平行平板流动腔是人们用以体外研究细胞在切应力作用下力学行为的主要工具之一。大多数研究者主要对定常层流情进行研究。本文通过对矩形平行平板流动腔内的层流脉动流进行详细分析,给出腔内速和腔室底部切应力的准确计算公式。当Ｗｏｍｅｒｓｌｅｙ数较小时,准确公式简化为准定常公式。数值计算结果表明,在脉动流条件下,对于人们常用的流动腔几何尺寸,准定常公式具有相当高的精度。这为在脉     

不同属性肋片对复杂方腔自然对流影响的数值分析

采用二阶全展开Euler-Taylor-Galerkin分裂步有限元方法,在指定的网格密度条件下,在流动对应的普朗特数取为0.71,雷诺数取为104的情况下,数值分析了热肋、冷肋、上绝热肋、下绝热肋等四种不同属性肋片对封闭方腔内典型自然对流流动的影响.计算结果表明,肋片的存在对封闭方腔内的自然对流及相应的传热效率具有较强的影响,对流流动结构以及平均Nusselt数随肋片的属性发生较大的改变.     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

长方腔自然对流第一次分岔突变现象的数值分析

在对网格密度作用进行详细分析的基础上,采用二阶全展开 Euler-Taylor-Galerkin分裂步有限元方法,对封闭水平矩形腔体内流体 (Pr=0.71)自然对流的第一次分岔过程进行了数值预报. 计算结果表明,第一次分岔相应的流动拓扑及临界Rayleigh ($Ra$)数随矩形腔体长宽比(W/B)取值的不同时会发生较大变化. 在所计算的长宽比取值范围内,封闭矩形腔内,流体自然对流第一次分岔拓扑的变化对应两种大的类型: 在较小的长宽比取值范围内(W/B\le 2.5),临界Ra数两侧,流动从单一涡核的定常流动突变成为具有不对称结构的定常双涡核运动, 在此范围内临界Ra数的取值随W/B取值的增加而减小;当对应长宽比取值2.6 \le W/B \le 4.0时,第一次分岔拓扑结构的变化呈现出更加复杂的特性,临界Ra数两侧流动从定常双涡核突变为定常非对称的三涡核流动,相应的临界Ra 数也随W/B的增加而减少. 而在区间[2.5,2.6]两端,临界Ra数的取值发生一次阶跃式突增,将该长宽比取值的区间定义为长方腔内该流体第一次分岔的突变区间.      

密度处理方法在圆柱腔体大温差自然对流数值模拟中的影响

采用Boussinesq近似、不可压-理想气体及密度线性差分三种密度处理方法对大温差(?Tmax=1000K)竖直圆柱腔体的定常自然对流进行了模拟.首先,通过模拟一封闭方腔流动并与文献对比,验证了数值模拟方法的有效性.在此基础上,进一步研究了圆柱腔体中表面曲率κ、瑞利数Ra及温差ΔT的影响,特别是这些因素影响下不同密度处理方法对热壁面平均努谢尔特数Nuave数值结果的影响.结果表明,表面曲率κ是影响圆柱腔体自然对流换热的重要因素,且随着κ的不断减小,其Nuave数会逐渐趋近于方腔对流;同时,在温差ΔT和表面曲率κ较大时,基于Boussinesq近似的方法对Nuave的预测存在较大偏差.     

格拉晓夫数Gr对侧向局部加热腔体内对流结构的影响

通过二维流体力学基本方程组数值模拟,研究了Pr=6.949时侧向局部加热腔体内不同格拉晓夫数Gr对对流结构的影响.发现随着格拉晓夫数Gr增加出现以下四种情况:准定常流→单局部周期滚动→双局部周期滚动→定常流.还分析了周期对流阶段的对流特性.随着格拉晓夫数Gr增大,对流卷随时间移动的周期减小.随Gr增大,在0≤≤Y 5.2和0.5≤≤Y 5.7范围内对流卷逐渐消失,在5.2≤≤Y5和5.7≤≤Y10范围内对流卷由单对流卷变为多对流卷.     

封闭方腔内空气和水自然对流每一次分岔数值预报

采用二阶全展开ETG分裂步有限元方法,通过对流动拓扑的详细分析,在排除网格密度影响的基础上,结合二分法给出封闭方腔内空气和水两种典型流体自然对流发生第一次分岔时的临界Rayleigh数。计算结果表明,该方法可用于进行不同Prandtl数条件下方腔内自然对流流动第一次分岔的数值预报,可作为后续各阶分岔及转捩数值预报研究的基础。在相应的条件下,封闭方腔内空气比水更容易发生分岔,且空气的流动结构相对于水表现出一定的倾斜性。     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

底部周期加热的局部对流斑图

运用Simple算法对二维流体力学基本方程组进行了数值模拟,探讨了普朗特数(Pr)为0.0272时矩形腔体底部周期加热对对流时空斑图的影响。当水平流动雷诺数(Re)为0时,发现了由正弦波周期加热引起的稳定的局部定常对流。当Re≠0时,由于正弦波周期加热与水平流动相互作用,获得了由正弦波周期加热和水平流动引起的局部行波对流。进一步比较和讨论了底部正弦波周期加热局部对流和混合流体Rayleigh-Benard局部对流的时空斑图,发现它们存在不同的机理。     

周期加热的Rayleigh-Benard对流时空结构

通过二维流体力学基本方程组的数值模拟,研究了普朗特数Pr=6.99时矩形渠槽周期加热对Rayleigh-Benard对流时空结构的影响．当水平流动强度Re=0时,发现稳定的由周期加热引起的局部定常对流．当Re比较小时,对流滚动抑制水平流动,获得了由周期加热引起的局部行波对流．当水平流动强度比较大时,由于周期加热与水平流动相互作用,水平流动抑制部分对流滚动,导致对流区域上游附近出现传导区域,对流区域减小,从而形成一种新的局部行波对流结构．并进一步讨论了Rayleigh-Benard对流时空结构的动力学特性．     

Nonlinear behaviour of the instability of steady natural convection in a vertical air layer

In this paper we deal with the flow of natural convection between two vertical planes with horizontal temperature gradient. We show that periodic steady flow occur when Rayleigh number increases. Critical values are obtained numerically and a nonlinear analysis is presented according to Center manifold method.     

Heat and mass transfer analogy for condensation of humid air in a vertical channel

This study examines energy transport associated with liquid film condensation in natural convection flows driven by differences in density due to temperature and concentration gradients. The condensation problem is based on the thin-film assumptions. The most common compositional gradient, which is encountered in humid air at ambient temperature is considered. A steady laminar Boussinesq flow of an ideal gas–vapor mixture is studied for the case of a vertical parallel plate channel. New correlations for the latent and sensible Nusselt numbers are established, and the heat and mass transfer analogy between the sensible Nusselt number and Sherwood number is demonstrated. Received on 15 November 1999     

A high‐order discontinuous Galerkin solver for low Mach number flows

In this work, we present a high‐order discontinuous Galerkin method (DGM) for simulating variable density flows at low Mach numbers. The corresponding low Mach number equations are an approximation of the compressible Navier–Stokes equations in the limit of zero Mach number. To the best of the authors'y knowledge, it is the first time that the DGM is applied to the low Mach number equations. The mixed‐order formulation is applied for spatial discretization. For steady cases, we apply the semi‐implicit method for pressure‐linked equation (SIMPLE) algorithm to solve the non‐linear system in a segregated manner. For unsteady cases, the solver is implicit in time using backward differentiation formulae, and the SIMPLE algorithm is applied to solve the non‐linear system in each time step. Numerical results for the following three test cases are shown: Couette flow with a vertical temperature gradient, natural convection in a square cavity, and unsteady natural convection in a tall cavity. Considering a fixed number of degrees of freedom, the results demonstrate the benefits of using higher approximation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

Magnetic control of natural convection in the horizontal Bridgman configuration: symmetric and non-symmetric cross sections

This study deals with the electromagnetic damping of free-convective flows in cavities such as those used in the crystal growth horizontal Bridgman configuration. The cavities are filled with a dilute electrically conducting alloy and are subjected to a horizontal temperature gradient. The flow is steady and laminar under an external, vertical, transversal and uniform magnetic field. Several cross sections of the cavities were investigated and can either be centro-symmetric or not. The governing equations for such problems are two coupled partial differential equations, for the velocity and the induced magnetic fields, coupled with a third integral equation for mass conservation. A finite element method has been developed, and the numerical results for the variation of the velocity and the induced magnetic field in terms of the Hartmann number show a considerable decrease in convection intensity as the Hartmann number increases. Results also reveal the presence of the well-known Hartmann and parallel layers. For non-centro-symmetric sections, results show the way the flow reorganises into two cells as the Hartmann number increases.     

Dual steady solutions in natural convection between horizontal concentric cylinders

Dual steady solutions in natural convection in an annulus between two horizontal concentric cylinders are numerically investigated for a fluid of Prandtl number 0.7. It is found that, when the Rayleigh number based on the gap width exceeds a certain critical value, dual steady two-dimensional (2-D) flows can be realized: one being the crescent-shaped eddy flow commonly observed and the other the flow consisting of two counter-rotating eddies and their mirror images. The critical Rayleigh number decreases as the inverse relative gap width increases.     

NUMERICAL INVESTIGATION ON FLOW STABILITY OF RAYLEIGH-B?NARD CONVECTION OF COLD WATER IN A RECTANGULAR CAVITY COOLED AND HEATED SYMMETRICALLY RELATIVE TO THE TEMPERATURE OF DENSITY MAXIMUM

为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.     

Double-diffusive natural convection in a porous enclosure partially heated from below and differentially salted

This paper reports numerical results of two-dimensional double-diffusive natural convection in a square porous cavity partially heated from below while its upper surface is cooled at a constant temperature. The vertical walls of the porous matrix are subjected to a horizontal concentration gradient. The parameters governing the problem are the thermal Rayleigh number (Ra=100 and 200), the Lewis number (Le=0.1, 1 and 10), the buoyancy ratio (−10N10) and the relative position of the heating element with respect to the vertical centerline of the cavity (δ=0 and 0.5). The effect of the governing parameters on fluid characteristics is analyzed. The multiplicity of solutions is explored and the existence of asymmetric bicellular flow is proved when the heated element is shifted towards a vertical boundary (δ=0.5). The solutal buoyancy forces induced by horizontal concentration gradient lead to the elimination of the multiplicity of solutions obtained in pure thermal convection when N reaches some threshold value which depends on Le and Ra.     

Natural convection in a reservoir sidearm subject to solar radiation: experimental observations

The natural convection in a reservoir sidearm induced by solar radiation is visualised using a shadowgraph technique. The flow visualisation reveals three stages of the flow development, namely an initial growth stage, a transitional stage and a quasi-steady stage. At the initial growth stage, a distinct thermal boundary layer grows rapidly along the sloping bottom. The transitional stage is characterised by the onset of convective instability in a form of rising plumes. At the quasi-steady state, the mean temperature across the enclosure increases steadily in time and the flow is sighted with quasi-regular presence of instabilities with reduced intensities. Received: 3 July 2001/Accepted: 10 December 2001     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.