Experimental evidence shows that plane Couette flow becomes unstable when the Reynolds number R reaches certain critical values. Linear stability theory does not predict these observations and has been unable to locate these instabilities. A Chebyshev/QR numerical technique is used to investigate much higher values of R than those previously tested. In particular, values of R up to 108 are confidently tested, whereas previously values of R up to only 2 × 104 have been considered.

The minimization of signal distortion was one approach applied successfully in the theory of optimum signal detection for arrays [3]. The processors considered operated on the input as received.

In some applications it is desirable to clip the received signal before processing and the problem of optimum processing of such clipped signals then arises. Several approaches to this problem are being studied, but the present paper is concerned with that based on minimum signal distortion.

The paper discusses solutions of period 4 for the difference equation

where k and m are real parameters, with k > 0. For given values of k and m there are at most three solutions with period 4 and equations are set up to determine the elements of these solutions and the stability of each solution. Only real solutions are considered. The procedure that is used to find these solutions allows unstable solutions to be identified as well as stable solutions.

In a previous paper, solutions of period 2 and period 3 were examined for this equation and there was evidence of anomalous behaviour in the way the stability intervals occurred. Some preliminary information about solutions of period 4 was mentioned in the discussion. The present paper provides more complete results, which confirm the anomalous behaviour and give a better idea of how the stability criterion changes for different families of solutions. These results are used to indicate the variety of behaviour that can be found for one-parameter systems by imposing suitable conditions on m and k.

Exact solutions are developed for instantaneous point sources subject to nonlinear diffusion and loss or gain proportional to nth power of concentration, with n > 1. The solutions for the loss give, at large times, power-law decrease to zero of slug central concentration and logarithmic increase of slug semi-width. Those for gain give concentration decreasing initially, going through a minimum, and then increasing, with blow-up to infinite concentration in finite time. Slug semi-width increases with time to a finite maximum in finite time at a blow-up. Taken in conjunction with previous studies, these new results provide an overall schema for instantaneous nonlinear diffusion point sources with nonlinear loss or gain for the total range n ≥ 0. Six distinct regimes of behaviour of slug semi-width and concentration are identified, depending on the range of n, 0 ≤ n < 1, n = 1, or n > 1. Three of them are for loss, and three for gain. The classical Barenblatt-Pattle nonlinear instantaneous point-source solutions with material concentration occupy a central place in the total schema.

Taylor's model of dispersion simply describes the long-term spread of material along a pipe, channel or river. However, often we need multi-mode models to resolve finer details in space and time. Here we construct zonal models of dispersion via the new principle of matching their long-term evolution with that of the original problem. Using centre manifold techniques this is done straightforwardly and systematically. Furthermore, this approach provides correct initial and boundary conditions for the zonal models. We expect the proposed principle of matched centre manifold evolution to be useful in a wide range of modelling problems.

The large-amplitude oscillations and buckling of an anisotropic cylindrical shell subjected to the initial inplane biaxial normal stresses have been analysed. The concept of anisotropy used by Lekhnitsky has been introduced into the field equations for cylindrical shells of isotropic material deduced by Donnell. The method of Galerkin and the method of successive approximation have been used to obtain the desired approximate solution. The expression for the critical loads for the buckling of anisotropic cylindrical shells has been obtained during intermediate stages of analysis. Some relevant frequency response graphs of the obtained solution are also presented. The minimum critical loads for various classes of anisotropy have also been given at the end of the discussion, to exhibit the effects of large deflections and imperfections on elastic buckling.

The scattering and diffraction of harmonic SH waves by an arbitrarily shaped alluvial valley in a layered material is considered. The problem is solved in terms of boundary integral equations which yield a numerical solution.

In an earlier paper [4], the author showed how Laplace transforms might be assigned to a class of superexponential functions for which the usual defining integral diverges. The present paper considers the case of the function exp(et), which arises in combinatorial contexts and whose Laplace transform may be assigned by means of an extension of techniques described in the previous paper.

A problem in combustion theory with tenperature-dependent conductivity is considered. It is shown that information regarding criticality dependence on data and parameters can be obtained from a transformed equation in which the conductivity is constant, while the nonlinear source term is modified. Some previous work can then be used in such a study.

In this note we consider various theoretical aspects of the problem of least-squares approximation subject to constraints on the range of the approximating polynomial. The problem is treated from an optimization theory viewpoint. Rice's parameter space procedure is discussed.

