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.     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

The Nature of Electronic Charge

Advances in gauge theories and unified theories have not thrown light on the meaning of electron. The problem of the origin of electronic charge is made precise, new insights gained from Weyl space are summarized, and the origin of charge in terms of fractional spin is suggested. A new perspective on the abelian Chern-Simons theory is presented to explain charge.     

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 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     

Homogeneous immunoassays

Applied Biochemistry and Biotechnology -     

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.     

磁电子学器件应用原理

本文介绍几种重要的磁电子器件的基本结构和工作原理,包括巨磁电阻与隧穿磁电阻传感器、巨磁电阻隔离器、巨磁电阻与隧穿磁电阻硬盘读出磁头、磁电阻随机存取存储器、自旋转移磁化反转与微波振荡器。自旋晶体管作为未来磁电子学或自旋电子学时代的基本元素,目前大都还处在概念型阶段,本文也将对几种自旋晶体管的大致原理作简要介绍。     

Lagrangean heuristics combined with reoptimization for the 0-1 bidimensional knapsack problem

First, this paper deals with lagrangean heuristics for the 0-1 bidimensional knapsack problem. A projected subgradient algorithm is performed for solving a lagrangean dual of the problem, to improve the convergence of the classical subgradient algorithm. Secondly, a local search is introduced to improve the lower bound on the value of the biknapsack produced by lagrangean heuristics. Thirdly, a variable fixing phase is embedded in the process. Finally, the sequence of 0-1 one-dimensional knapsack instances obtained from the algorithm are solved by using reoptimization techniques in order to reduce the total computational time effort. Computational results are presented.