曲率流的参数化有限元逼近

许多物理现象可以在数学上描述为受曲率驱动的自由界面运动,例如薄膜和泡沫的演变、晶体生长,等等.这些薄膜和界面的运动常依赖于其表面曲率,从而可以用相应的曲率流来描述,其相关自由界面问题的数值计算和误差分析一直是计算数学领域中的难点.参数化有限元法是曲率流的一类有效计算方法,已经能够成功模拟一些曲面在几类基本的曲率流下的演化过程.本文重点讨论曲率流的参数化有限元逼近,它的产生、发展和当前的一些挑战.     

6元素集合上T0拓扑总数的计算

文中证明了有限预序集与有限偏序集的一些性质,并基于有限集上的拓扑和其上预序的一一对应关系,利用这些性质通过对极小元和极大元个数进行分类讨论,以一种有别于计算机算法而又容易理解的计算方法得出6元素集合上的T₀拓扑总数为130023.     

Schrödinger方程的连续时空有限元方法

本文讨论Schr?dinger方程的连续时空有限元方法,通过引入相应的时空投影算子,利用实部虚部分离技巧,得到了变量u在时间节点处的L²范数,以及u和u_t的全局L²(H¹)和L²(L²)范数意义下的最优误差估计结果.该文的结论对进一步探索和设计Schr?dinger方程的数值算法是有益的.     

Randomized algorithms for the separation of point sets and for solving quadratic programs

A randomized algorithm for finding a hyperplane separating two finite point sets in the Euclidean space ^d and a randomized algorithm for solving linearly constrained general convex quadratic problems are proposed. The expected running time of the separating algorithm isO(dd! (m + n)), wherem andn are cardinalities of sets to be separated. The expected running time of the algorithm for solving quadratic problems isO(dd! s) wheres is the number of inequality constraints. These algorithms are based on the ideas of Seidel's linear programming algorithm [6]. They are closely related to algorithms of [8], [2], and [9] and belong to an abstract class of algorithms investigated in [1]. The algorithm for solving quadratic problems has some features of the one proposed in [7].This research was done when the author was supported by the Alexander von Humboldt Foundation, Germany.On leave from the Institute of Mathematics and Mechanics, Ural Department of the Russian Academy of Sciences, 620219 Ekaterinburg, S. Kovalevskaya str. 16, Russia.     

On the classification of geometric codes by polynomial functions

It is shown how to represent algebraically all functions that have a zero sum on all -dimensional subspaces ofPG(n,q) or ofAG(n,q). In this way one can calculate the dimensions of related codes, or one can represent interesting sets of points by functions.     

Fastiterative methods for least squares estimations

Least squares estimations have been used extensively in many applications, e.g. system identification and signal prediction. When the stochastic process is stationary, the least squares estimators can be found by solving a Toeplitz or near-Toeplitz matrix system depending on the knowledge of the data statistics. In this paper, we employ the preconditioned conjugate gradient method with circulant preconditioners to solve such systems. Our proposed circulant preconditioners are derived from the spectral property of the given stationary process. In the case where the spectral density functions() of the process is known, we prove that ifs() is a positive continuous function, then the spectrum of the preconditioned system will be clustered around 1 and the method converges superlinearly. However, if the statistics of the process is unknown, then we prove that with probability 1, the spectrum of the preconditioned system is still clustered around 1 provided that large data samples are taken. For finite impulse response (FIR) system identification problems, our numerical results show that annth order least squares estimator can usually be obtained inO(n logn) operations whenO(n) data samples are used. Finally, we remark that our algorithm can be modified to suit the applications of recursive least squares computations with the proper use of sliding window method arising in signal processing applications.Research supported in part by HKRGC grant no. 221600070, ONR contract no. N00014-90-J-1695 and DOE grant no. DE-FG03-87ER25037.     

An Arithmetic Reduction of Finite Rank 3 Geometries with Linear Spaces as Plane Residues and with Dual Linear Spaces as Point Residues

Let be a rank three incidence geometry of points, lines and planes whose planes are linear spaces and whose point residues are dual linear spaces (notice that we do not require anything on the line residues). We assume that the residual linear spaces of belong to a natural class of finite linear spaces, namely those linear spaces whose full automorphism group acts flag-transitively and whose orders are polynomial functions of some prime number. This class consists of six families of linear spaces. In the amalgamation of two such linear spaces imposes an equality on their orders leading, in particular, to a series of diophantine equations, the solutions of which provide a reduction theorem on the possible amalgams of linear spaces that can occur in .We prove that one of the following holds (up to a permutation of the words point and plane):A) the planes of and the dual of the point residues belong to the same family and have the same orders,B) the diagram of is in one of six families,C) the diagram of belongs to a list of seven sporadic cases.Finally, we consider the particular case where the line residues of are generalized digons.     

Unbiased Bayes estimates and improper priors

Given two random variables (X, Y) the condition of unbiasedness states that:E(X |Y=y)=y andE(Y |X=x)=x both almost surely (a.s.). If the prior onY is proper and has finite expectation or non-negative support, unbiasedness impliesX=Y a.s. This paper examines the implications of unbiasedness when the prior onY is improper. Since the improper case can be meaningfully analysed in a finitely additive framework, we revisit the whole issue of unbiasedness from this perspective. First we argue that a notion weaker than equality a.s., named coincidence, is more appropriate in a finitely additive setting. Next we discuss the meaning of unbiasedness from a Bayesian and fiducial perspective. We then show that unbiasedness and finite expectation ofY imply coincidence betweenX andY, while a weaker conclusion follows if the improper prior onY is only assumed to have positive support. We illustrate our approach throughout the paper by revisiting some examples discussed in the recent literature.This work was partially supported by C.N.R. grant N.80.02970.10 (G.C.) and by C.N.R. grant altri interventi (P.V.). A preliminary draft was written while the Authors were visiting the Department of Statistics at Carnegie Mellon University.The paper is the result of close cooperation between the two authors. However subsections 3.1 and 3.3 are mainly due to G.C. while subsection 3.2 and section 4 are mainly due to P.V.     

Surface tension and universality in the three-dimensional Ising model

The amplitude ₀ of the interfacial free energy per unit area (or surface tension) of the body-centered-cubic Ising model is found using a direct monte carlo simulation technique. The combination ²/k_BT_c, where is the correlation length, is shown to agree within the precision of the simulations with a previously reported estimate for the simple cubic lattice. Evidence is also presented for the universality of the finite-size scaling amplitude for the surface tension.     

Analysis of thermal behaviour of high power semiconductor laser arrays by means of the finite element method (FEM)

New results of steady-state two-dimensional finite-element computations of temperature distributions of high power semiconductor laser arrays are presented. The influence of different thermal loads on the 2D temperature distribution in AlGaAs/GaAs gain-guided laser arrays is investigated. TheFEM model is tested by comparing it with analytical solutions. For numerical convenience, the latter is rewritten in a novel form, which is free of overflow problems. The maximum temperatures calculated by both methods agree within 1%. Several factors determining the thermal resistance of the device are quantitatively examined: the ratio of light emitting to non-emitting areas along the active zone, the amount of Joule losses, the current spreading, the solder thickness, and voids in the solder. This yields design rules for optimum thermal performance.