Singular Boundary Research Articles

An efficient computational technique and its convergence analysis for a class of doubly singular boundary value problems

An efficient computational technique and its convergence analysis for a class of doubly singular boundary value problems

A new fourth-order compact finite difference method for solving Lane-Emden-Fowler type singular boundary value problems

We develop a novel fourth-order compact finite difference scheme to solve nonlinear singular ordinary differential equations. Such problems occur in many fields of science and engineering, such as studying the equilibrium of an isothermal gas sphere, reaction–diffusion in a spherical permeable catalyst, etc. These problems are challenging to solve because of their singularity or nonlinearity. By our proposed method, we can easily solve these complex problems without removing or modifying the singularity. To construct the new fourth-order compact difference method, Initially, we created a uniform mesh within the solution domain and developed a compact finite difference scheme. This scheme approximates the derivatives at the boundary nodal points to handle the problem’s singularity effectively. Employing a matrix analysis approach, we discussed the convergence analysis of the methods. To demonstrate its efficacy, we apply our approach to solve various real-life problems from the literature. The new method offers high-order accuracy with minimal grid points and provides better numerical results than the nonstandard finite difference method and exponential compact finite difference method.

An efficient fourth-order convergent scheme based on half-step spline function for two-point mixed boundary value problems

The present work aims to introduce a novel numerical method for solving second-order two-point mixed boundary value problems. Mixed boundary value problems occur in various scientific fields, including quantum mechanics, fluid dynamics, and chemical reactor theory. The proposed method provides a highly accurate, fourth-order convergent numerical solution, achieved by implementing a compact finite difference method based on non-polynomial spline in tension approximations. Notably, the discretisation involves the use of half-step grid points, eliminating the need to modify the method at singularities while dealing with singular boundary value problems. The article also includes a comprehensive convergence analysis that demonstrates the theoretical order of convergence of the method. Additionally, the results obtained from solving six diverse mixed boundary value problems, including Burger's equation and Bratu-type equation, are compared to showcase the superiority and effectiveness of the proposed method.

Classical echoes of quantum boundary conditions

Abstract We consider a non-relativistic particle in a one-dimensional box with all possible quantum boundary conditions that make the kinetic-energy operator self-adjoint. We determine the Wigner functions of the corresponding eigenfunctions and analyze in detail their classical limit, governed by their behavior in the high-energy regime. We show that the quantum boundary conditions split into two classes: all local and regular boundary conditions collapse to the same classical boundary condition, while a dependence on singular non-local boundary conditions persists in the classical limit.

ANN based optimization of nano-beam oscillations with intermolecular forces and geometric nonlinearity

In this study, we investigate the effect of Van der Waals and Casimir forces on the mathematical model of nano-electromechanical systems (NEMS) such as nano-beam actuators that contain cantilever and double cantilever beams. The singular nonlinear boundary value problem governing the beam-type actuators, including geometric nonlinearity is solved by using an intelligent strength of feedforward artificial neural networks (ANNs) and hybridization of optimization algorithms such as arithmetic optimization algorithm (AOA) and active set algorithm (ASA). The proposed ANN-AOA-AS algorithm is employed to quantify the effect of changes in applied voltage, dispersion forces, geometric nonlinearity parameters, and initial axial strain on the deflection of the beam. Furthermore, to validate the results obtained by the proposed algorithm, statistical analyses are conducted to compare the approximate solutions with state-of-the-art methodologies available in the latest literature. In addition, performance indicators are defined such as mean square error (MSE), Nash–Sutcliffe efficiency (NSE), mean absolute deviations (MAD), root mean square error (RMSE), and Error in Nash–Sutcliffe efficiency (ENSE) to study the accuracy and efficiency of the solutions. The results show that these indicators’ mean percentage values lie around 10−4 to 10−6 which reflects the perfect modeling of the approximate solutions.

Regularity for one-phase Bernoulli problems with discontinuous weights and applications

We study a one-phase Bernoulli free boundary problem with weight function admitting a discontinuity along a smooth jump interface. In any dimension N ≥ 2 N\ge 2 , we show the C 1 , α C^{1, \alpha } regularity of the free boundary outside of a singular set of Hausdorff dimension at most N − 3 N-3 . In particular, we prove that the free boundaries are C 1 , α C^{1, \alpha } regular in dimension N = 2 N=2 , while in dimension N = 3 N=3 the singular set can contain at most a finite number of points. We use this result to construct singular free boundaries in dimension N = 2 N=2 , which are minimizing for one-phase functionals with weight functions in L ∞ L^\infty that are arbitrarily close to a positive constant.

Construction of mathematical models of thermal conductivity for modern electronic devices with elements of a layered structure

