Numerical Evaluation of Fourier Integrals

A method is developed for evaluating Fourier integrals of theform A() = ¹_–1f(x) e^fax dx, 0. The method consists of expanding the function f in a seriesof Chebyshev polynomials and expressing the integral A() asa series of the Bessel functionsJ_r+(), r= 0, 1, 2,.... A partialsum A_N() of the series provides an approximant to A(). The principalfeature of the method is that one set of N+1 evaluations off(x) suffices for the calculation of A_N() for all , and alsothe truncation error A()–A_N() is essentially independentof . Numerical tests show that the method is accurate, economicaland reliable. An application to the inversion of Fourier andLaplace transforms is briefly described.     

Computer extended series solution for unsteady flow in a contracting or expanding pipe

Unsteady flow in a semi-infinite contracting or expanding pipeis reinvestigated using long series analysis. The proposed seriesmethod is useful in analysing the problem for a moderately largeconstant ( = aa/, where a = a(t), the radius of the pipe isa function of time, a(t) is the velocity of the wall, and iskinematic viscosity). For positive values of (expansion ofthe pipe) accuracy of the series representing shear stress andpressure gradient is increased from = 2.89 to = 6.0 by extractingthe singularity followed by completion of the series. For negativevalues of (contraction of the pipe), we revert the series whichresults into the increase of the region of validity of the transposedseries from = -25.0 to = -2.89. Later we use Padé approximantsfor summing them. Also, the asymptotic solution for large valuesof is obtained and it agrees closely with pure numerical valuesof shear stress at the wall and pressure gradient.     

Elliptic Problems Involving Phase Boundaries Satisfying a Curvature Condition

Two theorems related to equilibrium free-boundary problems arepresented. One arises as a time-independent solution to thephase-field equations. The other is the relevant time-independentproblem for the Stefan model, modified for the surface tensioneffect. It also serves as a preliminary result for the phase-fieldformulation. Under appropriate conditions, we prove that, givenan appropriate positive constant and a smooth function u: R;,where is an annular domain in R², there exists a curve suchthat u(x)=—K(x) for all x , where K is the curvature.Using this result, we prove the existence of solutions to O=²+ ?(—³) + 2u that have a transition layer behaviour (from=—1 to =+1) for small and make the transition on thecurve . This proves there exist solutions to the phase fieldmodel that satisfy a Gibbs-Thompson relation.     

On hearing the shape of a bounded domain with Robin boundary conditions

The asymptotic expansions of the trace of the heat kernel (t)= [sum ]_j=1 exp(-t_j) for small positive t, where {_j} _j=1 arethe eigenvalues of the negative Laplacian -_n = -[sum ]ⁿ_k=1 (/x^k)²in Rⁿ (n = 2 or 3), are studied for a general multiply connectedbounded domain which is surrounded by simply connected boundeddomains _i with smooth boundaries _i (i = 1,...,m), where smoothfunctions Y_i (i = 1,...,m) are assuming the Robin boundary conditions(n_i + Y_i) = 0 on _i. Here /n_i denote differentiations along theinward-pointing normals to _i (i = 1,...,m). Some applicationsof an ideal gas enclosed in the multiply connected bounded containerwith Neumann or Robin boundary conditions are given.     

The Model Plasma Problem: Existence of a Solution Branch

We are interested in the model plasma problem –u = u⁺in ,u = –d on , _au⁺ dx=j where is a bounded domain in with boundary ; here, j isa given positive number, the function u and the positive number are the unknowns of the problem, and d is a real parameter.Using a variant of the implicit function theorem, we can provethe existence of a global solution branch parametrized by d.The method has the advantage that it can be used for analysingthe approximation of the above problem by a finite-element method.     

The large deformation of nonlinearly elastic rings in a two-dimensional compressible flow

The nonlinear nonlocal system of the equilibrium equations ofan elastic ring under the action of an external two-dimensionaluniformly subsonic potential barotropic steady-state gas flowis considered. The configurations of the elastic ring are identifiedby a pair of functions (, ). The simple curve represents theshape of the ring and the real-valued function identifies theorientation of the material sections of the ring. The pressurefield on the ring depends nonlocally on , and on two parametersU and P which represent the pressure and the velocity at infinity.The system is shown to be equivalent to a fixed-point problem,which is then treated with continuation methods. It is shownthat the solution branch ensuing from certain equilibrium states((₀, ₀), 0, P₀) in the solution-parameter space of ((₀, ₀),0, P₀) either approaches the boundary of the admissible ((,), U,p)'s in a well-defined sense, or is unbounded, or is homotopicallynontrivial in the sense that there exists a continuous map from the branch to a two-dimensional sphere which is not homotopicin the sphere to a constant, while restricted to the branchminus ((₀, ₀), 0, P₀) is homotopic to a constant in the sphere.Furthermore, by fixing the pressure parameter at P₀ and by consideringthe one-parameter problem in ((, ), U), the following holds.Every hyperplane in the solution-parameter space of the ((,), U)'s which contains the equilibrium state ((₀, ₀), 0) anddoes not include a welldetermined one-dimensional subspace intersectsthe solution branch above at a point different from ((₀, ₀),0).     

