Fractal dynamics and computational analysis of local fractional Poisson equations arising in electrostatics

Jagdev Singh ,Hassan Kamil Jassim ,Devendra Kumar and Ved Prakash Dubey

1 Department of Mathematics,JECRC University,Jaipur,Rajasthan,India

2 Department of Mathematics,Kyung Hee University,26 Kyungheedae-ro,Dongdaemun-gu,Seoul,02447,Republic of Korea

3 Department of Computer Science and Mathematics,Lebanese American University,Beirut,Lebanon

4 Department of Mathematics,Faculty of Education for Pure Sciences,University of Thi-Qar,Nasiriyah,Iraq

5 Department of Mathematics,University of Rajasthan,Jaipur-302004,Rajasthan,India

6 Faculty of Mathematical and Statistical Sciences,Shri Ramswaroop Memorial University,Barabanki-225003,Uttar Pradesh,India

Abstract In this paper,the local fractional natural decomposition method(LFNDM)is used for solving a local fractional Poisson equation.The local fractional Poisson equation plays a significant role in the study of a potential field due to a fixed electric charge or mass density distribution.Numerical examples with computer simulations are presented in this paper.The obtained results show that LFNDM is effective and convenient for application.

Keywords: poisson equation,local fractional natural transform,adomian decomposition method,local fractional derivative,electrostatics,fractal media

1.Introduction

The Poisson equation (PE) plays a key role in the field of electrostatics,where it is solved to determine electric potential from a provided charge distribution.The Poisson equation is linear in potential and is the source term used to stipulate the object’s static electricity.Columb’s law and Gauss’s theorem derive the PE.The solution of the Poisson equation is actually a potential field subjected to a provided mass density distribution or electric charge and further,the determined potential field computes gravitational or electrostatic field.The Poisson equation models the phenomena of intersecting interfaces and electrodynamics [1,2].

The field of local fractional calculus (LFC) investigates the characteristics of physical models appearing in a fractal space.Fractals appear as random geometrical structures that do not show any change during amplification of their shapes.The use of fractal theory can be found in various areas of physical sciences along with important applications in electrostatics,quantum mechanics and high energy physics.As compared to classical derivatives,a local fractional derivative(LFD) provides more precise estimates of performance measures.Fractals elucidate those images that cannot be undertaken by Euclidean geometry and handle those objects that deal with dimensions of real order.The nature of fractionalorder modeling is nonlocal and therefore is irrelevant to deal with characteristics of local scaling phenomena or LF derivability.The LFD operator (LFDO) with fractal order is an efficient instrument for modeling physical phenomena and provides physical insights along with geometric observations.The motivational purpose of LFC is to explore the differential properties of extremely irregular &nondifferentiable functions.These LFDOs appeared as a generalized form of classical derivatives to fractal-order conserving local features of derivatives to explore local scaling properties of nowhere differentiable and extremely irregular functions [3].These reasons and features prompted motivation for the modeling of the Poisson equation with LFDO.The above-discussed features of LFDs show the significance of the chosen local fractional Poisson equations (LFPEs) and that obviously,the local fractional modeling of the Poisson equation is far better and superior compared to integer-and fractional-order modeling of the Poisson equation.

When physical variables in Poisson’s model are nondifferentiable functions defined on Cantor sets,the classical conservation law doesn’t fit and so the integer-order Poisson’s model becomes irrelevant in this sense.Therefore to deal with this difficulty,Chen et al [4] suggested the Poisson equation model with LFDO arising in electrostatistics within the LF conservation laws in the domain of LFC [5–7].The Poisson equation with LFDO was presented in [4] as follows:

subject to the initial and boundary conditions:

whereφ(μ,τ) is an unknown local fractional continuous nondifferentiable function,g(μ) is the given function,and the LFDO ofφ(μ) of orderεatμ=μ0is defined as

Many analytical and numerical techniques have been employed to obtain approximate solutions of local fractional partial differential equations(LFPDEs)for example,the local fractional function decomposition method [8,9],the local fractional Adomian decomposition technique (LFADT)[9–11],the local fractional series expansion technique [12,13],the local fractional Laplace transform approach[14],the local fractional variational iteration approach(LFVIA) [15–19],the local fractional reduced differential transform scheme (LFRDTS) [19],the local fractional homotopy analysis Sumudu transform method [20],the local fractional differential transform scheme [21,22],the local fractional Laplace VIA [23–28],the local fractional Laplace decomposition technique [29],the local fractional homotopy analysis scheme [30],the local fractional Laplace homotopy perturbation technique (LFLHPT) [31,32] and the local fractional natural homotopy perturbation technique [33].Recently,Dubey et al [34] and Kumar et al [35] presented fractal dynamics of LFPDEs occurring in physical sciences.Moreover,Dubey et al [36,37] also investigated the local fractional Tricomi equation and local fractional Klein–Gordon models in a fractal media using hybrid local fractional schemes.Recently,Alqhtani et al [38] discussed spatiotemporal chaos in spatially extended fractional dynamical systems.Moreover,Alqhtani et al [39] studied the chaotic Lorenz system.Srivastava et al[40]analyzed fractal–fractional Kuramoto–Sivashinsky and Korteweg–de Vries equations.

Recently,the LFPEs were studied by several authors using LFVIA [4],LFLHPT and LFRDTS [41].Moreover,Singh et al[42]and Li et al[43]also studied and investigated the LFPEs in fractal media.The purpose of this paper is to introduce a new method to obtain the analytical approximate solutions to the Poisson equation with LFDO.In this paper,we apply the local fractional natural decomposition method(LFNDM) for the solution of LFPE.The LFNDM actually appears as a coupling of the LFADT [9,44] and local fractional natural transform (LFNT) [33].The LFNDM was introduced earlier in [45].The main focus of the paper is to illustrate the implementation of LFNDM for different forms of LFPE and numerical simulations with the help of 3D graphic visuals on the cantor set.The LFPE explores the nature of the potential field in a fractal domain in view of nondifferentiable functions where a free charge occurs.The graphical analysis of the solution of LFPE provides significant physical characteristics of the LFPE in a fractal medium.

