In this paper, we extend the results in the literature for boundary layer flow over a horizontal plate, by considering the buoyancy force term in the momentum equation. Using a similarity transformation, we transform the partial differential equations of the problem into coupled nonlinear ordinary differential equations. We first analyse several special cases dealing with the properties of the exact and approximate solutions. Then, for the general problem, we construct series solutions for arbitrary values of the physical parameters. Furthermore, we obtain numerical solutions for several sets of values of the parameters. The numerical results thus obtained are presented through graphs and tables and the effects of the physical parameters on the flow and heat transfer characteristics are discussed. The results obtained reveal many interesting behaviours that warrant further study of the equations related to non-Newtonian fluid phenomena, especially the shear-thinning phenomena. Shear thinning reduces the wall shear stress.

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 consider fully three-dimensional time-dependent outflow from a source into a surrounding fluid of different density. The source is distributed over a sphere of finite radius. The nonlinear problem is formulated using a spectral approach in which two streamfunctions and the density are represented as a Fourier-type series with time-dependent coefficients that must be calculated. Linearized theories are also discussed and an approximate stability condition for early stages in the outflow is derived. Nonlinear solutions are presented and different outflow shapes adopted by the fluid interface are investigated.

The recent work of Cheng and Stokes on the processing of clipped signals from two or three receivers is extended and generalised by removing a number of restrictions. In particular, there is no restriction on the number of receivers and the restrictions on the statistical properties of the signal and noise processes have been considerably relaxed.

Mathematically—Plackett's result is used to expand the orthant proabilities involved in increasing powers of the input signal to noise ratio.

The explicit inverse and determinant of a class of matrices is given. The class is the Hadamard product of two already known classes. Its elements are defined by 3n − 1 parameters, analytical expressions of which compose the Hessenberg form inverse. These expressions enable a recursive formula to be obtained, which gives the inverse in O(n2) multiplications/divisions and O(n) additions/subtractions.

We shall study continuous-time Markov chains on the nonnegative integers which are both irreducible and transient, and which exhibit discernible stationarity before drift to infinity “sets in”. We will show how this ‘quasi’ stationary behaviour can be modelled using a limiting conditional distribution: specifically, the limiting state probabilities conditional on not having left 0 for the last time. By way of a dual chain, obtained by killing the original process on last exit from 0, we invoke the theory of quasistationarity for absorbing Markov chains. We prove that the conditioned state probabilities of the original chain are equal to the state probabilities of its dual conditioned on non-absorption, thus allowing to establish the simultaneous existence and then equivalence, of their limiting conditional distributions. Although a limiting conditional distribution for the dual chain is always quasistationary distribution in the usual sense, a similar statement is not possible for the original chain.

When two trains travel along the same track in the same direction, it is a common safety requirement that the trains must be separated by at least two signals. This means that there will always be at least one clear section of track between the two trains. If the safe-separation condition is violated, then the driver of the following train must adopt a revised strategy that will enable the train to stop at the next signal if necessary. One simple way to ensure safe separation is to define a prescribed set of latest allowed section exit times for the leading train and a corresponding prescribed set of earliest allowed section entry times for the following train. We will find strategies that minimize the total tractive energy required for both trains to complete their respective journeys within the overall allowed journey times and subject to the additional prescribed section clearance times. We assume that the drivers use a discrete control mechanism and show that the optimal driving strategy for each train is defined by a sequence of approximate speedholding phases at a uniquely defined optimal driving speed on each section and that the sequence of optimal driving speeds is a decreasing sequence for the leading train and an increasing sequence for the following train. We illustrate our results by finding optimal strategies and associated speed profiles for both trains in some elementary but realistic examples.

This paper examines the issue of derivative pricing within the framework of a fractional stochastic volatility model. We present a deterministic partial differential equation system to derive an approximate expression for the derivative price. The proposed approach allows for the stochastic volatility to be expressed as a composition of deterministic functions of time and a fractional Ornstein–Uhlenbeck process. We apply this method to the European option pricing under the fractional Stein–Stein volatility model, demonstrating its feasibility and reliability through numerical simulations. Our numerical simulations also illustrate the impact of the parameters in the fractional stochastic volatility model on the option price.

