Abstract In this research, the (2+1)‐dimensional normal biological population model, incorporating the Verhulst law for population growth, is employed to explore species population dynamics. Employing Lie symmetry analysis, we address a nonlinear degenerate parabolic partial differential equation, yielding much‐improved results. This analysis includes computing one‐dimensional optimal subalgebras, reduced ordinary differential equations, and obtaining invariant solutions with a visual depiction of the physical behavior of the Verhulst biological population model through symmetry group transformations. Additionally, the multiplier method leads to novel conservation laws and potential systems not locally connected to the governing partial differential equation (PDE). These findings have significant implications for understanding and controlling biological populations, offering insights for applications in ecology and the environment.

Abstract Forest fire spreading is a complex phenomenon characterized by a stochastic behavior. Nowadays, the enormous quantity of georeferenced data and the availability of powerful techniques for their analysis can provide a very careful picture of forest fires opening the way to more realistic models. We propose a stochastic spreading model continuous in space and time that has the potentiality to use such data in their full power. The state of the forest fire is described by the subprobability densities of the green trees and of the trees on fire that can be estimated thanks to data coming from satellites and earth detectors. The fire dynamics is encoded into a kernel function that can take into account wind conditions, land slope, spotting phenomena, and so on, bringing to a system of integrodifferential equations whose solutions provide the evolution in time of the subprobability densities. That makes the model complementary to models based on cellular automata that furnish single instantiations of the stochastic phenomenon. Moreover, stochastic models based on cellular automata can be derived from the present model by space and time discretization. Existence and uniqueness of the solutions is proved by using Banach's fixed‐point theorem. The asymptotic behavior of the model is analyzed as well. By specifying a particular structure for the kernel, we obtain numerical simulations of the fire spreading under different conditions. For example, in the case of a forest fire evolving toward a river, the simulations show that the probability density of the trees on fire is different from zero beyond the river due to the spotting phenomenon. The kernel could be slightly modified to include firefighters interventions and weather changes.

Abstract The aim of this paper is to present a survey and a detailed numerical study on a remarkable Boussinesq system describing weakly nonlinear, long surface water waves. In the one‐dimensional case, this system can be viewed as a dispersive perturbation of the hyperbolic Saint‐Venant (shallow water) system. The asymptotic stability of the solitary waves is numerically established. Blow‐up of solutions for initial data not satisfying the noncavitation condition as well as the appearance of dispersive shock waves are studied.

Abstract Reaction–diffusion equations (RDEs) model the spatiotemporal evolution of a density field according to diffusion and net local changes. Usually, the diffusivity is positive for all values of , which causes the density to disperse. However, RDEs with partially negative diffusivity can model aggregation, which is the preferred behavior in some circumstances. In this paper, we consider a nonlinear RDE with quadratic diffusivity that is negative for . We use a nonclassical symmetry to construct analytic receding time‐dependent, colliding wave, and receding traveling wave solutions. These solutions are multivalued, and we convert them to single‐valued solutions by inserting a shock. We examine properties of these analytic solutions including their Stefan‐like boundary condition, and perform a phase plane analysis. We also investigate the spectral stability of the and constant solutions, and prove for certain and that receding traveling waves are spectrally stable. In addition, we introduce a new shock condition where the diffusivity and flux are continuous across the shock. For diffusivity symmetric about the midpoint of its zeros, this condition recovers the well‐known equal‐area rule, but for nonsymmetric diffusivity it results in a different shock position.

Abstract In this paper, a compound equation which is a mix of the Wadati–Konno–Ichikawa (WKI) equation and the short‐pulse (SP) equation is first studied. By transforming both the independent and dependent variables in the equation, we introduce a novel hodograph transformation to convert the compound WKI–SP equation into the mKdV–SG (modified Korteweg–de Vries and sine‐Gordon) equation. The multiloop soliton solutions in the form of the parametric representation are found. It is shown that the ‐loop soliton solution may be decomposed exactly into separate soliton elements by using a Moloney–Hodnett‐type decomposition. By virtue of the decomposed soliton solutions, the asymptotic behaviors of and are investigated in detail. The corresponding phase shifts of each loop or antiloop soliton caused by its interaction with the other ones are calculated. Furthermore, a new hierarchy of WKI–SP‐type equations possessing multiloop soliton solutions is constructed. These deduced equations are all with time‐varying coefficients and the corresponding dispersion relation will have a time‐dependent velocity. The whole hierarchy of equations which include the WKI‐type equations, the SP‐type equations, and the compound generalized WKI–SP equations, are illustrated Lax integrable. The specific equation in the hierarchy is labeled as equation so that its Lax pairs can be directly written out with the help of and . A unified hodograph transformation is established to relate the compound WKI–SP hierarchy with the mKdV–SG hierarchy.

