Singular Perturbation Problems Research Articles

General-Kindred physics-informed neural network to the solutions of singularly perturbed differential equations

Physics-informed neural networks (PINNs) have become a promising research direction in the field of solving partial differential equations (PDEs). Dealing with singular perturbation problems continues to be a difficult challenge in the field of PINN. The solution of singular perturbation problems often exhibits sharp boundary layers and steep gradients, and traditional PINN cannot achieve approximation of boundary layers. In this manuscript, we propose the General-Kindred physics-informed neural network (GKPINN) for solving singular perturbation differential equations (SPDEs). This approach utilizes asymptotic analysis to acquire prior knowledge of the boundary layer from the equation and establishes a novel network to assist PINN in approximating the boundary layer. It is compared with traditional PINN by solving examples of one-dimensional, two-dimensional, and time-varying SPDE equations. The research findings underscore the exceptional performance of our novel approach, GKPINN, which delivers a remarkable enhancement in reducing the L2 error by two to four orders of magnitude compared to the established PINN methodology. This significant improvement is accompanied by a substantial acceleration in convergence rates, without compromising the high precision that is critical for our applications. Furthermore, GKPINN still performs well in extreme cases with perturbation parameters of 1×10−38, demonstrating its excellent generalization ability.

An Efficient Computational Algorithm for a Class of Nonlinear Singular Perturbation Problems of Convection Diffusion Type with Cauchy Data

An Efficient Computational Algorithm for a Class of Nonlinear Singular Perturbation Problems of Convection Diffusion Type with Cauchy Data

Rectifiability, finite Hausdorff measure, and compactness for non‐minimizing Bernoulli free boundaries

AbstractWhile there are numerous results on minimizers or stable solutions of the Bernoulli problem proving regularity of the free boundary and analyzing singularities, much less is known about critical points of the corresponding energy. Saddle points of the energy (or of closely related energies) and solutions of the corresponding time‐dependent problem occur naturally in applied problems such as water waves and combustion theory. For such critical points —which can be obtained as limits of classical solutions or limits of a singular perturbation problem—it has been open since (Weiss, 2003) whether the singular set can be large and what equation the measure satisfies, except for the case of two dimensions. In the present result we use recent techniques such as a frequency formula for the Bernoulli problem as well as the celebrated Naber–Valtorta procedure to answer this more than 20 year old question in an affirmative way: For a closed class we call variational solutions of the Bernoulli problem, we show that the topological free boundary (including degenerate singular points , at which as ) is countably ‐rectifiable and has locally finite ‐measure, and we identify the measure completely. This gives a more precise characterization of the free boundary of in arbitrary dimension than was previously available even in dimension two. We also show that limits of (not necessarily minimizing) classical solutions as well as limits of critical points of a singularly perturbed energy are variational solutions, so that the result above applies directly to all of them.

Rate of Convergence for First-Order Singular Perturbation Problems: Hamilton–Jacobi–Isaacs Equations and Mean Field Games of Acceleration

Rate of Convergence for First-Order Singular Perturbation Problems: Hamilton–Jacobi–Isaacs Equations and Mean Field Games of Acceleration

Theoretical and Numerical Analysis of Singular Perturbation Problems in Ordinary Differential Equations

Background: Singular perturbation troubles in everyday differential equations (ODEs) involve a small parameter ϵ that causes answers to show off behaviors on multiple scales, main to massive challenges in both theoretical and numerical evaluation. Objective: This paper at targets to comprehensively look at the theoretical and numerical techniques for reading and solving singular perturbation troubles in ODEs. Focus is located on know-how the behaviors precipitated through the small parameter and developing strong numerical techniques to accurately seize these behaviors. Methods: Theoretical approaches employed include matched asymptotic expansions and multiple scale analysis. These methods decompose the solution domain into regions with different scales, constructing and matching approximate solutions to ensure smooth transitions. Numerical techniques such as finite difference methods, finite element methods, and spectral methods are utilized, with particular emphasis on adaptive mesh refinement to handle boundary layers effectively Results: Theoretical evaluation demonstrates the effectiveness of matched asymptotic expansions and multiple scale analysis in presenting correct approximations and easy transitions among one of a kind solution regions. Numerical techniques showed high accuracy and performance, particularly whilst blended with adaptive techniques. Finite distinction strategies with non-uniform grids and adaptive mesh refinement, finite detail techniques with adaptive mesh strategies, and spectral strategies with special handling of boundary layers all proved successful. Stability and convergence analyses showed the reliability of those techniques. Conclusions: This comprehensive evaluation highlights the strengths and applicability of each theoretical and numerical strategies in tackling singular perturbation issues in ODEs. The mixed use of these techniques lets in for correct and green answers, presenting treasured equipment for applications in diverse clinical and engineering fields.

A parameter-uniform hybrid method for singularly perturbed parabolic 2D convection-diffusion-reaction problems

