Steady two-dimensional free surface flow past a semi-infinite flat plate is considered. The vorticity in the flow is assumed to be constant. For large values of the Froude number F, an analytical relation between F, the vorticity parameter ω and the steepness s of the waves in the far field is derived. In addition numerical solutions are calculated by a boundary integral equation method.

In this paper, a closed-form analytical solution for pricing convertible bonds on a single underlying asset with constant dividend yield is presented. A closed-form analytical formula has apparently never been found for American-style convertible bonds (CBs) of finite maturity time although there have been quite a few approximate solutions and numerical approaches proposed. The solution presented here is written in the form of a Taylor's series expansion, which contains infinitely many terms, and thus is completely analytical and in a closed form. Although it is only for the simplest CBs without call or put features, it is nevertheless the first closed-form solution that can be utilised to discuss convertibility analytically. The solution is based on the homotopy analysis method, with which the optimal converting price has been elegantly and temporarily removed in the solution process of each order, and consequently, the solution of a linear problem can be analytically worked out at each order, resulting in a completely analytical solution for the optimal converting price and the CBs' price.

A centre manifold or invariant manifold description of the evolution of a dynamical system provides a simplified view of the long term evolution of the system. In this work, I describe a procedure to estimate the appropriate starting position on the manifold which best matches an initial condition off the manifold. I apply the procedure to three examples: a simple dynamical system, a five-equation model of quasi-geostrophic flow, and shear dispersion in a channel. The analysis is also relevant to determining how best to account, within the invariant manifold description, for a small forcing in the full system.

The machinery of Lie theory (groups and algebras) is applied to the unsteady equations of motion of rotating fluid. A special-function type solution for the steady state is derived. It is then shown how the solution generates an infinite number of time-dependent solutions via three arbitrary functions of time. This algebraic structure also provides the mechanism to search for other solutions since its character is inferred from the basic equations.

Carlson has shown that if the predicted price in the linear cobweb model is taken as the average of all previous actual prices, then stability results independently of parameter values provided only that the demand–curve gradient is less than that of the supply curve. This result has subsequently been generalised by Manning and by Holmes and Manning. We investigate the robustness of their results.

The aggregation–decomposition method is used to derive a sufficient condition for the equi-ultimate boundedness of large-scale systems governed by nonlinear ordinary differential equations.

In this paper, a gradient method is developed for the optimal shape design in a nozzle problem described by variational inequalities. It is known that this method can be used for the optimal shape design for systems described by partial differential equations (Pironneau [6]); it is used here for differential inequalities by taking limits of the expression resulting from an approximations scheme. The computations are done by the finite element method; the gradient of the criteria as a function of the coordinates nodes is computed, and the performance criterion is then minimised by the gradient method.

Centre manifolds arise in a rational approach to the problem of forming low-dimensional models of dynamical systems with many degrees of freedom. When a dynamical system with a centre manifold is subject to a small forcing, F, there are two effects: to displace the centre manifold; and to alter the evolution thereon. We propose a formal scheme for calculating the centre manifold of such a forced dynamical system. Our formalism permits the calculation of these effects, with errors of order |F|2. We find that the displacement of the manifold allows a reparameterisation of its description, and we describe two “natural” ways in which this can be carried out. We give three examples: an introductory example; a five-mode model of the atmosphere to display the quasi-geostrophic approximation; and the forced Kuramoto-Sivashinsky equation.

In this paper we consider simplifying a model of the nitrogen cycle in Port Phillip Bay, Victoria, Australia. The approach taken is to aggregate state variables that are linearly related using a projection in state space. The technique involved is a modification of proper orthogonal decomposition and was developed so that a resulting simplified model retains an ecological interpretation. It can be applied automatically, and enables insights into the system to be gained that were not obvious beforehand. In the case of the Port Phillip Bay model, we find that the variables representing water and sand are unaffected by the remaining variables, while only variables on the same trophic level can be grouped together. The validity of the aggregation under several nutrient loads is also discussed.

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.

Many tens of serious incidents involving reactors occur in the developed countries each year. The disaster at the chemical plant in Bhopal, India in 1984 was particularly notable where a thermal runaway process led to more than 3000 tragic fatalities from the cloud of extremely toxic methyl isocyanate that boiled out of a storage tank. This signalled the design of special types of chemical reactors to reduce the risk of thermal runaway by planning (at the design stage) integral safety and thermal stabilization mechanisms. The Endex CSTR (continuously stirred tank reactor) proposed by Gray and Ball [3] involves a reactor in two parts with heat exchange allowed between them. The two parts of the reactor operate side by side in tandem, such that the thermal runaway of one part is offset by an endothermic reaction in the other reactor—hence the term ‘endex’.

It is found that the adiabatic endex system has a large region of parameter space where the operation can be made safe. However adiabatic conditions rely on the continuous supply of reactants to the endothermic side of the reactor, for operation of the system. The risks involved are such that it is always safer to operate batch reactors in a non-adiabatic mode. Thus we consider the limiting case of the approach to adiabatic conditions where although the mathematics produces no oscillatory causes for instability, yet there is a narrow but significant area where the stable solution branch is lost and consequently a persistent and unexpected region of instability in what otherwise appears to be a simple CSTR system.

A simple algebraic method is presented to determine the necessary condition for the existence of a Hamiltonian circuit in a directed graph of n vertices. A search procedure is then introduced to identify any or all of the existing Hamiltonian circuits. The procedure is based upon finding a set of edges which will then be candidates for being parts of circuits of length n at any vertex of the graph.

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.

The problem of transient two-dimensional transport by diffusion and advection of a decaying contaminant in two adjacent porous media is solved using a boundary-integral method. The method requires the construction of appropriate Green's functions. Application of Green's theorem in the plane then yields representations for the contaminant concentration in both regions in terms of an integral of the initial concentration over the region's interior and integrals along the boundaries of known quantities and the unknown interfacial flux between the two adjacent media. This flux is given by a first-kind integral equation, which can be solved numerically by a discretisation technique. Examples of contaminant transport in fractured porous media systems are presented.

With the help of the generalized Mayer theorem we obtain an improved inequality for free energies of model and approximating systems, where only “connected parts” over the approximating Hamiltonian are taken into account. For a concrete system we discuss the problems of convergence of appropriate series of “connected parts ”.

The likelihood ratio approach to the detection of small signals in the presence of noise is investigated in the case where the available data have been clipped. The statistic obtained is the ratio of orthant probabilities and appears intractable; accordingly an approximation to this statistic is developed by truncating an appropriate Taylor expansion. Approximations are obtained for the mean and variance of this modified statistic and compared with those obtained from computer simulations.

The temporal instability of a developing swirling incompressible jet is considered. The jet development (in the streamwise direction) is modelled by combining a near-field and far-field approximation to the jet velocity profile into a one parameter family of basic velocity fields. The single parameter in the jet velocity field then allows us to model the radial spreading of the jet and the decay of swirl observed experimentally. Two distinct modes of instability of this model profile are found. The first is that found from a stability analysis of a fully developed swirling jet in the far field whilst the second is relevant to a “top-hat” jet with an imposed rigid body rotation. We demonstrate that the effect of azimuthal swirl is to destabilise both modes of instability. Additionally our results suggest that the near-nozzle modes of instability will dominate; indeed the growth rates of these modes are significantly larger than those found from previous studies of a fully developed jet in the far-field region.

Using geometric quantization, and accepting the quantum Hamiltonian of previous authors, we propose some candidate formulae for the quantum operator of an observable which is a quadratic form in the momenta.