Many electricity market participants have a requirement to calculate the probabilistic risk measures, such as earnings at risk (EaR) and value at risk (VaR), for compliance reporting purposes. This requirement is currently hindered by the lack of analytical representations for forecasts of demand (load) and price curves; this motivates numerical simulation and models that need extensive calibration. In this paper, we derive an analytical representation of a state demand forecast which is the aggregated usage of all electricity consumers in a particular region (such as New South Wales or Victoria). We have used two probabilistic benchmarks from the Australian energy market operator as input, which are expressed as forecasted probability of exceedance.

Due to a number of considerations, including asymmetry of these quantiles with respect to the median, we have selected a series of truncated lognormal distributions with two parameters. The procedure of finding these parameters has been reduced to solving (for every half-hour) a single nonlinear equation. As a result, the two-year half-hourly forecast (expected curve) and demand volatility are found by explicit integration with the set of derived distributions. We have also tested an alternative method based on simplifying assumptions; using a nontruncated lognormal distribution, we found that under the test conditions this method produces an identical forward load and volatility curve.

Pontryagin's maximum principle is derived by elementary mathematical techniques. The conditions on the functions which enter are generally somewhat more stringent than in Pontryagin's derivation, but one (practically very awkward) condition of Pontryagin can be relaxed: continuity in the time variable can be replaced by a much weaker condition.

We discuss simulation of sensitivities or Greeks of multi-asset European style options under a special Lévy process model: that is, the subordinated Brownian motion model. The Malliavin calculus method combined with Monte Carlo and quasi-Monte Carlo methods is used in the simulations. Greeks are expressed in terms of the expectations of the option payoff functions multiplied by the weights involving Malliavin derivatives for multi-asset options. Numerical results show that the Malliavin calculus method is usually more efficient than the finite difference method for options with nonsmooth payoffs. The superiority of the former method over the latter is even more significant when both are combined with quasi-Monte Carlo methods.

The resolution of the identity formula for a localisation operator with two admissible wavelets on a separable and complex Hilbert space is given and the traces of these operators are computed.

In this paper a feasible direction method is presented to find all efficient extreme points for a special class of multiple objective linear fractional programming problems, when all denominators are equal. This method is based on the conjugate gradient projection method, so that we start with a feasible point and then a sequence of feasible directions towards all efficient adjacent extremes of the problem can be generated. Since methods based on vertex information may encounter difficulties as the problem size increases, we expect that this method will be less sensitive to problem size. A simple production example is given to illustrate this method.

Multiple integrals in ten or twenty variables are often needed by atomic, molecular and nuclear physicists, because of the large number of degrees of freedom in the quantum systems with which they must deal. In statistics too there is often a need to evaluate integrals with many degrees of freedom. It is in mathematical finance, however, that the most striking examples are seen, with claims of integrals being evaluated during recent years with many hundreds of variables.

Microbial competition for nutrients is a common phenomenon that occurs between species inhabiting the same environment. Bioreactors are often used for the study of microbial competition since the number and type of microbial species can be controlled, and the system can be isolated from other interactions that may occur between the competing species. A common type of competition is the so-called “pure and simple” competition, where the microbial populations interact in no other way except the competition for a single rate-limiting nutrient that affects their growth rates. The issue of whether pure and simple competition under time-invariant conditions can give rise to chaotic behaviour has been unresolved for decades. The third author recently showed, for the first time, that chaos can theoretically occur in these systems by analysing the dynamics of a model where both competing species grow following the biomass-dependent Contois model while the yield coefficients associated with the two species are substrate-dependent. In this paper we show that chaotic behaviour can occur in a much simpler model of pure and simple competition. We examine the case where only one species grows following the Contois model with variable yield coefficient while the other species is allowed to grow following the simple Monod model with constant yield. We show that while the static behaviour of the proposed model is quite simple, the dynamic behaviour is complex and involves period doubling culminating in chaos. The proposed model could serve as a basis to re-examine the importance of Contois kinetics in predicting complex behaviour in microbial competition.

