Conversion of Polynomials between Different Polynomial Bases

The purpose of this paper is to derive a recursive scheme forthe evaluation of the coefficients in the expansion , in terms of the coefficients in the expansion , where both q_k(x) and Q_k(x) are polynomials in xof degree k, and where both q_k(x) and Q_k{x} satisfy recursionformulae of the type satisfied by orthogonal polynomials. Thesets {Q_k(x)} and {q_k(x)} need not be orthogonal polynomials,though they usually are in the applications. An applicationis made to the evaluation of integrals with oscillatory andsingular integrands.     

Kazhdan-Lusztig Cells, q-Schur Algebras and James' Conjecture

We consider the Dipper–James q-Schur algebra S_q(n, r)_k,defined over a field k and with parameter q 0. An understandingof the representation theory of this algebra is of considerableinterest in the representation theory of finite groups of Lietype and quantum groups; see, for example, [6] and [11]. Itis known that S_q(n, r)_k is a semisimple algebra if q is nota root of unity. Much more interesting is the case when S_q(n,r)_k is not semisimple. Then we have a corresponding decompositionmatrix which records the multiplicities of the simple modulesin certain ‘standard modules’ (or ‘Weyl modules’).A major unsolved problem is the explicit determination of thesedecomposition matrices.     

Generating tuples of free products

Grushko's theorem [Mat. Sb. 8 (1940) 169–182] says thatany generating tuple (g₁, ..., g_m) of a free product H*K isNielsen-equivalent to a tuple (h₁, ..., h_l, k_l+1, ..., k_m) withh_i H and k_i K for all i. The h_i and k_i are clearly not unique.In this paper we address the extent of this non-uniqueness.     

The Fast Numerical Solution of Very Large Elliptic Difference Schemes

Let L_kv_k = g_k be a system of difference equations discretizingan elliptic boundary value problem. Assume the system to be"very large", that means that the number of unknowns exceedsthe capacity of storage. We present a method for solving theproblem with much less storage requirement. For two-dimensionalproblems the size of the needed storage decreases from O(h^–2)to (or even O(h^–5/4)). The computational work increasesonly by a factor about six. The technique can be generalizedto nonlinear problems. The algorithm is also useful for computerswith a small number of parallel processors.     

Strong Weighted Mean Summability and Kuttner's Theorem

In this paper we study sequence spaces that arise from the conceptof strong weighted mean summability. Let q = (q_n) be a sequenceof positive terms and set Q_n = ⁿ_k=1q_k. Then the weighted meanmatrix M_q = (a_nk) is defined by if kn, a_nk=0 if k>n. It is well known that M_q defines a regular summability methodif and only if Q_n. Passing to strong summability, we let 0<p<.Then , are the spaces of all sequences that are strongly M_q-summablewith index p to 0, strongly M_q-summable with index p and stronglyM_q-bounded with index p, respectively. The most important specialcase is obtained by taking M_q = C₁, the Cesàro matrix,which leads to the familiar sequence spaces w₀(p), w(p) and w(p), respectively, see [4, 21]. We remark that strong summabilitywas first studied by Hardy and Littlewood [8] in 1913 when theyapplied strong Cesàro summability of index 1 and 2 toFourier series; orthogonal series have remained the main areaof application for strong summability. See [32, 6] for furtherreferences. When we abstract from the needs of summability theory certainfeatures of the above sequence spaces become irrelevant; forinstance, the q_k simply constitute a diagonal transform. Hence,from a sequence space theoretic point of view we are led tostudy the spaces     

A new transform method for evolution partial differential equations

We introduce a new transform method for solving initial-boundary-valueproblems for linear evolution partial differential equationswith spatial derivatives of arbitrary order. This method isillustrated by solving several such problems on the half-line{t > 0, 0 < x < }, and on the quarter-plane {t >0, 0 < x_j < , j = 1, 2}. For equations in one space dimensionthis method constructs q(x, t) as an integral in the complexk-plane involving an x-transform of the initial condition anda t-transform of the boundary conditions. For equations in twospace dimensions it constructs q(x₁, x₂, t) as an integral inthe complex (k₁, k₂)-planes involving an (x₁, x₂)-transformof the initial condition, an (x₂, t)-transform of the boundaryconditions at x₁ = 0, and an (x₁, t)-transform of the boundaryconditions at x₂ = 0. This method is simple to implement andyet it yields integral representations which are particularlyconvenient for computing the long time asymptotics of the solution.     

