充满多孔介质的方腔内双扩散自然对流格子Boltzmann模拟

采用格子Boltzmann方法,对4个壁面均为低温、低浓度,内置高浓度发热圆的充满均匀多孔介质的方腔内双扩散自然对流现象进行了数值模拟研究.分析了Darcy(达西)数Da(10-4≤D_a≤10~(-2))和浮升力比B(-5.0≤B≤5.0)对内部发热圆表面平均Nusselt(努赛尔)数Nuav和平均Sherwood(舍伍德)数Shav的影响.模拟结果表明:除B=-1.0时,Nu_(av)和Sh_(av)随D_a的增加而增大;当-5.0B5.0,在D_a=10~(-4)时,Nu_(av)和Sh_(av)几乎不受B变化的影响;在Da=10~(-3)和Da=10~(-2)时,Nu_(av)和Sh_(av)随B的增加先减小后增大,在B=-1.0时取得最小值.     

基于Soret和Dufour效应的方腔内双扩散自然对流振荡特性研究

采用格子Boltzmann方法模拟研究了方腔内带Soret效应和Dufour效应的双扩散自然对流振荡特性.方腔内置高浓度发热圆且位于方腔中心,四周壁面均为低温低浓度.采用时间历程分析法和功率谱法分析了不同的浮升力比Br(2.0≤B_r≤10.0)、Soret数Sr(-0.6≤S_r≤0.0)和Dufour数Df(-0.6≤D_f≤0.0)下的方腔内部流动的振荡特性.研究结果表明:不考虑Soret和Dufour效应时,方腔内部流动呈现稳定状态,随着Df和Sr从0.0变化到-0.6,双扩散自然对流状态开始逐渐转变为周期性振荡和非周期性振荡,且振荡性随着浮升力比Br的增大而增强.     

竖直平板嵌入非Darcy多孔介质时,数值研究流过平板的不稳定MHD自然对流中的Dufour和Soret效应

就竖直平板嵌入非Darcy多孔介质中,导电流体流过平板时作不稳定的二维磁流体(MHD)双扩散对流,数值研究了Dufour和Soret效应对流动的影响.用Crank-Nicolson型的隐式有限差分法,按三对角矩阵处理,求解无量纲的非线性控制方程.详细地研究了问题中出现的各种参数对不稳定无量纲的速度、温度和浓度曲线的影响.进一步地,给出并分析了表面摩擦因数、Nus-selt数和Sherwood数随时间的变化.研究结果表明,不稳定速度、温度和浓度分布曲线,受Dufour和Soret的影响十分显著.随着Dufour数的增加或者Soret数的减小,表面摩擦因数和Sherwood数都在减小,而Nusselt数在增加.研究发现,当磁场参数增加时,边界层中的速度和温度在减小.     

非Newton流体在多孔介质中沿竖直面作自由对流时应力屈服和Soret效应对其传热/传质的影响

对饱和的非Newton幂律流体,流经多孔介质中竖直平板时的自由对流,在出现应力屈服和Soret效应时,研究化学反映对传热/传质的影响.用相似变换,将边界层控制方程及其边界条件转换为无量纲的形式,然后通过有限差分法求解该方程.给出并讨论了浓度曲线,以及本问题各种参数值时的Nusselt数和Sherwood数.发现化学反应参数γ、化学反应级m、Soret数Sr、浮力比N、Lewis数Le及无量纲流变参数Ω对流场有着显著的影响.     

倾斜层中的对流斑图及其临界条件

通过二维流体力学基本方程的数值模拟，探讨了Prandtl(普朗特)数Pr=6.99时，倾斜矩形腔体中的对流斑图和斑图转换的临界条件.根据倾角θ和相对Rayleigh(瑞利)数Ra_r的变化，倾斜矩形腔体中的对流斑图可以分为:单滚动圈对流斑图、充满腔体的多滚动圈对流斑图和过渡阶段的多滚动圈对流斑图.当θ一定时，随着Ra_r的减小，系统由充满腔体的多滚动圈对流斑图过渡到单滚动圈对流斑图.这时，对流振幅A和Nusselt(努塞尔)数Nu随着Ra_r的增加而增加.当Ra_r=9时，随着θ的增加，系统由充满腔体的多滚动圈对流斑图过渡到单滚动圈对流斑图，这时对流振幅A随着θ的增加而减小，Nusselt数Nu随着θ的增加而增加.在θ_c-Ra_r平面上对多滚动圈到单滚动圈对流斑图过渡的模拟结果表明， 在Ra_r=2时, 腔体中没有发现多滚动圈对流斑图.在Ra_r为2.5左右时，腔体中出现多滚动圈到单滚动圈对流斑图的过渡.当多滚动圈到单滚动圈对流斑图过渡的临界倾角θ_c<10°时，θ_{c随着Rar的减小而增加.当θc>10°时，θc随着Rar的增加而增加，在Rar≤5时，θc随着Rar的增加而迅速增加；当Rar>5时，θc随着Rar的增加而缓慢增加.θc与Rarθ的关系与Rar类似}     

