The effect of signals on stability, stable throughput region, and delay in a two-user slotted ALOHA-based random-access system with collisions is considered. This work gives rise to the development of random access G-networks, which can model security attacks, expiration of deadlines, or other malfunctions, and introduce load balancing among highly interacting queues. The users are equipped with infinite capacity buffers accepting external bursty arrivals. We consider both negative and triggering signals. Negative signals delete a packet from a user queue, while triggering signals cause the instantaneous transfer of packets among user queues. We obtain the exact stability region, and show that the stable throughput region is a subset of it. Moreover, we perform a compact mathematical analysis to obtain exact expressions for the queueing delay by solving a non-homogeneous Riemann boundary value problem. A computationally efficient way to obtain explicit bounds for the expected number of buffered packets at user queues is also presented. The theoretical findings are numerically evaluated and insights regarding the system performance are derived.

This paper presents a control strategy for solving the cable pseudo-drag problem of cable-driven parallel camera robots at high speeds. The control strategy belongs to a hybrid position/tension control method based on cable tension optimization. The cable catenary model and cable pseudo-drag problem are considered firstly. Then, the dynamic model of the cable-driven parallel camera robot is established. The cable tension optimization is proposed. And then a control strategy is put forward and its stability is proved. Simulation results of a four-cable camera robot are presented and discussed.

A nonlinear mathematical model for innovation diffusion is proposed. The system of ordinary differential equations incorporates variable external influences (the cumulative density of marketing efforts), variable internal influences (the cumulative density of word of mouth) and a logistically growing human population (the variable potential consumers). The change in population density is due to various demographic processes such as intrinsic growth rate, emigration, death rate etc. Thus the problem involves two dynamic variables viz. a non-adopter population density and an adopter population density. The model is analysed qualitatively using the stability theory of differential equations, with the help of the corresponding characteristic equation of the system. The interior equilibrium point can be stable for all time delays to a critical value, beyond which the system becomes unstable and a Hopf bifurcation occurs at a second critical value. Employing normal form theory and a centre manifold theorem applicable to functional differential equations, we derive some explicit formulas determining the stability, the direction and other properties of the bifurcating periodic solutions. Our numerical simulations show that the system behaviour can become extremely complicated as the time delay increases, with a stable interior equilibrium point leading to a limit cycle with one local maximum and minimum per cycle (Hopf bifurcation), then limit cycles with more local maxima and minima per cycle, and finally chaotic solutions.

The fractional derivatives include nonlocal information and thus their calculation requires huge storage and computational cost for long time simulations. We present an efficient and high-order accurate numerical formula to speed up the evaluation of the Caputo fractional derivative based on the L2-1 σ formula proposed in [A. Alikhanov, J. Comput. Phys., 280 (2015), pp. 424-438], and employing the sum-of-exponentials approximation to the kernel function appeared in the Caputo fractional derivative. Both theoretically and numerically, we prove that while applied to solving time fractional diffusion equations, our scheme not only has unconditional stability and high accuracy but also reduces the storage and computational cost.

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.

In this study, we present a phase-field model that describes the process of intercalation of Li ions into a layer of an amorphous solid such as amorphous silicon (a-Si). The governing equations couple a viscous Cahn–Hilliard-Reaction model with elasticity in the framework of the Cahn–Larché system. We discuss the parameter settings and flux conditions at the free boundary that lead to the formation of phase boundaries having a sharp gradient in lithium ion concentration between the initial state of the solid layer and the intercalated region. We carry out a matched asymptotic analysis to derive the corresponding sharp-interface model that also takes into account the dynamics of triple points where the sharp interface intersects the free boundary of the Si layer. We numerically compare the interface motion predicted by the sharp-interface model with the long-time dynamics of the phase-field model.

