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.     

关于混合型多项式Hammerstein积分方程正解迭代求法

本文给出了混合型多项式Ｈａｍｍｅｒｓｔｅｉｎ积分方程正解的迭代求法，并将所得结果应用到二阶非线性常微分方程的边值问题     

压力淬火过程中锗的相结构研究

 应用高压原位差热方法，直接测量了压力下锗的固化参数─—固化温度与过冷度。高压差热信号表明，当压力大于3 GPa时，锗在凝固过程中可能发生结构相变。X射线结构分析表明，在最终的样品中除GeⅠ相外，还形成GeⅢ相和GeⅣ相。     

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.     

SINGULAR PERTURBATION OF BOUNDARY VALUE PROBLEM FOR NONLINEAR HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS INVOLVING TWO SMALL PARAMETERS

The singularly perturbed boundary value problem for nonlinear higher order ordinary differential equation involving two small parameters has been considered. Under appropriate assumptions, for the three cases:ε/μ2→0(μ→0),μ2/ε→0 (ε→0) andε=μ2, the uniformly valid asymptotic solution is obtained by using the expansion method of two small parameters and the theory of differential inequality.     

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

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

Attractors for a Three-Dimensional Thermo-Mechanical Model of Shape Memory Alloys

Abstract In this note, we consider a Frémond model of shape memory alloys. Let us imagine a piece of a shape memory alloy which is fixed on one part of its boundary, and assume that forcing terms, e.g., heat sources and external stress on the remaining part of its boundary, converge to some time-independent functions, in appropriate senses, as time goes to infinity. Under the above assumption, we shall discuss the asymptotic stability for the dynamical system from the viewpoint of the global attractor. More precisely, we generalize the paper [12] dealing with the one-dimensional case. First, we show the existence of the global attractor for the limiting autonomous dynamical system; then we characterize the asymptotic stability for the non-autonomous case by the limiting global attractor. * Project supported by the MIUR-COFIN 2004 research program on “Mathematical Modelling and Analysis of Free Boundary Problems”.     

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.