The application of the implemented method is shown by using two different examples and obtained solutions have also been compared with solutions computed by other methods in previous works.The 3D figures have been constructed for solutions of LFPE using MATLAB.The 3D plots depict the fractal nature of the functionφ(μ,τ).The mathematical analysis shows that the implemented hybrid approach is beneficial to obtain the solutions for LFPEs.To make the convergence of LFNDM faster,the LFNT is selected.Two examples of LFPE are solved to demonstrate the application of LFNDM.The resulting solutions obtained from the applied scheme actually enter the picture as a special case of the classical Poisson equation model when the fractal orderε=1,which shows the convergence of fractal geometry to Euclidean geometry.The coupling of LFADT with LFNT provides faster computations compared to the LFADT.Furthermore,this merger reduces the computational procedure as compared to other conventional methods while still giving reliable results.This work examines two significant aspects of the LFNDM.One feature is linked to the easy decomposition of non-linear quantities in a simple form by adopting Adomian polynomials and the other feature is related to the delivery of closed-form solutions in a series form with faster convergence.The original contributions of the paper are concerned with solution and numerical simulations for given LFPEs on the Cantor set via the applied hybrid method.It is also observed that the achieved solutions exactly match with previously reported solutions computed in the recent past.

The novelty and original contributions of the paper are applications of LFNDM to LFPEs which were constructed in[45].Thus the new application along with computer simulations on the Cantorian set surely articulates the novelty of this work.The LFNDM provides fast convergence series solutions.The implemented method is time-saving,more trustable and proficient in comparison to other approaches.Moreover,computer-based simulations are also provided for the computed solutions of different forms of LFPE for the integer orderε=1.0and the fractal orderof an LFD by using MATLAB.The remaining sections of the present work are developed as follows: In section 2,some essential fundamentals of LFC are displayed.In section 3,the basic procedure of LFNDM is provided.Applications of the LFNDM are demonstrated in section 4.Section 5 deals with numerical simulations for LFPEs.The conclusion of the present paper is reported in section 6.

2.Fundamentals of the LFC and LFNT

Definition 1[46,47].The LF derivative ofφ(μ) of orderεat the pointμ0is defined as

Definition 2[47].The Mittag-Leffler function in a fractal space is defined by

Definition 3[33].The LFNT of a functionφ(τ) of order 0<ε≤1is stated as

Following (2.3),its inverse formula is expressed as

wheresεandvεare the local natural transform variables,andγis a real constant.

The LFNT of some special functions reported in[33]are given as follows:

3.Analysis of LFNDM

To describe outlines of the procedure,the following partial differential equation (PDE) with LFDOs is considered

Operating the LFNT on equation (3.1) and using the property of the LFNT,we get

Applying the inverse formula of LFNT on equation(3.2),we acquire

Now,we represent the solution in an infinite series in this way:

whereAnis the local fractional Adomian polynomial and can be calculated using the following formula:

Using of equations(3.4)and(3.5)in equation(3.3)yields the following result:

On comparing both the sides of equation (3.6),we obtain

The local fractional recursive relation in its general form is obtained as

Hence,the approximation solution of equation (3.1) is given by

4.Applications of the LFNDM

Example 1.Consider the following LFPE described in [41]

with initial conditions:

Using equations (3.8) and (4.1),the LF iteration algorithm is reported in the following form:

From equation (4.3),we can find the following components

and so on.

Therefore,the approximate solutionφ(μ,τ) of equation (4.1) is given by

The solution (4.4) is in good agreement with the solutions computed by LFLHPT and LFRDTS [41].

Example 2.We report the following LFPE described in [41]

Making use of equations (3.8) and (4.5),the local fractional iteration algorithm is generated in this way:

From equation (4.7),we can generate the components as follows.

and so on.

Therefore,the solutionφ(μ,τ) of equation(4.5)is given by

Figure 3.3D plot of φ (μ,τ)for Example 2 with respect toμ and τ for ε=1.0.

The acquired solution (4.8) is exactly the same as the solutions obtained by LFLHPT and LFRDTS [41].

5.Numerical results and analysis

In this section,the computer-based simulations for solutions of LFPEs obtained by LFNDM are presented.The numerical investigation of LFPEs considers distinct values ofε=Figures 1 and 2 elucidate the 3D variation of solutionφ(μ,τ)for Example 1 forε=1.0andrespectively.Similarly,figures 3 and 4,respectively,demonstrate the nature ofφ(μ,τ)for Example 2 forε=1.0 andFigures 2 and 4,respectively,depict the variation ofφ(μ,τ)in a fractal domain for Examples 1 and 2.The graphic visuals for the solutionφ(μ,τ) demonstrate a fractal nature and obviously depict the nondifferentiability of the function in certain phases.

Figure 4.3D plot ofφ (μ,τ)for Example 2 with respect toμ and τ for ε=log 2/log 3.

6.Conclusion

In this work,we have considered the LFPE with LFDOs.The method which is called the LFNDM has been applied successfully to attain the approximate solutions for LFPEs.The results are obtained in the closed form of an infinite power series.The examples show that the outcomes of LFNDM are in good agreement with the results obtained by LFLHPT and LFRDTS.This work also depicts that the applied method is systematic and can be helpful for solving nonlinear and linear LFPDEs with fractal order.In future studies,other fractal order physical models can also be solved by the applied technique to attain new results and conclusions.

ORCID iDs