In this paper we investigate the population dynamics of a species with age structure in the case where the diffusion and death rates of the matured population are both age-dependent. We develop a new application of the age-structure technique in terms of an integral equation. For unbounded spatial domains, we study the existence of travelling waves, whilst in bounded domains, we investigate the existence of positive steady-state solutions and their stability.

In this study, we examine a steady two-dimensional slow flow past a rigid cylinder coated with a thin layer of immiscible fluid. The Reynolds number for the external bulk flow is assumed small and flow within the film is driven by the action of the bulk fluid’s tangential viscous stress acting at the interface. Using double asymptotic expansions based on the bulk fluid’s Reynolds number and the aspect ratio of the film thickness to the cylinder’s radius, we derive the leading- and first-order equations governing the steady-state film dynamics, and obtain analytical solutions, in terms of the film thickness, for the bulk flow. We solve the governing film equations, finding that solutions feature a drained region. We briefly discuss the influence of the Capillary number and fluid viscosities, and conclude by showing how the presence of the film affects the drag on the film-coated cylinder.

A new framework for modelling the evolution of the thermal bar system in a lake is presented. The model assumes that the thermal bar is located between two regions: the deeper region, where spring warming leads to overturning of the entire water column, and the near shore shallower region, where a stable surface layer is established. In this model the thermal bar moves out slightly more quickly than predicted by a simple thermal balance. Also, the horizontal extent of the thermal bar region increases as it moves out from the shore.

The existence and multiplicity of solutions to a quasilinear, elliptic partial differential equation with singular non-linearity is analysed. The partial differential equation is a recently derived variant of a canonical model used in the modelling of micro-electromechanical systems. It is observed that the bifurcation curve of solutions terminates at single dead-end point, beyond which no classical solutions exist. A necessary condition for the existence of solutions is developed, revealing that this dead-end point corresponds to a blow-up in the solution's gradient at a point internal to the domain. By employing a novel asymptotic analysis in terms of two small parameters, an accurate characterization of this dead-end point is obtained. An arc length parameterization of the solution curve can be employed to continue solutions beyond the dead-end point; however, all extra solutions are found to be multi-valued. This analysis therefore suggests that the dead-end is a bifurcation point associated with the onset of multi-valued solutions for the system.

It is shown that the set of conservation laws for the nonlinear system of equations describing plane steady potential barotropic flow of gas is given by the set of conservation laws for the linear Chaplygin system. All the conservation laws of zero order for the Chaplygin system are found. These include both known and new nonlinear conservation laws. It is found that the number of conservation laws of the first order is not more than three, assuming that the laws do not depend on the velocity potential and are not non-obvious ones. The components of these conservation laws are quadratic with respect to the stream function and its derivatives. All the Chaplygin functions are found, for which the Chaplygin system has three non-obvious conservation laws of the first order that are independent of velocity potential. All such non-obvious first-order conservation laws are found.

We prove the existence of a one-parameter family of solutions of the porous medium equation in which the interface is a half line whose end point advances at a constant speed. Then we prove the stability of the solutions under a suitable class of perturbations. We discuss the relevance of these solutions to gravity-driven flows of thin films, and show that some solutions develop a very thin triangular plateau in the direction of propagation and that the angle of the plateau and its thickness are decreasing functions of the speed.

In this paper, a non-linear mathematical model is proposed and analysed to study the role of technology in combating social crimes in a dynamic population by considering immigration and emigration rates of susceptible population and criminals. The problem is modelled by considering five interacting variables, namely the density of susceptible population, the density of criminals, the density of removed (isolated) criminals, the density of crime burden and the level of technology used to control crime. The proposed model is analysed by using the stability theory of differential equation and simulation. The model analysis shows that the crime burden decreases considerably as the level of technology increases. It is noted that the crime in a society can be controlled almost completely if criminals from the general population are removed by intensive use of technology.