Implicit finite difference schemes for a linear model of well-reservoir coupling

Standard reservoir models usually consider wells as Dirac measuresover an interval length. Moreover, the well-reservoir couplingis taken into account under quite simplified assumptions. Mostrecently, however, attention has been drawn to the fact thatin some situations, such as those related to non-vertical wells,these simplifications do not allow us to model some relevantmechanisms of the coupled flow. Therefore, more complex alternativemodels have been proposed recently in the oil reservoir simulationliterature. A linearized version for the well-reservoir couplingcan be written, in an appropriate functional setting, in theform U'(t) + AU(t) = F(t). In this work we discuss implicitin time discretizations of this equation, of the form { Uⁿ⁺¹_h - Uⁿ_h÷+A_hUⁿ⁺¹_h = Fⁿ⁼¹_h, U⁰_h=U_0.h We propose two different approximations, corresponding to first-and second-order spatial truncation errors, and we establishthe convergence of both approximations.     

A multigrid-type method for thin plate spline interpolation on a circle

We consider the problem of thin plate spline interpolation ton equally spaced points on a circle, where the number of datapoints is sufficiently large for work of O(n³ to be unacceptable.We develop an iterative multigrid-type method, each iterationcomprising ngrid stages, and n being an integer multiple of2^ngrid–1. We let the first grid, V₁ be the full set ofdata points, V say, and each subsequent (coarser) grid, V_k,k=2, 3,...,ngrid, contain exactly half of the data points ofthe preceding (finer) grid, these data points being equallyspaced. At each stage of the iteration, we correct our current approximationto the thin plate spline interpolant by an estimate of the interpolantto the current residuals on V_k, where the correction is constructedfrom Lagrange functions of interpolation on small local subsetsof p data points in V_k. When the coarsest grid is reached, however,then the interpolation problem is solved exactly on its q=n/2^ngrid–1points. The iterative process continues until the maximum residualdoes not exceed a specified tolerance. Each iteration has the effect of premultiplying the vector ofresiduals by an n x n matrix R, and thus convergence will dependupon the spectral radius, (R), of this matrix. We investigatethe dependence of the spectral radius on the values of n, p,and q. In all the cases we have considered, we find (R) <<1, and thus rapid convergence is assured.     

SCALED ASYMPTOTICS FOR SOME q-SERIES

This work, investigates the asymptotics for Euler’s q-exponentialE_q(z), Ramanujan’s function A_q(z), Jackson’s q-Besselfunction J_v⁽²⁾ (z; q), the Stieltjes–Wigert orthogonalpolynomials S_n(x; q) and q-Laguerre polynomials L_n⁽⁾ (x; q)as q approaches 1.     

Mixed hp-DGFEM for incompressible flows II: Geometric edge meshes

We consider the Stokes problem of incompressible fluid flowin three-dimensional polyhedral domains discretized on hexahedralmeshes with hp-discontinuous Galerkin finite elements of typeQ_k for the velocity and Q_k–1 for the pressure. We provethat these elements are inf-sup stable on geometric edge meshesthat are refined anisotropically and non-quasiuniformly towardsedges and corners. The discrete inf-sup constant is shown tobe independent of the aspect ratio of the anisotropic elementsand is of O(k^–3/2) in the polynomial degree k, as in thecase of conforming Q_k–Q_k–2 approximations on thesame meshes.     

Buckling of an axisymmetric vesicle under compression: the effects of resistance to shear

We consider the axisymmetric deformation of an initially spherical,porous vesicle with incompressible membrane having finite resistanceto in-plane shearing, as the vesicle is compressed between parallelplates. We adopt a thin-shell balance-of-forces formulationin which the mechanical properties of the membrane are describedby a single dimensionless parameter, C, which is the ratio ofthe membrane's resistance to shearing to its resistance to bending.This results in a novel free-boundary problem which we solvenumerically to obtain vesicle shapes as a function of plateseparation, h. For small deformations, the vesicle contactseach plate over a small circular area. At a critical value ofplate separation, h_TC, there is a transcritical bifurcationfrom which a new branch of solutions emerges, representing buckledvesicles which contact each plate along a circular curve. Forthe values of C investigated, we find that the transcriticalbifurcation is subcritical and that there is a further saddle-nodebifurcation (fold) along the branch of buckled solutions ath = h_SN (where h_SN > h_TC). The resulting bifurcation structureis commensurate with a hysteresis loop in which a sudden transitionfrom an unbuckled solution to a buckled one occurs as h is decreasedthrough h_TC and a further sudden transition, this time froma buckled solution to an unbuckled one, occurs as h is increasedthrough h_SN. We find that h_SN and h_TC increase with C, thatis, vesicles that resist shear are more prone to buckling.     

Multiple Blocking Sets and Arcs in Finite Planes

This paper contains two main results relating to the size ofa multiple blocking set in PG(2, q). The first gives a verygeneral lower bound, the second a much better lower bound forprime planes. The latter is used to consider maximum sizes of(k, n)-arcs in PG(2, 11) and PG(2, 13), some of which are determined.In addition, a summary is given of the value of m_n(2, q) forq 13.     

Strang-type preconditioners for systems of LMF-based ODE codes

We consider the solution of ordinary differential equations(ODEs) using boundary value methods. These methods require thesolution of one or more unsymmetric, large and sparse linearsystems. The GMRES method with the Strang-type block-circulantpreconditioner is proposed for solving these linear systems.We show that if an A_k1,k2 -stable boundary value method is usedfor an m-by-m system of ODEs, then our preconditioners are invertibleand all the eigenvalues of the preconditioned systems are 1except for at most 2m(k₁ + k₂) outliers. It follows that whenthe GMRES method is applied to solving the preconditioned systems,the method will converge in at most 2m(k₁ + k₂) + 1 iterations.Numerical results are given to illustrate the effectivenessof our methods. Received 8 October 1999. Accepted 30 May 2000.     

Homogenization of the Stefan Problem and Application to Magnetic Composite Media

The theory of homogenization (Bensoussan, Lions & Papanicolaou,1978) shows that u, the solution of the diffusion equation [with k(y) periodic in the space-variable y and q = cu a linearfunction of u] has a weak limit u for = 0. This theory allowsone to compute, for a given k, the conductivity tensor of ananisotropic but homogeneous medium in which, for unchanged initialand boundary conditions, u is the solution of the diffusionequation. We examine here the case where the relation between q and uis given by a maximal monotone graph (i.e. the Stefan problem),depending on the space variable in the same manner as k. Applicationsto eddy-current problems in magnetic composite media (steelcables, laminations) are suggested. A numerical example is given.     

Well matched systems

The product (3.10) on page 33 is incorrectly called a cartesianproduct on pages 33 and 35. This misnomer in effect amountsto a wrong definition. The product (3.10) should be definedso that the right-hand member of (3.10) is the set of all sums_f=1 f_j (not the set of all ordered q-tuples) such that f₁ F(m₁,d₁), ..., f_q F(m_q, d_q).     

Stability in the Sense of Bounded Average Power

An R^m-valued sequence (x_k): = (x_k : k = 1, 2, ...), e.g. generatedrecursively by x_k = f_k (x_k–_k, U_k), is called ‘averagepth power bounded’ if (1/K) is bounded uniformly in K= 1, 2,.... (The case p = 2 may correspond to ‘power’in the physical sense.) This is a notion of stability. Givenestimates of the form: f_k (x, u) < a x + ¶ _k conditionsare obtained on the coefficient sequence (a_k) and the inputestimates e_k:=¶_k (u_k) which ensure this form of stabilityfor the output (x_k). In particular, a condition (utilized inan application to adaptive control) is obtained which imposes(i) a bound b on (a_k) and a ‘sparsity measure’ m(K) on #{kK: a_k>} as K ( >1) (ii) average pth power boundednesson (e_k), and (iii) a growth condition on (e_k) related to b andm (•). This condition is sharp.     

De Rham Cohomology and Hodge Decomposition For Quantum Groups

Let be one of the N²-dimensionalbicovariant first order differential calculi for the quantumgroups GL_q(N), SL_q(N), SO_q(N), or Sp_q(N), where q is a transcendentalcomplex number and z is a regular parameter. It is shown thatthe de Rham cohomology of Woronowicz' external algebra coincides with the de Rham cohomologiesof its left-coinvariant, its right-coinvariant and its (two-sided)coinvariant subcomplexes. In the cases GL_q(N) and SL_q(N) thecohomology ring is isomorphic to the coinvariant external algebra and to the vector space of harmonic forms. We prove a Hodge decomposition theorem in thesecases. The main technical tool is the spectral decompositionof the quantum Laplace-Beltrami operator. 2000 MathematicalSubject Classification: 46L87, 58A12, 81R50.     

An explicit solution for an instantaneous two-phase Stefan problem with nonlinear thermal coefficients

We consider a nonlinear heat conduction problem for a semi-infinitematerial x > 0, with phase-change temperature T₁, an initialtemperature T₂ (> T₁) and a heat flux of the type q (t) =q₀/t imposed on the fixed face x = 0. We assume that the volumetricheat capacity and the thermal conductivity are particular nonlinearfunctions of the temperature in both solid and liquid phases. We determine necessary and/or sufficient conditions on the parametersof the problem in order to obtain the existence of an explicitsolution for an instantaneous nonlinear twophase Stefan problem(solidification process).     

Asymptotics for the Heat Content of a Planar Region with a Fractal Polygonal Boundary

Let k 3 be an integer. For 0<s<1, let D_s R² be the setthat is constructed iteratively as follows. Take a regular openk-gon with sides of unit length, attach regular open k-gonswith sides of length s to the middles of the edges, and so on.At each stage of the iteration the k-gons that are added area factor s smaller than the previous generation and are attachedto the outer edges of the family grown so far. The set D_s isdefined to be the interior of the closure of the union of allthe k-gons. It is easy to see that there must exist some s_k> 0 such that no k-gons overlap if and only if 0 < s s_k. We derive an explicit formula for s_k. The set D_s is open, bounded, connected and has a fractal polygonalboundary. Let denote the heat content of D_s at time t when D_s initially has temperature 0and D_s is kept at temperature 1. We derive the complete short-timeexpansion of up to terms that are exponentially small in 1/t. It turns out that there arethree regimes, corresponding to 0<s<1/(k–1), s=1/(k–1),and 1/(k–1)<s s_k. For s 1/(k–1) the expansionhas the form where p_s is a log (1/s²)-periodic function, d_s=log (k–1)/log(1/s) is a similarity dimension, A_s and B are constants relatedto the edges and vertices, respectively, of D_s, and r_s is anerror exponent. For s=1/(k–1), the t^1/2-term carries anadditional log t. 1991 Mathematics Subject Classification: 11D25,11G05, 14G05.     

Convergent and spurious solutions of nonlinear elliptic equations

In this paper we investigate finite element approximations ofnonlinear elliptic equations in three dimensions. By applyingand extending the results of Lopez-Marcos and Sanz-Serna, weprove that the finite element approximation on a mesh of sizeh, has a solution U_k which converges to an exact solution ofthe differential equation as h0. This solution is unique withina suitably defined stability ball Bh. For the particular nonlinearequation u + (u + u^p) we show that the size of B_h depends uponh only if p > 5 when it tends to zero as h 0. In this casewe prove the existence of spurious solutions V_h of the Galerkinapproximation which become unbounded in the maximum norm ash0. The stability ball B_h then acts to separate the convergentand the spurious solutions. We present the results of some numericalexperiments to substantiate our claims.