Pressure-correction projection finite element methods (FEMs) are proposed to solve nonstationary natural convection problems in this paper. The first-order and second-order backward difference formulas are applied for time derivative, the stability analysis and error estimates of the semi-discrete schemes are presented using energy method. Compared with characteristic variational multiscale FEM, pressure-correction projection FEMs are more efficient and unconditionally energy stable. Ample numerical results are presented to demonstrate the effectiveness of the pressure-correction projection FEMs for solving these problems.

The computational work and storage of numerically solving the time fractional PDEs are generally huge for the traditional direct methods since they require total memory and work, where NT and NS represent the total number of time steps and grid points in space, respectively. To overcome this difficulty, we present an efficient algorithm for the evaluation of the Caputo fractional derivative of order α∈(0,1). The algorithm is based on an efficient sum-of-exponentials (SOE) approximation for the kernel t –1–α on the interval [Δt, T] with a uniform absolute error ε. We give the theoretical analysis to show that the number of exponentials N exp needed is of order for T≫1 or for TH1 for fixed accuracy ε. The resulting algorithm requires only storage and work when numerically solving the time fractional PDEs. Furthermore, we also give the stability and error analysis of the new scheme, and present several numerical examples to demonstrate the performance of our scheme.

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. (2012 Phys. Rev. Lett. 108, 194102).

We consider a fully practical finite element approximation of the Cahn–Hilliard–Stokes system:

$$\begin{align*} \gamma \tfrac{\partial u}{\partial t} + \beta v \cdot \nabla u - \nabla \cdot \left( \nabla w \right) & = 0 \,, \quad w= -\gamma \Delta u + \gamma ^{-1} \Psi ' (u) - \tfrac12 \alpha c'(\cdot,u) | \nabla \phi |^2\,, \\ \nabla \cdot (c(\cdot,u) \nabla \phi) & = 0\,,\quad \begin{cases} -\Delta v + \nabla p = \varsigma w \nabla u, \\ \nabla \cdot v = 0, \end{cases} \end{align*}$$

$\mathbb{R}_{>0}$

$\mathbb{R}_{\geq0}$

In this paper, we study two PDEs that generalize the urban crime model proposed by Short et al. (2008 Math. Models Methods Appl. Sci. 18, 1249–1267). Our modifications are made under assumption of the spatial heterogeneity of both the near-repeat victimization effect and the dispersal strategy of criminal agents. We investigate pattern formations in the reaction–advection–diffusion systems with non-linear diffusion over multi-dimensional bounded domains subject to homogeneous Neumann boundary conditions. It is shown that the positive homogeneous steady state loses its stability as the intrinsic near-repeat victimization rate ε decreases and spatially inhomogeneous steady states emerge through bifurcation. Moreover, we find the wavemode selection mechanism through rigorous stability analysis of these non-trivial spatial patterns, which shows that the only stable pattern must have a wavenumber that maximizes the bifurcation value. Based on this wavemode selection mechanism, we will be able to predict the formation of stable aggregates of the house attractiveness and criminal population density, at least when the diffusion rate ε is around the principal bifurcation value. Our theoretical results also suggest that large domains support more stable aggregates than small domains. Finally, we perform extensive numerical simulations over 1D intervals and 2D squares to illustrate and verify our theoretical findings. Our numerics also demonstrate the formation of other interesting patterns in these models such as the merging of two interior spikes and the emerging of new spikes, etc. These non-trivial solutions can model the well-observed phenomenon of aggregation in urban criminal activities.

A finite difference method which is second-order accurate in time and in space is proposed for two-dimensional fractional percolation equations. Using the Fourier transform, a general approximation for the mixed fractional derivatives is analyzed. An approach based on the classical Crank-Nicolson scheme combined with the Richardson extrapolation is used to obtain temporally and spatially second-order accurate numerical estimates. Consistency, stability and convergence of the method are established. Numerical experiments illustrating the effectiveness of the theoretical analysis are provided.

We consider in this paper numerical approximation of the linear Fluid-Structure Interaction (FSI). We construct a new class of pressure-correction schemes for the linear FSI problem with a fixed interface, and prove rigorously that they are unconditionally stable. These schemes are computationally very efficient, as they lead to, at each time step, a coupled linear elliptic system for the velocity and displacement in the whole region and a discrete Poisson equation in the fluid region.

