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

On the behaviour of blow-up interfaces for an inhomogeneous filtration equation

We study the asymptotic behaviour of blow-up interfaces of thesolutions to the one-dimensional nonlinear filtration equationin inhomogeneous media where m>1 isa constant and (x) = |x|^– (for |x| 1, with > 2) isa bounded, positive, smooth, and symmetric function. The initialdata are assumed to be smooth, bounded, compactly supported,symmetric, and monotone. It is known that due to the fast decayof the density (x) as |x| the support of the solution increasesunboundedly in a finite time T. We prove that as tT^– theinterface behaves like O((T – t)^–b), where the exponentb > 0 (which depends on m and only) is given by a uniqueself-similar solution of the second kind satisfying the equation|x|^– u_t = (u^m)_xx. The corresponding rescaled profilesalso converge. We establish the stability of the self-similarsolution of the second kind for the exponential density (x)=e^–|x|for |x| 1. We give a formal asymptotic analysis of the blow-upbehaviour for the non-self-similar density (x) = e^–|x|2.Several exact self-similar solutions and their correspondingasymptotics are constructed.     

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.     

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.     

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.     

Experimental Observation of the Large-Amplitude Solutions of Duffing's and Related Equations

Experiments with a nonlinear electronic model show that certainsimple features of the solutions of where f(u) is an odd monotonic function of u for example u³,repeat in a regular pattern as either is decreased or U isincreased. For fixed U, the position of these features is periodicin 1/ and, when f(u) has the form u|u|^k–1 a quantitativerelation between the period in 1/ and U can be found. The occurrenceof large-amplitude chaotic solutions is found to depend notonly on the nonlinearity of f(u) for large U but also on itsbehaviour near u = 0. For the Duffing equation, which can bereduced to the range of parameters accessible to experiment is 0<1 and0<F5000.     

The Range of Possible Values of f(x)

The purpose of this paper is to explain how to compute the rangeof possible values of a function of one variable, f(x), givenvalues of the function at n distinct points x₁ < x₂ <... < x_M–1 < x_M, and given a finite bound on thekth derivative of f: ||f^(k)|| L, 1 k n.     

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     

Injection of Ideal Fluid from a Slot into a Stream: Two Free Boundaries

A fluid is injected from a slot into a stream of another fluid.In a simple model this leads to a two-phase two-free-boundaryproblem with the jump relation |u^–|² – |u⁺|² = on the free boundary {u=0}, and |u| = ¹ on the free boundary{u > – Q}, where u is the stream function and Q isthe flux of the injected fluid. Using the variational theoryof Alt, Caffarelli & Friedman, we prove existence of (,₁, u) such that there is a smooth fit for both free boundaries.     

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.     

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.     

Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: the scalar case

^** Email: Paul.Houston{at}mcs.le.ac.uk^*** Email: Janice.Robson{at}comlab.ox.ac.uk^**** Email: Endre.Suli{at}comlab.ox.ac.uk We develop a one-parameter family of hp-version discontinuousGalerkin finite element methods, parameterised by [–1,1], for the numerical solution of quasilinear elliptic equationsin divergence form on a bounded open set ^d, d 2. In particular,we consider the analysis of the family for the equation –·{µ(x, |u|)u} = f(x) subject to mixed Dirichlet–Neumannboundary conditions on . It is assumed that µ is a real-valuedfunction, µ C( x [0, )), and thereexist positive constants m_µ and M_µ such that m_µ(t– s) µ(x, t)t – µ(x, s)s M_µ(t– s) for t s 0 and all x . Using a result from the theory of monotone operators for any valueof [–1, 1], the corresponding method is shown to havea unique solution u_DG in the finite element space. If u C¹() H^k(), k 2, then with discontinuous piecewise polynomials ofdegree p 1, the error between u and u_DG, measured in the brokenH¹()-norm, is (h^s–1/p^k–3/2), where 1 s min {p+ 1, k}.     

Correction: on the Intersection Matrices of Curves lying on Irregular Algebraic Surfaces

