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.

This paper deals with the motion of a fluid in a closed loop under the effect of natural convection and a given external heat flux. More precisely, we show that the stationary solutions of a system describing the intermediate asymptotics of the previous problem are structurally linearly unstable.

Although in recent years bone piezoelectricity has been normally neglected, lately a new interest has appeared to show the importance of bone piezoelectricity in wet bone's complex response to loading. Here we numerically study a problem, including a strain-adaptive bone remodelling and the piezoelectricity. Its variational formulation leads to a coupled system composed of two linear variational equations for displacements and electric potential, and a parabolic variational inequality for the apparent density. Fully discrete approximations are now introduced by using the finite element method to approximate spatial variable and the explicit Euler scheme to discretise time derivatives. Some a priori error estimates are proved and the linear convergence of the algorithm is deduced under additional regularity conditions. Finally, some one- and two-dimensional numerical simulations are described to show the accuracy of the proposed algorithm and the behaviour of the solution.

The inspection game has been widely used as a model for crime deterrence. In this article, we shall present a broad review of the literature and an original result. The original result concerns the interaction between a single, long-run inspector and a population of short-run, boundedly rational would-be criminals. This approach gives rise to an optimal control problem that can be solved explicitly. We show that a forward-looking inspector obtains the same payoff he would obtain if he could commit to a probability of inspection. The consequences of this result for the economics of crime deterrence are also explored.

The linear Schrödinger equation with piecewise constant potential in one spatial dimension is a well-studied textbook problem. It is one of only a few solvable models in quantum mechanics and shares many qualitative features with physically important models. In examples such as ‘particle in a box’ and tunnelling, attention is restricted to the time-independent Schrödinger equation. This paper combines the unified transform method and recent insights for interface problems to present fully explicit solutions for the time-dependent problem.

The paper deals with the probability density function (PDF) of the concentration of a scalar within a turbulent flow. Following some comments about the overall structure of the PDF, and its approach to a limit at large times, attention focusses on the so-called small scale mixing term in the evolution equation for the PDF. This represents the effect of molecular diffusion in reducing concentration uctuations, eventually to zero. Arguments are presented which suggest that this quantity could, in certain circumstances, depend inversely upon the PDF, and a particular example of this leads to a new closure hypothesis. Consequences of this, especially similarity solutions, are explored for the case when the concentration field is statistically homogeneous.

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.

We consider the behaviour of a premixed flame anchored on a flat burner. For Lewis numbers L < L* < 1, stationary spatially periodic solutions corresponding to stationary cellular flames bifurcate from the basic solution which corresponds to a steady planar flame. We study the existence and stability of two-dimensional patterns which correspond to certain imposed symmetries by considering the evolution of N pairs of wave vectors, each of which is separated from the next by angle π/N. In the neighbourhood of the critical Lewis number L*, we derive evolution equations for the amplitudes corresponding to N = 2, which corresponds to square patterns, and N = 3, which corresponds to triangular or hexagonal patterns. We determine existence and stability results in terms of m∈(0, 1), the flow rate of the fuel, and K > 2/e, the scaled heat loss to the burner. Square patterns exist for L < L* and are stable for values of m and K above a stability boundary in the m−K plane, which has a maximum at K = K* ∼ 4.77, so that for K > K* square patterns are stable for all m. The stability of the square patterns does not vary with L. Hexagonal patterns exist for L < LH, where 1 > LH > L*. The size of the stability region increases with decreasing L < L*. For a range of values of L there is bistability, that is, for given parameter values rolls and hexagons are simultaneously stable, each with its own domain of attraction.

We present the first macroscopical model for charge transport in compound semiconductors to make use of analytic ellipsoidal approximations for the energy dispersion relationships in the neighbours of the lowest minima of the conduction bands. The model considers the main scattering mechanisms charges undergo in polar semiconductors, that is the acoustic, polar optical, intervalley non-polar optical phonon interactions and the ionized impurity scattering. Simulations are shown for the cases of bulk 4H and 6H-SiC.

