Semi-Analytical Solutions for Fuzzy Caputo–Fabrizio Fractional-Order Two-Dimensional Heat Equation

In the analysis in this article, we developed a scheme for the computation of a semi-analytical solution to a fuzzy fractional-order heat equation of two dimensions having some external diffusion source term. For this, we applied the Laplace transform along with decomposition techniques and the Adomian polynomial under the Caputo–Fabrizio fractional differential operator. Furthermore, for obtaining a semi-analytical series-type solution, the decomposition of the unknown quantity and its addition established the said solution. The obtained series solution was calculated and approached the approximate solution of the proposed equation. For the validation of our scheme, three different examples have been provided, and the solutions were calculated in fuzzy form. All the three illustrations simulated two different fractional orders between 0 and 1 for the upper and lower portions of the fuzzy solution. The said fractional operator is nonsingular and global due to the presence of the exponential function. It globalizes the dynamical behavior of the said equation, which is guaranteed for all types of fuzzy solution lying between 0 and 1 at any fractional order. The fuzziness is also included in the unknown quantity due to the fuzzy number providing the solution in fuzzy form, having upper and lower branches.


Introduction
The area of fractional calculus has attained considerable attention in the last three decades. Famous scientists have provided their contributions on this aspect by introducing different fractional operators in different articles. Modern calculus has provided more realistic results than classical calculus. It has described the dynamics of different realworld phenomena lying between two integers. Furthermore, the fractional operators have more degrees of freedom, and they generalize the integer differential operators. To date, various researchers have published more research articles, more books and different monographs that touch upon the said area. Podlubny provides a physical and geometrical explanation of the fractional-order derivatives [1]. An analysis of various dynamical systems in the sense of fractional-order operators can be seen in [2]. The applications of the said calculus in physics may be studied in [3]. Some fuzzy fractional-order linear and non-linear dynamical problems have been analyzed for semi-analytical solutions using fractional Sumudo transform [4,5]. Many types of publications have also been based on existence, uniqueness and numerical analysis under fractional-order concepts.
Modern calculus can be extended to various fields of physical and natural sciences. These operators may be applicable in the field of applied and pure mathematics. Now, we will apply these operators to those fields where fuzziness lies in the data. Uncertainty may be found in quantum mechanics, chemistry and different distributions of chromosomes in humans or other living things. The concepts of fuzziness were first introduced by Zadeh [6] in 1965 for sets by defining the membership function. The idea of uncertainty was then applied in other areas, i.e., fuzzy topology, the fixed-point theorem, fuzzy automata, control systems, etc. Zadeh, Chang and a few other scientists have applied this idea of fuzzy sets to define some membership functions and fuzzy control [7]. By applying the idea of fuzziness and its operators, various scholars have published basic fuzzy differential calculus [8]. Elementary fuzzy calculus and the study of fuzzy differential equations may be seen in [9,10]. From 2001, differential equations have been presented in fuzzy form by taking account of uncertainty in the initial data. Some scholars such as "Dobius" and "Prada" have defined the fundamental concepts of integral equations in fuzzy form [11]. We deal with all those problems having uncertainty or vagueness in their information using fuzzy differential and integral equations or systems. Therefore, plenty of research work on such types of fuzzy fractional differential equations (FDEs) has been presented by many researchers [12][13][14].
Semi-analytical investigations for the solution of fuzzy FDEs, important techniques and schemes have been introduced in the past. Among them are the well-known integral transforms of Fourier, Natural, Z, Laplace, Sumudu, etc. [4,15,16]. Some of the iterative and series solution methods such as the "homotopy and improved homotopy method" [17][18][19], Adomian decomposition along with Adomian polynomial and Laplace Adomian decomposition methods (LADM) [20,21], "Taylor's" series method, etc. have also been used for dealing with such problems [22,23]. However, to the best of our knowledge, the aforementioned techniques have not been properly used to study fuzzy fractional partial differential equations (FPDEs). The idea of a nonsingular fractional derivative is new and has very recently received attention. Here, we remark that, recently, some valuable work related to applications of the nonsingular kernel type derivatives CF and ABC have been considered; see [24][25][26][27][28][29].
As far as novelty is concerned, we consider an initial value heat equation with external source terms under the fuzzy Caputo fractional operator. We take the fuzziness in an unknown quantity and initial conditions to find its approximate solution in fuzzy form with two branches, because of the fuzzy number. Different scientists have worked on diffusion equations, of both integer and fractional order; some have worked on the fuzzy heat equation. We consider both fuzziness and fractional order for the analysis of the said equation having different external source terms [12][13][14][30][31][32]. Therefore, a fuzzy solution is required, as many quantities can be obtained in fuzzy format when there is some ambiguity regarding the exact quantity. Keep in mind the fact that the considered method has the ability to produce series-type solutions for such problems under consideration. Furthermore, this method decreases the error and avoids wasting time. Additionally, this method is free from any axillary parameters that run the whole system. The adopted method can be used as a robust tool to analyze fuzzy problems under the mentioned fractional operator. Therefore, here, in the current work, we compute a semi-analytical series solution or approximate solution using analytical techniques applying Laplace Adomian decomposition under the fuzzy Caputo-Fabrizio derivative concept for the proposed 2D heat equation given as where CF D θ τ represents the Caputo-Fabrizio fractional operator while Generally, the transfer of heat in a thin rod from one point to the other can be formulated by the boundary value problem of the one-dimensional heat equation. The two-dimensional heat equation represents the transfer of heat through an infinite thin sheet. Here, U represents the temperature of the body at any point in the thin sheet. The phenomena of heat transfer may be found in many diffusion problems. Therefore, the investigation of partial differential equations such as the heat equation has much more application in real life [33,34]. The analysis for semi-analytical solutions for the problems of electrical circuits [35] has been performed using numerical approximation. The use of different analytical techniques such as Laplace transform may be found in [36]. The technique of double Laplace transform for waves and Laplace equations was conducted in [37]. The phenomenon of heat conduction was discussed in [38]. The semi-analytical solution, via the Laplace Adomian decomposition method, is provided in [39,40]. The semi-analytical solution of the 2D heat equation was found using double Laplace transform without an external diffusion term F in [37]. We analyze the 2D heat equation by applying the external term under the fuzzy concept. Different examples along with a numerical simulation were used to perform the verification of the theoretical results.

