Next Article in Journal / Special Issue
On the Existence and Ulam Stability of BVP within Kernel Fractional Time
Previous Article in Journal
An Analysis and Global Identification of Smoothless Variable Order of a Fractional Stochastic Differential Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Fuzzy Fractional Partial Differential Equations in the Realm of Fractal Calculus with Local Fractional Derivatives

1
School of Mathematical Sciences, Zhejiang Normal University, Jinhua 321004, China
2
School of Computer Science and Technology, Zhejiang Normal University, Jinhua 321004, China
3
Faculty of Mathematical Sciences and Informatics, University of Khartoum, Khartoum 11111, Sudan
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2023, 7(12), 851; https://doi.org/10.3390/fractalfract7120851
Submission received: 28 September 2023 / Revised: 22 October 2023 / Accepted: 6 November 2023 / Published: 29 November 2023

Abstract

:
In this study, local fuzzy fractional partial differential equations (LFFPDEs) are considered using a hybrid local fuzzy fractional approach. Fractal model behavior can be represented using fuzzy partial differential equations (PDEs) with local fractional derivatives. The current methods are hybrids of the local fuzzy fractional integral transform and the local fuzzy fractional homotopy perturbation method (LFFHPM), the local fuzzy fractional Sumudu decomposition method (LFFSDM) in the sense of local fuzzy fractional derivatives, and the local fuzzy fractional Sumudu variational iteration method (LFFSVIM); these are applied when solving LFFPDEs. The working procedure shows how effective solutions for specific LFFPDEs can be obtained using the applied approaches. Moreover, we present a comparison of the local fuzzy fractional Laplace variational iteration method (LFFLIM), the local fuzzy fractional series expansion method (LFFSEM), the local fuzzy fractional variation iteration method (LFFVIM), and the local fuzzy fractional Adomian decomposition method (LFFADM), which are applied to obtain fuzzy fractional diffusion and wave equations on Cantor sets. To demonstrate the effectiveness of the used techniques, some examples are given. The results demonstrate the major advantages of the approaches, which are equally efficient and simple to use in order to solve fuzzy differential equations with local fractional derivatives.

1. Introduction

The concept of fuzzy theory [1] was originally introduced by Zadeh in 1965; it has been considered extensively from several different aspects of the theory and its applications, such as linguistic information systems and approximate reasoning in [2,3,4], fuzzy decision making and fuzzy logic in [5,6], fuzzy analysis in [7], and fuzzy topology in [8,9].
Many mathematicians are interested in fuzzy fractional differential equations (FFDEs) and fuzzy fractional calculus since these theories are helpful in determining uncertainty influenced by ambiguity and inaccuracy. The concepts of Riemann–Liouville, Caputo–Hadamard, Caputo–Katugampola, Caputo–Atangana–Baleanu, Caputo–Fabrizio derivatives, and Caputo fuzzy fractional integrals and/or derivative operators have been applied in the majority of articles presented to date on this topic (see [10,11,12,13,14,15,16,17,18,19,20,21,22]).
Local fractional calculus has garnered the interest of many scientists due to its efficient treatment of non-differentiable functions. Local fractional calculus (LFC) represents differential and integral functions on fractal sets in their generalized form. Recently, physicists and engineers have appreciated LFC in addition to mathematicians. Fractional calculus is a generalized version of traditional calculus that examines real-order integrals and derivatives. Particularly when inherent constraints on a system influence its dynamics, these types of fractional derivative operators model a wide range of real-world phenomena more accurately. The literature contains a multitude of formulations pertaining to fractional integrals and derivatives. Due to the nonlocal nature of these formulations, they cannot be used to investigate the characteristics of local scaling phenomena or local fractional differentiability. Kolwankar and Gangal [23] introduced the notion of local fractional derivatives a decade ago; these derivatives show several features of integer-order derivatives. However, they lose the memory properties implicit to the fractional-order derivatives. The local differential operators with fractal order have been applied as effective and powerful mathematical tools to simulate complex real-world problems, communicating physical observations and geometrical clarifications. The purpose of LFC is to investigate the differential properties of fractal objects and nowhere differentiable functions. The fractal properties are typically observed in aquifers, porous media, turbulence, and other phenomena.
The local fractional derivative and calculus theory have recently been introduced by Yang et al. [24]. This is based on fractal geometry and is the best method for the characterization of the non-differentiable function defined on Cantor sets by Yang and Hua [25]. Golmankhaneh et al. [26] used fractional calculus in generalized Newtonian mechanics, Maxwell’s equations, and Hamiltonian mechanics. This then led to the emergence of ordinary differential equations and PDEs relating to this new concept, which became known as local fractional differential equations (LFDEs) and local fractional PDEs. This prompted some researchers to use the aforementioned methods to solve this new type of equation, including the local fractional ADM (LFADM) by Yang et al. [27], the local fractional homotopy perturbation method (LFHPM) [28], the local fractional variational iteration method [25], the local fractional variational iteration transform method [29], the local fractional Laplace decomposition method [30], the Yang–Machado–Baleanu–Cattain wave method, the fractional sech function method [31], and the fractional sech function method [32].
In 1993, Watugala [33] introduced the Sumudu transformation, which was thought to be the first of its kind in this field. This transformation was used to solve problems in control engineering. For this transformation, we devised a difficult inverse formula [34]. Asiru employed the Sumudu transformation to solve integral equations and systems of discrete dynamic equations in his studies [35,36]. Belgacem and Karaballi [37] presented a comprehensive list of functions and properties of the Sumudu transformation.
Recently, scientists and engineers have applied the ADM in order to solve linear, nonlinear differential, and integral problems ([38,39,40,41,42]). Moreover, some researchers have considered the ADM with the Sumudu transform. Trushit and Ramakanta [43] used the Adomian decomposition Sumudu transform method (ADSTM). We investigated the fluctuations in the temperature distribution, efficiency, and efficacy of porous fins for different fractional orders, porous parameters, and convection parameters. Saadeh et al. [44] presented the modified double ARA–Sumudu decomposition method for the obtained PDEs.
J.H. He invented the VIM, which was successfully used to solve autonomous ODEs [45,46]. It has been proven that this methodology works well when dealing with a variety of issues. Similar to the way in which the Shehu transform method is used to modify this method, the updated methodology is known as the variational iteration transform method (VITM). This technique has been used to solve many problems; see [40,47,48]. Recently, many authors have presented the VIM with the Sumudu transform. Prakash et al. [49] proposed a new computational method that was used to solve the numerically nonlinear time-fractional Zakharov–Kuznetsov (FZK) equation in two dimensions. Anac et al. [50] presented approximations of solutions to random component time-fractional PDEs with Caputo derivatives using the new Sumudu transform iterative method (SUTIM).
The Laplace variational iteration method (LVIM) is a combination of the Laplace transform and VIM. Many scientific projects have considered this technique. Bhargava et al. [51] used the LVIM to obtain series solutions to fractional-order heat equations that appear in many engineering applications. Nadeem et al. [52] presented a new amendment to the LVIM for the solution of fourth-order parabolic PDEs with variable coefficients.
In 2015, Yang et al. [28] and Zhang et al. [53] proposed the LFHPM. To further obtain the solution for the local fractional Tricomi equation (LFTE), Singh et al. [54] proposed the local fractional homotopy perturbation Sumudu transform method (LFHPSTM). Furthermore, Zhao et al. [55] used LFHPSTM to look into the local fractional heat conduction equations in fractal media. Others [56,57,58] have conducted work based on the application of LFHPSTM. The present work focuses on using the local fuzzy fractional homotopy perturbation Sumudu transform method (LFFHPSTM) to solve a few LFFPDEs in the fractal domain. Numerical simulations based on computers have been carried out in order to clarify the basic structure of the physical models expressed by LFFPDEs. The local fractional Sumudu transform (LFST) and the LFHPM as discussed in [59] are combined to generate the LFHPSTM.
The contributions and novelty of this work can be summarized as follows:
  • The LFFHPSTM is clearly advantageous in comparison to the local fuzzy fractional homotopy perturbation method and LFADM as it combines two powerful mathematical tools to solve nonlinear LFFDEs. The combination of LFFHPM with the local fuzzy fractional Sumudu transform (LFFST) produces faster computations than LFFHPM; hence, this combination saves time. Furthermore, without requiring linearization, perturbation, or discretization, this method provides a convergent series solution. The LFHPSTM does not involve rounding errors and so consumes less time. Further, by using He’s polynomials to solve nonlinear terms, this approach can handle nonlinear LFDEs. The novelty and uniqueness of this work arise from the fact that the hybrid method used has never been applied to the studied LFFPDEs in the recent past. The technique has two features: the first is that it decomposes nonlinear terms using He’s polynomials, and the second is that it produces closed-form series solutions with fast convergence. Furthermore, the LFFHPSTM does not require the solution of complex Adomian polynomials. These attributes are decisive in the selection of this method to solve LFFPDEs.
  • For the purpose of solving linear LFFPDEs, the combination of the ADM and the Sumudu transform method in the context of the local fractional derivative proves to be rather successful. Utilizing the series form of the solution proposed by the algorithm, rapid convergence will occur towards the exact solution. It is evident from the findings that the LFFSDM produces very accurate solutions with a minimal number of iterations. Therefore, given the effectiveness and flexibility of the application, as demonstrated by the provided examples, the work concludes that the LFFSDM can be used for additional linear LFFPDEs of higher order.
  • The coupling of the VIM and the Sumudu transform method in the sense of local fractional derivatives has been shown to be highly effective in solving linear and nonlinear LFFPDEs. The local fuzzy fractional Sumudu variational iteration method (LFFSVIM) is a user-friendly solution for such problems. This method combines two potent techniques to obtain exact or approximate solutions to linear–nonlinear LFFPDEs. The modified LFFSVIM is an alternative algorithm to solve linear–nonlinear LFFPDEs.
  • The fuzzy diffusion and wave equations defined on Cantor sets under fractal conditions are solved using the LFFLVIM. The method is found to be both useful and efficient in analytical applications. The LFFSEM, LFFVIM, and LFFADM all provide solutions to the same set of problems. All four approaches yield similar results; hence, the LFFLVIM is used as a viable alternative to the more standard technique of obtaining approximate solutions to linear–nonlinear fuzzy fractional differential equations.
This work is structured as follows. Section 2 is dedicated to providing the necessary notations and fundamental definitions. In Section 3, we present some LFFPDEs that are treated using a hybrid local fuzzy fractional approach. In Section 4, we propose the comparison of some analytical techniques applied to obtain fuzzy fractional diffusion and wave equations on the Cantor set. Finally, the conclusions end the work.

2. Preliminaries

The set of fuzzy numbers can be denoted by E 1 , whereas normal, fuzzy convex, upper semi-continuous and compactly supported fuzzy sets can be defined on the real line. For 0 < λ 1 , set [ ψ ] λ = { ϑ R | ψ ( ϑ ) λ } , and [ ψ ] 0 = c l { ϑ R | ψ ( ϑ ) > 0 } . We explain [ ψ ] ϱ = [ ψ ̲ ϱ , ψ ¯ ϱ ] ; consequently, if ψ E 1 , the ϱ -level set [ ψ ] ϱ is a closed interval for all ϱ [ 0 , 1 ] , refs. [7]. Suppose that ψ , E 1 and k R , and the addition and scalar multiplication are defined by
  • [ ψ + ] ϱ = [ ψ ] ϱ + [ ] ϱ ,
  • [ k ψ ] ϱ = k [ ψ ] ϱ .
The triangular fuzzy number defined as a fuzzy set in E 1 , is determined by ψ = ( a , b , c ) R and a b c such that ψ ̲ ϱ = a + ( b a ) ϱ and ψ ¯ ϱ = c ( c b ) ϱ are the endpoints of ϱ -level sets for all ϱ [ 0 , 1 ] . The support of fuzzy number ψ is given as
sup p ( ψ ) = c l { ϑ R | ψ ( ϑ ) > 0 } ,
where c l is the closure of set { ϑ R | ψ ( ϑ ) > 0 } .
Definition 1
([60]). Let us consider ψ , E 1 . If there exists a w E 1 such that ψ = + w , then w is called the Hukuhara difference of ψ , , and it is denoted using ψ . Note that
ψ ψ + ( 1 ) .
The Hausdorff distance between two fuzzy numbers is defined as D : E 1 × E 1 R + { 0 }
D ( ψ , ) = sup ϱ [ 0 , 1 ] max { | ψ ̲ ϱ ̲ ϱ | , | ψ ¯ ϱ ¯ ϱ | } ,
where the ϱ -level sets of ψ and are [ ψ ] ϱ = ( ψ ̲ ϱ , ψ ¯ ϱ ) and [ ] ϱ = ( ̲ ϱ , ¯ ϱ ) , respectively. It is easy to note that D is a metric in E 1 and ( D , E 1 ) is a complete metric space [61]. The ϱ -level set of fuzzy functions ψ : A R E 1 can be represented by [ ψ ( μ ) ] ϱ = [ ψ ̲ ( μ , ϱ ) , ψ ¯ ( μ , ϱ ) ] , μ A R and 0 ϱ 1 .
Definition 2
([62,63]). For arbitrary fuzzy numbers ψ , E 1 , ψ = [ ψ ̲ ϱ , ψ ¯ ϱ ] , = [ ̲ ϱ , ¯ ϱ ] , the quantity D ( ψ , ) = sup ϱ [ 0 , 1 ] max { | ψ ̲ ϱ ̲ ϱ | , | ψ ¯ ϱ ¯ ϱ | } is the distance between ψ and ˜
  • ( E 1 , D ) denote a complete metric space,
  • D ( ψ w , w ) = D ( ψ , ) , ψ , , w E 1 ,
  • D ( ψ , w e ) D ( ψ , w ) + D ( , e ) , ψ , , w , e E 1 ,
  • D ( ψ , 0 ˜ ) D ( ψ , 0 ˜ ) + D ( , 0 ˜ ) , ψ , E 1 ,
  • D ( ρ ψ , ρ ) = | ρ | D ( ψ , ) , ψ , E 1 , ρ R ,
  • D ( ρ 1 ψ , ρ 2 ψ ) = | ρ 1 ρ 2 | D ( ψ , 0 ˜ ) , ψ E 1 , ρ 1 , ρ 2 R , with ρ 1 · ρ 2 0 .
