This special issue of the European Journal of Applied Mathematics features research articles on networks, one of the most exciting areas of applied mathematics.

In this paper, we study a quasi-static frictional contact problem for a viscoelastic body with damage effect inside the body as well as normal compliance condition and multi-valued friction law on the contact boundary. The considered friction law generalizes Coulomb friction condition into multi-valued setting. The variational–hemi-variational formulation of the problem is derived and arguments of fixed point theory and surjectivity results for pseudo-monotone operators are applied, in order to prove the existence and uniqueness of solution.

In this paper, we show that there exists a critical number that stabilises the Reissner–Mindlin–Timoshenko system with frictional dissipation acting only on the equation for the transverse displacement. We identify that the Reissner–Mindlin–Timoshenko system has two speed characteristics v12 := K/ρ1 and v22 := D/ρ2 and we show that the system is exponentially stable if only if

\begin{equation*} v_{1}^{2}=v_{2}^{2}. \end{equation*}

In the general case, we prove that the system is polynomially stable with optimal decay rate. Numerical experiments using finite differences are given to confirm our analytical results. Our numerical results are qualitatively in agreement with the corresponding results from dynamical in infinite dimensional.

This special issue of the European journal of applied mathematics features research articles that involve the application of mathematical methodologies to the modelling of a broad range of problems related to crime and security. Some specific topics in this issue include recent developments in mathematical models of residential burglary, a dynamical model for the spatial spread of riots initiated by some triggering event, the analysis and development of game-theoretic models of crime and conflict, the study of statistically based models of insurgent activity and terrorism using real-world data sets, models for the optimal strategy of police deployment under realistic constraints, and a model of cyber crime as related to the study of spiking behaviour in social network cyberspace communications. Overall, the mathematical and computational methodologies employed in these studies are as diverse as the specific applications themselves and the scales (spatial or otherwise) to which they are applied. These methodologies range from statistical and stochastic methods based on maximum likelihood methods, Bayesian equilibria, regression analysis, self-excited Hawkes point processes, agent-based random walk models on networks, to more traditional applied mathematical methods such as dynamical systems and stability theory, the theory of Nash equilibria, rigorous methods in partial differential equations and travelling wave theory, and asymptotic methods that exploit disparate space and time scales.

In this paper, we consider the penalty method to solve the unilateral contact with friction between an electro-elastic body and a conductive foundation. Mathematical properties, such as the existence of a solution to the penalty problem and its convergence to the solution of the original problem, are reported. Then, we present a finite elements approximation for the penalised problem and prove its convergence. Finally, we propose an iterative method to solve the resulting finite element system and establish its convergence.

A spike solution is constructed on the infinite line for a sub-diffusive version of the Gierer–Meinhardt reaction – diffusion model. A non-local eigenvalue problem governs the spike's stability and is explicitly solvable for a certain choice of the kinetic parameters. Its solution generalises former results for the Gierer–Meinhardt model with regular diffusion, and the normal and anomalous systems' properties are juxtaposed. It is shown that a Hopf bifurcation occurs in the sub-diffusive system for larger values of the time parameter τo as compared to the normal counterpart, rendering the anomalous system more stable. Asymptotic solutions are obtained near important values of the diffusion anomaly index γ and collectively shown to be valid over most of the applicable range 0 < γ < 1. A bifurcation delay scenario is described for the sub-diffusive system, and the WKB exponent is computed analytically.

We examine the game theoretic properties of a model of crime first introduced by Short et al. (2010 Phys. Rev. E82, 066114) as the SBD Adversarial Game. We identify the rationalizable strategies and one-shot equilibria under multiple equilibrium refinements. We further show that SBD's main result about the effectiveness of defecting-punishers (“Informants”) in driving the system to evolve to the cooperative equilibrium under an imitation dynamic generalizes to a best response dynamic, though only under certain parameter regimes. The nature of this strategy's role, however, differs significantly between the two dynamics: in the SBD imitation dynamic, Informants are sufficient but not necessary to achieve the cooperative equilibrium, while under the best response dynamic, Informants are necessary but not sufficient for convergence to cooperation. Since a policy of simply converting citizens to Informants will not guarantee success under best response dynamics, we identify alternative strategies that may help the system reach cooperation in this case, e.g., the use of moderate but not too severe punishments on criminals.

