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

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

Poisson方程热源识别反问题

探讨了半带状区域上二维Poisson方程只含有一个空间变量的热源识别反问题.这类问题是不适定的,即问题的解(如果存在的话)不连续依赖于测量数据.利用Carasso-Tikhonov正则化方法,得到了问题的一个正则近似解,并且给出了正则解和精确解之间具有Holder型误差估计.数值实验表明Carasso-Tikhonov正则化方法对于这种热源识别是非常有效的.     

求解二维不可压缩Navier-Stokes方程的混合型微分求积法

微分求积法（ＤＱＭ）能以较少的网格点求得微分方程的高精度数值解,但采用单纯的微分求积法求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程时,只能对低雷诺数流动获得较好的数值解,当雷诺数较高时会导致数值解不收敛·　为此,提出了一种微分求积法与迎风差分法混合求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程的预估＿校正数值格式,用伪时间相关算法以较少的网格点获得了较高雷诺数流动的数值解·　作为算例,对１∶１ 和１∶２ 驱动方腔内的流动进行了计算,得到了较好的数值结果·     

一类热传导方程的反问题

在本篇文章中,主要研究的是用伴随问题方法解决热传导方程反问题中的系数识别问题。     

基于遗传算法的最佳摄动量法在反问题中的应用

抛物型方程参数的反演在工程中有重要的应用价值,针对最佳摄动量法对初始模型依赖性的严重不足,给出遗传算法对最佳摄动量法的改进的新算法并进行了数值模拟．从模拟结果可以看到模拟结果与估算值曲线图基本吻合,体现了该方法的有效性和高精度性,且新算法弥补了对初始数据严重依赖的不足以及保证了区域收敛的全局性．     

求解全局优化问题的三角进化算法

This paper presents a hybrid heuristic-triangle evolution （TE） for global optimization. It is a real coded evolutionary algorithm. As in differential evolution （DE）, TE targets each individual in current population and attempts to replace it by a new better individual. However, the way of generating new individuals is different. TE generates new individuals in a Nelder- Mead way, while the simplices used in TE is 1 or 2 dimensional. The proposed algorithm is very easy to use and efficient for global optimization problems with continuous variables. Moreover, it requires only one （explicit） control parameter. Numerical results show that the new algorithm is comparable with DE for low dimensional problems but it outperforms DE for high dimensional problems.     

改进的最佳摄动量法及在反问题中的应用

给出了改进的最佳摄动量法,并应用在双曲型方程参数反演问题的求解中.由遗传算法借助交叉和变异算子控制全局搜索来获得参数的初始迭代值,代入最佳摄动量法求解出稳定的高精度数值解.     

求解二维Poisson方程的重心插值配点法

采用重心Lagrange插值配点法计算了二维Poisson方程.采用重心Lagrange插值法构造近似函数,由配点法离散Poisson方程及其边界条件.数值算例表明方法具有理论简单、计算精度高的特点.     

求解二维时谐Maxwell方程的一种混合有限元新格式

本文,我们提出一种新的求解二维时谐Maxwell方程的H~1-协调节点连续混合有限元格式.由于加上若干稳定化项和投影项,得到的混合变分形式是稳定的.我们证明了双线性形式满足连续性,K_h-强制性和Inf-Sup条件,因此,解是存在唯一的.此外,我们也给出了拟优的误差估计和相应的收敛阶.     

求解二次Lagrangian有限元方程的代数两水平方法

基于矩阵图集的粗化算法,构造一种新的插值算子,提出了瀑布型代数两重网格法;然后结合部分几何信息,提出了求解二次Lagrangian有限元方程的代数两水平方法.数值实验表明该算法稳健性强、计算量更少.     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.     

费马q-差分微分方程整函数解的增长性研究

本文研究了费马q-差分微分方程的整函数解的相关问题.利用经典和差分的Nevanlinna理论和函数方程理论的研究方法,获得了q-差分微分方程整函数解增长性的几个结果.     

一种基于新型插值单元的稳态传热边界元法

为了提高边界元法在求解稳态热问题时的计算精度,通过使用一种新型单元插值方法(称为扩展单元插值法),实现对稳态传热问题的求解。扩展单元是在传统不连续单元的边界配置虚拟节点,把原非连续单元变成高阶的连续单元,并将其作为新型的插值单元。利用虚拟节点和内部源节点构造出的插值函数,可以精确插值边界上的连续和不连续物理场,插值精度要比原始不连续单元高两阶。另外,边界积分方程只在传统的不连续单元的内部节点处建立,只包含内部源节点的自由度,而虚拟节点的自由度可通过与内部源节点之间的关系消除掉,因此最终系统方程的求解规模不会增加。这种新型的插值单元继承了传统连续和不连续单元的优点,克服了它们的缺点。数值结果表明,此种单元插值方法用于求解稳态传热问题时可获得较高的计算精度和收敛性。

     

非线性问题的插值摄动解法

本文用插值摄动法[1]求解几个非线性问题.算例表明,本文方法有很好的精度.     

The Iyengar inequality,revisited

Detailed analysis shows that the famous Iyengar inequality actually says that the Trapezoidal formula is a central algorithm for approximating integrals over an appropriate interval for the class of functions whose derivatives are bounded by a positive number K in L _∞-sense. The inherent nonlinearity from central algorithms reflects the importance of the Iyengar inequality and thus makes familiar linear methods malfunction when one tries to generalize it. It is shown that the generalization depends on a nonlinear system of equations satisfied by a set of free nodes of a perfect spline. Explicit constructions are obtained in the spirit of the Iyengar inequality for the class of functions whose rth (r≤4) derivatives are bounded by a positive number K in L _∞-sense because a closed solution to the nonlinear system is only available for r≤4. Connections with computational mathematics, especially with best interpolation and best quadrature, are discussed. Numerical experiments are also included. AMS subject classification (2000)  65D30, 41A17     

边界层型问题的插值摄动解法

本文在文[1]的基础上用插值摄动法研究了最高阶导数乘以小参数的二阶常微分方程的定解问题。算例表明,本文方法计算过程简单,其精度甚至比多重尺度法的一级近似结果的精度还稍高一些。     

A New Third S.N.Bernstein Interpolation Polynomial

§１．　ＩｎｔｒｏｄｕｃｔｉｏｎＷｅｋｎｏｗｎｔｈａｔｉｓｎｏｔｐｏｓｓｉｂｌｅｔｈａｔｔｈｅＬａｇｒａｎｇｅｉｎｔｅｒｐｏｌａｔｉｏｎｐｏｌｙｎｏｍｉａｌｓｃｏｎｖｅｒｇｅｕｎｉ－ｆｏｒｍｌｙｔｏｆ（ｘ）ｆｏｒａｌｌｃｏｎｔｉｎｕｏｕｓｆｕｎｃｔｉｏｎｓｆ（ｘ）ｏｎ［－１，１］．Ａｓａｒｅｓｕｌｔ，ｏｎｅｃａｎｕｓｅｖａｒｉｏｕｓｗａｙｓｔｏｃｈａｎｇｅｔｈｅＬａｇｒａｎｇｅｉｎｔｅｒｐｏｌａｔｉｏｎｐｏｌｙｎｏｍｉａｌｓｕｃｈｔｈａｔｔｈｅｙｕｎｉｆｏｒｍｌｙｃｏｎｖｅｒｇｅｔｏａｌｌｃｏｎ…