Search

In a previous paper (Feng et al., J. Differential Equations (2023)), we studied a seasonally interactive model between closed seasons and open seasons with Michaelis-Menten type harvesting, in which we assumed that the harvesting quantity was relatively large (0 < κ = l G c E < 1) and obtained a length threshold of the closed season T ¯ ⁎ depending on the harvesting parameter and the seasonal fluctuation period. It was shown that the origin is globally asymptotically stable if and only if T ¯ ≤ T ¯ ⁎ , and there exists a unique globally asymptotically stable T -periodic solution if and only if T ¯ > T ¯ ⁎. In this paper, we continue to investigate the periodic dynamics of this model when the harvesting quantity is relatively small (κ = l G c E ≥ 1). By finding another smaller length threshold T ¯ ⁎ ∈ (0 , T ¯ ⁎) , we determine the number of periodic solutions and study their stability, which imply the occurrence of bifurcation of periodic solutions. Our results demonstrate that designing the closed harvesting season properly can prevent the species from extinction. Moreover, by comparing with the continuous harvesting model and combining the results in our previous paper and this article, we provide a complete understanding on the global dynamics for the periodic switching model. Some numerical simulations are also carried out to illustrate the theoretical results. [ABSTRACT FROM AUTHOR]

The aim of this paper is to obtain the sharp estimate for the lowest positive eigenvalue for the Camassa-Holm equation y ″ = 1 4 y + λ m (t) y , with the Neumann-Dirichlet boundary conditions, where potential m admits to change sign. We first study the optimal lower bound for the smallest positive eigenvalue in the measure differential equations. Then based on the relationship between the minimization problem of the smallest positive eigenvalue for the ODE and the one for the MDE, we find the explicit optimal lower bound of the smallest positive eigenvalue for this indefinite Camassa-Holm equation. [ABSTRACT FROM AUTHOR]

In this paper, we construct higher-order modular differential equations for elliptic modular forms, holomorphic Jacobi forms, and skew-holomorphic Jacobi forms by using the Rankin-Cohen brackets extended to the Eisenstein series of weight two. Rankin-Cohen brackets are typical tools in the construction of modular forms by differential operators. We introduce the extended Rankin-Cohen brackets for holomorphic or skew-holomorphic Jacobi forms paired with elliptic modular forms. They provide a unified method of construction of higher-order modular differential equations for the modular forms above. In particular, with these extended Rankin-Cohen brackets, we can reproduce the modular differential equations of previous studies. [ABSTRACT FROM AUTHOR]

In this paper, we study the existence, multiplicity and dynamics of positive periodic solutions to a generalized Liénard equation with repulsive singularities. The Ambrosetti-Prodi type result is proved in the absence of the so-called anticoercivity condition. Furthermore, with s as a parameter, under some conditions on the function h , it has been shown that for any M > 1 there exists s M ∈ R such that the equation x ″ + f (x) x ′ + h (t , x) = s has two positive T -periodic solutions u 1 (⋅ ; s) and u 2 (⋅ ; s) satisfying min ⁡ { u 1 (t ; s) : t ∈ [ 0 , T ] } > M and min ⁡ { u 2 (t ; s) : t ∈ [ 0 , T ] } < 1 / M for every s < s M. As a by-product of the property, we obtain sufficient conditions to guarantee the existence of positive T -periodic solutions of indefinite differential equations. • Establish a refinement result of the Ambrosetti-Prodi type in the absence of anti-coercivity condition. • Investigate the dynamic behaviors of positive periodic solutions to a generalized Liénard equation. • Provide a new approach to study the existence of positive periodic solutions for indefinite differential equations. • Find a relationship between the strict lower functions and Leray-Schauder degree. [ABSTRACT FROM AUTHOR]

The main aim of this article is to establish an effective criterion for a positive parametric topological entropy to differential equations with multivalued nonautonomous impulses on tori. The crucial tool for this goal consists in developing an appropriate Ivanov-like inequality for the lower estimate of a new kind of entropy by means of the logarithm of asymptotic Nielsen numbers for the compositions of a one-parameter family of multivalued admissible maps. This inequality is then applied to impulsive differential equations on tori via the associated Poincaré translation operators along their trajectories. The obtained results generalize in a significant way those in our recent papers, especially in the one with the same title (whence the indication by II), into a multivalued setting. Another nontrivial generalization can be regarded with respect to their simple reduction into a "nonparametric" case. [ABSTRACT FROM AUTHOR]

In this paper, we study exponentially fast mixing induced/enhanced by gravity and stochastic boundary in the kinetic theory of Vlasov equations. We consider the Vlasov equations with and without a vertical magnetic field inside a horizontally-periodic 3D half-space equipped with a non-isothermal diffusive reflection boundary condition of bounded continuous boundary temperature at the bottom. We construct both stationary solutions and global-in-time dynamical solutions in L ∞. We prove that moments of a dynamical fluctuation around the steady solutions decay exponentially fast in L ∞. As a key of this proof, we establish a uniform bound of so-called residual measures independently of the bouncing number of stochastic characteristics, by constructing a continuous stationary outgoing boundary flux which is strictly positive almost everywhere. [ABSTRACT FROM AUTHOR]

This paper is concerned with the asymptotic expansion of solutions to the initial-value problem of u ″ (t) + A u (t) + b (t) u ′ (t) = 0 in a Hilbert space with a nonnegative selfadjoint operator A and a coefficient b (t) ∼ (1 + t) − β (− 1 < β < 1). In the case b (t) ≡ 1 , it is known that the higher-order asymptotic profiles are determined via a family of first-order differential equations of the form v ′ (t) + A v (t) = F n (t) (Sobajima (2021) [10]). For the time-dependent case, it is only known that the asymptotic behavior of such a solution is given by the one of b (t) v ′ (t) + A v (t) = 0. The result of this paper is to find the equations for all higher-order asymptotic profiles. It is worth noticing that the equation for n -th order profile u ˜ n is given via v ′ (t) + m n (t) A v (t) = F n (t) which coefficient m n (time-scale) differs each other. [ABSTRACT FROM AUTHOR]

This paper is devoted to the smooth and stationary Wong-Zakai approximations for a class of rough differential equations driven by a geometric fractional Brownian rough path ω with Hurst index H ∈ (1 3 , 1 2 ]. We first construct the approximation ω δ of ω by probabilistic arguments, and then use the rough path theory to obtain the Wong-Zakai approximation for the solution on any finite interval. Finally, both the original system and the approximative system generate a continuous random dynamical systems φ and φ δ. As a consequence of the Wong-Zakai approximation of the solution, φ δ converges to φ as δ → 0. [ABSTRACT FROM AUTHOR]

