As a possible model for fluid turbulence, a Reiner–Rivlin-type equation is used to study Poiseuille–Couette flow of a viscous fluid in a rotating cylindrical pipe. The equations of motion are derived in cylindrical coordinates, and small-amplitude perturbations are considered in full generality, involving all three velocity components. A new matrix-based numerical technique is proposed for the linearized problem, from which the stability is determined using a generalized eigenvalue approach. New results are obtained in this cylindrical geometry, which confirm and generalize the predictions of previous recent studies. A possible mechanism for the transition to turbulent flow is discussed.

Safescape is a Western Australian company that has recently developed a device for improved safety in open-pit mines. Serious accidents can occur when large trucks veer off the roads running around the edge of the mine. The conventional technique to mitigate the risk is to pile waste rock to form a so-called bund on the edge of the road. This method is not fail safe though as vehicles can, and do, drive completely over the bund. In this paper, we describe a new device that consists of a row of filled polyethylene shell units which are linked together and sit on the road side of the rock bund. The vertical front face of the edge protector prevents out of control dump trucks from climbing over the bund and into the pit, so that they push against the barriers and heave the broken rock behind the bund. The models developed here suggest that the primary resistance to an impacting truck is provided by the large heaving force with the barrier simply facilitating this process. The theory indicates that the total resistance is independent of truck speed, meaning that simple barrier pushing experiments are sufficient to validate the analysis. The conclusions of the theory and field tests suggest that in a worst-case scenario involving the normal impact of a 500 tonne filled dump truck, the barriers and bund move a few metres before coming to rest.

Understanding how seasonal patterns change from year to year is important for the management of infectious disease epidemics. Here, we present a mathematical formalization of the application of complex demodulation, which has previously only been applied in an exploratory manner in the context of infectious diseases. This method extracts the changing amplitude and phase from seasonal data, allowing comparisons between the size and timing of yearly epidemics. We first validate the method using synthetic data that displays the key features of epidemic data. In particular, we analyse both annual and biennial synthetic data, and explore the effect of delayed epidemics on the extracted amplitude and phase. We then demonstrate the usefulness of complex demodulation using national notification data for influenza in Australia. This method clearly highlights the higher number of notifications and the early peak of the influenza pandemic in 2009. We also identify that epidemics that peaked later than usual generally followed larger epidemics and involved fewer overall notifications. Our analysis establishes a role for complex demodulation in the study of seasonal epidemiological events.

An initial-boundary value problem for a time-fractional diffusion equation is discretized in space, using continuous piecewise-linear finite elements on a domain with a re-entrant corner. Known error bounds for the case of a convex domain break down, because the associated Poisson equation is no longer $H^{2}$-regular. In particular, the method is no longer second-order accurate if quasi-uniform triangulations are used. We prove that a suitable local mesh refinement about the re-entrant corner restores second-order convergence. In this way, we generalize known results for the classical heat equation.

Pricing variance swaps have become a popular subject recently, and most research of this type come under Heston’s two-factor model. This paper is an extension of some recent research which used the dimension-reduction technique based on the Heston model. A new closed-form pricing formula focusing on a log-return variance swap is presented here, under the assumption that the underlying asset prices can be described by a mean-reverting Gaussian volatility model (Ornstein–Uhlenbeck process). Numerical tests in two respects using the Monte Carlo (MC) simulation are included. Moreover, we discuss a procedure of solving a quadratic differential equation with one variable. Our method can avoid the previously encountered limitations, but requires more time for calculation than other recent analytical discrete models.

We extend the evolution equation for weak nonlinear gravity–capillary waves by including fifth-order nonlinear terms. Stability properties of a uniform Stokes gravity–capillary wave train is studied using the evolution equation obtained here. The region of stability in the perturbed wave-number plane determined by the fifth-order evolution equation is compared with that determined by third- and fourth-order evolution equations. We find that if the wave number of longitudinal perturbations exceeds a certain critical value, a uniform gravity–capillary wave train becomes unstable. This critical value increases as the wave steepness increases.

We present a new iterative model order reduction method for large-scale linear time-invariant dynamical systems, based on a combined singular value decomposition–adaptive-order rational Arnoldi (SVD-AORA) approach. This method is an extension of the SVD-rational Krylov method. It is based on two-sided projections: the SVD side depends on the observability Gramian by the resolution of the Lyapunov equation, and the Krylov side is generated by the adaptive-order rational Arnoldi based on moment matching. The use of the SVD provides stability for the reduced system, and the use of the AORA method provides numerical efficiency and a relative lower computation complexity. The reduced model obtained is asymptotically stable and minimizes the error ($H_{2}$ and $H_{\infty }$) between the original and the reduced system. Two examples are given to study the performance of the proposed approach.