In this paper, an adaptive force reflection scheme is proposed for bilateral teleoperation. In order to achieve an ideal telepresence performance while keeping the system stable, the force reflection algorithm needs to take both the human force and the contact force into consideration. An observer based on the feature of the human operator is designed to estimate the force applied on the master device. The reflected force is calculated by performing the orthogonal decomposition of the contact force, and is adjusted adaptively according to the estimated human force. The direction of the reflected force becomes a key consideration in the design process, so the proposed approach has an advantage in the guiding contact task. Based on the small gain theorem, the master device with the force reflection scheme is proved to be input-to-output stable, and the derivation for stability criterion of the closed-loop teleoperation system is also given. The results of simulations and experiments on a 6-degree of freedom teleoperation system demonstrate the effectiveness of the proposed scheme.

In this paper we explore the eco-evolutionary dynamics of a predator-prey model, where the prey population is structured according to a certain life history trait. The trait distribution within the prey population is the result of interplay between genetic inheritance and mutation, as well as selectivity in the consumption of prey by the predator. The evolutionary processes are considered to take place on the same time scale as ecological dynamics, i.e. we consider the evolution to be rapid. Previously published results show that population structuring and rapid evolution in such predator-prey system can stabilize an otherwise globally unstable dynamics even with an unlimited carrying capacity of prey. However, those findings were only based on direct numerical simulation of equations and obtained for particular parameterizations of model functions, which obviously calls into question the correctness and generality of the previous results. The main objective of the current study is to treat the model analytically and consider various parameterizations of predator selectivity and inheritance kernel. We investigate the existence of a coexistence stationary state in the model and carry out stability analysis of this state. We derive expressions for the Hopf bifurcation curve which can be used for constructing bifurcation diagrams in the parameter space without the need for a direct numerical simulation of the underlying integro-differential equations. We analytically show the possibility of stabilization of a globally unstable predator-prey system with prey structuring. We prove that the coexistence stationary state is stable when the saturation in the predation term is low. Finally, for a class of kernels describing genetic inheritance and mutation we show that stability of the predator-prey interaction will require a selectivity of predation according to the life trait.

In this paper, a new discontinuous Galerkin method is developed for the parabolic equation with jump coefficients satisfying the continuous flow condition. Theoretical analysis shows that this method is L2 stable. When the finite element space consists of interpolative polynomials of degrees k, the convergent rate of the semi-discrete discontinuous Galerkin scheme has an order of . Numerical examples for both 1-dimensional and 2-dimensional problems demonstrate the validity of the new method.

This paper presents new proofs of four stability properties (semiglobal strict passivity, semiglobal asymptotic stability, semiglobal input-to-state stability, and semiglobal uniform ultimate boundedness with an arbitrarily reducible ultimate bound) of a rigid-link manipulator under proportional-integral-derivative (PID) position control. The proofs employ a strict Lyapunov function and a novel parameterization to provide four inequality conditions for the stability properties. In those inequalities, arithmetic operations on physical quantities are physically consistent if the joints are all revolute or all prismatic. A gain selection procedure is presented by which the ultimate bounds of velocity error, position error, and its integral can be independently designed.

This paper presents an adaptive force/position controller for a parallel robot executing constrained motions. This controller, based on an MLPNN (or Multi-Layer Perceptron Neural Network), does not require the inverse dynamic model of the robot to derive the control law. A neural identification of the dynamic model of the robot is proposed to determine the principal components of the MLPNN input vector. The latter is used to compensate the dynamic effects arising from the robot–environment interaction and its parameters are adjusted according to an adaptation law based on the Lyapunov-analysis methodology. The proposed controller is evaluated experimentally on the C5 parallel robot. This method is capable of tracking accurately the force/position trajectories and its stability robustness is proved.