We deal with the following fractional Choquard equation (− Δ) s u + V (x) u = (I μ ⁎ | u | 2 μ , s ⁎ ) | u | 2 μ , s ⁎ − 2 u , x ∈ R N , where I μ (x) is the Riesz potential, s ∈ (0 , 1) , 2 s < N ≠ 4 s , 0 < μ < min ⁡ { N , 4 s } and 2 μ , s ⁎ = 2 N − μ N − 2 s is the fractional critical Hardy-Littlewood-Sobolev exponent. By combining variational methods and the Brouwer degree theory, we investigate the existence and multiplicity of positive bound solutions to this equation when V (x) is a positive potential bounded from below. The results obtained in this paper extend and improve some recent works in the case where the coefficient V (x) vanishes at infinity. [ABSTRACT FROM AUTHOR]

Exact solutions of the SEIR epidemic model are derived, and various properties of solutions are obtained directly from the exact solution. In this paper Abel differential equations play an important role in establishing the exact solution of SEIR differential system, in particular the number of infected individuals can be represented in a simple form by using a positive solution of an Abel differential equation. It is shown that the parametric form of the exact solution satisfies some linear differential system including a positive solution of an Abel differential equation. [ABSTRACT FROM AUTHOR]

This paper addresses the problem of computing valuations of the Dieudonné determinants of matrices over discrete valuation skew fields (DVSFs). Under a reasonable computational model, we propose two algorithms for a class of DVSFs, called split. Our algorithms are extensions of the combinatorial relaxation of Murota (1995) and the matrix expansion by Moriyama and Murota (2013) , both of which are based on combinatorial optimization. While our algorithms require an upper bound on the output, we give an estimation of the bound for skew polynomial matrices and show that the estimation is valid only for skew polynomial matrices. We consider two applications of this problem. The first one is the noncommutative weighted Edmonds' problem (nc-WEP), which is to compute the degree of the Dieudonné determinants of matrices having noncommutative symbols. We show that the presented algorithms reduce the nc-WEP to the unweighted problem in polynomial time. In particular, we show that the nc-WEP over the rational field is solvable in time polynomial in the input bit-length. We also present an application to analyses of degrees of freedom of linear time-varying systems by establishing formulas on the solution spaces of linear differential/difference equations. [ABSTRACT FROM AUTHOR]

We study the existence and uniqueness of solution of a evolutionary partial differential equation originating from the continuum limit of a coupled process of totally asymmetric simple exclusion process (TASEP) and Langmuir kinetics (LK). In the fields of physics and biology, the TASEP-LK coupled process has been extensively studied by Monte Carlo simulations, numerical computations, and detailed experiments. However, no rigorous mathematical analysis so far has been given for the corresponding differential equations, especially the existence and uniqueness of their solutions. In this paper, we prove the existence of the C ∞ [ 0 , 1 ] steady-state solution by the method of upper and lower solution, and the uniqueness in both W 1 , 2 (0 , 1) and L ∞ (0 , 1) by a generalized maximum principle. We further prove the global existence and uniqueness of the time-dependent solution in C ([ 0 , 1 ] × [ 0 , + ∞)) ∩ C 2 , 1 ([ 0 , 1 ] × (0 , + ∞)) , which, for any continuous initial value, converges to the steady-state solution in C [ 0 , 1 ] (global attractivity). Our results support the numerical calculations and Monte Carlo simulations, and provide theoretical foundations for the TASEP-LK coupled process, especially the most important phase diagram of particle density along the travel track under different model parameters, which is difficult because the boundary layers (at one or both boundaries) and domain wall (separating high and low particle densities) may appear as the length of the travel track tends to infinity. The methods used in this paper may be instructive for studies of the more general cases of the TASEP-LK process, such as the one with multiple travel tracks and/or multiple particle species. [ABSTRACT FROM AUTHOR]

This paper is concerned with the invariant measures of non-autonomous stochastic reaction-diffusion equations on unbounded domains. We firstly investigate the existence of invariant measures for random dynamical systems over two parametric spaces as well as periodic invariant measures, and then discuss the convergence of (periodic) invariant measures with respect to some system parameters. Random Liouville type equation is also derived for invariant measures associated with non-autonomous random differential equations. Based on such abstract results, we prove the existence of (periodic) invariant measures for non-autonomous stochastic reaction-diffusion equation on unbounded domains, which are supported on the pullback attractors. And we further show every limit of a sequence of (periodic) invariant measures must be the (periodic) invariant measure of the corresponding limiting equation when the noise intensity varies in finite interval. Moreover, we prove that the invariant measures of the underlying equation satisfy a stochastic Liouville type equation. [ABSTRACT FROM AUTHOR]

It is shown in this paper that suitable weak solutions to the 6D steady incompressible Navier-Stokes are Hölder continuous at 0 provided that ∫ B 1 | u (x) | 3 d x + ∫ B 1 | f (x) | q d x or ∫ B 1 | ∇ u (x) | 2 d x + ∫ B 1 | ∇ u (x) | 2 d x (∫ B 1 | u (x) | d x) 2 + ∫ B 1 | f (x) | q d x with q > 3 is sufficiently small, which implies that the 2D Hausdorff measure of the set of singular points is zero. For the boundary case, we also obtain that 0 is regular provided that ∫ B 1 + | u (x) | 3 d x + ∫ B 1 + | f (x) | 3 d x or ∫ B 1 + | ∇ u (x) | 2 d x + ∫ B 1 + | f (x) | 3 d x is sufficiently small. These results improve previous regularity theorems by Dong-Strain ([8] , Indiana Univ. Math. J., 2012), Dong-Gu ([7] , J. Funct. Anal., 2014), and Liu-Wang ([27] , J. Differential Equations, 2018), where either the smallness of the pressure or the smallness of the scaling invariant quantities on all balls is necessary. [ABSTRACT FROM AUTHOR]

This paper aims to present a comprehensive framework for coupling tumor-bone remodeling processes in a 2-dimensional geometry. This is achieved by introducing a bio-inspired damage that represents the growing tumor, which subsequently affects the main populations involved in the remodeling process, namely, osteoclasts, osteoblasts, and bone tissue. The model is constructed using a set of differential equations based on the Komarova's and Ayati's models, modified to incorporate the bio-inspired damage that may result in tumor mass formation. Three distinct models were developed. The first two models are based on the Komarova's governing equations, with one demonstrating an osteolytic behavior and the second one an osteoblastic model. The third model is a variation of Ayati's model, where the bio-inspired damage is induced through the paracrine and autocrine parameters, exhibiting an osteolytic behavior. The obtained results are consistent with existing literature, leading us to believe that our in-silico experiments will serve as a cornerstone for paving the way towards targeted interventions and personalized treatment strategies, ultimately improving the quality of life for those affected by these conditions. [ABSTRACT FROM AUTHOR]