Let us consider the definition of the Hukuhara difference (H-difference) in [64]. The Hukuhara H-difference is proposed as a set w for which ψ g H ˜ = w ψ = w . The H-difference is unique; however, it does not always exist (a necessary condition for ψ g H to exist is that ψ contains a translation { c } o f ).
Definition 3
([63,64]). The generalized Hukuhara difference between two fuzzy numbers ψ , E 1 is given as
ψ g H = w ( i ) ψ = w , o r ( i i ) = ψ ( w ) .
In terms of the ϱ-level, we obtain [ ψ g H ] = [ min { ψ ϱ ̲ ϱ , ψ ¯ ϱ ¯ ϱ } , max { ψ ̲ ϱ ̲ ϱ , ψ ¯ ϱ ¯ ϱ } ] , and if the H-difference exists, then ψ = ψ g H ; the conditions for the existence of w = ψ g H E 1 are
C a s e ( i ) w ̲ ϱ = ψ ̲ ϱ ̲ ϱ a n d w ¯ ϱ = ψ ¯ ϱ ¯ ϱ , ϱ [ 0 , 1 ] , w i t h w ̲ ϱ i n c r e a s i n g , w ¯ ϱ d e c r e a s i n g , w ̲ ϱ w ¯ ϱ .
C a s e ( i i ) w ̲ ϱ = ψ ¯ ϱ ¯ ϱ a n d w ¯ ϱ = ψ ̲ ϱ ̲ ϱ , ϱ [ 0 , 1 ] , w i t h w ̲ ϱ i n c r e a s i n g , w ¯ ϱ d e c r e a s i n g , w ̲ ϱ w ¯ ϱ .
It is simple to demonstrate that (i) and (ii) are both true provided that w is a crisp number.
Definition 4
([65]). Assume that ψ ( ϑ , θ ) : D E 1 and ( ϑ 0 , θ ) D . We say that ψ is strongly generalized Hukuhara differentiable on ( ϑ 0 , θ ) (GH-differentiable for short) if there exists an element ψ ϑ | ( ϑ 0 , θ ) E 1 such that
(i) 
for all > 0 sufficiently small, ψ ( ϑ 0 + , θ ) H ψ ( ϑ 0 , θ ) , ψ ( ϑ 0 , t ) H ψ ( ϑ 0 , θ ) and the limits (in the metric D)
lim 0 + ψ ( ϑ 0 + , θ ) H ψ ( ϑ 0 , θ ) = lim 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) = ψ ϑ | ( ϑ 0 , θ ) ,
or
(ii) 
for all > 0 sufficiently small, ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 + , θ ) , ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) and the limits
lim 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 + , θ ) = lim 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) = ψ ϑ | ( ϑ 0 , θ ) ,
or
(iii) 
for all > 0 sufficiently small, ψ ( ϑ 0 + , θ ) H ψ ( ϑ 0 , θ ) , ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) and the limits
lim 0 + ψ ( ϑ 0 + , θ ) H ψ ( ϑ 0 , θ ) = lim 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) = ψ x | ( ϑ 0 , θ ) ,
or
(iv) 
for all > 0 sufficiently small, ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 + , θ ) , ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) and the limits
lim 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 + , θ ) = lim h 0 + ψ ( ϑ 0 , θ ) H ψ ( ϑ 0 , θ ) = ψ ϑ | ( ϑ 0 , θ ) .
Definition 5
([64]). Assume that ψ ( ϑ , θ ) : D E 1 is a function and set ψ ( ϑ , θ ) = [ ( ψ ̲ ( ϑ , θ ; ϱ ) , ψ ¯ ( ϑ , θ ; ϱ ) ] for each ϱ [ 0 , 1 ] . Then,
(1) If ψ is gH-differentiable in the first form (i), then ψ ̲ ( ϑ , θ ; ϱ ) and ψ ¯ ( ϑ , θ ; ϱ ) are differentiable functions and
ψ ϑ ϱ = ψ ̲ ( ϑ , θ ; ϱ ) ϑ , ψ ¯ ( ϑ , θ ; ϱ ) μ ,
(2) If ψ is gH-differentiable in the second form (ii), then ψ ¯ ( ϑ , θ ; ϱ ) and ψ ̲ ( ϑ , θ ; ϱ ) are differentiable functions and
ψ ϑ ϱ = ψ ¯ ( ϑ , θ ; ϱ ) ϑ , ψ ̲ ( ϑ , θ ; ϱ ) ϑ .
Definition 6
([66]). A fuzzy-valued function ψ of two variables is a rule that assigns to every ordered pair of real numbers, ( ϑ , θ ) , in a set D , a unique fuzzy number denoted by ψ ( ϑ , θ ) . The set D is the domain of ψ and its range is the set of values that ψ takes on, which is { ψ ( ϑ , θ ) | ( ϑ , θ ) D } . The parametric representation of the fuzzy-valued function ψ : D E 1 is expressed by ψ ( ϑ , θ ; σ ) = [ ψ ̲ ( ϑ , θ ; ϱ ) , ψ ¯ ( ϑ , θ ; ϱ ) ] , for any ( ϑ , θ ) D and ϱ [ 0 , 1 ] .
Definition 7
([66]). A fuzzy-valued function ψ : [ , ] E 1 is said to be continuous at θ 0 [ , ] if, for every positive ε, there is δ > 0 such that D ψ ( θ ) , ψ θ 0 < ε whenever θ θ 0 < δ . If ψ is continuous for each θ [ , ] , then we say that ψ is fuzzy continuous on [ , ] .
Next, we regard C F [ , ] as the space of all continuous fuzzy-valued functions on [ , ] , and we recall the space of all Lebesgue integrable fuzzy functions on the bounded interval [ , ] R by L F [ , ] , refs. [67].
Definition 8
([67]). Assume that ψ ( μ ) C F [ , ] L F [ , ] . The fuzzy Riemann–Liouville integral of fuzzy function ψ is given as
I a + α ψ ( μ ) = 1 Γ ( α ) μ ψ ( ι ) d ι ( μ ι ) 1 α , μ > , 0 < α 1 .
Let the ϱ-level expression of a fuzzy function ψ as f ( μ ; ϱ ) = [ ψ ̲ ( μ ; ϱ ) , ψ ¯ ( μ ; ϱ ) ] , for 0 ϱ 1 .
Definition 9
([60,67]). Let ψ C F [ , ] L F [ , ] be fuzzy function and 0 < α 1 . Then, ψ is said to be Caputo gH-differentiable at μ when
C D α ψ ( μ ; ϱ ) = 1 Γ ( 1 α ) μ 0 μ ( μ ι ) α ψ ( ι ; ϱ ) d ι .
Note that, later, we indicate C D 0 α ψ ( ι ; ϱ ) using C D α ψ ( ι ; ϱ ) .
Theorem 1
([60]). Let ψ C F [ , ] L F [ , ] , μ 0 ( a , b ) and 0 < α 1 . Then,
(i) 
if ψ is an (i)-differentiable fuzzy function, then
i C D μ 0 α ψ ( μ ; ϱ ) = C D ϑ 0 α ψ ̲ ( μ ; ϱ ) , C D μ 0 α ψ ¯ ( μ ; ϱ ) , 0 ϱ 1 ,
(ii) 
if ψ is an (ii)-differentiable fuzzy function, then
i i C D μ 0 α ψ ( μ ; ϱ ) = C D μ 0 α ψ ¯ ( μ ; ϱ ) , C D μ 0 α ψ ̲ ( μ ; ϱ ) , 0 ϱ 1 .

Local Fractional Calculus

We present the local fractional calculus, refs. [24,27,68], as follows.
Definition 10.
The real-valued function W ( ϑ ) is known local fractional continuous at ϑ = ϑ 0 and is displayed by lim ϑ = ϑ 0 W ( ϑ ) = W ( ϑ 0 ) if there exists a relation.
At ϑ = ϑ 0 , the real-valued function W ( ϑ ) is known as local fractional continuous and is shown by lim ϑ = ϑ 0 W ( ϑ ) = W ( ϑ 0 ) if ∃, if there is a relation.
| W ( ϑ ) W ( ϑ 0 ) | < δ β , 0 < β 1 with | ϑ ϑ 0 | < ε , for ε , δ R . In the same way, W ( ϑ ) is referred to as local fractional continuous on ( γ 1 , γ 2 ) and is symbolized as W ( ϑ ) C β ( γ 1 , γ 2 ) , provided that | W ( ϑ ) W ( ϑ 0 ) | < δ β holds for W ( γ 1 , γ 2 ) .
Definition 11.
Assume that ( γ 1 , γ 2 ) , is the interval and Δ θ = max { Δ θ 0 , Δ θ 1 , Δ θ 2 , } is a partition of ( γ 1 , γ 2 ) , with ( θ ε , θ ε + 1 ) , ε = 0 , 1 , 2 , , N 1 , θ 0 = γ 1 , θ N = γ 2 with Δ θ ε = θ ε + 1 θ ε . Then, the local fractional integral of W ( ϑ ) C β ( γ 1 , γ 2 ) is defined as
γ 1 I γ 2 β W ( ϑ ) = 1 Γ ( 1 + β ) γ 1 γ 2 W ( θ ) ( d θ ) β = 1 Γ ( 1 + β ) lim Δ θ 0 ε = 0 N 1 W ( θ ε ) ( Δ θ ε ) β .
Definition 12.
The Mittage–Leffler function in fractal media is represented as
E β ( ϑ β ) = m = 0 ϑ m β Γ ( 1 + m β ) , 0 < β 1 .
Definition 13.
In fractal media, the trigonometric function is stated as
sin β ( ϑ β ) = m = 0 ( 1 ) m ϑ ( 2 m + 1 ) β Γ ( 1 + ( 2 m + 1 ) β ) cos β ( ϑ β ) = m = 0 ( 1 ) m ϑ 2 m β Γ ( 1 + 2 m β ) , ϑ R , 0 < β 1 .
Definition 14.
Assume that W ( ϑ ) satisfies the condition of local fractional continuity according to Definition 10. Hence, we have the inverse formula of the local fractional integral described in Definition 11 as
D ϑ β W ( ϑ 0 ) = d β W ( ϑ 0 ) d ϑ β = W β ( ϑ 0 ) = lim ϑ ϑ 0 Δ β ( W ( ϑ ) W ( ϑ 0 ) ) ( ϑ ϑ 0 ) β , ϑ ( γ 1 , γ 2 ) ,
where Δ β ( W ( ϑ ) W ( ϑ 0 ) ) Γ ( β + 1 ) ( W ( ϑ ) W ( ϑ 0 ) ) , and D ϑ β W ( ϑ 0 ) denote the local fractional derivative of W ( ϑ ) C β ( γ 1 , γ 2 ) of order β at ϑ = ϑ 0 .
The local fractional partial derivative of W ( ϑ , ξ ) C β ( γ 1 , γ 2 ) of order β was provided by Yang [24,68] as follows:
β ξ β W ( ϑ , ξ ) = Δ β ( W ( ϑ , ξ ) W ( ϑ , ξ 0 ) ) ( ξ ξ 0 ) β , ϑ ( γ 1 , γ 2 ) ,
where Δ β ( W ( ϑ , θ ) W ( ϑ , θ 0 ) Γ ( β + 1 ) ( W ( ϑ , θ ) W ( ϑ , θ 0 ) ) .
This study makes use of the formulas for local fractional derivatives and local fractional integrals of a few spatial functionals that are listed in [24,68]: D ϑ β a W ( ϑ ) = a D ϑ β W ( ϑ ) , where a is constant
D ϑ β ϑ m β Γ ( 1 + m β ) = ϑ ( m 1 ) β Γ ( 1 + ( m 1 ) β ) , I ϑ β ϑ m β Γ ( 1 + m β ) = ϑ ( m + 1 ) β Γ ( 1 + ( m + 1 ) β ) , m N ,
where ϑ β specifies the Cantor function.

3. Local Fuzzy Fractional Partial Differential Equations

In this section, the hybrid fuzzy fractional approach is applied to some local fuzzy fractional partial differential equations.

3.1. Local Fractional Sumudu Transform

We introduce the following definitions of the local fractional Sumudu transform (LFST), refs. [59,69].
Definition 15.
The LFST of a function W ( θ ) is given as
L F S β { W ( θ ) } = W ¯ β ( z ) = 1 Γ ( 1 + β ) 0 E β ( z β θ β ) W ( θ ) z β ( d θ ) β , 0 < β 1 .
Applying the formula above, the inverse LFST is defined as
L F S β 1 { W ¯ β ( z ) } = W ( θ ) , 0 < β 1 .
Remark 1.
L F S β θ β Γ ( 1 + β ) = z β .
Theorem 2.
(Linearity). If L F S β { W ( θ ) } = W ¯ β ( z ) and L F S β { V ( θ ) } = V ¯ β ( z ) , we get
L F S β { W ( θ ) + V ( θ ) } = W ¯ β ( z ) + V ¯ β ( z ) , 0 < β 1 .
Theorem 3.
(Local fractional convolution). If L F S β { W ( θ ) } = W ¯ β ( z ) and L F S β { V ( θ ) } = V ¯ β ( z ) , we obtain
L F S β { W ( θ ) V ( θ ) } = z β W ¯ β ( z ) V ¯ β ( z ) with W ( θ ) V ( θ ) = 1 Γ ( 1 + β ) 0 W ( y ) V ( θ y ) ( d θ ) β .
Theorem 4.
If L F S β { W ( θ ) } = W ¯ β ( z ) , then
L F S β d β W ( θ ) d θ β = W ¯ β ( z ) W ( 0 ) z β , L F S β d m β W ( θ ) d θ m β = 1 z m β W ¯ β ( z ) k = 0 m 1 z k β W k β ( 0 ) .
For m = 2 , we obtain
L F S β d 2 β W ( θ ) d θ 2 β = 1 z 2 β W ¯ β ( z ) W ( 0 ) z β W β ( 0 ) ,
where W ¯ β ( z ) = L F S β { W ( θ ) } .
Theorem 5.
If L F S β { W ( θ ) } = W ¯ β ( z ) , we obtain L F S β 0 I θ ( β ) W ( θ ) = z β W ¯ β ( z ) .
Theorem 6.
The LFST of some special functions
L F S β { a } = a , L F S β { E β ( a θ β ) } = 1 1 a z β ,
where a is a constant.
L F S β θ β Γ ( 1 + β ) = z β , L F S β k = 0 a k θ β k = k = 0 Γ ( 1 + k β ) a k z β k , L F S β { cos β ( a θ β ) } = 1 1 + a 2 z 2 β , L F S β { sin β ( a θ β ) } = a z β 1 + a 2 z 2 β , L F S β { sinh β ( a θ β ) } = a z β 1 a 2 z 2 β , L F S β { cosh β ( a θ β ) } = a z β 1 a 2 z 2 β .

3.2. Local Fuzzy Fractional Homotopy Perturbation Sumudu Transform Method

We provide the fundamental plan for the LFFHPSTM. The following LFFPDEs are used to demonstrate the basic steps of the LFFHPSTM.
L β W ( ϑ , θ ) U β W ( ϑ , θ ) R β W ( ϑ , θ ) = h ( ϑ , θ ) , 0 < ϑ < 1 , 0 < θ < 1 ,
where L β = k β θ k β is called the linear local fuzzy fractional differential operator (LFFDO) of order k β , U β is called the LFFDO of order less than L β , R β indicates the nonlinear local differential operator in ϑ and θ are variables, W ( ϑ , θ ) specifies the unknown local fractional continuous function, and h ( ϑ , θ ) is called the nowhere differentiable source term.
Furthermore, the fundamental strategy of LFFHPSTM proposes the use of LFFST in Equation (16)
L F F S β L β W ( ϑ , θ ) L F F S β U β W ( ϑ , θ ) L F F S β R β W ( ϑ , θ ) = L F F S β h ( ϑ , θ ) .
Using the property of LFFST, we have
W β ( ϑ , z ) = W β ( ϑ , 0 ) z β W β ( ϑ , 0 ) z 2 β W 2 β ( ϑ , 0 ) z ( k 1 ) β W ( k 1 ) β ( ϑ , 0 ) z k β L F F S β U β W ( ϑ , θ ) z k β L F F S β R β W ( ϑ , θ ) z k β L F F S β R β h ( ϑ , θ ) ,
where W ¯ β ( ϑ , z ) = L F F S β W ( ϑ , θ ) , or
W ( ϑ , θ ) = m = 0 k 1 z m β m β W ( ϑ , 0 ) θ m β z k β L F F S β h ( ϑ , θ ) z k β L F F S β U β ( ϑ , θ ) z k β L F F S β R β W ( ϑ , θ ) .
Using Equation (19) and the reverse of LFFST, we obtain
W ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β L F F S β h ( ϑ , θ ) L F F S β 1 z k β L F F S β U β ( ϑ , θ ) z k β L F F S β R β W ( ϑ , θ ) .
The basic steps include the following. The expansion of W ( ϑ , θ ) takes the place of a power series of ρ as
W ( ϑ , θ ) = n = 0 ρ n β W n ( ϑ , θ ) ,
and the nonlinear component is presented as
R β W ( ϑ , θ ) = n = 0 ρ n β H n ( W ) ,
where ρ ( ϑ , θ ) is an embedding variable and H n ( W ) is a local fractional He’s polynomial represented as
H n ( W 0 , W 1 , , W n ) = 1 Γ ( n β + 1 ) n β ρ n β R β i = 0 ρ β i W i ρ = 0 , n = 0 , 1 , 2 ,
Applying the values of W ( ϑ , θ ) and R β W ( ϑ , θ ) in Equation (20) yields the following homotopy equation:
m = 0 k 1 ρ n β W n ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β L F F S β h ( ϑ , θ ) ρ β L F F S β 1 z k β L F F S β U β n = 0 ρ n β W n ( ϑ , θ ) z k β L F F S β n = 0 ρ n β H n ( W ) .
Furthermore, comparing the coefficients of the same powers of ρ yields
ρ 0 : W 0 ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β L F F S β h ( ϑ , θ ) ρ 1 β : W 1 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 0 ( ϑ , θ ) ) H 0 ( W ) ρ 2 β : W 2 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 1 ( ϑ , θ ) H 1 ( W ) ρ 3 β : W 3 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 2 ( ϑ , θ ) ) H 2 ( W ) .
Finally, the local fractional series solution of Equation (16) can be written as
W ( ϑ , θ ) = lim N n = 0 N W n ( ϑ , θ ) .