This work considers the problem of binary classification: given training data x1, . . ., xn from a certain population, together with associated labels y1,. . ., yn ∈ {0,1}, determine the best label for an element x not among the training data. More specifically, this work considers a variant of the regularized empirical risk functional which is defined intrinsically to the observed data and does not depend on the underlying population. Tools from modern analysis are used to obtain a concise proof of asymptotic consistency as regularization parameters are taken to zero at rates related to the size of the sample. These analytical tools give a new framework for understanding overfitting and underfitting, and rigorously connect the notion of overfitting with a loss of compactness.

Linear stability theory is developed for an activator–inhibitor model where fractional derivative operators of generally different exponents act both on diffusion and reaction terms. It is shown that in the short wave limit the growth rate is a power law of the wave number with decoupled time scales for distinct anomaly exponents of the different species. With equal anomaly exponents an exact formula for the anomalous critical value of reactants diffusion coefficients' ratio is obtained.

We analyze a model for non-isothermal superconductivity, derived independently by G. Maugin, and K. Miya and S. A. Zhou. The model is described by a parabolic system based on the Time-Dependent Ginzburg–Landau (TDGL) equation, the Maxwell equations, and an energy equation such that the Clausius–Duhem inequality holds. The principal unknown fields are the complex valued Ginzburg–Landau order parameter $\psi$, the magnetic vector potential $A$, and the temperature $T$. A significant feature for this model is that it accounts for the interchange of thermal and electro-magnetic energies through Joule heating. The sensitivity of superconducting materials to thermal variations is a major obstacle to their applications. In practice one sees that time varying currents and magnetic fields generate thermal energy, producing ‘hot spots’ and increasing the temperature. This can result in suppressing superconductivity, i.e. $\psi$ tending to $0$. Our principal result is that we exhibit this phenomenon. We prove that if the electro-magnetic energy of the superconductor is sufficiently large at time $t\,{=}\,0$ then $T$ will rise beyond the critical temperature, $T_c$, and $\psi\to 0$ as $t\to\infty$. Earlier analytic work investigating electro-dynamic effects in superconductivity centered on an isothermal model consisting of the TDGL equation and Maxwell equations. This setting takes into account just electro-magnetic effects and produces markedly different evolutions.

We study the effect of stochastic volatility on option prices. In the fast mean-reversion model for stochastic volatility of [5], we show that there is a full asymptotic expansion for the option price, centered at the Black-Scholes price. We show how to callibrate the first two terms in the expansion with the implied volatility surface. We show, however, that this price does not converge in a strong sense to Black-Scholes as the mean-reversion rate increases.

We present an unsteady Hele–Shaw model of the fluid–fluid displacements that take place during primary cementing of an oil well, focusing on the case where one Herschel–Bulkley fluid displaces another along a long uniform section of the annulus. Such unsteady models consist of an advection equation for a fluid concentration field coupled to a third-order non-linear PDE (Partial differential equation) for the stream function, with a free boundary at the boundary of regions of stagnant fluid. These models, although complex, are necessary for the study of interfacial instability and the effects of flow pulsation, and remain considerably simpler and more efficient than computationally solving three-dimensional Navier–Stokes type models. Using methods from gradient flows, we demonstrate that our unsteady evolution equation for the stream function has a unique solution. The solution is continuous with respect to variations in the model physical data and will decay exponentially to a steady-state distribution if the data do not change with time. In the event that density differences between the fluids are small and that the fluids have a yield stress, then if the flow rate is decreased suddenly to zero, the stream function (hence velocity) decays to zero in a finite time. We verify these decay properties, using a numerical solution. We then use the numerical solution to study the effects of pulsating the flow rate on a typical displacement.

We study the effect of population mobility on the transmission dynamics of infectious diseases by considering a susceptible-exposed-infectious-recovered (SEIR) epidemic model with graph Laplacian diffusion, that is, on a weighted network. First, we establish the existence and uniqueness of solutions to the SEIR model defined on a weighed graph. Then by constructing Liapunov functions, we show that the disease-free equilibrium is globally asymptotically stable if the basic reproduction number is less than unity and the endemic equilibrium is globally asymptotically stable if the basic reproduction number is greater than unity. Finally, we apply our generalized weighed graph to Watts–Strogatz network and carry out numerical simulations, which demonstrate that degrees of nodes determine peak numbers of the infectious population as well as the time to reach these peaks. It also indicates that the network has an impact on the transient dynamical behaviour of the epidemic transmission.