The solution of the singular perturbation problems (SPP) of convection-diffusion-reaction type may exhibit regular and corner layers in a rectangular domain. In this work, we construct and analyze a parameter-uniform operator-splitting alternating direction implicit (ADI) scheme to efficiently solve a two-dimensional parabolic singularly perturbed problem with two positive parameters. The proposed model is a combination of the backward-Euler method defined on a uniform mesh in time and a hybrid method in space defined on a special Shishkin mesh. The analysis is presented on a layer adapted piecewise-uniform Shishkin mesh. The developed numerical method is proved to be first-order convergent in time and almost second-order convergent in space. The numerical experiments are performed to validate the theoretical convergence results and illustrate the efficiency of the current strategy.

A singular perturbation result for a class of periodic-parabolic BVPs

Abstract In this article, we obtain a very sharp version of some singular perturbation results going back to Dancer and Hess [Behaviour of a semilinear periodic-parabolic problem when a parameter is small, Lecture Notes in Mathematics, Vol. 1450, Springer-Verlag, Berlin, 1990, pp. 12–19] and Daners and López-Gómez [The singular perturbation problem for the periodic-parabolic logistic equation with indefinite weight functions, J. Dynam. Differential Equations 6 (1994), 659–670] valid for a general class of semilinear periodic-parabolic problems of logistic type under general boundary conditions of mixed type. The results of Dancer and Hess [Behaviour of a semilinear periodic-parabolic problem when a parameter is small, Lecture Notes in Mathematics, Vol. 1450, Springer-Verlag, Berlin, 1990, pp. 12–19] and [The singular perturbation problem for the periodic-parabolic logistic equation with indefinite weight functions, J. Dynam. Differential Equations 6 (1994), 659–670] were found, respectively, for Neumann and Dirichlet boundary conditions with L = − Δ {\mathfrak{L}}=-\Delta . In this article, L {\mathfrak{L}} stands for a general second-order elliptic operator.

A Posteriori Error Analysis of Defect Correction Method for Singular Perturbation Problems With Discontinuous Coefficient and Point Source

Abstract The paper presents a defect correction method to solve singularly perturbed problems with discontinuous coefficient and point source. The method combines an inexpensive, lower-order stable, upwind difference scheme and a higher-order, less stable central difference scheme over a layer-adapted mesh. The mesh is designed so that most mesh points remain in the regions with rapid transitions. A posteriori error analysis is presented. The proposed numerical method is analyzed for consistency, stability, and convergence. The error estimates of the proposed numerical method satisfy parameter-uniform second-order convergence on the layer-adapted grid. The convergence obtained is optimal because it is free from any logarithmic term. The numerical analysis confirms the theoretical error analysis.

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.

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.

Lowest-degree robust finite element schemes for inhomogeneous bi-Laplace problems

In this paper, we study the numerical method for the bi-Laplace problems with inhomogeneous coefficients; particularly, we propose finite element schemes on rectangular grids respectively for an inhomogeneous fourth-order elliptic singular perturbation problem and for the Helmholtz transmission eigenvalue problem. The new methods use the reduced rectangle Morley (RRM for short) element space with piecewise quadratic polynomials, which are of the lowest degree possible. For the finite element space, a discrete analogue of an equality by Grisvard is proved for the stability issue and a locally-averaged interpolation operator is constructed for the approximation issue. Optimal convergence rates of the schemes are proved, and numerical experiments are given to verify the theoretical analysis.

A Singular Tempered Sub-Diffusion Fractional Equation with Changing-Sign Perturbation

In this paper, we establish some new results on the existence of positive solutions for a singular tempered sub-diffusion fractional equation involving a changing-sign perturbation and a lower-order sub-diffusion term of the unknown function. By employing multiple transformations, we transform the changing-sign singular perturbation problem to a positive problem, then establish some sufficient conditions for the existence of positive solutions of the problem. The asymptotic properties of solutions are also derived. In deriving the results, we only require that the singular perturbation term satisfies the Carathéodory condition, which means that the disturbance influence is significant and may even achieve negative infinity near some time singular points.

The virtual element method with interior penalty for the fourth-order singular perturbation problem

We present the virtual element method with interior penalty to solve a fourth-order singular perturbation problem. In order to estimate the nonconformity error, the degrees of freedom on edges are changed to the moments of functions in the interior penalty scheme. To do this, we design a special H1-type projection that can be uniquely determined by the new degrees of freedom. With the help of the H1-type projection, the local discrete space is defined by imposing some restrictions on the local space of H2-conforming VE, such that it has the local H2 regularity. Then for the numerical method, we derive the error estimates under the condition of enough smooth solution and the uniform error estimates for cases with boundary layers. Finally, we display some numerical examples. From that, we see that the interior penalty method gives a better performance on the convergence than the theoretical predictions.

Dynamics of a Class of Prey–Predator Models with Singular Perturbation and Distributed Delay