Convergence Analysis

We consider the Banach space C [ 0 , T ] of all continuous fuzzy functions on [ 0 , T ] with the supremum norm. Throughout this section, we consider W ( ϑ , θ ) , W n ( ϑ , θ ) C [ 0 , T ] .
Theorem 7.
The solution obtained by the LFFHPSTM of fuzzy partial differential Equation (16) has a unique solution, whenever 0 < ζ < 1 .
Proof. 
The solution of Equation (16) is of the form W ( ϑ , θ ) = n = 0 ρ n β W n ( ϑ , θ ) . Here,
W ̲ ( ϑ , θ ; ϱ ) = k = 0 n 1 θ k k ! W ̲ ( k ) ( ϑ , 0 ) ( ϱ ) + L F F S β 1 w n β L F F S β h ̲ ( ϑ , θ ; ϱ ) U β W ̲ ( ϑ , θ ; ϱ ) R β W ̲ ( ϑ , θ ; ϱ ) , W ¯ ( ϑ , θ ; ϱ ) = k = 0 n 1 θ k k ! W ¯ ( k ) ( ϑ , 0 ) ( ϱ ) + L F F S β 1 w n β L F F S β h ¯ ( ϑ , θ ; ϱ ) U β W ¯ ( ϑ , θ ; ϱ ) R β W ¯ ( ϑ , θ ; ϱ ) ,
where ϱ [ 0 , 1 ] .
Assume that W ( ϑ , θ ) and Y ( ϑ , θ ) are the distinct solutions of Equation (16); then,
| W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | = | L F F S β 1 [ w n β L F F S β [ U ( W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) + R W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) ) ] ] | ,
| W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | = | L F F S β 1 [ w n β L F F S β [ U ( W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) + R W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) ) ] ] | .
Applying the convolution theorem,
| W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | 0 θ ( | [ U ( W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) ) | + | R W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) ] | ) ( θ η ) n n ! d η 0 θ ( μ | [ ( W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) ) | + δ | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) ] | ) ( θ η ) n n ! d η ,
and
| W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | 0 θ ( | [ U ( W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) ) | + | R W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) ] | ) ( θ η ) n n ! d η 0 θ ( μ | [ ( W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) ) | + δ | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) ] | ) ( θ η ) n n ! d η .
U is a bounded operator,
| U W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | μ | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | | U W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | μ | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) |
and R satisfies the Lipschitz condition with δ > 0 such that
| R W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | δ | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | | R W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | δ | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) |
| W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | 0 θ ( μ + δ ) | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | ( θ η ) n n ! d η | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | 0 θ ( μ + δ ) | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | ( θ η ) n n ! d η .
Using the value theorem of integral calculus,
| W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | ( μ + δ ) | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | MT , | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | ( μ + δ ) | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | MT ,
where M = max ( θ η ) n and θ [ 0 , T ] .
Hence,
| W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | ζ , | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | ζ ,
where ζ = ( n + δ ) MT . Then,
( 1 ζ ) | W ̲ ( ϑ , θ ; ϱ ) Y ̲ ( ϑ , θ ; ϱ ) | 0 , ( 1 ζ ) | W ¯ ( ϑ , θ ; ϱ ) Y ¯ ( ϑ , θ ; ϱ ) | 0 ,
implies W ( ϑ , θ ) = Y ( ϑ , θ ) whenever 0 < ζ < 1 .
Theorem 8.
Assume n = 0 W n W if ζ ( 0 , 1 ) such that W n ζ W n 1 n N .
Proof. 
For the convergence of sequence { S n ( ϑ , θ ) } of the partial sums of the series (26), we prove that { S n ( ϑ , θ ) } is a Cauchy sequence in ( C [ 0 , T ] , ) . As
S ̲ n + 1 ( ϑ , θ ; ϱ ) S ̲ n ( ϑ , θ ; ϱ ) = W ̲ n + 1 ( ϑ , θ ; ϱ ) ζ W ̲ n ( ϑ , θ ; ϱ ) ζ 2 W ̲ n 1 ( ϑ , θ ; ϱ ) ζ n + 1 W ̲ 0 ( ϑ , θ ; ϱ )
and
S ¯ n + 1 ( ϑ , θ ; ϱ ) S ¯ n ( ϑ , θ ; ϱ ) = W ¯ n + 1 ( ϑ , θ ; ϱ ) ζ W ¯ n ( ϑ , θ ; ϱ ) ζ 2 W ¯ n 1 ( ϑ , θ ; ϱ ) ζ n + 1 W ¯ 0 ( ϑ , θ ; ϱ )
Hence,
S ̲ n ( ϑ , θ ; ϱ ) S ̲ m ( ϑ , θ ; ϱ ) = i = m + 1 n W ̲ i ( ϑ , θ ; ϱ ) i = m + 1 n W ̲ i ( ϑ , θ ; ϱ ) ζ m + 1 i = 0 n m 1 ζ i W ̲ 0 ( ϑ , θ ; ϱ ) = ζ m + 1 ( 1 ζ n m ) 1 ζ W ̲ 0 ( ϑ , θ ; ϱ ) , n , m N ,
and
S ¯ n ( ϑ , θ ; ϱ ) S ¯ m ( ϑ , θ ; ϱ ) = i = m + 1 n W ¯ i ( ϑ , θ ; ϱ ) i = m + 1 n W ¯ i ( ϑ , θ ; ϱ ) ζ m + 1 i = 0 n m 1 ζ i W ¯ 0 ( ϑ , θ ; ϱ ) = ζ m + 1 ( 1 ζ n m ) 1 ζ W ¯ 0 ( ϑ , θ ; ϱ ) , n , m N .
since 0 < ζ < 1 . Hence,
S ̲ n ( ϑ , θ ; ϱ ) S ̲ m ( ϑ , θ ; ϱ ) ζ m + 1 1 ζ W ̲ 0 ( ϑ , θ ; ϱ ) , S ¯ n ( ϑ , θ ; ϱ ) S ¯ m ( ϑ , θ ; ϱ ) ζ m + 1 1 ζ W ¯ 0 ( ϑ , θ ; ϱ ) .
Moreover, W 0 ( ϑ , θ ) is bounded; therefore, S n + 1 ( ϑ , θ ) S n ( ϑ , θ ) 0 as m , n . Thus, { S n ( ϑ , θ ) } is a Cauchy sequence in C [ 0 , T ] , and  n = 0 W n ( ϑ , θ ) is convergent. □

3.3. Local Fuzzy Fractional Sumudu Decomposition Method

In this section, we present the fuzzy linear operator with a local fractional derivative as follows,
L β U ( ϑ , θ ) R β U ( ϑ , θ ) = g ( ϑ , θ )
where L β = m β t m β ( m N * ) represents the linear local fuzzy fractional derivative operator of order m β , R β represents the linear local fuzzy fractional derivative operator of order less than L β , and g ( ϑ , θ ) is a non-differentiable source term. Taking the local fractional Sumudu transform (denoted by L F F S β ) on both sides of Equation (33), we obtain
L F F S β L β U ( ϑ , θ ) L F F S β R β U ( ϑ , θ ) = L F F S β g ( ϑ , θ ) .
Applying the prosperity of the local fractional Sumudu transform,
U ( ϑ , θ ) = k = 0 m 1 u k v k v U ( r , 0 ) θ k v θ k v Γ ( 1 + k v ) + L F F S β 1 u m v ( L F F S [ g ( ϑ , θ ) ] ) L F F S β 1 u m v ( L F F S [ R β U ( ϑ , θ ) ] ) .
The unknown function U is then divided into infinite series by
W ( ϑ , θ ) = n = 0 W n ( ϑ , θ ) .
Substituting (36) into (35), we obtain
n = 0 W n ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β ( L F F S β g ( ϑ , θ ) ) L F F S β 1 z k β L F F S β R β n = 0 W n ( ϑ , θ ) .
On comparing Equation (37), we obtain
W 0 ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β ( L F F S β g ( ϑ , θ ) ) , W 1 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 0 ( ϑ , θ ) ) , W 2 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 1 ( ϑ , θ ) , W 3 ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W 2 ( ϑ , θ ) ) , .
The general form of the local fractional recursive relation is
W 0 ( ϑ , θ ) = m = 0 k 1 θ m β Γ ( 1 + m β ) m β W ( ϑ , 0 ) θ m β L F F S β 1 z k β ( L F F S β g ( ϑ , θ ) ) , W n ( ϑ , θ ) = L F F S β 1 z k β L F F S β U β ( W n 1 ( ϑ , θ ) ) ,
where 0 < v 1 , n N * and m = 1 , 2 , 3 ,

3.4. Local Fuzzy Fractional Sumudu Variational Iteration Method

We consider the following nonlinear fuzzy operator with a local fractional derivative as
L β W ( ϑ , θ ) U β W ( ϑ , θ ) R β W ( ϑ , θ ) = h ( ϑ , θ ) ,
where L β = 2 β θ 2 β specifies the linear local fuzzy fractional differential operator of order 2 β , U β is called the linear local fuzzy fractional derivative operator of order less than in L β , N β is represent the nonlinear local fractional operator, and h ( ϑ , θ ) is a non-differentiable source term. Using the LFFST of Equation (40), we obtain
L F F S β L β W ( ϑ , θ ) L F F S β U β W ( ϑ , θ ) L F F S β R β W ( ϑ , θ ) = L F F S β h ( ϑ , θ ) .
Applying the property of the LFFST, we obtain
L F F S β W ( ϑ , θ ) = W ( ϑ , 0 ) W θ ( β ) ( ϑ , 0 ) u β u 2 β L F F S β h ( ϑ , θ ) u 2 β L F F S β N β W ( ϑ , θ ) R β W ( ϑ , θ ) .
Using the inverse LFFST of Equation (42) yields
W ( ϑ , θ ) = W ( ϑ , 0 ) W θ ( β ) ( ϑ , 0 ) θ β Γ ( 1 + β ) L F F S β 1 u 2 β L F F S β W ( ϑ , θ ) N β W ( ϑ , θ ) R β W ( ϑ , θ ) .
Taking β θ β of Equation (43), we obtain
β W ( ϑ , θ ) θ β W θ ( β ) ( ϑ , 0 ) L F F S β 1 u 2 β L F F S β N β W n ( ϑ , θ ) R β W n ( ϑ , θ ) h ( ϑ , θ ) = 0 .
The following limit is used to determine the solution
W ( ϑ , θ ) = lim ϑ W n ( ϑ , θ ) .

3.5. Examples

