Finite element approximation of the Neumann eigenvalue problem in domains with multiple cracks

^** Email: belhach{at}poncelet.univ-metz.fr^*** Email: bucur{at}math.univ-metz.fr^**** Email: jmse{at}math.univ-metz.fr We study the Neumann–Laplacian eigenvalue problem in domainswith multiple cracks. We derive a mixed variational formulationwhich holds on the whole geometric domain (including the cracks)and implements efficient finite-element discretizations forthe computation of eigenvalues. Optimal error estimates aregiven and several numerical examples are presented, confirmingthe efficiency of the method. As applications, we numericallyinvestigate the behaviour of the low eigenvalues in domainswith a large number of cracks.     

Error estimates for stochastic differential games: the adverse stopping case

^** Email: frederic.bonnans{at}inria.fr^*** Email: stefania.maroso{at}inria.fr^**** Email: zidani{at}ensta.fr We obtain error bounds for monotone approximation schemes ofa particular Isaacs equation. This is an extension of the theoryfor estimating errors for the Hamilton–Jacobi–Bellmanequation. To obtain the upper error bound, we consider the ‘Krylovregularization’ of the Isaacs equation to build an approximatesub-solution of the scheme. To get the lower error bound, weextend the method of Barles & Jakobsen (2005, SIAM J. Numer.Anal.) which consists in introducing a switching system whosesolutions are local super-solutions of the Isaacs equation.     

A cell-centred finite-volume approximation for anisotropic diffusion operators on unstructured meshes in any space dimension

^** Email: eymard{at}math.univ-mlv.fr^*** Email: gallouet{at}cmi.univ-mrs.fr^**** Corresponding author. Email: herbin{at}cmi.univ-mrs.fr Finite-volume methods for problems involving second-order operatorswith full diffusion matrix can be used thanks to the definitionof a discrete gradient for piecewise constant functions on unstructuredmeshes satisfying an orthogonality condition. This discretegradient is shown to satisfy a strong convergence property forthe interpolation of regular functions, and a weak one for functionsbounded in a discrete H¹-norm. To highlight the importance ofboth properties, the convergence of the finite-volume schemefor a homogeneous Dirichlet problem with full diffusion matrixis proven, and an error estimate is provided. Numerical testsshow the actual accuracy of the method.     

On a Monte Carlo method for neutron transport criticality computations

^** Email: maire{at}univ-tln.fr^*** Email: denis.talay{at}sophia.inria.fr We give a stochastic representation of the principal eigenvalueof some homogeneous neutron transport operators. Our constructionis based upon the Feynman–Kac formula for integral transportequations, and uses probabilistic techniques only. We developa Monte Carlo method for criticality computations. We numericallytest this method on various homogeneous and inhomogeneous problems,and compare our results with those obtained by standard methods.     

A finite-volume scheme for dynamic reliability models

^** Email: christiane.cocozza{at}univ-mlv.fr^*** Email: robert.eymard{at}univ-mlv.fr^**** Email: sophie.mercier{at}univ-mlv.fr In a model arising in the dynamic reliability study of a system,the probability of the state of the system is completely describedby the Chapman–Kolmogorov equations, which are scalarlinear hyperbolic partial differential equations coupled bytheir right-hand side, the solution of which are probabilitymeasures. We propose in this paper a finite-volume scheme toapproximate these measures. We show, thanks to the proof ofthe tightness of the approximate solution, that the conservationof the probability mass leads to a compactness property. Theconvergence of the scheme is then obtained in the space of continuousfunctions with respect to the time variable, valued in the setof probability measures on [graphic: see PDF] . We finally show on a numerical example the accuracy and efficiencyof the approximation method.     

Finite time stability of differential inclusions

^** Email: emmanuel.moulay{at}ec-lille.fr^*** Email: wilfrid.perruquetti{at}ec-lille.fr In this paper, the problem of finite time stability is investigatedfor differential inclusion. Two sufficient conditions for finitetime stability, using a smooth Lyapunov function and a non-smoothone, are established. Then, the same idea is used to give twonecessary conditions. Examples are developed using the conceptof Krasovskii solutions for differential equation with discontinuousright-hand sides.     

