Time-Space Fractional Heat Equation in the Unit Disk

and Applied Analysis 3 (a) (b) Figure 1: Distributions of the heat in the unit disk. It is well known that, if the function f(z) given by (10) is in the class S, then 󵄨󵄨󵄨an 󵄨󵄨󵄨 ≤ n, n ∈ N \ {1} . (11) Equality holds for the Koebe function f (z) = z (1 − z)2 , z ∈ U. (12) Moreover, if the function f(z) given by (10) is in the classC, then 󵄨󵄨󵄨an 󵄨󵄨󵄨 ≤ 1, n ∈ N. (13) Equality holds for the function f (z) = z 1 − z , z ∈ U. (14) A function f ∈ A is called starlike of order μ if it satisfies the following inequality: R{ zf󸀠 (z) f (z) } > μ, (z ∈ U) (15) for some 0 ≤ μ < 1. We denoted this class S(μ). A function f ∈ A is called convex of order μ if it satisfies the following inequality: R{ zf󸀠󸀠 (z) f (z) + 1} > μ, (z ∈ U) (16) for some 0 ≤ μ < 1. We denoted this class C(μ). Note that f ∈ C(μ) if and only if zf󸀠 ∈ S(μ). For α j ∈ C, j = 1, . . . , l and β j ∈ C \ {0, −1, −2, . . .}, j = 1, . . . , m, the generalized hypergeometric function l F m (α 1 , . . . , α l ; β 1 , . . . , β m ; z) is defined by the infinite series l F m (α 1 , . . . , α l ; β 1 , . . . , β m ; z) := ∞


Introduction
Fractional calculus (integrals and derivatives) of any positive order can be considered a branch of mathematical physics which concerns with differential equations, integral equations, and integrodifferential equations, where integrals are of convolution form with weakly singular kernels of power law type.It has prevailed more and more interest in applications in several fields of applied sciences.
Fractional differential equations (real and complex) are viewed as models for nonlinear differential equations.Varieties of them play important roles and tools not only in mathematics but also in physics, dynamical systems, control systems, and engineering to create the mathematical modeling of many physical phenomena.Furthermore, they are employed in social science such as food supplement, climate, and economics.One of these equations is the heat equation of fractional order.Over the last two decades, many authors studied this equation in fractional type.Recently, Vázquez et al. analyzed a set of fractional generalized heat equations [1], Karatay and Bayramoglu considered the numerical solution of a time-fractional heat equation, which is obtained from the standard diffusion equation by replacing the first-order time derivative with the Riemann-Liouville fractional derivative [2], Hu et al. established a version of the Feynman-Kac formula for the multidimensional stochastic heat equation with a multiplicative fractional Brownian sheet [3], Khan and Gondal designed a reliable recipe of homotopy analysis method and Laplace decomposition method, namely, homotopy analysis transform method to solve fuzzy fractional heat and wave equations [4], and Povstenko investigated Caputo time-fractional heat equation where solution is obtained by applying the Laplace and finite Hankel integral transforms [5].
In this paper, we shall study a solution of the timespace fractional heat equation in the unit disk.The fractional time is taken in the sense of the Riemann-Liouville operator while the fractional space is assumed in the Srivastava-Owa operator.Here we shall apply some properties of the univalent functions to determine the upper bound of this solution.The maximal solution is illustrated in terms of the generalized hypergeometric functions.

Fractional Calculus
The idea of the fractional calculus (i.e., calculus of integrals and derivatives of any arbitrary real or complex order) was found over 300 years ago.Abel in 1823 scrutinized the generalized tautochrone problem and for the first time applied fractional calculus techniques in a physical problem.