In this part, we present the application of LFFHPSTM for local fuzzy fractional partial differential equations.
Example 1.
We consider the following local fuzzy fractional differential equation on Cantor:
β W ( ϑ , θ ) θ β 2 β W ( ϑ , θ ) ϑ 2 β W ( ϑ , θ ) = 0 , θ > 0 , ϑ R , 0 < β 1 ,
with the initial condition
W ( ϑ , 0 ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) ,
where n = 1 , 2 , 3 ,
1. Local fuzzy fractional homotopy perturbation Sumudu transform method
Using Equation (46) and the LFFST operator LFFS, we obtain
L F F S β β W ( ϑ , θ ) θ β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β = L F F S β ( W ( ϑ , θ ) ) .
Applying the LFFST for the LFFD formula to Equation (48) yields
L F F S β 2 β W ( ϑ , θ ) ϑ 2 β 1 z β L F F S β ( W ( ϑ , θ ) ) W ( ϑ , 0 ) = L F F S β ( W ( ϑ , θ ) ) .
Rearranging the terms in Equation (48) yields
L F F S β ( W ( ϑ , θ ) ) = W ( ϑ , 0 ) z β L F F S β ( W ( ϑ , θ ) ) z β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
Further simplification on account of Equation (47) reduces Equation (50) to
L F F S β ( W ( ϑ , θ ) ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) z β L F F S β ( W ( ϑ , θ ) ) z β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
Using the inverse LFFST of Equation (51), we obtain
W ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) L F F S β 1 z β L F F S β ( W ( ϑ , θ ) ) L F F S β 1 z β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
The LFHPM suggests the following homotopy equation:
n = 0 ρ n β W n ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) ρ β L F F S β 1 z β L F F S β n = 0 ρ n β W n ( ϑ , θ ) ρ β L F F S β 1 z β L F F S β 2 β n = 0 ρ n β W n ( ϑ , θ ) ϑ 2 β .
Comparing the like powers of ρ , we obtain
ρ 0 : W 0 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) ρ 1 β : W 1 ( ϑ , θ ) = L F F S β 1 z β L F F S β W 0 ( ϑ , θ ) L F F S β 1 z β L F F S β 2 β W 0 ( ϑ , θ ) ϑ 2 β ρ 2 β : W 2 ( ϑ , θ ) = L F F S β 1 z β L F F S β W 1 ( ϑ , θ ) L F F S β 1 z β L F F S β 2 β W 1 ( ϑ , θ ) ϑ 2 β ρ 3 β : W 3 ( ϑ , θ ) = L F F S β 1 z β L F F S β W 2 ( ϑ , θ ) L F F S β 1 z β L F F S β 2 β W 2 ( ϑ , θ ) ϑ 2 β ρ n β : W n ( ϑ , θ ) = L F F S β 1 z β L F F S β W n 1 ( ϑ , θ ) L F F S β 1 z β L F F S β 2 β W n 1 ( ϑ , θ ) ϑ 2 β .
After simplification, we get
W 0 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) W 1 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 2 θ β Γ ( 1 + β ) W 2 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 4 θ 2 β Γ ( 1 + 2 β ) W 3 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 8 θ 3 β Γ ( 1 + 3 β ) W n ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 2 θ β n Γ ( 1 + n β ) .
Consequently, the solution of Equation (46) is as follows:
W ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 1 + 2 θ β Γ ( 1 + β ) + 4 θ 2 β Γ ( 1 + 2 β ) + 8 θ 3 β Γ ( 1 + 3 β ) + + 2 θ β n Γ ( 1 + n β ) + = sin β ( ϑ β ) n = 0 ( 2 θ β ) n Γ ( 1 + n β ) ,
or
W ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) E β ( 2 θ β ) .
2. Local fuzzy fractional Sumudu decomposition method
The subsequent approximations from (39) and (46) are
U 0 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) U n ( ϑ , θ ) = L F F S β 1 u β L F F S β 2 β U n 1 ( ϑ , θ ) ϑ 2 β U n 1 ( ϑ , θ ) , n 1 .
The subsequent formula (57) yields
U 0 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) U 1 ( ϑ , θ ) = L F F S β 1 u β L F F S β 2 β U 0 ( ϑ , θ ) ϑ 2 β U 0 ( ϑ , θ ) , U 2 ( ϑ , θ ) = L F F S β 1 u β L F F S β 2 β U 1 ( ϑ , θ ) ϑ 2 β U 1 ( ϑ , θ ) , U 3 ( ϑ , θ ) = L F F S β 1 u β L F F S β 2 β U 2 ( ϑ , θ ) ϑ 2 β U 2 ( ϑ , θ ) ,
The first terms of LFFSDM are given by the equations (58) as
U 0 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) U 1 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 2 θ v Γ ( 1 + 2 v ) , U 2 ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 4 θ 2 Γ ( 1 + 2 v ) ,
The local fuzzy fractional series form is
W ( ϑ , θ ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) 1 + 2 θ β Γ ( 1 + β ) + 4 θ 2 β Γ ( 1 + 2 β ) + 8 θ 3 β Γ ( 1 + 3 β ) + + 2 θ β n Γ ( 1 + n β ) + .
Thus, we can obtain the exact solution as follows:
W ( ϑ , θ ) = sin β ( ϑ β ) n = 0 ( 2 θ β ) n Γ ( 1 + n β ) = ( 2.5 + 0.5 ϱ ) n , ( 3.5 0.5 ϱ ) n sin β ( ϑ β ) E β ( 2 θ β ) .
In Figure 1, we plot 2 D graphs of the exact and approximation solutions for local fuzzy fractional differential equations (LFFDEs) on Cantor set. Figure 1a shows that for θ = 45 , ϑ = 2 , β = 1 , n = 2 , the LFFDEs become bounded and closed. In addition, the W ̲ sign shows increasing functions and W ¯ denotes decreasing functions on the ϱ -level set of W . To reflect the concept of the ϱ -level set, Figure 1a illustrates that the ϱ -level set of LFFDEs is bounded and closed for θ = 45 . Similarly, in Figure 1b, we can observe the same explanation of ϱ -level set closedness and boundedness for Example 1.
Example 2.
We consider the following local fuzzy fractional partial differential equation with a local fuzzy fractional differential operator:
4 β W ( ϑ , θ ) θ 4 β 2 β 3 β W ( ϑ , θ ) ϑ 3 β = 1 , θ > 0 , ϑ R , 0 < β 1 ,
with the initial conditions
W ( ϑ , 0 ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β , β W ( ϑ , 0 ) θ β = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] , 2 β W ( ϑ , 0 ) θ 2 β = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] 2 2 β E β ( 2 ϑ ) β , 3 β W ( ϑ , 0 ) θ 3 β = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] ,
where n = 1 , 2 , 3 ,
1. Local fuzzy fractional homotopy perturbation Sumudu transform method.
Using the LFFST for the LFFD formula in Equation (61), we obtain
L F F S β 4 β W ( ϑ , θ ) θ 4 β L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β = 1 .
Applying the LFFST for the LFD formula in Equation (63), we obtain
1 z 4 β L F F S β W ( ϑ , θ ) W ( ϑ , 0 ) z β W ( β ) ( ϑ , 0 ) z 2 β W 2 β ( ϑ , 0 ) z ( 3 β ) W ( 3 β ) ( ϑ , 0 ) L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β = 1 .
After the rearrangement of the terms in Equation (64), we obtain
L F F S β W ( ϑ , θ ) = W ( ϑ , 0 ) W ( ϑ , 0 ) z β W ( β ) ( ϑ , 0 ) z 2 β W 2 β ( ϑ , 0 ) z ( 3 β ) W ( 3 β ) ( ϑ , 0 ) z 4 β L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β z 4 β .
The initial conditions (62) can transform Equation (68) as
L F F S β W ( ϑ , θ ) = E β ( 2 ϑ ) β z 2 β 2 2 β E β ( 2 ϑ ) β z 4 β L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β z 4 β .
Taking the inverse LFFST on Equation (66), we obtain
W ( ϑ , θ ) = L F F S β 1 E β ( 2 ϑ ) β L F F S β 1 z 2 β 2 2 β E β ( 2 ϑ ) β L F F S β 1 z 4 β L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β L F F S β 1 z 4 β .
After simplification, Equation (67) yields
W ( ϑ , θ ) = E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) L F F S β 1 z 4 β L F F S β 2 β 3 β W ( ϑ , θ ) ϑ 3 β .
The LFFHPM proposes homotopy expression creation as
n = 0 ρ n β W n ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) ρ β L F F S β 1 z 4 β 2 β L F F S β 3 β n = 0 ρ n β W n ( ϑ , θ ) ϑ 3 β .
The following components for the series solution are calculated by comparing the powers of ρ as
ρ 0 : W 0 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) , ρ 1 β : W 1 ( ϑ , θ ) = L F F S β 1 z 4 β 2 β L F F S β 3 β ( W 0 ( ϑ , θ ) ) ϑ 3 β , ρ 2 β : W 2 ( ϑ , θ ) = L F F S β 1 z 4 β 2 β L F F S β 3 β ( W 1 ( ϑ , θ ) ) ϑ 3 β , ρ 3 β : W 3 ( ϑ , θ ) = L F F S β 1 z 4 β 2 β L F F S β 3 β ( W 2 ( ϑ , θ ) ) ϑ 3 β , .
After simplification, we obtain
W 0 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) , W 1 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 4 β θ 4 β Γ ( 1 + 4 β ) 2 6 β θ 6 β Γ ( 1 + 6 β ) W 2 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 8 β θ 4 β Γ ( 1 + 8 β ) 2 10 β θ 10 β Γ ( 1 + 10 β ) W 3 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 12 β θ 12 β Γ ( 1 + 12 β ) 2 14 β θ 14 β Γ ( 1 + 14 β ) .
Hence, we can obtain the solutions of Equation (61) as follows:
W ( ϑ , θ ) = n = 0 W n ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) 2 4 β θ 4 β Γ ( 1 + 4 β ) 2 6 β θ 6 β Γ ( 1 + 6 β ) 2 8 β θ 8 β Γ ( 1 + 8 β ) 2 10 β θ 10 β Γ ( 1 + 10 β ) 2 12 β θ 12 β Γ ( 1 + 12 β ) 2 14 β θ 14 β Γ ( 1 + 14 β ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β 1 ( 2 θ ) 2 β Γ ( 1 + 2 β ) ( 2 θ ) 4 β Γ ( 1 + 4 β ) ( 2 θ ) 6 β Γ ( 1 + 6 β ) ( 2 θ ) 8 β Γ ( 1 + 8 β ) ( 2 θ ) 10 β Γ ( 1 + 10 β ) ( 2 θ ) 12 β Γ ( 1 + 12 β ) ( 2 θ ) 14 β Γ ( 1 + 14 β ) ( 1 ) n ( 2 θ ) 2 n β Γ ( 1 + 2 n β ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] ( θ ) 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β n = 0 ( 1 ) n ( 2 θ ) 2 n β Γ ( 1 + 2 n β )
or
W ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β cos ( 2 ϑ ) β .
2. Local fuzzy fractional Sumudu decomposition method.
Utilizing equations (39) and (61), we can obtain the subsequent mathematical expression
U 0 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 2 β E β ( ( 2 ϑ ) β ) θ 2 β Γ ( 1 + 2 ϑ ) θ 4 β Γ ( 1 + 4 β ) U n ( ϑ , θ ) = L F F S v 1 u 4 β L F F S β 2 β 3 β U n 1 ( ϑ , θ ) ϑ 3 β , n 1 .
Applying the formula (74), successive approximations are obtained as follows:
U 0 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) , U 1 ( ϑ , θ ) = L F F S β 1 u 4 β L F F S β 2 β 3 β U 0 ( ϑ , θ ) ϑ 3 β , U 2 ( ϑ , θ ) = L F F S β 1 u 4 β L F F S β 2 β 3 β U 1 ( ϑ , θ ) ϑ 3 β , U 3 ( ϑ , θ ) = L F F S β 1 u 4 β L F F S β 2 β 3 β U 2 ( ϑ , θ ) ϑ 3 β , .
The first terms of LFFSDM have the following form:
W 0 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 1 2 2 β θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) W 1 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 4 β θ 4 β Γ ( 1 + 2 β ) 2 6 β θ 6 β Γ ( 1 + 6 β ) W 2 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 8 β θ 8 β Γ ( 1 + 8 β ) 2 10 β θ 10 β Γ ( 1 + 10 β ) W 3 ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] E β ( 2 ϑ ) β 2 12 β θ 12 β Γ ( 1 + 12 β ) 2 14 β θ 14 β Γ ( 1 + 14 β ) .
The local fractional series form can be expressed as
W ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) + E β ( 2 ϑ ) β 1 ( 2 θ ) 2 β Γ ( 1 + 2 β ) ( 2 θ ) 4 β Γ ( 1 + 4 β ) ( 2 θ ) 6 β Γ ( 1 + 6 β ) ( 2 θ ) 8 β Γ ( 1 + 8 β ) ( 2 θ ) 10 β Γ ( 1 + 10 β ) ( 1 ) n ( 2 θ ) 2 n β Γ ( 1 + 2 n β ) .
Therefore, the exact solution of (61), as determined by LFFSDM, is given by
W ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β n = 0 ( 1 ) n ( 2 θ ) 2 n β Γ ( 1 + 2 n β ) .
Therefore, the exact solution can be achieved as
W ( ϑ , θ ) = [ ( 3.5 + 0.5 ϱ ) n , ( 4.5 0.5 ϱ ) n ] θ 4 β Γ ( 1 + 4 β ) E β ( 2 ϑ ) β cos ( 2 θ ) β .
Table 1 shows the error term between exact and approximate solutions of Example 2 for ϱ [ 0 , 1 ] .
In Figure 2, we plot 2 D graphs of the exact and approximation solutions for local fuzzy fractional partial differential equations with local fuzzy fractional differential operators. Figure 2a shows that for θ = 30 , ϑ = 3 , β = 1 , n = 1 , the LFFPDEs become bounded and closed. In addition, the W ̲ sign shows increasing functions and W ¯ denotes decreasing functions on the ϱ -level set of W . To reflect the concept of the ϱ -level set, Figure 2a illustrates that the ϱ -level set of LFFPDEs is bounded and closed for θ = 30 . Similarly, in Figure 2b, we can observe the same explanation of ϱ -level set closedness and boundedness for Example 2.
Example 3.
We consider the following local fuzzy fractional partial differential equation
2 β W ( ϑ , θ ) θ 2 β 2 β W ( ϑ , θ ) ϑ 2 β 3 W ( ϑ , θ ) = 0 , θ > 0 , ϑ R , 0 < β 1 ,
with the initial condition
W ( ϑ , 0 ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) , β W ( ϑ , 0 ) θ β = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n 2 sin β ( ϑ β )
where n = 1 , 2 , 3 ,
1. Local fuzzy fractional homotopy perturbation Sumudu transform method
Applying the LFFST of Equation (79), we obtain the following:
L F F S β 2 β W ( ϑ , θ ) θ 2 β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β = 3 L F F S β W ( ϑ , θ ) .
Taking the formula of LFFST for the LFFDs in Equation (81), we obtain
L F F S β 2 β W ( ϑ , θ ) ϑ 2 β 1 z 2 β L F F S β W ( ϑ , θ ) W ( ϑ , 0 ) z β W ( β ) ( ϑ , 0 ) = 3 L F F S β W ( ϑ , θ ) .
Rearranging the terms given in Equation (82) yields
L F F S β W ( ϑ , θ ) = W ( ϑ , 0 ) z β W ( β ) ( ϑ , 0 ) 3 z 2 β L F F S β W ( ϑ , θ ) z 2 β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
Moreover, Equation (83) is simplified for the initial condition (80) in the following manner:
L F F S β W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 2 z β sin β ( ϑ β ) 3 z 2 β L F F S β W ( ϑ , θ ) z 2 β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
We apply the inverse of LFFST on Equation (83) to obtain
W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 2 L F F S β 1 z β sin β ( ϑ β ) 3 L F F S β 1 z 2 β L F F S β W ( ϑ , θ ) L F F S β 1 z 2 β L F F S β 2 β W ( ϑ , θ ) ϑ 2 β .
The homotopy expression can be established using the fundamental method of LFFHPM:
n = 0 ρ n β W n ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 2 sin β ( ϑ β ) θ β Γ ( 1 + β ) 3 ρ β L F F S β 1 z 2 β L F F S β n = 0 ρ n β W n ( ϑ , θ ) ρ β L F F S β 1 z 2 β L F F S β 2 β n = 0 ρ n β W n ( ϑ , θ ) ϑ 2 β .
The comparison of the similar powers of ρ yields the subsequent components of the series solution:
ρ 0 : W 0 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 2 sin β ( ϑ β ) θ β Γ ( 1 + β ) ρ 1 β : W 1 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n 3 L F F S β 1 z 2 β L F F S β W 0 ( ϑ , θ ) L F F S β 1 z 2 β L F F S β 2 β W 0 ( ϑ , θ ) ϑ 2 β ρ 2 β : W 2 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n 3 L F F S β 1 z 2 β L F F S β W 1 ( ϑ , θ ) L F F S β 1 z 2 β L F F S β 2 β W 1 ( ϑ , θ ) ϑ 2 β ρ n β : W n ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n 3 L F F S β 1 z 2 β L F F S β W n 1 ( ϑ , θ ) L F F S β 1 z 2 β L F F S β 2 β W n 1 ( ϑ , θ ) ϑ 2 β .
Following the reduction technique, we have
W 0 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 2 sin β ( ϑ β ) θ β Γ ( 1 + β ) W 1 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 4 θ 2 β Γ ( 1 + 2 β ) 8 θ 3 β Γ ( 1 + 3 β ) W 2 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 16 θ 4 β Γ ( 1 + 4 β ) 32 θ 5 β Γ ( 1 + 5 β ) W 3 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 64 θ 6 β Γ ( 1 + 6 β ) 128 θ 7 β Γ ( 1 + 7 β ) W n ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) ( 2 θ β ) 2 n Γ ( 1 + 2 n β ) ( 2 θ β ) ( 2 n + 1 ) Γ ( 1 + ( 2 n + 1 ) β ) .
Therefore, the solution to Equation (79) is expressed as
W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 1 2 θ β Γ ( 1 + β ) 4 θ 2 β Γ ( 1 + 2 β ) 8 θ 3 β Γ ( 1 + 3 β ) 16 θ 4 β Γ ( 1 + 4 β ) 32 θ 5 β Γ ( 1 + 5 β ) ( 2 θ β ) n Γ ( 1 + n β ) .
Thus, we achieve our desired solution as
W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) E β 2 θ β .
2. Local fuzzy fractional Sumudu variational iteration method
The formula for consecutive approximations given in (45) and (79) is
W n + 1 ( ϑ , θ ) = W n ( ϑ , θ ) 0 I t ( β ) 2 β W n ( ϑ , θ ) θ 2 β β W ( ϑ , 0 ) θ β S β 1 u 2 β S β 2 β W n ( ϑ , θ ) ϑ 2 β 3 W n ( ϑ , θ ) .
The initial conditions (80), the successive formula (89), and the result are as follows:
W 0 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) + 2 sin β ( ϑ β ) θ β Γ ( 1 + β ) W 1 ( ϑ , θ ) = W 0 ( ϑ , θ ) 0 I θ ( β ) 2 β W 0 ( ϑ , θ ) θ 2 β β W ( ϑ , 0 ) θ β S β 1 u 2 β S β 2 β W 0 ( ϑ , θ ) ϑ 2 β 3 W 0 ( ϑ , θ ) W 2 ( ϑ , θ ) = W 1 ( ϑ , θ ) 0 I θ ( β ) 2 β W 1 ( ϑ , θ ) θ 2 β β W ( ϑ , 0 ) θ β S β 1 u 2 β S β 2 β W 1 ( ϑ , θ ) ϑ 2 β 3 W 1 ( ϑ , θ ) .
According to the above equations, the first of LFFSVIM can be obtained as
W 0 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 1 2 θ β Γ ( 1 + β ) , W 1 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 1 2 θ β Γ ( 1 + β ) 4 θ β Γ ( 1 + 2 β ) 8 θ β Γ ( 1 + 3 β ) , W 2 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 1 2 θ β Γ ( 1 + β ) 4 θ β Γ ( 1 + 2 β ) 8 θ β Γ ( 1 + 3 β ) 16 θ 4 β Γ ( 1 + 4 β ) 32 θ 5 β Γ ( 1 + 5 β ) W 3 ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) 1 2 θ β Γ ( 1 + β ) 4 θ β Γ ( 1 + 2 β ) 8 θ β Γ ( 1 + 3 β ) 16 θ 4 β Γ ( 1 + 4 β ) 32 θ 5 β Γ ( 1 + 5 β ) 64 θ 6 β Γ ( 1 + 6 β ) 128 θ 7 β Γ ( 1 + 7 β ) .
Then, the non-differentiable solution of (38) has the form
W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) lim n 1 2 θ β Γ ( 1 + β ) 4 θ β Γ ( 1 + 2 β ) 8 θ β Γ ( 1 + 3 β ) 16 θ 4 β Γ ( 1 + 4 β ) + 32 θ 5 β Γ ( 1 + 5 β ) 64 θ 6 β Γ ( 1 + 6 β ) 128 θ 7 β Γ ( 1 + 7 β ) ( 2 θ β ) n Γ ( 1 + n β ) .
Finally, the exact solution can be obtained as
W ( ϑ , θ ) = ( 5.5 + 0.5 ϱ ) n , ( 6.5 0.5 ϱ ) n sin β ( ϑ β ) E β 2 θ β .