The famous J.C.P. Miller formula provides a recurrence algorithm for the composition B a ∘ f , where B a is the formal binomial series and f is a formal power series, however it requires that f has to be a nonunit. In this paper we provide the general J.C.P. Miller formula which eliminates the requirement of nonunitness of f and, instead, we establish a necessary and sufficient condition for the existence of the composition B a ∘ f. We also provide the general J.C.P. Miller recurrence algorithm for computing the coefficients of that composition, if B a ∘ f is well defined, obviously. Our generalizations cover both the case in which f is a one–variable formal power series and the case in which f is a multivariable formal power series. In the central part of this article we state, using some combinatorial techniques, the explicit form of the general J.C.P. Miller formula for one-variable case. As applications of these results we provide an explicit formula for the inverses of polynomials and formal power series for which the inverses exist, obviously. We also use our results to investigation of approximate solution to a differential equation which cannot be solved in an explicit way. [ABSTRACT FROM AUTHOR]

The paper aims to study the geometry and physics of the Raychaudhuri equation (RE) in the background of a homogeneous and anisotropic space–time described by Kantowski–Sachs (KS) metric. Role of anisotropy/shear in the context of convergence and possible avoidance of singularity has been analyzed subject to a physically motivated constraint. Moreover, using a suitable transformation the first order RE has been converted to a second order differential equation analogous to a Harmonic Oscillator and criterion for convergence has been shown to be associated with the time varying frequency of the Oscillator. Finally, the paper points out a geometric and physical notion of anisotropy along with their corresponding behavior towards convergence. [ABSTRACT FROM AUTHOR]

In this paper the half-linear differential equation with one-dimensional p -Laplacian (p = α + 1), (| u ′ | α sgn u ′) ′ = α (λ α + 1 + b (t)) | u | α sgn u , t ≥ t 0 , is considered, where α and λ are positive constants. It is proved that if the function b (t) is absolutely integrable on [ t 0 , ∞) , then the above equation has two nonoscillatory solutions u + (t) and u − (t) such that u ± (t) ∼ c e ± λ t and u ± ′ (t) ∼ ± c λ e ± λ t (t → ∞), and moreover, any nontrivial solution u (t) of the above equation satisfies either u (t) ∼ c e λ t , u ′ (t) ∼ c λ e λ t (t → ∞) or u (t) ∼ c e − λ t , u ′ (t) ∼ − c λ e − λ t (t → ∞). Here, c is a nonzero constant. A weaker result is also shown under the weaker condition that b (t) is conditionally integrable on [ t 0 , ∞). [ABSTRACT FROM AUTHOR]

• The normal form formula of double-Hopf bifurcation is generated for general differential equations with the nonlocality. • The formula is applicable to functional differential equations and partial differential equations with or without the nonlocality. • The algorithm is favourable to reveal complex dynamics, e.g. periodic, quasi-periodic and multi-periodic oscillations. • Stable spatially nonhomogeneous periodic and quasi-periodic solutions are found in a nonlocal predator-prey model. This paper is concerned with the calculation of the normal form on the center manifold up to the third-order term at the double-Hopf singularity for general partial functional differential equations with the nonlocal effect. We derive explicit formulas of the normal form, that can be applied for both functional differential equations and partial differential equations with or without nonlocal effects in a bounded spatial domain. This provides an effective tool to establish existences of multi-periodic and quasi-periodic oscillations of a double-Hopf singularity for such equations. As an example, the Holling-Tanner predator-prey model with the nonlocal intraspecific competition of prey is considered. It turns out that double-Hopf bifurcation will occur because of the nonlocal interaction. Many spatio-temporal dynamics are found, including stable spatially homogeneous or nonhomogeneous periodic solutions, and stable spatially nonhomogeneous quasi-periodic solutions. [ABSTRACT FROM AUTHOR]

In this paper we associate to a linear differential equation with coefficients in the field of Laurent formal power series a new geometric object, a spectrum in the sense of Berkovich. We compute this spectrum and show that it contains interesting information about the equation. [ABSTRACT FROM AUTHOR]

In this paper we present a number of examples of exact solutions for the Friedmann cosmological equation for metric F (R) gravity model. Emphasis was placed on the possibility of obtaining exact time dependences of the main cosmological physical quantities: scale factor, scalar curvature, Hubble rate and function F (R). For this purpose an ansatz was used to reduce the Friedmann equation to an ordinary differential equation for function F = F (H 2). This made it possible to obtain a number of exact solutions, both already known and new. • Exact solutions for the Friedmann cosmological equation for metric F (R) gravity. • Friedmann–Lemaître–Robertson–Walker (FLRW) metric interval. • An ansatz used to reduce the Friedmann equation to an ordinary differential equation. • Cosmological interpretation of solution obtained. [ABSTRACT FROM AUTHOR]

Recently, there has been a growing interest in the development of intelligent structural damage identification methods that utilize machine learning techniques. However, many existing machine learning-based damage identification methods rely on labeled datasets, which can be large-scale and therefore costly, time-consuming, and/or unrealistic to obtain in practice. To address this issue, a baseline-free damage identification method for plate structures is proposed in this paper that utilizes physics-informed neural networks (PINNs). PINNs are a unique type of neural network that embed physical laws, frequently represented as differential equations, into their loss function. The method assumes that a plate structure in its pristine state is an isotropic and homogeneous thin plate, which can be well modeled by the Kirchhoff–Love plate theory. The proposed damage identification method aims to isolate and intensify local anomalies caused by damage in flexural guided wavefields of damaged thin plates. To accomplish this goal for such a damaged plate, a PINN is constructed and trained using data from a measured flexural guided wavefield. Such a trained PINN is referred to as a pseudo-pristine model, as it can generate a wavefield approximating that of an isotropic and homogeneous thin plate. This unique functionality arises from penalizing the governing equation of the plate in the loss function of the PINN and the fact that the governing equation does not consider the existence of the damage. Any local anomalies in the measured wavefield can be isolated by comparing it with the wavefield generated by the pseudo-pristine model, and they are then intensified using the Teager energy operator. An accumulative damage index is formulated and the damage can be identified within neighborhoods with high index values. The effectiveness and noise-robustness of the proposed method are demonstrated through numerical and experimental examples. Parameter studies are conducted to investigate the robustness of the PINNs with different hyper-parameters. The results show that the damage can be accurately identified by the proposed method under the different hyper-parameters. • PINNs applied to isotropic and homogeneous thin plates with inverse problems. • Baseline-free damage identification method using PINNs is proposed. • Pseudo-pristine models generate wavefields approximating those of pristine plates. • Local anomalies in measured wavefields are isolated by pseudo-pristine models. • Teager energy of isolated local anomalies used for damage identification. [ABSTRACT FROM AUTHOR]

