This paper addresses the effect of general Lewis number and heat losses on the calculation of combustion wave speeds using an asymptotic technique based on the ratio of activation energy to heat release being considered large. As heat loss is increased twin flame speeds emerge (as in the classical large activation energy analysis) with an extinction heat loss. Formulae for the non-adiabatic wave speed and extinction heat loss are found which apply over a wider range of activation energies (because of the nature of the asymptotics) and these are explored for moderate and large Lewis number cases—the latter representing the combustion wave progress in a solid. Some of the oscillatory instabilities are investigated numerically for the case of a reactive solid.

In this paper, we study the asymptotic behavior of an SIRS epidemic model with a time delay in the recovered class and a nonlinear incidence rate. A conjecture of Hethcote et al. [5] on the global stability of the disease-free equilibrium is solved. Moreover, we analyse the model when the contact number takes its threshold value. We show that solutions tend to either the disease-free equilibrium or to a unique positive endemic equilibrium, and there is no periodic solution.

The present paper proves criteria for oscillation of the soultions of functional differential equations of the type

where λ, τ > 0.

Closed-form analytical expressions are derived for the reflection and transmission coefficients for the problem of scattering of surface water waves by a sharp discontinuity in the surface-boundary-conditions, for the case of deep water. The method involves the use of the Havelock-type expansion of the velocity potential along with an analysis to solve a Carleman-type singular integral equation over a semi-infinite range. This method of solution is an alternative to the Wiener-Hopf technique used previously.

Here we consider a particular class of stochastic geometric programs in which the randomness occurs in the decision variables. Specifically we analyse a program in which we specify a joint normal probability for the dicision variables and require the constraint set to be satisfied in the chance constrained mode. A numerical example is given to illustrate the approach.

A monoenergetic point-source solution of the steady-state cosmic-ray equation of transport for cosmic-rays in the interplanetary region in which monoenergetic particles are released isotropically and continuously from a fixed heliocentric position is derived by Lie Theory. A spherically-symmetric model of the propagation region is assumed incorporating anisotropic diffusion with a diffusion tensor symmetric about the radial direction, and the solar wind velocity is radial and of constant speed V. Because of the point release the solution is non-spherically-symmetric.

In this paper, we use an ordinary differential equation approach to study the existence of similarity solutions for the equation u1 = Δ(uα) + θu–β in Rn × (0, ∞) where β > 0, θ ∈ [0, 1}, and n ≥ 1. This includes the slow diffusion equation when α > = 1, and the standard heat equation when α = 1, and the fast diffusion equation when 0 < α < 1. We prove that there are forward self-similar solutions for this equation with initial data of the form c|x|p, where p = 2/(α + β) if θ = 1; p ≥ 0 and 2 + (1 – α)p > 0 if θ = 0, for some positive constant c.

We have finally obtained for each of the 6 Painlevés an expression of z, w, w′ that behaves as 1/(z − Z0) + O(1) at each kind of movable singular point. This expression is polynomial in w′ (at most quadratic), and rational in w and z. After it is integrated and exponentiated it yields a function that has a simple zero at each of the singular points.

We refine Shannon's inequality, in its discrete and integral forms, by presenting upper estimates of the difference between its two sides.

Under the assumptions that the spatial variable is one dimensional and the distributed delay kernel is the general Gamma distributed delay kernel, when the average delay is small, the existence of travelling wave solutions for the population genetics model with distributed delay is obtained by using the linear chain trick and geometric singular perturbation theory. On the other hand, for the population genetics model with small discrete delay, the existence of travelling wave solutions is obtained by employing a technique which is based on a result concerning the existence of the inertial manifold for small discrete delay equations.

We investigate oscillatory behaviour in the famous Belousov-Zhabotinskii chemical reaction, as described by the simple two-variable Oregonator model. It is shown that oscillations are possible only in certain parameter regions. Numerical results are presented, and the presence of fold bifurcations discussed.

In this paper we have introduced extensions νυ(α, x; b) and Γν(α, x; b) of the generalized Gamma functions γ(α x; b) and Γ(α, x; b) considered recently by Chaudhry and Zubair. These extensions are found useful in the representations of the Laplace and K-transforms of a class of functions. We have also defined a generalization of the inverse Gaussian distribution. The cumulative and the reliability functions of the generalized inverse Gaussian distribution are expressed in terms of these functions. Some useful properties of the functions are also discussed.

The paper extends earlier work by using the factorisation method to discuss solutions of period four for the difference equation

This equation was suggested by R. M. May as a simple mathematical model for the effect of frequency-dependent selection in genetics. It is shown that for a given value of the parameter, a, the identification of solutions of period four can be reduced to finding real roots for a polynomial equation of degree eight. The appropriate values of xn follow from a quartic equation. By splitting up the problem in this way it becomes relatively straightforward to determine the critical values of a at which the various solutions of period four first appear and to discuss the stability of these solutions. Intervals of stability are tabulated in the paper.

Newton's method is applied to an operator that satisfies stronger conditions than those of Kantorovich. Convergence and error estimates are compared in the two situations. As an application, we obtain information on the existence and uniqueness of a solution for differential and integral equations.

Viscous fluid is squeezed out from a shrinking (or expanding) tube whose radius varies with time as (1 – βt)½. The full Navier–Stokes equations reduce to a non-linear ordinary differential equation governed by a non-dimensional parameter S representing the relative importance of unsteadiness to viscosity. This paper studies the analytic solutions for large | S | through the method of matched asymptotic expansions. A simple numerical scheme for integration is presented. It is found that boundary layers exist near the walls for large | S |. In addition, flow reversals and oscillations of the velocity profile occur for large negative S (fast expansion of the tube).

The Maxwell-Dirac equations model an electron in an electromagnetic field. The two equations are coupled via the Dirac current which acts as a source in the Maxwell equation, resulting in a nonlinear system of partial differential equations (PDE's). Well-behaved solutions, within reasonable Sobolev spaces, have been shown to exist globally as recently as 1997 [12]. Exact solutions have not been found—except in some simple cases.

We have shown analytically in [6, 18] that any spherical solution surrounds a Coulomb field and any cylindrical solution surrounds a central charged wire; and in [3] and [19] that in any stationary case, the surrounding electron field must be equal and opposite to the central (external) field. Here we extend the numerical solutions in [6] to a family of orbits all of which are well-behaved numerical solutions satisfying the analytic results in [6] and [11]. These solutions die off exponentially with increasing distance from the central axis of symmetry. The results in [18] can be extended in the same way. A third case is included, with dependence on z only yielding a related fourth-order ordinary differential equation (ODE) [3].