二维瞬态热传导问题的无单元Galerkin法分析

采用无单元Galerkin(element-free Galerkin,EFG)法求解具有混合边界条件的二维瞬态热传导问题.首先采用二阶向后微分公式离散热传导方程的时间变量,将该问题转化为与时间无关的混合边值问题;然后采用罚函数法处理Dirichlet边界条件,建立了二维瞬态热传导问题的无单元Galerkin法;最后基于移动最小二乘近似的误差结果,详细推导了无单元Galerkin法求解二维瞬态热传导问题的误差估计公式.给出的数值算例表明计算结果与解析解或已有数值解吻合较好,该方法具有较高的计算精度和较好的收敛性.     

An Element-Free Galerkin Method for Time-Fractional Diffusion-Wave Equations北大核心CSCD

利用无单元Galerkin法,对Caputo意义下的时间分数阶扩散波方程进行了数值求解和相应误差理论分析。首先用L1逼近公式离散该方程中的时间变量,将时间分数阶扩散波方程转化成与时间无关的整数阶微分方程;然后采用罚函数方法处理Dirichlet边界条件,并利用无单元Galerkin法离散整数阶微分方程;最后推导该方程无单元Galerkin法的误差估计公式。数值算例证明了该方法的精度和效果。     

对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

一类三维逆时热传导问题的数值求解

该文考虑一类三维逆时热传导问题的数值解法.基于有限差分时间离散,并结合伽辽金(Galerkin)方法对空间进行有限元离散,导出刚度矩阵及载荷向量,对热传导问题进行数值求解.针对反问题,利用分离变量法建立T时刻温度场与初始温度场之间的对应关系,给出了反演公式,并在一定先验假设条件下证明了反问题的局部稳定性.为克服反问题求...     

三维热传导方程的高精度有限差分方法

提出了数值求解三维热传导方程的一个四阶精度的有限差分格式,首先对三个空间方向上的二阶导数项,采用四次样条函数来近似,从而得到半离散的常微分方程.然后利用常微分方程的解析解表达式,时间矩阵利用Padé近似,得到时间和空间均为四阶精度的差分格式.最后利用方法计算了两个数值算例,并与文献中结果进行了对比,从而验证了高精度格式的性能.     

有限差分-配点法求解二维Burgers方程

将重心插值配点法结合Crank-Nicolson差分格式来求解Burgers方程.首先,利用Hopf-Cole变换将Burgers方程转化为线性热传导方程;空间方向采用重心插值配点法进行离散,时间方向采用Crank-Nicolson格式离散,导出对应的线性代数方程组,并对此计算格式进行相容性分析;最后,通过数值算例验证此计算格式具有高精度和有效性.     

非定常 N-S 方程的非线性 Galerkin 算法及其误差估计

本文给出了二维非定常N-S方程的三种数值格式,其中空间变量用谱非线性Galerkin算法进行离散,时间变量用有限差分离散,并研究了这些格式数值解的逼近精度．最后,给出了部分数值计算结果．     

求解一维侧边热传导方程反问题的正则化方法

研究了一维侧边热传导方程反问题.在求解一维侧边热传导方程的基础上,利用数值积分法进行离散化处理,然后引入正则化方法,采用偏差原理确定正则化参数,从而得到一维侧边热传导方程反问题的数值解.数值模拟结果表明,给出的正则化方法对于求解一维侧边热传导方程反问题是可行有效的.     

边界节点法计算二维瞬态热传导问题

采用边界节点法(BKM)结合双重互易法(DRM)求解二维瞬态热传导问题.采用差分格式处理时间变量,可将原瞬态热传导方程转化为一系列非齐次修正的Helmholtz方程.随后,方程的解可分为特解和齐次解两部分计算,引入双重互易法在区域内部配点求解方程的特解,采用边界节点法仅需边界配点求解方程的齐次解.给出的数值算例显示该方法计算精度高,适用性好,具有很好的稳定性和收敛性,适合求解瞬态热传导问题.     

无限域上一维热传导方程的Gauss-Laguerre数值解

针对无限域上一维热传导方程的解析解为反常积分形式,直接计算往往比较困难.首先采用Fourier变换给出问题解析解,其次结合解析解的形式和无限域上Gauss型数值积分法精度高的优点,将半无限域上的一维热传导方程问题利用Gauss-Laguerre数值积分计算数值解,对无限域上的一维热传导方程的解析解转化为半无限域上的形式后用Gauss-Laguerre数值积分计算.实验结果表明,本文给出的数值解方法具有很高的精度.     