• We proposed variable linear transformation improved physics-informed neural networks (VLT-PINNs) to solve thin-layer flow forward and inverse problems. • Three principles for determining the VLT parameters were identified for three types of variables in PINNs respectively. • The VLT-PINNs with the principle-suggested parameters show excellent performance over the ten test cases. • Only under the constraints of the VLT principles can we obtain satisfactory results, while traditional linear transformation may yield totally wrong results. • VLT method is an attempt to combine the three advantages of accuracy, universality and simplicity. Physics-informed neural networks (PINNs) have attracted wide attention due to their ability to seamlessly embed the learning process with physical laws and their considerable success in solving forward and inverse differential equation (DE) problems. While most studies are improving the learning process and network architecture of PINNs, less attention has been paid to the modification of the DE system, which may play an important role in addressing some limitations of PINNs. One of the simplest modifications that can be implemented to all DE systems is the variable linear transformation (VLT). Therefore, in this work, we propose the VLT-PINNs that solve the DE systems of the linear-transformed variables instead of the original ones. To clearly illustrate the importance of prior knowledge in determining the VLT parameters, we choose the thin-layer flow problems as our focus. Ten related cases were tested, including the jet flows, wake flows, mixing layers, boundary layers and Kovasznay flows. Based on the principle of normalization and for a better match of the DE system to the preference of NNs, we identify three principles for determining the VLT parameters: magnitude normalization for dependent variables (principle 1), local normalization for independent variables (principle 2), and appropriate scaling for physics-related parameters in inverse problems (principle 3). The VLT-PINNs with the VLT parameters suggested by the proposed principles show excellent performance over all the test cases, while the results are quite poor with the VLT parameters suggested by traditional linear transformations, such as nondimensionalization and global normalization. Comparison studies also show that only under the constraints of the VLT principles can we obtain satisfactory results. Besides, we find tanh is more appropriate as the activation function than sin for thin-layer flow problems, from both posteriori results and priori analyses with physical intuition. We highlight that our VLT method is an attempt to combine the three advantages of accuracy, universality and simplicity, and hope that it can provide new insights into the better integration of prior knowledge, physical intuition and the nature of NNs. The code for this paper is available on https://github.com/CAME-THU/VLT-PINN. [ABSTRACT FROM AUTHOR]

Projection methods are applicable in many fields. It is a natural and practical approach to devise the invariant-preserving schemes for conservative systems. The idea is to project the solution of any underlying numerical scheme onto the manifold determined by the invariant, and this process will be referred to as the projection integrator. Generally, the projection integrator chooses the gradient of invariant as its projection direction and has the same order as the underlying method. In this paper, we propose a different projection direction to construct a new projection integrator whose order is higher than the underlying method. According to this novel direction, we further summarize high-order projection integrators with superconvergence and rigorously prove the truncation error by utilizing the linear integral method as a central tool. Apart from the invariant-preserving property, symmetry is an important geometric property for reversible differential equations. The design and analysis of another high-order projection integrators with symmetry and superconvergence are also presented in this paper. Numerical experiments are provided to verify our theoretical results and illustrate that our proposed projection integrators have superior behaviors in a long time numerical simulation. [ABSTRACT FROM AUTHOR]

The aim of this paper is to obtain results on topological properties of flows for nonnegative time in the framework of generalized ODEs. We define the concept of generalized semiflow and we present some recursive properties as minimality and recurrence. Using correspondence theorems, we translate our results to measure differential equations which are known to be special cases of generalized ODEs. Some examples are presented to illustrate the theory. [ABSTRACT FROM AUTHOR]

In [3] André showed that the minimal differential equations of Э-functions and E -functions have a basis of holomorphic solutions at every point except for zero and infinity. In addition, in [5] he observed that if f 1 (z) is an E -function with rational coefficients such that f 1 (1) = 0 , then f 1 (z) / (z − 1) is also an E -function. With this additional result André derived transcendence results for the values of E -functions. These results were further applied by Beukers [7] to obtain a strengthened Siegel-Shidlovskii theorem for E -functions. The arguments of André and Beukers should aid in obtaining a Siegel-Shidlovskii type theorem for Э-functions provided that one can show that if f − 1 (z) is a Э-function with rational coefficients such that the 1-summation F − 1 (z) vanishes at z = 1 , then f − 1 (z) / (z − 1) is also a Э-function. In this paper, we investigate this problem and derive several criteria to check its validity. [ABSTRACT FROM AUTHOR]

This paper is devoted to studying the averaging principle for a fast-slow system of rough differential equations driven by mixed fractional Brownian rough path. The fast component is driven by Brownian motion, while the slow component is driven by fractional Brownian motion with Hurst index H (1 / 3 < H ≤ 1 / 2). Combining the fractional calculus approach to rough path theory and Khasminskii's classical time discretization method, we prove that the slow component strongly converges to the solution of the corresponding averaged equation in the L 1 -sense. The averaging principle for a fast-slow system in the framework of rough path theory seems new. [ABSTRACT FROM AUTHOR]

In this paper, we are interested in investigating stability results for generalized ordinary differential equations (generalized ODEs in short), and their applications to measure differential equations and dynamic equations on time scales. First, we establish stability, asymptotic and exponential stability for the trivial solution of generalized ODEs. Secondly, we use the well known correspondence between solutions of generalized ODEs and solutions of measure differential equations, obtaining analogues results for the last equations. Finally, we apply some of these results for dynamic equations on time scales. [ABSTRACT FROM AUTHOR]

This paper considers compound operators on a general vector space. Explicit representations for these operators are obtained in various contexts with emphasis on the compound differential equations corresponding to a given linear differential equation. Applications to estimation of eigenvalues of self-adjoint boundary value problems and of the codimension of the stable manifold of a linear differential equation in a Banach space are given. [ABSTRACT FROM AUTHOR]