The Helmholtz equation in a locally perturbed half-plane with passive boundary

^** Email: mduran{at}ing.puc.cl^*** Email: ignacio.muga{at}ucv.cl^**** Email: nedelec{at}cmapx.polytechnique.fr In this article, we study the existence and uniqueness of outgoingsolutions for the Helmholtz equation in locally perturbed half-planeswith passive boundary. We establish an explicit outgoing radiationcondition which is somewhat different from the usual Sommerfeld'sone due to the appearance of surface waves. We work with thehelp of Fourier analysis and a half-plane Green's function framework.This is an extended and detailed version of the previous articleDurán et al. (2005, The Helmholtz equation with impedancein a half-plane. C. R. Acad. Sci. Paris, Ser. I, 340, 483–488).     

A priori and a posteriori analysis of non-conforming finite elements with face penalty for advection-diffusion equations

^** Corresponding author. Email: l.elalaoui{at}imperial.ac.uk^*** Email: ern{at}cermics.enpc.fr^**** Email: erik.burman{at}epfl.ch We analyse a non-conforming finite-element method to approximateadvection–diffusion–reaction equations. The methodis stabilized by penalizing the jumps of the solution and thoseof its advective derivative across mesh interfaces. The a priorierror analysis leads to (quasi-)optimal estimates in the meshsize (sub-optimal by order in the L²-norm and optimal in thebroken graph norm for quasi-uniform meshes) keeping the Pécletnumber fixed. Then, we investigate a residual a posteriori errorestimator for the method. The estimator is semi-robust in thesense that it yields lower and upper bounds of the error whichdiffer by a factor equal at most to the square root of the Pécletnumber. Finally, to illustrate the theory we present numericalresults including adaptively generated meshes.     

Electrical coupling in proton exchange membrane fuel cell stacks: mathematical and computational modelling

^** Corresponding author. Email: wetton{at}math.ubc.ca^*** Email: Peter.Berg{at}uoit.ca^**** Email: caglara{at}uwgb.edu^***** Email: kpromisl{at}math.msu.edu^****** Email: jean.st-pierre{at}ballard.com A mathematical model describing the effects of electrical couplingof proton exchange membrane unit fuel cells through shared bipolarplates is developed. Here, the unit cells are described by simple,steady-state, 1D models appropriate for straight reactant gaschannel designs. A linear asymptotic version of the model isused to give analytic insight into the effect of the coupling,including estimates of the extent of the coupling in terms ofthe number of adjacent cells affected. An efficient numericalmethod is developed to solve the non-linear coupled system.Numerical results showing the effects on stack voltage due toa single cell with anomalous oxidant flow rate are given. Theeffects on stack performance due to end plate effects are alsogiven. It is shown that electrical coupling has a significanteffect on fuel cell performance.     

Post-processing with computable error bounds for the finite element approximation of a nonlinear heat conduction problem

Email: ain{at}mcs.le.ac.uk Email: D.Kelly{at}unsw.edu.au^* Email: I.Sloan{at}unsw.edu.au^** Email: swang{at}cs.curtin.edu.au It is shown how the finite element approximation of a nonlinearheat conduction problem may be post-processed to yield enhancedapproximations to the solution and the flux at any point inthe domain. Sharp computable bounds on the accuracy of the post-processedapproximations are derived. A criterion is identified for guidingadaptive refinements of the finite element discretization. Anumerical example is given illustrating the theoretical results.     

An LMI-based approach for robust stabilization of time delay systems containing saturating actuators

