瞬态热传导问题的复变量重构核粒子法

在重构核粒子法的基础上,引入复变量,讨论了复变量重构核粒子法.复变量重构核粒子法的优点是在构造形函数时采用一维基函数建立二维问题的修正函数.然后,将复变量重构核粒子法应用于瞬态热传导问题的求解,结合瞬态热传导问题的Galerkin积分弱形式,采用罚函数法引入本质边界条件,建立了瞬态热传导问题的复变量重构核粒子法,推导了相应的计算公式.与传统的重构核粒子法相比,复变量重构核粒子法具有计算量小、精度高的优点.最后通过数值算例证明了该方法的有效性. 关键词： 重构核粒子法 复变量重构核粒子法 修正函数 瞬态热传导问题     

光滑粒子动力学SPH方法应力不稳定性的一种改进方案

光滑粒子动力学方法是一种拉格朗日型无网格粒子方法,在模拟大变形和自由表面流方面具有特殊的优势,已经在工程和科学领域得到了广泛的应用．然而,长期以来,传统光滑粒子动力学方法一直受到应力不稳定性的困扰,从而限制了它的进一步发展和应用．应力不稳定性的根本原因在于应力状态与核函数的不匹配：负压状态下粒子间产生吸引力,吸引力随着粒子间距的减小而增大,导致拉伸不稳定性;正压状态下粒子间产生排斥力,排斥力随着粒子间距的减小而先增大后减小,导致压缩不稳定性．本文通过改进光滑粒子动力学方法的核函数和离散格式,使得无论在正压还是负压状态下粒子间的作用力恒为排斥力,且排斥力随着粒子间距的减小而增大,从而防止粒子聚集等现象,解决应力不稳定问题．分别使用改进前后的光滑粒子动力学方法模拟两个典型的应力不稳定算例,结果表明本文的改进方法能够有效地消除应力不稳定性．     

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.     

气粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟

在气粒两相流动问题中,颗粒间以及气体与颗粒间的传热问题不可忽略.光滑离散颗粒流体动力学(SDPH)模型作为一种新的求解气粒两相流动问题的方法,已经成功应用于模拟风沙运动等问题.在此基础上,提出了SDPH方法的热传导模型,模拟了气粒两相流动问题中的热传导过程以及颗粒蒸发过程.首先引入各相的能量方程,利用有限差分与光滑粒子流体动力学一阶导数相结合的方法,处理各相内部热传导项中的二阶导数问题,基于气粒两相间温度差及对流换热系数计算颗粒与气体间的热传导量,推导得到了含热传导模型的气粒两相流SDPH计算方程组,模拟计算了圆盘形颗粒团算例及鼓泡流化床内部热传导算例,并与双流体模型计算结果进行对比,结果基本符合;其次利用离散液滴模型中的颗粒蒸发传质传热定律计算颗粒的蒸发过程,数值模拟了颗粒射流蒸发过程,并与离散颗粒模型结果进行对比,两者符合得较好,验证了该方法的准确性及实用性.     

气体热导率实验中的传热机制

分析了从低气压到常压的实验数据,并进行了理论计算,得出了实验中气体传热机制的演变过程.     

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比

为了对比研究弱可压光滑粒子动力学(WCSPH)方法和不可压光滑粒子动 力学(ISPH)方法在模拟封闭方腔自然对流问题时的特性, 采用粒子位移技术有效地解决了高瑞利数条件下, 拉格朗日型SPH方法模拟封闭方腔自然对流时流体域内的粒子聚集和空穴问题, 将拉格朗日型SPH 方法求解封闭方腔自然对流问题的最高瑞利数提高到了10⁶; 进而通过对比瑞利数分别为10⁴, 10⁵, 10⁶的条件下, 采用拉格朗日型WCSPH、 拉格朗日型ISPH、欧拉型ISPH三种SPH方法模拟得到的封闭方腔速度分布云图、 温度分布云图、壁面努赛尔特数分布曲线和平均努塞尔特数, 分析了三种SPH方法在模拟封闭方腔自然对流时的精度、稳定性和计算效率. 结果表明: 在低瑞利数条件下, 以上三种SPH方法都可以较好地模拟此问题, 在高瑞利数条件下, 欧拉型ISPH方法的模拟结果最为精确; 拉格朗日型WCSPH方法模拟所得结果比拉格朗日型ISPH方法模拟所得结果稍好些. 关键词： 光滑粒子动力学 不可压光滑粒子动力学 粒子位移技术 自然对流     

非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟

为准确、有效地模拟非等温非牛顿黏性流体的流动问题,本文基于一种不含核导数计算的核梯度修正格式和不可压缩条件给出了一种改进光滑粒子动力学(SPH)离散格式,它较传统SPH离散格式具有较高精度和较好稳定性.同时,为准确地描述温度场的演化过程,建立了非牛顿黏性的SPH温度离散模型.通过对等温Poiseuille流、喷射流和非等温Couette流、4:1收缩流进行模拟,并与其他数值结果作对比,分别验证了改进SPH方法模拟非牛顿黏性流动问题的可靠性和提出的SPH温度离散模型求解非等温流动问题的有效性和准确性.随后,运用改进SPH方法结合SPH温度离散模型对环形腔和C形腔内非等温非牛顿黏性流体的充模过程进行了试探性模拟研究,分析了数值模拟的收敛性,讨论了不同位置处热流参数对温度和流动的影响.     

