{"title":"Preface: Special issue on continuation methods and applications","authors":"B. Krauskopf, H. Osinga","doi":"10.3934/jcd.2022015","DOIUrl":"https://doi.org/10.3934/jcd.2022015","url":null,"abstract":"<jats:p xml:lang=\"fr\" />","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"19 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82909485","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Smooth complete Lyapunov functions for multifunctions","authors":"S. Suhr","doi":"10.3934/jcd.2022025","DOIUrl":"https://doi.org/10.3934/jcd.2022025","url":null,"abstract":"<p style='text-indent:20px;'>We prove existence of smooth complete Lyapunov functions for convex semicontinuous multifunctions.</p>","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"13 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78265567","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Emergence of quasiperiodic regimes in a neutral delay model of flute-like instruments: Influence of the detuning between resonance frequencies","authors":"S. Terrien, C. Vergez, B. Fabre, P. de la Cuadra","doi":"10.3934/jcd.2022011","DOIUrl":"https://doi.org/10.3934/jcd.2022011","url":null,"abstract":"Musical instruments display a wealth of dynamics, from equilibria (where no sound is produced) to a wide diversity of periodic and non-periodic sound regimes. We focus here on two types of flute-like instruments, namely a recorder and a pre-hispanic Chilean flute. A recent experimental study showed that they both produce quasiperiodic sound regimes which are avoided or played on purpose depending on the instrument. We investigate the generic model of sound production in flute-like musical instruments, a system of neutral delay-differential equations. Using time-domain simulations, we show that it produces stable quasiperiodic oscillations in good agreement with experimental observations. A numerical bifurcation analysis is performed, where both the delay time (related to a control parameter) and the detuning between the resonance frequencies of the instrument – a key parameter for instrument makers – are considered as bifurcation parameters. This demonstrates that the large detuning that is characteristic of prehispanic Chilean flutes plays a crucial role in the emergence of stable quasiperiodic oscillations.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"26 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88327269","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Monotonicity and discretization of Hammerstein integrodifference equations","authors":"Magdalena Nockowska-Rosiak, C. Pötzsche","doi":"10.3934/jcd.2022023","DOIUrl":"https://doi.org/10.3934/jcd.2022023","url":null,"abstract":"The paper provides sufficient conditions for monotonicity, subhomogeneity and concavity of vector-valued Hammerstein integral operators over compact domains, as well as for the persistence of these properties under numerical discretizations of degenerate kernel type. This has immediate consequences on the dynamics of Hammerstein integrodifference equations and allows to deduce a local-global stability principle.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"199 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86833897","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Construction of mean-square Lyapunov-basins for random ordinary differential equations","authors":"Florian Rupp","doi":"10.3934/jcd.2022024","DOIUrl":"https://doi.org/10.3934/jcd.2022024","url":null,"abstract":"We propose a straightforward basin search algorithm to determine a suitably large level set of the mean-square Lyapunov-function that corresponds to the linearization about an path-wise equilibrium solution of a random ordinary differential equation (RODE). Noise intensity plays a crucial role for how similar the behavior of solutions of RODEs is compared to the corresponding deterministic system. In this regards, the basin search algorithm also allows to numerically estimate up to which noise intensities linearized mean-square asymptotic stability remains.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"30 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88764458","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
E. Celledoni, C. Evripidou, D. McLaren, B. Owren, G. Quispel, B. Tapley
{"title":"Detecting and determining preserved measures and integrals of birational maps","authors":"E. Celledoni, C. Evripidou, D. McLaren, B. Owren, G. Quispel, B. Tapley","doi":"10.3934/jcd.2022014","DOIUrl":"https://doi.org/10.3934/jcd.2022014","url":null,"abstract":"In this paper we use the method of discrete Darboux polynomials to calculate preserved measures and integrals of rational maps. The approach is based on the use of cofactors and Darboux polynomials and relies on the use of symbolic algebra tools. Given sufficient computing power, most, if not all, rational preserved integrals can be found (and even some non-rational ones). We show, in a number of examples, how it is possible to use this method to both determine and detect preserved measures and integrals of the considered rational maps, thus lending weight to a previous ansatz [2]. Many of the examples arise from the Kahan-Hirota-Kimura discretization of completely integrable systems of ordinary differential equations.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"104 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74484878","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Modified refinement algorithm to construct Lyapunov functions using meshless collocation","authors":"N. Mohammed, P. Giesl","doi":"10.3934/jcd.2022022","DOIUrl":"https://doi.org/10.3934/jcd.2022022","url":null,"abstract":"Lyapunov functions are functions with negative derivative along solutions of a given ordinary differential equation. Moreover, sublevel sets of a Lyapunov function are subsets of the domain of attraction of the equilibrium. One of the numerical construction methods for Lyapunov functions uses meshless collocation with radial basis functions.Recently, this method was combined with a grid refinement algorithm (GRA) to reduce the number of collocation points needed to construct Lyapunov functions. However, depending on the choice of the initial set of collocation point, the algorithm can terminate, failing to compute a Lyapunov function. In this paper, we propose a modified grid refinement algorithm (MGRA), which overcomes these shortcomings by adding appropriate collocation points using a clustering algorithm. The modified algorithm is applied to two- and three-dimensional examples.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"15 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74609607","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Low-rank kernel approximation of Lyapunov functions using neural networks","authors":"K. Webster","doi":"10.3934/jcd.2022026","DOIUrl":"https://doi.org/10.3934/jcd.2022026","url":null,"abstract":"We study the use of single hidden layer neural networks for the approximation of Lyapunov functions in autonomous ordinary differential equations. In particular, we focus on the connection between this approach and that of the meshless collocation method using reproducing kernel Hilbert spaces. It is shown that under certain conditions, an optimised neural network is functionally equivalent to the RKHS generalised interpolant solution corresponding to a kernel function that is implicitly defined by the neural network. We demonstrate convergence of the neural network approximation using several numerical examples, and compare with approximations obtained by the meshless collocation method. Finally, motivated by our theoretical and numerical findings, we propose a new iterative algorithm for the approximation of Lyapunov functions using single hidden layer neural networks.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"13 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78822206","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A novel moving orthonormal coordinate-based approach for region of attraction analysis of limit cycles","authors":"Eva Ahbe, A. Iannelli, Roy S. Smith","doi":"10.3934/jcd.2022016","DOIUrl":"https://doi.org/10.3934/jcd.2022016","url":null,"abstract":"The paper proposes a Lyapunov theory-based method to compute inner estimates of the region of attraction (ROA) of stable limit cycles. The approach is based on a transformation of the system to transverse coordinates, defined on a moving orthonormal coordinate system (MOC) for which a novel construction is presented. The proposed center point MOC (cp-MOC) is associated with a user-defined center point and provides flexibility to the construction of the transverse coordinates. In particular, compared to the standard approach based on hyperplanes orthogonal to the flow, the new construction allows the analyst to obtain larger regions of the state space where the well-definedness property of the transformation is satisfied. This has important benefits when using transverse coordinates to compute inner estimates of the ROA. To demonstrate these improvements, a sum-of-squares optimization-based formulation is proposed for computing inner estimates of the ROA of limit cycles for polynomial dynamics described in transverse coordinates. Different algorithmic options are explored, taking into account computational and accuracy aspects. Results are shown for three different systems exhibiting increasing complexity. The presented algorithms are extensively compared, and the newly cp-MOC is shown to markedly outperform existing approaches.","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"75 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77335582","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Rotation invariant patterns for a nonlinear Laplace-Beltrami equation: A Taylor-Chebyshev series approach","authors":"J. B. Berg, Gabriel William Duchesne, J. Lessard","doi":"10.3934/jcd.2022005","DOIUrl":"https://doi.org/10.3934/jcd.2022005","url":null,"abstract":"<p style='text-indent:20px;'>In this paper, we introduce a rigorous computational approach to prove existence of rotation invariant patterns for a nonlinear Laplace-Beltrami equation posed on the 2-sphere. After changing to spherical coordinates, the problem becomes a singular second order boundary value problem (BVP) on the interval <inline-formula><tex-math id=\"M1\">begin{document}$ (0,frac{pi}{2}] $end{document}</tex-math></inline-formula> with a <i>removable</i> singularity at zero. The singularity is removed by solving the equation with Taylor series on <inline-formula><tex-math id=\"M2\">begin{document}$ (0,delta] $end{document}</tex-math></inline-formula> (with <inline-formula><tex-math id=\"M3\">begin{document}$ delta $end{document}</tex-math></inline-formula> small) while a Chebyshev series expansion is used to solve the problem on <inline-formula><tex-math id=\"M4\">begin{document}$ [delta,frac{pi}{2}] $end{document}</tex-math></inline-formula>. The two setups are incorporated in a larger zero-finding problem of the form <inline-formula><tex-math id=\"M5\">begin{document}$ F(a) = 0 $end{document}</tex-math></inline-formula> with <inline-formula><tex-math id=\"M6\">begin{document}$ a $end{document}</tex-math></inline-formula> containing the coefficients of the Taylor and Chebyshev series. The problem <inline-formula><tex-math id=\"M7\">begin{document}$ F = 0 $end{document}</tex-math></inline-formula> is solved rigorously using a Newton-Kantorovich argument.</p>","PeriodicalId":37526,"journal":{"name":"Journal of Computational Dynamics","volume":"97 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89109291","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}