Within the framework of variational modelling we derive a one-phase moving boundary problem describing the motion of a semipermeable membrane enclosing a viscous liquid, driven by osmotic pressure and surface tension of the membrane. For this problem we prove the existence of classical solutions for a short-time.

We study the behaviour of network diffusions based on the PageRank random walk from a set of seed nodes. These diffusions are known to reveal small, localized clusters (or communities), and also large macro-scale clusters by varying a parameter that has a dual-interpretation as an accuracy bound and as a regularization level. We propose a new method that quickly approximates the result of the diffusion for all values of this parameter. Our method efficiently generates an approximate solution path or regularization path associated with a PageRank diffusion, and it reveals cluster structures at multiple size-scales between small and large. We formally prove a runtime bound on this method that is independent of the size of the network, and we investigate multiple optimizations to our method that can be more practical in some settings. We demonstrate that these methods identify refined clustering structure on a number of real-world networks with up to 2 billion edges.

The sine-Gordon (SG) partial differential equation (PDE) with an arbitrary perturbation is initially considered. Using the method of Kuzmak–Luke, we investigate the conditions, which the perturbation must satisfy, for a breather solution to be a valid leading-order asymptotic approximation to the perturbed problem. We analyse the cases of both stationary and moving breathers. As examples, we consider perturbing terms which include typical linear damping, periodic sinusoidal driving, and dispersion. The motivation for this study is that the mathematical modelling of physical systems often leads to the discrete SG system of ordinary differential equations, which are then approximated in the long wavelength limit by the continuous SG PDE. Such limits typically produce fourth-order spatial derivatives as correction terms. The new results show that the stationary breather solution is a consistent solution of both the quasi-continuum SG equation and the forced/damped SG system. However, the moving breather is only a consistent solution of the quasi-continuum SG equation and not the damped SG system.

The discrete non-linear Schrödinger equation is one of the most important inherently discrete models, having a crucial role in the modelling of a great variety of phenomena, ranging from solid-state and condensed-matter physics to biology. In this paper, a class of discrete non-linear Schrödinger equations are considered. Using critical point theory in combination with periodic approximations, we establish some new sufficient conditions on the existence results for solitons of the equation. The classical Ambrosetti–Rabinowitz superlinear condition is improved.

In this paper, we propose a new algebraic method to construct non-autonomous discrete integrable systems. The method starts from constructing generalizations of convergence acceleration algorithms related to discrete integrable systems. Then the non-autonomous version of the corresponding integrable systems are derived. The molecule solutions of the systems are also obtained. As an example of the application of the method, we propose a generalization of the multistep ϵ-algorithm, and then derive a non-autonomous discrete extended Lotka–Volterra equation. Since the convergence acceleration algorithm from the lattice Boussinesq equation is just a particular case of the multistep ϵ-algorithm, we have therefore arrived at a generalization of this algorithm. Finally, numerical experiments on the new algorithm are presented.

This paper deals with a fully-coupled thermoelastic problem, in a heterogeneous medium, arising from the metallurgical industry. The aim is to prove regularity properties of the solution with respect to space and time. Regularity in space is obtained by means of regularity properties for elliptic operators. In order to prove regularity in time, a mathematical induction technique, together with an existence and uniqueness result for this type of problems, is applied.

We consider a storage allocation model with a finite number of storage spaces. There are m primary spaces that are ranked {1,2,. . .,m} and R secondary spaces ranked {m + 1, m + 2,. . .,m + R}. Items arrive according to a Poisson process, occupy a space for a random exponentially distributed time, and an arriving item takes the lowest ranked available space. Letting N1 and N2 denote the numbers of occupied primary and secondary spaces, we study the joint distribution Prob[N1 = k, N2 = r] in the steady state. The joint process (N1, N2) behaves as a random walk in a lattice rectangle. We shall obtain explicit expressions for the distribution of (N1, N2), as well as the marginal distribution of N2. We also give some numerical studies to illustrate the qualitative behaviors of the distribution(s). The main contribution is to study the effects of a finite secondary capacity R, whereas previous studies had R = ∞.