基于修正并行光滑粒子动力学方法三维变系数瞬态热传导问题的模拟

为了解决传统光滑粒子动力学(SPH)方法求解三维变系数瞬态热传导方程时出现的精度低、稳定性差和计算效率低的问题,本文首先基于Taylor展开思想拓展一阶对称SPH方法到三维热传导问题的模拟,其次引入稳定化处理的迎风思想,最后基于相邻粒子标记和MPI并行技术,结合边界处理方法得到一种能够准确、高效地求解三维变系数瞬态热传导问题的修正并行SPH方法.通过对带有Direclet和Newmann边界条件的常/变系数三维热传导方程进行模拟,并与解析解进行对比,对提出的方法的精度、收敛性及计算效率进行了分析;随后,运用提出的修正并行SPH方法对三维功能梯度材料中温度变化进行了模拟预测,并与其他数值结果做对比,准确地展现了功能梯度材料中温度变化过程.     

燃料抛撒成雾及其燃烧爆炸的光滑离散颗粒流体动力学方法数值模拟研究

燃料在炸药爆炸驱动下形成燃料空气爆炸云团, 进而引燃爆炸, 对目标造成毁伤. 本文在前期提出的光滑离散颗粒流体动力学方法(SDPH)的基础上, 引入描述炸药由爆轰到膨胀整个过程的Jones-Wilkins-Lee状态方程及描述气体快速燃烧过程的EBU-Arrhenius燃烧模型, 建立了求解战斗部起爆、燃料抛撒和燃料二次引燃爆炸问题的新型SDPH方法. 设计了圆环形燃料颗粒在炸药爆炸驱动下运动抛撒的算例进行数值验证, 结果与理论相符; 对燃料空气炸药(FAE)云雾的形成和发展过程进行了数值模拟, 分析了云雾的形态, 并与实验结果进行对比, 符合较好, 同时分析了不同起爆方式对云雾团成型的影响; 最后, 在云雾团成型的基础上, 引入蒸发燃烧模型对FAE的燃烧爆炸过程进行了模拟研究. 结果表明, 本文建立的数学模型和计算方法可以较好的模拟燃料空气炸药抛撒成雾及云雾燃烧爆炸过程, 为该类武器装备的设计研究提供了较好的数值方法.     

The shape gradient of the least-squares objective functional in optimal shape design problems of radiative heat transfer

Optimal shape design problems of steady-state radiative heat transfer are considered. The optimal shape design problem (in the three-dimensional space) is formulated as an inverse one, i.e., in the form of an operator equation of the first kind with respect to a surface to be optimized. The operator equation is reduced to a minimization problem via a least-squares objective functional. The minimization problem has to be solved numerically. Gradient minimization methods need the gradient of a functional to be minimized. In this paper the shape gradient of the least-squares objective functional is derived with the help of the shape sensitivity analysis and adjoint problem method. In practice a surface to be optimized may be (or, most likely, is to be) given in a parametric form by a finite number of parameters. In this case the objective functional is, in fact, a function in a finite-dimensional space and the shape gradient becomes an ordinary gradient. The gradient of the objective functional, in the case that the surface to be optimized is given in a finite-parametric form, is derived from the shape gradient. A particular case, that a surface to be optimized is a “two-dimensional” polyhedral one, is considered. The technique, developed in the paper, is applied to a synthetic problem of designing a “two-dimensional” radiant enclosure.     

弹性壳结构静力与动力分析的光滑粒子法

本文通过采用移动最小二乘函数作为近似函数 和完全拉格朗日方程作为近似方程来改善光滑粒子法的稳定性和数值精度; 在此基础上, 提出了壳结构静力分析的光滑粒子法, 并完善了壳结构动力分析方法; 最后, 采用国际公认的壳结构的标准测试模型对静力和动力问题分别进行了验证, 所得结果与已有数据吻合良好, 证明了本文数值模型的有效性和可靠性, 为光滑粒子法进一步在裂纹、破碎等非线性壳结构中的应用提供参考. 关键词： 弹性壳 静力与动力分析 光滑粒子法 完备性和稳定性     

Analysis of frost visualization over a fin and tube heat exchanger by natural convection

In this article, studies are made on frost formation and flow over a fin and tube heat exchanger due to natural convection for various conditions of relative humidity, air ambient temperature, and mean refrigerant temperature. The results include frost deposition, steps of frost formation, and its effect on heat transfer rate for different conditions. The results show that frost is formed only on the tip of the fins with higher thickness from top to bottom due to small distance between the fins. Frost causes to trap the air which increases the thermal resistance and reduces heat transfer in the system.     

Stability of cylindrical conducting fluids with heat and mass transfer in longitudinal periodic electric field

