JET: Preparing the future in fusion

JET (Joint European Torus) is the largest tokamak in the world and the only fusion facility able to operate with Tritium, the fusion fuel, and Beryllium, the ITER first wall material. JET also features the most complete remote handling equipment for invessel maintenance. As a multinational research center, JET provides logistic experience in preparing for operation of the global facility, tokamak ITER.Experiments on JET are focused on ITER-relevant studies, in particular on detailing the operational scenarios (EL My H-modes and advanced regimes), on enhancing the heating systems, on developing diagnostics for burning plasmas etc. Pioneering real-time control techniques have been implemented that maximize performance and minimize internal disturbances of JET plasmas. In helium plasmas, ion cyclotron heating (ICRH) created fast α-particles, mimicking their populations in future burning plasmas. The recent successful Trace Tritium campaign provided important new data on fuel transport. Current enhancements on JET include a new ITER-like ELM-resilient high power ICRH antenna (7 MW) and over twenty new diagnostics that will further extend the JET scientific capabilities and push the facility even closer to the ITER parameters.A special mention is given to the involvement of the fusion experts from Association EURATOM-IPP.CR, who have been actively participating in the collective use of JET facility for more than three years.     

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.