In this paper we aim to find the highest number of critical periods in a class of planar systems of polynomial differential equations for fixed degree having a center. We fix our attention to lower bounds of local criticality for low degree planar polynomial centers. The main technique is the study of perturbations of reversible holomorphic (isochronous) centers, inside the reversible centers class. More concretely, we study the Taylor developments of the period constants with respect to the perturbation parameters. First, we see that there are systems of degree 3 ≤ n ≤ 16 for which up to first order at least (n 2 + n − 4) / 2 critical periods bifurcate from the center. Second, we improve this number for centers with degree from 3 to 9. In particular, we obtain 6 and 10 critical periods for cubic and quartic degree systems, respectively. [ABSTRACT FROM AUTHOR]

Orthogonal and biorthogonal (multi)wavelets on the real line have been extensively studied and employed in applications with success. On the other hand, a lot of problems in applications such as images and solutions of differential equations are defined on bounded intervals or domains. Therefore, it is important in both theory and application to construct all possible wavelets on intervals with some desired properties from (bi)orthogonal (multi)wavelets on the real line. Then wavelets on rectangular domains such as [ 0 , 1 ] d can be obtained through tensor product. Vanishing moments of compactly supported wavelets are the key property for sparse wavelet representations and are closely linked to polynomial reproduction of their underlying refinable (vector) functions. Boundary wavelets with low order vanishing moments often lead to undesired boundary artifacts as well as reduced sparsity and approximation orders near boundaries in applications. Scalar orthogonal wavelets and spline biorthogonal wavelets on the interval [ 0 , 1 ] have been extensively studied in the literature. Though multiwavelets enjoy some desired properties over scalar wavelets such as high vanishing moments and relatively short support, except a few concrete examples, there is currently no systematic method for constructing (bi)orthogonal multiwavelets on bounded intervals. In contrast to current literature on constructing particular wavelets on intervals from special (bi)orthogonal (multi)wavelets, from any arbitrarily given compactly supported (bi)orthogonal multiwavelet on the real line, in this paper we propose two different approaches to construct/derive all possible locally supported (bi)orthogonal (multi)wavelets on [ 0 , ∞) or [ 0 , 1 ] with or without prescribed vanishing moments, polynomial reproduction, and/or homogeneous boundary conditions. The first approach generalizes the classical approach from scalar wavelets to multiwavelets, while the second approach is direct without explicitly involving any dual refinable functions and dual multiwavelets. We shall also address wavelets on intervals satisfying general homogeneous boundary conditions. Though constructing orthogonal (multi)wavelets on intervals is much easier than their biorthogonal counterparts, we show that some boundary orthogonal wavelets cannot have any vanishing moments if these orthogonal (multi)wavelets on intervals satisfy the homogeneous Dirichlet boundary condition. In comparison with the classical approach, our proposed direct approach makes the construction of all possible locally supported (multi)wavelets on intervals easy. Seven examples of orthogonal and biorthogonal multiwavelets on the interval [ 0 , 1 ] will be provided to illustrate our construction approaches and proposed algorithms. [ABSTRACT FROM AUTHOR]

In this paper we consider a second order nonlinear hyperbolic equation of Monge-Ampère type arizing in the Light Curve Inversion problem [7,8,4] when dealing with objects with regions with negative gaussian curvature. The above works imply that such regions are of practical importance in the Light Curve Inversion problem. It is shown here that in such regions existence and uniqueness of solutions to the original equation can be obtained by solving a system of four first order quasi-linear differential equations supplemented by suitable initial and boundary conditions. Furthermore, these solutions can be found by integrating a system of ordinary differential equations. The numerical methods in [14] can be applied to actually find these solutions. We present here the completed part of the overall project, which establishes the necessary basis for efficient numerical solution. It is expected that the part of the project dealing with numerics will be addressed in a separate paper. This paper is written in a form accessible to physicists, engineers and applied mathematicians working on the Light Curve Inversion problem and design of reflectors. [ABSTRACT FROM AUTHOR]

Mathematical models play an increasing role in understanding and predicting machining processes, in particular milling. However, despite the considerable efforts that have been dedicated to this problem, a majority of milling models still rely on simplifying assumptions to calculate the chip thickness. In this paper, the chip thickness is determined without these simplifications, based on a surface function that describes the milled surface and on information about the workpiece boundary. By combining the partial differential equation (PDE) governing the evolution of this surface function with the ordinary differential equations (ODE) governing the tool/machine dynamics, a mixed PDE–ODE formulation is proposed to describe the dynamics of the milling process. The coupled system of differential equations is solved using an algorithm that combines finite difference (ODE) and finite volume (PDE) methods. A case study is presented to compare the proposed approach with the classical delay differential equations (DDE) model formulation for milling processes based on a simplified chip thickness model. The PDE–ODE formulation represents an explicit mathematical model for milling process dynamics; it yields a theoretically exact chip thickness and offers a means to assess the validity of models based on DDE formulation. Moreover, the proposed formulation is capable of simulating transient tool behaviors when the tool is milling the outer region of the workpiece, which is in general neglected by the DDE-based models. • Accurate mathematical description of the milling process. • Evolution of the machined surface around the tool is described. • Chip thickness model affects the accuracy of milling simulation in certain scenarios. • Transient tool behaviors are captured while milling the outer part of the workpiece. [ABSTRACT FROM AUTHOR]

• Dynamical model used for long-distance vibration conveyors is proposed. • The four exciters split-driving synchronization theory scheme is constructed. • Synchronization, stability and Sommerfeld effect are analyzed. • Machine functions in the different resonant regions are discussed. • Experimental prototype is designed to verify the correctness of the theory. The previous studies on the vibrating system with multiple exciters less considered the Sommerfeld effect, which were mainly focused on synchronization and stability problems of the system in the whole resonant region. A dynamical model is proposed in this paper to study the synchronization, stability and the Sommerfeld effect of four exciters in the multi-resonant types vibrating system (MRTVS). In order to improve the power of the system, the vibrating system with two rigid frames is driven by four especially distributed exciters. The motion differential equations of the system are given by Lagrange's equations, and the dimensionless coupling equations and synchronization criterion are constructed as well by the average method. The theory condition for stability of four exciters in synchronous states is derived from the Hamilton principle. The coupling dynamic characteristics, stability and Sommerfeld effect of the system are discussed numerically. The whole resonant region of the system is divided into three segments by natural frequencies, and the synchronous and stable states in the different resonant types are analyzed respectively. The diversity phenomenon of the nonlinear system is observed, in other words, the system exists multiple stable equilibrium solutions at a certain particular resonant region. Then the Sommerfeld effect around the natural frequencies (NFs) is further revealed. The correctness of theoretical analyses is examined by simulation and the experiment. It's shown that the reasonable working point in engineering should be selected in the sub-resonant region with respect to the natural frequency of the main vibrating system. The present work can provide theoretical guidance for designing some new types of vibrating machines with high driving power. [ABSTRACT FROM AUTHOR]

