A technique to construct symmetric variable-stepsize linear multistep methods for second-order systems

Some previous works show that symmetric fixed- and variable-stepsize linear multistep methods for second-order systems which do not have any parasitic root in their first characteristic polynomial give rise to a slow error growth with time when integrating reversible systems. In this paper, we give a technique to construct variable-stepsize symmetric methods from their fixed-stepsize counterparts, in such a way that the former have the same order as the latter. The order and symmetry of the integrators obtained is proved independently of the order of the underlying fixed-stepsize integrators. As this technique looks for efficiency, we concentrate on explicit linear multistep methods, which just make one function evaluation per step, and we offer some numerical comparisons with other one-step adaptive methods which also show a good long-term behaviour.

     

Heuristic Rejection in Interval Global Optimization

Based on the investigation carried out in Ref. 1, this paper incorporates new studies about the properties of inclusion functions on subintervals while a branch-and-bound algorithm is solving global optimization problems. It is found that the relative place of the global minimum value within the inclusion function value of the objective function at the current interval indicates mostly whether the given interval is close to a minimizer point. This information is used in a heuristic interval rejection rule that can save a considerable amount of computation. Illustrative examples are discussed and an extended numerical study shows the advantages of the new approach.     

Measurement of the dynamic elastic constants of short isotropic cylinders

A method based on a single test is proposed to characterize the elasticity of an isotropic homogeneous material in the shape of a cylinder of any slenderness (length-diameter) ratio. Firstly, the Rayleigh-Ritz method is used to determine the natural frequencies of the cylinders vibrating axisymmetrically. The study is focused on cylindrical samples with diameter and length of similar magnitude so that the shear modulus and the Poisson ratio can be calculated simultaneously. Subsequently, the theoretical results for cylinders of slenderness ratio between 0.1 and 3 are analyzed in order to obtain the data required to determine the elastic constants from one of the two lowest measured natural frequencies and their quotient. The analysis of the results demonstrates that any slenderness ratio is useful in the calculation of the elastic constants, although in some cases the third natural frequency should be used. Furthermore, the influence of the length-diameter quotient on the sensitivity of the method is analyzed by evaluating the systematic uncertainties for both dynamic elastic constants. Finally, the method is experimentally tested by characterizing two steel cylinders with slenderness ratios 0.1 and 1, respectively. The results demonstrate that uncertainties for both Poisson ratio and the shear modulus are smaller when the slenderness ratio is 1.     

Magnetic phases in Mn<Subscript>1</Subscript>_<Subscript>x</Subscript>Fe<Subscript>x</Subscript>WO<Subscript>4</Subscript> studied by neutron powder diffraction

The magnetic structures of Mn_1-xFe_xWO₄ with x = 0.0, 0.16, 0.21, 0.225, 0.232, 0.24, 0.27, 0.29, and 1.0 were refined from neutron powder diffraction data. The magnetic phase diagram could be completed in the coexistence range of different magnetic structures up to x = 0.29. For the magnetic state at 1.5 K a commensurate antiferromagnetic structure with a propagation vector = (±1/4, 1/2, 1/2) was found for x ⩽ 0.22 while the magnetic spins order with = (1/2, 0, 0) for x ≥ 0.22. In the latter phase, additionally, weak magnetic reflections indexed to an incommensurate ordering with = (- 0.214, 1/2, 0.457) occur in the diffraction pattern up to x = 0.29 indicating the occurence of a reentrant phase. For 0.12 ⩽ x ⩽ 0.29 the low temperature phases are separated from a magnetic high temperature phase showing only magnetic reflections indexed to a spin arrangement with = (1/2, 0, 0). The magnetic phase diagram is discussed qualitatively considering random superexchange between the statistically distributed Mn²⁺- and Fe²⁺-ions in the coexistence range 0.12 ⩽ x ⩽ 0.29 of different magnetic structures related to those of pure MnWO₄ and FeWO₄. Received 9 October 2002 Published online 14 March 2003     

Asymptotically balanced schemes for non-homogeneous hyperbolic systems – application to the Shallow Water equations

In this work we introduce a class of balanced numerical schemes, up to second order, for the solution of general non-homogeneous hyperbolic systems of conservation laws. We give a general technique to build such schemes. We also prove that they balance up to second order a large class of steady solutions in the whole domain but some subset whose measure tends to zero as the grid size decreases to zero. We finally present an application to Shallow Water equations that exhibit the good performances of some of the schemes introduced. To cite this article: T. Chacón Rebollo et al., C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

Equivalence between mixed finite element and multi-point finite volume methods

We consider the lowest-order Raviart–Thomas mixed finite element method for elliptic problems on simplicial meshes in two or three space dimensions. This method produces saddle-point type problems for scalar and flux unknowns. We show how to easily eliminate the flux unknowns, which implies an equivalence between this method and a particular multi-point finite volume scheme, without any approximate numerical integration. We describe the stencil of the final matrix and give sufficient conditions for its symmetry and positive definiteness. We present a numerical example illustrating the performance of the proposed method. To cite this article: M. Vohralík, C. R. Acad. Sci. Paris, Ser. I 339 (2004).     

Hysteresis and Eddy-current Relaxation in SiFe Materials

The rate-dependent effects in metallic ferromagnetic materials with magnetization processes attributed to domain-walls motion were studied. The experimental method for eddy-current relaxation time determination was proposed. Some discrepancies between the data and theoretical predictions were observed and discussed. The origin of these discrepancies is accounted for by the decrease of viscous-type relaxation process as the rate of change of averaged magnetization is approaching zero.     

The Gross-Neveu model at finite temperature at next-to-leading order in the 1/N expansion

We present new results on the Gross-Neveu model at finite temperature and at next-to-leading order in the 1/N expansion. In particular, a new expression is obtained for the effective potential which is explicitly invariant under renormalization group transformations. The model is used as a playground to investigate various features of field theory at finite temperature. For example we verify that, as expected from general arguments, the cancellation of ultraviolet divergences takes place at finite temperature without the need for introducing counterterms beyond those of zero temperature. As well known, the discrete chiral symmetry of the (1+1)-dimensional model is spontaneously broken at zero temperature and restored, in leading order, at some temperature T_c; we find that the 1/N approximation breaks down for temperatures below T_c: as the temperature increases, the fluctuations become eventually too large to be treated as corrections, and a Landau pole invalidates the calculation of the effective potential in the vicinity of its minimum. Beyond T_c, the 1/N expansion becomes again regular: it predicts that in leading order the system behaves as a free gas of massless fermions and that, at the next-to-leading order, it remains weakly interacting. In the limit of large temperature, the pressure coincides with that given by perturbation theory with a coupling constant defined at a scale of the order of the temperature, as expected from asymptotic freedom.     

A theoretical analysis of the TiO2/Sn doped (1 1 0) surface properties

We have used the periodic quantum-mechanical method with density functional theory at the B3LYP level in order to study TiO₂/Sn doped (1 1 0) surfaces and have investigated the structural, electronic and energy band properties of these oxides. Our calculated relaxation directions for TiO₂ is the experimental one and is also in agreement with other theoretical results. We also observe for the doped systems relaxation of lattice positions of the atoms. Modification of Sn, O and Ti charges depend on the planes and positions of the substituted atoms. Doping can modify the Fermi levels, energy gaps as well as the localization and composition of both valence and conduction band main components. Doping can also modify the chemical, electronic and optical properties of these oxides surfaces increasing their suitability for use as gas sensors and optoelectronic devices.