A CQM-based BEM for transient heat conduction problems in homogeneous materials and FGMs

A method for numerical solution of time-domain boundary integral formulations of transient problems governed by the heat equation is presented. The heat conduction problem is analyzed considering homogeneous and non-homogeneous media. In the case of the non-homogeneous media, the conductor material is assumed to be a functionally graded material, i.e., the material properties vary spatially according to known smooth functions. For some specific spatial variations of the material properties, the fundamental solution and the boundary integral equation of the problem are obtained thanks to a change of variables that transforms the original problem to the standard heat conduction problem for homogeneous materials. For the treatment of time-dependent terms, the convolution quadrature method is adopted to approximate numerically the integral equation of the time-domain boundary element method. In the case that the responses are required at a large number of interior points, the convolution performed to calculate them is very time consuming. It is shown that the discrete convolution of the proposed formulation can be computed by means of the fast Fourier transform technique, which considerably reduces the computational complexity. Results for some transient heat conduction examples are presented to validate the numerical techniques studied.     

Study of convergence of the projection-difference method for hyperbolic equations

We consider the Cauchy problem for an abstract quasilinear hyperbolic equation with variable operator coefficients and a nonsmooth but Bochner integrable free term in a Hilbert space. Under study is the scheme for approximate solution of this problem which is a combination of the Galerkin scheme in space variables and the three-layer difference scheme with time weights. We establish an a priori energy error estimate without any special conditions on the projection subspaces. We give a concrete form of this estimate in the case when discretization in the space variables is carried out by the finite element method (for a partial differential equation) and by the Galerkin method in Mikhlin form.     

The Method of Fundamental Solution for a Radially Symmetric Heat Conduction Problem with Variable Coefficient

We consider an inverse heat conduction problem with variable coefficient on an annulus domain. In many practice applications, we cannot know the initial temperature during heat process, therefore we consider a non-characteristic Cauchy problem for the heat equation. The method of fundamental solutions is applied to solve this problem. Due to ill-posedness of this problem, we first discretize the problem and then regularize it in the form of discrete equation. Numerical tests are conducted for showing the effectiveness of the proposed method.     

Multiscale computation method for parabolic problems of composite materials

In this paper, we consider the initial-boundary value problem of parabolic type equation with rapidly oscillating coefficients in both time and space. A multiscale asymptotic expansion of solution for this kind of problem is presented. The full discrete finite element method for computing above problem is introduced. This method can apply to heat conduction analysis of composite materials. The main advantages of this method are that it can greatly save computer memory and CPU time, and it has good precision at the same time. Finally numerical results show that the method presented in this paper is effective and reliable.     

Linear barycentric rational collocation method for solving heat conduction equation

The linear barycentric rational collocation method for solving heat conduction equation is presented. The matrix form of discrete heat conduction equation by collocation method is also obtained. With the help of convergence rate of the barycentric interpolation, the convergence rate of linear barycentric rational collocation method for solving heat conduction equation is proved. At last, several numerical examples are provided to validate the theoretical analysis.     

Numerical analysis of a continuous Galerkin method for damped sine-Gordon equation

In this article, we discuss the numerical solution for the two-dimensional (2-D) damped sine-Gordon equation by using a space–time continuous Galerkin method. This method allows variable time steps and space mesh structures and its discrete scheme has good stability which are necessary for adaptive computations on unstructured grids. Meanwhile, it can easily get the higher-order accuracy in both space and time directions. The existence and uniqueness to the numerical solution are strictly proved and a priori error estimate in maximum-norm is given without any space–time grid conditions attached. Also, we prove that if the mesh in each time level is generated in a reasonable way, we can get the optimal order of convergence in both temporal and spatial variables. Finally, the convergence rates are presented and analyzed by some numerical experiments to illustrate the validity of the scheme.     

Poisson方程的混合无网格方法(英文)

以Poisson方程的混合变分形式为基础,采用移动最小二乘方法建立插值形函数空间,给出了Poisson方程的混合无网格方法,理论上证明了Poisson方程混合无网格解的存在唯一性,并给出了误差估计.本质边界条件的处理采用Lagrange乘子法.数值算例表明,在应用相同阶次的基函数条件下,利用混合无网格方法求解Poisson方程所得的解的梯度值优于传统的无网格方法及有限元法.     

三维热传导型半导体问题的特征有限元方法和分析

本文研究三维热传导型半导体瞬态问题的特征有限元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出Ｇａｌｅｒｋｉｎ逼近;对电子,空穴浓度方程采用特征有限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近．应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。