#### Abstract

We will study a maximal solution of the time-space fractional heat equation in complex domain. 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 employ some properties of the univalent functions in the unit disk to determine the upper bound of this solution. The maximal solution is illustrated in terms of the generalized hypergeometric functions.

#### 1. 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 time-space 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.

#### 2. 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.

##### 2.1. 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 1. *The fractional (arbitrary) order integral of the function of order is defined by
When , we write , where denoted the convolution product (see [8]), , and , , and as , where is the delta function.

*Definition 2. *The fractional (arbitrary) order derivative of the function of order is defined by

*Remark 3. *From Definitions 1 and 2, , we have
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 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 containing the origin, and the multiplicity of is removed by requiring to be real when . Furthermore, for , 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 containing the origin, and the multiplicity of is removed by requiring to be real when .

*Remark 6. *Consider,
Recently, the Srivastava-Owa operators are generalized into two fractional parameters [10]. Moreover, stability of admissible functions in the unit disk is defined by using these fractional operators [11–13]. Other studies involving these operators can be found in [8, 14–16].

Note that the Srivastava-Owa operators are the complex version (in a simply-connected region) of the Riemann-Liouville operators.

#### 3. 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 is said to be starlike with respect to the point if the line segment joining to all points lies in . A function which is analytic and univalent in the unit disk , , and maps onto a starlike domain with respect to the origin, is called starlike function.

Let denote the class of functions normalized by Also, let and denote the subclasses of 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 , then Equality holds for the Koebe function Moreover, if the function given by (10) is in the class , then Equality holds for the function A function is called starlike of order if it satisfies the following inequality: for some . We denoted this class .

A function is called convex of order if it satisfies the following inequality: for some . We denoted this class . Note that if and only if .

For , and , the generalized hypergeometric function is defined by the infinite series where is the Pochhammer symbol defined by

#### 4. 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)).

**(a)**

**(b)**

We have the following results.

Theorem 7. *Let be univalent function in the unit disk for all , and let , . Then
**
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.

Theorem 8. *Let be convex function in the unit disk for all , and let be defined as in Theorem 7. Then
*

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

Theorem 9. *Let be univalent function in the unit disk for all . Then (19) has a maximal solution in the domain of the form
*

*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).

Theorem 10. *Let be convex function in the unit disk for all . Then (19) has a maximal solution in the domain of the form
*

#### 5. 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 well-known 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 , and for Figures 3(a)–3(d), respectively.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

#### Acknowledgments

The authors would like to thank the reviewers for their comments on earlier versions of this paper. This research has been funded by university of Malaya, under Grant no. (HIR/132).