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.     

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     

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.     

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

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

Exponential stability and maximal attractors for a one-dimensional nonlinear thermoviscoelasticity

This paper is concerned with the global existence, exponentialstability of solutions and associated nonlinear C₀-semigroupas well as the existence of maximal attractors in Hⁱ (i = 1,2, 4) for a nonlinear one-dimensional thermoviscoelasticitydescribing a kind of solid-like material. Some new ideas andmore delicated estimates are employed to prove the global existenceand exponential stability of solutions. The important featurefor the existence of maximal attractors in Hⁱ₊ (i = 1, 2, 4)is that the metric spaces H¹₊, H²₊ and H⁴₊ we work with arethree incomplete metric spaces, as can be seen from the physicalconstraints, i.e. > 0 and u > 0, with and u being absolutetemperature and deformation gradient (strain). For any positiveparameters ₁, ₂, ..., ₅ verifying some conditions, a sequenceof closed subspaces Hⁱ Hⁱ+ (i = 1, 2, 4) is found, and theexistence of maximal attractors in Hⁱ (i = 1, 2, 4) is established.     

Great time delay in a system with material cycling and delayed biomass growth

In this paper, we consider the effects of time delay on thelinear stability of the positive equilibrium of a reaction–diffusionmodel with time delayed nutrient recycling. We take the twodelay kernels as f(s) = (s), g(s) = p(0) + (1 – p)(r),s (0, +), r [0, +), 0 < p < 1, where is the Dirac deltafunction. When both N₁ and N₂ are space-independent, we canobserve that the increasing time delay destabilizes the positiveequilibrium, and stability-switches from stability to instabilityto stability occur. That is, for some special time delay kernels,the equilibrium is stable for small time delay, bifurcates towardsstable oscillations when the time delay is increased, then regainsstability for larger values. When both N₁ and N₂ are space-dependent,we are concerned with the effects of the delay kernels as wellas diffusion on the dynamics of the system.     

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.     

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 Structure of the Solution Set for Periodic Oscillations in a Suspension Bridge Model

Numerical results are reported for the computation of periodicsolution paths for a suspension bridge model represented bythe equation. u_n + EIu_xxxx + u_t, + Ku⁺ = W(^x) + sin x sin µt. with hinged-end boundary conditions, as the forcing amplitude and frequency µ are varied. The term Ku⁺ models the factthat there is restoring force due to the cables only when theyare being stretched. It is found that an S-shaped curve is obtainedwhen the displacement amplitudes are plotted against the forcingamplitudes for some frequency regimes. As a two-parameter problem,it appears that the solution set resembles a cusp-like surfacewith the singular point near linear resonance. While the effectof strengthening the cable (i.e. increasing K) will enhancethe occurrence of the multiple solutions, the effect of increasingthe damping coefficient gives the opposite effect.     

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.     

The Decay of Eddy-currents in Thin Sheets and Related Water-wave Problems

The decay of the eddy-currents that are induced in a thin, uniform,imperfectly-conducting sheet by switching off the source ofan external magnetic field is investigated. For the two-dimensionalproblem of an infinite strip the (non-dimensional) decay constants_n and eddy-current distributions i_n(x) are the eigenvalues andeigenfunctions of the integral equation with the constraint. For the circular disc the corresponding equation is where and K and E are complete elliptic integrals. For both problemsthe initial eddy-currents have inverse-square-root singularitiesat the edges but during their decay the eddy currents are finiteat the edges and the normal magnetic fields have logarithmicsingularities there. Numerical results are given for variousinitial-value problems. The eddy current problems are closely related to water-waveproblems in which there is a strip-shaped or circular aperturein a horizontal rigid dock. If _n and _n are the decay constantsand magnetic scalar potentials for the strip and _n and _n theangular frequencies and velocity potentials for the normal modesin the strip-shaped aperture, then _n =_n² and _n and _n are thereal and imaginary parts respectively of a holomorphic function.The velocities in the normal modes are deduced from the solutionof the eddy-current problem and are found to agree with resultsgiven in Miles (1972). For circular geometries the eigenvaluesand eigenfunctions of the axisymmetric eddy-current problemare the same as those of the water-wave problem that has angularvariation eⁱ; where (, , z) are cylindrical polar co-ordinateslocated at the centre of the basin.     

Convergence of Pad? Kernel Approximants to the Delayed Dirac Function

We derive from the solutions of Kummer's equation an expressionfor the exponential function, which when restricted to integerparameters gives the [M, N] Pad? approximation and its error.Laplace inversion, using Bromwich's integral, then yields anapproximation of the form to ( – t) the delayed Dirac kernel, in which the constants{_{i, v}, K_i,v: i = 1, 2 , ...,N, = M – N + 1} are thoseof the partial fraction decomposition of the [M, N] Pad? approximation.These constants also occur in direct quadrature formulae toinvert Laplace transforms. Finally, we show that the kernelapproximants converge for > t.     

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 effects of an electric field on an autocatalytic ionic reaction in a system with high ionic strength

