Finite-element methods for singularly perturbed high-order elliptic two-point boundary value problems. II: convection-diffusion-type problems

e-mail address: stynes{at}bureau.ucc.ie We consider singularly perturbed high-order elliptic two-pointboundary value problems of convection-diffusion type. Undersuitable hypotheses, the coercivity of the associated bilinearform is proved and a representation result for the solutionsof such problems is given. A family of Galerkin finite-elementmethods based on piecewise polynomial test/trial functions ona Shishkin mesh is constructed and proved to be convergent,uniformly in the perturbation parameter, in energy and Wnorms. Numerical results are presentedfor a second-order problem and fourth-order problems.     

Convergence of Methods for the Numerical Solution of the Korteweg--de Vries Equation

We prove that a family of methods for the numerical solutionof the Korteweg—de Vries equation is convergent. Thisfamily includes as particular cases some known finite differenceand finite element schemes. It is also found that the stabilityproperties of the methods vary significantly with the treatmentof the non-linear term. On leave from Shanghai University of Science and Technology,Shanghai, China.     

Long-time behaviour of arbitrary order continuous time Galerkin schemes for some one-dimensional phase transition problems

The long-time behaviour of continuous time Galerkin (CTG) approximationsof some well-known one-dimensional non-linear evolution problemswhich model phase transitions are analyzed. These numericalschemes are fully discrete and of arbitrary order. Partially supported by the Army Research Office through grant28535-MA. Partially supported by the ONR Contract No N00014-90-J-1238.     

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.     

Estimating vaccine coverage by using computer algebra

Email: altmannd{at}rki.deEmail: altmann{at}mathematik.hu-berlin.de The approach of N Gay for estimating the coverage of a multivalentvaccine from antibody prevalence data in certain age cohortsis complemented by using computer aided elimination theory ofvariables. Hereby, Gay's usage of numerical approximation canbe replaced by exact formulae which are surprisingly nice, too.     

A Simple Proof for a Spectral Factorization Theorem

It is shown using simple methods that the transform Z(z), z C, of the coefficient sequence of the Wold decomposition ofany full-rank wide-sense stationary purely non-deterministicstochastic process satisfies (i) Z(z) H² (D) and (ii) Z^–1(z) H(D). Further it is shown that all spectral factors satisfying(i) and (ii) are equal up to right multiplication by orthogonalmatrices, and that among these the normalized (Z(0) =I) spectralfactors are equal to the transform of the Wold decomposition.An elementary proof of Youla's Theorem is then given togetherwith a simple proof that the rows of a Cholesky factor of abanded block Toeplitz matrix converge to the coefficients ofa stable matrix polynomial. Computer and Automation Institute of the Hungarian Academyof Sciences, Budapes, Hungary.     

Numerical stabilization of polynomial and matrix

Email: jdhan{at}sia.cn Email: zhjiang{at}sia.cn Corresponding author Email: nyy{at}sia.cn Received on March 6, 2006; Accepted on September 4, 2006 In this paper, the conception of numerical stabilization, whichis related to mantissa digits of computer and dimensions ofsystem, is described; and several strategies for the numericalstabilization of polynomial and matrix are presented.     

On the Structure of the Moving Finite-element Equations

Address from 1st April 1985, School of Mathematics, Universityof Bristol, University Walk, Bristol BS8 1TW. The morning finite-element method for evolutionary partial differentialequations leads to a coupled non-linear system of ordinary differentialequations in time, with a coefficien matrix A, say, for thetime derivaties, We show for linear elements in any number ofdimensions, A can be written in the form M^TCM, where the matrixC depends solely on the mesh geometry and the matrix M on thegradient of the section, As a simple consequence we show thatA is singular only in the cases (i) element degeneracy () and (ii) collinearity of nodes (M not out of fullrank). We give constructions for the inversion of A in all cases. In one dimension, if A is non-singular, it has a simple explicitinverse. If A is singular we replace it by reduced matrix A^*.It can be shown that every case the spectral radius of the Jacobiiteration matrix ia ?and that A or A^* can be efficiently invertedby conjugate gradient methods. Finally, we discuss the applicability of these arguments tosystem of equations in any number of dimensions.     

Necessary conditions for infinite-horizon Volterra optimal control problems with time delay

Email: csfvega{at}dm.uba.ar Received on August 17, 2006; Accepted on September 8, 2007 Necessary conditions are proved for optimal control problemsinvolving an infinite horizon and terminal conditions at infinitywhose states are governed by Volterra integral equations withnon-linear time delay.     

Error estimates for the finite-element solution of an elliptic singularly perturbed problem

e-mail:ang{at}am.uni-erlangen.de The paper deals with an upwind discretization of finite-volume-typefor singularly perturbed elliptic boundary value problems. Stability,inverse monotonicity and convergence are considered. The mainobjective is to pursue the dependence of the error estimateon the perturbation parameter.     

Controllability and observability of impulsive hybrid dynamic systems

Corresponding author, Email: cbsoh{at}auto.eee.ntu.ac.sg This paper considers the controllability and observability ofhybrid dynamic systems with impulsive switching of the otherwisecontinuously evolving state. It derives necessary and sufficientconditions for impulsive hybrid dynamic systems to be controllableand observable.     

Structure of solution manifolds in a strongly coupled elliptic system