The Magnetogasdynamic Boundary Layer on a Thermally Insulated Flat Plate

This paper examines the effect of compressibility on the flowin the boundary layer on a semi-infinite, thermally insulatedflat plate placed at zero incidence to a uniform stream of electricallyconducting gas, with an aligned magnetic field at large distancesfrom the plate. The present discussion is limited to small values of the conductivityparameter = 4µv, and the Prandtl number is taken to beunity. The latter assumption permits a simplification of theanalysis, and the former allows the dependence of the flow onthe parameters ß = µH²/4U² and M = U/cto beadequately illustrated without excessive computation. A seriessolution valid for small values of the conductivity parameterand for Mach numbers not too large is derived. Values of ß = 0.3 and 0.5, = 0.01 and 0.1 are consideredand for those values the skin friction decreases with increasingMach number, similar to the case when ß = 0. The analysissuggests that for larger values of ß the skin frictionmight even increase with the Mach number initially. This iscertainly the case with the tangential component of the magneticfield, which for ß = 0.5 exhibits a maximum at approximatelyunit Mach number. The reason for this behaviour lies in thefact that, in view of the temperature changes taking place inthe flow, the electrical conductivity and thereby the localvalue of can change by more than an order of magnitude. Thishas the effect of giving results which are akin to those forarbitrary large in incompressible flow even though the valueof based on the main stream gas properties remains low.     

An Asymptotic Analysis of the Delayed Logistic Equation when the Delay is Large

We show how to construct an asymptotic solution to the delayedlogistic equation = y(1–y₁), corresponding to the asymptoticlimit . The results of the analysis are compared with a numericalcomputation, and found to be comparatively accurate for >2. Since the approach adopted is novel, we comment on some featureswhich may be relevant in other problems.     

An analogue of Berry's conjecture for the phase in fractal obstacle scattering

In this paper, the authors consider the high-frequency asymptoticsof the phase s() of acoustic waves scattered by an obstacleRⁿ with fractal boundary. Under certain conditions, it is provedthat if is –Minkowski measurable with –Minkowskimeasure µ then there exists a positive constant C_n, dependingonlyon n and such that where     

A Criticality Criterion for Thermal Explosions in Well-stirred Systems with Reactant Consumption

The autonomous differential equations for the temperature andreactant consumption in a first-order well-stirred exothermicreaction are considered. An examination of the phase-plane solutionsallows the qualitative behaviour of the Semenov number as afunction of maximum temperature rise ^* to be established. Inthe limit of infinite adiabatic temperature rise (B) and zeroactivation energy parameter ( = 0), the relationship between and stationary temperature _s is known to be e_¹ = _s. Criticalityarises at the maximum of (_s) and leads to the critical Semenovvalues (_s)_cr = 1, _cr = e^–1. For sufficiently large B,it is shown that the (^*) curve has a bifurcation at ^* = 1, withthe upper branch monotonically increasing and the lower branchmonotonically decreasing for ^* > 1. In the limit B thesebecome respectively the straight line = e^–1, _s 1 andthe unstable branch of = _se^–1, _s 1 and the unstablebranch of = _s e^–. Criticality for finite B is definedas occurring at the bifurcation, namely ^*_cr = 1, with _cr(B)the value of at this point. Values of these Semenoy numbersare obtainable from the numerical calculations of Boddingtonet al. [Proc. R. Soc. Lond. (1983), 390, 13–30]. The newcriterion is applied to an approximate phase-plane solution.The corresponding critical parameter is found to be _cr = e^–1[1+B^–(2–e^–1)+O(B^–1)].     

Geometry of Critical Loci

Let :(Z,z)(U,0) be the germ of a finite (that is, proper with finite fibres)complex analytic morphism from a complex analytic normal surfaceonto an open neighbourhood U of the origin 0 in the complexplane C². Let u and v be coordinates of C² defined on U. Weshall call the triple (, u, v) the initial data. Let stand for the discriminant locus of the germ , that is,the image by of the critical locus of . Let ()_A be the branches of the discriminant locus at O whichare not the coordinate axes. For each A, we define a rational number d by where I(–, –) denotes the intersection number at0 of complex analytic curves in C². The set of rational numbersd, for A, is a finite subset D of the set of rational numbersQ. We shall call D the set of discriminantal ratios of the initialdata (, u, v). The interesting situation is when one of thetwo coordinates (u, v) is tangent to some branch of , otherwiseD = {1}. The definition of D depends not only on the choiceof the two coordinates, but also on their ordering. In this paper we prove that the set D is a topological invariantof the initial data (, u, v) (in a sense that we shall definebelow) and we give several ways to compute it. These resultsare first steps in the understanding of the geometry of thediscriminant locus. We shall also see the relation with thegeometry of the critical locus.     

