The Mesa Problem: Diffusion Patterns for ut = {nabla} {middle dot} (um{nabla}u) as m -> +{infty}

In this paper we consider the limit m+ of solutions of the porous-mediumequation u_t = · (u^mu) (xR^N), with N > 1. We conjecturethat, for initial data with a unique maximum, the evolutionis characterized by the onset of a ‘mesa’ region,in which the solution is nearly spatially independent, surroundedby a region in which u is nearly equal to its initial value.The transition between these regions occurs near a surface whichis identified with the free boundary in a certain Stefan problemwhich can be studied using variational inequalities. Moreover,singular-perturbation theory can be used to describe the structureof the transition region.     

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.     

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

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.     

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

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.     

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.     

Blow up of solutions of a generalized Boussinesq equation

Consider the Cauchy problem u_tt = (f(u))_xx + u_xxtt x R, t 0, u(x,0) = u₀(x), u_t(x,0) = u₁(x),7rcub; where f : R R C, f(0) = ). After treatment of the local existenceproblem, we show the blow up of the solution of the equation(1) under the folowing assumptions. Let > 0 be real, such that 2(l + 2)F(u) uf(u), (v₀, Pv₀)_l² + _- F(u₀)dx < 0 where P = 1 - d²/dx², F'(s), and v₀ is given by u₁(x,0) = (v₀(x))_x. Then we focus on various perturbations of the question. We alsostudy the vectorial case in the same way, and finally we giveexamples.     

Exact Remainders for Asymptotic Expansions of Fractional Integrals

This is a continuation of work begun in an earlier paper inwhich we used the theory of distributions to derive explicitexpressions for the remainder terms associated with the asymptoticexpansions of the Stieltjes transform. In this paper similarresults are obtained for the fractional integral of order definedby 1f(x)=1/f()^x_o(x-)^x-1 f(t)dt, >. Heref(t) is a locally integrable function on [0, ) and satisfies f(t) a_st-5–0(ó >0), s=0 as      

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.     

A Comparison Between Chebyshev and Ultraspherical Expansions

This paper examines the norms of the projections from C[–1,1] to P_M[–1, 1] which are obtained by truncating Chebyshevand ultraspherical expansions after n terms. The norms are evaluatednumerically for n = 1, 2, 3 , ... , 10 and –0.1 5 wherethe ultraspherical weight function is (1 – x²)^x–.In particular the data show that the norm of the Chebyshev projectionis not the smallest in this range for 1 n 10.     

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.     

Optimal Parameter Choice in the {theta}-Method for Parabolic Equations

The general first-order method, known as the -method, is appliedto the semi-discrete form of a parabolic equation. It is shownthat to every required local accuracy there corresponds a valueof the parameter that is optimal in the sense of allowing thelargest step for which the error remains bounded below . Anasymptotic formula for in terms of is obtained, showing thatthe maximum step-size for the optimal -method is more than twiceas large as that for the Crank-Nicolson method. A numericalexample is given, showing good agreement between theory andpractice.     

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 .     

Self-similar solutions to the isothermal compressible Navier-Stokes equations

^** Email: guo_zhenhua{at}iapcm.ac.cn^*** Email: jiang{at}iapcm.ac.cn We investigate the self-similar solutions to the isothermalcompressible Navier–Stokes equations. The aim of thispaper is to show that there exist neither forward nor backwardself-similar solutions with finite total energy. This generalizesthe results for the incompressible case in Neas, J., Rika, M.& verák, V. (1996, On Leray's self-similar solutionsof the Navier-Stokes equations. Acta. Math., 176, 283–294),and is consistent with the (unproved) existence of regular solutionsglobally in time for the compressible Navier–Stokes equations.     

Spreading speed and travelling wave solutions of a partially sedentary population

In this paper, we extend the population genetics model of Weinberger(1978, Asymptotic behavior of a model in population genetics.Nonlinear Partial Differential Equations and Applications (J.Chadam ed.). Lecture Notes in Mathematics, vol. 648. New York:Springer, pp. 47–98.) to the case where a fraction ofthe population does not migrate after the selection process.Mathematically, we study the asymptotic behaviour of solutionsto the recursion u_n+1 = Q_g[u_n], where In the above definition of Q_g, K is a probabilitydensity function and f behaves qualitatively like the Beverton–Holtfunction. Under some appropriate conditions on K and f, we showthat for each unit vector R^d, there exists a c^*_g() which hasan explicit formula and is the spreading speed of Q_g in thedirection . We also show that for each c c^*_g(), there existsa travelling wave solution in the direction which is continuousif gf '(0) 1.     

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.