In this work, we analyse a broad class of generalized Poisson–Boltzmann equations and reveal a common mathematical structure. In the limit of a wide electrode, we show that a broad class of generalized Poisson–Boltzmann equations admits a reduction that affords an explicit connection between the functional form of the corresponding free energy and the associated differential capacitance data. We exploit the relation to we show that differential capacitance curves generically undergo an inflection transition with increasing salt concentration, shifting from a local minimum near the point of zero charge for dilute solutions to a local maximum point near the point of zero charge for concentrated solutions. In addition, we develop a robust numerical method for solving generalized Poisson–Boltzmann equations which is easily applicable to the broad class of generalized Poisson–Boltzmann equations with very few code adjustments required for each model

We introduce several novel and computationally efficient methods for detecting “core–periphery structure” in networks. Core–periphery structure is a type of mesoscale structure that consists of densely connected core vertices and sparsely connected peripheral vertices. Core vertices tend to be well-connected both among themselves and to peripheral vertices, which tend not to be well-connected to other vertices. Our first method, which is based on transportation in networks, aggregates information from many geodesic paths in a network and yields a score for each vertex that reflects the likelihood that that vertex is a core vertex. Our second method is based on a low-rank approximation of a network's adjacency matrix, which we express as a perturbation of a tensor-product matrix. Our third approach uses the bottom eigenvector of the random-walk Laplacian to infer a coreness score and a classification into core and peripheral vertices. We also design an objective function to (1) help classify vertices into core or peripheral vertices and (2) provide a goodness-of-fit criterion for classifications into core versus peripheral vertices. To examine the performance of our methods, we apply our algorithms to both synthetically generated networks and a variety of networks constructed from real-world data sets.

Recently, many security-related problems have gained increasing attention from a quantitative perspective. In this paper, we propose a game-theoretical approach to model the interaction between police forces and delinquents in public places. In the well-known Stackelberg game, a leader is faced with only one follower. However, in our application, the police are simultaneously faced with many offenders, who may be organized or act independently of each other. This application motivates the development of two games: a classical leader-follower interaction between police and organized criminals on the one hand and a novel approach between the leader and selfishly acting offenders on the other. It is of special interest that the effect of crime displacement under police surveillance be anticipated by the proposed models. Results using data from a simulated environment emphasise how these models can provide decision support for policing outperforming traditional strategies.

In this paper, we give exact solutions for the convective viscous Cahn--Hilliard equation. This equation with a general symmetric double-well potential and Burgers-type convective term was introduced by T. P. Witelski (1996 Studies in Applied Mathematics96, 277–300) to study the joint effects of nonlinear convection and viscosity. We consider this equation with a polynomial, generally asymmetric potential. We also consider both Burgers-type and cubic convective terms. We obtained exact travelling-wave solutions for both cases. For the former case, with an additional constraint on nonlinearity and viscosity, we also obtained an exact two-wave solution.

The central theme of complex systems research is to understand the emergent macroscopic properties of a system from the interplay of its microscopic constituents. The emergence of macroscopic properties is often intimately related to the structure of the microscopic interactions. Here, we present an analytical approach for deriving necessary conditions that an interaction network has to obey in order to support a given type of macroscopic behaviour. The approach is based on a graphical notation, which allows rewriting Jacobi's signature criterion in an interpretable form and which can be applied to many systems of symmetrically coupled units. The derived conditions pertain to structures on all scales, ranging from individual nodes to the interaction network as a whole. For the purpose of illustration, we consider the example of synchronization, specifically the (heterogeneous) Kuramoto model and an adaptive variant. The results complete and extend the previous analysis of Do et al. (2012Phys. Rev. Lett.108, 194102).

We present local existence and uniqueness results for the following 2 + 1 diffusive–dispersive equation due to P. Hall arising in modelling of river braiding:

$$\begin{equation*} u_{yyt} - \gamma u_{xxx} -\alpha u_{yyyy} - \beta u_{yy} + \left ( u^2 \right )_{xyy} = 0 \end{equation*}$$