A new complex variable meshless method for transient heat conduction problems

<正>In this paper,based on the improved complex variable moving least-square(ICVMLS) approximation,a new complex variable meshless method(CVMM) for two-dimensional(2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations,and the essential boundary conditions are imposed by the penalty method.As the transient heat conduction problems are related to time,the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization.Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained.In order to demonstrate the applicability of the proposed method,numerical examples are given to show the high convergence rate,good accuracy,and high efficiency of the CVMM presented in this paper.     

A complex variable meshless local Petrov-Galerkin method for transient heat conduction problems

On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.     

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

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

An improved local radial point interpolation method for transient heat conduction analysis

The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.     

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.     

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

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

FGMs稳态热传导分析的重心Lagrange插值配点法

提出一种求解二维功能梯度材料（FGMs）稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明：该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析.     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

亚表面异质缺陷对功能梯度材料表面温度场的影响

基于双曲型热传导方程,采用镜像法和波函数展开法,求解了含亚表面异质圆柱缺陷的半无限功能梯度材料的表面温度场,给出了功能梯度材料中热波散射的一般解.分析了亚表面异质圆柱缺陷的几何参数(如埋藏深度)和热物理参数(如导热系数、热扩散长度、热扩散率及热弛豫时间等)对功能梯度材料表面温度场的影响.温度波由调制的超短脉冲激光在功能梯度材料表面激发,异质圆柱缺陷表面的边界条件为导热边界.研究结果可望为功能梯度材料的红外热波无损检测、导热反问题提供计算方法和参考数据.     

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

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

A general modelling method for functionally graded materials with an arbitrarily oriented crack

In order to analytically solve crack problems regarding functionally graded materials (FGMs), some ideal assumptions are often made. They are: (1) the properties of FGMs are usually assumed to be described by very particular functions; (2) the crack is assumed to be vertical to (or parallel to) the gradient of FGMs. However, these assumptions may not be practical for actual FGMs. Since the controlling differential equations with general mechanical properties are very difficult to solve and the arbitrarily oriented crack causes great trouble in the analytical procedure, a general piecewise-exponential model (GPE model) is proposed to investigate the fracture behaviour of FGMs with general mechanical properties and an arbitrarily oriented crack. “General mechanical properties” means that the mechanical properties in the GPE model are not required to be very particularly pre-defined functions but arbitrary functions determined by fitting the experimental results of FGMs. The studied FGMs are divided into some sub-layers with each layer’s properties varying exponentially so that the general mechanical properties can be approximated by a series of exponential functions and hence the stresses and displacements of each layer which may contain a mixed-mode crack can be solved analytically. By use of integral transform methods, principle of superposition, residual theorem and theory of singular integral equations, the mixed-mode crack problem can be turned into solving a group of singular integral equations from which mixed-mode stress intensity factors (SIFs) can be obtained. Finally, the influences of the nonhomogeneous and geometric parameters on the mixed-mode SIFs are analysed.     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

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.     

A model of heat conduction

We construct a model of a chain of atoms coupled at its ends to two reservoirs at different temperatures. In a weak coupling limit the atoms obey a stochastic evolution law and have an equilibrium state with a uniform temperature gradient along the chain.     

On the entropy production for the Davies model of heat conduction

The formula for the entropy production in open quantum systems is examined for the Davies model of heat conduction.This work is supported by Polish Ministry of Higher Education Science and Technology, project MRI 7.     

导热几何形状反演的变分原理及边界条件的确立

本文采用变域变分原理,建立了导热几何形状反演问题的变分原理,同时获得了该问题所需满足的边界条件和附加条件.该变分原理能将未知形状的几何变量及控制方程结合在一个变分泛函中,使得数学描述简洁、紧凑,且几何变量及控制方程的求解能耦合地进行.介绍了运用该变分原理并结合有限元方法进行数值计算的方法. 关键词： 几何形状反演 变分原理 有限元 导热