Singular Perturbation Research Articles

Superconvergence analysis of the conforming discontinuous Galerkin method on a Bakhvalov-type mesh for singularly perturbed reaction–diffusion equation

The conforming discontinuous Galerkin (CDG) method maximizes the utilization of all degrees of freedom of the discontinuous Pk polynomial to achieve a convergence rate two orders higher than its counterpart conforming finite element method employing continuous Pk element. Despite this superiority, there is little theory of the CDG methods for singular perturbation problems. In this paper, superconvergence of the CDG method is studied on a Bakhvalov-type mesh for a singularly perturbed reaction–diffusion problem. For this goal, a pre-existing least squares method has been utilized to ensure better approximation properties of the projection. On the basis of that, we derive superconvergence results for the CDG finite element solution in the energy norm and L2-norm and obtain uniform convergence of the CDG method for the first time.

Existence of weak solutions for a class of (p(u),q(u))-Laplacian problems

This paper investigates the existence of weak solutions for (p,q)-Laplacian problems where p and q depend on the unknown solution. We focus on the case where p and q are local quantities. By means of a singular perturbation technique, we prove the existence of weak solutions for certain Dirichlet problem.

Stability Analysis of Some Types of Singularly Perturbed Time-Delay Differential Systems: Symmetric Matrix Riccati Equation Approach

Several types of linear and nonlinear singularly perturbed time-delay differential systems are considered. Asymptotic stability of the linear systems and asymptotic stability of the trivial solution of the nonlinear systems, valid for any sufficiently small value of the parameter of singular perturbation, are analyzed. For the stability analysis in the linear case, a partial exact slow–fast decomposition of the original system and an application of the Symmetric Matrix Riccati Equation method are proposed. Such an analysis yields parameter-free conditions, providing the asymptotic stability of the considered linear singularly perturbed time-delay differential systems for any sufficiently small value of the parameter of singular perturbation. Using the asymptotic stability results for the considered linear systems and the method of asymptotic stability in the first approximation, parameter-free conditions, guaranteeing the asymptotic stability of the trivial solution to the considered nonlinear systems for any sufficiently small value of the parameter of singular perturbation, are derived. Illustrative examples are presented.

On the Scaling of the Cubic-to-Tetragonal Phase Transformation with Displacement Boundary Conditions

We provide (upper and lower) scaling bounds for a singular perturbation model for the cubic-to-tetragonal phase transformation with (partial) displacement boundary data. We illustrate that the order of lamination of the affine displacement data determines the complexity of the microstructure. As in (Rüland and Tribuzio in ESAIM Control Optim. Calc. Var. 29:68, 2023) we heavily exploit careful Fourier space localization methods in distinguishing between the different lamination orders in the data.

Nonlinear Observer-Based Visual Servoing and Vibration Control of Flexible Robotic Manipulators With a Fixed Camera.

This article studies the visual servoing and vibration suppression control for flexible manipulators when the system states are unmeasurable and only the image feedback is available. The dynamic equations of flexible manipulators are decomposed into the slow and fast subsystems based on the singular perturbation theory. The nonlinear observers based on the state transformation using the Lie derivatives are proposed to estimate the unmeasurable system states and unknown camera intrinsic parameters at the same time. Then, the image-based controllers utilizing the estimated states are, respectively, designed in the slow and fast subsystems to regulate the image positions of feature points and suppress the vibration of flexible manipulators simultaneously. In the proposed approach, only the visual feedback is required to generate the control input for flexible manipulators, which simplifies the controller implementation. The stability of the proposed control scheme is proved based on the Lyapunov theory. Finally, experimental results on a flexible single-link manipulator are provided to demonstrate the effectiveness of the proposed control approach.

Edge Private Graph Neural Networks with Singular Value Perturbation

