The Z-transform and linear multistep stability

^** Email: mapjjc{at}maths.bath.ac.uk^*** Corresponding author. Email: ath{at}maths.bath.ac.uk^**** Email: hl{at}maths.bath.ac.uk This paper makes systematic use of control-theoretic methodssuch as the -transform, small-gain theorems and frequency-domainstability criteria in the analysis of the stability behaviourof linear multistep methods. Some of the results in Nevanlinna'swork are recovered and a number of new boundedness and asymptoticproperties of solutions of numerical schemes are obtained. Inparticular, we give a careful and detailed analysis of the nonlinearstability properties of strictly zero-stable methods.     

An a posteriori error indicator for discontinuous Galerkin discretizations of H(curl)-elliptic partial differential equations

^** Email: paul.houston{at}nottingham.ac.uk^*** Corresponding author. Email: ilaria.perugia{at}unipv.it^**** Email: schoetzau{at}math.ubc.ca We introduce a residual-based a posteriori error indicator fordiscontinuous Galerkin discretizations of H(curl; )-ellipticboundary value problems that arise in eddy current models. Weshow that the indicator is both reliable and efficient withrespect to the approximation error measured in terms of a naturalenergy norm. We validate the performance of the indicator withinan adaptive mesh refinement procedure and show its asymptoticexactness for a range of test problems.     

A unified framework for the numerical solution of general quadratic matrix equations

^** Email: vassilios.tsachouridis{at}ieee.org^*** Email: N.karcanias{at}city.ac.uk^**** Email: ixp{at}le.ac.uk Algebraic quadratic equations are special cases of a singlegeneralized algebraic quadratic matrix equation (GQME). Thispaper focuses on the numerical solution of the GQME using probability-1homotopy methods. A synoptic review of these methods and theirapplication to algebraic matrix equations is provided as background.A large variety of analysis and design problems in systems andcontrol are reported as special cases of the presented frameworkand some of them are illustrated via numerical examples fromthe literature.     

Alternate strip-based substructuring algorithms for elliptic PDEs in two dimensions

^** Email: angela.mihai{at}strath.ac.uk^*** Email: alan.craig{at}durham.ac.uk The alternate strip-based substructuring algorithms are efficientpreconditioning techniques for the discrete systems which arisefrom the finite-element approximation of symmetric ellipticboundary-value problems in 2D Euclidean spaces. The new approachis based on alternate decomposition of the given domain intoa finite number of strips. Each strip is a union of non-overlappingsubdomains and the global interface between subdomains is partitionedas a union of edges between strips and edges between subdomainsthat belong to the same strip. Both scalability and efficiencyare achieved by alternating the direction of the strips. Thisapproach generates algorithms in two stages and allows the useof a two-grid V cycle. Numerical estimates illustrate the behaviourof the new domain decomposition techniques.     

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

A Lobatto interpolation grid over the triangle

^** Email: m.blyth{at}uea.ac.uk^*** Email: cpozrikidis{at}ucsd.edu A sequence of increasingly refined interpolation grids overthe triangle is proposed, with the goal of achieving uniformconvergence and ensuring high interpolation accuracy. The numberof interpolation nodes, N, corresponds to a complete mth-orderpolynomial expansion with respect to the triangle barycentriccoordinates, which arises by the horizontal truncation of thePascal triangle. The proposed grid is generated by deployingLobatto interpolation nodes along the three edges of the triangle,and then computing interior nodes by averaged intersectionsto achieve three-fold rotational symmetry. Numerical computationsshow that the Lebesgue constant and interpolation accuracy ofthe proposed grid compares favorably with those of the best-knowngrids consisting of the Fekete points. Integration weights correspondingto the set of Lobatto triangle base points are tabulated.     

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.     

An explicit expression for the surface-impedance tensor of a compressible monoclinic material in a state of plane strain

