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

: 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 speciﬁc 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 efﬁcient and simple to use in order to solve fuzzy differential equations with local fractional derivatives.


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].
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.

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 [7].Suppose that ψ, h ∈ E 1 and k ∈ R, and the addition and scalar multiplication are defined by The triangular fuzzy number defined as a fuzzy set in E where cl is the closure of set {ϑ ∈ R|ψ(ϑ) > 0}.
Definition 1 ([60]).Let us consider ψ, h ∈ E 1 .If there exists a w ∈ E 1 such that ψ = h + w, then w is called the Hukuhara difference of ψ, h, and it is denoted using ψ h.Note that The Hausdorff distance between two fuzzy numbers is defined as where the -level sets of ψ and h are [ψ] = (ψ , ψ ) and [h] = (h , h ), 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 ψ : Definition 2 ([62,63]).For arbitrary fuzzy numbers ψ, 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 ψ gH h = w ⇔ ψ = h ⊕ w.The H-difference is unique; however, it does not always exist (a necessary condition for ψ gH h to exist is that ψ contains a translation {c} ⊕ h of h).Definition 3 ([63,64]).The generalized Hukuhara difference between two fuzzy numbers ψ, h ∈ E 1 is given as In terms of the -level, we obtain [ and if the H-difference exists, then ψ h = ψ gH h; the conditions for the existence of w = ψ gH h ∈ E 1 are Case (i) with w increasing, w decreasing, w ≤ w . (2) with w increasing, w decreasing, w ≤ w . ( It is simple to demonstrate that (i) and (ii) are both true provided that w is a crisp number.
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 In the same way, W(ϑ) is referred to as local fractional continuous on (γ 1 , γ 2 ) and is symbolized as Definition 12.The Mittage-Leffler function in fractal media is represented as Definition 13.In fractal media, the trigonometric function is stated as 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 where The local fractional partial derivative of W(ϑ, ξ) ∈ C β (γ 1 , γ 2 ) of order β was provided by Yang [24,68] as follows: where 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]: where ϑ β specifies the Cantor function.

Local Fuzzy Fractional Partial Differential Equations
In this section, the hybrid fuzzy fractional approach is applied to some local fuzzy fractional partial differential equations.

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 Applying the formula above, the inverse LFST is defined as For m = 2, we obtain where Theorem 6.The LFST of some special functions where a is a constant.

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.
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) (17)   Using the property of LFFST, we have where Using Equation ( 19) and the reverse of LFFST, we obtain The basic steps include the following.The expansion of W(ϑ, θ) takes the place of a power series of ρ as and the nonlinear component is presented as where ρ(ϑ, θ) is an embedding variable and H n (W) is a local fractional He's polynomial represented as Applying the values of W(ϑ, θ) and R β W(ϑ, θ) in Equation (20) yields the following homotopy equation: Furthermore, comparing the coefficients of the same powers of ρ yields Finally, the local fractional series solution of Equation ( 16) can be written as 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 ].Proof.The solution of Equation ( 16) Assume that W(ϑ, θ) and Y(ϑ, θ) are the distinct solutions of Equation ( 16); then, Applying the convolution theorem, and U is a bounded operator, and R satisfies the Lipschitz condition with δ > 0 such that Using the value theorem of integral calculus, 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 and Hence, and

Local Fuzzy Fractional Sumudu Decomposition Method
In this section, we present the fuzzy linear operator with a local fractional derivative as follows, 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 LFFS β ) on both sides of Equation (33), we obtain Applying the prosperity of the local fractional Sumudu transform, The unknown function U is then divided into infinite series by Substituting ( 36) into (35), we obtain On comparing Equation (37), we obtain The general form of the local fractional recursive relation is where 0 < v ≤ 1, n ∈ N * and m = 1, 2, 3, . . .

Local Fuzzy Fractional Sumudu Variational Iteration Method
We consider the following nonlinear fuzzy operator with a local fractional derivative as 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 (41)   Applying the property of the LFFST, we obtain Using the inverse LFFST of Equation ( 42) yields Taking ∂ β ∂θ β of Equation ( 43), we obtain The following limit is used to determine the solution

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: with the initial condition where n = 1, 2, 3, . . .

Local fuzzy fractional homotopy perturbation Sumudu transform method
Using Equation (46) and the LFFST operator LFFS, we obtain Applying the LFFST for the LFFD formula to Equation (48) yields Rearranging the terms in Equation (48) yields Further simplification on account of Equation (47) reduces Equation (50) to Using the inverse LFFST of Equation ( 51), we obtain The LFHPM suggests the following homotopy equation: Comparing the like powers of ρ, we obtain After simplification, we get Consequently, the solution of Equation ( 46) is as follows:

