基于中心差分的对流扩散方程四阶紧凑格式

在经典中心差分格式的基础上,提出对流扩散方程的四阶紧凑差分格式。具体方法是,先就一维情形,将中心差分格式改造为不受网格Reynolds数限制的恒稳二阶格式,再在不增加相关网格点的前提下,通过格式中对流系数和源项的摄动处理,使稳格式的精度提高至四阶。本文并作一、二、三维流动模型方程及高Rayleigh数自然对流传热问题的数值求解,例示本文格式的优良性态。     

含不同形状内热源的圆管内纳米流体自然对流及换热数值研究

采用格子玻尔兹曼方法对有三种恒温热源（圆形、三角形、方形）参与的圆管内纳米流体（铜-水）自然对流进行数值研究。主要研究瑞利（Ra）数,纳米颗粒体积分数以及热源几何形状等控制参数对纳米流体的流动与传热的影响。结果发现纳米颗粒体积分数的增加有利于强化传热,且在Ra数较小时,平均努塞尔（Nu）数增加的幅度要优于Ra数较大的情况。在所研究的控制参数范围内,方形热源的平均Nu数最大。根据数值结果给出不同热源表面的平均Nu数、纳米颗粒体积分数、Ra数三者之间的函数关系式,该函数关系可为此类工程的设计提供理论指导。     

二维半线性扩散反应方程的高精度全隐格式及其多重网格方法

针对二维非定常半线性扩散反应方程，空间导数项采用四阶紧致差分公式离散，时间导数项采用四阶向后Euler公式进行离散，提出一种无条件稳定的高精度五层全隐格式.格式截断误差为O（τ⁴+τ²h²+h⁴），即时间和空间均具有四阶精度.对于第一、二、三时间层采用Crank-Nicolson方法进行离散，并采用Richardson外推公式将启动层时间精度外推到四阶.建立适用于该格式的多重网格方法，加快在每个时间层上迭代求解代数方程组的收敛速度，提高计算效率.最后通过数值实验验证格式的精确性和稳定性以及多重网格方法的高效性.     

滑移边界对二维Rayleigh-Bénard热对流流态和传热特性的影响

采用并行直接数值模拟（PDM-DNS）计算无滑移和滑移边界二维Rayleigh-Bénard热对流.与无滑移边界形成的随机羽流运动的湍流热对流不同,滑移边界热对流最终形成湍流特征消失,且温度仅分布于四壁的一个大尺度环流的流动形态.平均场近底板的温度分布特性,无滑移边界逐渐变化而滑移边界出现过冲现象.宽高比Γ=1时,Nusselt数（Nu）随Rayleigh数（Ra）的变化具有相同标度指数,Nu~Ra^0.3.滑移边界热对流具有传热增强作用.滑移边界热对流Nu随Γ变化明显,并分为两个阶段,在Γ=0.5时出现Nu_max≈250,是无滑移边界热流Nu的5倍.     

三维定常问题的高阶紧致差分格式向非定常问题的推广

将已经建立的求解三维定常对流扩散方程的高阶紧致差分格式直接推广到三维非定常对流扩散方程的数值求解,时间导数项利用二阶向后欧拉差分公式,所得到的高阶隐式紧致差分格式时间为二阶精度,空间为四阶精度,并且是无条件稳定的.数值实验结果验证了本文方法的精确性和稳健性.     

摄动有限体积法重构近似高精度的意义

研讨有限体积(FV)方法重构近似高精度的作用问题.FV方法中积分近似采用中点规则为二阶精度时,重构近似高精度(精度高于二阶)的意义和作用是一个有争议的问题.利用数值摄动技术[1,2]构造了标量输运方程的积分近似为二阶精度、重构近似为任意阶精度的迎风型和中心型摄动有限体积(PFV)格式.迎风PFV格式无条件满足对流有界准则(CBC),中心型PFV格式为正型格式,两者均不会产生数值振荡解.利用PFV格式求解模型方程的数值结果表明:与一阶迎风和二阶中心格式相比,PFV格式精度高、对解的间断分辨率高、稳定性好、雷诺数的适用范围大,数值地"证实"重构近似高精度和PFV格式的实际意义和好处.     