Graph neural networks (GNNs) play a key role in learning representations from graph-structured data and are demonstrated to be useful in many applications. However, the GNN training pipeline has been shown to be vulnerable to node feature leakage and edge extraction attacks. This paper investigates a scenario where an attacker aims to recover private edge information from a trained GNN model. Previous studies have employed differential privacy (DP) to add noise directly to the adjacency matrix or a compact graph representation. The added perturbations cause the graph structure to be substantially morphed, reducing the model utility. We propose a new privacy-preserving GNN training algorithm, Eclipse, that maintains good model utility while providing strong privacy protection on edges. Eclipse is based on two key observations. First, adjacency matrices in graph structures exhibit low-rank behavior. Thus, Eclipse trains GNNs with a low-rank format of the graph via singular values decomposition (SVD), rather than the original graph. Using the low-rank format, Eclipse preserves the primary graph topology and removes the remaining residual edges. Eclipse adds noise to the low-rank singular values instead of the entire graph, thereby preserving the graph privacy while still maintaining enough of the graph structure to maintain model utility. We theoretically show Eclipse provide formal DP guarantee on edges. Experiments on benchmark graph datasets show that Eclipse achieves significantly better privacy-utility tradeoff compared to existing privacy-preserving GNN training methods. In particular, under strong privacy constraints (𝜖 < 4), Eclipse shows significant gains in the model utility by up to 46%. We further demonstrate that Eclipse also has better resilience against common edge attacks (e.g., LPA), lowering the attack AUC by up to 5% compared to other state-of-the-art baselines.

Predicting transient dynamics in a model of reed musical instrument with slowly time-varying control parameter.

When playing a self-sustained reed instrument (such as the clarinet), initial acoustical transients (at the beginning of a note) are known to be of crucial importance. Nevertheless, they have been mostly overlooked in the literature on musical instruments. We investigate here the dynamic behavior of a simple model of reed instrument with a time-varying blowing pressure accounting for attack transients performed by the musician. In practice, this means studying a one-dimensional non-autonomous dynamical system obtained by slowly varying in time the bifurcation parameter (the blowing pressure) of the corresponding autonomous systems, i.e., whose bifurcation parameter is constant. In this context, the study focuses on the case for which the time-varying blowing pressure crosses the bistability domain (with the coexistence of a periodic solution and an equilibrium) of the corresponding autonomous model. Considering the time-varying blowing pressure as a new (slow) state variable, the considered non-autonomous one-dimensional system becomes an autonomous two-dimensional fast-slow system. In the bistability domain, the latter has attracting manifolds associated with two stable branches of the bifurcation diagram of the system with constant parameter. In the framework of the geometric singular perturbation theory, we show that a single solution of the two-dimensional fast-slow system can be used to describe the global system behavior. Indeed, this allows us to determine, depending on the initial conditions and rate of change of the blowing pressure, which manifold is approached when the bistability domain is crossed and to predict whether a sound is produced during transient as a function of the musician's control.

Exponential stabilization of linear systems with feedback optimization and its application to microgrids control

This paper proposes a feedback-optimization-based control method for linear time-invariant systems, which is aimed to exponentially stabilize the system and, meanwhile, drive the system output to an approximate solution of an optimization problem with convex set constraints and affine inequality constraints. To ensure the exponential stability of the closed-loop system, the original optimization problem is first reformulated into a counterpart that has only convex set constraints. It is shown that the optimal solution of the new optimization problem is an approximate optimal solution of the original problem. Then, based on this new optimization problem, the projected primal–dual gradient dynamics algorithm is used to design the controller. By using the singular perturbation method, sufficient conditions are provided to ensure the exponential stability of the closed-loop system. The proposed method is also applied to microgrid control.

Robust mixed finite element methods for a quad-curl singular perturbation problem

