底部加热长方体腔内自然对流的非线性特性

采用SIMPLE算法,QUICK差分格式,对底部加热三维长方体腔内空气的自然对流进行了数值模拟。根据模拟结果,探讨了方腔内流体流动与换热的静态分岔与振荡等非线性现象。数值结果显示,在固定的几何尺寸和不同Ra的情况下,当初始场不同时,会出现若干不同的解,即存在解的静态分岔;在固定的几何尺寸和相同的初始场情况下,低Ra时流动和换热处于稳态,当Ra超过某一临界值时,流动和换热就会随时间振荡,并通过倍周期分岔过渡到混沌;当方腔的几何尺寸不同时,分岔点的特征值Ra也发生变化。     

三维紧致修正方法和圆筒内自然对流数值模拟

本文发展了差分法求解流动与换热问题的三维非均分网格7点紧致格式,并利用延迟修正方法将其与SIMPLE算法相结合形成了一种三维紧致修正方法。利用所发展的紧致修正方法对圆筒内同心开缝圆筒的三维自然对流与换热问题进行了数值模拟,所获得的数值结果与实验结果一致。采用Richardson方法证实所发展的三维紧致修正方法大约具有4阶精度。进一步的数值计算表明,在特征参数Ra数大于一定值时,由圆筒内同心开缝圆筒的三维自然对流问题简化的二维模型数值结果与实验结果逐渐加大,用三维模型才能得到比较可靠的结果。     

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

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

基于SIMPLE的紧致方法

通过泰勒展式系数匹配的方法发展了基于非等距网格的有限容积紧致格式,采用延迟修正的方法建立了基于SIMPLE的紧致方法,,该方法能够得到高精度的数值解,增加迭代求解代数方程组的稳定性。对底部加热的方腔内自然对流换热问题进行数值模拟,结果表明,紧致方法比二阶中心差分方法具有更高的精度。     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

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

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

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

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

具有隔板的平行通道内空气混合对流换热数值模拟

本文以锅炉干排渣装置为背景,对抽象的理论模型具有隔板的平行通道内空气混合对流换热进行了数值模拟.数值计算表明,在Re＞1000时应采用非稳态数学模型进行数值模拟;在Re＞500时,自然对流机制对流动和换热的影响基本可以忽略.数值计算给出了不同Re时的进出口无量纲压差、局部的Nux和平均Nu以及流线图.这些结果可为深入研究干排渣装置中流动和换热特性提供参考.     

振荡流共轭换热格子——Boltzmann模拟

振荡流共轭换热现象广泛存在于热声热机等工程应用中.基于双分布格子-Boltzmann模型,对平行平板间振荡流共轭换热进行了数值模拟.通过假定共轭界面处流体和固体的未知内能分布函数均为对应的平衡态滑移修正格式,提出了一种处理共轭换热边界的新方法.模拟结果表明,该方法可以保证共轭界面上温度连续和热流连续.分析了不同流体与固体导热系数比情况下振荡流共轭换热的速度场、温度场以及热流分布的特点.     

方形空间内混合对流换热的数值研究

以建筑物内人工环境控制为应用背景,对有对称空气射流的方形空间内混合对流换热进行了数值模拟,探讨了这种具有对称结构的混合对流换热解的分岔问题。数值结果表明,Ｒｅｙｎｏｌｄｓ数、强制通风气流的射流角度、以及方形空间的宽高比都会影响解的分岔。当Ｒｅ数超过某一临界值时,会出现非对称数值解。宽高比减小,出现非对称解的临界Ｒｅ数也随之减小。Ｒｅ数、宽高比一定,仅当通风气流的射流角度在某个范围内时,能够得到非对称的数值解。     

Transition to chaos in lid-driven square cavity flow

To date, there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity, due to the difficulties in theoretical analysis and numerical simulations. In this paper, we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us. We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of (13944.7021,13946.5333) by the method of bisection. Through Fourier analysis, it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval. Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram, Kolmogorov entropy and maximal Lyapunov exponent. The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.     

The second Hopf bifurcation in lid-driven square cavity

To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs.     

三维方腔温盐双扩散的格子Boltzmann方法数值模拟

