Accuracy Of Approximate Solution Research Articles

A novel multilevel finite element method for a generalized nonlinear Schrödinger equation

In this article, we focus on an efficient multilevel finite element method to solve the time-dependent nonlinear Schrödinger equation which is one of the most important equations of mathematical physics. For the time derivative, we adopt implicit schemes including the backward Euler method and the Crank–Nicolson method. Based on these stable implicit schemes, the proposed method requires solving a nonlinear elliptic problem at each time step. For these nonlinear elliptic equations, a multilevel mesh sequence is constructed. At each mesh level, we first derive a rough approximation by correcting the approximation of the previous mesh level in a special correction subspace. The correction subspace is composed of a coarse finite element space and an additional approximate solution derived from the previous mesh level. Next, we only need to solve a linearized elliptic equation by inserting the rough approximation into the nonlinear term. Then, we derive an accurate approximate solution by performing the aforementioned solving process on the multilevel mesh sequence until we reach the final mesh level. Owing to the special construct of the correction subspace, we derive a multilevel finite element method to solve the nonlinear Schrödinger equation for the first time, and meanwhile we also derive an optimal error estimate with linear computational complexity. Additionally, unlike the existing multilevel methods for nonlinear problems, that typically require bounded second-order derivatives of the nonlinear terms, the nonlinear term in our study requires only one-order derivatives. Numerical results are provided to support our theoretical analysis and demonstrate the efficiency of the presented method.

Novel Approach by Shifted Fibonacci Polynomials for Solving the Fractional Burgers Equation

This paper analyzes a novel use of the shifted Fibonacci polynomials (SFPs) to treat the time-fractional Burgers equation (TFBE). We first develop the fundamental formulas of these polynomials, which include their power series representation and the inversion formula. We establish other new formulas for the SFPs, including integer and fractional derivatives, in order to design the collocation approach for treating the TFBE. These derivative formulas serve as tools that aid in constructing the operational metrics for the integer and fractional derivatives of the SFPs. We use these matrices to transform the problem and its underlying conditions into a system of nonlinear equations that can be treated numerically. An error analysis is analyzed in detail. We also present three illustrative numerical examples and comparisons to test our proposed algorithm. These results showed that the proposed algorithm is advantageous since highly accurate approximate solutions can be obtained by choosing a few terms of retained modes of SFPs.

An Accurate Calculation Method on Blasingame Production Decline Model of Horizontal Well with Dumbbell-like Hydraulic Fracture in Tight Gas Reservoirs

Blasingame production decline is an effective method to obtain permeability and single-well controlled reserves. The accurate Blasingame production decline curve needs an accurate wellbore pressure approximate solution of the real-time domain. Therefore, the aim of this study is to present a simple and accurate wellbore pressure approximate solution and Blasingame production decline curves calculation method of a multi-stage fractured horizontal well (MFHW) with complex fractures. A semi-analytical model of MFHWs in circle-closed reservoirs is presented. The wellbore pressure and dimensionless pseudo-steady productivity index JDpss (1/bDpss) are verified with a numerical solution. The comparison result reaches a good match. Wellbore pressure and Blasingame production decline curves are used to analyze parameter sensitivity. Results show that when the crossflow from matrix to natural fracture appears after the pseudo-state flow regime, the value of the inter-porosity coefficient has an obvious influence on the pressure approximate solution of the pseudo-steady flow regime in naturally fractured gas reservoirs. The effects of relevant parameters on wellbore pressure and the Blasingame decline curve are also analyzed. The method of pseudo-steady productivity index JDpss can applied to all well and reservoir models.

Semi-analytic PINN methods for boundary layer problems in a rectangular domain

