A Monte Carlo method for high dimensional integration

Summary A new method for the numerical integration of very high dimensional functions is introduced and implemented based on the Metropolis' Monte Carlo algorithm. The logarithm of the high dimensional integral is reduced to a 1-dimensional integration of a certain statistical function with respect to a scale parameter over the range of the unit interval. The improvement in accuracy is found to be substantial comparing to the conventional crude Monte Carlo integration. Several numerical demonstrations are made, and variability of the estimates are shown.     

Well-posed absorbing layer for hyperbolic problems

Summary. The perfectly matched layer (PML) is an efficient tool to simulate propagation phenomena in free space on unbounded domain. In this paper we consider a new type of absorbing layer for Maxwell's equations and the linearized Euler equations which is also valid for several classes of first order hyperbolic systems. The definition of this layer appears as a slight modification of the PML technique. We show that the associated Cauchy problem is well-posed in suitable spaces. This theory is finally illustrated by some numerical results. It must be underlined that the discretization of this layer leads to a new discretization of the classical PML formulation. Received May 5, 2000 / Published online November 15, 2001     

Holographic interferometry as electrochemical emission spectroscopy of a carbon steel in 5–20 ppm TROS C-70 inhibited seawater

In the present investigation, holographic interferometry was utilized for the first time to determine the rate change of the number of the fringe evolutions during the corrosion test of a carbon steel in blank seawater and seawater with different concentrations of a corrosion inhibitor. In other words, the anodic dissolution behaviors (corrosion) of the carbon steel were determined simultaneously by holographic interferometry, as an electromagnetic method, and by the electrochemical impedance (EI) spectroscopy, as an electronic method. So, the abrupt rate change of the number of the fringe evolutions during corrosion tests, EI spectroscopy, of the carbon steel is called electrochemical emission spectroscopy. The corrosion process of the steel samples was carried out in blank seawater and seawater with different concentrations, 5–20 ppm, of TROS C-70 corrosion inhibitor using the EI spectroscopy method, at room temperature. The electrochemical-emission spectra of the carbon steel in different solutions represent a detailed picture of the rate change of the anodic dissolution of the steel throughout the corrosion processes. Furthermore, the optical interferometry data of the carbon steel were compared to the data, which was obtained from the EI. spectroscopy. Consequently, holographic interferometry is found very useful for monitoring the anodic dissolution behaviors of metals, in which the number of the fringe evolutions of the steel samples can be determined in situ.     

Structural and lithium intercalation studies of Mn(0.5−x)CaxTi2(PO4)3 phases (0≤x≤0.50)

Materials from the Mn_(0.5−x)Ca_xTi₂(PO₄)₃ (0≤x≤0.50) solid solution were obtained by solid-state reaction in air at 1000 °C. Selected compositions were investigated by powder X-ray diffraction analysis, ³¹P nuclear magnetic resonance (NMR) spectroscopy and electrochemical lithium intercalation. The structure of all samples determined by Rietveld analysis is of the Nasicon type with the R space group. Mn²⁺/Ca²⁺ ions occupy only the M1 sites in the Ti₂(PO₄)₃ framework. The divalent cations are ordered in one of two M1 sites, except for the Mn_0.50Ti₂(PO₄)₃ phase, where a small departure from the ideal order is observed by XRD and ³¹P MAS NMR. The electrochemical behaviour of Mn_0.50Ti₂(PO₄)₃ and Mn_(0.5−x)Ca_xTi₂(PO₄)₃ phases was characterised in Li cells. Two Li ions can be inserted without altering the Ti₂(PO₄)₃ framework. In the 0≤y≤2 range, the OCV curves of Li//Li_yMn_0.50Ti₂(PO₄)₃ cells show two main potential plateaus at 2.90 and 2.50–2.30 V. Comparison between the OCV curves of Li//Li_(1+y)Ti₂(PO₄)₃ and Li//Li_yMn_0.50Ti₂(PO₄)₃ shows that the intercalation occurs first in the unoccupied M1 site of Mn_0.50Ti₂(PO₄)₃ at 2.90 V and then, for compositions y>0.50, at the M2 site (2.50–2.30 V voltage range). The effect of calcium substitution in Mn_0.50Ti₂(PO₄)₃ on the lithium intercalation is also discussed from a structural and kinetic viewpoint. In all systems, the lithium intercalation is associated with a redistribution of the divalent cation over all M1 sites. In the case of Mn_0.50Ti₂(PO₄)₃, the stability of Mn²⁺ either in an octahedral or tetrahedral environment facilitates cationic migration.     

