While the quasilinear isothermal Euler equations are an excellent model for gas pipeline flow, the operation of the pipeline flow with high pressure and small Mach numbers allows us to obtain approximate solutions by a simpler semilinear model. We provide a derivation of the semilinear model that shows that the semilinear model is valid for sufficiently low Mach numbers and sufficiently high pressures. We prove an existence result for continuous solutions of the semilinear model that takes into account lower and upper bounds for the pressure and an upper bound for the magnitude of the Mach number of the gas flow. These state constraints are important both in the operation of gas pipelines and to guarantee that the solution remains in the set where the model is physically valid. We show the constrained exact boundary controllability of the system with the same pressure and Mach number constraints.

In the current work, we study a stochastic parabolic problem. The presented problem is motivated by the study of an idealised electrically actuated MEMS (Micro-Electro-Mechanical System) device in the case of random fluctuations of the potential difference, a parameter that actually controls the operation of MEMS device. We first present the construction of the mathematical model, and then, we deduce some local existence results. Next for some particular versions of the model, relevant to various boundary conditions, we derive quenching results as well as estimations of the probability for such singularity to occur. Additional numerical study of the problem in one dimension follows, which also allows the further investigation the problem with respect to its quenching behaviour.

Biofilms, and collections of embedded microbial communities, present structural heterogeneities with functional consequences for important processes, such as transport. The origin of such structures has been unclear. Here, we propose that they can arise as a consequence of diffusive transport limitation. To illustrate, a model allowing internal heterogeneity is developed. Linear analysis is applied to a simplified version of the model suggesting that heterogeneity forms on (or below) the active layer length, a length scale that may not be suitable for homogenization, with non-trivial implications for system scale properties such as reduction in system-wide diffusive transport efficiency. Numerics suggest that the simplified model provides useful insight into behaviour of the full model. We then show examples based on microcolony formation in host domains and argue that internal heterogeneity can be related to community function.

In this contribution, we present a modelling and simulation framework for parametrised lithium-ion battery cells. We first derive a continuum model for a rather general intercalation battery cell on the basis of non-equilibrium thermodynamics. In order to efficiently evaluate the resulting parameterised non-linear system of partial differential equations, the reduced basis method is employed. The reduced basis method is a model order reduction technique on the basis of an incremental hierarchical approximate proper orthogonal decomposition approach and empirical operator interpolation. The modelling framework is particularly well suited to investigate and quantify degradation effects of battery cells. Several numerical experiments are given to demonstrate the scope and efficiency of the modelling framework.

Droplet profiles under the influence of van der Waals forces

Euro. Jnl. of Applied Maths.12, 367–393

ELI MINKOV and AMY NOVICK-COHEN

The intrusion of a constant density fluid at the interface of a two-layer fluid is considered. Numerical solutions are computed for a model of a steady intrusion resulting from flow down a bank and across a broad lake or reservoir. The incoming fluid is homogeneous and spreads across the lake at its level of neutral buoyancy. Solutions are obtained for a range of different inflow angles, flow rate and density differences. Except in extreme cases, the nature of the solution is predicted quite well by linear theory, with the wavelength at any Froude number given by a dispersion relation and wave steepness determined largely by entry angle. However, some extreme solutions with rounded meandering flows and non-unique solutions in the parameter space are also obtained.

We investigate the Brusselator reaction–diffusion equations with periodic boundary conditions. We consider the range of values of the parameters used by Kuramoto in his study of chaotic concentration waves. We determine numerically the bifurcation diagram of the long-time travelling and standing wave solutions using a highly accurate Fourier pseudo-spectral method. For moderate values of the bifurcation parameter, we have found a sequence of instabilities leading either to periodic and quasiperiodic standing waves, or to chaotic regimes. However, for large values of the control parameter, we have found only uniform time-periodic solutions or time-periodic travelling wave solutions. Our numerical study motivates a new asymptotic analysis of the Brusselator equations for large values of the control parameter and small diffusion coefficients. This analysis explains the numerical predictions. The chaotic regime is limited to moderate values of the control parameter and periodic solutions are the only solutions for large values of the control parameter. We identify the stabilizing mechanism as the relaxation oscillations which appear when the control parameter is large. Our asymptotic result on the stability of periodic solutions is then generalized to a class of two-variable reaction-diffusion equations.

Amortization of a shock in an electro-rheological shock absorber is carried out in the motion of a piston in an electrorheological fluid. The drag force acting on the piston is regulated by varying the voltage applied to electrodes. A model of an electrorheological shock absorber is constructed. A problem on shock absorber reduces to the solution of a coupled problem for motion equation of the piston and non-linear equations of fluid flow in an unknown domain that varies with the time. A method of semi-discretization for approximate solution of the coupled problem is considered. Results on the existence and on the uniqueness of the solution of the coupled problem are obtained. Convergence of approximate solutions to the exact solution is proved. Numerical simulation of the operation of the shock absorber is performed.

Uniqueness and solvability theorems are proved for the two-dimensional Neumann–Kelvin problem in the case when a body is totally submerged in a subcritical stream of finite depth fluid. A version of source method is developed to find a solution. The Green's identity coupling the solution with a solution of the problem with opposite stream direction is used to prove that the solution is unique.

We study the propagation of a crack in critical equilibrium for a brittle material in a Mode III field. The energy variations for small virtual extensions of the crack are handled in a novel way: the amount of energy released is written as a functional over a family of univalent functions on the upper half plane. Classical techniques developed in connection to the Bieberbach Conjecture are used to quantify the energy-shape relationship. By means of a special family of trial paths generated by the so-called Löwner equation we impose a stability condition on the field which derives in a local crack propagation criterion. We called this the anti-symmetry principle, being closely related to the well known symmetry principle for the in-plane fields.