Singularly perturbed boundary value problems pose a significant challenge for their numerical approximations because of the presence of sharp boundary layers. These sharp boundary layers are responsible for the stiffness of solutions, which leads to large computational errors, if not properly handled. It is well-known that the classical numerical methods as well as the Physics-Informed Neural Networks (PINNs) require some special treatments near the boundary, e.g., using extensive mesh refinements or finer collocation points, in order to obtain an accurate approximate solution especially inside of the stiff boundary layer. In this article, we modify the PINNs and construct our new semi-analytic SL-PINNs suitable for singularly perturbed boundary value problems. Performing the boundary layer analysis, we first find the corrector functions describing the singular behavior of the stiff solutions inside boundary layers. Then we obtain the SL-PINN approximations of the singularly perturbed problems by embedding the explicit correctors in the structure of PINNs or by training the correctors together with the PINN approximations. Our numerical experiments confirm that our new SL-PINN methods produce stable and accurate approximations for stiff solutions.

A unified theory of the self-similar supersonic Marshak wave problem

We present a systematic study of the similarity solutions for the Marshak wave problem in the local thermodynamic equilibrium (LTE) diffusion approximation and in the supersonic regime. Self-similar solutions exist for a temporal power law surface temperature drive and a material model with power law temperature dependent opacity and energy density. The properties of the solutions in both linear and nonlinear conduction regimes are studied as a function of the temporal drive, opacity, and energy density exponents. We show that there exists a range of the temporal exponent for which the total energy in the system decreases, and the solution has a local maxima. For nonlinear conduction, we specify the conditions on the opacity and energy density exponents under which the heat front is linear or even flat and does possess its common sharp characteristic; this characteristic is independent of the drive exponent. We specify the values of the temporal exponents for which analytical solutions exist and employ the Hammer–Rosen perturbation theory to obtain highly accurate approximate solutions, which are parameterized using only two numerically fitted quantities. The solutions are used to construct a set of benchmarks for supersonic LTE radiative heat transfer, including some with unusual and interesting properties such as local maxima and non-sharp fronts. The solutions are compared in detail to implicit Monte Carlo and discrete-ordinate transport simulations as well gray diffusion simulations, showing a good agreement, which highlights their usefulness as a verification test problem for radiative transfer simulations.

Explicit and approximate series solutions for nonlinear fractional wave-like differential equations with variable coefficients

This work investigates the exact and accurate approximate series solutions of the strongly nonlinear wave-like differential equations of fractional order with variable coefficients utilizing the Laplace residual power series technique in the 2D, and 3D-space. The posed model is a notable model in mathematical physics, engineering physics, and biophysics, to characterize the variations in laser light intensity, and fluid particle velocity distributions in turbulent flows. In the meaning of Caputo fractional derivative, the posed models are solved analytically via simulation of the generalized Taylor formula in the Laplace space. Different from the other analytic methods, an easy recurrence formula for finding the Laplace power series coefficients is derived. To validate the proposed approach, it is applied to three interesting problems. The approximate analytical solutions of studied fractional problems are expressed in fast-converging Caputo-fractional Maclaurin formulas. The results are numerically and graphically analyzed to demonstrate the efficiency and applicability of the technique, as well as to examine the effects of partial arrangement on the behavior of the solutions. Ultimately, the results and numerical simulations indicate that the used technique is an effective approach and accurate tool for solving fractional evolution models and other fractional models arising in nonlinear phenomena.

A novel asymptotically consistent approximation for integral evaporation from a spherical cap droplet

The total evaporation rate due to a volatile capillarity-dominated droplet diffusively evaporating into the surrounding gas is a critically important quantity in industrial and engineering applications such as Q/OLED screen manufacturing. However, the analytical expression in terms of integrals in toroidal coordinates can be unwieldy in applications, as well as expensive to compute. Therefore, simple yet highly accurate approximate solutions are frequently used in practical settings. Herein we present a new approximate form that is both accurate and fast to compute, but also retains the correct asymptotic behaviour in the key physical regimes, namely hydrophilic and superhydrophobic substrates, and a hemispherical droplet. We illustrate this by comparison to several previous approximations and, in particular, illustrate its use in calculating droplet lifetimes, as well as approximating the local evaporative flux.