Herman Rings and Arnold Disks

For (,a) C* x C, let f_,a be the rational map defined by f_,a(z)= z² (az+1)/(z+a). If R/Z is a Brjuno number, we let D bethe set of parameters (,a) such that f_,a has a fixed Hermanring with rotation number (we consider that (e²ⁱ,0) D). Resultsobtained by McMullen and Sullivan imply that, for any g D, theconnected component of D(C* x (C/{0,1})) that contains g isisomorphic to a punctured disk. We show that there is a holomorphic injection F:DD such thatF(0) = (e²ⁱ ,0) and , where r is the conformal radius at 0 of the Siegel disk of the quadraticpolynomial z e²ⁱ z(1+z). As a consequence, we show that for a (0,1/3), if f_l,a has afixed Herman ring with rotation number and if m_a is the modulusof the Herman ring, then, as a0, we have e m_a=(r/a) + O(a). We finally explain how to adapt the results to the complex standardfamily z e^(a/2)(z-1/z).     

Sub- and Superadditive Properties of Fejer's sine Polynomial

Let be Fejér's sine polynomial. We prove the following statements.

1. The inequality holds for all x, y (0, ) with x + y < if and only if 0 and + rß 1.
2. The converse of the above inequality is valid for allx, y (0, ) with x + y < if and only if 0 and + rß 1.
3. For all n N and x, y [0, ] we have . Both bounds are best possible.

2000 Mathematics Subject Classification 42A05, 26D05 (primary),39B62 (secondary).     

The Critical Conditions in a Reactive Slab with Slowly Varying Surface Temperature

Liouville's non-linear partial differential equation is consideredfor an infinite rectangular strip domain with a slowly varyingboundary condition. The equation describes a layer of chemicallyreactive material under conditions where the resistance to surfaceheat transfer is negligible and the ambient temperature variesslowly along the surface. Symmetrical heating by a zero orderexothermic reaction is assumed. If is a small dimensionlesstemperature difference between regions where the surface temperatureis effectively constant, a perturbation series solution in may be determined provided the Frank-Kamenetskii parameter satisfies _c(). It is shown that a plausible value for thecritical parameter is _c() = _c(0) e^–e,where _c(0) = 0.878.The corresponding critical temperature distribution is shownto have a dependence on different from that for subcriticalcases.     

On the Slow Rotation of a Sphere about a Diameter Parallel to a Nearby Plane Wall