Delayed production can substantially alter the qualitative behaviour of feedback systems. Motivated by stochastic mechanisms in gene expression, we consider a protein molecule which is produced in randomly timed bursts, requires an exponentially distributed time to activate and then partakes in positive regulation of its burst frequency. Asymptotically analysing the underlying master equation in the large-delay regime, we provide tractable approximations to time-dependent probability distributions of molecular copy numbers. Importantly, the presented analysis demonstrates that positive feedback systems with large production delays can constitute a stable toggle switch even if they operate with low copy numbers of active molecules.

The quasi-steady power-law Stokes flow of a mixture of incompressible fluids with shear-dependent viscosity is studied. The fluids are immiscible and have constant densities. Existence results are presented for both the no-slip and the no-stick boundary value conditions. Use is made of Schauder's fixed-point theorem, compactness arguments, and DiPerna–Lions renormalized solutions.

An infinite regular two-dimensional network is composed of unit resistances joining adjacent nodes. What is the resistance of the whole network measured between two different nodes? Three cases are considered, in which the meshes of the network are square, hexagonal or triangular. Some extensions of the problem are also treated. Asymptotic approximations are given for large distances.

This paper considers the role of costless decisions relating to the extraction of a non-renewable resource in the presence of uncertainty. We begin by deriving a size scale of the extractable resource, above which the solution to the valuation and optimal control strategy can be described by analytic solutions; we produce solutions for a general form of operating cost function. Below this critical resource size level the valuation and optimal control strategy must be solved by numerical means; we present a robust numerical algorithm that can solve such a class of problem. We also allow for the embedding of an irreversible investment decision (abandonment) into the optimisation. Finally, we conduct experimentation for each of these two approaches (analytical and numerical), and show how they are consistent with one another when used appropriately. The extensions of this paper's techniques to renewable resources are explored.

We consider a processor-sharing storage allocation model, which has m primary holding spaces and infinitely many secondary ones, and a single processor servicing the stored items. All of the spaces are numbered and ordered. An arriving customer takes the lowest available space. Dynamic storage allocation and the fragmentation of computer memory are well-known applications of this model. We define the traffic intensity ρ to be λ/μ, where λ is the customers' arrival rate and μ is the service rate of the processor. We study the joint probability distribution of the numbers of occupied primary and secondary spaces. We study the problem in two asymptotic limits: (1) m → ∞ with a fixed ρ < 1, and (2) ρ ↑ 1, m → ∞ with m(1-ρ) = O(1). The asymptotics yield insight into how many secondary spaces tend to be needed, and into the sample paths leading to the occupation of the two types of spaces. We show that the asymptotics lead to accurate numerical approximations.

We consider periodic solutions which bifurcate from equilibria in simple population models which incorporate a state-dependent time delay of the discrete kind. The delay is a function of the current size of the population. Solutions near equilibria are constructed using perturbation methods to determine the sub/supercriticality of the bifurcation and hence their stability. The stability of the bifurcating solutions depends on the qualitative form of the delay function. This is in contrast to the stability of an equilibrium, which is determined purely by the actual value of this function at the equilibrium.

The Stefan problem is coupled with a spatially inhomogeneous and anisotropic Gibbs–Thomson condition at the phase boundary. We show the long-time existence of weak solutions for the non-degenerate Stefan problem with a spatially inhomogeneous and anisotropic Gibbs–Thomson law and a conditional existence result for the corresponding degenerate Stefan problem. To this end, approximate solutions are constructed by means of variational problems for energy functionals with spatially inhomogeneous and anisotropic interfacial energy. By passing to the limit, we establish solutions of the Stefan problem with a spatially inhomogeneous and anisotropic Gibbs–Thomson law in a weak generalised BV-formulation.

A free boundary problem arising in a model for inviscid, incompressible shallow water entry at small deadrise angles is derived and analysed. The relationship between this novel free boundary problem and the well-known viscous squeeze film problem is described. An inverse method is used to construct explicit solutions for certain body profiles and to find criteria under which the splash sheet can ‘split’. A variational inequality formulation, conservation of certain generalized moments and the Schwarz function formulation are introduced.

We discuss the fourth-order thin film equation with a stable second-order diffusion term, in the context of a standard free-boundary problem with zero height, zero contact angle and zero-flux conditions imposed at an interface. For the first critical exponent where N ≥ 1 is the space dimension, there are continuous sets (branches) of source-type very singular self-similar solutions of the form For p ≠ p0, the set of very singular self-similar solutions is shown to be finite and consists of a countable family of branches (in the parameter p) of similarity profiles that originate at a sequence of critical exponents {pl, l ≥ 0}. At p = pl, these branches appear via a non-linear bifurcation mechanism from a countable set of second kind similarity solutions of the pure thin film equation Such solutions are detected by a combination of linear and non-linear ‘Hermitian spectral theory’, which allows the application of an analytical n-branching approach. In order to connect with the Cauchy problem in Part I, we identify the cauchy problem solutions as suitable ‘limits’ of the free-boundary problem solutions.

In a recent paper, a new three-parameter class of Abel type equations, so-called AIR, all of whose members can be mapped into Riccati equations, is shown. Most of the Abel equations with solution presented in the literature belong to the AIR class. Three canonical forms were shown to generate this class, according to the roots of a cubic. In this paper, a connection between those canonical forms and the differential equations for the hypergeometric functions $_2{\rm F}_1$, $_1{\rm F}_1$ and $_0{\rm F}_1$ is unveiled. This connection provides a closed form $_p{\rm F}_q$ solution for all Abel equations of the AIR class.