In this paper we study the steady-state solutions of a reaction-diffusionmodel, the Selkov scheme for glycolysis, under homogeneous Dirichletboundary conditions. Near to thermodynamic equilibrium, thestructure and stability of solutions are fully described. Abifurcation analysis is carried out, using the size of the regionin which the reaction takes place and one diffusion coefficientas main bifurcation parameters. The analysis helps us to understandthe nature of the bifurcation points, and determines the shapesand stability of the bifurcating manifolds in the neighbourhoodof the constant state. Local convergence of spectral methodsis shown, and some global pictures are calculated using path-followingtechniques. The framework we use can be applied to a wide varietyof reaction-diffusion systems. ^*Permanent address: Dpto. de Matemtica Aplicada, Fac. de CenciasQumicas, Universidad Complutense, 28040-Madrid, Spain. Current address: Department of Mathematics, Paisley Collegeof Technology, High Street, Paisley, Renfrewshire PA1 2BE, UK. Permanent address: Dpto. de Matemtica, Fac. de Matemtca, Astronomiay Fisica, Universidad National de Crdoba, 5000-Crdoba, R.Argentina.     

Estimates of the region of attraction of continuous-time cascade systems

Xiao-Song Yang Department of Mathematics, Huazhong University of Science and Technology Wuhan, Hubei 430074, People's Republic of China Huimin Li Departments of Mathematics and Control and Engineering, Huazhong University of Science and Technology Wuhan, Hubei 430074, People's Republic of China Corresponding author. Email: yangxs{at}cqupt.edu.cn Received on May 1, 2006; Accepted on October 6, 2006 In this paper, we first present some new sufficient conditionsfor global asymptotic stability of continuous-time cascade systems.Then, we give some detailed results on the region of attractionof continuous-time cascade systems. Some examples are presentedto illustrate these results.     

An analysis of the coupling of finite-element and Nystrm methods in acoustic scattering

In this paper we analyze the use of a combined Nystrm and finite-elementprocedure for approximating the scattered time-harmonic acousticfield produced when an incident wave interacts with a boundedinhomogeneity. The coupling technique uses a smooth artificialboundary and explicitly decouples the integral equation andfinite-element computations. We prove convergence of the method.One highlight of this analysis is that we prove Sobolev spaceerror estimates for the Nystrm scheme (which is usually analyzedin Hlder spaces). We also present some numerical results. E-mail: kirsch{at}am.uni-erlangen.de E-mail: monk{at}math.udel.edu     

The Convergence of Operator Approximations at Turning Points

Present address: Department of Mathematics, University of Reading, Reading RG6 2AX. We consider the convergence of solution curves of approximationsto parameter-dependent operator equations of the form G(, x)= 0. Provided G_x(, x) remains non-singular this problem is cateredfor by a simple extension to standard theory. In this paper,however, attention is concentrated on solution curves throughcertain singular points (⁰, x⁰), and the main result is thatconvergence depends on consistency and stability results forthe linear eigenvalue problem G_x(⁰, x⁰)0 = 0.     

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.     

Computation and parametrization of periodic and connecting orbits

Email: g.moore{at}ma.ic.ac.uk or na.gmoore New algorithms for the computation of periodic and connectingorbits of autonomous differential equations are developed. Previousmethods have parametrized these curves by the independent timevariable, but our procedure is based on the canonical arclengthparametrization.     

Time optimal control of retarded parabolic systems

Email: skrakowi{at}usk.pk.edu.pl Received on June 15, 2005; Accepted on July 12, 2006 In this paper, the time-optimal control problem for parabolicsystems in which time lags appear in the integral form bothin the state equation and in the Neumann boundary conditionis presented. The particular properties of the optimal controlare discussed.     

Therapy burden, drug resistance, and optimal treatment regimen for cancer chemotherapy

Email: boldrini{at}ime.unicamp.brAuthor to whom correspondence should be addressed. Email: michel{at}Incc.br Three nonlinear models of tumour cell growth under continuousdelivery of cycle nonspecific anticancer agents are studied.A dynamical optimization problem with the objective of minimizingthe final level of tumour cells is posed for these mathematicalsetups. The simplest setup does not possess toxicity constraints,whereas the other setups contain a dynamical equation describingthe therapy burden as a toxicity criterion. In addition, thethird setting contains the dynamics of drug resistant cells.A discussion conceming the optimal strategies of the respectivemodels is performed.     

Multiple criteria models for evaluation of competitive bids

Email: slliu{at}mail.east.net.cnAuthor for correspondence Email: mskklai{at}cityu.edu.hkEmail: sywang{at}iss02.iss.ac.cn In this paper, two general bidding systems are proposed to aidan owner to select one contractor based on multiple-attributedecision making models. They are named the general multiple-attributelower-bidder system and the general multiple-attribute average-bid system and are an extension (and refreshment) of the multiparameterbidding system and the average-bid method, respectively. Thepreference of the owner over the criteria (or attributes) relevantto construction work award is incorporated into the two newsystems if necessary and required by the owner. The values ofthe attributes (excluding the cost) required to be determinedby the owner in the multiparameter bidding system are not necessarilydetermined in the general multiple-attribute lower-bidder system.This reduces the burden to the owner. The two new bidding systemsare demonstrated via an example, and are compared with the multiparameterbidding system and the average-bid method. The comparisons showthat the preference of the owner in regard to the attributesrelevant to selection of contractors has a significant effecton bid evaluation. The two new systems preserve the form ofthe competitive bidding concept and can easily be applied topractice.