A Computational Approach for the Calculation of Temperature Distribution in Casting-Mould Heterogeneous System with Fractional Order

The purpose of this paper is to investigate the approximate solution of the casting-mould heterogeneous system with Caputo derivative under the homotopy idea. The symmetry design of the system contains the integer partial differential equations and the fractional-order partial differential equations. We apply Yang transform homotopy perturbation method (𝒴T-HPM) to find the approximate solution of temperature distribution in the casting-mould heterogeneous system. The 𝒴T-HPM is a combined form of Yang transform (𝒴T) and the homotopy perturbation method (HPM) using He's polynomials. Some examples are provided to demonstrate the superiority of the suggested technique. The significant findings reveal that 𝒴T-HPM minimizes the enormous without imposing any assumptions. Due to its powerful and robust support for nonlinear problems, this approach presents a remarkable appearance in the functional studies of fractal calculus.


Introduction
A differential problem of symmetry is a modification that generates the differential equation continuously in such a way that these symmetries can help to achieve the solution of the differential equation. Solving these equations is sometimes easier than solving the original differential equations.
In the past few decades, nonlinear fractional differential equations (FDEs) in mathematical physics have been contending for a prominent role in a variety of fields, including biological research, applied science, signal processing, control theory, finance, and fractal dynamics [1][2][3]. Debnath [4] presents some recent applications of fractional calculus and obtained the numerical computation of fractional derivatives and integrals. Heydari et al. [5] applied the Legendre polynomials to obtain the numerical solution of nonlinear fractal-fractional optimal control problems. Wang and Wang [6] employed a semi-inverse method to obtain the fractal variational principles for two different types of discontinuous plasma physics. FDEs are the generalized forms of the integer-order differential equations but some nonlinear mathematical models of integer-order derivatives do not implement well in most of the circumstances [7][8][9]. This is because integer-order derivatives are limited operators that are unsuitable for infinite variance, whereas the fractional-order derivatives are global to account for neighborhood dominance.
Various types of differential equations with the fractional derivative can be used to precisely characterize many proceedings of physics and engineering. Khan et al. [10] constructed an operator using the Caputo fractional differentiation to validate the performance of this approach. The challenge of discovering approximate and exact solutions to FDEs is much critical. The homotopy perturbation technique (HPM) [11] is a well-known method for obtaining series solutions to a variety of linear and nonlinear differential equations of arbitrary order. Many powerful and efficient strategies have been proposed such as Laplace homotopy perturbation method [12], weighted least squares method [13], iterative method [14], homotopy perturbation Sumudu transform method [15], Elzaki transform decomposition approach [16], Laplace decomposition method [17], and natural homotopy transform method [18] with a logic sensitivity function and small diffusivity. Grzymkowski et al. [19] employed HPM whereas Tripathi and Mishra [20] adopted HPM together with the Laplace transform to determine the temperature distribution in the casting-mould heterogeneous system as a continuous function, which is particularly useful for analyzing the mould. Vanani et al. [21] used a weighted approach based on HPM to solve the heat equation in the cast-mould heterogeneous domain. Later, this proposed approach has also been examined in more than one spatial dimension, indicating that this method has a broader application in nonlinear PDE systems [22,23]. This study is particularly powerful for fractal theory and fractal calculus, and it can be seen as dependable in getting analytical solutions and suitable for other nonlinear issues [24][25][26].
This study presents the idea of YT-HPM to obtain the solution of casting-mould heterogeneous system with fractional order in Caputo sense. Yang transform coupled with the homotopy perturbation method presents the results in the form of series and this series approaches to the exact solution very rapidly. The quality of the current method is appropriate to provide the analytical results to the given examples. This study is summarized as follows: in Section 2, we start with some primary definitions in Caputo sense. In Sections 3 and 4, we formulate the problem for the implementation of YT-HPM. In Section 6, we apply this scheme to two numerical problems to show its capability and efficiency. Results and discussion with concluding remarks are given in Sections 7 and 8.