A method using a matched asymptotic expansions technique ispresented for obtaining the Stokes flow solution for a rigidspherical body of radius a rotating uniformly about a diameterparallel to a fixed plane wall when the minimum clearance ais very much smaller than a. An inner solution is constructedwhich is valid for the region in the neighbourhood of the nearestpoints of the sphere and the wall where the flow is stronglysheared with large velocity gradients and pressure; in thisregion the leading term of the asymptotic expansion of the solutionsatisfies the equations of lubrication theory. A matching outersolution is constructed which is valid in the remainder of thefluid where the flow is weakly sheared and it is possible toassume = 0. The forces and couples acting on the sphere andthe wall are shown to be of the form (₀+₁) log +ß₀+0(,where ₀, ₁ and ß₀ are constants which have been determinedexplicitly. By use of these results it is shown that the problemwhen the sphere rolls on the wall is not well posed.     

On a Penalty-Perturbation Theory for Plate Problems

A penalty-perturbation method previously proposed by Westbrook(J. Inst. Maths Applics (1974) 14, 79–82) for the solutionsof static bending problems for elastic plates is analysed here.The method replaces the single fourth-order biharmonic equationby a system of three second-order equations which is "singularly"perturbed with respect to a small penalty parameter . The existenceof solutions of the perturbed problem for each > 0 is establishedand the behaviour of these solutions as 0 0 is studied. Inparticular, the results show that while these solutions arecontinuous in at = 0, analyticity in at = 0 is lost exceptin special cases.     

Thermodynamic Inequalities in Continuum Mechanics

The constitutive relations for the transport of heat, stress,electric charge, etc., in a continuum must be chosen so thatthe second law of thermodynamics is not violated; the constraintstake the form of inequalities, typically requiring the entropygenerated within a material element to be non-negative. Thepaper is concerned with this concept—its history, thephysical principles on which it depends, how to apply it whensecond-order or non-linear effects are important and how itis widely misused in modern continuum mechanics. The history is reduced to the contributions of five leadingthermodynamicists—Clausius, Maxwell, Gibbs, Boltzmannand Duhem. The object here was to try to discover which formof the inequality one should regard as being fundamental. Oneimportant conclusion is that entropy S must be defined simultaneouslywith the identification of the inequality, and that in generalthis cannot be done until the constitutive equations are known.The empirical element enters with the notion of irreversibility,which is given a precise meaning with the aid of the motionreversed parity (x), a variable x having = +1 or = –1if, when time and motions are reversed, x x or x –x.The macroscopic parity of x, _*(x), is obtained by first replacingx by the constitutive equation for x. The entropy production rate has both irreversible (f) and reversible(r) parts. It is shown that the reciprocal relations followfrom the requirement that the macroscopic parity of (i) mustbe +1. Continuum thermodynamics is based on various principles extractedfrom theory developed for uniform systems, the example chosento illustrate the ideas being the simple monatomic gas. Second-orderconstitutive relations are introduced, and the expressions forentropy and its production rate per unit volume, , obtained.It is shown that the stability condition 0 cannot, in general,be satisfied merely by imposing constraints on the constitutiverelations. To second-order = ₁ + ₂, where ₁ is the usual bilinearform, and the terms in ₂ have an additional derivative. Thesecond-order term ₂ can have both signs, and is not dissipative.The relation between this fact and the frame-dependence of constitutiverelations is explained. The final section illustrates the errors frequently found inthe thermodynamic arguments appearing in books and papers onrational continuum mechanics. The principle of these is that 0 is interpreted as being a constraint on the constitutiverelations alone. Another is the idea that the balance equationscan be set aside as constraints by regarding them as mere definitionsof a heat source and a body force, an error based partly onthe misconception that constitutive relations should be frame-indifferent.Finally, an inequality due to Glansdorff & Prigogine isexamined and found to be in error.     

The Isolated Points of G{rho} and the w*-Strongly Exposed Points of P{rho}(G)0

Let G be a separable locally compact group and let be its dualspace with Fell's topology. It is well known that the set P(G)of continuous positive-definite functions on G can be identifiedwith the set of positive linear functionals on the group C^*-algebraC^*(G). We show that if is discrete in , then there exists anonzero positive-definite function associated with such that is a w^*-strongly exposed point of P(G)₀, where P(G)₀={f P(G):f(e)1. Conversely, if some nonzero positive-definite function associatedwith is a w^*-strongly exposed point of P(G)₀, then is isolatedin . Consequently, G is compact if and only if, for every ,there exists a nonzero positive-definite function associatedwith that is a w^*-strongly exposed point of P(G)₀. If, in addition,G is unimodular and , then is isolated in if and only if somenonzero positive-definite function associated with is a w^*-stronglyexposed point of P(G)₀, where is the left regular representationof G and is the reduced dual space of G. We prove that if B(G)has the Radon–Nikodym property, then the set of isolatedpoints of (so square-integrable if G is unimodular) is densein . It is also proved that if G is a separable SIN-group, thenG is amenable if and only if there exists a closed point in. In particular, for a countable discrete non-amenable groupG (for example the free group F₂ on two generators), there isno closed point in its reduced dual space .     

Application of JMN Approximants to Numerical Initial-value Problems in Linear Differential-algebraic Systems

For f: [0,) R, the J_MN approximant of f(t) is where M and N are non-negative integers and ß_k, _i,K_iare defined constants. Under appropriate conditions on f andprovided Re(_i) > 0 The approximants are the bases of recursions for numerical initial-valueproblems in linear differential-algebraic systems with constantcoefficients. The recursions are stable when N M N–2.Each step of a recursion involves mainly the solution of N/2uncoupled algebraic systems.     

Calculs De Facteurs Epsilon De Paires Pour GLn Sur Un Corps Local, I

Soient F un corps commutatif localement compact non archimédienet un caractère additif non trivial de F. Soient unereprésentation du groupe de Weil–Deligne de F,et sa contragrédiente. Nous calculons le facteur (, , ). De manière analogue, nous calculons le facteur (x, , ) pour toute représentationadmissible irréductible de GL_n(F). En conséquence,si F est de caractéristique nulle et si et se correspondentpar la correspondance de Langlands construite par M. Harris,ou celle construite par les auteurs, alors les facteurs (, , s) et (x, , s) sont égaux pour tout nombre complexe s. Let F be a non-Archimedean local field and a non-trivial additivecharacter of F. Let be a representation of the Weil–Delignegroup of F and its contragredient representation. We compute (, , ). Analogously, we compute (x, , ) for all irreducible admissible representations of GL_n(F).Consequently, if F has characteristic zero, and , correspondvia the Langlands correspondence established by M. Harris orthe correspondence constructed by the authors, then we have(, , s) = (x, , s) for all sC. 1991 Mathematics Subject Classification22E50.