This paper considers a heat conduction process for isotropic layered media with internal thermal heating. As a result of the heterogeneity of environments, significant temperature gradients arise as a result of the thermal load. In order to establish the temperature regimes for the effective operation of electronic devices, linear and non-linear mathematical models for determining the temperature field have been constructed, which would make it possible to further analyze the temperature regimes in these heat-active environments. The coefficient of thermal conductivity for the above structures was represented as a whole using asymmetric unit functions. As a result, the conditions of ideal thermal contact were automatically fulfilled on the surfaces of the conjugation of the layers. This leads to solving one heat conduction equation with discontinuous and singular coefficients and boundary conditions at the boundary surfaces of the medium. For linearization of nonlinear boundary value problems, linearizing functions were introduced. Analytical solutions to both linear and nonlinear boundary value problems were derived in a closed form. For heat-sensitive environments, as an example, the linear dependence of the coefficient of thermal conductivity of structural materials on temperature, which is often observed when solving many practical problems, was chosen. As a result, analytical relations for determining the temperature distribution in these environments were obtained. Based on this, a numerical experiment was performed, and it was geometrically represented depending on the spatial coordinates. This proves that the constructed linear and nonlinear mathematical models testify to their adequacy to the real physical process. They make it possible to analyze heat-active media regarding their thermal resistance. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual nodes and their elements but also the entire structure

A highly accurate and efficient Genocchi‐based spectral technique applied to singular fractional order boundary value problems

This article focuses on an efficient and highly accurate approximate solver for a class of generalized singular boundary value problems (SBVPs) having nonlinearity and with two‐term fractional derivatives. The involved fractional derivative operators are given in the form of Liouville–Caputo. The developed algorithm for solving the generalized SBVPs consists of two main stages. The first stage is devoted to an iterative quasilinearization method (QLM) to conquer the (strong) nonlinearity of the governing SBVPs. Secondly, we employ the generalized Genocchi polynomials (GGPs) to treat the resulting sequence of linearized SBVPs numerically. An upper error estimate for the Genocchi series solution in the norm is obtained via a rigorous error analysis. The main benefit of the presented QLM‐GGPs procedure is that the required number of iteration in the first stage is within a few steps, and an accurate polynomial solution is obtained through computer implementations in the second stage. Three widely applicable test cases are investigated to observe the efficacy as well as the high‐order accuracy of the QLM‐GGPs algorithm. The comparable accuracy and robustness of the presented algorithm are validated by doing comparisons with the results of some well‐established available computational methods. It is apparently shown that the QLM‐GGPs algorithm provides a promising tool to solve strongly nonlinear SBVPs with two‐term fractional derivatives accurately and efficiently.

Nonlinear Contractions Employing Digraphs and Comparison Functions with an Application to Singular Fractional Differential Equations

After the initiation of Jachymski’s contraction principle via digraph, the area of metric fixed point theory has attracted much attention. A number of outcomes on fixed points in the context of graph metric space employing various types of contractions have been investigated. The aim of this paper is to investigate some fixed point theorems for a class of nonlinear contractions in a metric space endued with a transitive digraph. The outcomes presented herewith improve, extend and enrich several existing results. Employing our findings, we describe the existence and uniqueness of a singular fractional boundary value problem.

A fast semi-analytical meshless method in two-dimensions

This paper introduces a fast direct algorithm for the singular boundary method (SBM) in two-dimensional (2D) problems, utilizing the hierarchical off-diagonal low-rank (HODLR) matrix concept as the foundation of the fast direct solver. The HODLR matrix is constructed by hierarchically partitioning the coefficient matrix into blocks using a binary tree, with all off-diagonal blocks at each level being represented as low-rank factors. Furthermore, The Sherman-Morrison-Woodbury formula can be used to efficiently compute the inverse of the HODLR matrix. The numerical experiment results demonstrate that the new fast solver can significantly improve the computational efficiency of the SBM while maintaining the same level of precision.

Spectral Solutions of Specific Singular Differential Equations Using A Unified Spectral Galerkin-Collocation Algorithm

This paper presents a novel numerical approach to addressing three types of high-order singular boundary value problems. We introduce and consider three modified Chebyshev polynomials (CPs) of the third kind as proposed basis functions for these problems. We develop new derivative operational matrices for the three modified CPs of the third kind by deriving formulas for their first derivatives. Our approach follows a unified method for numerically handling singular differential equations (DEs). To transform these equations into algebraic systems suitable for numerical treatment, we employ the collocation method in combination with the introduced operational matrices of derivatives of the modified CPs of the third kind. We address the convergence examination for the three expansions in a unified manner. We present numerous numerical examples to demonstrate the accuracy and efficiency of our unified numerical approach.