基于树网格的格子Boltzmann方法以及曲线边界的处理

将格子玻尔兹曼方法(LBM)和BFECC技术相结合,在四叉树网格上实现,并用于二维不可压缩流动的求解.二维四叉树(三维八叉树)网格是一种非均匀网格,利用树结构进行数据存储,具有非常高的检索效率.BFECC是一种数值技巧,可以把求解输运方程的奇数阶格式通过误差修正提高一阶精度.基于四叉树网格的LBM,引入BFECC后,能有效减小输运过程中非均匀网格插值产生的误差,在总体上实现二阶精度,和经典的LBM方法相当.最后给出二维顶盖驱动方腔流动算例和二维圆柱绕流算例,从中可以看出格式的有效性.     

Poisson方程的高精度有限差分格式

本文提出数值求解Poisson方程的含选择因子α的预示校正差分格式,它具有四阶精度。第一种格式处理Drichlet边界条件的Poisson方程,它包括Bramble的差分格式和林群、吕涛提出的差分外推格式。第二种格式处理Neumann边界条件的Poisson方程。对于工程计算常用的粗糙网络,作者建议采用α≲0.5的预示校正差分格式。     

瞬态热传导问题的一阶对称SPH方法模拟

为提高传统光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法模拟瞬态热传导问题的精度和稳定性,本文提出了一种一阶对称光滑粒子动力学(first order symmetric SPH, FO-SSPH)方法.该方法将具有二阶热传导方程分解成两个一阶偏微分方程,然后基于梯度离散和Taylor级数展开思想,对一阶核梯度形式进行修正,并将得到的局部矩阵对称化.数值结果表明:与传统SPH方法相比,FO-SSPH方法精度高、数值稳定性好; 该方法能较准确地直接施加混合边值 关键词： 瞬态热传导 光滑粒子动力学 非线性     

过渡区微尺度流动的有效黏性多松弛系数格子Boltzmann模拟

为提高采用二维九速离散速度模型的格子Boltzmann方法 (LBM)模拟微尺度流动中非线性现象的精度和效率,引入Dongari等提出的有效平均分子自由程对黏性进行修正(Dongari N,Zhang Y H,Reese J M2011 J.Fluids Eng.133 071101);并针对以往研究微尺度流动时采用边界处理格式含有离散误差的问题,采用多松弛系数格子Boltzmann方法结合二阶滑移边界条件,对微尺度Couette流动和周期性Poiseuille流动进行模拟,并将速度分布以及质量流量等模拟结果与直接模拟蒙特卡罗方法模拟数据、线性Boltzmann方程的数值解以及现有的LBM模型模拟结果进行对比.结果表明,相对于现有的LBM模型,引入新的修正函数所建立的有效黏性多松弛系数LBM模型有效提高了LBM模拟过渡区的微尺度流动中的非线性现象的能力.     

Entropy generation minimization and statistical declaration with probable error for skin friction coefficient and Nusselt number

Main emphasis of present work is to analyze the novel feature of entropy generation in MHD nanomaterial flow between two rotating disks. Heat transfer process is explored in the presence of Joule heating and thermal radiation. Tiwari–Das nanofluid model is employed in mathematical modeling. Aluminum oxide and copper water nanoparticles are accounted. Statistical declaration and probable error for problem accuracy are computed. Total entropy generation subject to Bejan number is scrutinized. Suitable variables are utilized to transform nonlinear PDEs to ordinary ones. Convergent series solutions are computed. Zeroth and mth order problems are discussed for stability analysis. The impact of physical flow variables like Reynolds number, magnetic parameter, porosity parameter, stretching parameter, rotational parameter, radiation parameter, Eckert number, suction injection parameter, Brinkman number and temperature ratio parameter on velocities, temperature, total entropy generation and Bejan number are examined and discussed through graphs. Velocity and thermal gradients at the surface of disks are computed.     

