Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

Numerical integration of ordinary differential equations on manifolds

Summary This paper is concerned with the problem of developing numerical integration algorithms for differential equations that, when viewed as equations in some Euclidean space, naturally evolve on some embedded submanifold. It is desired to construct algorithms whose iterates also evolve on the same manifold. These algorithms can therefore be viewed as integrating ordinary differential equations on manifolds. The basic method “decouples” the computation of flows on the submanifold from the numerical integration process. It is shown that two classes of single-step and multistep algorithms can be posed and analyzed theoretically, using the concept of “freezing” the coefficients of differential operators obtained from the defining vector field. Explicit third-order algorithms are derived, with additional equations augmenting those of their classical counterparts, obtained from “obstructions” defined by nonvanishing Lie brackets.     

Local and Parallel Finite Element Algorithms for Eigenvalue Problems

Abstract Some new local and parallel finite element algorithms are proposed and analyzed in this paper foreigenvalue problems.With these algorithms, the solution of an eigenvalue problem on a fine grid is reduced tothe solution of an eigenvalue problem on a relatively coarse grid together with solutions of some linear algebraicsystems on fine grid by using some local and parallel procedure.A theoretical tool for analyzing these algorithmsis some local error estimate that is also obtained in this paper for finite element approximations of eigenvectorson general shape-regular grids.     

Estimation of multivariate binary density using orthogonal functions

In this paper, the authors studied certain properties of the estimate of Liang and Krishnaiah (1985, J. Multivariate Anal. 16, 162–172) for multivariate binary density. An alternative shrinkage estimate is also obtained. The above results are generalized to general orthonormal systems.     

四角晶相HfO2(001)表面原子和电子结构研究

采用基于第一性原理的密度泛函理论研究了四角晶相二氧化铪(t-HfO₂)体相及 其(001)表面的原子几何与电子结构.理论计算结果表明，t-HfO₂(001)表面不会 产生重构现象.与体相电子结构相比, t-HfO₂(001)表面态密度明显高于体相态 密度.其次，表面原子的态密度更靠近费米能级(E_F)，价带往低能量处移动，并 有表面态产生.计算结果表明了t-HfO₂表面禁带宽度明显低于体相的禁带宽度. t-HfO₂(001)的表面态产生以及表面禁带宽度减小是由于Hf原子与O原子的配位 数减少，表面原子周围的环境发生变化而引起的. 关键词： 密度泛函理论 2(001)')" href="#">t-HfO₂(001) 表面电子结构     

Structural and thermodynamic properties of AlB2 compound

We employ a first-principles plane wave method with the relativistic analytic pseudopotential of Hartwigsen, Goedecker and Hutter (HGH) scheme in the frame of DFT to calculate the equilibrium lattice parameters and the thermodynamic properties of AlB₂ compound with hcp structure. The obtained lattice parameters are in good agreement with the available experimental data and those calculated by others. Through the quasi-harmonic Debye model, obtained successfully are the dependences of the normalized lattice parameters a/a₀ and c/c₀ on pressure P, the normalized primitive cell volume V/V₀ on pressure P, the variation of the thermal expansion α with pressure P and temperature T, as well as the Debye temperature \Theta_D and the heat capacity C_V on pressure P and temperature T.     

水斗非定常自由水膜流三维贴体数值模拟

本研究采用水斗三维非正交贴体坐标系进行了非定常自由水膜流动的数值解析。对不规则水斗内表面采用三维非正交贴体坐标系下离散点进行拟合,推导了曲面离散点的法向矢量和曲面微元面高斯曲率、平均曲率等几何特征量的计算公式,进而导出流体粒子在运动方向上曲率计算式。在水斗三维贴体坐标系中,还推导了流体粒子在水斗曲面上的运动控制方程。最后对某水轮机水斗内表面非定常自由水膜流进行了数值模拟,得到其非定常水膜流态分布。     

超强脉冲激光在低密度等离子体中的相对论自导引效应

分析了相对论条件下激光超短脉冲在等离子体中的传输特性 ,在傍轴近似和慢变振幅近似条件下 ,推导了折射率、电子密度、静电场以及电子空腔尺度的表达式。当激光功率超过产生自导引阈值功率时 ,激光束斑沿着传输光轴方向振荡。在有质动力产生的压力非常强时 ,聚焦光束中央部分的电子被全部排开形成电子空腔。给出了电子空腔的尺寸以及在出现电子空腔时的处理方法。在超过形成电子空腔的阈值功率 (Pc≈ 2 .5TW )时 ,空腔的尺度几乎与激光功率无关 ,这意味着电子空腔阻止了激光脉冲的进一步聚焦。