On the regularity of solutions to a class of nonlinear Volterra integral equations with singularities

We study the smoothness properties of solutions to nonlinear Volterra integral equations of the second kind on a bounded interval [0,b]. The kernel of the integral operator of the underlying equation may have a diagonal singularity and a boundary singularity. Information about them is given through certain estimates. To characterize the regularity of solutions of such equations we show that the solution belongs to an appropriately weighted space of smooth functions on (0,b], with possible singularities of the derivatives of the solution at the left endpoint of the interval [0,b].

Analytical method for systems of nonlinear singular boundary value problems

The standard Lane–Emden equations model several physical phenomena such as isotropic continuous media, thermal behaviour of a spherical cloud of gas, and isothermal gas spheres. Systems of Lane–Emden-type equations appear in the modelling of the concentration of carbon substrate and oxygen, catalytic diffusion reactions, the steady state concentration of carbon dioxide, dusty fluid models, and pattern formation. In solving singular boundary value problems, one faces challenges resulting from the divergence of the associated variable coefficients at the singular points. This research article explores the power series approach to analytically approximate solutions to a class of systems of strongly nonlinear singular boundary value problems of Lane–Emden-type. The nonlinear terms in the proposed problems are transformed into power series using the generalised Cauchy product before establishing explicit recursion formulae for the expansion coefficients of the system of series solutions. The initial conditions required in the proposed boundary value problems are assumed and determined from a set of nonlinear algebraic equations resulting from the given right boundary conditions. Three special systems of nonlinear singular boundary value problems of Lane–Emden type are presented to demonstrate the proposed method’s reliability, effectiveness, and accuracy. The obtained approximate solutions are compared with the exact solutions (where they are available), or with other existing results (where the exact solution is not readily available). The series solution of the first example has a slow convergent rate, while the results of the other two examples are in excellent agreement with the exact solutions and other published results.

Bifurcation analysis in the regularized four-sided cavity flow: Equilibrium states in a D2-symmetric fluid system

In this study, we investigate a handful of equilibrium states and their hydrodynamic stability for the incompressible flow in a square cavity driven by the tangential simultaneous motion of all of its lids at the same speed. This problem has been investigated in the recent past, albeit only partially and always disregarding the singular boundary conditions at the corners. The implications of the system symmetries have not been rigorously addressed, either. In our study, we introduce a regularized version of the boundary conditions to avoid the corner singularities, in a way that preserves the main features of the original setup. In addition, the Navier–Stokes equations are discretized by means of highly accurate Chebyshev spectral methods that provide exponential convergence of all computed flows and consistently eliminate any potential source of structural instability of the bifurcation scenario. We employ Newton–Krylov solvers, implemented within continuation algorithms, to accurately compute equilibrium solutions. Linear stability analysis of both the primary symmetric base flow and the secondary asymmetric states uncovers new branches of fully asymmetric steady states. The analysis has allowed identification of six previously undetected bifurcations, all of which are associated with the disruption of either the rotational invariance or the reflection symmetry. Some of these bifurcations have been found to be quite clustered in some regions of the parameter space, which points at the underlying action of higher codimension mechanisms. Notably, all the bifurcations reported occur within the range of low to moderate Reynolds numbers, making the regularized four-sided lid-driven cavity flow a reliable benchmark for assessing different numerical schemes in the context of Navier–Stokes equivariant bifurcation theory.

On a localization-in-frequency approach for a class of elliptic problems with singular boundary data

We consider a class of nonhomogeneous elliptic equations in the half-space with critical singular boundary potentials and nonlinear fractional derivative terms. The forcing terms are considered on the boundary and can be taken as singular measure. Employing a functional setting and approach based on localization-in-frequency and Littlewood–Paley decomposition, we obtain results on solvability, regularity, and symmetry of solutions.

Structural Shape Optimization Based on Multi-Patch Weakly Singular IGABEM and Particle Swarm Optimization Algorithm in Two-Dimensional Elastostatics

In this paper, a multi-patch weakly singular isogeometric boundary element method (WSIGABEM) for two-dimensional elastostatics is proposed. Since the method is based on the weakly singular boundary integral equation, quadrature techniques, dedicated to the weakly singular and regular integrals, are applied in the method. A new formula for the generation of collocation points is suggested to take full advantage of the multi-patch technique. The generated collocation points are essentially inside the patches without any correction. If the boundary conditions are assumed to be continuous in every patch, no collocation point lies on the discontinuous boundaries, thus simplifying the implementation. The multi-patch WSIGABEM is verified by simple examples with analytical solutions. The features of the present multi-patch WSIGABEM are investigated by comparison with the traditional IGABEM. Furthermore, the combination of the present multi-patch WSIGABEM and the particle swarm optimization algorithm results in a shape optimization method in two-dimensional elastostatics. By changing some specific control points and their weights, the shape optimizations of the fillet corner, the spanner, and the arch bridge are verified to be effective.