Spectral collocation with generalized Laguerre operational matrix for numerical solutions of fractional electrical circuit models

In this paper, we introduce a pioneering numerical technique that combines generalized Laguerre polynomials with an operational matrix of fractional integration to address fractional models in electrical circuits. Specifically focusing on Resistor-Inductor ($RL$), Resistor-Capacitor ($RC$), Resonant (Inductor-Capacitor) ($LC$), and Resistor-Inductor-Capacitor ($RLC$) circuits within the framework of the Caputo derivative, our approach aims to enhance the accuracy of numerical solutions. We meticulously construct an operational matrix of fractional integration tailored to the generalized Laguerre basis vector, facilitating a transformation of the original fractional differential equations into a system of linear algebraic equations. By solving this system, we obtain a highly accurate approximate solution for the electrical circuit model under consideration. To validate the precision of our proposed method, we conduct a thorough comparative analysis, benchmarking our results against alternative numerical techniques reported in the literature and exact solutions where available. The numerical examples presented in our study substantiate the superior accuracy and reliability of our generalized Laguerre-enhanced operational matrix collocation method in effectively solving fractional electrical circuit models.

Approximate and Parametric Solutions to SIR Epidemic Model

This article provides a detailed exploration of the SIR epidemic model, starting with its meticulous formulation. The study employs a novel approach called the upper and lower bounds technique to approximate the solution to the SIR model, providing insights into the dynamic interplay between susceptible S, infected I, and recovered R populations. A new parametric solution to this model has been presented. Applying the Adomian decomposition method (ADM) allows for the attaining of highly accurate approximate solutions in the context of the SIR epidemic model. To validate the accuracy and robustness of the proposed approach, a numerical exploration is conducted, considering a diverse range of experimental parameters. This numerical analysis provides valuable insights into the sensitivity and responsiveness of the SIR epidemic model under varying conditions, contributing to the broader understanding of infectious disease dynamics. The interplay between theoretical formulation and numerical exploration establishes a comprehensive framework for studying the SIR model, with implications for refining our ability to predict and manage the spread of infectious diseases.

A New Method for Solving a Neutral Functional-Differential Equation with Proportional Delays

This study presents and implements a new hybrid technique that combines the Sawi transform (ST) and Homotopy perturbation method (HPM) to solve neutral functional-differential equations with proportional delays. Some of the important properties of the method are established and validated. We start the method by first applying ST to obtain the recurrence relation. We, next, implement HPM to find convergent series solutions of the recurrence relation. The series is free of assumptions and restrictions, highlighting its adaptability and robustness. Moreover, the convergence of the method is established through convincing proof. To demonstrate its effectiveness and applicability, we provide five examples. The method yields accurate approximate solutions, or in some cases exact solutions, with a few number of iterations, reinforcing its reliability and validity. Moreover, the performance of the method is compared with some available methods and demonstrates its superiority and efficiency.

Debonding of an elastic layer with a cavity from a rigid substrate caused by rotation of a bonded rigid cylinder

Debonding of an elastic layer with a circular cylindrical cavity 0 ≤ ζ ≤ h , ρ ≥ 1 , from a rigid substrate under action of a rigid cylinder is the object of this study. The annular debonding zone 1 ≤ ρ < 1 + η is caused by rotation of a cylinder bonded to the cavity surface. The problem is reformulated as dual integral equations with Weber integral transforms kernels. A Volterra operator transforming a Weber transforms kernel into a Bessel function of the first kind and Hankel integral transforms allow us to reduce dual equations to a Fredholm integral equation of the second kind and then, by some transformation, to another Fredholm integral equation which is more suitable for approximate methods. As h ≥ 0 . 5 and 0 < η ≤ 1 , a highly accurate analytic approximate solution of the problem is suggested. The asymptotic solution of the problem is obtained as the width of debonding zone is very small while the thickness is not small. When the thickness is small, the Fredholm integral equations are computationally inefficient. A new method based on an operator transforming a Bessel function of the first kind into the kernel of Mehler–Fock integral transforms enabled us to convert one of the above-mentioned Fredholm equations into an equivalent Fredholm integral equation of the second kind that is effective for a small thickness. The asymptotic solution of the problem is obtained when both the layer thickness and the debonding zone width are small. Accurate mathematical methods, in particular investigations, and transformations of equations, developed in this study can be interesting to researchers employing dual integral equations technique in problems of mechanics and mathematical physics with mixed boundary conditions.

