Asymptotic Methods Applied to Contact Problems

Asymptotic methods for contact problems are expounded. Some typical integral equations are considered     

An implicit characteristic-flux-averaging method for the euler equations for real gases

A formulation of an implicit characteristic-flux-averaging method for the unsteady Euler equations with real gas effects is presented. Incorporation of a real gas into a general equation of state is achieved by considering the pressure as a function of density and specific internal energy. The Ricmann solver as well as the flux-split algorithm are modified by introducing the pressure derivatives with respect to density and internal energy. Expressions for calculating the values of the flow variables for a real gas at the cell faces are derived. The Jacobian matrices and the eigenvectors are defined for a general equation of state. The solution of the system of equations is obtained by using a mesh-sequencing method for acceleration of the convergence. Finally, a test case for a simple form of equation of state displays the differences from the corresponding solution for an ideal gas.     

On‐line Determination of Particle Size Distributions in Flowing Dispersions

Emulsions are of great importance to industry. They are involved in many engineering operations, including chemical reactions, extraction, emulsification and suspension polymerization, etc. However, an important problem for these processes is how to control the size distribution of the dispersed phase. Indeed, off‐line analysis of the emulsion may generate uncertainties due to sampling and dilution of the product, which are likely to change the dispersion state and physico‐chemical properties. In this work, an on‐line optical method is proposed to characterize dispersed media in real flowing conditions. This method is based on the time‐analysis of back‐scattered light fluctuations. The present paper deals with the development of this method and its application to dispersions of alumina in water. The results obtained with the on‐line optical method are compared with those acquired by classical laser light scattering and microscopy.     

Ultrathin Fe film on Cu(0 0 1): Exchange splitting of image states from first principles

The theoretical investigation of the image states in front of an ultrathin iron film grown on copper has been performed by means of the embedding method and a recently developed procedure for the inclusion of the image potential tail in a first principle calculation. From the electronic response to an applied electric field, the image plane position has been evaluated. This also allows one to obtain useful information about the spin dependent screening properties of the system. Exchange splitting, effective mass, and lifetime of such surface states result in good agreement with recently performed two-photon photoemission experiments [see A.B. Schmidt, M. Pickel, M. Wiemhöfer, M. Donath, M. Weinelt, Phys. Rev. Lett. 95 (2005) 107402].     

Pattern Search Method for Discrete <Emphasis Type="Italic">L</Emphasis><Subscript>1</Subscript>–Approximation

We propose a pattern search method to solve a classical nonsmooth optimization problem. In a deep analogy with pattern search methods for linear constrained optimization, the set of search directions at each iteration is defined in such a way that it conforms to the local geometry of the set of points of nondifferentiability near the current iterate. This is crucial to ensure convergence. The approach presented here can be extended to wider classes of nonsmooth optimization problems. Numerical experiments seem to be encouraging. This work was supported by M.U.R.S.T., Rome, Italy.     

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.     

粒子鉴别中dE/dx和TOF的权重

在北京谱仪数据分析中,粒子鉴别使用了dE/dx和TOF的信息.在利用dE/dx和TOF联合鉴别粒子时,通常的做法是把两者的χ²值等权重相加(χ²方法).在不同的动量区间,dE/dx和TOF对粒子的鉴别能力不同.本文给出两者联合鉴别时不同动量区间的权重,并构造一个新的线性变量来进行粒子鉴别     

Charge density fluctuation of low frequency in a dusty plasma

The charge density fluctuation of low frequency in a dusty plasma, which is derived from the longitudinal dielectric permittivity of the dusty plasma, has been studied by kinetic theory. The results show that theP value, which describes the relative charge density on the dust in the plasma, and the charging frequency of a dust particle Ω _c , which describes the ratio of charge changing of the dust particles, determine the character of the charge density fluctuation of low frequency. For a dusty plasma ofP≪1, when the charging frequency Ω _c , is much smaller than the dusty plasma frequency ω_d, there is a strong charge density fluctuation which is of character of dust acoustic eigenwave. For a dusty plasma ofP≫1, when the frequency Ω _c , is much larger than ω _d there are weaker fluctuations with a wide spectrum. The results have been applied to the ionosphere and the range of radius and density of dust particles is found, where a strong charge density fluctuation of low frequency should exist.     

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.