Robust mixed finite element methods are developed for a quad-curl singular perturbation problem. Lower order H(gradcurl)-nonconforming but H(curl)-conforming finite elements are constructed, which are extended to nonconforming finite element Stokes complexes and the associated commutative diagrams. Then H(gradcurl)-nonconforming finite elements are employed to discretize the quad-curl singular perturbation problem, which possess the sharp and uniform error estimates with respect to the perturbation parameter. The Nitsche’s technique is exploited to achieve the optimal convergence rate in the case of the boundary layers. Numerical results are provided to verify the theoretical convergence rates. In addition, the regularity of the quad-curl singular perturbation problem is established.

A dynamical systems approach to WKB-methods: The simple turning point

In this paper, we revisit the classical linear turning point problem for the second order differential equation ϵ2x″+μ(t)x=0 with μ(0)=0,μ′(0)≠0 for 0<ϵ≪1. Written as a first order system, t=0 therefore corresponds to a turning point connecting hyperbolic and elliptic regimes. Our main result is that we provide an alternative approach to WBK that is based upon dynamical systems theory, including GSPT and blowup, and we bridge – perhaps for the first time – hyperbolic and elliptic theories of slow-fast systems. As an advantage, we only require finite smoothness of μ. The approach we develop will be useful in other singular perturbation problems with hyperbolic–to–elliptic turning points.

Output-feedback path-following control of underactuated AUVs via singular perturbation and interconnected-system technique

This paper focuses on the output-feedback control for path-following of underactuated autonomous underwater vehicles subject to multiple uncertainties and unmeasured velocities. First, a novel extended state observer is proposed to estimate the mismatched lumped disturbance and recover the unmeasured velocities. Based on this premise, to overcome the limitation of relying solely on the accurate kinematic model, a disturbance observer-based stabilizing controller is developed. The difference in bandwidths between the observer and the vehicle dynamics allows for a mathematical setup amenable to standard singular perturbation theory. In the fast mode, a kinematic observer is designed to reject system uncertainty caused by unknown attack angular velocity and prohibitive path-tangential angular velocity, using a novel physical perspective. In the slow mode, an interconnected-system control law is proposed by integrating the backstepping technique with the time scale decomposition method. Furthermore, the stability of the overall closed-loop system is established. Finally, simulation results are presented to demonstrate the effectiveness of the proposed method for path-following of underactuated autonomous underwater vehicles in the vertical plane.

On inherent hyperelastic crease

We present analyses of the inherent hyperelastic crease that entails an unstable degenerative singular perturbation field over a uniformly compressed hyperelastic half-space. The analyses reveal asymptotic far-field characteristics of the singular field that bring out admissible incremental elastic deformation mechanisms of the inherent creasing instability typically observed at a compressive strain of ∼0.354 in neo-Hookean solids. The admissible singular perturbation field has an asymptotic leading-order incremental-stress singularity decaying 1/r2 as the distance from the crease tip r→∞. In general, asymptotic incremental-deformation fields of 1/r2 stress singularity can be introduced by surface flaws. We found that the singular field creates a far-field eigenmode of three energy-release angular sectors separated by two energy-elevating sectors of incremental deformation. The far-field eigenmode braces the near-tip energy-release field of the surface flaw against the transition to an unstable self-similar expansion field of creasing. Two asymptotic far-field parameters can characterize the braced-incremental-deformation (bid) and crease fields. One is the Burgers vector of a projected subsurface dislocation of the asymptotically singular far field. The other is a dimensionless shape factor representing the ratio of the subsurface depth of the dislocation to the Burgers vector. Our analyses show that the shape factor gauges the transition from the stable bid field to the unstable inherent crease field at the crease limit point as the compression increases. Two tangential-manifold conservation integrals are revealed to identify the two parameters. At the crease-limit point, matched asymptotes lead to the ratio between the two separate scaling parameters of the asymptotic solutions of the crease-tip folding field and the leading-order far field. To this end, we introduced a novel finite element method (FEM) for simulating the crease field nearly in a hyperelastic half-space with a finite-size simulation domain. Furthermore, we uncovered with the Gent model (1996) that strain-stiffening in hyperelasticity alters the dependence of the shape factor on the compressive strain, raising crease resistance. The new findings in hyperelastic crease mechanisms will help study ruga mechanics of self-organization and design soft-material structures and skin conditions for high crease resistance.

