Hybridization and Postprocessing Techniques for Mixed Eigenfunctions

We introduce hybridization and postprocessing techniques for the Raviart–Thomas approximation of second-order elliptic eigenvalue problems. Hybridization reduces the Raviart–Thomas approximation to a condensed eigenproblem. The condensed eigenproblem is nonlinear, but smaller than the original mixed approximation. We derive multiple iterative algorithms for solving the condensed eigenproblem and examine their interrelationships and convergence rates. An element-by-element postprocessing technique to improve accuracy of computed eigenfunctions is also presented. We prove that a projection of the error in the eigenspace approximation by the mixed method (of any order) superconverges and that the postprocessed eigenfunction approximations converge faster for smooth eigenfunctions. Numerical experiments using a square and an L-shaped domain illustrate the theoretical results

Reconstructing initial data using observers: error analysis of the semi-discrete and fully discrete approximations

A new iterative algorithm for solving initial data inverse problems from partial observations has been recently proposed in Ramdani et al. (Automatica 46(10), 1616-1625, 2010 ). Based on the concept of observers (also called Luenberger observers), this algorithm covers a large class of abstract evolution PDE's. In this paper, we are concerned with the convergence analysis of this algorithm. More precisely, we provide a complete numerical analysis for semi-discrete (in space) and fully discrete approximations derived using finite elements in space and an implicit Euler method in time. The analysis is carried out for abstract Schrödinger and wave conservative systems with bounded observation (locally distributed)

HDG-NEFEM with Degree Adaptivity for Stokes Flows

This paper presents the first degree adaptive procedure able to directly use the geometry given by a CAD model. The technique uses a hybridisable discontinuous Galerkin discretisation combined with a NURBS-enhanced rationale, completely removing the uncertainty induced by a polynomial approximation of curved boundaries that is common within an isoparametric approach. The technique is compared against two strategies to perform degree adaptivity currently in use. This paper demonstrates, for the first time, that the most extended technique for degree adaptivity can easily lead to a non-reliable error estimator if no communication with CAD software is introduced whereas if the communication with the CAD is done, it results in a substantial computing time. The proposed technique encapsulates the CAD model in the simulation and is able to produce reliable error estimators irrespectively of the initial mesh used to start the adaptive process. Several numerical examples confirm the findings and demonstrate the superiority of the proposed technique. The paper also proposes a novel idea to test the implementation of high-order solvers where different degrees of approximation are used in different elements

Successive Measurements of Cosmic-Ray Antiproton Spectrum in a Positive Phase of the Solar Cycle

The energy spectrum of cosmic-ray antiprotons has been measured by BESS successively in 1993, 1995, 1997 and 1998. In total, 848 antiprotons were clearly identified in energy range 0.18 to 4.20 GeV. From these successive measurements of the antiproton spectrum at various solar activity, we discuss about the effect of the solar modulation and the origin of cosmic-ray antiprotons. Measured antiproton ratios were nearly identical during this period, and were consistent with a prediction taking the charge dependent solar modulation into account.Comment: 15 pages, 5 figure

Systems of Differential Algebraic Equations in Computational Electromagnetics

Starting from space-discretisation of Maxwell's equations, various classical formulations are proposed for the simulation of electromagnetic fields. They differ in the phenomena considered as well as in the variables chosen for discretisation. This contribution presents a literature survey of the most common approximations and formulations with a focus on their structural properties. The differential-algebraic character is discussed and quantified by the differential index concept

Article electronically published on November 16, 2007 COUPLING OF GENERAL LAGRANGIAN SYSTEMS

Abstract. This work is devoted to the coupling of two fluid models, such as two Euler systems in Lagrangian coordinates, at a fixed interface. We define coupling conditions which can be expressed in terms of continuity of some well chosen variables and then solve the coupled Riemann problem. In the present setting where the interface is characteristic, a particular choice of dependent variables which are transmitted ensures the overall conservativity. 1