A posteriori error estimation with finite element methods of lines for one-dimensional parabolic systems

Summary Consider the solution of one-dimensional linear initial-boundary value problems by a finite element method of lines using a piecewiseP ^th -degree polynomial basis. A posteriori estimates of the discretization error are obtained as the solutions of either local parabolic or local elliptic finite element problems using piecewise polynomial corrections of degreep+1 that vanish at element ends. Error estimates computed in this manner are shown to converge in energy under mesh refinement to the exact finite element discretization error. Computational results indicate that the error estimates are robust over a wide range of mesh spacings and polynomial degrees and are, furthermore, applicable in situations that are not supported by the analysis.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR 90-0194; by the U.S. Army Research Office under Contract Number DAAL03-91-G-0215; and by the National Science Foundation under Institutional Infrastructure Grant Number CDA-8805910     

On the Moments of Overflow and Freed Carried Traffic for the GI/M/C/0 System

In circuit-switched networks call streams are characterized by their mean and peakedness (two-moment method). The GI/M/C/0 system is used to model a single link, where the GI-stream is determined by fitting moments appropriately. For the moments of the overflow traffic of a GI/M/C/0 system there are efficient numerical algorithms available. However, for the moments of the freed carried traffic, defined as the moments of a virtual link of infinite capacity to which the process of calls accepted by the link (carried arrival process) is virtually directed and where the virtual calls get fresh exponential i.i.d. holding times, only complex numerical algorithms are available. This is the reason why the concept of the freed carried traffic is not used. The main result of this paper is a numerically stable and efficient algorithm for computing the moments of freed carried traffic, in particular an explicit formula for its peakedness. This result offers a unified handling of both overflow and carried traffics in networks. Furthermore, some refined characteristics for the overflow and freed carried streams are derived.     

Point approximation of a space-homogeneous transport equation

Summary We present an approximation method of a space-homogeneous transport equation which we prove is convergent. The method is very promising for numerical computation. Comparison of a numerical computation with an exact solution is given for the Master equation.     

C72,La2@C72几何结构和电子性质的计算研究

采用密度泛涵理论(density functional theory,DFT)中的广义梯度近似(generally gradient approximation,GGA)对富勒烯C72和内掺金属La富勒烯La2@C72三种同分异构体的几何结构和电子结构进行研究.发现在C72的三种同分异构体中,满足独立五边形规则(isolated-pentagon-rule,IPR)的C72(D6d)结构最为稳定;在La2@C72三种同分异构体中,有着两对两两相邻五元环(twofused-pentagon)的La2@C72(#10611)结构最为稳定,而满足IPR的La2@C72(D6d)的结构变成了最不稳定结构.从能级图和态密度图分析得知,笼子稳定性的变化与La原子的原子轨道与C72原子轨道之间的杂化有关.Mulliken电荷分析得知,La2@C72(#10611)的两个La原子共转移了约3个电子给C72,并且,它们几乎分布在整个C笼上,形成的电子结构为La23+@C3-72.净自旋分析得知,La2@C72(#10611)中La原子磁性完全淬灭.     

Optimal tristance anticodes in certain graphs

For z₁,z₂,z₃∈Zⁿ, the tristance d₃(z₁,z₂,z₃) is a generalization of the L₁-distance on Zⁿ to a quantity that reflects the relative dispersion of three points rather than two. A tristance anticodeA_d of diameter d is a subset of Zⁿ with the property that d₃(z₁,z₂,z₃)?d for all z₁,z₂,z₃∈A_d. An anticode is optimal if it has the largest possible cardinality for its diameter d. We determine the cardinality and completely classify the optimal tristance anticodes in Z² for all diameters d?1. We then generalize this result to two related distance models: a different distance structure on Z² where d(z₁,z₂)=1 if z₁,z₂ are adjacent either horizontally, vertically, or diagonally, and the distance structure obtained when Z² is replaced by the hexagonal lattice A₂. We also investigate optimal tristance anticodes in Z³ and optimal quadristance anticodes in Z², and provide bounds on their cardinality. We conclude with a brief discussion of the applications of our results to multi-dimensional interleaving schemes and to connectivity loci in the game of Go.     

Angle theorems for the Lagrangian mean curvature flow

We prove that symplectic maps between Riemann surfaces L, M of constant, nonpositive and equal curvature converge to minimal symplectic maps, if the Lagrangian angle for the corresponding Lagrangian submanifold in the cross product space satisfies . If one considers a 4-dimensional K?hler-Einstein manifold of nonpositive scalar curvature that admits two complex structures J, K which commute and assumes that is a compact oriented Lagrangian submanifold w.r.t. J such that the K?hler form w.r.t.K restricted to L is positive and , then L converges under the mean curvature flow to a minimal Lagrangian submanifold which is calibrated w.r.t. . Received: 11 April 2001 / Published online: 29 April 2002     

Absolutely indecomposable symmetric matrices

Let A be a symmetric matrix of size n×n with entries in some (commutative) field K. We study the possibility of decomposing A into two blocks by conjugation by an orthogonal matrix T∈Mat_n(K). We say that A is absolutely indecomposable if it is indecomposable over every extension of the base field. If K is formally real then every symmetric matrix A diagonalizes orthogonally over the real closure of K. Assume that K is a not formally real and of level s. We prove that in Mat_n(K) there exist symmetric, absolutely indecomposable matrices iff n is congruent to 0, 1 or −1 modulo 2s.     

Global stability of the Armstrong-Frederick model with periodic biaxial inputs

The paper is concerned with the study of plasticity models described by differential equations with stop and play operators. We suggest sufficient conditions for the global stability of a unique periodic solution for the scalar models and for the vector models with biaxial inputs of a particular form, namely the sum of a uniaxial function and a constant term. For another class of simple biaxial inputs, we present an example of the existence of unstable periodic solutions. The paper was written during the research stay of D. Rachinskii at the Technical University Munich supported by the research fellowship from the Alexander von Humboldt Foundation. His work was partially supported by the Russian Science Support Foundation, Russian Foundation for Basic Research (Grant No. 01-01-00146, 03-01-00258), and the Grants of the President of Russia (Grant No. MD-87.2003.01, NS-1532.2003.1). The support is gratefully acknowledged.