基于对流传热模型的LD端面泵浦圆形激光晶体的热效应研究

通过对半导体端面泵浦棒状Nd∶YAG晶体的热效应进行了理论分析,研究了端面抽运圆形截面激光晶体内部温度场,建立了符合条件的激光晶体热模型。考虑晶体侧面与冷却液之间的对流传热,以及晶体端面与外界非绝热边界条件,从而建立更为合理的边界条件,得出更符合实际的晶体的温度分布场。研究结果表明,考虑端面的对流传热后,计算的晶体中心温度降低,而相应的热焦距稍有增加;空气传热系数增加时,晶体中心温度明显降低,热焦距显著增加,减弱了晶体热效应。     

General Solutions for Hydromagnetic Free Convection Flow over an Infinite Plate with Newtonian Heating,Mass Diffusion and Chemical Reaction

The problem of hydromagnetic free convection flow over a moving infinite vertical plate with Newtonian heating, mass diffusion and chemical reaction in the presence of a heat source is completely solved. Radiative and porous effects are not taken into consideration but they can be immediately included by a simple rescaling of Prandtl number and magnetic parameter. Exact general solutions for the dimensionless velocity and concentration fields and the corresponding Sherwood number and skin friction coefficient are determined under integral form in terms of error function or complementary error function of Gauss. They satisfy all imposed initial and boundary conditions and can generate exact solutions for any problem with technical relevance of this type. As an interesting completion, uncommon in the literature, the differential equations which describe the thermal, concentration and momentum boundary layer, as well as the exact expressions for the thicknesses of thermal, concentration or velocity boundary layers were determined.Numerical results have shown that the thermal boundary layer thickness decreases for increasing values of Prandtl number and the concentration boundary layer thickness is decreasing with Schmidt number. Finally, for illustration,three special cases are considered and the influence of physical parameters on some fundamental motions is graphically underlined and discussed. The required time to reach the flow according with post-transient solution(the steady-state),for cosine/sine oscillating concentrations on the boundary is graphically determined. It is found that, the presence of destructive chemical reaction improves this time for increasing values of chemical reaction parameter.     

An optimization method for the problems of thermal cloaking of material bodies

Inverse heat-transfer problems related to constructing special thermal devices such as cloaking shells, thermal-illusion or thermal-camouflage devices, and heat-flux concentrators are studied. The heatdiffusion equation with a variable heat-conductivity coefficient is used as the initial heat-transfer model. An optimization method is used to reduce the above inverse problems to the respective control problem. The solvability of the above control problem is proved, an optimality system that describes necessary extremum conditions is derived, and a numerical algorithm for solving the control problem is proposed.     

Multiphysical simulation analysis of the dislocation structure in germanium single crystals

To grow high-quality germanium crystals is one of the most important problems of growth industry. The dislocation density is an important parameter of the quality of single crystals. The dislocation densities in germanium crystals 100 mm in diameter, which have various shapes of the side surface and are grown by the Czochralski technique, are experimentally measured. The crystal growth is numerically simulated using heat-transfer and hydrodynamics models and the Alexander–Haasen dislocation model in terms of the CGSim software package. A comparison of the experimental and calculated dislocation densities shows that the dislocation model can be applied to study lattice defects in germanium crystals and to improve their quality.     

The Reynolds analogy and a new formulation of the temperature-defect law for a turbulent boundary layer on a plate

A rational asymptotic theory describing the dynamic and thermal turbulent boundary layer on a plate at zero pressure gradient is proposed. The fact that the flow depends on a finite number of governing parameters makes it possible to formulate algebraic closure conditions, which relate the turbulent shear stress and heat flux to mean velocity and temperature gradients. As a result of an exact asymptotic solution of the boundary-layer equations, the known laws of the wall for the velocity and temperature and the velocity and temperature defect laws as well as the expression for the skin-friction coefficient, the Stanton number, and the Reynolds-analogy factor are obtained. The latter implies two new formulations of the temperature-defect law, one of which is completely similar to the velocity-defect law and does not contain the Stanton number and the turbulent Prandtl number, and the other does not contain the skin-friction coefficient. A heat-transfer law that relates only thermal quantities is also obtained. The conclusions of the theory agree well with experimental data.     