Preliminary Concepts
In this segment, we demonstrate some fundamental properties of fractional calculus along with Yang transform, which help to construct the idea of YT-HPM. Definition 1. The fractional-order derivative in Caputo sense is given as [27] Definition 2. Recently, Yang [28] introduced the Yang-Laplace transform that if ΨðηÞ is a function, then YT can be written as where Y −1 is the inverse Yang operator.
Definition 5. The Yang transform for fractional-order derivatives is defined as [29]

Remarks
The YT of some helpful expressions are as follows:

The Description of the Problem
In this segment, we formulate the casting-mould system to analyze the temperature distribution. Let us consider, two regions, Ψðθ, ηÞ indicating for casting and Φðθ, ηÞ for mould on the boundary of the problem as shown in Figure 1, such that with the boundaries on these domains δ i , i = 1, 2, 3, 4, 5 are distributed as These functions satisfy the heat conduction equation where ∂ α /∂η α is the derivative of functions Φðθ, ηÞ and Ψð θ, ηÞ order α in Caputo sense, a and b are the thermal diffusivity, Ψ and Φ represent the temperature, and η and θ refer to the time and spatial, respectively [30]. These boundaries satisfy the following initial and boundary conditions: The selection of these boundary conditions is an important task for the determination of the casting-mould problem.

Idea of YT-HPM
In this part, we will demonstrate the concept of YHPTM. Let us assume fractional-order PDE such as where R and N are linear and nonlinear differential operators, respectively, and gðθ, ηÞ is called the source function.
Applying the YT to Equation (11), By using inverse Y, However, HPM is stated as where p is the homotopy parameter and The following strategy can be operated to acquire He's polynomials: With the help of Equations (15) and (16), we can get Equation (14) such as We can get the following terms by evaluating the p components Thus, we can summarize the set of Equation (19) in the

Numerical Examples
Case I: let us consider Thus, system of Equation (9) becomes Now, taking YT and using its property definition, we get Thus, inverse YT takes place as Using the initial condition Equation (21) into Equation (24), we get Applying HPM to get with He's polynomials, we get Start with the initial condition to get the following iteration in the form of series

Computational and Mathematical Methods in Medicine
Consequently, the series may be demonstrated as: which can be written as follows For α = 1, the above equations may reduce to the classical casting system Case II: let us consider again Thus, system of Equation (9) becomes According to YT-HPM, we get

Computational and Mathematical Methods in Medicine
Start with the initial condition to get the following iteration in the form of series Consequently, the series may be demonstrated as: which can be written as follows: For α = 1, the above equations may reduce to the classical casting system  Table 1 represents the absolute error of the casting system, and Table 2 represents absolute error of the mould system. In case II, Figures 6(a)-6(d) indicate the surface solution of casting system with fractional order α = 0:25, α = 0:50, α = 0:75, and α = 1 at θ = 2 and η = 2. Figure 7 indicates the graphical results of Ψðθ, ηÞ for different values of α at θ = 2 and η = 0:1. However, the graphical results for the mould system remain same because only changing the initial condition in casting system is studied to show the performance of this approach in this case. Table 3 represents the absolute error of casting system.
These graphical results reveal that they are virtually similar and validate towards the exact solutions, which encourages us to interpret the physical behavior of the coupled system. The solutions results are demonstrated in both 2D and 3D to realize the physical description of the coupled system.

Conclusion
In this survey, we successfully utilized YT-HPM to investigate the approximate solution of the casting-mould heterogeneous system with the Gerasimov-Caputo derivative. This approach does not involve any hypothesis and restriction of variables to ruin the nature of the problems in the recurrence relation. Two examples are tested to verify the excellent performance of this hybrid scheme. It is seen that YT-HPM has less computational effort which shows that the solution of the system of PDEs has a fast rate of convergence. We performed all the calculations with the help of Wolfram Mathematica software 11.0.1. The graphical representation of surface solution and plot distributing validate that YT-HPM results are very precise and effective which demonstrates that this approach is very simple and straightforward for other nonlinear evolution problems with fractal derivatives in the future demands.

Data Availability
All the data are available within the article.