Email: elhaous_fati{at}yahoo.fr Corresponding author. Email: elh_tissir{at}yahoo.fr Received on September 8, 2005; Accepted on July 24, 2006 This paper deals with the problem of robust stabilization foruncertain systems with input saturation and time delay in thestate. The parameter uncertainties are time-varying and unknownbut are norm bounded. Sufficient conditions obtained via a linearmatrix inequality formulation are stated to guarantee the localstabilization. The method of synthesis consists in determiningsimultaneously a state feedback control law and an associateddomain of safe admissible states for which the stability ofthe closed-loop system is guaranteed when control saturationseffectively occur. Numerical examples are used to demonstratethe effectiveness of the proposed design technique.     

Reduced-basis output bound methods for parabolic problems

^** Email: rovas{at}uiuc.edu^*** Email: luc_machiels{at}mckinsey.com^**** Corresponding author. Email: maday{at}ann.jussieu.fr In this paper, we extend reduced-basis output bound methodsdeveloped earlier for elliptic problems, to problems describedby ‘parameterized parabolic’ partial differentialequations. The essential new ingredient and the novelty of thispaper consist in the presence of time in the formulation andsolution of the problem. First, without assuming a time discretization,a reduced-basis procedure is presented to ‘efficiently’compute accurate approximations to the solution of the parabolicproblem and ‘relevant’ outputs of interest. In addition,we develop an error estimation procedure to ‘a posteriorivalidate’ the accuracy of our output predictions. Second,using the discontinuous Galerkin method for the temporal discretization,the reduced-basis method and the output bound procedure areanalysed for the semi-discrete case. In both cases the reduced-basisis constructed by taking ‘snapshots’ of the solutionboth in time and in the parameters: in that sense the methodis close to Proper Orthogonal Decomposition (POD).     

A coupled multivalued model for ice streams and its numerical simulation

^** Email: nati{at}dma.uvigo.es^*** Email: durany{at}dma.uvigo.es^**** Email: anaisabel.munoz{at}urjc.es^***** Email: emanuele.schiavi{at}urjc.es^****** Email: carlosv{at}udc.es This paper deals with the numerical solution of a non-linearmodel describing a free-boundary problem arising in modern glaciology.Considering a shallow, viscous ice sheet flow along a soft,deformable bed, a coupled non-linear system of differentialequations can be obtained. Particularly, an obstacle problemis then deduced and solved in the framework of its complementarityformulation. We present the numerical solution of the resultingmultivalued system modelling the ice sheet non-Newtonian dynamicsdriven by the underlying drainage system. Our numerical resultsshow the existence of fast ice streams when positive wave-likeinitial conditions are considered. The solutions are numericallycomputed with a decoupling iterative method and finite-elementtechnique. A duality algorithm and a projected Gauss–Seidelmethod are the alternatives used to cope with the resultingvariational inequality while the explicit treatment, Newtonmethod or a duality method are proposed to deal with the non-linearsource term. Finally, the numerical solutions are physicallyinterpreted and some comparisons among the numerical methodsare then discussed.     

Solving parallel machines scheduling problems with sequence-dependent setup times using variable neighbourhood search

^** Corresponding author. Email: mahdi{at}dcc.ufmg.br^*** Email: martin{at}dcc.ufmg.br^**** Email: mateus{at}dcc.ufmg.br^***** Email: pardalos{at}ufl.edu Variable neighbourhood search (VNS) is a modern metaheuristicbased on systematic changes of the neighbourhood structure withina search to solve optimization problems. The aim of this paperis to propose and analyse a VNS algorithm to solve schedulingproblems with parallel machines and sequence-dependent setuptimes, which are of great importance on the industrial context.Three versions of a greedy randomized adaptive search procedurealgorithm are used to compare with the proposed VNS algorithmto highlight its advantages in terms of generality, qualityand speed for large instances.     

A descent modified Polak-Ribiere-Polyak conjugate gradient method and its global convergence