Local fuzzy fractional Sumudu decomposition method
The subsequent approximations from ( 39) and ( 46) are The subsequent formula (57) yields The first terms of LFFSDM are given by the equations ( 58) as The local fuzzy fractional series form is Thus, we can obtain the exact solution as follows: In Figure 1, we plot 2D 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: with the initial conditions W(ϑ, 0) = [(3.5 + 0.5 ) n , (4.5 − 0.5 where n = 1, 2, 3, . . .

Local fuzzy fractional homotopy perturbation Sumudu transform method.
Using the LFFST for the LFFD formula in Equation ( 61), we obtain Applying the LFFST for the LFD formula in Equation ( 63), we obtain After the rearrangement of the terms in Equation ( 64), we obtain The initial conditions (62) can transform Equation (68) as Taking the inverse LFFST on Equation (66), we obtain After simplification, Equation (67) yields The LFFHPM proposes homotopy expression creation as The following components for the series solution are calculated by comparing the powers of ρ as , After simplification, we obtain Hence, we can obtain the solutions of Equation ( 61) as follows:

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 Applying the formula (74), successive approximations are obtained as follows: The first terms of LFFSDM have the following form: The local fractional series form can be expressed as Therefore, the exact solution of ( 61), as determined by LFFSDM, is given by Therefore, the exact solution can be achieved as Table 1 shows the error term between exact and approximate solutions of Example 2 for ∈ [0, 1].
Table 1.For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.5.In Figure 2, we plot 2D 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

Local fuzzy fractional homotopy perturbation Sumudu transform method
Applying the LFFST of Equation (79), we obtain the following: Taking the formula of LFFST for the LFFDs in Equation ( 81), we obtain Rearranging the terms given in Equation (82) yields Moreover, Equation ( 83) is simplified for the initial condition (80) in the following manner: We apply the inverse of LFFST on Equation (83) to obtain (W(ϑ, θ)) = [(5.5 + 0.5 ) n , (6.5 − 0.5 The homotopy expression can be established using the fundamental method of LFFHPM: The comparison of the similar powers of ρ yields the subsequent components of the series solution: Following the reduction technique, we have Therefore, the solution to Equation ( 79) is expressed as W(ϑ, θ) = [(5.5 + 0.5 ) n , (6.5 − 0.5 .

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.

Local Fuzzy Fractional Laplace Variational Iteration Method
We consider the following local fuzzy fractional partial differential equations, 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 where ) 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 Applying the Yang-Laplace transform on the value of (93), yields, in particular Assume the local fuzzy fractional variation of Equation ( 95), as expressed by We apply the computation of (96) to obtain Thus, from (96), we get where Consequently, we obtain Using the Yang-Laplace transform's inverse, we have and, as a result, we get Consequently, we have the subsequent iteration algorithm, where the first value can be written as Hence, the solution of the above equation (92) in terms of a local fuzzy fractional series can be expressed as follows:

Local Fuzzy Fractional Series Expansion Method
We consider the following local fractional differential equation 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 ϑ where T i (θ) and X i (ϑ) are local fuzzy fractional continuous functions.Moreover, there exists a non-differential series term where P i represents a coefficient.We can represent the solution in the form presented below: Then, following (109), we obtain Hence, Considering (111), we have Consequently, we can derive a recursion from (112) and when n = 1, the following relation is obtained: For n = 2, we get the following relation: The solution of (106) can be obtained through the recursion process, whereas convergence can be achieved using

Local Fuzzy Fractional Variational Iteration and Decomposition Methods
We study the following nonlinear local fuzzy fractional differential equations to demonstrate two analysis techniques: β represents linear local fuzzy fractional operators with n = 1, 2 and R β represents linear local fuzzy fractional operators of order less than L (n) β .

Local Fuzzy Fractional Variational Iteration Method
In this sub-section, we present the local fuzzy fractional variational iteration technique as follows where ψ n is a confined local fuzzy fractional variation, i.e., δ β ψ n = 0.For n = 0, we have and, therefore, this iteration can be defined as Thus, the solution can be in the following form:

. Local Fuzzy Fractional Decomposition Method
The local fuzzy fractional differential operator L β in Equation ( 118) has order 2β, and we define The n-fold fuzzy fractional integral operator is defined as and we have Hence, where r(ϑ) is identified by the initial conditions of the fractal.Consequently, the iterative formula is with ψ 0 (ϑ) = r(ϑ).Therefore, assuming that n ≥ 0, the recurrence connectionthat we have achieved is The solution can be formulated as

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 subject to the initial condition where n = 1, 2, 3, • • • .

Local fuzzy fractional Laplace variational iteration method.
Using relation (103), we obtain The initial value can be expressed as We obtain the initial approximation Therefore, The second approximation is presented as = [(0.25 + 0.25 ) n , (0.75 + 0.25 = [(0.25 + 0.25 ) n , (0.75 + 0.25 Therefore, Therefore, the local fuzzy fractional series solution is

Local fuzzy fractional series expansion method
Following (114), we have the following recursive formula X 0 (ϑ) = [(0.25 + 0.25 ) n , (0.75 + 0.25 which gives us X 0 (ϑ) = [(0.25 + 0.25 ) n , (0.75 + 0.25  In Figure 3, we plot 2D 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 with the initial condition where n = 1, 2, 3, • • • .

Local fuzzy fractional Laplace variational iteration method
Applying relation (103), the iterative relation is structured as follows: According to (104), the initial value is Thus, we obtain The second approximation is written as Thus, we have which gives us the local fuzzy fractional series solutions Hence, we find our desired solution as

Local fuzzy fractional series expansion method
Following (114), we have We apply the recursive formula (152) to get As a result of these recursive calculations, one can obtain Table 5 shows the error term between exact and approximate solutions of Example 5. for β = 0.50, 1.
Table 5.For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50.In Figure 4, we plot 2D 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

Local fuzzy fractional Laplace variational iteration method
Using relation (103), we structure the iterative relation as According to (104), the initial value is Hence, we obtain the first approximation The second approximation can be written as We obtain The local fuzzy fractional series solutions are    In Figure 5, we plot 2D 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.

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.

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

Theorem 7 .
The solution obtained by the LFFHPSTM of fuzzy partial differential Equation (16) has a unique solution, whenever 0 < ζ < 1.
[67]nition 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}.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].

Table 2 .
For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.25.

Table 3 .
For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.50.

Table 4 .
For the absolute error between exact solutions (E-S) at β = 1 and approximate solutions (A-S) at β = 0.75.