Evaluation of the spatial resolution of passive acoustic thermal tomography

Possible spatial resolution of a passive acoustic thermal tomograph is evaluated in numerical experiments. A criterion for the evaluation of this parameter is proposed. The criterion is based on the formation of a dip between two reconstructed temperature peaks with increasing distance between the thermal sources in the region under investigation. The spatial resolution depends on the positions of a pair of temperature peaks relative to the body surface and on the algorithm used for the reconstruction of the temperature distribution. An algorithm version providing an adequate reconstruction of the heights of the distribution peaks is suggested. It takes into account a priori information on the ultrasonic absorption coefficient and also the characteristics of the heat-transfer processes and the presence of blood circulation in a human body. In the case of lateral positions of the pair of temperature peaks relative to the surface, the spatial resolution is ≈1.7 cm. The prospects for the improvement of the spatial resolution of a passive acoustic thermal tomographs are discussed.     

Stability of a spherical flame ball in a porous medium

Gaseous flame balls and their stability to symmetric disturbances are studied numerically and asymptotically, for large activation temperature, within a porous medium that serves only to exchange heat with the gas. Heat losses to a distant ambient environment, affecting only the gas, are taken to be radiative in nature and are represented using two alternative models. One of these treats the heat loss as being constant in the burnt gases and linearizes the radiative law in the unburnt gas (as has been studied elsewhere without the presence of a solid). The other does not distinguish between burnt and unburnt gas and is a continuous dimensionless form of Stefan's law, having a linear part that dominates close to ambient temperatures and a fourth power that dominates at higher temperatures.

Numerical results are found to require unusually large activation temperatures in order to approach the asymptotic results. The latter involve two branches of solution, a smaller and a larger flame ball, provided heat losses are not too high. The two radiative heat loss models give completely analogous steady asymptotic solutions, to leading order, that are also unaffected by the presence of the solid which therefore only influences their stability. For moderate values of the dimensionless heat-transfer time between the solid and gas all flame balls are unstable for Lewis numbers greater than unity. At Lewis numbers less than unity, part of the branch of larger flame balls becomes stable, solutions with the continuous radiative law being stable over a narrower range of parameters. In both cases, for moderate heat-transfer times, the stable region is increased by the heat capacity of the solid in a way that amounts, simply, to decreasing an effective Lewis number for determining stability, just as if the heat-transfer time was zero.     

Thermal explosion in a stirred medium

The critical conditions of a thermal explosion in a cylindrical chemical reactor with several symmetrically positioned agitators for creating a forced convection of the reaction mixture are examined. The analysis is carried out for an arbitrary Peclet number on the assumption that the flow of the fluid is laminar. The critical value of the thermal explosion parameter (Frank-Kamenetskii parameter) is determined as a function of the number of agitators, Peclet number, and the distances between the axes of the agitators and the reactor axis. In contrast to the previous studies, the dependence of the thermal explosion parameter on the rate of stirring of the medium is considered. Calculations show that, as the number of agitators and mixing rate increases, the thermal explosion probability decreases, so that the critical value of the thermal explosion parameter can be several times greater than its classical value. It is demonstrated that, depending on the values of the parameters, the problem can have four stationary solutions, two of which are stable.     

真空下气固界面的传热

已知对于真空下的传热稀薄气体起着非常重要的作用。该文给出了稀薄气体在不同克努曾数区域下的气 -固界面间的传热表达式 ,这些表达式中包含了气体反射系数和热适应系数。文中讨论了平行平板、同轴圆柱表面和同心圆球表面的情况 ,给出了不同传热区域热流依照克努曾数倒数关系变化的数值计算结果