In this paper, optimal consumption and investment decisions are studied for an investor who has available a bank account and a stock whose price is a log normal diffusion. The bank pays at an interest rate r(t) for any deposit, and vice takes at a larger rate r′(t) for any loan. Optimal strategies are obtained via Hamilton-Jacobi-Bellman (HJB) equation which is derived from dynamic programming principle. For the specific HARA case, we get the optimal consumption and optimal investment explicitly, which coincides with the classical one under the condition r′(t) ≡ r(t)

We construct a new stochastic interest rate model with two stochastic factors, by introducing a stochastic long-run equilibrium level into the Vasicek interest rate model which follows another Ornstein–Uhlenbeck process. With the interest rate under the Black–Scholes model being assumed to follow the newly proposed model, a closed-form representation of European option prices is successfully presented, when the analytical characteristic function of the underlying log-price under a forward measure is derived. To assess the model performance, a preliminary empirical study is conducted using S&P 500 index and its options, with the Vasicek model and an alternative two-factor Vasicek model taken as benchmarks.

The paper is mainly concerned with the difference equation

where k and m are parameters, with k > 0. This equation arises from a method proposed for solving a cubic equation by iteration and represents a standardised form of the general problem. In using the above equation it is essential to know when the iteration process converges and this is discussed by means of the usual stability criterion. Critical values are obtained for the occurrence of solutions with period two and period three and the stability of these solutions is also examined. This was done by considering the changes as k increases, for a give value of m, which makes it effectively a one-parameter problem, and it turns out that the change with k can differ strongly from the usual behaviour for a one-parameter difference equation. For m = 2, for example it appears that the usual picture of stable 2-cycle solutions giving way to stable 4-cycle solutions is valid for smaller values of k but the situation is recersed for larger values of k where stable 4-cycle solutions precede stable 2-cycle solutions. Similar anomalies arise for the 3-cycle solutions.

We propose a Legendre–Laguerre spectral approximation to price the European and double barrier options in the time-fractional framework. By choosing an appropriate basis function, the spectral discretization is used for the approximation of the spatial derivatives of the time-fractional Black–Scholes equation. For the time discretization, we consider the popular $L1$ finite difference approximation, which converges with order $\mathcal {O}((\Delta \tau )^{2-\alpha })$ for functions which are twice continuously differentiable. However, when using the $L1$ scheme for problems with nonsmooth initial data, only the first-order accuracy in time is achieved. This low-order accuracy is also observed when solving the time-fractional Black–Scholes European and barrier option pricing problems for which the payoffs are all nonsmooth. To increase the temporal convergence rate, we therefore consider a Richardson extrapolation method, which when combined with the spectral approximation in space, exhibits higher order convergence such that high accuracies over the whole discretization grid are obtained. Compared with the traditional finite difference scheme, numerical examples clearly indicate that the spectral approximation converges exponentially over a small number of grid points. Also, as demonstrated, such high accuracies can be achieved in much fewer time steps using the extrapolation approach.

The motion of a two-dimensional bubble rising at a constant velocity U in an inclined tube of width H is considered. The bubble extends downwards without limit, and is bounded on the right by a wall of the tube, and on the left by a free surface. The same flow configuration describes also a jet emerging from a nozzle and falling down along an inclined wall. The acceleration of gravity g and the surface tension T are included in the free surface condition. The problem is characterized by the Froude number the angle β between the left wall and the horizontal, and the angle γ between the free surface and the right wall at the separation point. Numerical solutions are obtained via series truncation for all values of 0 < β < π. The results extend previous calculations of Vanden-Broeck [12–14] for β = π/2 and of Couët and Strumolo [3] for 0 < β < π/2. It is found that the behavior of the solutions depends on whether 0 < β 2π/3 or 2π/3 ≤ β < π. When T = 0, it is shown that there is a critical value F of Froude number for each 0 < β 2π/3 such that solutions with γ = 0, π/3 and π - β occur for F > Fc F = Fc and F < Fc respectively, and that all solutions are characterized by γ = 0 for 2π/3 ≤ β < π. When a small amount of surface tension T is included in the free surface condition, it is found that for each 0 < β < π there exists an infinite discrete set of values of F for which γ = π - β. A particular value F* of the Froude number for which T = 0 and γ = π - β is selected by taking the limit as T approaches zero. The numerical values of F* and the corresponding free surface profiles are found to be in good agreement with experimental data for bubbles rising in an inclined tube when 0 < β < π/2.