热动力学问题的数值计算

 对含热传导的流体动力学方程组，用有限元方法进行数值求解。采用傅里叶热传导计算热流、用热流连续条件计算单元间接触面的温度、用三角形传热法计算体单元表面的热流，考虑各向同性弹塑性流体材料模型、三项式物态方程和导热系数与状态的相关性，给出了傅里叶热传导、接触传热、热应力应变效应、以及混合物冲击压缩特性等算例。对混合物冲击温度的数值模拟表明，小颗粒混合物在冲击压缩过程中，颗粒间的温度有差别、稍有波动，并随时间趋向于一致，以至热平衡。     

Constraint incorporation in optimization

Numerical methods for solving constrained optimization problems need to incorporate the constraints in a manner that satisfies essentially competing interests; the incorporation needs to be simple enough that the solution method is tractable, yet complex enough to ensure the validity of the ultimate solution. We introduce a framework for constraint incorporation that identifies a minimal acceptable level of complexity and defines two basic types of constraint incorporation which (with combinations) cover nearly all popular numerical methods for constrained optimization, including trust region methods, penalty methods, barrier methods, penalty-multiplier methods, and sequential quadratic programming methods. The broad application of our framework relies on addition and chain rules for constraint incorporation which we develop here.     

沪深股市回报率、波动率和交易量关系的实证研究

本文首先对回报率与交易量之间的关系进行了研究,发现并不存在非对称的数量关系,但存在双向的葛兰杰因果关系;同时将交易量对波动率的解释能力进行了研究,发现在沪市交易量对波动率具有解释力,而在深市交易量对波动率没有解释力。     

On the Patterson–Sullivan Measure for Geometrically Finite Groups Acting on Complex or Quaternionic Hyperbolic Space

The goal of this paper is to provide a tool, the Global Measure Formula, that will facilitate the study of the limit set of discrete geometrically finite groups of isometries of the rank one symmetric spaces. We consider the shadow of a ball from a fixed reference point onto the boundary, and prove a formula that describes the measure of the shadow in terms of the center of the shadowed ball, generalizing a result from real hyperbolic geometry.     

Field-of-values analysis of preconditioned iterative methods for nonsymmetric elliptic problems

Summary. The convergence rate of Krylov subspace methods for the solution of nonsymmetric systems of linear equations, such as GMRES or FOM, is studied. Bounds on the convergence rate are presented which are based on the smallest real part of the field of values of the coefficient matrix and of its inverse. Estimates for these quantities are available during the iteration from the underlying Arnoldi process. It is shown how these bounds can be used to study the convergence properties, in particular, the dependence on the mesh-size and on the size of the skew-symmetric part, for preconditioners for finite element discretizations of nonsymmetric elliptic boundary value problems. This is illustrated for the hierarchical basis and multilevel preconditioners which constitute popular preconditioning strategies for such problems. Received May 3, 1996     

The influence of temperature gradients on the kinetic boundary layer problem for a condensing droplet

We extend an earlier method for solving kinetic boundary layer problems to the case of particles moving in aspatially inhomogeneous background. The method is developed for a gas mixture containing a supersaturated vapor and a light carrier gas from which a small droplet condenses. The release of heat of condensation causes a temperature difference between droplet and gas in the quasistationary state; the kinetic equation describing the vapor is the stationary Klein-Kramers equation for Brownian particles diffusing in a temperature gradient. By means of an expansion in Burnett functions, this equation is transformed into a set of coupled algebrodifferential equations. By numerical integration we construct fundamental solutions of this equation that are subsequently combined linearly to fulfill appropriate mesoscopic boundary conditions for particles leaving the droplet surface. In view of the intrinsic numerical instability of the system of equations, a novel procedure is developed to remove the admixture of fast growing solutions to the solutions of interest. The procedure is tested for a few model problems and then applied to a slightly simplified condensation problem with parameters corresponding to the condensation of mercury in a background of neon. The effects of thermal gradients and thermodiffusion on the growth rate of the droplet are small (of the order of 1%), but well outside of the margin of error of the method.     

Sparse-grid finite-volume multigrid for 3D-problems

We introduce a multigrid algorithm for the solution of a second order elliptic equation in three dimensions. For the approximation of the solution we use a partially ordered hierarchy of finite-volume discretisations. We show that there is a relation with semicoarsening and approximation by more-dimensional Haar wavelets. By taking a proper subset of all possible meshes in the hierarchy, a sparse grid finite-volume discretisation can be constructed.The multigrid algorithm consists of a simple damped point-Jacobi relaxation as the smoothing procedure, while the coarse grid correction is made by interpolation from several coarser grid levels.The combination of sparse grids and multigrid with semi-coarsening leads to a relatively small number of degrees of freedom,N, to obtain an accurate approximation, together with anO(N) method for the solution. The algorithm is symmetric with respect to the three coordinate directions and it is fit for combination with adaptive techniques.To analyse the convergence of the multigrid algorithm we develop the necessary Fourier analysis tools. All techniques, designed for 3D-problems, can also be applied for the 2D case, and — for simplicity — we apply the tools to study the convergence behaviour for the anisotropic Poisson equation for this 2D case.     

The many facets of linear programming

We examine the history of linear programming from computational, geometric, and complexity points of view, looking at simplex, ellipsoid, interior-point, and other methods. Received: June 22, 2000 / Accepted: April 4, 2001?Published online October 2, 2001     

On the rigidity of sphericalt-designs that are orbits of finite reflection groups

The concept of rigid sphericalt-designs was introduced by Bannai. He conjectured that there is a functionf(t, d) such that ifX is a sphericalt design in thed-dimensional Euclidean space so that |X|>f(t, d), theX is non-rigid. Furthermore, he asked to find examples of rigid but not tight sperical designs. In the present article we shall investigate the case whenX is an orbit of a finite reflection group and prove thatX is rigid iff tight for the groupsA _n ,B _n ,C _n ,D _n ,E ₆,E ₇,F ₄,I ₃.Research was partially supported by Hungarian National Research fund Grant No. 4267.