We provide a qualitative analysis of a system of nonlinear differential equations that model the spread of alcoholism through a population. Alcoholism is viewed as an infectious disease and the model treats it within a sir framework. The model exhibits two generic types of steady-state diagram. The first of these is qualitatively the same as the steady-state diagram in the standard sir model. The second exhibits a backwards transcritical bifurcation. As a consequence of this, there is a region of bistability in which a population of problem drinkers can be sustained, even when the reproduction number is less than one. We obtain a succinct formula for this scenario when the transition between these two cases occurs.

We investigate the approximate controllability of a size- and space-structured population model, for which the control function acts on a subdomain and corresponds to the migration of individuals. We establish the main result via the unique continuation property of the adjoint system. The desired controller is the minimizer of an infinite-dimensional optimization problem.

Conditions are given for a Ck map T to be a Newton map, that is, the map associated with a differentiable real-valued function via Newton's method. For finitely differentiable maps and functions, these conditions are only necessary, but in the smooth case, that is, for k = ∞, they are also sufficient. The characterisation rests upon the structure of the fixed point set of T and the value of the derivative T′ there, and it is best possible as is demonstrated through examples.

The boundary integral equation method is obtained by expressing a solution to a particular partial differential equation in terms of an integral taken round the boundary of the region under consideration. Various methods exist for the numerical solution of this integral equation and the purpose of this note is to outline an improvement to one of these procedures.

A model governing the combustion of a material is considered. The model consists of two non-linear coupled parabolic equations with initial and boundary conditions. An approximation for the rate of reactant consumption is made to enable the system to the treated by laplace transform. Three simple geometries are considered; namely, an infinite slab, an infinite circular and a sphere. The results obtained are then compared with numerical solutions for spme specific values of the parameters. There is good agreement over time duration for which numerical work was performed.

We consider the stability of high Reynolds number flow past a heated, curved wall. The influence of both buoyancy and curvature, with the appropriate sense, can render a flow unstable to longitudinal vortices. However, conversely each mechanism can make a flow more stable; as with a stable stratification or a convex curvature. This is partially due to their influence on the basic flow and also due to additional terms in the stability equations. In fact the presence of buoyancy in combination with an appropriate local wall gradient can actually increase the wall shear and these effects can lead to supervelocities and the promotion of a wall jet. This leads to the interesting discovery that the flow can be unstable for both concave and convex curvatures. Furthermore, it is possible to observe sustained vortex growth in stably stratified boundary layers over convexly curved walls. The evolution of the modes is considered in both the linear and nonlinear régimes.

We study the existence of the invariant region for the equations of nonisentropic gas dynamics. We obtain the mean-integral of the conserved quantity after making an intensive study of the Riemann problem. Using the extremum principle and the Lagrangian multiplier method, we prove that the one-dimensional equations of nonisentropic gas dynamics for an ideal gas possess a unique invariant region. However, the invariant region is not bounded.

A three stage procedure for the analysis and least-cost design of looped water distribution networks is considered in this paper. The first stage detects spanning trees and identifies the true global optimum for the system. The second stage determines hydraulically feasible pipe flows for the network by the numerical solution of a set of non-linear simultaneous equations and shows that these solutions are contained within closed convex polygonal regions in the solution space bounded by singularities resulting from zero flows in individual pipes. Ideal pipe diameters, consistent with the pipe flows and the constant velocity constraint adopted to prevent the system degenerating into a branched network, are selected and costed. It is found that the most favourable optimum is in the vicinity of a vertex in the solution space corresponding to the minimum spanning tree. In the third stage, commercial pipes are specified and the design finalised. Upper bound formulae for the number of spanning trees and hydraulically feasible solutions in a network have also been proposed. The treatment of large networks by a heuristic procedure is described which is shown to result in significant economies compared with designs obtained by non-linear programming.

Elastic behaviour of a nonhomogeneous transversely isotropic half-space is studied under the action of a smooth rigid axisymmetric indentor. Hankel transforms of different orders have been used. It is observed that in contrast to a homogeneous medium, the pressure distribution in the contact region in a nonhomogeneous medium is not directly available, rather it is obtainable from the solution of a Fredholm integral equation. The integral equation is solved for a flat-ended punch and paraboloidal indentations for various values of the nonhomogeneity parameter, and the effects of nonhomogeneity in elastic behaviour on stresses have been shown graphically. The results of the associated homogeneous case are readily available from the results of the present study.