The effects of an applied electric field on an ionic autocatalyticreaction with a quadratic rate law are considered, where thereacting species, A⁺ and B⁺, are present in a system which alsoincludes non-reacting species C^- and D⁺. The conditions areestablished under which the general terms which describe theelectric field effects in the reaction-diffusion equations canbe simplified to those used in previous studies, where theseeffects are modelled by linear advection terms. The resultingequations are then studied in detail by first obtaining conditionsfor the existence of travelling waves of permanent form. Thisdiscussion shows that _B, the ratio of the diffusion coefficientsof B⁺ and A⁺, is a critical parameter, with different formsof behaviour arising for _B < 1 and _B > 1. This analysisis augmented by obtaining solutions valid for large times andlarge values of (the dirnensionless applied field). Numericalsolutions of initial-value problems are obtained for a rangeof values of and _B, guided by and interpreted through the analysispreviously obtained. These numerical integrations show the formationof reaction fronts, with the possibility of greatly increasedreaction rates caused by the applied electric field, as wellas propagating electrophoretic fronts in B⁺ being formed incases where a reaction front is also initiated. There is alsothe possibility of separate electrophoretic fronts in A⁺ andB⁺ being formed, which become increasingly separated as timeincreases with the reaction being completely inhibited.     

A Phase-Change Model with a Zone of Coexistence of Phases

A mathematical model for change of phase is presented, accountingfor the presence of regions in which liquid and solid coexist.The basic variables are temperature and solid fraction v. Westart from a relationship of the type =(v), supposed valid inthermodynamical equilibrium. Then for dynamical processes weintroduce a perturbation causing v to be less than its equilibriumvalue in any solidification process. This solid fraction deficiencyis governed by an ordinary differential equation containing_t, in the forcing term. The heat-balance equation is in turncoupled to the ordinary differential equation through the termv_t, ( is latent heat). Some existence and uniqueness resultsare proved and some monotonicity properties are described forpure melting or pure solidification processes.     

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.     

Analysis of Certain Slowly Converging Series

Slowly converging series frequently arise in the solution ofdifferential equations. The slow convergence is commonly dueto the presence of some kind of singularity. This paper demonstratesthat, if sufficient information is available as to the natureof the singularity, it is possible to accelerate massively theconvergence by extracting most singular terms. The prime illustration is a separating laminar boundary layer,where there is a weak singularity at separation if the pressureis independently prescribed. This singularity was extensivelystudied by Goldstein (1948), Stewartson (1958) and Terrill (1960),who showed that the leading terms in an expansion about thesingularity depend upon a single parameter ₁, with later termsdepending upon parameters ₅, ₉, etc. In a particular case, first studied by Howarth (1938), a computersolution by Leigh (1955) gave ₁ = 0.492, and more recently Davies& Walker (1977) gave the (presumably) more accurate ₁ =0.502. The present work uses knowledge of the singularity toaccelerate the convergence of a series expansion (in powersof distance from the leading edge) and yields ₁ = 0.50287;though the last figure may not be reliable. Further, thoughthe computer solutions are not sufficiently accurate to estimate₅, the present method yields ₅ =– 2.0(6). Other detailsof the flow are predicted, and compared favourably with thecomputer solutions. The above major claims assume that the position of the singularityis known precisely. Additional calculations, where the positionof the singularity is first to be estimated, still yield excellentresults, except very close to the singularity. Even here, theyare remarkably accurate.     

Tensor Products of C*-Algebras Over Abelian Subalgebras

Suppose that A is a C*-algebra and C is a unital abelian C*-subalgebrawhich is isomorphic to a unital subalgebra of the centre ofM(A), the multiplier algebra of A. Letting = , so that we maywrite C = C(), we call A a C()-algebra (following Blanchard[7]). Suppose that B is another C()-algebra, then we form A_CB, the algebraic tensor product of A with B over C as follows:A B is the algebraic tensor product over C, I_C = {ⁿ_i–1(f_i 1–1f_i)x|f_iC, xAB} is the ideal in AB generated by f1–1f|fC,and A _CB = AB/I_C. Then A_CB is an involutive algebra over C,and we shall be interested in deciding when A_CB is a pre-C*-algebra;that is, when is there a C*-norm on A_C B? There is a C*-semi-norm,which we denote by ||·||_C-min, which is minimal in thesense that it is dominated by any semi-norm whose kernel containsthe kernel of ||·||_C-min. Moreover, if A _C B has a C*-norm,then ||·||_C-min is a C*-norm on A_C B. The problem isto decide when ||·||_C-min is a norm. It was shown byBlanchard [7, Proposition 3.1] that when A and B are continuousfields and C is separable, then ||·||_C-min is a norm.In this paper we show that ||·||_C-min is a norm whenC is a von Neumann algebra, and then we examine some consequences.     

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.