在有一级化学反应时粘弹性流体流经无限竖直多孔平板时的边界层流动

在有一级化学反应时,研究不可压缩的粘弹性流体,在竖直多孔连续运动平板上的不稳定自然对流.控制方程用隐式有限差分法进行数值求解.与解析解的结果比较,证明所选用的数值方法有效.详细图示了速度分布的数值结果.研究了粘弹性参数、无量纲化学反应参数和平板运动速度,对稳定的速度分布、与时间相关的摩擦因数、Nusselt数和Sherwood数的影响.     

一类带约束的二维弱奇异积分方程的解*

本文找出二维弱奇异第一类积分方程作用着约束方程的解p.p=p(r,θ)={2/[π2k(φ₀]}√F(r,θ)-c_*(0≤r≤r_*)其中是(s,φ)原点在M(r,θ)的局部极坐标,(r,θ)是原点在O(0,0)的总体极坐标:k和F是给出的连续函数;φ₀是一常数;F(_*,θ)=c_*(常数)是研究域Q的边界围线。所用方法可推广到三维情形。     

二阶奇异非线性微分方程边值问题的正解

分别在0≤f₀+<M₁,m₁<f_∞-≤∞和0≤f_∞+<M₁,m₁<f₀-≤∞的情形下研究了非线性奇异边值问题u″+g(t)f(u)=0,0<t<1,αu(0)-βu′(0)=0,γu(1)+δu′(1)=0正解的存在性,其中f₀+=₀f(u)/u,f_∞-=_∞f(u)/u,f₀-=₀f(u)/u,f_∞+=_∞f(u)/u,g在区间[0,1]的端点可以具有奇性。     

一种新的混合流体对流竖向镜面对称对传波斑图

利用SIMPLE算法对混合流体对流的流体力学基本方程组进行了数值模拟,在混合流体分离比ψ=－0.6和矩形腔体长高比Γ=20的情况下,首次发现了一种新的竖向镜面对称对传波斑图,并初步探讨了它的动力学特性.竖向镜面对称对传波斑图的中心为驻波,随着时间的发展驻波的波长伸长.当波长增加到某个临界值时,一个滚动分裂成两个滚动,在这两个滚动之间产生一个具有180°相位差的新滚动.位于中心线上的滚动只有相位的突变及其波长的压缩或者伸长,没有对流滚动的移动,在它的两侧是向左右传播的对流滚动.驻波两次相位突变形成一个周期,驻波周期随着相对Rayleigh(瑞利)数Ra_r的增加而增加.这种对流结构存在于相对Rayleigh数Ra_r∈(3.6,4.3]的范围.当相对Rayleigh数Ra_r≤3.6时,系统出现具有缺陷的行波斑图;当Ra_r>4.3时系统过渡到行波斑图.说明竖向镜面对称对传波斑图是存在于具有缺陷的行波斑图和行波斑图之间的一种稳定的对流斑图.     

色散方程的四点显式差分格式

本文对色散方程u_t=au>_xxx构造了一类高稳定性的、在中间层涉及四个网格点的三层显式差分格式,其局部截断误差为O(τ+h),其稳定条件为|R|=|α|τ/h3≤0.25至|R|≤10,它们较大地改善了同类格式的稳定条件|R|≤0.25[1].     

Structural stability for the Brinkman equations of flow in double diffusive convection

This paper continues the investigation of structural stability for the Brinkman equations modeling the double diffusive convection for flow in a porous medium. It supplements earlier results of Straughan and Hutter [B. Straughan, K. Hutter, A priori bounds and structural stability for double diffusive convection incorporating the Soret effect, Proc. R. Soc. Lond. Ser. A 455 (1999) 767-777].     

Exact analytical solution for a thermal boundary layer in a saturated porous medium

The forced convection thermal boundary layer in a porous medium as an analytically tractable special case of a mixed convection problem is considered. It is shown that some general features of the mixed convection solutions reported recently by other authors [B. Brighi, J.-D. Hoernel, On the concave and convex solutions of mixed convection boundary layer approximation in a porous medium, Appl. Math. Lett. (published online, 2005); M. Guedda, Multiple solutions of mixed convection boundary layer approximations in a porous medium, Appl. Math. Lett. (published online, 2005)] can already be recovered from this exactly solvable case.     

Attractors for autonomous double‐diffusive convection systems based on Brinkman–Forchheimer equations

In this paper, we consider the existence of global attractor and exponential attractor for some dynamical system generated by nonlinear parabolic equations in bounded domains with the dimension N≤4 which describe double‐diffusive convection phenomena in a porous medium. We deal with both of homogeneous Dirichlet and Neumann boundary condition cases. Especially, when Neumann condition is imposed, we need some assumptions and restrictions for the external forces and the average of initial data, since the mass conservation law holds. Copyright © 2015 John Wiley & Sons, Ltd.     

An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq fluid, which is heated and salted from below in the presence of Soret coefficient is studied analytically using both linear and nonlinear stability analyses. The normal mode technique is used in the linear stability analysis while a weak nonlinear analysis based on a minimal representation of double Fourier series method is used in the nonlinear analysis. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and frequency of oscillations are obtained analytically using linear theory. The effect of anisotropy parameters, solute Rayleigh number, Soret parameter and Lewis number on the stationary, oscillatory, finite amplitude convection and heat and mass transfer are shown graphically.     

侧向加热腔体中的多圈型对流斑图

基于流体力学方程组的数值模拟,研究了倾角θ=90°时侧向加热的大高宽比腔体中的对流斑图.对于Prandtl数Pr=6.99的流体,在相对Rayleigh数2≤Ra r≤25的范围内,腔体中发生的是单圈型对流斑图.对于Pr=0.0272的流体,取Ra r=13.9,随着计算时间的发展,腔体中由最初的单圈型对流斑图过渡到多圈型对流斑图,这是出现在侧向加热大高宽比腔体中的新型对流斑图.对不同Ra r情况的计算结果表明,Ra r对对流斑图的形成存在明显的影响.当Ra r≤4.4时是单圈型对流滚动;当Ra r=8.9~11.1时是过渡状态;当Ra r≥13.9时是多圈型对流滚动.对流最大振幅和Nusselt数Nu随着相对Rayleigh数的增加而增加.该对流斑图与Pr=6.99时对流斑图的比较说明,对流斑图的形成依赖于Prandtl数.     

Unsteady convective diffusion of a solute in a Hagen-Poiseuille flow through a tube with permeable wall

Influence of Interphase Mass Transfer (IMT) on the unsteady convective diffusion in a fluid flow through a tube surrounded by a porous medium is examined against the background of no IMT. The three coefficients namely exchange coefficient, convection coefficient, and dispersion coefficient are evaluated asymptotically at large-time. The exchange coefficient exists due to IMT. All-time analysis is made analytically when there is no IMT. The mean concentration distribution is measured at a point inside and outside the slug. The peak of mean concentration is higher than that of pure convection and it is further enhanced with increase of porous parameter. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermohaline convection with cross-diffusion in an anisotropic porous medium

Using normal mode technique it has been shown that (i) values of the anisotropy parameter are important in deciding the mode of convection in a doubly diffusive fluid saturating a porous medium. (ii) Depending on the values of the Soret and Dufour parameters, an increase in anisotropy parameter either promotes or inhibits instability, (iii) cross-diffusion induces instability even in a potentially stable set-up and (iv) for certain values of the Dufour and Soret parameters there is a discontinuity in the critical thermal Rayleigh number, which disappears if the porous medium has horizontal isotropy.     

Onset of triply diffusive convection in a Maxwell fluid saturated porous layer

The linear stability of triply diffusive convection in a binary Maxwell fluid saturated porous layer is investigated. Applying the normal mode method theory, the criterion for the onset of stationary and oscillatory convection is obtained. The modified Darcy–Maxwell model is used as the analysis model, this allows us to make a thorough investigation of the processes of viscoelasticity and diffusions that causes the convection to set in through oscillatory rather than stationary. The effects of Vadasz number, normalized porosity parameter, relaxation parameter, Lewis number and solute Rayleigh number on the system are presented numerically and graphically.     

On the concave and convex solutions of a mixed convection boundary layer approximation in a porous medium

In this paper we investigate the similarity solutions of a plane mixed convection boundary layer flow near a semi-vertical plate, with a prescribed power law function of the distance from the leading edge for the temperature, that is embedded in a porous medium. We show the existence and uniqueness of convex and concave solutions for positive values of the power law exponent.     

Magnetohydrodynamic non-Darcy mixed convection heat transfer from a vertical heated plate embedded in a porous medium with variable porosity

A numerical model is developed to study magnetohydrodynamics (MHD) mixed convection from a heated vertical plate embedded in a Newtonian fluid saturated sparsely packed porous medium by considering the variation of permeability, porosity and thermal conductivity. The boundary layer flow in the porous medium is governed by Forchheimer–Brinkman extended Darcy model. The conservation equations that govern the problem are reduced to a system of non-linear ordinary differential equations by using similarity transformations. Because of non-linearity, the governing equations are solved numerically. The effects of magnetic field on velocity and temperature distributions are studied in detail by considering uniform permeability (UP) and variable permeability (VP) of the porous medium and the results are discussed graphically. Besides, skin friction and Nusselt number are also computed for various physical parameters governing the problem under consideration. It is found that the inertial parameter has a significant influence in increasing the flow field and the rate of heat transfer for variable permeability case. The important finding of the present work is that the magnetic field has considerable effects on the boundary layer velocity and on the rate of heat transfer for variable permeability of the porous medium. Further, the results obtained under the limiting conditions were found to be in good agreement with the existing ones.