The Riemann-Liouville
Operators.This section concerns with some preliminaries and notations regarding the Riemann-Liouville operators.The Riemann-Liouville fractional derivative strongly poses the physical interpretation of the initial conditions required for the initial value problems involving fractional differential equations.Moreover, this operator possesses advantages of fast convergence, higher stability, and higher accuracy to derive different types of numerical algorithms [6,7].
Definition 2. The fractional (arbitrary) order derivative of the function  of order 0 ≤  < 1 is defined by Remark 3. From Definitions 1 and 2,  = 0, we have + ,  > −1;  > 0. ( The Leibniz rule is where 2.2.The Srivastava-Owa Operators.In [9], Srivastava and Owa defined and studied fractional operators (derivative and integral) in the complex -plane C for analytic functions.
Definition 4. The fractional derivative of order  is defined, for a function (), by where the function () is analytic in simply connected region of the complex -plane C containing the origin, and the multiplicity of ( − ) − is removed by requiring log( − ) to be real when ( − ) > 0. Furthermore, for  ≤  <  + 1, the fractional differential operator is defined as Definition 5.The fractional integral of order  is defined, for a function (), by where the function () is analytic in simply-connected region of the complex z-plane (C) containing the origin, and the multiplicity of ( − ) −1 is removed by requiring log( − ) to be real when ( − ) > 0.
Note that the Srivastava-Owa operators are the complex version (in a simply-connected region) of the Riemann-Liouville operators.

The Class of Univalent Functions
One of the major branches of complex analysis is univalent function theory: the study of one-to-one analytic functions.A domain  of the complex plane is said to be convex if and only if the line segment joining any two points of  lies entirely in .An analytic, univalent function  in the unit disk  mapping the unit disk onto some convex domain is called a convex function.Moreover, A set  ⊂ C is said to be starlike with respect to the point  0 ∈  if the line segment joining  0 to all points  ∈  lies in .A function () which is analytic and univalent in the unit disk , (0) = 0, and maps  onto a starlike domain with respect to the origin, is called starlike function.
Let A denote the class of functions () normalized by Also, let S and C denote the subclasses of A consisting of functions which are, respectively, univalent and convex in .It is well known that, if the function () given by ( 10) is in the class S, then Equality holds for the Koebe function Moreover, if the function () given by ( 10) is in the class C, Equality holds for the function A function  ∈ A is called starlike of order  if it satisfies the following inequality: for some 0 ≤  < 1.We denoted this class S * ().
A function  ∈ A is called convex of order  if it satisfies the following inequality: for some 0 ≤  < 1.We denoted this class C().Note that  ∈ C() if and only if   ∈ S * ().

Time-Space Fractional Heat Equation
In this section, we will express the maximal solution (a solution which has no extension is called a maximal solution) for the problem and  :  ×  → .The source of the heat is in the unit disk in neighborhood of the origin.Therefore, the distribution of the heat is either in homogeneous state; in this case we have a convex form (Figure 1(a)); or it is nonhomogeneous then we obtain the starlike (univalent) form (Figure 1(b)).
We have the following results.
where the equality holds true for the Koebe function.
Proof.Since  is univalent, then it can be expressed by In view of Remark 6, we impose that Applying (11) implies This completes the proof.Proof.Since  is convex, then it yields Employing (13) implies This completes the proof.
We proceed to determine the maximal solution of (19) using Theorems 7 and 8. (29) Proof.By applying the upper bound of the operator    (, ) in (19), Theorem 7 implies Operating (30) by   and applying Remark 3, we receive Hence this comletes the proof.
Note that any extension of the solution  in (31) should be out the unit disk .Therefore, we obtain some kind of stable solutions.Moreover, the real case can be founded in [1].
In the same manner of Theorem 9, we may obtain the following result which describes the maximal solution for the convex function (homogeneous case).(32)

Conclusion and Discussion
We introduced a method depending on some properties of the geometric function theory (univalent, starlike, and convex) for obtaining the maximal solution of heat equation of fractional order in a complex domain.This solution was obtained for two cases: homogeneous and nonhomogeneous distributions of the source.We realized that this solution can be considered as a positive, bounded, and stable solution in the unit disk.We may apply this new method on wellknown diffusion equations of fractional order such as the fractional wave equation in complex domain.Furthermore, since these solutions are determined in terms of the hypergeometric function and its generalization, then this leads to have a convergence solution.Note that the hypergeometric function involves the Mittag-Leffler function functions.Figure 2 shows nonhomogeneous distribution (starlike), while Figure 3 shows the homogeneous distribution (convex).The values of (, ) are taken as (0.5, 0.5), (0.19, 0.19), (0.78, 0.78), and (0.99, 0.99) for Figures 3(a)-3(d), respectively.

Figure 1 :
Figure 1: Distributions of the heat in the unit disk.