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.     

利用杨辉三角形对称性推导高阶运动微分方程

利用Mathematica数学软件计算函数r=r(q(t),t)各变量之间偏导和高阶导数的关系，发现具有杨辉三角形对称性.结合杨辉三角形的对称性规律和牛顿第二定律推导出了高阶运动微分方程，并讨论了理想约束系统下的高阶运动微分方程. 关键词： 杨辉三角形 牛顿第二定律 高阶运动微分方程 高阶力变率 高阶速度能量 理想约束     

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.     

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”.     

Geometry of Differential Equations: A Concise Introduction

A short introduction to geometrical theory of nonlinear differential equations is given to provide a unified overview to the collection 'Symmetries of differential equations and related topics'.     

亚纯函数及其微分多项式的唯一性

设f是非常数亚纯函数,g是f的线性微分多项式．a和b是f的两个不同的小函数．本文证明如果f和g几乎CM分担a和b,则f≡g;此外,若f是非常数整函数,且f和f^(k)(k≥1)IM分担a和b,b-a≠Pe^λz,测f≡g．     

Stable multi-sets

In this paper we introduce a generalization of stable sets: stable multi-sets. A stable multi-set is an assignment of integers to the vertices of a graph, such that specified bounds on vertices and edges are not exceeded. In case all vertex and edge bounds equal one, stable multi-sets are equivalent to stable sets.  For the stable multi-set problem, we derive reduction rules and study the associated polytope. We state necessary and sufficient conditions for the extreme points of the linear relaxation to be integer. These conditions generalize the conditions for the stable set polytope. Moreover, the classes of odd cycle and clique inequalities for stable sets are generalized to stable multi-sets and conditions for them to be facet defining are determined.  The study of stable multi-sets is initiated by optimization problems in the field of telecommunication networks. Stable multi-sets emerge as an important substructure in the design of optical networks. Received: February 14, 2001/Revised version: September 7, 2001     

方程w"-w+f(t,w)=O的Dirichlet边值问题的正解存在性与多解性

考察了下列常微分方程的Dirichlet边值问题的正解[w″(t)-w(t) f(t,w(t))=0,0≤t≤1 w(0)=w(1)=0建立了n正解的存在性，其中n是一个任意的自然数。     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

(e, 3e) test on e-e correlations in helium

The angular variations of the five-fold differential cross section obtained by using different wave functions of helium are compared with experimental data. It is found that in the coplanar geometry two kinematical arrangements, (i) equal energy sharing between the two ejected electrons with one of them ejected along the momentum transfer direction and the other along varying direction and (ii) the Bethe ridge condition with fixed sum of ejected electron energies and varying angle between them, are very sensitive to e-e correlations contained in the target wave function. This comparison has been used to show that open-shell class of wave functions better incorporate e-e correlations than the closed-shell class.