A novel system to study the effect of an axial periodic electric field on the stability of a system of cylinders of conducting fluids in the presence of heat and mass transfer is investigated. The stability of a cylindrical interface between the vapor and liquid phases of a fluid is studied when the vapor is hotter than the liquid and the two phases are enclosed between two cylindrical surfaces coaxial with the interface. The linear dispersion relation is found to be of damped Mathieu-type equation with real coefficients. The method of multiple time scales is used to obtain approximate solution and analyze the stability criteria for both the nonresonant and resonant cases. The stability of the system is also discussed analytically and numerically for such cases. It is found that both the heat and mass transfer coefficient and the dimensions of the system have destabilizing influences on the considered system, while azimuthal wavenumber is found to have a stabilizing effect. The dual role of the electric field frequency is also observed on the stability of the system depending on the electrical conductivities values. Finally the behaviour of the resonance points corresponding to the effects of each of the above physical parameters are determined, and a comparison between the obtained results with the corresponding results in the case of a constant applied electric field is achieved.     

竖直恒温面自然对流边界层不稳定性研究

在无扰动、随机式扰动以及正弦式扰动下,通过对竖直恒温面处状态Ra为1.328×10^9、Pr为6.24的自然对流进行模拟,探索了热边界层的不稳定性和共振强化自然对流换热。结果表明:(1)竖直自然对流边界层上游位置的随机式扰动对热边界层的影响主要体现在稳定阶段;(2)该状态下的竖直自然对流边界层的特征频率为15 067,且相比于无扰动状态,频率为15 067的正弦式扰动能在竖直恒温面处提高5.15%的换热量;(3)在竖直自然对流边界层上游位置加入特征频率的正弦式扰动,竖直恒温面处的局部努塞尔数Nu均出现明显波动,且波动随着边界层高度的增加而增大。     

Soret and Dufour effects on MHD non-Darcian mixed convection heat and mass transfer over a stretching sheet with non-uniform heat source/sink

An analysis has been carried out to study the effects of thermal-diffusion and diffusion-thermo on non-Darcian mixed convection heat and mass transfer of an incompressible, electrically conducting fluid over a stretching sheet embedded in a porous medium in the presence of an external magnetic field and non-uniform heat source/sink. Similarity transformations are used to convert highly non-linear partial differential equations into ordinary differential equations. Similarity equations are then solved numerically using shooting algorithm with Runge-Kutta-Fehlberg scheme over the entire range of physical parameters. The effects of various physical parameters on the dimensionless velocity, temperature and concentration profiles are depicted graphically. Present results are compared with previously published work on various special cases of the problem and the results are found to be in very good agreement. Numerical results for local skin-friction, local Nusselt number and local Sherwood number are tabulated for different physical parameters.     

Radiative and conductive heat transfer in a nongrey semitransparent medium. Application to fire protection curtains

This paper deals with heat transfer in nongrey media which scatter, absorb and emit radiation. Considering a two dimensional geometry, radiative and conductive phenomena through the medium have been taken into account. The radiative part of the problem was solved using the discrete ordinate method with classical S_n quadratures. The absorption and scattering coefficients involved in the radiative transfer equation (RTE) were obtained from the Mie theory. Conduction inside the medium was linked to the RTE through the energy conservation. Validation of the model has been achieved with several simulation of water spray curtains used as fire protection walls.     

Interaction between coherent structures and surface temperature and its effect on ground heat flux in an unstably stratified boundary layer

Surface layer plumes, thermals, downdrafts and roll vortices are the most prominent coherent structures in an unstably stratified boundary layer. They contribute most of the temperature and vertical velocity variance, and their time scales increase with height. The effects of these multi-scale structures (surface layer plumes scale with surface layer depth, thermals scale with boundary layer height and the resulting roll vortices scale with convective time scale) on the surface temperature and ground heat flux were studied using turbulence measurements throughout the atmospheric boundary layer and the surface temperature measurements from an infrared camera. Plumes and thermals imprint on the surface temperature as warm structures and downdrafts imprint as cold structures. The air temperature trace shows a ramp-like pattern, with small ramps overlaid on a large ramp very close to the surface; on the other hand, surface temperature gradually increases and decreases. Turbulent heat flux and ground heat flux show similar patterns, with the former lagging the latter. The maximum values of turbulent heat flux and ground heat flux are 4 and 1.2 times the respective mean values during the ejection event. Surface temperature fluctuations follow a similar power-law exponent relationship with stability as suggested by surface layer similarity theory.     

Transient radiation and conduction heat transfer inside a plane-parallel participating gray medium with boundaries having different reflecting characteristics

A hybrid ray-tracing method is developed for the solution to the radiative transfer in a plane-parallel participating medium having one specular surface and another diffuse surface. By this method, radiative transfer coefficients (RTCs) for specular–diffuse (S–D) surfaces are deduced. The medium surfaces are considered to be semitransparent. The effects of convection–radiation parameter, conduction–radiation parameter and refractive index on transient coupled heat transfer are investigated. Results show that the temperature curves of the medium having S–D surfaces is higher than those of the medium having S–S surfaces (two specular surfaces); the total heat flux at steady state for the S–D surfaces is lower than that for the S–S surfaces.