Numerical solution of unsteady elastic equations with C-Bézier basis functions

&lt;abstract&gt;&lt;p&gt;In this paper, the finite element method is applied to solve the unsteady elastic equations, C-Bézier basis functions are used to construct the shape function spaces, the semi-discrete scheme of the unsteady elastic equations is obtained by Galerkin finite element method and then the fully discretized Galerkin method is obtained by further discretizing the time variable with $ \theta $-scheme finite difference. Furthermore, for several numerical examples, the accuracy of approximate solutions are improved by 1–3 order-of magnitudes compared with the Lagrange basis function in $ L^\infty $ norm, $ L^2 $ norm and $ H^1 $ semi-norm, and the numerical examples show that the method proposed possesses a faster convergence rate. It is fully demonstrated that the C-Bézier basis functions have a better approximation effect in simulating unsteady elastic equations.&lt;/p&gt;&lt;/abstract&gt;

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;

A new approach in handling one-dimensional time-fractional Schrödinger equations

&lt;abstract&gt; &lt;p&gt;Our aim of this paper was to present the accurate analytical approximate series solutions to the time-fractional Schrödinger equations via the Caputo fractional operator using the Laplace residual power series technique. Furthermore, three important and interesting applications were given, tested, and compared with four well-known methods (Adomian decomposition, homotopy perturbation, homotopy analysis, and variational iteration methods) to show that the proposed technique was simple, accurate, efficient, and applicable. When there was a pattern between the terms of the series, we could obtain the exact solutions; otherwise, we provided the approximate series solutions. Finally, graphical results were presented and analyzed. Mathematica software was used to calculate numerical and symbolic quantities.&lt;/p&gt; &lt;/abstract&gt;

Collocation method with Lagrange polynomials for variable‐order time‐fractional advection–diffusion problems

In this study, we proposed a numerical technique to solve a class of variable‐order time‐fractional advection–diffusion equations (VOTFADEs) by applying an operational matrix of differentiation based on fractional‐order Lagrange polynomials (FOLPs). The variable‐order fractional derivative is assumed to be Caputo's derivative. Using the operational matrix and collocation method, the advection–diffusion equation can be reduced to an algebraic system of equations that can be solved using Newton's iterative method. Error analysis also has been carried out for the proposed method. The current approach is simple to use and computer oriented and provides highly accurate approximate solutions. The effectiveness and accuracy of the proposed method are demonstrated using a few numerical examples.

Mixed Fractional-Order and High-Order Adaptive Image Denoising Algorithm Based on Weight Selection Function

In this paper, a mixed-order image denoising algorithm containing fractional-order and high-order regularization terms is proposed, which effectively suppresses the staircase effect generated by the TV model and its variants while better preserving the edges and details of the image. Adding different regularization penalties in different regions is fundamental to improving the denoising performance of the model. Therefore, a weight selection function is designed using the structure tensor to achieve a more effective selection of regularization terms in different regions. In each iteration, the regularization parameters are adaptively adjusted according to the Morozov discrepancy principle to promote the performance of the algorithm. Based on the primal–dual theory, the original algorithm is improved by using the predictor–corrector scheme to obtain a more accurate approximate solution while ensuring the convergence of the algorithm. The effectiveness of the proposed algorithm is demonstrated through simulation experiments.

Modified Regularized Long Wave Denkleminin Nümerik Çözümü İçin Kuintik Trigonometrik B-spline Algoritması