To study electrostatic actuation, researchers commonly use a setup proposed by G. I. Taylor in [Proc. R. Soc. Lond. Ser. A, 306 (1968), pp. 423–434]. It consists of soap film held at a distance h above a rigid plate so that when a voltage difference is applied between the two components, the top film deflects towards the bottom plate. The most striking feature of this system is when the voltage difference exceeds a critical value V*, the electrostatic forces dominate the surface forces and the soap film gets ‘pulled-into’ or collapses onto the bottom plate. This so-called ‘pull-in’ instability is a ubiquitous feature of electrostatic actuation and as a result, has been the subject of many studies. Recently, Siddique et al. [J. Electrostatics, 69 (2011), pp. 1–6] measured the value of V* as a function of the separation distance and found that the standard prediction breaks down as h increases. Here, we continue the work done in [N. D. Brubaker and J. A. Pelesko, European J. Appl. Math., 22 (2011), pp. 455–470] by investigating the cause of this discrepancy. Specifically, we model the effect of gravity on the generalized version of Taylor's model and study whether it provides the proper correction to the predicted value of V*. In doing so, we derive two nonlinear eigenvalue value problems and investigate their solutions sets.

In this paper, we analytically investigate multi-component Cahn–Hilliard and Allen–Cahn systems which are coupled with elasticity and uni-directional damage processes. The free energy of the system is of the form ∫Ω½Γ∇c : ∇c + ½|∇z|2+Wch(c)+Wel(e,c,z)dx with a polynomial or logarithmic chemical energy density Wch, an inhomogeneous elastic energy density Wel and a quadratic structure of the gradient of damage variable z. For the corresponding elastic Cahn–Hilliard and Allen–Cahn systems coupled with uni-directional damage processes, we present an appropriate notion of weak solutions and prove existence results based on certain regularization methods and a higher integrability result for strain e.

In this paper we consider a two-phase model describing the growth of avascular solid tumors when taking into account the effects of cell-to-cell adhesion and taxis due to nutrient. The tumor is surrounded by healthy tissue which is the source of nutrient for tumor cells. In a three-dimensional context, we prove that the mathematical formulation corresponds to a well-posed problem, and find radially symmetric steady-state solutions of the problem. They appear in the regime where the rate of cell apoptosis to cell proliferation is less than the far field nutrient concentration. Furthermore, we study the stability properties of those radially symmetric equilibria and find, depending on the biophysical parameters involved in the problem, both stable and unstable regimes for tumor growth.

In this paper, a bilayer model is derived to simulate the evolution of a thin film flow over water. This model is derived from the incompressible Navier–Stokes equations together with suitable boundary conditions including friction and capillary effects. The derivation is based on the different properties of the fluids; thus, we perform a multiscale analysis in space and time, and a different asymptotic analysis to derive a system coupling two different models: the Reynolds lubrication equation for the upper layer and the shallow water model for the lower one. We prove that the model verifies a dissipative entropy inequality up to a second-order term. Moreover, we propose a correction of the model – by taking into account the second-order extension for the pressure – that admits an exact dissipative entropy inequality. Two numerical tests are presented. In the first test, we compare the numerical results with the viscous bilayer shallow water model proposed in Narbona-Reina et al. (Comput. Model. Eng. Sci., 2009, Vol. 43, pp. 27–71). In the second test, the objective is to show some of the characteristic situations that can be studied with the proposed model. We simulate a problem of pollutant dispersion near the coast. For this test, the influence of the friction coefficient on the coastal area affected by the pollutant is studied.

We prove an upper bound for the convergence rate of the homogenization limit ε → 0 for a linear transmission problem for a advection–diffusion(–reaction) system posed in areas with low and high diffusivity, where ε is a suitable scale parameter. In this way we rigorously justify the formal homogenization asymptotics obtained in [37] (van Noorden, T. and Muntean, A. (2011) Homogenization of a locally-periodic medium with areas of low and high diffusivity. Eur. J. Appl. Math. 22, 493–516). We do this by providing a corrector estimate. The main ingredients for the proof of the correctors include integral estimates for rapidly oscillating functions with prescribed average, properties of the macroscopic reconstruction operators, energy bounds, and extra two-scale regularity estimates. The whole procedure essentially relies on a good understanding of the analysis of the limit two-scale problem.

We investigate the motion of two immiscible fluids in a porous medium described by a two-phase flow system. In the capillary pressure relation, we include static and dynamic hysteresis. The model is well established in the context of the Richards equation, which is obtained by assuming a constant pressure for one of the two phases. We derive an existence result for this hysteresis two-phase model for non-degenerate permeability and capillary pressure curves. A discretization scheme is introduced and numerical results for fingering experiments are obtained. The main analytical tool is a compactness result for two variables that are coupled by a hysteresis relation.