A 2.5D hybrid SBM-MFS methodology for elastic wave propagation problems

This paper proposes a novel hybrid methodology that combines the singular boundary method (SBM) and the method of fundamental solutions (MFS) for the computational simulation of elastic wave propagation. Particularly, the methodology aims to address radiation or scattering problems of longitudinally invariant systems in the wavenumber–frequency domain involving boundaries with intricate geometries. The approach uses the SBM to deal with the complex parts of these geometries and the MFS for the smooth ones. The method is studied in the framework of three case studies involving longitudinally infinite cavities in a homogeneous full-space with circular, square-shaped and five-cusped hypocycloid cross-sections. These three examples are selected to assess the accuracy and robustness of the hybrid SBM-MFS approach in comparison with alternative modelling strategies. The comparisons show that the proposed method inherits the accuracy of the MFS while keeping the robustness of the SBM when dealing with complex geometries. The method is found to be computationally more efficient than the SBM or the boundary element method (BEM). Moreover, the hybrid approach naturally mitigates the effect of fictitious eigenfrequencies, a feature that neither conventional versions of the SBM nor the BEM have.

A novel hybrid SBM-MFS methodology for acoustic wave propagation problems

In this paper, a novel hybrid meshless approach that combines the singular boundary method (SBM) and the method of fundamental solutions (MFS) to deal with two-dimensional (2D) exterior acoustic wave propagation problems is proposed and studied. The methodology is particularly devised to solve problems with complex boundary geometries containing geometric singularities such as corners and sharp edges. It employs the SBM to model intricate segments of these geometries and the MFS for the smooth ones. The proposed hybrid SBM-MFS method is studied in a 2D context in the framework of three benchmark examples involving acoustic radiation problems of circular-, square- and L-shaped objects in a full-space acoustic medium. In addition, the applicability of the proposed hybrid SBM-MFS methodology to predict the acoustic performance of a T-shaped thin barrier is also investigated. These examples are specifically designed to assess the feasibility, validity and accuracy of the hybrid SBM-MFS approach in comparison with the available analytical solutions and alternative numerical strategies such as the MFS, the SBM and the boundary element method (BEM). Numerical simulations demonstrate that the proposed method can match the level of accuracy of the MFS while keeps the robustness of the SBM when dealing with complex geometries, overcoming one the most important drawbacks of the traditional MFS. Moreover, the proposed hybrid SBM-MFS naturally avoids the non-uniqueness problem arising at the fictitious eigenfrequencies associated with the corresponding interior problems, a feature that neither the SBM nor the BEM inherently possesses.

Nonlinear coupled multi-physics field analysis of permanent magnet synchronous motor combining finite element analysis with optimized meshless method

Aiming at significant disadvantages of finite element analysis (FEA) performed in coupled multiphysics field calculation, including high computing cost and strong dependency on meshing quality. Taking a high-speed permanent magnet synchronous motor (P MSM) as an example, an approach combining FEA with optimized meshless method is proposed to solve the heat transfer problem. Firstly, electromagnetic losses under rated and 1.2 times overloaded operation conditions are calculated using FEA respectively. Secondly, three-dimensional FEA flow models considering the influence of viscosity varying with temperature are established. Thirdly, convective heat transfer coefficients are calculated as the boundary condition, after which the electromagnetic losses as heat sources are mapped to the threedimensional steady-state temperature field calculations using FEA and conventional meshless method respectively. Applying the above two methods, a singular boundary method (SBM) is employed to replace the temporal derivative term of the transient calculation equation. Finally, to mitigate the illconditioned matrix appearing in the conventional meshless method coupled with FEA, an optimized meshless method is proposed to recalculate the heat transfer equation. After that, the calculation results of the above methods are compared with experimental results. This research is beneficial for the multiphysical coupling analysis of electrical equipment with complex structures considering nonlinear factors.

Singular solutions for space-time fractional equations in a bounded domain

This paper is devoted to describing a linear diffusion problem involving fractional-in-time derivatives and self-adjoint integro-differential space operators posed in bounded domains. One main concern of our paper is to deal with singular boundary data which are typical of fractional diffusion operators in space, and the other one is the consideration of the fractional-in-time Caputo and Riemann–Liouville derivatives in a unified way. We first construct classical solutions of our problems using the spectral theory and discussing the corresponding fractional-in-time ordinary differential equations. We take advantage of the duality between these fractional-in-time derivatives to introduce the notion of weak-dual solution for weighted-integrable data. As the main result of the paper, we prove the well-posedness of the initial and boundary-value problems in this sense.