4. Fuzzy Diffusion and Wave Equations on Cantor Sets within Local Fractional Operators

In this section, we present the fuzzy fractional diffusion and wave equations on Cantor sets with local fractional operators by using some techniques.

4.1. Local Fuzzy Fractional Laplace Variational Iteration Method

We consider the following local fuzzy fractional partial differential equations,
L β ψ ( ϑ , θ ) R β ψ ( ϑ , θ ) = f ( ϑ , θ ) ,
where the linear local fractional operator L β and the linear local fractional operator R β , which has an order less than L β , are defined. Additionally, f ( ϑ , θ ) is a non-differentiable source term.
The correction functional for (92) is constructed as follows, based on the rule of the local fractional variational iteration technique
ψ n + 1 ( ϑ ) = ψ n ( ϑ ) 0 I ϑ ( β ) λ ( ϑ θ ) β Γ ( 1 + β ) L β ψ n ( θ ) R β ψ n ( θ ) f ( θ ) ,
where λ ( ϑ θ ) β Γ ( 1 + β ) is a fractal Lagrange multiplier and L β in (92) are k β times LFFPDEs. The initial value problem of (92) can be approached by beginning with
ψ 0 ( ϑ ) = ψ ( 0 ) ϑ β Γ ( 1 + β ) u ( β ) ( 0 ) ϑ ( k 1 ) β Γ [ 1 + ( k 1 ) β ] ψ ( ( k 1 ) β ) ( 0 ) .
Applying the Yang–Laplace transform on the value of (93), yields, in particular
L β ψ n + 1 ( ϑ ) = L β ψ n ( ϑ ) L β 0 I ϑ ( β ) λ ( ϑ θ ) β Γ ( 1 + β ) L β ψ n ( θ ) R β ψ n ( θ ) f ( θ ) .
Assume the local fuzzy fractional variation of Equation (95), as expressed by
δ β L β ψ n + 1 ( ϑ ) = δ β L β ψ n ( ϑ ) δ β L β λ ( ϑ ) β Γ ( 1 + β ) L β L β ψ n ( ϑ ) R β ψ ˜ n ( ϑ ) f ( ϑ ) .
We apply the computation of (96) to obtain
δ β Ł β ψ n + 1 ( ϑ ) = δ β Ł β ψ n ( ϑ ) Ł β λ ( ϑ ) β Γ ( 1 + β ) δ β Ł β L β ψ n ( ϑ ) .
Thus, from (96), we get
1 + Ł β λ ( ϑ ) β Γ ( 1 + β ) s k β = 0 ,
where
δ β Ł β L β ψ n ( ϑ ) = δ β s k β Ł β ψ n ( ϑ ) s ( k 1 ) β ψ n ( 0 ) ψ n ( ( k 1 ) β ) ( 0 ) = s k β δ β Ł β ψ n ( ϑ ) .
Consequently, we obtain
Ł β λ ( ϑ ) β Γ ( 1 + β ) s k β = 1 s k β .
Using the Yang–Laplace transform’s inverse, we have
λ ( ϑ ) β Γ ( 1 + β ) = Ł β 1 1 s k β = ( ϑ ) ( k 1 ) β Γ [ 1 + ( k 1 ) β ] , k N ,
and, as a result, we get
Ł β ψ n + 1 ( ϑ ) = Ł β ψ n ( ϑ ) Ł β 0 I ϑ ( β ) ( ϑ θ ) ( k 1 ) β Γ [ 1 + ( k 1 ) β ] L β ψ n ( θ ) R β ψ n ( θ ) f ( ϑ ) .
Consequently, we have the subsequent iteration algorithm,
Ł β ψ n + 1 ( ϑ ) = Ł β ψ n ( ϑ ) Ł β ( ϑ ) ( k 1 ) β Γ [ 1 + ( k 1 ) β ] Ł β L β ψ n ( ϑ ) R β ψ n ( ϑ ) f ( ϑ ) ,
where the first value can be written as
ψ 0 ( ϑ ) = Ł β 1 s ( k 1 ) β ψ ( 0 ) s ( k 2 ) β ψ ( β ) ( 0 ) ψ n ( ( k 1 ) β ) ( 0 ) s k β = ψ ( 0 ) ϑ β Γ ( 1 + β ) ψ ( β ) ( 0 ) ϑ 2 β Γ ( 1 + 2 β ) ψ ( 2 β ) ( 0 ) ϑ ( k 1 ) β Γ [ 1 + ( k 1 ) β ] ψ ( ( k 1 ) β ) ( 0 ) .
Hence, the solution of the above Equation (92) in terms of a local fuzzy fractional series can be expressed as follows:
ψ ( ϑ , θ ) = lim n Ł β 1 Ł β ψ n ( ϑ , θ ) .

4.2. Local Fuzzy Fractional Series Expansion Method

We consider the following local fractional differential equation
ψ θ n β ( ϑ , θ ) = L β ψ ( ϑ , θ ) ,
where L stands for a linear local operator with regard to ϑ , n { 1 , 2 } . According to the findings of [68,70], there are multi-term separated functions of independent variables θ and ϑ
ψ ( ϑ , θ ) = i = 0 T i ( θ ) X i ( ϑ ) ,
where T i ( θ ) and X i ( ϑ ) are local fuzzy fractional continuous functions. Moreover, there exists a non-differential series term
T i ( θ ) = P i θ i β Γ ( 1 + i β ) ,
where P i represents a coefficient. We can represent the solution in the form presented below:
ψ ( ϑ , θ ) = i = 0 P i θ i β Γ ( 1 + i β ) X i ( ϑ ) .
Then, following (109), we obtain
ψ ( ϑ , θ ) = i = 0 θ i β Γ ( 1 + i β ) X i ( ϑ ) .
Hence,
ψ θ n β = i = 0 1 Γ ( 1 + i β ) θ i β X i + 1 ( ϑ ) = i = 0 1 Γ ( 1 + i β ) θ i β X i + n ( ϑ ) , L β ψ = L β i = 0 θ i β Γ ( 1 + i β ) X i ( ϑ ) = i = 0 θ i β Γ ( 1 + i β ) L β X i ( ϑ ) .
Considering (111), we have
i = 0 1 Γ ( 1 + i β ) θ i β X i + n ( ϑ ) = i = 0 θ i β Γ ( 1 + i β ) ( L β , X i ) ( ϑ ) .
Consequently, we can derive a recursion from (112)
X i + n ( ϑ ) = ( L β , X i ) ( ϑ ) ,
and when n = 1 , the following relation is obtained:
X i + 1 ( ϑ ) = ( L β , X i ) ( ϑ ) .
For n = 2 , we get the following relation:
X i + 2 ( ϑ ) = ( L β , X i ) ( ϑ ) .
The solution of (106) can be obtained through the recursion process,
ψ ( ϑ , θ ) = i = 0 θ i β Γ ( 1 + i β ) X i ( ϑ ) ,
whereas convergence can be achieved using
lim n θ i β Γ ( 1 + i β ) X i ( ϑ ) = 0 .

4.3. Local Fuzzy Fractional Variational Iteration and Decomposition Methods

We study the following nonlinear local fuzzy fractional differential equations to demonstrate two analysis techniques:
L β ( n ) ψ R β ψ = 0 ,
where L β ( n ) represents linear local fuzzy fractional operators with n = 1 , 2 and R β represents linear local fuzzy fractional operators of order less than L β ( n ) .

4.3.1. Local Fuzzy Fractional Variational Iteration Method