^** Email: y.fu{at}keele.ac.uk^*** Email: d.w.brookes{at}keele.ac.uk We present an explicit expression for the surface-impedancetensor associated with a compressible monoclinic elastic materialin a state of plane strain. Among a wide range of applications,such an explicit expression can, for instance, be used to writedown an explicit secular equation for the surface or interfacialwave speed, and an explicit bifurcation condition for surfaceor interfacial wrinkling. A previously known quartic secularequation for the surface-wave speed is recovered as a by-product.     

The multifacility maximin planar location problem with facility interaction

^** Email: s.salhi{at}kent.ac.uk Two branch-and-bound algorithms are proposed to optimally solvethe maximin formulation for locating p facilities in the plane.Tight upper and lower bounds are constructed and suitable methodsof guiding the search developed. To enhance the method, efficientmeasures for identifying specific squares for subdivision aresuggested. The proposed algorithms are evaluated on a set ofrandomly generated problems of up to five facilities and 120nodes.     

Optimal allocation of reliability tasks to mitigate faults during system development

^** Email: walter.johnston{at}baesystems.com^*** Email: j.quigley{at}strath.ac.uk^**** Email: lesley.walls{at}strath.ac.uk This paper considers reliability planning for a concept designfor a new system where a portfolio of possible reliability developmenttasks exists; the goal is to find a selection and sequence oftasks to achieve reliability targets subject to time constraintsat minimal cost. This is non-trivial given that each task potentiallycan expose several different weaknesses and each weakness potentiallycan be exposed by several different tasks. We use a Bayesianpoint process model to estimate the system reliability. Theprior distribution maps to a fault register and relates directlyto a set of potential engineering modifications.The likely impactof each task can be assessed using the point process model.An integer programming approach is used to sequence and scheduletasks under the constraint that contractual reliability requirementsmust be met. An illustrative example is provided and an extensionto system availability is proposed.     

On the formation of acceleration and reaction-diffusion wavefronts in autocatalytic-type reaction-diffusion systems

^** Email: d.j.needham{at}reading.ac.uk We consider generalization of the theory for the evolution ofreaction–diffusion and accelerating wavefronts in KPP-typesystems as developed in Needham (2004, Proc. R. Soc. Lond. A,460, 1921–1934) (DN). These generalizations allow forthe removal of a number of technical restrictions imposed inthe paper of DN.     

Solving the continuous space p-centre problem: planning application issues

^** Email: wei.97{at}osu.edu^*** Email: murray.308{at}osu.edu^**** Email: xiao.37{at}osu.edu The Voronoi diagram heuristic has been proposed for solvingthe p-centre problem in continuous space. However, importantassumptions underlie this heuristic and may be problematic forpractical applications. These simplifying assumptions includeuniformly distributed demand, representing a region as a rectangle;analysis of a simple Voronoi polygon in solving associated one-centreproblems and no restrictions on potential facility locations.In this paper, we explore the complexity of solving the continuousspace p-centre problem in location planning. Considering theissue of solution space feasibility, we present a spatiallyrestricted version of this problem and propose methods for solvingit heuristically. Theoretical and empirical results are provided.     

An adaptive boundary element method for the exterior Stokes problem in three dimensions

^** Email: vjervin{at}clemson.edu^*** Email: norbert.heuer{at}brunel.ac.uk We present an adaptive refinement strategy for the h-versionof the boundary element method with weakly singular operatorson surfaces. The model problem deals with the exterior Stokesproblem, and thus considers vector functions. Our error indicatorsare computed by local projections onto 1D subspaces definedby mesh refinement. These indicators measure the error separatelyfor the vector components and allow for component-independentadaption. Assuming a saturation condition, the indicators giverise to an efficient and reliable error estimator. Also we describehow to deal with meshes containing quadrilaterals which arenot shape regular. The theoretical results are underlined bynumerical experiments. To justify the saturation assumption,in the Appendix we prove optimal lower a priori error estimatesfor edge singularities on uniform and graded meshes.     

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.     

Improving the accuracy of the mini-element approximation to Navier-Stokes equations