Canard Cycles and Their Cyclicity of a Fast–Slow Leslie–Gower Predator–Prey Model with Allee Effect

We study canard cycles and their cyclicity of a fast–slow Leslie–Gower predator–prey system with Allee effect. More specifically, we find necessary and sufficient conditions of the exact number (zero, one or two) of positive equilibria of the slow–fast system and its location (or their locations), and then we further completely determine its (or their) dynamics under explicit conditions. Besides, by geometric singular perturbation theory and the slow–fast normal form, we find explicit sufficient conditions to characterize singular Hopf bifurcation and canard explosion of the system. Additionally, the cyclicity of canard cycles is completely solved, and of particular interest is that we show the existence and uniqueness of a canard cycle, whose cyclicity is at most two, under corresponding precise explicit conditions.

Mode stability and shallow quasinormal modes of Kerr–de Sitter black holes away from extremality

A Kerr–de Sitter black hole is a solution (M,g_{\Lambda,\mathfrak{m},\mathfrak{a}}) of the Einstein vacuum equations with cosmological constant \Lambda&gt;0 . It describes a black hole with mass \mathfrak{m}&gt;0 and specific angular momentum \mathfrak{a}\in\mathbb{R} . We show that for any \varepsilon&gt;0 there exists \delta&gt;0 so that mode stability holds for the linear scalar wave equation \Box_{g_{\Lambda,\mathfrak{m},\mathfrak{a}}}\phi=0 when |\mathfrak{a}/\mathfrak{m}|\in[0,1-\varepsilon] and \Lambda\mathfrak{m}^{2}&lt;\delta . In fact, we show that all quasinormal modes \sigma in any fixed half-space \operatorname{Im}\sigma&gt;-C\sqrt{\Lambda} are equal to 0 or -i\sqrt{\Lambda/3}(n+o(1)) , n\in\mathbb{N} , as \Lambda\mathfrak{m}^{2}\searrow 0 . We give an analogous description of quasinormal modes for the Klein–Gordon equation. We regard a Kerr–de Sitter black hole with small \Lambda\mathfrak{m}^{2} as a singular perturbation either of a Kerr black hole with the same angular momentum-to-mass ratio, or of de Sitter spacetime without any black hole present. We use the mode stability of subextremal Kerr black holes, proved by Whiting and Shlapentokh-Rothman, as a black box; the quasinormal modes described by our main result are perturbations of those of de Sitter space. Our proof is based on careful uniform a priori estimates, in a variety of asymptotic regimes, for the spectral family and its de Sitter and Kerr model problems in the singular limit \Lambda\mathfrak{m}^{2}\searrow 0 .

The soliton and periodic solutions for the perturbed sine–cosine–Gordon equation

In the paper, the existence of soliton and periodic solutions is studied for the perturbed sine–cosine–Gordon equation. The corresponding traveling wave system is converted to a regularly Hamilton system via using bifurcation theory of differential equations and the geometric singular perturbation theory. Then by using Melnikov's method and symbolic computation, periodic wave solutions, solitary solutions, and kink and antikink solutions are obtained. The wave speed selection principles are given, respectively. Numerical simulations are performed to illustrate our theoretical analysis.

Numerical method for second order singularly perturbed delay differential equations with fractional order in time via fitted computational method

