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.     

一类非线性抛物型变分不等方程的全局吸引子及弱近似惯性流形

对由一类非线性抛物型变分不等方程所描述的无穷维动力系统，给出了存在全局吸引子及弱近似惯性流形的充分条件．     

Studies on basicity of alkali metal catron exchanged β and X zeolites by pyrrole-IR spectroscopy

The strength of basic sites has been measured by pyrrole-IR on alkali metal cation exchanged β and X zeolites, as well as NaOH loaded Naβ. The influence of cation type and the structure of zeolites on their basicity has been studied. The acidic and basic properties of the samples were investigated by NH3-TPD and isopropanol reaction. It was shown that the strength of basic sites on samples could be characterized by the shift of vNH band in the pyrrole-IR spectra. The framework oxygen charges were calculated from the Sanderson electronegativity. The changes in basic properties with various alkali metal cation are consistent with the changes of local oxygen charges of the zeolite framework.     

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.     

A likelihood-based method for haplotype association studies of case-control data with genotyping uncertainty

This paper discusses the associations between traits and haplotypes based on Fl (fluorescent intensity) data sets. We consider a clustering algorithm based on mixtures of t distributions to obtain all possible genotypes of each individual (i.e. "GenoSpec-trum"). We then propose a likelihood-based approach that incorporates the genotyping uncertainty to assessing the associations between traits and haplotypes through a haplotype-based logistic regression model. Simulation studies show that our likelihood-based method can reduce the impact induced by genotyping errors.     

VRPTW的扰动恢复及其TABU SEARCH算法

本文对带时间窗的车辆路线安排扰动恢复问题进行了讨论,分析了各种可能的扰动:增加减少客户,时间窗、客户需求及路线可行性的扰动,构造了扰动模型.利用禁忌搜索算法对问题进行求解,同时通过对模型参数重新设置,得到了多个满足要求的不同的解,这样使解更具有实际可行性和有效性.     

Simulation of a High Harmonic Gyrotron with Axis-Encircling Electron Beam and Permanent Magnet

A computer code for simulation of beam-field interaction in a resonant cavity of a gyrotron has been developed. It is based on a self-consistent, time-independent, single-mode physical model. The code has been applied to the analysis of the electrodynamical system of a novel high-harmonic gyrotron with axis-encircling beam and a permanent magnet. In this paper both the physical model and numerical techniques used are outlined. Some results of the numerical experiments are presented and discussed.     

B_2分子X~3Σ_g~-和A~3Σ_u~-态的势能函数

使用SAC/SAC-CI方法,利用D95(d),6-311g**以及cc-PVTZ等基组,对B2分子的基态(X3Σg-)和第一激发态(A3Σu-)的平衡结构和谐振频率进行了优化计算.通过对3个基组的计算结果的比较,得出了D95(d)基组为3个基组中的最优基组的结论;使用D95(d)基组,利用SAC的GSUM(Group Sum of Operators)方法对基态(X3Σg-),SAC-CI的GSUM方法对激发态(A3Σu-)进行单点能扫描计算,用正规方程组拟合Murrell-Sorbie函数,得到了相应电子态的完整势能函数;从得到的势能函数计算了与基态(X3Σg-)和第一激发态(A3Σu-)相对应的光谱常数(Be,αe,ωe和ωeχe),结果与实验数据吻合.