In this sub-section, we present the local fuzzy fractional variational iteration technique as follows
ψ n + 1 ( θ ) = ψ n ( θ ) 1 Γ ( 1 + β ) 0 θ λ β Γ ( 1 + β ) L β ( n ) ψ n ( s ) R β ψ n ( s ) ( d s ) β ,
where ψ n is a confined local fuzzy fractional variation, i.e., δ β ψ n = 0 . For n = 0 , we have
λ β = ( s θ ) β Γ ( 1 + β ) ,
and, therefore, this iteration can be defined as
ψ n + 1 ( θ ) = ψ n ( θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) L β ( n ) ψ n ( s ) R β ψ n ( s ) ( d s ) β .
Thus, the solution can be in the following form:
ψ ( ϑ ) = lim n ψ n ( ϑ ) .

4.3.2. Local Fuzzy Fractional Decomposition Method

The local fuzzy fractional differential operator L β ( n ) in Equation (118) has order 2 β , and we define
L β ( n ) = L ϑ ϑ ( 2 β ) = 2 β ϑ 2 β , R β ψ ( θ ) = β ϑ β ψ ( θ ) f ( θ ) .
The n-fold fuzzy fractional integral operator is defined as
L β ( 2 ) m ( s ) = 0 I ϑ β 0 I ϑ β m ( s )
and we have
L β ( 2 ) L β ( 2 ) ψ ( s ) = L β ( 2 ) R β ψ ( s ) .
Hence,
ψ ( s ) = r ( ϑ ) L β ( 2 ) R β ψ ( s ) ,
where r ( ϑ ) is identified by the initial conditions of the fractal. Consequently, the iterative formula is
ψ ( ϑ ) = ψ 0 ( ϑ ) L β ( 2 ) R β ψ ( s ) ,
with ψ 0 ( ϑ ) = r ( ϑ ) . Therefore, assuming that n 0 , the recurrence connectionthat we have achieved is
ψ n + 1 ( ϑ ) = L β ( 2 ) R β ψ n ( s ) , ψ 0 ( ϑ ) = r ( ϑ ) .
The solution can be formulated as
ψ ( ϑ ) = lim n ϕ n ( ϑ ) = lim n n = 0 ψ n ( ϑ ) .

4.4. Applications

In this section, we give examples of fuzzy fractional diffusion and wave equations on Cantor sets to demonstrate the efficiency of the used techniques.
Example 4.
We consider the following fuzzy fractional diffusion equation on a Cantor set
β ψ ( ϑ , θ ) θ β 2 β ψ ( ϑ , θ ) ϑ 2 β = 0 , 0 < β 1 ,
subject to the initial condition
ψ ( ϑ , 0 ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) ,
where n = 1 , 2 , 3 , .
1. Local fuzzy fractional Laplace variational iteration method.
Using relation (103), we obtain
Ł β ψ n + 1 ( ϑ , θ ) = Ł β ψ n ( ϑ , θ ) Ł β 1 Ł β β ψ ( ϑ , θ ) θ β 2 β ψ ( ϑ , θ ) ϑ 2 β = Ł β ψ n ( ϑ , θ ) 1 s β s β Ł β ψ n ( ϑ , θ ) ψ n ( ϑ , 0 ) 2 β Ł β ψ n ( ϑ , θ ) ϑ 2 β = 1 s β ψ n ( ϑ , 0 ) 1 s β 2 β Ł β ψ n ( ϑ , θ ) ϑ 2 β .
The initial value can be expressed as
ψ 0 ( ϑ , θ ) = ψ ( ϑ , 0 ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) .
We obtain the initial approximation
Ł β ψ 1 ( ϑ , θ ) = 1 s β ψ 0 ( ϑ , 0 ) 1 s β 2 β Ł β ψ 0 ( ϑ , θ ) ϑ 2 β , = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n 1 s β ϑ β Γ ( 1 + β ) 1 s β 2 β ϑ 2 β Ł β ϑ β Γ ( 1 + β ) , = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n 1 s β ϑ β Γ ( 1 + β ) .
Therefore,
ψ 1 ( ϑ , θ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n Ł β 1 1 s β ϑ β Γ ( 1 + β ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) .
The second approximation is presented as
Ł β ψ 2 ( ϑ , θ ) = 1 s β ψ 1 ( ϑ , 0 ) 1 s β 2 β Ł β ψ 1 ( ϑ , θ ) ϑ 2 β , = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n 1 s β ϑ β Γ ( 1 + β ) 1 s β 2 β ϑ 2 β Ł β ϑ β Γ ( 1 + β ) , = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n 1 s β ϑ β Γ ( 1 + β ) .
Therefore,
ψ 2 ( ϑ , θ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n Ł β 1 1 s β ϑ β Γ ( 1 + β ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) .
Therefore, the local fuzzy fractional series solution is
ψ ( ϑ , θ ) = lim n Ł β 1 Ł β ψ n ( ϑ , θ ) , = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n lim n Ł β 1 1 s β ϑ β Γ ( 1 + β ) .
Thus, we achieve the required solution as
ψ ( ϑ , θ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n lim n ϑ β Γ ( 1 + β ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) .
2. Local fuzzy fractional series expansion method
Following (114), we have the following recursive formula
X i + 1 ( ϑ ) = 2 β X i ( ϑ ) ϑ 2 β X 0 ( ϑ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) ,
which gives us
X 0 ( ϑ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) , X 1 ( ϑ ) = 0 , X 2 ( ϑ ) = 0 , .
Finally, using Equation (141), we obtain
ψ ( ϑ , θ ) = ( 0.25 + 0.25 ϱ ) n , ( 0.75 + 0.25 ϱ ) n ϑ β Γ ( 1 + β ) .
Table 2, Table 3 and Table 4 show the error term between exact and approximate solutions of Example 4 for different values of β = 0.25, 0.50, 0.75.
In Figure 3, we plot 2 D graphs of the approximation solutions for fuzzy fractional diffusion equations (FFDEs) on the Cantor set. Figure 3a shows that for θ = 1 , ϑ = 3 , n = 1 , β = 0.25 , 0.50 , 0.75 , the FFDEs becomes bounded and closed. Moreover, the ψ ̲ sign shows increasing functions and ψ ¯ denotes decreasing functions on the ϱ -level set of ψ . To reflect the concept of the ϱ -level set, Figure 3a illustrates that the ϱ -level set of the FFDEs is bounded and closed for θ = 1 . Similarly, in Figure 3b,c, we can observe the same explanation of the ϱ -level set closedness and boundedness for Example 4.
Example 5.
We consider the following fuzzy fractional diffusion equation on the Cantor set
β ψ ( ϑ , θ ) θ β ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ ( ϑ , θ ) ϑ 2 β = 0 , 0 < β 1 ,
with the initial condition
ψ ( ϑ , 0 ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) ,
where n = 1 , 2 , 3 , .
1. Local fuzzy fractional Laplace variational iteration method
Applying relation (103), the iterative relation is structured as follows:
Ł β ψ n + 1 ( ϑ , θ ) = Ł β ψ n ( ϑ , θ ) Ł β 1 Ł β β ψ ( ϑ , θ ) θ β ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ ( ϑ , θ ) ϑ 2 β = Ł β ψ n ( ϑ , θ ) 1 s β s β Ł β ψ n ( ϑ , θ ) ψ n ( ϑ , 0 ) ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ ( ϑ , θ ) ϑ 2 β = 1 s β ψ n ( ϑ , 0 ) 1 s β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ n ( ϑ , θ ) ϑ 2 β .
According to (104), the initial value is
ψ 0 ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ψ ( ϑ , 0 ) ϑ β Γ ( 1 + β ) , ψ ( β ) ( ϑ , 0 ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) .
Thus, we obtain
Ł β ψ 1 ( ϑ , θ ) = 1 s β ψ 0 ( ϑ , 0 ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ 0 ( ϑ , θ ) ϑ 2 β , = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n 1 s β ϑ 2 β Γ ( 1 + 2 β ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ϑ 2 β ϑ 2 β Γ ( 1 + 2 β ) , = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 3 β .
Therefore,
ψ 1 ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 2 β = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) .
The second approximation is written as
Ł β ψ 2 ( ϑ , θ ) = 1 s β ψ 1 ( ϑ , 0 ) 1 s β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ 1 ( ϑ , θ ) ϑ 2 β , = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n 1 s β ϑ 2 β Γ ( 1 + 2 β ) 1 s β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ϑ 2 β ϑ 2 β Γ ( 1 + 2 β ) 1 θ β Γ ( 1 + β ) , = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 2 β 1 s 3 β .
Thus, we have
ψ 1 ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 2 β 1 s 3 β = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ β Γ ( 1 + β ) t 2 β Γ ( 1 + 2 β ) ,
which gives us the local fuzzy fractional series solutions
ψ ( ϑ , θ ) = lim n Ł β 1 β ψ n ( ϑ , θ ) , = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n lim n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) k = 0 n 1 s ( n + 1 ) β = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n lim n ϑ 2 β Γ ( 1 + 2 β ) k = 0 n θ k β Γ ( 1 + 2 β ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) k = 0 θ k β Γ ( 1 + 2 β ) .
Hence, we find our desired solution as
ψ ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) E β ( θ β ) .
2. Local fuzzy fractional series expansion method
Following (114), we have
X i + 1 ( ϑ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 2 β X i ( ϑ ) ϑ 2 β , X 0 ( ϑ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) ,
We apply the recursive formula (152) to get
X 0 ( ϑ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) , X 1 ( ϑ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) , X 2 ( ϑ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) , .
As a result of these recursive calculations, one can obtain
ψ ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) k = 0 θ k β Γ ( 1 + 2 β ) , ψ ( ϑ , θ ) = ( 0.75 + 0.25 ϱ ) n , ( 1.25 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) E β ( θ β ) .
Table 5 shows the error term between exact and approximate solutions of Example 5. for β = 0.50 , 1 .
In Figure 4, we plot 2 D graphs of the exact and approximation solutions of local fuzzy fractional diffusion equation on Cantor set. Figure 4a shows that for θ = 8 , ϑ = 1 , β = 1 , n = 2 , the local fuzzy fractional diffusion equation on Cantor set become bounded and closed. In addition, the ψ ̲ sign shows increasing functions and ψ ¯ denotes decreasing functions on the ϱ -level set of ψ . To reflect the concept of the ϱ -level set, Figure 4a illustrates that the ϱ -level set of local fuzzy fractional diffusion equation on Cantor set is bounded and closed for θ = 8 . Similarly, in Figure 4b, we can observe the same explanation of ϱ -level set closedness and boundedness for Example 5.
Example 6.
We consider the following fuzzy fractional wave equation on a Cantor set
2 β ψ ( ϑ , θ ) θ 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ ( ϑ , θ ) ϑ 2 β = 0 , 1 < β 2 ,
subject to the initial value condition
ψ ( ϑ , 0 ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) , β ψ ( ϑ , 0 ) θ β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ,
where n = 1 , 2 , 3 , .
1. Local fuzzy fractional Laplace variational iteration method
Using relation (103), we structure the iterative relation as
Ł β ψ n + 1 ( ϑ , θ ) = Ł β ψ n ( ϑ , θ ) Ł β θ β Γ ( 1 + β ) Ł β 2 β ψ ( ϑ , θ ) θ 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ ( ϑ , θ ) ϑ 2 β = Ł β ψ n ( ϑ , θ ) 1 s 2 β s 2 β Ł β ψ n ( ϑ , θ ) s β ψ n ( ϑ , 0 ) ψ n ( β ) ( ϑ , 0 ) ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ n ( ϑ , θ ) ϑ 2 β = 1 s β ψ n ( ϑ , 0 ) 1 s 2 β ψ n ( β ) ( ϑ , 0 ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ n ( ϑ , θ ) ϑ 2 β .
According to (104), the initial value is
ψ 0 ( ϑ , θ ) = ψ ( ϑ , θ ) ψ ( β ) ( ϑ , 0 ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) .
Hence, we obtain the first approximation
Ł β ψ 1 ( ϑ , θ ) = 1 s β ψ 0 ( ϑ , 0 ) 1 s 2 β ψ 0 ( β ) ( ϑ , 0 ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ 0 ( ϑ , θ ) ϑ 2 β , = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n 1 s β ϑ 2 β Γ ( 1 + 2 β ) 1 s β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ϑ 2 β ϑ 2 β Γ ( 1 + 2 β ) , = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 3 β .
Thus,
ψ 1 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 3 β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) .
The second approximation can be written as
Ł β ψ 2 ( ϑ , θ ) = 1 s β ψ 1 ( ϑ , 0 ) 1 s 2 β ψ 1 ( β ) ( ϑ , 0 ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β Ł β ψ 1 ( ϑ , θ ) ϑ 2 β , = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n 1 s β ϑ 2 β Γ ( 1 + 2 β ) 1 s 2 β ϑ 2 β Γ ( 1 + 2 β ) 2 β ϑ 2 β Ł β ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) , = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 2 β ) 1 s β 1 s 3 β 1 s 5 β .
We obtain
ψ 2 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) 1 s β 1 s 3 β 1 s 5 β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) .
The local fuzzy fractional series solutions are
ψ ( ϑ , θ ) = lim n Ł β 1 β ψ n ( ϑ , θ ) , = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n lim n Ł β 1 ϑ 2 β Γ ( 1 + 2 β ) k = 0 n 1 s ( 2 n + 1 ) β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n lim n ϑ 2 β Γ ( 1 + 2 β ) k = 0 n θ 2 k β Γ ( 1 + 2 k β ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) k = 0 θ k β Γ ( 1 + k β ) .
Thus, we can obtain the solutions as
ψ ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) cosh β ( θ β ) .
2. Local fuzzy fractional variation iteration method
Using (121), we have the iterative formula
ψ n + 1 ( ϑ , θ ) = ψ n ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) 2 β ψ n ( ϑ , θ ) s 2 β ( d s ) β 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ n ( ϑ , θ ) ϑ 2 β ( d s ) β ,
where the initial condition is defined as
ψ 0 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) .
After computing (156), we have
ψ 1 ( ϑ , θ ) = ψ 0 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) 2 β ψ 0 ( ϑ , θ ) s 2 β ( d s ) β 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ 0 ( ϑ , θ ) ϑ 2 β ( d s ) β = ψ 0 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) ( d s ) β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) ,
ψ 2 ( ϑ , θ ) = ψ 1 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) 2 β ψ 1 ( ϑ , θ ) s 2 β ( d s ) β 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ 1 ( ϑ , θ ) ϑ 2 β ( d s ) β = ψ 1 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) s 2 β Γ ( 1 + 2 β ) ( d s ) β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) 1 θ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) ,
ψ 3 ( ϑ , θ ) = ψ 2 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) 2 β ψ 2 ( ϑ , θ ) s 2 β ( d s ) β 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ 2 ( ϑ , θ ) ϑ 2 β ( d s ) β = ψ 2 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) ( d s ) β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) i = 0 3 θ 2 i β Γ ( 1 + 2 i β ) ,
ψ n ( ϑ , θ ) = ψ n 1 ( ϑ , θ ) 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) 2 β ψ n 1 ( ϑ , θ ) s 2 β ( d s ) β 1 Γ ( 1 + β ) 0 θ ( s θ ) β Γ ( 1 + β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ n 1 ( ϑ , θ ) ϑ 2 β ( d s ) β = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) i = 0 n θ 2 i β Γ ( 1 + 2 i β ) .
Consequently, we obtain the solution as follows:
ψ ( ϑ , θ ) = lim n ψ n ( ϑ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) cosh β ( θ β ) .
3. Local fuzzy fractional decomposition method
From (128), we obtain
ψ n + 1 ( ϑ , θ ) = 0 I θ ( β ) 0 I θ ( β ) ϑ 2 β Γ ( 1 + 2 β ) 2 β ψ n ( ϑ , θ ) ϑ 2 β , ψ 0 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) .
From (167), the components are
ψ 0 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) , ψ 1 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 2 β Γ ( 1 + 2 β ) , ψ 2 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 4 β Γ ( 1 + 4 β ) , ψ 3 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 3 β Γ ( 1 + 3 β ) , ψ 4 ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 8 β Γ ( 1 + 8 β ) , ψ κ ( ϑ , θ ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) θ 2 κ β Γ ( 1 + 2 κ β ) .
Thus, the exact solution can be defined as
ψ ( ϑ , θ ) = lim n κ = 0 ψ κ ( ϑ , θ ) = lim n κ = 0 ϑ 2 β Γ ( 1 + 2 β ) θ 2 κ β Γ ( 1 + 2 n β ) = ( 2.25 + 0.25 ϱ ) n , ( 2.75 0.25 ϱ ) n ϑ 2 β Γ ( 1 + 2 β ) cosh β ( ϑ β ) ,
where
cosh β ( ϑ β ) = κ = 0 θ 2 n β Γ ( 1 + 2 n β ) .
In Figure 5, we plot 2 D graphs of the approximation solutions for fuzzy fractional wave equations on Cantor set. Figure 5a shows that for θ = 5 , ϑ = 2 , n = 3 , β = 1.25 , 1.50 , 1.75 , the fuzzy fractional wave equation becomes bounded and closed. Moreover, the ψ ̲ sign shows increasing functions and ψ ¯ denotes decreasing functions on the ϱ -level set of ψ . To reflect the concept of the ϱ -level set, Figure 5a illustrates that the ϱ -level set of the fuzzy fractional wave equation is bounded and closed for θ = 5 . Similarly, in Figure 5b,c, we can observe the same explanation of the ϱ -level set closedness and boundedness for Example 6.