The aim of this paper is to develop an analytical approach to obtain the sharp estimates for the lowest positive periodic eigenvalue and all Dirichlet eigenvalues of a general Sturm-Liouville problem y ″ = q (t) y + λ m (t) y , where q is a nonnegative potential and another potential m admits to change sign. A typical example of such problems is the well-known Camassa-Holm equations with indefinite potentials, which corresponds to the case q (t) ≡ 1 4. It is shown that the solution of the minimization problems of the lowest positive periodic eigenvalues and Dirichlet eigenvalues will lead to more general distributions of potentials which have no densities with respect to the Lebesgue measure. As a result, it is very natural to choose the general setting of the measure differential equations d y • = y (t) d μ (t) + λ y d ν (t) , to understand the eigenvalues and their minimization, where μ and ν are two suitable measures. The variational characterization of lowest positive eigenvalues, together with a strong continuous dependence of eigenvalues on the potentials, will play crucial roles in our analysis. Different from the periodic case, we are able to obtain the optimal bounds for higher Dirichlet eigenvalues. [ABSTRACT FROM AUTHOR]

Summation-by-parts (SBP) operators allow us to systematically develop energy-stable and high-order accurate numerical methods for time-dependent differential equations. Until recently, the main idea behind existing SBP operators was that polynomials can accurately approximate the solution, and SBP operators should thus be exact for them. However, polynomials do not provide the best approximation for some problems, with other approximation spaces being more appropriate. We recently addressed this issue and developed a theory for one-dimensional SBP operators based on general function spaces, coined function-space SBP (FSBP) operators. In this paper, we extend the theory of FSBP operators to multiple dimensions. We focus on their existence, connection to quadratures, construction, and mimetic properties. A more exhaustive numerical demonstration of multi-dimensional FSBP (MFSBP) operators and their application will be provided in future works. Similar to the one-dimensional case, we demonstrate that most of the established results for polynomial-based multi-dimensional SBP (MSBP) operators carry over to the more general class of MFSBP operators. Our findings imply that the concept of SBP operators can be applied to a significantly larger class of methods than is currently done. This can increase the accuracy of the numerical solutions and/or provide stability to the methods. • Multi-dimensional summation-by-parts operators for general (non-polynomial) function spaces. • Existence proof that connects them to quadratures. • Construction procedure for general multi-dimensional domains/elements/blocks. • Proof of discrete conservation and energy-boundedness. • Numerical comparison demonstrating that non-polynomial function spaces can be of advantage. [ABSTRACT FROM AUTHOR]

In this paper, we propose a cusp-capturing physics-informed neural network (PINN) to solve discontinuous-coefficient elliptic interface problems whose solution is continuous but has discontinuous first derivatives on the interface. To find such a solution using neural network representation, we introduce a cusp-enforced level set function as an additional feature input to the network to retain the inherent solution properties; that is, capturing the solution cusps (where the derivatives are discontinuous) sharply. In addition, the proposed neural network has the advantage of being mesh-free, so it can easily handle problems in irregular domains. We train the network using the physics-informed framework in which the loss function comprises the residual of the differential equation together with certain interface and boundary conditions. We conduct a series of numerical experiments to demonstrate the effectiveness of the cusp-capturing technique and the accuracy of the present network model. Numerical results show that even using a one-hidden-layer (shallow) network with a moderate number of neurons and sufficient training data points, the present network model can achieve prediction accuracy comparable with traditional methods. Besides, if the solution is discontinuous across the interface, we can simply incorporate an additional supervised learning task for solution jump approximation into the present network without much difficulty. • Cusp-capturing physics informed neural network for elliptic interface problems. • Our network can present continuous solutions that inherently have discontinuous first derivatives on interfaces. • A mesh-free approach for solving PDEs with interfaces is presented. • Our network can achieve high prediction accuracy with large contrast discontinuous coefficients. [ABSTRACT FROM AUTHOR]

The paper deals with a suitable transformation related to the metric scalar of the hyper-surface so that the Raychaudhuri Equation (RE) can be written as a second order nonlinear differential equation. A first integral of this second order differential equation gives a possible analytic solution of the RE. Also, it is shown that construction of a Lagrangian (and hence a Hamiltonian) is possible, from which the RE can be derived. Wheeler-Dewitt equation has been formulated in canonical quantization scheme and norm of its solution (wave function of the universe) is shown to affect the singularity analysis in the quantum regime for any spatially homogeneous and isotropic cosmology. Finally Bohmian trajectories are formulated with causal interpretation and these quantum trajectories unlike classical geodesics obliterate the initial big-bang singularity when the quantum potential is included. [ABSTRACT FROM AUTHOR]

Aiming at the problem that the dynamic displacement at critical positions of the cable cannot be measured due to economic or difficult installation, this paper proposes a method to invert the dynamic displacement of the cable in the whole field by using the dynamic displacement of a single measuring point. The method can realize the inversion of the lateral dynamic displacement response (L-DDR) and angular dynamic displacement response (A-DDR). Firstly, the exact mode shape functions of the L-DDR and A-DDR, considering bending stiffness, sag, inclination angle, and boundary conditions, are derived from the cable vibration control differential equation. The exact mode shape functions are a closed-form solution; thus, it can accurately calculate the cable mode shape at any frequency in real time. Then, an inversion method of the dynamic displacement of the cable in the whole field is proposed based on the amplitude ratio of the exact mode shape. This method can be used to invert the dynamic displacement of the cable in the frequency domain and time domain. Finally, through a numerical simulation and an actual cable test, it is verified that the accuracy of the proposed method is satisfactory. The cable dynamic displacement components at modal and external load frequencies are the main contributors to the dynamic displacement of the cable. [ABSTRACT FROM AUTHOR]