Limits of parallelism in explicit ODE methods

Numerical methods for ordinary initial value problems that do not depend on special properties of the system are usually found in the class of linear multistage multivalue methods, first formulated by J.C. Butcher. Among these the explicit methods are easiest to implement. For these reasons there has been considerable research activity devoted to generating methods of this class which utilize independent function evaluations that can be performed in parallel. Each such group of concurrent function evaluations can be regarded as a stage of the method. However, it turns out that parallelism affords only limited opportunity for reducing the computing time with such methods. This is most evident for the simple linear homogeneous constant-coefficient test problem, whose solution is essentially a matter of approximating the exponential by an algebraic function. For a given number of stages and a given number of saved values, parallelism offers a somewhat enlarged set of algebraic functions from which to choose. However, there is absolutely no benefit in having the degree of parallelism (number of processors) exceed the number of saved values of the method. Thus, in particular, parallel one-step methods offer no speedup over serial one-step methods for the standard linear test problem. Although the implication of this result for general nonlinear problems is unclear, there are indications that dramatic speedups are not possible in general. Also given are some results relevant to the construction of methods.Work supported in part by National Science Foundation grants DMS 89 11410 and DMS 90 15533 and by US Department of Energy grant DOE DEFG02-87ER25026. Work of the second author was completed while at the University of Illinois.     

中心—焦点区分问题的I l'yashenko算法

对二维平面系统的二维系统的中心焦点区分问题 ,I l'yashenko曾建议一个算法 ,本文给出此法的详细证明 .据此 ,我们讨论了区分问题在 Arnold意义下的代数可解性与不可解性     

Multiderivations of Coxeter arrangements

Let V be an ℓ-dimensional Euclidean space. Let G⊂O(V) be a finite irreducible orthogonal reflection group. Let ? be the corresponding Coxeter arrangement. Let S be the algebra of polynomial functions on V. For H∈? choose α _H ∈V ^* such that H=ker(α _H ). For each nonnegative integer m, define the derivation module D ^(m) (?)={θ∈Der _S |θ(α _H )∈Sα ^m _H }. The module is known to be a free S-module of rank ℓ by K. Saito (1975) for m=1 and L. Solomon-H. Terao (1998) for m=2. The main result of this paper is that this is the case for all m. Moreover we explicitly construct a basis for D ^(m) (?). Their degrees are all equal to mh/2 (when m is even) or are equal to ((m−1)h/2)+m _i (1≤i≤ℓ) (when m is odd). Here m ₁≤···≤m _ℓ are the exponents of G and h=m _ℓ+1 is the Coxeter number. The construction heavily uses the primitive derivation D which plays a central role in the theory of flat generators by K. Saito (or equivalently the Frobenius manifold structure for the orbit space of G). Some new results concerning the primitive derivation D are obtained in the course of proof of the main result. Oblatum 27-XI-2001 & 4-XII-2001?Published online: 18 February 2002     

On the existence of nontrivial extremal metrics on complete noncompact surfaces

In this paper, we consider the 2-dimensional local Calabi flow on a complete noncompact surface . Then, based on the Harnack-type estimate, we show the long-time existence and asymptotic convergence of a subsequence of solutions of such a flow on with and bounded from above by a negative constant on a ball. For its applications, this will lead to the existence of extremal metrics on a complete noncompact surface of finite topological type. In particular, there exists an extremal metric of nonconstant Gaussian curvature on or Received: 21 June 2001 / 18 January 2002 / Published online: 27 June 2002 Research supported in part by NSC and NCTS.     

Multispectral imaging of Mars: ISM results

Summary Multispectral imaging of planets opens a new class of studies of atmospheric and surface signatures. With ISM experiment on board Phobos, the first IR multispectral imaging of Mars from space was achieved. In viem of the Future Omega experiment on board Mars 94, we present here a sample of data and results (and problems) which can be achieved by studying spectral reflectance. Paper presented at the V Cosmic Physics National Conference, S. Miniato, November 27–30, 1990.     

Numerical and theoretical study of a Dual Mesh Method using finite volume schemes for two phase flow problems in porous media

Summary. In this paper we are interested in two phase flow problems in porous media. We use a Dual Mesh Method to discretize this problem with finite volume schemes. In a simplified case (elliptic - hyperbolic system) we prove the convergence of approximate solutions to the exact solutions. We use the Dual Mesh Method in physically complex problems (heterogeneous cases with non constant total mobility). We validate numerically the Dual Mesh Method on practical examples by computing error estimates for different test-cases. Received March 21, 1997 / Revised version received October 13, 1997