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

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

The complex variable reproducing kernel particle method for two-dimensional elastodynamics

On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.     

Analysis of variable coefficient advection-diffusion problems via complex variable reproducing kernel particle method

The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method （RKPM） and the element- free Galerkin （EFG） method.     

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.     

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.     

弹性力学的复变量重构核粒子法

在重构核粒子法的基础上,提出了复变量重构核粒子法.复变量重构核粒子法的优点是采用一维基函数建立二维问题的修正函数.然后,将复变量重构核粒子法应用于弹性力学,提出了弹性力学的复变量重构核粒子法,并推导了相关公式.与传统的重构核粒子法相比,复变量重构核粒子法具有计算量小、效率高的优点.最后给出了数值算例证明了该方法的有效性. 关键词： 重构核粒子法 复变量重构核粒子法 弹性力学 无网格方法     

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.     

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

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

弹性力学的插值型重构核粒子法

采用具有离散点插值特性的重构核粒子法形函数, 较精确地重构弹性体 变形的位移试函数, 再与弹性力学的最小势能原理相结合, 形成新的分析弹性力 学平面问题的插值型重构核粒子法. 由于插值型重构核粒子法形函数具有点插值特性和不低于核函数 的高阶光滑性, 因而既克服了多数无网格方法处理本质边界条件的困难, 也保证了较高的数值精度. 与早期的无网格方法相比, 本方法具有精度高、解题规模较小、可直接施加边界条件等优点. 通过对典型弹性力学问题数值模拟, 验证了所提方法的有效性和正确性.     

An interpolating reproducing kernel particle method for two-dimensional scatter points

An interpolating reproducing kernel particle method for two-dimensional(2D) scatter points is introduced. It eliminates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating reproducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.     

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

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

弹性力学的重构核粒子边界无单元法

将重构核粒子法(RKPM)和边界积分方程方法结合，提出了一种新的边界积分方程无网格方法——重构核粒子边界无单元法(RKP-BEFM).对弹性力学问题，推导了其重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核粒子边界无单元法离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位移和应力积分公式.重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点的精确值，所以本方法具有更高的精度.最后给出了数值算例，验证了本方法的有效性和正确性. 关键词： 重构核粒子法 弹性力学 边界无单元法     

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

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

Improved kernel gradient free-smoothed particle hydrodynamics and its applications to heat transfer problems

Kernel gradient free-smoothed particle hydrodynamics(KGF-SPH) is a modified smoothed particle hydrodynamics(SPH) method which has higher precision than the conventional SPH.However,the Laplacian in KGF-SPH is approximated by the two-pass model which increases computational cost.A new kind of discretization scheme for the Laplacian is proposed in this paper,then a method with higher precision and better stability,called Improved KGF-SPH,is developed by modifying KGF-SPH with this new Laplacian model.One-dimensional(1D) and two-dimensional(2D) heat conduction problems are used to test the precision and stability of the Improved KGF-SPH.The numerical results demonstrate that the Improved KGF-SPH is more accurate than SPH,and stabler than KGF-SPH.Natural convection in a closed square cavity at different Rayleigh numbers are modeled by the Improved KGF-SPH with shifting particle position,and the Improved KGF-SPH results are presented in comparison with those of SPH and finite volume method(FVM).The numerical results demonstrate that the Improved KGF-SPH is a more accurate method to study and model the heat transfer problems.     

The complex variable meshless local Petrov—Galerkin method of solving two-dimensional potential problems

Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square(MLS) approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin(MLPG) method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.     

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.     

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

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

New complex variable meshless method for advection-diffusion problems

In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.     

黏弹性问题的改进的复变量无单元Galerkin方法

基于改进的复变量移动最小二乘法,提出了二维黏弹性问题的改进的复变量无单元Galerkin方法.采用改进的复变量移动最小二乘法建立形函数,根据Galerkin积分弱形式建立求解方程,并用罚函数法施加本质边界条件,推导了二维黏弹性问题的改进的复变量无单元Galerkin方法的计算公式.最后,通过实际算例,将计算结果与复变量无单元Galerkin方法及有限元法的结果进行了对比,说明了本文方法具有更高的计算精度和计算效率.