We consider an anisotropic diffusion equation of the form ut = ∇∇(D(x)u) in two dimensions, which arises in various applications, including the modelling of wolf movement along seismic lines and the invasive spread of certain brain tumours along white matter neural fibre tracts. We consider a degenerate case, where the diffusion tensor D(x) has a zero-eigenvalue for certain values of x. Based on a regularisation procedure and various pointwise and integral a priori estimates, we establish the global existence of very weak solutions to the degenerate limit problem. Moreover, we show that in the large time limit these solutions approach profiles that exhibit a Dirac-type mass concentration phenomenon on the boundary of the region in which diffusion is degenerate, which is quite surprising for a linear diffusion equation. The results are illustrated by numerical examples.

In a one-dimensional domain, the stability of localized spike patterns is analysed for two closely related singularly perturbed reaction–diffusion (RD) systems with Brusselator kinetics. For the first system, where there is no influx of the inhibitor on the domain boundary, asymptotic analysis is used to derive a non-local eigenvalue problem (NLEP), whose spectrum determines the linear stability of a multi-spike steady-state solution. Similar to previous NLEP stability analyses of spike patterns for other RD systems, such as the Gierer–Meinhardt and Gray–Scott models, a multi-spike steady-state solution can become unstable to either a competition or an oscillatory instability depending on the parameter regime. An explicit result for the threshold value for the initiation of a competition instability, which triggers the annihilation of spikes in a multi-spike pattern, is derived. Alternatively, in the parameter regime when a Hopf bifurcation occurs, it is shown from a numerical study of the NLEP that an asynchronous, rather than synchronous, oscillatory instability of the spike amplitudes can be the dominant instability. The existence of robust asynchronous temporal oscillations of the spike amplitudes has not been predicted from NLEP stability studies of other RD systems. For the second system, where there is an influx of inhibitor from the domain boundaries, an NLEP stability analysis of a quasi-steady-state two-spike pattern reveals the possibility of dynamic bifurcations leading to either a competition or an oscillatory instability of the spike amplitudes depending on the parameter regime. It is shown that the novel asynchronous oscillatory instability mode can again be the dominant instability. For both Brusselator systems, the detailed stability results from NLEP theory are confirmed by rather extensive numerical computations of the full partial differential equations system.

The subcritical flow of a stream over a bottom obstruction or depression is considered with particular interest in obtaining solutions with no downstream waves. In the linearised problem this can always be achieved by superposition of multiple obstructions, but it is not clear whether this is possible in a full nonlinear problem. Solutions computed here indicate that there is an effective nonlinear superposition principle at work as no special shape modifications were required to obtain wave-cancelling solutions. Waveless solutions corresponding to one or more trapped waves are computed at a range of different Froude numbers and are shown to provide a rather elaborate mosaic of solution curves in parameter space when both negative and positive obstruction heights are included.

In this paper, we consider rational or multi-cut solutions to Hele-Shaw flow driven by suction. A lower bound of the distance between the moving boundary and the suction point is given and this lower bound enables us to find necessary conditions for rational or multi-cut solutions to blow up due to the boundary touching the sink.

In this paper, we study singularly perturbed Filippov systems. More specifically, our main question is to know how the dynamics of Filippov systems is affected by singular perturbations. We extend the Fenichel theory developed in Fenichel (J. Differ. Equ., 1979, Vol. 31, pp. 53–98) to these systems. In addition, the study of non-smooth constrained systems is considered.

A free boundary problem, which comes from the model of the perpetual American call options with utility functions in financial market, is investigated. It is a degenerative parabolic free boundary problem and is studied by the line method. The existence, regularity and uniqueness of the solution as well as some properties of the free boundary are established.

We use the heat equation as an illustrative example to show that the unified method introduced by one of the authors can be employed for constructing analytical solutions for linear evolution partial differential equations in one spatial dimension involving non-separable boundary conditions as well as non-local constraints. Furthermore, we show that for the particular case in which the boundary conditions become separable, the unified method provides an easier way for constructing the relevant classical spectral representations avoiding the classical spectral analysis approach. We note that the unified method always yields integral expressions which, in contrast to the series or integral expressions obtained by the standard transform methods, are uniformly convergent at the boundary. Thus, even for the cases that the standard transform methods can be implemented, the unified method provides alternative solution expressions which have advantages for both numerical and asymptotic considerations. The former advantage is illustrated by providing the numerical evaluation of typical boundary value problems.