Amaranth, with its high genetic variability, holds promise for global food security, income generation and climate resilience. Developing stable, high-yielding genotypes is essential for sustainable production. In this study, stability analysis was conducted on five Amaranth accessions over two seasons at three Malawian sites. Significant trait variations, including grain yield, plant height and leaf characteristics, underscored the dynamic nature of Amaranth cultivation. Notably, LL-BH-04 consistently exhibited superior grain yield, while others showed variable performance, highlighting the importance of stability analysis. Employing the Eberhart and Russell model, stable accessions in leaf and grain yield were identified. Additionally, AMMI (Additive Main effects and Multiplicative Interaction) biplot analysis revealed genetic diversity and stability patterns, aiding resilient cultivar selection. Consequently, LL-BH-04 and PE-UP-BH-01, identified as stable genotypes, were recommended for release, thereby enhancing agricultural sustainability and food security. These findings emphasize the need for site-specific breeding evaluations for sustainable productivity and underscore the importance of selecting stable cultivars to address agricultural challenges. LL-BH-04 and PE-LO-BH-01 were proposed for release to boost Amaranth production in Malawi, serving as the foundation for tailored breeding efforts aimed at improving productivity and resilience. This study contributes valuable insights into the stability and performance of Amaranth cultivars, offering guidance for sustainable crop production and variety development strategies.

This paper presents some of the first results of global linear stability analyses performed using a bespoke eigensolver that has recently been implemented in the next generation flow solver framework CODA. The eigensolver benefits from the automatic differentiation capability of CODA that allows computation of the exact product of the Jacobian matrix with an arbitrary complex vector. It implements the Krylov–Schur algorithm for solving the eigenvalue problem. The bespoke tool has been validated for the case of laminar flow past a circular cylinder with numerical results computed using the TAU code and those reported in the literature. It has been applied with both second-order finite volume and high-order discontinuous Galerkin schemes for the case of laminar flow past a square cylinder. It has been demonstrated that using high-order schemes on coarser grids leads to well-converged eigenmodes with a shorter computation time compared to using second-order schemes on finer grids.

In this study, a non-linear deterministic model for the transmission dynamics of skin sores (impetigo) disease is developed and analysed by the help of stability of differential equations. Some basic properties of the model including existence and positivity as well as boundedness of the solutions of the model are investigated. The disease-free and endemic equilibrium were investigated, as well as the basic reproduction number, R0, also calculated using the next-generation matrix approach. When R0 < 1, the model's stability analysis reveals that the system is asymptotically stable at disease-free critical point globally as well as locally. If R0 > 1, the system is asymptotically stable at disease-endemic equilibrium both locally and globally. The long-term behaviour of the skin sores model's steady-state solution in a population is investigated using numerical simulations of the model.

In this paper, we present a systematic approach to improve the understanding of stability and robustness of stability against the external disturbances of a passive biped walker. First, a multi-objective, multi-modal particle swarm optimization (MOMM-PSO) algorithm was employed to suggest the appropriate initial conditions for a given biped walker model to be stable. The MOMM-PSO with ring topology and special crowding distance (SCD) used in this study can find multiple local minima under multiple objective functions by limiting each agent’s search area properly without determining a large number of parameters. Second, the robustness of stability under external disturbances was studied, considering an impact in the angular displacement sampled from the probabilistic distribution. The proposed systematic approach based on MOMM-PSO can find multiple initial conditions that lead the biped walker in the periodic gait, which could not be found by heuristic approaches in previous literature. In addition, the results from the proposed study showed that the robustness of stability might change depending on the location on a limit cycle where immediate angular displacement perturbation occurs. The observations of this study imply that the symmetry of the stable region about the limit cycle will break depending on the accelerating direction of inertia. We believe that the systematic approach developed in this study significantly increased the efficiency of finding the appropriate initial conditions of a given biped walker and the understanding of robustness of stability under the unexpected external disturbance. Furthermore, a novel methodology proposed for biped walkers in the present study may expand our understanding of human locomotion, which in turn may suggest clinical strategies for gait rehabilitation and help develop gait rehabilitation robotics.

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 Nexp 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. (2012Phys. 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.