Singular Perturbation Problems Research Articles

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 analysed 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.

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.

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

<abstract> <p>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.</p> </abstract>

ASYMPTOTIC EXPANSION OF THE SOLUTION FOR SINGULARPERTURBED LINEAR IMPULSIVE SYSTEMS

In this study, a singularly perturbed linear impulsive system with singularly perturbed impulses is considered. Many books discuss different types of singular perturbation problems. In the present work, an impulse system is considered in which a small parameter is introduced into the impulse equation. This is the main novelty of our study, since other works [25] have only considered a small parameter in the differential equation. A necessary condition is also established to prevent the impulse function from bloating as the parameter approaches zero. As a result, the notion of singularity for discontinuous dynamics is greatly extended. An asymptotic expansion of the solution of a singularly perturbed initial problem with an arbitrary degree of accuracy for a small parameter is constructed. A theorem for estimating the residual term of the asymptotic expansion is formulated, which estimates the difference between the exact solution and its approximation. The results extend those of [32], which formulates an analogue of Tikhonov’s limit transition theorem. The theoretical results are confirmed by a modelling example.

A singular perturbation problem for a nonlinear Schrödinger system with three wave interaction

A singular perturbation problem for a nonlinear Schrödinger system with three wave interaction

A Numerical scheme to Solve Boundary Value Problems Involving Singular Perturbation

The Wang-Ball polynomials operational matrices of the derivatives are used in this study to solve singular perturbed second-order differential equations (SPSODEs) with boundary conditions. Using the matrix of Wang-Ball polynomials, the main singular perturbation problem is converted into linear algebraic equation systems. The coefficients of the required approximate solution are obtained from the solution of this system. The residual correction approach was also used to improve an error, and the results were compared to other reported numerical methods. Several examples are used to illustrate both the reliability and usefulness of the Wang-Ball operational matrices. The Wang Ball approach has the ability to improve the outcomes by minimizing the degree of error between approximate and exact solutions. The Wang-Ball series has shown its usefulness in solving any real-life scenario model as first- or second-order differential equations (DEs).

An asymptotic-preserving finite element method for a forth order singular perturbation problem with boundary layers

We propose an asymptotic-preserving finite element method for a fourth order singular perturbation problem, which completely preserves the asymptotic transition of the underlying partial differential equation. As a representative, we analyze the NZT element and a linear convergence rate is proved for the solution with sharp boundary layer. Numerical examples in two and three dimensions are consistent with the theoretical prediction.

Couette-Poiseuille flow in a fluid overlying an anisotropic porous layer

In this study, Couette-Poiseuille flow inside a channel partially filled with a porous medium is treated. Considering the applicability to the practical scenarios, we assume that the porous medium is anisotropic in nature. More precisely, the model under steady-state assumption is described by the coupling of Brinkman-Forchheimer and Stokes equations. We develop the existence and uniqueness results corresponding to this two-dimensional Stokes-Brinkman-Forchheimer coupled system using the Browder-Minty theorem. For this we restrict the permeability to a diagonal matrix. The coupled system involves a non-linear term due to the Brinkman-Forchheimer equation. We propose an algorithm based on shooting and searching techniques to develop a numerical solution for this system, however, for an arbitrary anisotropic permeability matrix. The robustness of the presented numerical solution is established by comparing it with asymptotic solutions corresponding to small and large Darcy numbers. For the case of a large Darcy number, the problem of interest is a regular perturbation while the same is a singular perturbation problem for the case of a small Darcy number, which is dealt with using Prandtl's matching principle. With the velocity distribution determined, we have presented a comparative analysis of the shear stress distribution at the liquid-porous interface and at the bottom plate, which is of interest to analyze the stress distribution on the arterial wall. The asymptotic results are found to be in good agreement with those obtained numerically.

Potentiality of the HOFiD_bvp code in solving different kind of second order boundary value problems

The aim of this paper is to present the potentiality of the code HOFiD_bvp based on high order finite difference, following the approach of the boundary value methods and the upwind methods. Strategies, as variable stepsize and variable order, that helps to handle with solutions of singular perturbation problems, so as the continuation strategy applied for solving nonlinear problems, will be discussed. Now, the method is suitable to solve second order scalar boundary value problems, not only singularly perturbed with one and two parameters, but also problems with discontinuous source terms. Concerning that point some numerical tests will be taken in consideration in order to show the potentiality of the code.

A robust higher-order finite difference technique for a time-fractional singularly perturbed problem

A higher-order finite difference method is developed to solve the variable coefficients convection–diffusion singularly perturbed problems (SPPs) involving fractional-order time derivative with the order α∈(0,1). The solution to this problem class possesses a typical weak singularity at the initial time t=0 and an exponential boundary layer at the right lateral surface as the perturbation parameter ɛ→0. Alikhanov’s L2−1σ approximation is applied in the temporal direction on a suitable graded mesh, and the spatial variable is discretized on a piecewise uniform Shishkin mesh using a combination of midpoint upwind and central finite difference operators. Stability estimates and the convergence analysis of the fully discrete scheme are provided. It is shown that the fully discrete scheme is uniformly convergent with a rate of O(M−p+N−2(logN)2), where p=min{2,rα}, r is the graded mesh parameter and M,N are the number of mesh points in the time and space direction, respectively. Two numerical examples are taken in counter to confirm the sharpness of the theoretical estimates.

A Parameter-Robust Numerical Method Based on Defect-Correction Technique for Parabolic Singular Perturbation Problems with Discontinuous Convection Coefficient and Source

The paper presents a uniformly convergent finite difference method based on the defect correction technique to solve parabolic singular perturbation problems with discontinuous convection coefficient and source. The solution to the problem exhibits interior layer across the discontinuity and demonstrates turning point behavior. The simultaneous presence of perturbation parameters and discontinuity makes the problem stiff. A higher-order method is developed using an implicit difference scheme in time on a uniform mesh and a combination of the upwind difference method and the central difference scheme over an adaptive mesh in space. The method involves iteratively solving increasingly accurate discrete problems by computing and using the defect to correct the approximate solution. Parameter uniform error estimates show uniform convergence of first-order in time and second-order in space. Numerical experiments confirm the accuracy of the proposed scheme and support the theoretical analysis.

A singular perturbation problem for mean field games of acceleration: application to mean field games of control

The singular perturbation of mean field game systems arising from minimization problems with control of acceleration is addressed, that is, we analyze the behavior of solutions as the acceleration costs vanishes. In this setting, the Hamiltonian fails to be strictly convex and coercive w.r.t. the momentum variable and, so, the classical results for Tonelli Hamiltonian systems cannot be applied. However, we show that the limit system is of MFG type in two different cases: we first study the convergence to the classical MFG system and, then, by a finer analysis of the Euler–Lagrange flow associated with the control of acceleration, we prove the convergence to a class of MFG systems, known as, MFG of control.