In this paper we examine the movement of hard contact lenses on the eye. In so doing, we take into account hydrodynamic forces underneath the lens, as well as surface tension forces at the lens periphery. This involves solving for the free surface of the tear film away from the lens in order to determine the magnitudes of the pressure and surface tension forces on the lens. The analysis, which assumes quasi-steady motion, is carried out in both two and three dimensions.

Liquid crystal elastomers present features not found in ordinary elastic materials, such as semi-soft elasticity and the related stripe domain phenomenon. In this paper, the two-dimensional Bladon–Terentjev–Warner model and the one-constant Oseen–Frank energy expression are combined to study the liquid crystal elastomer. We also impose two material constraints, the incompressibility of the elastomer and the unit director norm of the liquid crystal. We prove existence of minimiser of the energy for the proposed model. Next we formulate the discrete model, and also prove that it possesses a minimiser of the energy. The inf-sup values of the discrete linearised system are then related to the smallest singular values of certain matrices. Next the existence and uniqueness of the Lagrange multipliers associated with the two material constraints are proved under the assumption that the inf-sup conditions hold. Finally numerical simulations of the clamped-pulling experiment are presented for elastomer samples with aspect ratio 1 or 3. The semi-soft elasticity is successfully recovered in both cases. The stripe domain phenomenon, however, is not observed, which might be due to the relative coarse mesh employed in the numerical experiment. Possible improvements are discussed that might lead to the recovery of the stripe domain phenomenon.

An analysis is undertaken of the formation and stability of localised patterns in a 1D Schanckenberg model, with source terms in both the activator and inhibitor fields. The aim is to illustrate the connection between semi-strong asymptotic analysis and the theory of localised pattern formation within a pinning region created by a subcritical Turing bifurcation. A two-parameter bifurcation diagram of homogeneous, periodic and localised patterns is obtained numerically. A natural asymptotic scaling for semi-strong interaction theory is found where an activator source term \[a = O(\varepsilon )\] and the inhibitor source \[b = O({\varepsilon ^2})\], with ε2 being the diffusion ratio. The theory predicts a fold of spike solutions leading to onset of localised patterns upon increase of b from zero. Non-local eigenvalue arguments show that both branches emanating from the fold are unstable, with the higher intensity branch becoming stable through a Hopf bifurcation as b increases beyond the \[O(\varepsilon )\] regime. All analytical results are found to agree with numerics. In particular, the asymptotic expression for the fold is found to be accurate beyond its region of validity, and its extension into the pinning region is found to form the low b boundary of the so-called homoclinic snaking region. Further numerical results point to both sub and supercritical Hopf bifurcation and novel spikeinsertion dynamics.

Semi-supervised and unsupervised machine learning methods often rely on graphs to model data, prompting research on how theoretical properties of operators on graphs are leveraged in learning problems. While most of the existing literature focuses on undirected graphs, directed graphs are very important in practice, giving models for physical, biological or transportation networks, among many other applications. In this paper, we propose a new framework for rigorously studying continuum limits of learning algorithms on directed graphs. We use the new framework to study the PageRank algorithm and show how it can be interpreted as a numerical scheme on a directed graph involving a type of normalised graph Laplacian. We show that the corresponding continuum limit problem, which is taken as the number of webpages grows to infinity, is a second-order, possibly degenerate, elliptic equation that contains reaction, diffusion and advection terms. We prove that the numerical scheme is consistent and stable and compute explicit rates of convergence of the discrete solution to the solution of the continuum limit partial differential equation. We give applications to proving stability and asymptotic regularity of the PageRank vector. Finally, we illustrate our results with numerical experiments and explore an application to data depth.

We present an adaptation of the Monge–Ampère (MA) lattice basis reduction scheme to the MA equation with second boundary value condition, provided the target is a convex set. This yields a fast adaptive method to numerically solve the optimal transport (OT) problem between two absolutely continuous measures, the second of which has convex support. The proposed numerical method actually captures a specific Brenier solution which is minimal in some sense. We prove the convergence of the method as the grid step size vanishes and show with numerical experiments that it is able to reproduce subtle properties of the OT problem.