^** Email: zl606{at}tom.com^*** Corresponding author. Email: weijunzhou{at}126.com^**** Email: dhli{at}hnu.cn In this paper, we propose a modified Polak–Ribière–Polyak(PRP) conjugate gradient method. An attractive property of theproposed method is that the direction generated by the methodis always a descent direction for the objective function. Thisproperty is independent of the line search used. Moreover, ifexact line search is used, the method reduces to the ordinaryPRP method. Under appropriate conditions, we show that the modifiedPRP method with Armijo-type line search is globally convergent.We also present extensive preliminary numerical experimentsto show the efficiency of the proposed method.     

On the numerical solution of involutive ordinary differential systems: enhanced linear algebra

^** Email: jukka.tuomela{at}joensuu.fi^*** Corresponding author. Email: arponen{at}maths.warwick.ac.uk^**** Email: villesamuli.normi{at}joensuu.fi We analyse some Runge–Kutta type methods for computing1D integral manifolds, i.e. solutions to ordin-ary differentialequations and differential-algebraic equations. We show thatwe can compute the solutions which respect all the constraintsof the problem reliably and reasonably quickly. Moreover, weshow that the so-called impasse points are regular points inour approach and hence require no special attention.     

A stabilized finite-element method for the Stokes problem including element and edge residuals

^** Email: raraya{at}ing-mat.udec.cl^*** Corresponding author. Email: gbarrene{at}ing-mat.udec.cl^**** Email: valentin{at}lncc.br A new stabilized finite-element method is presented for theStokes problem. The method is of a Douglas–Wang type,and includes a positive jump term controlling the residual ofthe Cauchy stress tensor on the internal edges of the triangulation.A priori error estimates are obtained in the natural norms ofthe unknowns and an a posteriori error estimator is proposed,analysed and tested through numerical experiments.     

On a modified version of ILDM approach: asymptotic analysis based on integral manifolds

^** Email: vbykov{at}cs.bgu.ac.il^*** Email: goldfarb{at}cs.bgu.ac.il^**** Email: vladimir{at}bgumail.bgu.ac.il^***** Email: umaas{at}itt.mach.uni-karlsruhe.de Using the method of integral (invariant) manifolds, the intrinsiclow-dimensional manifolds (ILDM) method is analysed. This isa method for identifying invariant manifolds of a system's slowdynamics and has proven to be an efficient tool in modellingof laminar and turbulent combustion. It allows treating multi-scalesystems by revealing their hidden hierarchy and decomposingthe system dynamics into fast and slow motions. The performedanalysis shows that the original ILDM technique can be interpretedas one of the many possible realizations of the general framework,which is based on a special transformation of the original coordinatesin the state space. A modification of the ILDM is proposed basedon a new definition of the transformation matrix. The proposednumerical procedure is demonstrated on linear examples and highlynon-linear test problems of mathematical theory of combustionand demonstrates in some cases better performance with respectto the existing one.     

Supraconvergence of a finite difference scheme for solutions in Hs(0, L)

^** Email: silvia{at}mat.uc.pt^*** Email: ferreira{at}mat.uc.pt^**** Email: grigo{at}math.tu-berlin.de In this paper we study the convergence of a centred finite differencescheme on a non-uniform mesh for a 1D elliptic problem subjectto general boundary conditions. On a non-uniform mesh, the schemeis, in general, only first-order consistent. Nevertheless, weprove for s (1/2, 2] order O(h^s)-convergence of solution andgradient if the exact solution is in the Sobolev space H^1+s(0,L), i.e. the so-called supraconvergence of the method. It isshown that the scheme is equivalent to a fully discrete linearfinite-element method and the obtained convergence order isthen a superconvergence result for the gradient. Numerical examplesillustrate the performance of the method and support the convergenceresult.     

A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations

^** Email: jingtang{at}lsec.cc.ac.cn^*** Email: hermann{at}math.mun.ca In this paper we establish a posteriori error estimates forthe discontinuous Galerkin (DG) method applied to linear, semilinearand non-standard (non-linear) Volterra integro-differentialequations. We also present an analysis of the DG method withquadrature for the memory term. Numerical experiments basedon three integro-differential equations are used to illustratevarious aspects of the error analysis.