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.     

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.     

Implementation of the finite field perturbation method in the CRYSTAL program for calculating the dielectric constant of periodic systems

The finite field approach has been implemented in the periodic ab initio CRYSTAL program and been used for calculating the dielectric constants of crystalline LiF and MgO (FCC structure) and BeO (wurtzite structure). To maintain the periodicity along the applied field direction, a "sawtooth" potential is used in conjunction with a supercell scheme. Supercells four to five times longer than the primitive cell in the direction of the applied field provide well-converged results. The influence of the computational parameters is discussed. An alternative scheme has also been implemented, for inner check, that consists of applying a static electric field to a slab of increasing thickness in the direction orthogonal to the surface; the dielectric response at the center of the slab is shown to converge rapidly to the bulk value evaluated with the sawtooth field. The method is accurate and permits the determination of nonlinear corrections to the dielectric constant. When used in conjunction with the local density approximation (LDA) scheme, it provides for the dielectric constant of the three above-mentioned compounds values close to those recently obtained with a time-dependent density functional theory approach.     

On the modelling of solid state reactions.Synthesis of YAG

There is a model of yttrium aluminium garnet (YAG) synthesis presented in this article. The developed model is based on nonlinear reaction–diffusion partial differential equations. The solution was carried out numerically using finite difference techniques. We got dependability curves for diffusion and reaction rates and offered possible method to localize values of diffusion and reaction rate constants precisely enough.AMS subject classification: 35K57, 65M06     

The hole-pressure problem

There is a need to unify present hypotheses of the nature and role of the hole-pressure,p _e , and thus provide consolidation on which to base future research and understanding. This paper is intended to meet this need. Attention is directed towards the calculation ofp _e from the velocity and stress fields for viscoelastic fluids flowingacross rectangular holes. The constitutive models used are the Newtonian, Second-order and Maxwell models, for values of Reynolds number up to 10 and Weissenberg number up to 0.1.The numerical complications involved are studied through an investigation of the constituent parts ofp _e . Verification of present theory is then sought, from which justification may be derived for the estimation of elasticity fromp _e measurements. Attention is directed towards the predictions of Higashitani and Pritchard and the extension to the Tanner and Pipkin theory for Second-order fluids. The effects of variation of geometric dimensions and flow type uponp _e are also discussed.