• This paper explains that the dislocation of the cross-sectional centroids is the main reason for the longitudinal–transverse coupled vibration of the variable-height asymmetric beam. • The main mechanical characteristics of variable-height asymmetric beams are analyzed in detail. • The dynamic model of the variable-height asymmetric beam is proposed and verified by experiments. • This paper demonstrates that a single excitation source can excite vibrations in two directions of a variable-height asymmetric beam. In general, the cross-sectional centroids of the variable-height asymmetric beam change in the height direction. The Timoshenko beam theory does not consider the influence of the varying cross-sectional centroids and cannot accurately reveal the vibration of variable-height asymmetric beams. Considering that the cross-sectional centroids are dislocated, we propose a novel dynamic model for the variable-height asymmetric beam element. The proposed dynamic model is based on the rigid-plane assumption and can reflect its main mechanical characteristics. The non-collinear axial forces produce the complementary moment. The rotational radiuses of the two ends of longitudinal fibers are not equal. After retaining the main mechanical characteristics, the variable-height element is roughly treated as a constant-height element to simplify the governing differential equations. A bonded-type piezoelectric transducer with variable-height asymmetric beams is taken as an example to verify the developed dynamic model. The theoretical and experimental resonant frequencies are 39338 Hz and 38704 Hz, respectively, with a difference of 634 Hz and a relative error of 1.638%. The calculated displacement vibration shapes are consistent with the experimental results. Finally, the modeling differences between variable-height symmetric and asymmetric beams are clarified from three aspects. [ABSTRACT FROM AUTHOR]

In this paper, a closed-form solution of the moving harmonic load problem for continuous Euler–Bernoulli beam systems is presented. The generality of the boundary conditions is taken into account by solving the characteristic equation of the system, obtaining its natural frequencies and mode shapes. The undetermined coefficient method is applied to solve the governing differential equation of motion, determining the base functions of the solution space of the problem. For vertical vibrations, three numerical examples of footbridges are presented. The main contribution of this paper is to provide the closed-form solution of the moving harmonic load problem applied to continuous footbridges including the phase angle in the load definition. In this way, it is possible to find the solution in the time domain of the harmonic component of any load spectra. • A closed-form solution for the moving harmonic load problem is presented. • The solution includes the phase angle in the load definition. • The solution can be applied to a general continuous beam with flexible supports. • The solution can be used in the time domain and for any load spectra. • Great agreement was found between FE-models and the closed-form solution. [ABSTRACT FROM AUTHOR]

• New high-order energy-preserving integrators are developed for Poisson systems. • Theoretical analysis including the existence and uniqueness, orders, implementation issues, and conservation of Casimir function is made. • The new integrators show much better numerical behaviour in long-term computation compared with the existing integrators. As is known, a system of differential equations possessing a first integral can be written in a non-canonical Hamiltonian system that often leads to a Poisson system under moderate conditions. Therefore, energy-preserving numerical integrators for such systems have become a subject of intensive investigation. In this paper, we design and analyse two classes of energy-preserving functionally-fitted integrators of arbitrary order for Poisson systems with highly oscillatory solutions using a new framework. We also show that both the order and the stage order of the proposed integrators may be affected by the used quadrature formula in practice. Furthermore, the existence and uniqueness of energy-preserving functionally-fitted integrators, their implementation issues, and the conservation of Casimir functions are investigated in detail. Finally, numerical experiments are carried out with the fourth-order and sixth-order energy-preserving integrators proposed in this paper, including the Lotka-Volterra system and the charged-particle dynamics in a strong magnetic field, and the numerical results demonstrate the remarkable accuracy and efficiency of our high-order energy-preserving integrators compared with the other energy-preserving methods in the literature. [ABSTRACT FROM AUTHOR]

The concept of physics-informed neural networks has become a useful tool for solving differential equations due to its flexibility. A few approaches use this concept to solve the eikonal equation that describes the first-arrival traveltimes of waves propagating in smooth heterogeneous velocity models. However, the challenge of the eikonal is exacerbated by the velocity models producing caustics, resulting in instabilities and deterioration of accuracy due to the non-smooth solution behavior. In this paper, we revisit the problem of solving the eikonal equation using neural networks to tackle caustic pathologies. We introduce the novel Neural Eikonal Solver (NES) for solving the isotropic eikonal equation in two formulations: the one-point problem is for a fixed source location; the two-point problem is for an arbitrary source-receiver pair. We present several techniques which provide stability in the case of caustics: improved factorization; non-symmetric loss function based on Hamiltonian; gaussian activation; symmetrization. In our tests, NES showed the relative mean-absolute error of 0.2-0.4% from the second-order factored Fast Marching Method with a similar inference time, and outperformed existing neural-network solvers giving 10-60 times lower errors and 2-30 times faster training. The one-point NES provides the most accurate solution, while the two-point NES gives slightly lower accuracy but an extremely compact representation with all spatial derivatives. It can be useful in many seismic problems: massive computations of traveltimes for millions of source-receiver pairs in Kirchhoff migration; modeling of ray amplitudes using spatial derivatives; traveltime tomography; earthquake localization; ray multipathing analysis. • Eikonal equation describes wave propagation in heterogeneous models. • Physics-informed neural networks can approximate the eikonal solution. • Physics-based modifications of the neural-network architecture enhance the accuracy. • Neural network compresses the solution, useful in large computational problems. • Neural network gives spatial derivatives, useful for calculation of wave amplitudes. [ABSTRACT FROM AUTHOR]

In this paper a comprehensive and structured analysis of the stabilizing properties linked to the combination of different chatter suppression techniques in milling was carried out. Specifically, the spindle speed variation, the stiffness variation and the usage of tools with variable pitch and variable helixangles were considered. For this purpose, a unique, general and efficient frequency domain formulation (extended harmonic solution) was developed. It is suitable for dealing with differential delayed equations characterized by multiple and distributed time dependent delays and time periodic dynamics. It was demonstrated that the developed solution converges to the reference cases, which were studied in literature through semi-discretization and its extensions. Moreover, by selecting a proper number of considered harmonics, the computation time was reduced of one order of magnitude, in the face of a slight reduction of accuracy with respect to time domain approaches. The computed stability lobe diagrams demonstrated that the combination of multi-chatter suppression techniques allows exploiting the best features of each methodology. The simultaneous use of the spindle speed variation and the stiffness variation was studied for the first time. The absolute limit of stability increased up to 204%, much more than by adopting the techniques in a separate manner. If the enlargement of the stability region is considered as the reference target, the best results were granted by the combination of all the three strategies. Indeed, the improvement with respect to the use of only two techniques was at least 10%. • Stability of multiple point, distributed, time varying delay differential equation. • Frequency domain formulation for complex milling applications. • Investigation of novel chatter suppression technique combinations. • Improvement of chatter stability with respect to previous results in literature. [ABSTRACT FROM AUTHOR]