5. Conclusions

In this paper, the hybrid LFFHPSTM, LFFSDM, and LFFSVIM have been used to successfully obtain the solutions of LFFPDEs. Computer-based numerical simulations can clarify the fundamental properties of the physical models given in LFFPDEs. The results demonstrate the efficiency and simplicity of the used technique in identifying LFFPDEs. The originality of this work comes from the fact that the utilized methodology has not recently been used for LFFPDEs. The combination that is used has two features; He’s polynomials are used to decompose nonlinear terms in the first case, while fast convergent series solutions in closed form are produced in the second. Furthermore, the LFHPSTM does not necessitate the calculation of complex Adomian polynomials. The LFFHPM and local fuzzy fractional Sumudu transform method (LFFSTM) can be coupled more quickly and provide better mathematical computations than the LFFHPM. The process of determining a solution demonstrates the efficacy and precision of the proposed method. The LFFSDM offers the solution in the form of a series that, if one exists, quickly converges to the precise solution. It is evident from the findings that the LFFSDM produces very accurate solutions with a minimal number of iterations. Moreover, we have investigated the LFFLIM, LFFSEM, LFFVIM, and LFFADM applied to solving fuzzy diffusion and wave equations defined on Cantor sets with fractal conditions. We show that the LFFVIM utilizing the iteration of the correction local fractional function yields numerous successive approximations. However, the exact solution, which is a local fractional continuous function, is provided by the LFFADM, where these components are local fractional continuous functions. The techniques produce approximate solutions to linear and nonlinear fuzzy fractional differential equations.

Author Contributions

Conceptualization, M.O.; Formal analysis, M.O. and L.L.; Funding acquisition M.O.; Investigation, M.M., S.O.S. and L.L.; Methodology, M.O. and M.S.; Validation, M.O. and E.B.; writing original draft, M.O.; Writing—Review and editing M.O. and A.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by Zhejiang Normal University Research Fund under Grant ZC304022909.

Data Availability Statement

No data was used for the research in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zadeh, L.A. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef]
  2. Zadeh, L.A. The concept of a linguistic variable and its application to approximate reasoning-I. Inf. Sci. 1975, 8, 199–249. [Google Scholar] [CrossRef]
  3. Zadeh, L.A. The concept of a linguistic variable and its application to approximate reasoning-II. Inf. Sci. 1975, 8, 301–357. [Google Scholar] [CrossRef]
  4. Zadeh, L.A. The concept of a linguistic variable and its application to approximate reasoning-III. Inf. Sci. 1975, 9, 43–80. [Google Scholar] [CrossRef]
  5. Zadeh, L.A. Toward a generalized theory of uncertainty (GTU)—An outline. Inf. Sci. 2005, 172, 1–40. [Google Scholar] [CrossRef]
  6. Zadeh, L.A. Is there a need for fuzzy logici? Inf. Sci. 2008, 178, 2751–2779. [Google Scholar] [CrossRef]
  7. Negoita, C.V.; Ralescu, D. Applications of Fuzzy Sets to Systems Analysis; Wiley: New York, NY, USA, 1975. [Google Scholar]
  8. Aytar, S.; Pehlivan, S. Statistically monotonic and statistically bounded sequences of fuzzy numbers. Inf. Sci. 2006, 176, 734–744. [Google Scholar] [CrossRef]
  9. Aytar, S. Statistical limit points of sequences of fuzzy numbers. Inf. Sci. 2004, 165, 129–138. [Google Scholar] [CrossRef]
  10. Osman, M.; Xia, Y. Solving fuzzy fractional q-differential equations via fuzzy q-differential transform. J. Intell. Fuzzy Syst. 2023, 44, 2791–2846. [Google Scholar] [CrossRef]
  11. Allahviranloo, T.; Ghanbari, B. On the fuzzy fractional differential equation with interval Atangana-Baleanu fractional derivative approach. Chaos Solitons Fractals 2020, 130, 109397. [Google Scholar] [CrossRef]
  12. Allahviranloo, T.; Gouyandeh, Z.; Armand, A. Fuzzy fractional differential equations under generalized fuzzy Caputo derivative. J. Intell. Fuzzy Syst. 2014, 26, 1481–1490. [Google Scholar] [CrossRef]
  13. An, T.V.; Hoa, N.V. Fuzzy differential equations with Riemann-Liouville generalized fractional integrable impulses. Fuzzy Sets Syst. 2022, 429, 74–100. [Google Scholar] [CrossRef]
  14. Agarwal, R.P.; Arshad, S.; O’Regan, D.; Lupulescu, V. Fuzzy frac-tional integral equations under compactness type condition. Fract. Calc. Appl. Anal. 2012, 15, 572–590. [Google Scholar] [CrossRef]
  15. Alinezhad, M.; Allahviranloo, T. On the Solution of Fuzzy Fractional Optimal Control Problems with the Caputo Derivative. Inf. Sci. 2017, 421, 218–236. [Google Scholar] [CrossRef]
  16. Dai, R.; Chen, M. The structure stability of periodic solutions for first-order uncertain dynamical systems. Fuzzy Sets Syst. 2020, 400, 134–146. [Google Scholar] [CrossRef]
  17. Dong, N.P.; Long, H.V.; Giang, N.L. The fuzzy fractional SIQR model of com-puter virus propagation in wireless sensor network using Caputo Atangana-Baleanu derivatives. Fuzzy Sets Syst. 2021, 429, 28–59. [Google Scholar] [CrossRef]
  18. Long, H.V.; Son, N.T.K.; Tam, H.T.T.; Yao, J.C. Ulam stability for fractional partial integro-differential equation with uncertainty. Acta Math. Vietnam 2017, 42, 675–700. [Google Scholar] [CrossRef]
  19. Liu, Y.; Zhu, Y.; Lu, Z. On Caputo-Hadamard uncertain fractional differential equations. Chaos Solitons Fractals 2021, 146, 110894. [Google Scholar] [CrossRef]
  20. Mazandarani, M.; Kamyad, A.V. Modified fractional Euler method for solving fuzzy fractional initial value problem. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 12–21. [Google Scholar] [CrossRef]
  21. Mazandarani, M.; Najariyan, M. Type-2 fuzzy fractional derivatives. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2354–2372. [Google Scholar] [CrossRef]
  22. Osman, M.; Xia, Y. Solving fuzzy fractional differential equations with applications. Alex. Eng. J. 2023, 69, 529–559. [Google Scholar] [CrossRef]
  23. Kolwankar, K.M.; Gangal, A.D. Fractional differentiability of nowhere differentiable functions and dimension. Chaos 1996, 6, 505–513. [Google Scholar] [CrossRef] [PubMed]
  24. Yang, X.J. Advanced Local Fractional Calculus and Its Applications; World Science: New York, NY, USA, 2012. [Google Scholar]
  25. Yang, Y.J.; Hua, L.Q. Variational iteration transform method for fractional dif-ferential equations with local fractional derivative. Abs. Appl. Anal. 2014, 9, 760957. [Google Scholar]
  26. Golmankhaneh, A.K.; Yang, X.J.; Baleanu, D. Einstein field equations within local fractional calculus. Rom. J. Phys. 2015, 60, 22–31. [Google Scholar]
  27. Yang, X.J.; Baleanu, D.; Zhong, W.P. Approximate solutions for diffusion equations on cantor space-time. Proc. Rom. Acad. Ser. A 2013, 14, 127–133. [Google Scholar]
  28. Yang, X.J.; Srivastava, H.M.; Cattani, C. Local fractional homotopy perturbation method for solving fractal partial differential equations arising in mathematical physics. Rom. Rep. Phys. 2015, 67, 752–761. [Google Scholar]
  29. Yang, A.M.; Li, J.; Srivastava, H.M.; Xie, G.N.; Yang, X.J. Local fractional laplace variational iteration method for solving linear partial differential equations with local fractional derivative. Dis. Dyn. Nat. Soc. 2014, 8, 365981. [Google Scholar] [CrossRef]
  30. Jassim, H.K. Local fractional Laplace decomposition method for nonhomogeneous heat equations arising in fractal heat flow with local fractional derivative. Int. J. Adv. Appl. Math. Mech. 2015, 2, 1–7. [Google Scholar]
  31. Wang, K. Solitary wave dynamics of the Local fractional Bogoyavlen-sky-Konopelchenko model. Fractals 2023, 31, 2350054. [Google Scholar] [CrossRef]
  32. Wang, K. New solitary wave solutions for the fractional Jaulent-Miodek Hierarchy model. Fractals 2023, 31, 2350060. [Google Scholar]
  33. Watugala, G.K. Sumudu transform—A new integral transform to solve differential equations and control engineering problems. Int. J Math. Educ. Sci. Technol. 1998, 24, 35–43. [Google Scholar] [CrossRef]
  34. Weerakoon, S. Complex inversion formula for Sumudutransform. Int. J. Educ. Math. Sci. Technol. 1998, 29, 618–621. [Google Scholar]
  35. Asiru, M.A. Sumudu transform and solution of integralequations of convolution type. Int. Educ. Math. Sci. Technol. 2001, 32, 906–910. [Google Scholar] [CrossRef]
  36. Asiru, M.A. Applications of Sumudu transform to discretedynamic system. Int. J. Educ. Math. Sci. Technol. 2003, 34, 944–949. [Google Scholar] [CrossRef]
  37. Belgacem, F.M.; Karaballi, A.A. Sumudu transform fundamental properties investigation, applications. J. Appl. Math. Stoch. Anal. 2006, 2006, 91083. [Google Scholar] [CrossRef]
  38. Osman, M.; Xia, Y.; Omer, O.A.; Hamoud, A. On the fuzzy solution of linear-nonlinear partial differential equations. Mathematics 2022, 10, 2295. [Google Scholar] [CrossRef]
  39. Osman, M.; Gong, Z.T.; Mustafa, A.M.; Yong, H. Solving fuzzy (1 + n)-dimensional Burgers’ equation. Adv. Differ. Equ. 2021, 219, 219. [Google Scholar] [CrossRef]
  40. Osman, M.; Gong, Z.T.; Mustafa, A.M. Comparison of fuzzy Adomian decomposition method with fuzzy VIM for solving fuzzy heat-like and wave-like equations with variable coefficients. Adv. Differ. Equ. 2020, 327, 327. [Google Scholar] [CrossRef]
  41. Osman, M.; Almahi, A.; Omer, O.A.; Mustafa, A.M.; Altaie, S.A. Approximation solution for fuzzy fractional-order partial differential equations. Fractal Fract. 2022, 6, 646. [Google Scholar] [CrossRef]
  42. Osman, M.; Xia, Y.; Marwan, M.; Omer, O.A. Novel Approaches for solving fuzzy fractional partial differential equations. Fractal Fract. 2022, 6, 656. [Google Scholar] [CrossRef]
  43. Patela, T.; Meher, R. A Study on Temperature Distribution, Efficiency and Effectiveness of longitudinal porous fins by using Adomian Decomposition Sumudu Transform Method. Procedia Eng. 2015, 127, 751–758. [Google Scholar] [CrossRef]
  44. Saadeh, R.; Ahmed, S.A.; Qazza, A.; Elzaki, T.M. Adapting partial differential equations via the modified double ARA-Sumudu decomposition method. Partial. Differ. Equ. Appl. Math. 2023, 8, 100539. [Google Scholar] [CrossRef]
  45. He, J.H. Approximate solution of nonlinear differential equations with convolution product nonlinearities. Comput. Methods Appl. Mech. Eng. 1998, 167, 69–73. [Google Scholar] [CrossRef]
  46. He, J.H. Variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput. 2000, 114, 115–123. [Google Scholar] [CrossRef]
  47. Shah, N.A.; Dassios, I.; El-Zahar, E.R.; Chung, J.D.; Taherifar, S. The Variational It-eration Transform Method for Solving the Time-Fractional Fornberg-Whitham Equation and Comparison with Decomposition Transform Method. Mathematics 2021, 9, 141. [Google Scholar] [CrossRef]
  48. Singh, G.; Singh, I. Semi-analytical solutions of three-dimensional (3D) coupled Burgers’ equations by new Laplace variational iteration method. Part. Differ. Equ. Appl. Math. 2022, 6, 100438. [Google Scholar] [CrossRef]
  49. Prakasha, A.; Kumar, M.; Baleanu, D. A new iterative technique for a fractional model of nonlinear Zakharov-Kuznetsov equations via Sumudu transform. Appl. Math. Comput. 2018, 334, 30–40. [Google Scholar] [CrossRef]
  50. Anac, H.; Merdan, M.; Kesemen, T. Solving for the random component time-fractional partial diferential equations with the new Sumudu transform iterative method. SN Appl. Sci. 2020, 2, 1112. [Google Scholar] [CrossRef]
  51. Bhargava, A.; Jain, D.; Suthar, D.L. Applications of the Laplace variational iteration method to fractional heat like equations. Partial. Differ. Equ. Appl. Math. 2023, 8, 100540. [Google Scholar] [CrossRef]
  52. Nadeem, M.; Li, F.; Ahmad, H. Modified Laplace variational iteration method for solving fourth-order parabolic partial differential equation with variable coefficients. Comput. Math. Appl. 2019, 78, 2052–2062. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Cattani, C.; Yang, X.J. Local fractional homotopy perturbation method for solving non-homogeneous heat conduction equations in fractal domains. Entropy 2015, 17, 6753–6764. [Google Scholar] [CrossRef]
  54. Singh, J.; Kumar, D.; Nieto, J.J. A reliable algorithm for a local fractional tricomi equation arising in fractal transonic flow. Entropy 2016, 18, 206. [Google Scholar] [CrossRef]
  55. Zhao, D.; Singh, J.; Kumar, D.; Rathore, S.; Yang, X.J. An efficient computational technique for local fractional heat conduction equations in fractal media. J. Nonlinear Sci. Appl. 2017, 10, 1478–1486. [Google Scholar] [CrossRef]
  56. Kumar, D.; Tchier, F.; Singh, J.; Baleanu, D. An efficient computational technique for fractal vehicular traffic flow. Entropy 2018, 20, 259. [Google Scholar] [CrossRef]
  57. Prakash, A.; Kaur, H. An efficient hybrid computational technique for solving nonlinear local fractional partial differential equations arising in fractal media. Nonlinear Eng. 2018, 7, 229–235. [Google Scholar] [CrossRef]
  58. Dubey, V.P.; Singh, J.; Alshehri, A.M.; Dubey, S.; Kumar, D. A comparative analysis of two computational schemes for solving local fractional Laplace equations. Math. Meth. Appl. Sci. 2021, 44, 13540–13559. [Google Scholar] [CrossRef]
  59. Srivastava, H.M.; Golmankhaneh, A.K.; Baleanu, D.; Yang, X.J. Local fractional Sumudu transform with application to IVPs on Cantor sets. Abstr. Appl. Anal. 2014, 2014, 620529. [Google Scholar] [CrossRef]
  60. Rivaz, A.; Fard, O.S.; Bidgoli, T.A. Solving fuzzy fractional differential equations by generalized differential transform method. SeMA J. 2016, 73, 149–170. [Google Scholar] [CrossRef]
  61. Kaleva, O. Fuzzy differential equations. Fuzzy Sets Syst. 1987, 24, 301–317. [Google Scholar] [CrossRef]
  62. Wu, C.; Gong, Z. On Henstock integral of fuzzy-number-valued functions (1). Fuzzy Sets Syst. 2001, 120, 523–532. [Google Scholar] [CrossRef]
  63. Yang, H.; Gong, Z. I11-Posedness for fuzzy Fredholm integral equations of the first kind and regularization methods. Fuzzy Sets Syst. 2019, 358, 132–149. [Google Scholar] [CrossRef]
  64. Bede, B.; Stefanini, L. Generalized differentiability of fuzzy-valued functions. Fuzzy Sets Syst. 2013, 230, 119–141. [Google Scholar] [CrossRef]
  65. Bede, B.; Gal, S.G. Generalizations of the differentiable fuzzy number valued functions with applications to fuzzy differential equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar] [CrossRef]
  66. Allahviranloo, T.; Gouyandeha, Z.; Armanda, A.; Hasanoglub, A. On fuzzy solutions for heat equation based on generalized Hukuhara differentiability. Fuzzy Sets Syst. 2015, 265, 1–23. [Google Scholar] [CrossRef]
  67. Salahshour, S.; Allahviranloo, T.; Abbasbandy, S. Solving fuzzy fractional differential equa-tions by fuzzy Laplace transforms. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1372–1381. [Google Scholar] [CrossRef]
  68. Yang, X.J. Local Fractional Functional Analysis and Its Applications; Asian Academic: Hong Kong, China, 2011. [Google Scholar]
  69. Ziane, D.; Baleanu, D.; Belghaba, K.; Cherif, M. Local fractional Sumudu decomposition method for linear partial differential equations with local fractional derivative. J. King Saud. Univ. Sci. 2019, 31, 83–88. [Google Scholar] [CrossRef]
  70. Hu, M.S.; Agarwal, R.P.; Yang, X.J. Local fractional Fourier series with application to wave equation in fractal vibrating string. Abstr. Appl. Anal. 2012, 2012, 567401. [Google Scholar] [CrossRef]
