Study of oxygen vacancy ordering in mullite at high temperatures

High temperature X-ray diffraction and quenching experiments of mullite single crystals with Al₂O₃:SiO₂ ratio 2:1 have been performed to investigate the stability of the oxygen vacancy ordering close to the melting point of mullite. The experiments show that the structure of mullite exhibits an extremely stable, temperature-independent incommensurate modulation. Inspection of satellite reflections at different temperatures leads to the conclusion that the ordering scheme of oxygen vacancies after the crystallization of mullite persists to the melting point and does not show any disordering effects. The experimental results are in agreement with former theoretical calculations using a statistical mechanics approach which yield the critical temperature T_c > 3000°C.     

Numerical transformation methods: a constructive approach

Within group invariance theory we consider a constructive approach in order to define numerical transformation methods. These are initial-value methods for the solution of boundary value problems governed by ordinary differential equations. Here we consider the class of free boundary value problems governed by the most general second-order equation in normal form. For this class of problems the main theorem is concerned with the definition of an iterative transformation method. The definition of a noniterative method, applicable to a subclass of the original class of problems, follows as a corollary. Therefore, the proposed constructive approach allows us to establish a unifying framework for noniterative and iterative transformation methods.     

Parallel globally adaptive quadrature on the KSR-1

New algorithms for parallel one-dimensional globally adaptive quadrature are developed. The algorithms are implemented on a Kendall Square Research KSR-1 parallel computer and numerical results are presented. The most successful algorithm gives significant speedups on a range of hard problems, including ones with singular integrands. Both authors acknowledge the support of the EEC Esprit Basic Research Action Programme, Project 6634 (APPARC). The second author acknowledges the support of the NATO Collaborative Research Grant 920037.     

Monitoring cooperative equilibria in a stochastic differential game

This paper deals with a class of equilibria which are based on the use of memory strategies in the context of continuous-time stochastic differential games. In order to get interpretable results, we will focus the study on a stochastic differential game model of the exploitation of one species of fish by two competing fisheries. We explore the possibility of defining a so-called cooperative equilibrium, which will implement a fishing agreement. In order to obtain that equilibrium, one defines a monitoring variable and an associated retaliation scheme. Depending on the value of the monitoring variable, which provides some evidence of a deviation from the agreement, the probability increases that the mode of a game will change from a cooperative to a punitive one. Both the monitoring variable and the parameters of this jump process are design elements of the cooperative equilibrium. A cooperative equilibrium designed in this way is a solution concept for a switching diffusion game. We solve that game using the sufficient conditions for a feedback equilibrium which are given by a set of coupled HJB equations. A numerical analysis, approximating the solution of the HJB equations through an associated Markov game, enables us to show that there exist cooperative equilibria which dominate the classical feedback Nash equilibrium of the original diffusion game model.This research was supported by FNRS-Switzerland, NSERC-Canada, FCAR-Quebec.     

Numerical solution of eigenvalue problems for Hamiltonian systems

This paper discusses the numerical solution of eigenvalue problems for Hamiltonian systems of ordinary differential equations. Two new codes are presented which incorporate the algorithms described here; to the best of the author’s knowledge, these are the first codes capable of solving numerically such general eigenvalue problems. One of these implements a new new method of solving a differential equation whose solution is a unitary matrix. Both codes are fully documented and are written inPfort-verifiedFortran 77, and will be available in netlib/aicm/sl11f and netlib/aicm/sl12f.     

DQAINF: an algorithm for automatic integration of infinite oscillating tails

We describe an automatic quadrature routine which is specifically designed for real functions having a certain type of infinite oscillating tails. The algorithm is designed to integrate a vector function over an infinite interval. A FORTRAN implementation of the algorithm is included.The algorithm combines an adaptive subdivision strategy with extrapolation and requires that the decay of all the functions in the vector is the same. The algorithm is based on the assumption that the oscillating behavior is due to a periodic function with the property that it changes sign when evaluated at points of distance half a period. We assume that this period is known. The algorithm offers a choice of three different quasi-linear extrapolation procedures, namely the Euler transformation and two modifications of this transformation.This work was supported by The Norwegian Research Council for Science and the Humanities.     

A parallel method for finding the global minimum of univariate functions

We describe a new parallel method for solving global optimization problems. The formulation of the decision rules of this method is presented. We examine convergence conditions of the proposed algorithm and establish conditions which guarantee a considerable speedup with respect to the sequential version of the algorithm. We also present some numerical experiments executed on Alliant FX/80 for one class of multiextremal functions.The authors are greatly indebted to R. G. Strongin who stimulated the fulfillment of this research. They also would like to thank the anonymous referees for their useful suggestions.The research of the first author was partially supported by Grant 9494/NC/89 from the Italian Government under the Italian-Soviet Agreement about the Cultural and Scientific Exchange in 1990–1991. He thanks the Systems Department, University of Calabria, where he was a Visitor.     

On the distribution of the number of packets in the fluid flow approximation of packet arrival streams

Fluid flow approximations are widely used for approximating models of communication systems where packet arrival streams are generated in a regular manner over certain intervals (constant rate). The appropriate mathematical model for describing those bursty arrival streams in the fluid flow framework are the well-known Markov modulated rate processes (MMRP). The paper deals with the distribution of the numberN(t) of packets in the interval [0,t] of MMRP. For two-state MMRPs and their superpositions we derive formulas for the distribution ofN(t) and its density. Further we give asymptotic results. The presented numerical results and simulation studies illustrate the goodness of the fluid flow approximation and show that the proposed numerical algorithms work well even in the case of multiplexing a large number of burst silence sources.This work was partially supported by a grant from the Deutsche Bundespost TELEKOM.     

离散Lyapunov矩阵方程X—AXB=C的一种数值解法

本文提供了求解离散Lyapunov矩阵方程的一种数值解法。首先讨论了系数矩阵为三对角形矩阵的公式解法,然后通过相似变换,将该方法推广到一般情形。     

一次强沙尘暴长距离输送的数值模拟研究

利用HYSPLIT-4传输,扩散模式对2003年4月9日发生在我国西北的强沙尘暴进行长距离输送的数值模拟研究,并进行轨迹和浓度分析.研究结果表明:(1)2003年4月9日强沙尘暴过程的沙尘来源并非某单一地区;(2)此次沙尘暴长距离输送过程中,主要有两个高浓度中心向东传输.一个在北纬35度附近,主要在1 500～3000m高空传输,另一个在北纬40度附近,且主要在1 500m以下传输.后期两者合二为一.