This study introduces a new numerical algorithm for solving the modified regularized long-wave (MRLW) equation. To discretize the spatial variables and their derivatives, the collocation technique with quintic trigonometric B-spline functions is utilized and for the temporal derivative, the Adam's Moulton scheme is implemented. The performance and efficiency of the computational algorithm is tested on sample problem including the motion of single solitary wave. The error norm and three conservation constants are computed and compared with some of those available in the literature. The computed results verify that the suggested algorithm has the advantage in obtaining a highly accurate approximate solution of the MRLW equation as compared to the existing methods. The advantage of the method is that it is easy to implement and requires the low computational cost.

Analysis of the Vibro-Impact Nonlinear Damped and Forced Oscillator in the Dynamics of the Electromagnetic Actuation

In this work, the effect of vibro-impact nonlinear, forced, and damped oscillator on the dynamics of the electromagnetic actuation (EA) near primary resonance is studied. The vibro-impact regime is given by the presence of the Hertzian contact. The EA is supplied by a constant current generating a static force and by an actuation generating a fast alternative force. The deformations between the solids in contact are supposed to be elastic and the contact is maintained. In this study, a single degree of freedom nonlinear damped oscillator under a static normal load is considered. An analytical approximate solution of this problem is obtained using the Optimal Auxiliary Functions Method (OAFM). By means of some auxiliary functions and introducing so-called convergence-control parameters, a very accurate approximate solution of the governing equation can be obtained. We need only the first iteration for this technique, applying a rigorous mathematical procedure in finding the optimal values of the convergence-control parameters. Local stability by means of the Routh-Hurwitz criteria and global stability using the Lyapunov function are also studied. It should be emphasized that the amplitude of AC excitation voltage is not considered much lower than bias voltage (in contrast to other studies). Also, the Hertzian contact coupled with EA is analytically studied for the first time in the present work. The approximate analytical solution is determined with a high accuracy on two domains. Local stability is established in five cases with some cases depending on the trace of the Jacobian matrix and of the discriminant of the characteristic equation. In the study of global stability, the estimate parameters which are components of the Lyapunov function are given in a closed form and a graphical form and therefore the Lyapunov function is well-determined.

A Novel Solution Approach for Time-Fractional Hyperbolic Telegraph Differential Equation with Caputo Time Differentiation

In the current analysis, a specific efficient and applicable novel solution approach, based on a fractional power series technique and Laplace transform operator, is considered to predict certain accurate approximate solutions (ASs) for a time-fractional hyperbolic telegraph equation by aid of time-fractional derivatives in a Caputo sense. The solutions are obtained in a fractional Maclurian series formula by solving the original problem in the Laplace space aided by a limit concept having fewer small iterations than the classical fractional power series technique. To confirm applicability and feasibility of the proposed approach, three appropriate initial value problems are considered. Consequently, some simulations of gained outcomes are numerically and graphically implemented to support the effect of the fractional-order parameter on the geometric behavior of the obtained solutions. In addition, graphical representations are also fulfilled to verify the convergence analysis of the fractional series solutions of the classical solution. The proposed technique is therefore proposed to be a straightforward, accurate and powerful approach for handling varied time-fractional models in various physical phenomena.

An algebraically stabilized method for convection–diffusion–reaction problems with optimal experimental convergence rates on general meshes

Algebraically stabilized finite element discretizations of scalar steady-state convection–diffusion–reaction equations often provide accurate approximate solutions satisfying the discrete maximum principle (DMP). However, it was observed that a deterioration of the accuracy and convergence rates may occur for some problems if meshes without local symmetries are used. The paper investigates these phenomena both numerically and analytically and the findings are used to design a new algebraic stabilization called Symmetrized Monotone Upwind-type Algebraically Stabilized (SMUAS) method. It is proved that the SMUAS method is linearity preserving and satisfies the DMP on arbitrary simplicial meshes. Moreover, numerical results indicate that the SMUAS method leads to optimal convergence rates on general simplicial meshes.