Figure 1. The exact and approximate lower and upper solutions of (46) at (a) θ = 45 , ϑ = 2 , β = 1 , n = 2 (b) θ = 45 , ϑ = 2 , β = 0.5 , n = 2 , respectively. Moreover, (a,b) are 2D figures for the exact and approximate solutions of local fuzzy fractional differential equations on Cantor of W ( ϑ , θ ) in Example 1.
Figure 1. The exact and approximate lower and upper solutions of (46) at (a) θ = 45 , ϑ = 2 , β = 1 , n = 2 (b) θ = 45 , ϑ = 2 , β = 0.5 , n = 2 , respectively. Moreover, (a,b) are 2D figures for the exact and approximate solutions of local fuzzy fractional differential equations on Cantor of W ( ϑ , θ ) in Example 1.
Fractalfract 07 00851 g001
Figure 2. The exact and approximate lower and upper solutions of (61) at (a) θ = 30 , ϑ = 3 , β = 1 , n = 1 (b) θ = 30 , ϑ = 3 , β = 0.5 , n = 1 . Further, (a,b) are 2D figures for the exact and approximate solutions of local fuzzy fractional partial differential equations with local fuzzy fractional differential operators of W ( ϑ , θ ) in Example 2.
Figure 2. The exact and approximate lower and upper solutions of (61) at (a) θ = 30 , ϑ = 3 , β = 1 , n = 1 (b) θ = 30 , ϑ = 3 , β = 0.5 , n = 1 . Further, (a,b) are 2D figures for the exact and approximate solutions of local fuzzy fractional partial differential equations with local fuzzy fractional differential operators of W ( ϑ , θ ) in Example 2.
Fractalfract 07 00851 g002
Figure 3. The approximate lower and upper solutions of (130) at (ac) θ = 1 , ϑ = 3 , n = 1 , β = 0.25 , 0.50, 0.75, respectively. Moreover, (ac) are 2D figures for the approximate solutions of the fuzzy fractional diffusion equation on the Cantor set of ψ ( ϑ , θ ) in Example 4.
Figure 3. The approximate lower and upper solutions of (130) at (ac) θ = 1 , ϑ = 3 , n = 1 , β = 0.25 , 0.50, 0.75, respectively. Moreover, (ac) are 2D figures for the approximate solutions of the fuzzy fractional diffusion equation on the Cantor set of ψ ( ϑ , θ ) in Example 4.
Fractalfract 07 00851 g003
Figure 4. The exact and approximate lower and upper solutions of (143) at (a) θ = 8 , ϑ = 1 , β = 1 , n = 2 (b) θ = 8 , ϑ = 1 , β = 0.50 , n = 2 . Further, (a,b) 2D figures for the exact and approximate solutions of fuzzy fractional diffusion equation on Cantor set of ψ ( ϑ , θ ) in Example 5.
Figure 4. The exact and approximate lower and upper solutions of (143) at (a) θ = 8 , ϑ = 1 , β = 1 , n = 2 (b) θ = 8 , ϑ = 1 , β = 0.50 , n = 2 . Further, (a,b) 2D figures for the exact and approximate solutions of fuzzy fractional diffusion equation on Cantor set of ψ ( ϑ , θ ) in Example 5.
Fractalfract 07 00851 g004
Figure 5. The approximate lower and upper solutions of (167) at (ac) θ = 5 , ϑ = 2 , n = 3 , β = 1.25 , 1.50, 1.75, respectively. Moreover, (ac) are 2D figures for the approximate solutions of the fuzzy fractional wave equation on the Cantor set of ψ ( ϑ , θ ) in Example 6.
Figure 5. The approximate lower and upper solutions of (167) at (ac) θ = 5 , ϑ = 2 , n = 3 , β = 1.25 , 1.50, 1.75, respectively. Moreover, (ac) are 2D figures for the approximate solutions of the fuzzy fractional wave equation on the Cantor set of ψ ( ϑ , θ ) in Example 6.
Fractalfract 07 00851 g005
Table 1. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.5 .
Table 1. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.5 .
ϱ Lower E-SLower A-SLower ErrorUpper E-SUpper A-SUpper Error
033,369482.4232,88733370483.4232,887
0.133,369482.4732,88733,370483.3732,887
0.233,369482.5232,88733,370483.3232,887
0.333,369482.5732,88733,370483.2732,887
0.433,369482.6232,88733,370483.2232,887
0.533,370482.6732,88733,370483.1732,887
0.633,370482.7232,88733,370483.1232,887
0.733,370482.7732,88733,370483.0732,887
0.833,370482.8232,88733,370483.0232,887
0.933,370482.8732,88733,370482.9732,887
133,370482.9232,88733,370482.9232,887
Table 2. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.25 .
Table 2. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.25 .
ϱ Lower E-SLower A-SLower ErrorUpper E-SUpper A-SUpper Error
00.750.362990.387012.251.0891.161
0.10.8250.399290.425712.1751.05271.1223
0.20.90.435590.464412.11.01641.0836
0.30.9750.471890.503112.0250.980081.0449
0.41.050.508190.541811.950.943781.0062
0.51.1250.544490.580511.8750.907480.96752
0.61.20.580790.619211.80.871190.92881
0.71.2750.617090.657911.7250.834890.89011
0.81.350.653390.696611.650.798590.85141
0.91.4250.689690.735311.5750.762290.81271
11.50.725990.774011.50.725990.77401
Table 3. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50 .
Table 3. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50 .
ϱ Lower E-SLower A-SLower ErrorUpper E-SUpper A-SUpper Error
00.750.48860.26142.251.46580.78419
0.10.8250.537460.287542.1751.41690.75805
0.20.90.586320.313682.11.36810.73191
0.30.9750.635180.339822.0251.31920.70577
0.41.050.684040.365961.951.27040.67963
0.51.1250.73290.39211.8751.22150.65349
0.61.20.781760.418241.81.17260.62735
0.71.2750.830620.444381.7251.12380.60121
0.81.350.879480.470521.651.07490.57507
0.91.4250.928340.496661.5751.02610.54893
11.50.977210.522791.50.977210.52279
Table 4. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.75 .
Table 4. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.75 .
ϱ Lower E-SLower A-SLower ErrorUpper E-SUpper A-SUpper Error
00.750.620060.129942.251.86020.38981
0.10.8250.682070.142932.1751.79820.37682
0.20.90.744080.155922.11.73620.36382
0.30.9750.806080.168922.0251.67420.35083
0.41.050.868090.181911.951.61220.33784
0.51.1250.930090.194911.8751.55020.32484
0.61.20.99210.20791.81.48820.31185
0.71.2751.05410.220891.7251.42610.29885
0.81.351.11610.233891.651.36410.28586
0.91.4251.17810.246881.5751.30210.27287
11.51.24010.259871.51.24010.25987
Table 5. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50 .
Table 5. For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50 .
ϱ Lower E-SLower A-SLower ErrorUpper E-SUpper A-SUpper Error
0838.39161.01677.392328.9447.251881.6
0.1895.22171.92723.32236.7429.541807.1
0.2953.91183.19770.712146.3412.181734.1
0.31014.5194.82819.642057.8395.191662.6
0.41076.9206.81870.061971.2378.551592.6
0.51141.1219.159221886.4362.271524.1
0.61207.3231.85975.441803.5346.351457.1
0.71275.3244.911030.41722.4330.781391.7
0.81345.2258.331086.81643.3315.581327.7
0.91416.9272.11144.81565.9300.731265.2
11490.5286.241204.21490.5286.241204.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Osman, M.; Marwan, M.; Shah, S.O.; Loudahi, L.; Samar, M.; Bittaye, E.; Mohammed Mustafa, A. Local Fuzzy Fractional Partial Differential Equations in the Realm of Fractal Calculus with Local Fractional Derivatives. Fractal Fract. 2023, 7, 851. https://doi.org/10.3390/fractalfract7120851

AMA Style

Osman M, Marwan M, Shah SO, Loudahi L, Samar M, Bittaye E, Mohammed Mustafa A. Local Fuzzy Fractional Partial Differential Equations in the Realm of Fractal Calculus with Local Fractional Derivatives. Fractal and Fractional. 2023; 7(12):851. https://doi.org/10.3390/fractalfract7120851

Chicago/Turabian Style

Osman, Mawia, Muhammad Marwan, Syed Omar Shah, Lamia Loudahi, Mahvish Samar, Ebrima Bittaye, and Altyeb Mohammed Mustafa. 2023. "Local Fuzzy Fractional Partial Differential Equations in the Realm of Fractal Calculus with Local Fractional Derivatives" Fractal and Fractional 7, no. 12: 851. https://doi.org/10.3390/fractalfract7120851

APA Style

Osman, M., Marwan, M., Shah, S. O., Loudahi, L., Samar, M., Bittaye, E., & Mohammed Mustafa, A. (2023). Local Fuzzy Fractional Partial Differential Equations in the Realm of Fractal Calculus with Local Fractional Derivatives. Fractal and Fractional, 7(12), 851. https://doi.org/10.3390/fractalfract7120851

Article Metrics

Back to TopTop