Professor W. F. Hammond has kindly drawn my attention to a blunderin 4 of the above paper. He referred to the ( – 2r) xß submatrix D of the skew-symmetric matrix displayednear the top of page 181, of which it is asserted that it issquare and non-singular, and pointed out that, from the factthat the matrix of which D forms part is regular, it may onlybe deduced that the columns of D are linearly independent; thatis, it only follows that – 2r ß. The validity of the equation – 2r = ß is essentialto the succeeding argument and, fortunately, may be establishedby alternative means. Using the nomenclature of the paper, wehave on F the set ₁^*, ..., _2r^*, ₁^*, ..., _ß^* of independent3-cycles (independent because they cut independent 1-cycleson the curve C), which may be completed, to form a basis forsuch cycles on F, by a further set ₁', ..., _2q–2r–pof independent 3-cycles, each of which meets C in a cycle homologousto zero on C. The cycles ₁^*, ..., ^* are invariant cycles andare independent on F so that, if > 2r + ß, thereis a non-trivial linear combination ^* of these having zero intersectionon C with each of the cycles ₁^*, ..., _2r^*, ₁^*, ..., _ß^*.Thus we have. (^* ._k^*)c = 0 = (^* ._i^*)c i.e. (^* ._k^*) = 0 = (^* ._i^* on F (1 k 2r; 1 i ß). Furthermore, (_j . C) 0 on C and we have (^* ._j .C)_C = 0 i.e. (^* ._j) = 0 on F (1 j 2q – 2r – ß). It now follows that ^* 0 on F (for it has zero intersectionwith every member of a basic set of 3-cycles on F). But thiscondradicts the assumption that ^* is a non-trivial linear combinationof the independent cycles ₁^*, ...,^*; and hence < 2r + ß.     

Asymptotic Solution of Hyperbolic Systems with Constant Coefficients

An integral representation of the exact solution of the initialvalue problem for the hyperbolic equation of the form is derived. Here A^o, A^v, B, and Care constant m x m matrices, u(t, X; ) is an m-component columnvector, and is a positive parameter. Various conditions areimposed on the coefficient matrices that permit the applicationof the method of stationary phase in several variables to theintegral representation of the exact solution. The leading termof the asymptotic expansion as of the exact solution is obtainedfor several types of initial data and source functions whichdepend on the parameter .     

Free convection boundary layers on a vertical surface in a heat-generating porous medium

The natural convection boundary-layer flow on a solid verticalsurface with heat generated within the boundary layer at a rateproportional to (T – T)^p (p 1) is considered. The surfaceis held at the ambient temperature T except near the leadingedge where it is held at a temperature above ambient. The behaviourof the flow as it develops from the leading edge is examinedand is seen to become independent of the initial heat input;however, it does depend strongly on the exponent p. For 1 p 2, the local heating eventually dominates at large distancesand there is a convective flow driven by this mechanism. Forp 4, the local heating does not have a significant effect,the fluid temperature remains relatively small throughout andthe heat transfer dies out through a wall jet flow. For 2 <p < 4, the local heating has a significant effect at relativelysmall distances, with a thermal runaway developing at a finitedistance along the surface.     

Some Integrals of Gamma, Polygamma and Volterra Functions

Infinite integrals of the gamma and polygamma functions areexpressed in terms of the Volterra and related functions usingthe Laplace transformation method. A number of infinite, convolution,indefinite and definite integrals of v(x), v(x, ), µ(x,ß) and µ(x, ß, ) has been evaluated.A short table of the Volterra function v(x) in the range 0 x 10–0 is presented.     

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 Representation of Some Integers as a Subset Sum

Let A N. The cardinality (the sum of the elements) of A willbe denoted by |A| ((A)). Let m N and p be a prime. Let A {1, 2,...,p}. We prove thefollowing results. If |A| [(p+m–2)/m]+m, then for every integer x such that0 x p – 1, there is B A such that |B| = m and (B) x mod p. Moreover, the bound is attained. If |A| [(p+m–2)/m]+m!, then there is B A such that |B| 0 mod m and (B) = (m – 1)!p. If |A| [(p + 1)/3]+29, then for every even integer x such that4p s x p(p + 170)/48, there is S A such that x = (S). In particular,for every even integer a 2 such that p 192a – 170, thereare an integer j 0 and S A such that (S) = a^j+1.     

Majorations Uniformes de Normes D'Inverses Dans Les Algebres de Beurling

The Beurling algebras l¹(D,)(D=N,Z) that are semi-simple, withcompact Gelfand transform, are considered. The paper gives anecessary and sufficient condition (on ) such that l¹(D,) possessesa uniform quantitative version of Wiener's theorem in the sensethat there exists a function :]0,+[]0,+ such that, for everyinvertible element x in the unit ball of l¹(D,), we have ||x^–1||(r(x^–1)) r(x^–1) is the spectral radiusof x^–1.     

On Lp boundedness of wave operators for 4-dimensional Schrodinger operators with threshold singularities

Let H=–+V(x) be a Schrödinger operator on L²(R⁴),H₀=–. Assume that |V(x)|+| V(x)|C x ^– for some>8. Let be the wave operators. It is known that W_± extend to bounded operators in L^p(R⁴)for all 1p, if 0 is neither an eigenvalue nor a resonance ofH. We show that if 0 is an eigenvalue, but not a resonance ofH, then the W_± are still bounded in L^p(R⁴) for all psuch that 4/3<p<4.