Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Anomalous electric-field effect and glassy behaviour in granular aluminium thin films: electron glass?

We present a study of non-equilibrium phenomena observed in the electrical conductance of insulating granular aluminium thin films. An anomalous field effect and its slow relaxation are studied in some detail. The phenomenology is very similar to the one already observed in indium oxide. The origin of the phenomena is discussed. In granular systems, the present experiments can naturally be interpreted along two different lines. One relies on a slow polarisation in the dielectric surrounding the metallic islands. The other one relies on a purely electronic mechanism: the formation of an electron Coulomb glass in the granular metal. More selective experiments and/or quantitative predictions about the Coulomb glass properties are still needed to definitely distinguish between the two scenarios.     

Runge–Kutta Methods for Itô Stochastic Differential Equations with Scalar Noise

A general class of stochastic Runge–Kutta methods for Itô stochastic differential equation systems w.r.t. a one-dimensional Wiener process is introduced. The colored rooted tree analysis is applied to derive conditions for the coefficients of the stochastic Runge–Kutta method assuring convergence in the weak sense with a prescribed order. Some coefficients for new stochastic Runge–Kutta schemes of order two are calculated explicitly and a simulation study reveals their good performance.     

Second order Runge–Kutta methods for Stratonovich stochastic differential equations

The weak approximation of the solution of a system of Stratonovich stochastic differential equations with a m–dimensional Wiener process is studied. Therefore, a new class of stochastic Runge–Kutta methods is introduced. As the main novelty, the number of stages does not depend on the dimension m of the driving Wiener process which reduces the computational effort significantly. The colored rooted tree analysis due to the author is applied to determine order conditions for the new stochastic Runge–Kutta methods assuring convergence with order two in the weak sense. Further, some coefficients for second order stochastic Runge–Kutta schemes are calculated explicitly. AMS subject classification (2000)  65C30, 65L06, 60H35, 60H10     

Electron-confined LO-phonon scattering in GaAs-Al0.45Ga0.55As superlattice

We develop a theoretical model to the scattering time due to the electron-confined LO-phonon in GaAs-Al_xGa_1-xAs superlattice taking into account the sub-band parabolicity. Using the new analytic wave function of electron miniband conduction of superlattice and a reformulation slab model for the confined LO-phonon modes, an expression for the electron-confined LO-phonon scattering time is obtained. In solving numerically a partial differential equation for the phonon generation rate, our results show that forx = 0.45, the LO-phonon in superlattice changes from a bulk-like propagating mode to a confined mode. The dispersion of the relaxation time due to the emission of confined LO-phonons depends strongly on the total energy.     

A Distributional Approach to Multiple Stochastic Integrals and Transformations of the Poisson Measure

We study a class of ‘nonpoissonian’ transformations of the configuration space and the corresponding transformations of the Poisson measure. For some class of Poisson measures we find conditions which are sufficient for the transformed measure (which in general is nonpoissonian) to be absolutely continuous with respect to the initial Poisson measure and get the expression for the corresponding Radon–Nikodym derivative. To solve this problem we use a distributional approach to Poisson multiple stochastic integrals.     

A Relaxation Scheme for Solving the Boltzmann Equation Based on the Chapman-Enskog Expansion

Abstract In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jinand Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weaklyparabolic, has a linearly hyperbolic convection part, and is endowed with a generalized eotropy inequality. Itagrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system forthe Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme,and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and bythe extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments showthat the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equationobtained by the DSMC, for a range of Mach numbers for hypersonic flows, th     

Converged properties of clean metal surfaces by all-electron first-principles calculations

All-electron full-potential linearized augmented plane-wave calculations of the surface energy, work function, and interlayer spacings of close-packed metal surfaces are presented, in particular, for the free-electron-like metal surfaces, Mg(0 0 0 1) and Al(1 1 1), and for the transition metal surfaces, Ti(0 0 0 1), Cu(1 1 1), Pd(1 1 1), and Pt(1 1 1). We investigate the convergence of the surface energy as a function of the number of layers in the slab, using the Cu(1 1 1) surface as an example. The results show that the surface energy, as obtained using total energies of the slab and bulk from separate calculations, converges well with respect to the number of layers in the slab. Obviously, it is necessary that bulk and surface calculations are performed with the same high accuracy. Furthermore, we discuss the performance of the local-density and generalized gradient approximations for the exchange-correlation functional in describing the various surface properties.     

带ARMA(1,1)条件异方差相关的随机波动模型的MCMC算法

我们首先提出了一个带ARMA(1,1)条件异方差相关的随机波动模型,它是基本的随机波动模型的一个自然的推广.进一步,对于这一新模型,我们给出了一个马尔可夫链蒙特卡罗(M CM C)算法.最后,利用该模型的模拟数据,展示了M CM C算法在这种模型中的应用.