• Eigenproperties of a multi-cracked Euler-Bernoulli inextensible circular arch. • Cracks are accounted for in governing equation by means of the distribution theory. • Governing equation defined over a unique integration domain regardless n° cracks. • Integration via Laplace transform is adopted to infer closed form mode shapes. • Both validations and parametric analyses are presented. Despite the numerous explicit solutions of free vibration of arches with regular cross sections, in case of concentrated defects such as cracks, no procedure is available to analyse arch vibrations without sub-division of the integration domain. As a result, curved sub-elements comprised between crack and external constraints, or successive cracks, are considered. In this paper a distributional approach is adopted to provide a formulation of the free vibration differential governing equations of circular inextensible arches over a unique integration domain in presence of multiple concentrated open (non-breathing) cracks. Discontinuities due to the presence of an arbitrary number of cracks are modelled by means of Dirac's deltas. An integration procedure is devised to offer closed form solutions of the relevant vibration modes together with the relevant frequency determinantal equation. Natural frequencies and mode shapes of damaged arches with different damage and restraint configurations have been evaluated and compared with experimental results available in the literature as well as finite element numerical simulations. The presented closed form solutions are also employed for two parametric studies to evaluate the influence of an increasing number of along axis concentrated cracks as well as of the location of cracks along the arch span. [ABSTRACT FROM AUTHOR]

The propagation of sound waves in a horizontally stratified environment, a classic problem in ocean acoustics, has traditionally been calculated using normal modes. Most programs based on the normal mode model are discretized using the finite difference method (FDM). In this paper, a Legendre collocation method (LCM) based on domain decomposition is proposed to solve this problem. A set of collocation points cannot penetrate multiple layers of media, thus necessitating domain decomposition and the use of multiple sets of collocation points. The solution process of this method proceeds entirely in physical space, requiring that the original differential equation be strictly established at the collocation points; thus, a dense matrix eigenvalue system is formed, from which the solution for the horizontal wavenumbers and modes can be directly obtained. Numerical experiments are presented to demonstrate the validity and applicability of this method. A comparison with other methods shows that the LCM proposed in this article is more accurate than the FDM and offers approximately the same accuracy as, but a faster calculation speed than, other types of spectral methods. • This paper proposes a new method to calculate acoustic propagation in horizontally stratified marine environment. • The accuracy of the new method is greater than that of the classical finite difference method. • The new method is far faster than the Legendre-Galerkin and Chebyshev-Tau spectral methods. [ABSTRACT FROM AUTHOR]

The principal aim of this paper is to establish an iteration method on the space of resurgent functions. We discuss endless continuability of iterated convolution products of resurgent functions associated with rooted trees and derive their estimates. We show the resurgence of formal series solutions of nonlinear differential and difference equations by applying the estimates. [ABSTRACT FROM AUTHOR]

Mixed-precision algorithms combine low- and high-precision computations in order to benefit from the performance gains of reduced-precision without sacrificing accuracy. In this work, we design mixed-precision Runge–Kutta–Chebyshev (RKC) methods, where high precision is used for accuracy, and low precision for stability. Generally speaking, RKC methods are low-order explicit schemes with a stability domain growing quadratically with the number of function evaluations. For this reason, most of the computational effort is spent on stability rather than accuracy purposes. In this paper, we show that a naïve mixed-precision implementation of any Runge–Kutta scheme can harm the convergence order of the method and limit its accuracy, and we introduce a new class of mixed-precision RKC schemes that are instead unaffected by this limiting behavior. We present three mixed-precision schemes: a first- and a second-order RKC method, and a first-order multirate RKC scheme for multiscale problems. These schemes perform only the few function evaluations needed for accuracy (1 or 2 for first- and second-order methods respectively) in high precision, while the rest are performed in low precision. We prove that while these methods are essentially as cheap as their fully low-precision equivalent, they retain the stability and convergence order of their high-precision counterpart. Indeed, numerical experiments confirm that these schemes are as accurate as the corresponding high-precision method. • Mixed-precision algorithms can obtain high-accuracy results at a low-precision cost. • A naïve mixed-precision implementation can affect the convergence of ODE solvers. • A few high-precision function evaluations are sufficient to retain full accuracy. • Our mixed-precision Runge–Kutta–Chebyshev methods: fast and high-precision accurate. • We consider explicit stabilized methods for both single-scale and multirate problems. [ABSTRACT FROM AUTHOR]

• New analysis of the ligand-induced dimerisation process for VEGFR is presented. • A key idea in the mathematical model and analysis is kinetic binding cooperativity. • The nonlinear model is in contrast with the corresponding linear model for GPCRs. • Numerical and asymptotic analysis reveal insights into the underlying LID dynamics. • Our simulations show excellent fit to recent experimental data for VEGFR. The vascular endothelial growth factor (VEGF) receptor (VEGFR) system plays a role in cancer and many other diseases. It is widely accepted that VEGFR receptors dimerise in response to VEGF binding. However, analysis of these mechanisms and their implications for drug development still requires further exploration. In this paper, we present a mathematical model representing the binding of VEGF to VEGFR and the subsequent ligand-induced dimerisation. A key factor in this work is the qualitative and quantitative effect of binding cooperativity, which describes the effect that the binding of a ligand to a receptor has on the binding of that ligand to a second receptor, and the dimerisation of these receptors. We analyse the ordinary differential equation system at equilibrium, giving analytical solutions for the total amount of ligand bound. For time-course dynamics, we use numerical methods to explore possible behaviours under various parameter regimes, while perturbation analysis is used to understand the intricacies of these behaviours. Our simulation results show an excellent fit to experimental data, towards validating the model. [ABSTRACT FROM AUTHOR]

• On the stability boundary the mass-beam-foundation system is in the form of a steady state periodic motion. • The infinite beam is divided into three segments, in each of which the differential equation is solved analytically. • The point mass model overestimates the stability of the mass-beam-foundation system. The stability of a distributed mass of finite length travelling uniformly on an infinite beam resting on a visco-elastic foundation is considered in this paper. The focus is placed on the effects of the length of the distributed mass on the stability of the mass-beam-foundation system. It is argued that on the stability boundary, the mass-beam-foundation system is in the form of a steady state periodic motion. To seek the periodic solution, the infinite beam is divided into three segments, in each of which the differential equation is solved analytically. By comparing the stability boundary curves of the distributed mass of finite length with the one of infinite length, it is found that the stability boundary curves of finite length do not converge to the one of infinite length. It is believed that the course of the failure to converge towards the infinite length model is the differences in the boundary conditions at infinity. On the other hand, the stability boundary curves of finite length do converge to the point mass model when the mass length approaches zero. Compared to the point mass model, the distributed mass model is less stable. More precisely, for a specified system damping and total mass, the critical speed is in general much smaller if distributed model is adopted. The longer the mass length, the smaller the critical speed is. In short, the conventional point mass model overestimates the stability of the mass-beam-foundation system. [ABSTRACT FROM AUTHOR]