研究了温盐双扩散系统的多组分格子Boltzmann方法.通过对二维方腔的温盐双扩散系统的数值模拟，检验了方法的可行性及有效性，所得到的结果与差分法结果符合良好，继而将此方法推广到三维，建立了三维温盐双扩散系统的格子Boltzmann方法，对三维方腔双扩散问题进行了模拟和分析，并与差分法模拟的结果进行了比较，结果令人满意.最后，分析了格子Boltzmann方法在模拟双扩散对流问题时存在的局限性. 关键词： 格子Boltzmann方法 温盐双扩散 Boussinesq近似 数值模拟     

格子Boltzmann方法中三维运动边界的统一模型

格子Boltzmann方法(LBM)中边界条件的处理很复杂,在现有的边界条件处理方法中,动力学格式能够精确满足宏观边界条件,但由于要解一个不定方程,必须引入附加假设确保方程非奇异.作为动力学格式和反弹格式的一种扩展,提出一种处理三维任意速度运动边界的统一模型,其中人口速度和固体壁面速度是该模型的特殊情形.给出用于三维15速度的表达式.为了检验该模型,模拟对角顶盖驱动三维空腔流,并将结果与有限差分法计算的结果进行比较,说明所提出的统一模型是合理可行的.     

Fully HOC Scheme for Mixed Convection Flow in a Lid-Driven Cavity Filled with a Nanofluid

A fully higher-order compact (HOC) finite difference scheme on the 9-point two-dimensional (2D) stencil is formulated for solving the steady-state laminar mixed convection flow in a lid-driven inclined square enclosure filled with water-$Al_2O_3$ nanofluid. Two cases are considered depending on the direction of temperature gradient imposed (Case I, top and bottom; Case II, left and right). The developed equations are given in terms of the stream function-vorticity formulation and are non-dimensionalized and then solved numerically by a fourth-order accurate compact finite difference method. Unlike other compact solution procedure in literature for this physical configuration, the present method is fully compact and fully higher-order accurate. The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and averaged Nusselt number. Comparisons with previously published work are performed and found to be in excellent agreement. A parametric study is conducted and a set of graphical results is presented and discussed to elucidate that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

求解Navier-Stokes方程组的组合紧致迎风格式

给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟.     

基于非结构网格的不可压流动耦合求解器

采用非结构化网格有限容积法求解了不可压N-S方程组,对流项采用GAMMA格式,扩散项采用二阶中心差分格式建立离散方程,用SOAR算法处理压力与速度的耦合关系,得到了一种求解不可压N-S方程的非结构网格耦合求解器。通过方腔顶盖驱动流、后台阶绕流以及方腔自然对流等几个典型的算例,考察了求解器的计算精度及收敛特性,并与SIMPLE算法进行了比较,结果表明该求解器是有效可行的。     

A compact finite difference scheme for the fractional sub-diffusion equations

In this paper, a compact finite difference scheme for the fractional sub-diffusion equations is derived. After a transformation of the original problem, the L1 discretization is applied for the time-fractional part and fourth-order accuracy compact approximation for the second-order space derivative. The unique solvability of the difference solution is discussed. The stability and convergence of the finite difference scheme in maximum norm are proved using the energy method, where a new inner product is introduced for the theoretical analysis. The technique is quite novel and different from previous analytical methods. Finally, a numerical example is provided to show the effectiveness and accuracy of the method.     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

Multigrid method based on the transformation-free HOC scheme on nonuniform grids for 2D convection diffusion problems

In this paper, a multigrid method based on the high order compact (HOC) difference scheme on nonuniform grids, which has been proposed by Kalita et al. [J.C. Kalita, A.K. Dass, D.C. Dalal, A transformation-free HOC scheme for steady convection–diffusion on non-uniform grids, Int. J. Numer. Methods Fluids 44 (2004) 33–53], is proposed to solve the two-dimensional (2D) convection diffusion equation. The HOC scheme is not involved in any grid transformation to map the nonuniform grids to uniform grids, consequently, the multigrid method is brand-new for solving the discrete system arising from the difference equation on nonuniform grids. The corresponding multigrid projection and interpolation operators are constructed by the area ratio. Some boundary layer and local singularity problems are used to demonstrate the superiority of the present method. Numerical results show that the multigrid method with the HOC scheme on nonuniform grids almost gets as equally efficient convergence rate as on uniform grids and the computed solution on nonuniform grids retains fourth order accuracy while on uniform grids just gets very poor solution for very steep boundary layer or high local singularity problems. The present method is also applied to solve the 2D incompressible Navier–Stokes equations using the stream function–vorticity formulation and the numerical solutions of the lid-driven cavity flow problem are obtained and compared with solutions available in the literature.