Abstract In this paper, we study a predator–prey mite model of Leslie type with generalized Holling IV functional response. The model is shown to have very rich bifurcation dynamics, including subcritical and supercritical Hopf bifurcations, degenerate Hopf bifurcation, focus‐type and cusp‐type degenerate Bogdanov–Takens bifurcations of codimension 3, originating from a nilpotent focus or cusp of codimension 3 that acts as the organizing center for the bifurcation set. Coexistence of multiple steady states, multiple limit cycles, and homoclinic cycles is also found. Interestingly, the coexistence of two limit cycles is guaranteed by investigating generalized Hopf bifurcation and degenerate homoclinic bifurcation, and we also find that two generalized Hopf bifurcation points are connected by a saddle‐node bifurcation curve of limit cycles, which indicates the existence of global regime for two limit cycles. Our work extends some results in the literature.

Abstract In this paper, taking into account the inevitable impact of environmental perturbations on disease transmission, we primarily investigate a stochastic model for measles infection with nonlinear incidence. The transmission rate in this model follows a logarithmic normal distribution influenced by an Ornstein–Uhlenbeck (OU) process. To analyze the dynamic properties of the stochastic model, our first step is to establish the existence and uniqueness of a global solution for the stochastic equations. Next, by constructing appropriate Lyapunov functions and utilizing the ergodicity of the OU process, we establish sufficient conditions for the existence of a stationary distribution, indicating the prevalence of the disease. Furthermore, we provide sufficient conditions for disease elimination. These conditions are derived by considering the interplay between the model parameters and the stochastic dynamics. Finally, we validate the theoretical conclusions through numerical simulations, which allow us to assess the practical implications of the established conditions and observe the dynamics of the stochastic model in action. By combining theoretical analysis and numerical simulations, we gain a comprehensive understanding of the stochastic model's behavior, contributing to the broader understanding of measles transmission dynamics and the development of effective control strategies.

Abstract The sigmoid functional response in the predator–prey model was posed in 1977. But its dynamics has not been completely characterized. This paper completes the classification of the global dynamics for the classical predator–prey model with the sigmoid functional response, whose denominator has two different zeros. The dynamical phenomena we obtain here include global stability, the existence of the heteroclinic and homoclinic loops, the consecutive canard explosions via relaxation oscillation, and the canard explosion to a homoclinic loop among others. As we know, the last one is a new dynamical phenomenon, which has never been reported previously. In addition, with the help of geometric singular perturbation theory, we solve the problem of connection between stable and unstable manifolds from different singularities, which has not been well settled in the published literature.

Abstract Many reaction–diffusion models produce traveling wave solutions that can be interpreted as waves of invasion in biological scenarios such as wound healing or tumor growth. These partial differential equation models have since been adapted to describe the interactions between cells and extracellular matrix (ECM), using a variety of different underlying assumptions. In this work, we derive a system of reaction–diffusion equations, with cross‐species density‐dependent diffusion, by coarse‐graining an agent‐based, volume‐filling model of cell invasion into ECM. We study the resulting traveling wave solutions both numerically and analytically across various parameter regimes. Subsequently, we perform a systematic comparison between the behaviors observed in this model and those predicted by simpler models in the literature that do not take into account volume‐filling effects in the same way. Our study justifies the use of some of these simpler, more analytically tractable models in reproducing the qualitative properties of the solutions in some parameter regimes, but it also reveals some interesting properties arising from the introduction of cell and ECM volume‐filling effects, where standard model simplifications might not be appropriate.

Abstract We introduce a new family of multiple orthogonal polynomials satisfying orthogonality conditions with respect to two weights on the positive real line, with the gamma density and a density related to the exponential integral . We give explicit formulas for the type I functions and type II polynomials, their Mellin transform, Rodrigues formulas, hypergeometric series, and recurrence relations. We determine the asymptotic distribution of the (scaled) zeros of the type II multiple orthogonal polynomials and make a connection to random matrix theory. Finally, we also consider two related families of mixed‐type multiple orthogonal polynomials.