The numerical solution of time fractional parabolic differential equations with singular perturbations and delay is the subject of this article. An arbitrarily small perturbation parameter ɛ(0<ɛ<<1) is multiplied the highest order derivative term. The problem’s solution shows an exponential boundary layer on the right side of the spatial domain for small values of ɛ. The solution and its derivatives are discussed together with their properties and bounds. The Crank–Nicolson approach for time discretization and the fitted operator non-standard finite difference method for space discretization are used to solve the time fractional singularly perturbed time delay problem under consideration. The comparison principle and solution bound are used to examine and analyze the scheme’s stability. The scheme’s uniform convergence is examined and demonstrated. The proposed scheme has a uniform convergence order of Mx−1+(Δt)2−γ. Two numerical illustrations for various values of ɛ are taken into consideration to validate the theoretical analysis of the scheme. The constructed scheme provides a better order of convergence and more accurate results than methods found in the literature.

Nonstandard Nearly Exact Analysis of the FitzHugh–Nagumo Model

The FitzHugh–Nagumo model has been used empirically to model certain types of neuronal activities. It is also a non-linear dynamical system applicable to chemical kinetics, population dynamics, epidemiology and pattern formation. In the literature, many approaches have been proposed to study its dynamics. In this paper, initially, we have employed cutting-edge tools from discrete dynamics for discretization and fixed points. It has been proven that an exact discrete scheme exists for this paradigm. This project also considers the phase space and integral surfaces of these evolutionary equations. In addition, it carries out a thorough symmetry analysis of this reaction diffusion system to find equivalent systems. Moreover, steady-state solutions are obtained using ansatzes for traveling wave solutions. The existence of infinite traveling wave solutions has also been proven. Yet again, this investigation establishes the potential of symmetry methods to unravel non-linearity. Finally, singular perturbation theory has been employed to obtain analytical approximations and to study stability in different parameter regimes.

Identification and Control of Flexible Joint Robots Based on a Composite-Learning Optimal Bounded Ellipsoid Algorithm and Prescribe Performance Control Technique

This paper presents an indirect adaptive neural network (NN) control algorithm tailored for flexible joint robots (FJRs), aimed at achieving desired transient and steady-state performance. To simplify the controller design process, the original higher-order system is decomposed into two lower-order subsystems using the singular perturbation technique (SPT). NNs are then employed to reconstruct the aggregated uncertainties. An adaptive prescribed performance control (PPC) strategy and a continuous terminal sliding mode control strategy are introduced for the reduced slow subsystem and fast subsystem, respectively, to guarantee a specified convergence speed and steady-state accuracy for the closed-loop system. Additionally, a composite-learning optimal bounded ellipsoid algorithm (OBE)-based identification scheme is proposed to update the NN weights, where the tracking errors of the reduced slow and fast subsystems are integrated into the learning algorithm to enhance the identification and tracking performance. The stability of the closed-loop system is rigorously established using the Lyapunov approach. Simulations demonstrate the effectiveness of the proposed identification and control schemes.

Uniqueness of Single Peak Solutions for a Kirchhoff Equation

We deal with the following singular perturbation Kirchhoff equation: −ϵ2a+ϵb∫R3|∇u|2dyΔu+Q(y)u=|u|p−1u,u∈H1(R3), where constants a,b,ϵ&gt;0 and 1&lt;p&lt;5. In this paper, we prove the uniqueness of the concentrated solutions under some suitable assumptions on asymptotic behaviors of Q(y) and its first derivatives by using a type of Pohozaev identity for a small enough ϵ. To some extent, our result exhibits a new phenomenon for a kind of Q(x) which allows for different orders in different directions.

New solitary waves in a convecting fluid

It is critical to examine the effects of small perturbations on solitary wave solutions to nonlinear wave equations because of the existence of the unavoidable perturbations in real applications. In this work, we aim to study the solitary wave solutions for the dissipative perturbed mKdV equation which is proposed to model the evolution of shallow wave in a convecting fluid. By using the geometric singular perturbation theorem, the three-dimensional traveling wave system is reduced to a planar dynamical system with regular perturbation. By examining the homoclinic bifurcations via Melnikov method, we show that not only some solitary waves with particularly chosen wave speed persist but also some new type of solitary waves appear under small perturbation. The theoretical analysis results are illustrated by numerical simulations.