Background Materials
In the following part, we present some theorems, lemmas and definitions of classical calculus [1,2] and some well-known fractional operators from [11,41].

Definition 1. Take a continuous fuzzy operator
where is the space of the fuzzy operator along with continuity while L F [0, b] is the space of Lebesgue fuzzy integrable operators, respectively, then the non-integer-order Caputo-Fabrizio fuzzy integration is given as:

Definition 2. Similarly, for an operator
the fractional-order Caputo-Fabrizio (CF) differential operator in the fuzzy sense is given as here where the integration exists or converges and m = θ + 1. As θ lies in the interval (0, 1], m = 1. Definition 3. For x, either real or complex, the transform of Laplace for the fuzzy function F (x) is given as where 0 < β.
is called a fuzzy number if the following conditions hold: (i) k is continuous up to a peak value;

k is normal;
(iv) cl{y ∈ R, k(y) > 0} is bounded and continuous, where cl represents closure for the support of y.
We denote the set of fuzzy numbers collectively as E .
If k(r) = k(r) = r, then r is the crisp number.

Main Work
In this section, we investigate our proposed model for a semi-analytical solution. For this, we use the Laplace transform of the Caputo-Fabrizio fractional differential operator along with decomposition techniques for (1) as where θ is in (0, 1]; therefore, the Laplace transform of (7) is On using the initial condition, we obtain Decompose the solution as U(τ, Take parts of the solution by choice of comparison as Taking the inverse Laplace transform, we obtain n ≥ 0.
Thus, the solution becomes Equation (12) is the solution in series form, whose convergence was also derived in [42].

Examples
Now, in this part, we provide some illustrations of 2D fuzzy arbitrary-order heat equations with external source terms by applying the aforementioned method for obtaining the approximate results. Example 1. The first example is of a fuzzy heat equation of fractional order having the external source term F (τ, x 2 , x 1 ) = (x 2 + x 1 + 1) along with a fuzzy initial condition as in [12][13][14], as follows: Using the scheme of Equation (11), we obtain In a similar fashion, we can obtain the higher terms. The series solution is obtained using Equation (12); therefore, we write

Example 2.
Take another illustration of the fuzzy fractional heat equation having the external source term F (τ, x 2 , x 1 ) = τ 2 + x 2 + x 1 with an initial approximation as in [12][13][14] as follows: . Upon the application of the Laplace transform of the Caputo-Fabrizio derivative, one may obtain From (12), we can write U(τ, x 2 ,   Example 3. Let us consider the last example of a non-integer-order heat equation having uncertainty in the initial condition along with the source term F (τ, x 2 , x 1 ) = τ 4 + x 2 + x 1 like in [12][13][14] as follows: Using the scheme of Laplace transform and by the decomposition of U(τ, x 2 , x 1 ) into small terms, we obtain a hierarchy of terms as follows: Now, we can write all the decomposed terms into sum form as In Figure 5, we have drawn a 3D plot showing the upper and lower regions of the fuzzy solution for the fractional-order fuzzy heat equation at x 1 = 0.8, x 2 = 0.5, and at different fractional orders of θ. Figure 6 shows

Conclusions and Discussion
Finally, we conclude that a successful scheme for the computation of semi-analytical or approximate solutions has been used for two-dimensional fuzzy fractional-order diffusion or heat partial differential equations having diffusion terms as external terms. The techniques of Laplace Adomian decomposition (LADM) were applied for the analysis of a semi-analytical solution in the form of a series solution and may be used for different fractional-order problems. We also verified our general techniques of solution using three specific illustrations. All of the three examples were analyzed using the said techniques by taking account of the uncertainty in the initial values by multiplying the fuzzy number at two different fractional orders. The obtained results for all the three examples were simulated with two different formats of 2D and 3D graphs. The achieved results were also obtained in the form of uncertainty, having upper branches and lower branches of fuzzy solutions. We may also take different values of fractional orders to check and compare them with the integer order. By increasing the order of θ, we can converge to the value of integer order 1. For the future, the implementation of this scheme may be applied to different partial differential equations for series-type solutions, for both integer and non-integer orders in fuzzy form. Furthermore, the adopted techniques can be used as a robust tool to study fractional-order non-linear and linear dynamical systems. Data Availability Statement: This is not applicable for our research work.