^** Email: blanca{at}imati.cnr.it^*** Email: frutos{at}mac.cie.uva.es^**** Corresponding author. Email: julia.novo{at}uam.es A technique to improve the accuracy of the mini-element approximationto incompressible the Navier–Stokes equations is introduced.Once the mini-element approximation has been computed at a fixedtime, the linear part of this approximation is postprocessedby solving a discrete Stokes problem. The bubble functions neededto stabilize the approximation to the Navier–Stokes equationsare not used at the postprocessing step. This postprocessingprocedure allows us to increase by one unit (up to a logarithmicterm) the H¹ norm rate of convergence of the velocity and correspondinglythe L² norm of the pressure. An error analysis of the algorithmis performed.     

Non-variational approximation of discrete eigenvalues of self-adjoint operators

^** Email: L.Boulton{at}ma.hw.ac.uk We establish sufficient conditions for approximation of discreteeigenvalues of self-adjoint operators in the second-order projectionmethod suggested recently in Levitin & Shargorodsky (2004,Spectral pollution and second order relative spectra for self-adjointoperators. IMA J. Numer. Anal., 24, 393–416). We findfairly explicit estimates for the eigenvalue error and studyin detail two concrete model examples. Our results show thatsecond-order projection strategies not only are universallypollution free but also achieve approximation under naturalconditions on the discretising basis.     

Differentiated service delivery in cooperative IP-based broadcast and mobile telecommunications networks

Email: angelo.centonza{at}brunel.ac.uk Corresponding author. Email: thomas.owens{at}brunel.ac.uk Email: john.cosmas{at}brunel.ac.uk^¶ Email: y.h.song{at}brunel.ac.uk Received on 26 November 2005. Accepted on 20 March 2007. The evolution of wireless systems has led recently to the deploymentof cooperative network infrastructures where networks basedon different technologies cooperate together to offer innovativeservices which the networks individually could not offer. Cooperativenetwork infrastructures are hybrid systems. Examples of hybridsystems are already in use in areas as diverse power systemsand locomotion systems. It is very important to apply efficientsystem management techniques to hybrid systems that produceadvantageous business case scenarios for each participatingsystem and efficient use of the available resources. In thispaper, scenarios are investigated where cooperation betweenan integrated project (IP)-based broadcast network and an IP-basedmobile telecommunications network has the potential to allowservices to be provided to mobile users which would not be economicalto offer over either an IP-based broadcast network or an IP-basedmobile telecommunications network alone. The paper presentsa novel technique for determining which network to use to deliversuch services at a given point in time. The application of thistechnique in appropriate scenarios has the potential to generateadditional income for both the IP-based broadcast network operatorand the IP-based mobile telecommunications network operatorin an infrastructure where the operators cooperate to offerinnovative services. The paper explains the construction ofutility functions for cooperative IP-based broadcast and mobiletelecommunications networks. These utility functions are thenused to provide results that enable the efficiency of the managementof the network to be assessed in terms of the utility volumegenerated by the innovative services provided.     

The cyclic Barzilai--Borwein method for unconstrained optimization

^** Email: dyh{at}lsec.cc.ac.cn^*** Email: hager{at}math.ufl.edu^**** Email: klaus.schittkowski{at}uni-bayreuth.de^***** Email: hzhang{at}math.ufl.edu In the cyclic Barzilai–Borwein (CBB) method, the sameBarzilai–Borwein (BB) stepsize is reused for m consecutiveiterations. It is proved that CBB is locally linearly convergentat a local minimizer with positive definite Hessian. Numericalevidence indicates that when m > n/2 3, where n is the problemdimension, CBB is locally superlinearly convergent. In the specialcase m = 3 and n = 2, it is proved that the convergence rateis no better than linear, in general. An implementation of theCBB method, called adaptive cyclic Barzilai–Borwein (ACBB),combines a non-monotone line search and an adaptive choice forthe cycle length m. In numerical experiments using the CUTErtest problem library, ACBB performs better than the existingBB gradient algorithm, while it is competitive with the well-knownPRP+ conjugate gradient algorithm.     

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.