In this paper, two prey–predator models with distributed delays are presented based on the growth and loss rates of the predator, which are much smaller than that of the prey, leading to a singular perturbation problem. It is obtained that Hopf bifurcation can occur, where the coexistence equilibrium becomes unstable leading to a stable limit cycle. Subsequently, considering the perturbation parameter [Formula: see text], the fact that the solution crossing the transcritical point converges to a stable equilibrium is discussed for the model with Holling type I using the linear chain criterion, center-manifold reduction, the geometric singular perturbation theory and entry–exit function. The existence and uniqueness of relaxation oscillation cycle for the model with Holling type II are obtained. In addition, numerical simulations are provided to verify the analytical results.

Numerical analysis of singularly perturbed parabolic reaction diffusion differential difference equations

The authors present numerical analysis of singularly perturbed parabolic problems. The previous papers in this direction focussed on the convection-diffusion problems with regular type of boundary layers. It is very difficult and almost impossible as stated in the literature, e.g, in Chapter 14 of ‘Fitted numerical methods for singular perturbation problems: error estimates in the maximum norm for linear problems in one and two dimensions’ (1996) by Miller, John JH and O'Riordan, Eugene and Shishkin, Grigorii I, to capture singularities due to parabolic layers and develop a robust scheme for reaction-diffusion problems with general values of space shifts. The work is in progress now and this is the first paper in this direction. In this work, we are investigating such problems and develop a robust numerical scheme using Mickens's non-standard finite difference scheme and special mesh. Furthermore, numerical examples have been presented to verify the theoretical findings.

Asymptotic study of critical wave fronts for parameter-dependent Born–Infeld models: physically predicted behaviors and new phenomena

In this paper, we study some parameter-dependent reaction-diffusion models governed by the Born–Infeld (or Minkowski) operator. In dependence on two parameters , related to the field strength and to the diffusivity, we investigate the limit critical speed for traveling fronts, together with the limit behavior of the associated critical profiles. As the two main results, on the one hand we rigorously show that for arbitrarily large electric fields the critical speed and the critical profile converge to the ones of the linear diffusion problem, agreeing with the well-known physical prediction that, in this case, Born’s law and Maxwell’s law coincide. Such a result is accompanied by a counterpart for arbitrarily small electric fields and a complete analysis for vanishing/large diffusion. On the other hand, we prove the onset of a new and unexpected phenomenon for the singular perturbation problem: the critical speed of propagation does not converge to zero and, for KPP or combustion-type reactions, the limit front profile becomes sharp on one side only. In fact, this latter coincides with the C 1-gluing of a piecewise linear profile having slope equal to 1 when nonconstant, and of a regular inviscid profile propagating at the limit speed. The gluing point and the value of the limit speed are determined explicitly. The outcomes of the whole discussion, based on a careful analysis of the first-order reduction associated with the original equation, show substantial differences and peculiarities of the Born–Infeld operator with respect to linear or saturating diffusions, and are supported by several numerical experiments.

Multi-Scale-Matching neural networks for thin plate bending problem

Physics-informed neural networks are a useful machine learning method for solving differential equations, but encounter challenges in effectively learning thin boundary layers within singular perturbation problems. To resolve this issue, multi-scale-matching neural networks are proposed to solve the singular perturbation problems. Inspired by matched asymptotic expansions, the solution is decomposed into inner solutions for small scales and outer solutions for large scales, corresponding to boundary layers and outer regions, respectively. Moreover, to conform neural networks, we introduce exponential stretched variables in the boundary layers to avoid semi-infinite region problems. Numerical results for the thin plate problem validate the proposed method.

Alternative Methods of Regular and Singular Perturbation Problems

Alternative Methods of Regular and Singular Perturbation Problems

Piecewise approximate analytical solutions of high-order reaction-diffusion singular perturbation problems with boundary and interior layers

&lt;abstract&gt; &lt;p&gt;This work aims to present a reliable algorithm that can effectively generate accurate piecewise approximate analytical solutions for third- and fourth-order reaction-diffusion singular perturbation problems. These problems involve a discontinuous source term and exhibit both interior and boundary layers. The original problem was transformed into a system of coupled differential equations that are weakly interconnected. A zero-order asymptotic approximate solution was then provided, with known asymptotic analytical solutions for the boundary and interior layers, while the outer region solution was obtained analytically using an enhanced residual power series approach. This approach combined the standard residual power series method with the Padé approximation to yield a piecewise approximate analytical solution. It satisfies the continuity and smoothness conditions and offers higher accuracy than the standard residual power series method and other numerical methods like finite difference, finite element, hybrid difference scheme, and Schwarz method. The algorithm also provides error estimates, and numerical examples are included to demonstrate the high accuracy, low computational cost, and effectiveness of the method within a new asymptotic semi-analytical numerical framewor.&lt;/p&gt; &lt;/abstract&gt;