Next Article in Journal
Creep Properties of Cylinder Metal Rubber under Static Compression at Elevated Temperatures
Next Article in Special Issue
Modified Exp-Function Method to Find Exact Solutions of Microtubules Nonlinear Dynamics Models
Previous Article in Journal
Compound-Combination Synchronization for Fractional Hyperchaotic Models with Different Orders
Previous Article in Special Issue
Asymptotic Constancy for the Solutions of Caputo Fractional Differential Equations with Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mild Solutions for the Time-Fractional Navier-Stokes Equations with MHD Effects

by
Kinda Abuasbeh
1,*,
Ramsha Shafqat
2,*,
Azmat Ullah Khan Niazi
2,* and
Muath Awadalla
1
1
Department of Mathematics and Statistics, College of Science, King Faisal University, Hafuf 31982, Al Ahsa, Saudi Arabia
2
Department of Mathematics and Statistics, The University of Lahore, Sargodha 40100, Pakistan
*
Authors to whom correspondence should be addressed.
Symmetry 2023, 15(2), 280; https://doi.org/10.3390/sym15020280
Submission received: 22 December 2022 / Revised: 9 January 2023 / Accepted: 16 January 2023 / Published: 19 January 2023
(This article belongs to the Special Issue Differential/Difference Equations and Its Application)

Abstract

:
Recently, various techniques and methods have been employed by mathematicians to solve specific types of fractional differential equations (FDEs) with symmetric properties. The study focuses on Navier-Stokes equations (NSEs) that involve MHD effects with time-fractional derivatives (FDs). The (NSEs) with time-FDs of order β ( 0 , 1 ) are investigated. To facilitate anomalous diffusion in fractal media, mild solutions and Mittag-Leffler functions are used. In H δ , r , the existence, and uniqueness of local and global mild solutions are proved, as well as the symmetric structure created. Moderate local solutions are provided in J r . Moreover, the regularity and existence of classical solutions to the equations in J r . are established and presented.

1. Introduction

The fractional calculus branch of mathematical analysis, cited in references [1,2], is concerned with the analysis of a large number of interpretations of real or complex number powers as defined by the differentiation operator D and the integration operator J, as well as the development of a calculus for these operators. In fluid mechanics, the Navier-Stokes equation is a partial differential equation that is used to simulate the flow of incompressible fluids [3]. The model is an extended version of the one developed by Swiss mathematician Leonhard Euler in the 18th century to describe the motion of incompressible fluids, as shown in [4]. For incompressible fluid, the Navier-Stokes equation (NSE) is
u ς + ( u . ) u μ 2 u = 1 ρ P + f .
The Navier-Stokes equations (NSE) are a family of equations that fundamentally represent how a fluid flows through its environment, since (NSE) explains the movement of an incompressible fluid. For Newtonian fluid flows, which reflect the conservation of momentum and mass [5,6], covering both the lubrication and greasing of ball bearings and large-scale atmospheric movements, the advantage of using (NSE) is that the unsteady (NSE) are directly solved by them, and they are also able to resolve the smallest eddies and temporal scales of turbulence in the flow. Furthermore, they can offer all the data for each instantaneous flow in the flow field. We see that this system has so many occurrences that the existence, regularity and boundary conditions must be explained using the full power of mathematical theories [7,8]. It is interesting to note that Leray conducted an early work that revealed that a boundary value problem for time-dependent (NSE) has an interesting, excellent solution on particular time intervals if the data are properly smooth [9]. Liquid metals and space plasmas are only two examples of the numerous physical substances that fall under the category of Magneto-hydrodynamics (MHD). MHD, also referred to as magneto-fluid dynamics or hydromagnetics, is the examination of the magnetic properties and “behaviour” of electrically conducting fluids. These “magneto-fluids” include, to name a few, plasmas, liquid metals, salt water, and electrolytes, according to Reference [10]. The words “magneto-hydrodynamics” are derived from the words magneto, which means magnetic field, hydro, which means water, and dynamics, which means motion.Hannes Alfven founded the domain of MHD, for in which he was awarded the 1970 Nobel Prize in Physics. MHD equations are a combination of NSE, which describe the motion of fluids and Maxwell’s equations for electromagnetics (electric field and magnetic field); these system of equations are coupled together to form the hydro-magnetic or magneto-hydrodynamic system [11]. In this research, applications of MHD in the medical sciences are categorized into four groups. These areas include MHD applications in simple flow, peristaltic flow, pulsatile flow and drag delivery. These groups’ respective numerical studies are examined and reported separately. It is really worth stating the significance of Leray’s pioneering work in establishing whether the solution of (NS) decays to zero in L 2 as time reaches infinity. His work inspired numerous scholars to examine this topic, and there is now a large and strong body of literature on the subject; we may list a few examples. As a result, we may argue that the decay aspects of this issue are well recognized. As a result, it has piqued the curiosity of researchers during the last few years. The NSE are regarded as crucial mathematical tools for better understanding a variety of real-world problems in disciplines such as thermo-hydraulics, aeronautical sciences, meteorology, the petroleum industry, plasma physics and others. These equations provide a natural characterization of the interaction of a viscous liquid with rigid bodies [12]. Qayyum et al. [13], Rehman et al. [14] and Saeed et al. [15] worked on MHD. Niazi et al. [16], Shafqat et al. [17], Alnahdi [18], Khan [19] and Abuasbeh et al. [20] investigated the existence and uniqueness of the fractional evolution equations. Symmetry analysis is a powerful tool for understanding partial differential equations, especially when working with equations derived from accounting-related mathematical ideas. The secret of nature is symmetry, despite the fact that it is lacking from the majority of natural observations. It is preferable to hide symmetry when unanticipated symmetry-breaking occurs. Finite and infinitesimal symmetry can be divided into two distinct categories, Finite and infinitesimal. Finite symmetry can either be discrete or continuous, while infinitesimal symmetry is always continuous. Discrete finite symmetry refers to symmetry that exists in a finite number of elements, while continuous finite symmetry is a symmetry that is present across all elements in a given system. In the past few decades, fractional calculus has become increasingly important in mathematics. While space is a continuous transformation, natural symmetries such as parity and temporal inversion are discrete. These two types of symmetry provide different perspectives on the same phenomenon and can be used to gain insight into a variety of mathematical problems. Differential equations of fractional order are more appropriate for a variety of physical problems than equations of integer order. FDEs are extensively used in many fields of science, as electrical circuits, viscoelasticity materials, neural networks, engineering, chemistry, control theory, biology, mechanics and physics [21,22,23]. We draw attention to the fact that over the last three years, FDEs have significantly evolved and are a useful tool for describing certain materials and processes [24,25]. In this paper, we used the previously mentioned explanation. By considering the smooth boundary Ω = K in R n for ( n 3 ) , we investigate the time-fractional NSE:
ς β u ν Δ u + ( u . ) u = p + ( σ B 0 2 u ρ ) , ς > 0 , . u = 0 , u | K = 0 , u ( x , ς ) = a x c o s β + b ς s i n β ,
where the CFD of order β ( 0 , 1 ) is denoted by ς β . Motion is magneto-hydrodynamics (MHD), σ is electrical conductivity, B 0 is the magnetic field and ρ is constant due to incompressible fluid. Here, the velocity field at a point x K and time ς > 0 is denoted by u = ( u 1 ( ς , x ) , u 2 ( ς , x ) , u 3 ( ς , x ) , . . . , u n ( ς , x ) ) , p = p ( ς , x ) denotes the pressure term, the kinematics viscosity shown by the symbol ν , ς represents the time, and a = a(x) defines the initial velocity. We set the smooth boundary to be K . This model was modified by replacing the first time derivative with a fractional derivative of order β , where 0 < β 1 . There was the first time a fractional derivative of this order was in the model. The flow of fluid is detected at an angle of 45 , then we have
u ( 0 , x ) = a x c o s 45 = a x 2 ,
ς β u ν Δ u + ( u . ) u = P + ( σ B 0 2 u ρ ) , ς > 0 , . u = 0 , u | K = 0 , u ( 0 , x ) = a x 2 .
Hereby, we apply the Helmholts leray projector P L on Equation (1) to convert the NSE into a time-fractional model.
The operator ν P Δ with Dirichlet boundary conditions is simply the same operator A, just like in the divergence-free function space under consideration. Next, we write (1) in its abstract sense, which is
C D ζ β u ( ς ) = A u + F ( u , u ) + P L ( σ B 0 2 u ρ ) , ς > 0 , u ( 0 ) = a x 2 ,
where F ( u , ν ) = P L ( u . ) ν . If the Stokes function A and the Helmholtz–Leray projection P L seem alike, then the Equation (2) has a similar solution as that of Equation (1). For convenience, we simply write P instead of P L . The goal of this paper is to demonstrate the presence and distinction of moderate global and local problem solutions of (2) in H δ , r . Additionally, we demonstrate the regularity findings, which indicate that there is just one classical solution if P ( σ B 0 2 u ρ ) is Hölder-continuous. In order for A u and C D ς β u ( ς ) to be Hölder-continuous in J r , u ( ς ) has to be such solution.
The essential idea behind MHD is that magnetic fields have the ability to induce currents in conductive fluids that are in motion, which in turn produces forces on the fluid and modifies the magnetic field. The Navier-Stokes equations for fluid dynamics and Maxwell’s equations for electromagnetism combine to provide the set of equations that describe MHD. It is necessary to simultaneously solve these differential equations, either analytically or numerically. The flow of conductive fluid is influenced by a magnetic field. A current that travels through the magnet at a 45 angle as it goes down the screw and enters the magnetic field.

2. Preliminaries

In this section, we set the representations, definitions and introductory information that are used throughout the research [26]. Consider Ω = K = { ( x 1 , x 2 , . . . x n ) : x n > 0 } to be an open subset of R n , where n 3 . Let 1 < r < , then the bounded Hodge projection to ( L r ( K ) ) n on projection P, whose range is the closure of
C σ ( K ) : = { u ( C ( K ) ) n : . u = 0 , u h a s c o m p a c t s u b s p a c e i n K } ,
and its null space is given by the closure of
{ u ( C ( K ) ) n : u = η , η C ( K ) } .
For clear and simple notations, let J r := C σ ( K ) ¯ | . | r , which is a closed subspace of ( L r ( K ) ) n . Additionally, ( W m , r ( K ) ) n is a Sobolev space with the norm | . | m , r . A = ν P Δ stands for the Stokes operator in J r with the domain is D r ( A ) = D r ( Δ ) J r ; we have
D r ( Δ ) = { u ( Z 2 , r ( K ) ) n : u | K = 0 } .
The closed linear operator -A generates the bounded analytic semigroup { e ζ A } on J r . To present our results, we present the traditional power spaces which is associated with -A. For δ > 0 and u J r , describe
A δ u = 1 Γ ( δ ) 0 ς δ 1 e ς A u d ς .
Then, A δ is a one-to-one operator on J r . A δ is the inverse of A δ , and set H δ , r for the range of A δ supplemented by norm for δ > 0 ,
| u | H δ , r = | A δ u | r .
e ς A is extended (or restricted) to a bound analytic semigroup on H δ , r is simple. Let X be a Banach space and J be an R interval. The family of continuous X-valued functions is given by C ( J , X ) . C ϑ ( J , X ) is the family of all Hölder-continuous functions considering exponent ϑ for 0 < ϑ < 1 . Let β ( 0 , 1 ] and ν : [ 0 , ) X . The fractional integral of order β for a function ν with a lower limit of zero is defined as
I ς β ν ( ς ) = 0 ς g β ( ς s ) ν ( s ) d s , ς > 0 .
Let the right-hand side be point-wise defined the interval [ 0 , ) , where g a tends for Riemann Liouville kernel,
g β ( ς ) = ς β 1 Γ ( β ) , ς > 0 .
Furthermore, the CFD operator of order β is denoted by C D ς β ; and defined as
C D ς β ν ( ς ) = d d ς [ I ς 1 β ( ν ( ς ) ν ( 0 ) ) ] = d d ς 0 ς g 1 β ( ς s ) ( ν ( ς ) ν ( 0 ) ) d s , ς > 0 .
Generally, for u : [ 0 , ] × R n R n . The expression for the CFD of the function u with respect to time is
ς β u ( ς , x ) = ς 0 ς g 1 β ( ς s ) ( u ( ς , x ) u ( 0 , x ) ) d s , ς > 0 .
Let us look have s look at Mittag-Leffler functions
E β ( ς β A ) = 0 M ¯ β ( s ) e s ς β A d s , E β , β ( ς β A ) = 0 β s M ¯ β ( s ) e s ς β A d s ,
here M ¯ β ( λ ) stands for Mainardi Wright function, which is defined by
M ¯ β ( λ ) = k = 0 λ m m ! Γ ( 1 β ( 1 + m ) ) .
Proposition 1. 
Here are some properties of the Mittag-Leffler special functions
(i)
E β , β ( ς β A ) = 1 2 π i Γ θ E β , β ( ς β μ ) ( μ I + A ) 1 d μ ;
(ii)
A α E β , β ( ς β A ) = 1 2 π i Γ θ μ α E β , β ( ς β μ ) ( μ I + A ) 1 d μ .
Proof. 
(i) By having 0 β s M ¯ β ( s ) e s ς d s = E β , β ( ς ) with Fubini theorem, we introduce
E β , β ( ς β A ) = 0 β s M ¯ β ( s ) e s ς β A d s = 1 2 π i 0 β s M ¯ β ( s ) Γ θ e μ s ς β ( μ I + A ) 1 d μ d s = 1 2 π i Γ θ E β , β ( ς β μ ) ( μ I + A ) 1 d μ
with the appropriate Γ θ integral route.
(ii) In the same manner,
A α E β , β ( ς β A ) = 0 β s M ¯ β ( s ) A α e s ς β A d s = 1 2 π i 0 β s M ¯ β ( s ) Γ θ μ α e μ s ς β ( μ I + A ) 1 d μ d s = 1 2 π i Γ θ μ α E β , β ( μ ς β ) ( μ I + A ) 1 d μ .
The results are gotten are like. □
Lemma 1. 
([27]). The operators E β ( ς β A ) and E β , β ( ς β A ) are continuous for ς > 0 , in uniform operator topology. Moreover, they are uniformly continuous on [ r , ) for r > 0 .
Lemma 2. 
([28]).
Take 0 < β < 1 , which implies
(i)
u X , lim ς 0 + E β ( ς β A ) u = u ;
(ii)
u D ( A ) and ς > 0 , C D ς β E β ( ς β A ) u = A E β ( ς β A ) u ;
(iii)
u X , E β ( ς β A ) u = ς β 1 A E β , β ( ς β A ) u ;
(iv)
ς > 0 , E β ( ς β A ) u = I ς 1 β ς β 1 E β , β ( ς β A ) u .
We introduce the next lemma for the function h : [ 0 , ) X , before presenting the notion of a mild solution of (2). For this, see [6].
Lemma 3. 
C D ς β u ( ς ) = A u + F ( u , u ) + P ( σ B 0 2 u ρ ) , ς > 0 , u ( 0 ) = a x 2 ,
satisfying solution is
u ( ς ) = a ς 2 + 1 Γ ( β ) 0 ς ( ς s ) β 1 A u ( s ) + F ( u ( s ) , u ( s ) ) + P ( σ B 0 2 u ρ ) d s , f o r ς 0 u ( ς ) = a ς 2 + 1 Γ ( β ) 0 ς ( ς s ) β 1 ( A u ( s ) ) d s + 1 Γ ( β ) 0 ς ( ς s ) β 1 F ( u ( s ) , u ( s ) ) d s + 1 Γ ( β ) 0 ς ( ς s ) β 1 P ( σ B 0 2 u ρ ) d s .
Taking Laplace on both sides
u ( λ ) = a 2 λ 2 + 1 λ β [ A u ( λ ) ] + 1 λ β F ( u , u ( λ ) ) + 1 λ β ( P σ B 0 2 ρ u ( λ ) ) .
Multiplying both sides by λ β
λ β u ( λ ) = a 2 λ β 2 + [ A u ( λ ) ] + F ( u , u ( λ ) ) + ( P σ B 0 2 ρ u ( λ ) ) ( λ β + A ) u ( λ ) = a 2 λ β 2 + F ( u , u ( λ ) ) + ( P σ B 0 2 ρ u ( λ ) ) u ( λ ) = ( λ β + A ) 1 a 2 λ β 2 + ( λ β + A ) 1 F ( u , u ( λ ) ) + ( λ β + A ) 1 ( P σ B 0 2 ρ u ( λ ) ) .
Taking Laplace inverse for both sides
u ( ς ) = a 2 0 ς E β ( ( ς s ) β A ) d s + 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s + 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P σ B 0 2 ρ u ( s ) d s .
Definition 1. 
A function u : [ 0 , ) H δ , r o r ( J r ) is termed the global mild solution of problem (2) in H δ , r , if u C ( [ 0 , ) , H δ , r ) and ς [ 0 , ) ,
u ( ς ) = a 2 0 ς E β ( ( ς s ) β A ) d s + 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s + 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P σ B 0 2 ρ u ( s ) d s .
Definition 2. 
Let 0 < ˜ < . If u C ( [ 0 , ˜ ] , H δ , r ) or C ( [ 0 , ˜ ] , H δ , r ) and u and satisfy (4) for ς [ 0 , ˜ ] , then the function u : [ 0 , ˜ ] H β , r or ( J r ) is called a local mild solution of problem (2) in H β , r or ( J r ) .
Consider the three operators for case ξ , η , ϕ ,
ξ ( ς ) = a 2 0 ς E β ( ( ς s ) β A ) d s , η ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( ς ) ) d s , ϕ ( u , ν ) ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , ν ( s ) ) ) d s .
Definition 3. 
A non-negative measurable function f defined on a measurable set E is integrable if, E f < .
Definition 4. 
Let g be integrable over E , and consider < f n > be a sequence of a measurable function with | f n | g on E and lim n f n = f , i.e,
lim n E f n = E f .
It is clear that every function f n is integrable on E and, furthermore, it follows that lim n f n = f on E and | f n | g on E that | f | g and, hence, f is integrable on E .
Lemma 4. 
Let ( X , | | . | | X ) denote the Banach space, and let a bilinear operator be defined as G : X × X X and a positive real number L such that
| | G ( u , ν ) | | X L | | u | | X | | ν | | X , u , ν X .
Then, for any u 0 X with | | u 0 | | X < 1 4 L , there is just a unique solution u X to the equation u = u 0 + G ( u , u ) .

3. Global and Local Existence in H δ , r

For the existence and unique property of a mild solution to the situation (2) in H δ , r , we provide adequate conditions. For this, we suppose that:
(e) For ς > 0 , ( P σ B 0 2 ρ u ) is continuous and | P ( σ B 0 2 ρ u ( s ) ) | r = 0 ( ς β ( 1 δ ) for 0 < δ < 1 as ς 0 .
Lemma 5. 
Let 1 < r < and δ 1 δ 2 . The existence of constant is C = C ( δ 1 , δ 2 ) such that
| e ς A ν | H δ 2 , r C ς ( δ 2 δ 1 ) | ν | H δ 1 , r , ς > 0 ,
for ν H δ 1 , r . Furthermore,
lim t 0 ς ( δ 2 δ 1 ) | e ς A ν | H δ 2 , r = 0 .
We now examine a fundamental lemma that enables us to demonstrate the final major theorems of this section.
Lemma 6. 
Let 1 < r < and also δ 1 δ 2 . Then, for any ˜ > 0 , there is constant C 1 = C 1 ( δ 1 , δ 2 ) > 0 such that
| E β ( ς β A ) | H δ 2 , r C 1 ς β ( δ 2 δ 1 ) | ν | H δ 1 , r a n d | E β , β ( ς β A ) | H δ 2 , r C 1 ς β ( δ 2 δ 1 ) | ν | H δ 1 , r ,
for all ν H δ 1 , r and ς ( 0 , ˜ ] . Furthermore, lim ς 0 ς β ( δ 2 δ 1 ) | E β ( ς β A ) ν | H δ 2 , r = 0 .
Proof. 
Let ν H δ 1 , r . By the previous Lemma 5, we find that
| E β ( ς β A ) ν | H δ 2 , r 0 M ¯ β ( s ) | e s ς β A ν | H δ 2 , r d s C 0 M ¯ β ( s ) s ( δ 2 δ 1 ) d s ς β ( δ 2 δ 1 ) | ν | H δ 1 , r C 1 ς β ( δ 2 δ 1 ) | ν | H δ 1 , r .
Additionally, the dominated convergence theorem of Lebesgue demonstrates that
lim ς 0 ς β ( δ 2 δ 1 ) | E β ( ς β A ) ν | H δ 2 , r 0 M ¯ β ( s ) lim ς 0 ς β ( δ 2 δ 1 ) | e s ς β A ν | H β 2 , r d s = 0 .
Similarly,
| E β , β ( ς β A ) ν | H δ 2 , r 0 β s M ¯ β ( s ) | e s ς β A ν | H δ 2 , r d s β C 0 M ¯ β ( s ) s 1 ( δ 2 δ 1 ) d s ς β ( δ 2 δ 1 ) | ν | H δ 1 , r C 1 ς β ( δ 2 δ 1 ) | ν | H δ 1 , r ,
where constant C 1 = C 1 ( β , δ 1 , δ 2 ) is
C 1 C m a x Γ ( 1 δ 2 + δ 1 ) Γ ( 1 + β ( δ 1 δ 2 ) ) , β Γ ( 2 δ 2 + δ 1 ) Γ ( 1 + β ( 1 + δ 1 δ 2 ) ) .

4. Global Existence in H δ , r

Now, a portion of the above section of this article deals with existence of a global mild solution of problem (2) in H δ , r , and we let
M ( ς ) = sup s ( 0 , ς ] s β ( 1 δ ) | P ( σ B 0 2 ρ u ( s ) ) | r , B 1 = C 1 m a x B ( β ( 1 δ ) ) , 1 ( β ( 1 δ ) ) , B ( β ( 1 α ) , 1 β ( 1 δ ) ) , L M C 1 m a x B ( β ( 1 δ ) ) , 1 2 β ( α δ ) ) , B ( β ( 1 α ) , 1 2 β ( α δ ) ) ,
assuming M is provided afterward.
Theorem 1. 
Let 1 < r < , 0 < δ < 1 and (e) hold. For all a ς 2 H δ , r , suppose that
C 1 | a ς 2 | H δ , r + B 1 M < 1 4 L .
The above-mentioned M is defined as M : = sup s ( 0 , ) s β ( 1 δ ) P ( σ B 0 2 u ( s ) ρ ) . If n 2 r 1 2 < δ , then there is a α > m a x { δ , 1 2 } and a function which is unique u : [ 0 , ) H δ , r helps to satisfy:
(a)
The function u : [ 0 , ) H δ , r i s c o n t i n u o u s a n d u ( 0 ) = a ς 2 ;
(b)
The function u : ( 0 , ) H α , r i s c o n t i n u o u s a n d l i m ς 0 ς β ( α δ ) | u ( ς ) | H α , r = 0 ;
(c)
u s a t i s f y ( 4 ) f o r ς [ 0 , ) .
Proof. 
Take α = 1 + δ 2 . Here, X is a space containing all the well-defined curves u : ( 0 , ) H δ , r , also X = X [ ] , and the X term is a complete metric space, and it is nonempty:
(i)
The continuous and bounded function is u : [ 0 , ) H δ , r ;
(ii)
Additionally, the function is continuous and bounded u : ( 0 , ) H α , r ; moreover,
lim ς 0 ς β ( α δ ) | u ( ς ) | H α , r = 0 ;
having a fundamental norm
| | u | | X = m a x sup ς 0 | u ( ς ) | H δ , r , sup ς 0 ς β ( α δ ) | u ( ς ) | H α , r .
Since it is clear the mapping F : H α , r × H α , r J r is a well-defined, bounded as well as bilinear mapping because of a Weissler argument, ∃ M in such a way that for u , ν H α , r ,
| F ( u , ν ) | r M | u | H α , ς | ν | H α , r , | F ( u , u ) F ( ν , ν ) | r M ( | u | H α , r + | ν | H α , r ) | u ν | H α , r .
Step I
Let us suppose that u , ν X . Here, the operator ϕ ( u ( ς ) , ν ( ς ) ) C ( [ 0 , ) , H δ , r ) and also the operator ϕ ( u ( ς ) , ν ( ς ) ) C ( ( 0 , ) , H α , r ) . Consider ς < ς 0 for completely arbitrary, ς 0 0 is fixed and ε > 0 is very small (the following situation is related). There are
| ϕ ( u ( ς ) ν ( ς ) ) ϕ ( u ( ς 0 ) ν ( ς 0 ) ) | H δ , r ς 0 ς ( ς s ) β 1 | E β , β ( ( ς s ) β A ) F ( u ( s ) ν ( s ) | H δ , r d s + 0 ς 0 | ( ( ς s ) β 1 ( ς 0 s ) β 1 ) E β , β ( ( ς s ) β A ) F ( u ( s ) ν ( s ) | H δ , r d s + 0 ς 0 ε ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) F ( u ( s ) , ν ( s ) | H δ , r d s + ς 0 ε ς 0 ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) F ( u ( s ) , ν ( s ) | H δ , r d s : = G 11 ( ς ) + G 22 ( ς ) + G 33 ( ς ) + G 44 ( ς ) .
Each of these four terms is estimated independently. For G 11 ( ς ) , in the light of the above Lemma 6, we find
G 11 ( ς ) C 1 ς 0 ς ( ς s ) β ( 1 δ ) 1 | F ( u ( s ) , ν ( s ) | r d s M C 1 ς 0 ς ( ς s ) β ( 1 δ ) 1 | u ( s ) | H α , r | ν ( s ) | H α , r d s M C 1 ς 0 ς ( ς s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s sup s [ 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r | ν ( s ) | H α , r = M C 1 ς 0 / ς 1 ( 1 s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s sup s [ 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r | ν ( s ) | H α , r .
There exists δ ˜ > 0 , and δ ˜ is very small for 0 < ς ς 0 < δ ˜ from definition of the β function
ς 0 / ς 1 ( 1 s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s 0 .
Consequently, G 11 ( ς ) approaches to 0 as ς ς 0 approaches to 0. For G 12 ( ς ) ,
G 12 ( ς ) C 1 0 ς 0 ( ς 0 s ) β 1 ( ς s ) β 1 ( ς s ) β δ | F ( u ( s ) , ν ( s ) | r d s M C 1 0 ς 0 ( ς 0 s ) β 1 ( ς s ) β 1 ( ς s ) β δ s 2 β ( α δ ) d s sup s [ 0 , ς 0 ] s 2 β ( α δ ) | ( u ( s ) | H α , r | ν ( s ) | H α , r .
It is interesting to note
0 ς 0 | ( ς 0 s ) β 1 ( ς s ) β 1 | ( ς s ) β δ s 2 β ( α δ ) d s 0 ς 0 ( ς s ) β 1 ( ς s ) β δ s 2 β ( α δ ) d s + 0 ς 0 ( ς 0 s ) β 1 ( ς s ) β δ s 2 β ( α δ ) d s 2 0 ς 0 ( ς 0 s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s = 2 B ( β ( 1 δ ) , 1 2 β ( α δ ) ) .
So, by theorem (LDC), there are
0 ς 0 ( ς 0 s ) β 1 ( ς s ) β 1 ( ς s ) β δ s 2 β ( α δ ) d s 0 a s ς ς 0 .
We conclude the limiting value of G 12 ( ς ) is equal to zero as ς ς 0 . Now, we move towards G 13 ( ς ) ,
G 13 ( ς ) 0 ς 0 ε ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) + E β , β ( ( ς 0 s ) β A ) F ( u ( s ) , ν ( s ) | H δ , r d s 0 ς 0 ε ( ς 0 s ) β 1 ( ς s ) α β + ( ς 0 s ) α β | F ( u ( s ) , ν ( s ) | H δ , r d s 2 M C 1 0 ς 0 ε ( ς 0 s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s sup s [ 0 , ς 0 ] s 2 β ( α δ ) | u ( s ) | H α , r | ν ( s ) | H α , r ,
by using the (LDC) theorem one more time, and the operator E β , β ( ς β A ) is uniform continuous by Lemma 1, which shows
lim ς ς 0 G 13 ( ς ) = 0 ς 0 ε ( ς 0 s ) β 1 lim ς ς 0 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) F ( u ( s ) , ν ( s ) ) | H δ , r d s = 0 .
For G 14 ( ς ) , from calculations, we find conclusions that
G 14 ( ς ) ς 0 ε ς 0 ( ς 0 s ) β 1 ( ς s ) β δ + ( ς 0 s ) β δ | F ( u ( s ) , ν ( s ) | r d s 2 M C 1 ς 0 ε ς 0 ( t 0 s ) β ( 1 δ ) 1 s 2 β ( α δ ) sup s [ ς 0 ε , ς 0 ] s 2 β ( α δ ) | ( u ( s ) | H δ , r | ν ( s ) | H δ , r 0 , a s ε 0
From the characteristics of β -function, we find that
| ϕ ( u ( ς ) , ν ( ς ) ) ϕ ( u ( ς 0 ) , ν ( ς 0 ) ) | H δ , r 0 a s ς ς 0 .
The continuous operator ϕ ( u , ν ) is calculated in C ( ( 0 , ) , H α , r ) . The same conversation as before follows. So, we skip the explanation.
Step II
This must prove that ϕ : X × X X is a bilinear, as well as continuous, operator. By Lemma 6, we have
| ϕ ( u ( ς ) , ν ( ς ) ) | H δ , r | 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , ν ( s ) d s | H α ς , r C 1 0 ς ( ς s ) β ( 1 δ ) 1 | F ( u ( s ) , ν ( s ) | r d s M C 1 0 ς ( ς s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s sup s [ 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r | ν ( s ) | H δ , r = M C 1 B ( β ( 1 δ ) , 1 2 β ( α δ ) ) | | u | | X | | ν | | X .
In addition to
| ϕ ( u ( ς ) , ν ( ς ) ) | H α , r | 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , ν ( s ) ) d s | H α , r C 1 0 ς ( ς s ) β ( 1 α ) 1 | F ( u ( s ) , ν ( s ) | r d s M C 1 0 ς ( ς s ) β ( 1 α ) 1 s 2 β ( α δ ) d s sup s [ 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r | ν ( s ) | H α , r = M C 1 ς β ( α δ ) B ( β ( 1 α ) , 1 2 β ( α δ ) ) | | u | | X | | ν | | X .
Hence,
sup ς [ 0 , ) ς β ( α δ ) | ϕ ( u ( ς ) , ν ( ς ) ) | H α , r M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) | | u | | X | | ν | | X .
To be more accurate,
lim ς ς 0 ς β ( α δ ) | ϕ ( u ( ς ) , ν ( ς ) ) | H α , r = 0 .
Thus, ϕ ( u , ν ) belongs to X , and | | ϕ ( u ( ς ) , ν ( ς ) ) | | X L | | u | | X | | ν | | X .
Step III
Let 0 ς ς 0 . Since
| η ( ς ) η ( ς 0 ) | H δ , r ς 0 ς ( ς s ) β 1 | E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H δ , r d s + 0 ς 0 ( ( ς 0 s ) β 1 ( ς s ) ) β 1 | E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H δ , r d s + 0 ς 0 ε ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H δ , r d s + ς 0 ε ς 0 ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H δ , r d s C 1 ς 0 ς ( ς s ) β ( 1 δ ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s + C 1 0 ς 0 ( ς 0 s ) β 1 ( ς s ) β 1 ( ς s ) β δ | P ( σ B 0 2 ρ u ( s ) ) | r d s + C 1 0 ς 0 ε ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H δ , r d s
+ 2 C 1 ς 0 ε ς 0 ( ς 0 s ) β ( 1 δ ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s M ( ς ) C 1 ς 0 ς ( ς s ) β ( 1 δ ) 1 s β ( 1 δ ) d s + M ( ς ) C 1 0 ς 0 ( ( ς 0 s ) β 1 ( ς s ) β 1 ) ( ς s ) β δ s β ( 1 δ ) d s + M ( ς ) C 1 0 ς 0 ε ( ς 0 s ) β 1 | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) s β ( 1 δ ) d s + 2 M ( ς ) C 1 ς 0 ε ς 0 ( ς 0 s ) β ( 1 δ ) 1 s β ( 1 δ ) d s .
From the result of Lemma 1 as well as the property of β function, the 1st, 2nd, 3rd and final integral approaches to 0 as ς ς 0 as ε tends to 0, which suggests
| η ( ς ) η ( ς 0 ) | H δ , r 0 a s ς ς 0 .
We evaluated that η ( ς ) is continuous in H α , r , which is implied by the similar prior explanation.
| η ( ς ) | H δ , r | 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) d s | H δ , r C 1 0 ς ( ς s ) β ( 1 δ ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s M ( ς ) C 1 0 ς ( ς s ) β ( 1 δ ) 1 s β ( 1 δ ) d s = M ( ς ) C 1 B ( β ( 1 δ ) , 1 β ( 1 δ ) ) ,
and also
| η ( ς ) | H α , r | 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) d s | H α , r C 1 0 ς ( ς s ) β ( 1 α ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s M ( ς ) C 1 0 ς ( ς s ) β ( 1 α ) 1 s β ( 1 δ ) d s = ς β ( α δ ) C 1 M ( ς ) B ( β ( 1 α ) , 1 β ( 1 δ ) ) .
Specifically,
ς β ( α δ ) | η ( ς ) | H α , r M ( ς ) C 1 B ( β ( 1 α ) , 1 β ( 1 δ ) ) 0 , a s ς 0 .
As we know, if ς 0 , then M ( ς ) 0 , as a result of an assumption (e). This guarantees that η ( ς ) X and | | η ( ς ) | | B 1 M . For a ς 2 H δ , r . From the statement of Lemma 1, it is simple to see this
E β ( ς β A ) a 2 C ( [ 0 , ) , H δ , r ) E β ( ς β A ) a 2 C ( ( 0 , ) , H α , r ) .
Therefore, also
0 ς E β ( ( ς s ) β A ) a 2 d s C ( [ 0 , ) , H δ , r ) 0 ς E β ( ( ς s ) β A ) a 2 d s C ( ( 0 , ) , H α , r ) .
Lemma 6 suggests that ς ( 0 , ˜ ] ,
0 ς E β ( ( ς s ) β A ) a 2 d s X , ς β ( α δ ) 0 ς E β ( ( ς s ) β A ) a 2 d s C ( ( 0 , ) , H α , r ) , | | 0 ς E β ( ( ς s ) β A ) a 2 d s | | X 0 ς | | E β ( ( ς s ) β A ) a 2 d s | | X , 0 t C 1 | a 2 | H δ , r d s , = C 1 | a ς 2 | H δ , r .
With the help of Equation (5), the inequality gives
| | 0 ς E β ( ( ς s ) β A ) a 2 d s + η ( ς ) | | X | | 0 ς E β ( ( ς s ) β A ) a 2 d s | | X + | | η ( ς ) | | X C 1 | a ς 2 | H δ , r + B 1 M 1 4 L
which is showing the result that F has a unique and special fixed point.
Step IV
For the purpose of demonstrating that u ( ς ) a ς 2 in H δ , r by assigning ς 0 . We must demonstrate this as
lim ς 0 ξ ( ς ) = lim ς 0 a 2 0 ς E β ( ς β A ) d s = 0 , lim ς 0 η ( ς ) = lim ς 0 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P σ B 0 2 ρ u ( s ) d s = 0 , lim ς 0 ϕ ( u , ν ) ( ς ) = lim ς 0 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s = 0
in H δ , r . It is understood that lim ς ς 0 η ( ς ) = 0 and lim ς ς 0 M ( ς ) = 0 with Equation (7). Additionally,
| 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s | H δ , r C 1 0 ς ( ς s ) β ( 1 δ ) 1 | F ( u ( s ) , u ( s ) ) | r d s M C 1 0 ς ( ς s ) β ( 1 δ ) 1 | u ( s ) | H α , r 2 d s M C 1 0 ς ( ς s ) β ( 1 δ ) 1 | u ( s ) | H α , r 2 d s M C 1 0 ς ( ς s ) β ( 1 δ ) 1 s 2 β ( α δ ) d s sup s ( 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r 2 = M C 1 B ( β ( 1 δ ) , 1 2 β ( α δ ) ) sup s ( 0 , ς ] s 2 β ( α δ ) | u ( s ) | H α , r 2 0 a s ς 0 .

5. Local Existence in H δ , r

This section is further divided into the H δ , r local mild solution to problem (2).
Theorem 2. 
Let 1 < r < , 0 < δ < 1 and (e) exists. Let us suppose that
n 2 r 1 2 < δ .
The existence of a function α > m a x δ , 1 2 knowing that for all a ς 2 H δ , r ˜ > 0 and a distinctive continuous function u : [ 0 , ˜ ] H δ , r such that
(i)
A continuous mapping in H δ , r is defined as u : [ 0 , ˜ ] H δ , r with u ( 0 ) = a t 2 ;
(ii)
A continuous mapping in H α , r is defined as u : ( 0 , ˜ ] H α , r with limiting value of function lim ς 0 ς β ( α δ ) | u ( ς ) | H α , r = 0 ;
(iii)
u holds (4) for ς [ 0 , ˜ ] .
Proof. 
Take α = 1 + δ 2 . Additionally, value a ς 2 H δ , r . Let us define X ˜ as the space of curves in such a way that u : ( 0 , ˜ ] H δ , r , moreover X ˜ = X [ ˜ ] :
( a )
A continuous mapping is defined as u : [ 0 , ˜ ] H δ , r ;
( b )
A continuous mapping with limiting value of function is defined as u : ( 0 , ˜ ] H α , r , with lim ς 0 ς β ( α δ ) | u ( ς ) | H α , r = 0 ,
with a norm defined by
| | u | | X = sup ς [ 0 , ˜ ] ς β ( α δ ) | u ( ς ) | H α , r .
From the proof of Theorem 1, we notice that the operator ϕ is continuous and it is a linear map ϕ : X × X X , and η ( ς ) X . From Lemma 1, we can claim ∀ ς ( 0 , ˜ ] ,
E β ( ς β A ) a 2 C ( [ 0 , ˜ ] , H δ , r ) , E β ( ς β A ) a 2 C ( [ 0 , ˜ ] , H α , r ) , 0 ς E β ( ( ς s ) β A ) a 2 d s C ( [ 0 , ˜ , H δ , r ) , 0 ς E β ( ( ς s ) β A ) a 2 d s C ( ( 0 , ˜ , H α , r ) .
Hence, from the previous Lemma 6, this yields that
0 ς E β ( ( ς s ) β A ) a 2 d s X , ς β ( α δ ) 0 ς E β ( ( ς s ) β A ) a 2 d s C ( [ 0 , ˜ ] , H α , r ) .
| | 0 ς E β ( ( ς s ) β A ) a 2 d s | | X 0 ς | | E β ( ( ς s ) β A ) a 2 d s | | X , 0 ς C 1 | a 2 | H δ , r d s , = C 1 | a ς 2 | H δ , r .
With the help of Equation (5), the inequality gives us
| | 0 ς E β ( ( ς s ) β A ) a 2 d s + η ( ς ) | | X [ ˜ ] | | 0 ς E β ( ( ς s ) β A ) a 2 d s | | X [ ˜ ] + | | η ( ς ) | | X [ ˜ ] C 1 | a ς 2 | H δ , r + B 1 M 1 4 L ,
which gives the result that F has a unique fixed point due to Lemma 4. □

6. Existence Locally in J r

We are using the Iteration method in this part intended for thinking of the local existence of a mild solution to problem (2) in J r . Let α = 1 + δ 2 .
Theorem 3. 
Suppose 1 < r < , 0 < δ < 1 and (e) holds. Let us take value in H δ , r
a ς 2 H δ , r w i t h n 2 r 1 2 < δ .
So, the mild solution of (2) is a unique answer for a ς 2 H δ , r in J r . Additionally, ς β ( α δ ) A α u ( ς ) is bounded as ς 0 . Moreover, u and A α u are both continuous function in [ 0 , ˜ ] and ( 0 , ˜ ] , respectively.
Proof. 
Step I
κ ˜ ( ς ) : = sup s ( 0 , ς ] s β ( α δ ) | A α u ( s ) | r
also
ψ ˜ ( ς ) : = ϕ ( u , u ) ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s .
The instant results from Step II in Theorem 1 reveal that A is continuous functions ψ ˜ ( ς ) , and A α ψ ˜ ( ς ) exists in [ 0 , ˜ ] and ( 0 , ˜ ] , respectively, with value
| A α ψ ˜ ( ς ) | r | 0 ς A α ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s | r C 1 0 ς A α ( ς s ) β ( 1 α ) 1 | F ( u ( s ) , u ( s ) ) | r d s M C 1 0 ς A α ( ς s ) β ( 1 α ) 1 s 2 β ( α δ ) sup s ( 0 , ς ] s 2 β ( α δ ) | u ( s ) | r | u ( s ) | r d s = M C 1 0 ς ( ς s ) β ( 1 α ) 1 s 2 β ( α δ ) d s sup s ( 0 , ς ] A α s 2 β ( α δ ) | ( u ( s ) | 2 | A α ψ ˜ ( ς ) | r M C 1 ς β ( α δ ) B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ 2 ( ς ) .
We also take into account the integral η ( ς ) . Given that ( e ) is true,
| P ( σ B 0 2 ρ u ( s ) ) | r M ( ς ) s β ( 1 δ )
the above inequality with a continuous function M ( ς ) holds. Theorem 1’s third stage reveals that A α η ( ς ) is continuous in ( 0 , ˜ ] , with results
| A α η ( ς ) | r = | 0 ς ( ς s ) β 1 A α E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) | H α , r d s C 1 0 ς ( ς s ) β ( 1 α ) 1 | A α P ( σ B 0 2 ρ u ( s ) ) | r d s C 1 M ( ς ) 0 ς ( ς s ) β ( 1 α ) 1 s β ( 1 δ ) d s = ς β ( α δ ) C 1 M ( ς ) B ( β ( 1 δ ) , 1 β ( 1 δ ) ) | A α η ( ς ) | r ς β ( α δ ) C 1 M ( ς ) B ( β ( 1 δ ) , 1 β ( 1 δ ) ) .
For | P ( σ B 0 2 ρ u ( ς ) ) | r = 0 ( ς β ( 1 δ ) ) as ς tends to zero, and M ( ς ) = 0 . The above Equation (10) concludes that | A α η ( ς ) | r = 0 ( ς β ( α δ ) ) , ς approaches to 0. We establish the continuity of η in J r . In actuality, we take 0 < ς 0 < ς < ˜ , resulting in
| η ( ς ) η ( ς 0 ) | r C 3 ς 0 ς ( ς s ) β 1 | P ( σ B 0 2 u ( s ) ) ρ | r d s + C 3 0 ς 0 ( ( ς 0 s ) β 1 ( ς s ) ) β 1 | P ( σ B 0 2 u ( s ) ) ρ | r d s + C 3 0 ς 0 ε ( ς 0 s ) β 1 | | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) | | | P ( σ B 0 2 u ( s ) ) ρ | r d s + 2 C 3 ς 0 ε ς 0 ( ς 0 s ) β 1 | P ( σ B 0 2 u ( s ) ) ρ | r d s C 3 M ( ς ) ς 0 ς ( ς s ) β 1 s β ( 1 δ ) d s + C 3 M ( ς ) 0 ς 0 ( ( ς s ) β 1 ( ς 0 s ) β 1 ) s β ( 1 δ ) d s + C 3 M ( ς ) 0 ς 0 ε ( ς 0 s ) β 1 s β ( 1 δ ) d s sup s [ 0 , ς ε ] | | E β , β ( ( ς s ) β A ) E β , β ( ( ς 0 s ) β A ) | | + 2 C 3 M ( ς ) ς 0 ε ς 0 ( ς 0 s ) β 1 s β ( 1 δ ) d s 0 , a s ς ς 0 .
Additionally, we think about the function 0 ς E β ( ( ς s ) β A ) a 2 d s . It is clear from the Lemma 6 that
| 0 ς A α E β ( ( ς s ) β A ) a 2 d s | r 0 ς | A α E β ( ( ς s ) β A ) a 2 | r d s C 1 ς β ( α δ ) A δ | 0 ς a 2 d s | r = C 1 ς β ( α δ ) | a 2 ς | H δ , r lim ς 0 ς β ( α δ ) | 0 ς A α E β ( ( ς s ) β A ) a 2 d s | r , = lim ς 0 ς β ( α δ ) | 0 ς E β ( ( ς s ) β A ) a 2 d s | H α , r = 0 .
Step II
Using the successive approximation method, we now arrive at the following solution:
u 0 ( ς ) = 0 ς E β ( ( ς s ) β A ) a 2 d s + η ( ς ) , u n + 1 ( ς ) = u 0 ( ς ) + ϕ ( u n , u n ) ( ς ) , n = 0 , 1 , 2 , . . . .
Using the results above, we have that κ ˜ n ( ς ) : = sup s ( 0 , ς ] s β ( α δ ) | A α u n ( s ) | r are continuous and increasing in [ 0 , ˜ ] with value κ ˜ n ( 0 ) = 0 . However, given (9) and (10), and that the inequality is satisfied by κ ˜ n ( ς ) ,
κ ˜ n + 1 ( ς ) κ ˜ 0 ( ς ) + M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ n 2 ( ς ) .
For κ ˜ 0 ( 0 ) = 0 , select ˜ > 0 ,
4 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ 0 ( ˜ ) < 1 .
Therefore, it is assured that the sequence κ ˜ n ( ˜ ) is constrained by the fundamental consideration (12), i.e.,
κ ˜ n ( ˜ ) ρ ˜ ( ˜ ) , n = 0 , 1 , 2 , . . . .
There are
ρ ˜ ( ς ) = 1 1 4 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ 0 ( ς ) 2 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) .
It is true that κ ˜ n ( ς ) ρ ˜ ( ς ) holds for any value of ς ( 0 , ˜ ] , similarly to how we say ρ ˜ ( ς ) 2 κ ˜ 0 ( ς ) .
Let us think about the equality
z n + 1 ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) [ F ( u n + 1 ( s ) , u n + 1 ( s ) F ( u n ( s ) , u n ( s ) ) ] d s ,
for ς ( 0 , ˜ ] and z n = u n + 1 u n , n = 0 , 1 , 2 , . . . by writing
Z n ( ς ) : = sup s ( 0 , ς ] s β ( α δ ) | A α z n ( s ) | r .
In light of (6), there are
| F ( u n + 1 ( s ) , u n + 1 ( s ) F ( u n ( s ) , u n ( s ) ) | r M | u n + 1 | H α , r + | u n | H α , r | u n + 1 u n | = M | u n + 1 | H α , r + | u n | H α , r A α z n sup s ( 0 , ς ] s β ( α δ ) s β ( α δ ) = M | A α u n + 1 | r + | A α u n | r z n sup s ( 0 , ς ] s β ( α δ ) s β ( α δ ) = M ( sup s ( 0 , ς ] s β ( α δ ) | A α u n + 1 | r + sup s ( 0 , ς ] s β ( α δ ) | A α u n | r ) z n s β ( α δ ) M κ ˜ n + 1 + κ ˜ n z n s β ( α δ ) sup s ( 0 , ς ] s β ( α δ ) s β ( α δ ) A α M κ ˜ n + 1 + κ ˜ n sup s ( 0 , ς ] s β ( α δ ) | A α z n ( s ) | r s 2 β ( α δ ) | F ( u n + 1 ( s ) , u n + 1 ( s ) F ( u n ( s ) , u n ( s ) ) | r M κ ˜ n + 1 + κ ˜ n Z n ( s ) s 2 β ( α δ ) .
It is implied by step II in Theorem 1 that
ς β ( α δ ) | A α z n + 1 ( ς ) | r 2 M C 1 B ( β ( 1 α ) , 1 β ( 1 δ ) ) ρ ˜ ( ˜ ) Z n ( ς ) .
Such inequality results in
Z n + 1 ( ˜ ) 2 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) ρ ˜ ( ˜ ) Z n ( ˜ ) 4 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ 0 ( ˜ ) Z n ( ˜ ) .
As per (13) with (14), it is indeed obvious that
lim n 0 Z n + 1 ( ˜ ) Z n ( ˜ ) < 4 M C 1 B ( β ( 1 α ) , 1 2 β ( α δ ) ) κ ˜ 0 ( ˜ ) < 1 .
Thus, convergence of the series Σ n = 0 Z n ( ˜ ) implies the uniform convergence of series Σ n = 0 ς β ( α δ ) A α z n ( ς ) for ς ( 0 , ˜ ] ; the uniform convergence of sequence { ς β ( α δ ) A α u n ( ς ) } holds in ( 0 , ˜ ] . So,
lim n u n ( ς ) = u ( ς ) D ( A α )
and
lim n ς β ( α δ ) A α u n ( ς ) = ς β ( α δ ) A α u ( ς ) u n i f o r m l y .
From the boundedness theorem, “A function f continuous on a bounded and closed interval is necessarily a bounded function”. So, both A α and A α are bound and closed, respectively. The function κ ˜ ( ς ) = sup s ( 0 , ς ] s β ( α δ ) | A α u ( s ) | r satisfies
κ ˜ ( ς ) ρ ˜ ( ς ) 2 κ ˜ 0 ( ς ) , ς ( 0 , ς ] ,
as well as
ς n = sup s ( 0 , ˜ ] s 2 β ( α δ ) | F ( u n ( s ) , u n ( s ) F ( u ( s ) , u ( s ) ) | r M ( κ ˜ n ( ˜ ) + κ ˜ ( ˜ ) ) sup s ( 0 , ˜ ] s β ( α δ ) | A α ( u n ( s ) u ( s ) ) | r 0 , a s n .
To finalize this step, it is important to verify that u is a suitable solution to the issue (2) in the range of [ 0 , ˜ ] .
| ϕ ( u n , u n ) ( ς ) ϕ ( u , u ) ( ς ) | r 0 ς ( ς s ) β 1 ς n s 2 β ( α δ ) d s = ς β δ ς n 0 , ( n ) .
To put it another way, we get ϕ ( u n , u n ) ( ς ) ϕ ( u , u ) ( ς ) . Limits are taken on both sides of (10), and we conclude
u ( ς ) = u 0 ( ς ) + ϕ ( u , u ) ( ς ) .
Suppose u ( 0 ) = a ς 2 ; what we learn is that (16) is true for both ς [ 0 , ˜ ] and u C ( [ 0 , ˜ ] , J r ) . Additionally, the continuity of A α u ( ς ) in ( 0 , ˜ ] is derived from the uniform convergence of ς β ( α δ ) A α u n ( ς ) to ς β ( α δ ) A α u ( ς ) . We conclude that | A α u ( ς ) | r = 0 ( ς β ( α δ ) ) is clear from (15) and κ ˜ 0 ( ς ) = 0 .
Step III
We illustrate the difference between “mild solutions.” Assume that the issue has mild solutions in u and ν (2). Suppose z = u ν ; consider the inequality.
z ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) [ F ( u ( s ) , u ( s ) F ( ν ( s ) , ν ( s ) ) ] d s .
Expressing the solution
κ ˜ ( ς ) : = m a x { sup s ( 0 , ς ] s β ( α δ ) | A α ( u ( s ) ) | r , sup s ( 0 , ς ] s β ( α δ ) | A α ( ν ( s ) ) | r } .
From (5) and Lemma 6, we have the inequality
| A α z ( ς ) | r M C 1 κ ˜ ( ς ) 0 ς ( ς s ) β ( 1 α ) 1 s β ( α δ ) | A α z ( s ) | r d s .
The Gronwall inequality shows that ς ( 0 , ˜ ] , A α z ( ς ) = 0 , which indicates it for ς [ 0 , ˜ ] , z ( ς ) = u ( ς ) ν ( ς ) 0 . As a result, the mild solution is unique. □

7. Regularity

In this Section (2), the regularity of “the solution u that” resolves the issue is studied. In this essay, we’ll assume: ( e 1 ) with an exponent ϑ ( 0 , β ( 1 α ) ) , P σ B 0 2 u ( ς ) ρ is Hölder-continuous, that is,
| f ( x ) f ( y ) | | | x y | | β , | P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( s ) ρ ) | r L | ς s | ϑ , 0 < ς , s ˜ .
Definition 5. 
An expression u : [ 0 , ˜ ] J r , . If u C ( [ 0 , ˜ ] , J r ) with C D ς β u ( ς ) C ( ( 0 , ˜ ] , J r ) , which accepts values in D ( A ) and solves (2) for every ς ( 0 , ˜ ] , then J r is referred to as the classical solution of (2).
Lemma 7. 
Assume ( e 1 ) be satisfied. If
η 1 ( ς ) : = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) d s , ς ( 0 , ˜ ] ,
so, η 1 ( ς ) D ( A ) and, also, A η 1 ( ς ) C ϑ ( [ 0 , ˜ ] , J r ) .
Proof. 
To be fixed ς ( 0 , ˜ ] . Let us think about
( ς s ) β 1 | A E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r .
Lemma 6 with ( e 1 ) , give us
( ς s ) β 1 | A E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r C 1 ( ς s ) β 1 ( ς s ) β | P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r = C 1 ( ς s ) 1 | P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r C 1 L ( ς s ) ϑ 1 L 1 ( [ 0 , ˜ ] , J r ) .
Then,
| A η 1 ( ς ) | r 0 ς ( ς s ) β 1 | A E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r d s C 1 L 0 ς ( ς s ) ϑ 1 d s C 1 L ς ϑ ϑ < .
Since A is closed, we can write η 1 ( ς ) D ( A ) .
It is necessary to demonstrate that A η 1 ( ς ) is Hölder-continuous because
d d ς ς β 1 E β , β ( μ ς β ) = ς β 2 E β , β 1 ( μ ς β ) .
Then,
d d ς ς β 1 A E β , β ( ς β A ) = 1 2 π i Γ θ ς β 2 E β , β 1 ( μ ς β ) A ( μ I + A ) 1 d μ = 1 2 π i Γ θ ς β 2 E β , β 1 ( μ ς β ) d μ 1 2 π i Γ θ ς β 2 μ E β , β 1 ( μ ς β ) ( μ I + A ) 1 d μ .
Put
μ ς β = ψ ς β d μ = d ψ d μ = d ψ ς β = 1 2 π i Γ θ ς β 2 E β , β 1 ( ψ ) 1 ς β d ψ 1 2 π i Γ θ ς β 2 E β , β 1 ( ψ ) ψ ς β ψ ς β I + A 1 1 ς β d ψ ,
in the light of
| | ( μ I + A ) 1 | | C | μ | ,
we derive that
| | d d ς ς β 1 A E β , β ( ς β A ) | | C β ς 2 , 0 < ς < ˜ .
From MVT, we derive that ∀ 0 < s < ς ˜ , and we obtain
| | ς β 1 A E β , β ( ς β A ) s β 1 A E β , β ( s β A ) | | = | | s ς d d τ τ β 1 A E β , β ( τ β A ) d τ | | s ς | | d d τ τ β 1 A E β , β ( τ β A ) | | d τ C β s ς τ 2 d τ = C β ( s 1 ς 1 ) .
Take h > 0 such that 0 < ς < ς + h ˜ , then
A η 1 ( ς + h ) A η 1 ( ς ) = 0 ς ( ( ς + h s ) β 1 A E β , β ( ( ς + h s ) β A ) ( ς s ) β 1 A E β , β ( ( ς s ) β A ) ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) d s + 0 ς ( ς + h s ) β 1 A E β , β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) d s + ς ς + h ( ς + h s ) β 1 A E β , β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) d s : = I 1 ( ς ) + I 2 ( ς ) + I 3 ( ς ) .
The three major terms are discussed here one by one. We have Equation (18) with ( e 1 ) for I 1 ( ς )
| I 1 ( ς ) | r 0 ς | | ( ς + h s ) β 1 A E β , β ( ( ς + h s ) β A ) ( ς s ) β 1 A E β , β ( ( ς s ) β A ) | | | P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) | r d s C β L h 0 ς ( ς + h s ) 1 ( ς s ) ϑ 1 d s C β L h 0 ς ( h + s ) 1 ( ς s ) ϑ 1 d s C β L 0 ς h s + h s ϑ 1 d s + C β L h h s s + h s ϑ 1 d s C β L h ϑ .
For solving I 2 ( ς ) , Lemma 6 and ( e 1 ) are used here, so we have
| I 2 ( ς ) | r 0 ς ( ς + h s ) β 1 | A E β , β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) | r d s C 1 0 ς ( ς + h s ) 1 | P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) | r d s C 1 L h ϑ 0 ς ( ς + h s ) 1 d s = C 1 L [ l n ( h ) l n ( ς + h ) ] h ϑ .
Moreover, for solving I 3 ( ς ) , Lemma 6 and ( e 1 ) are used here, so we have
| I 3 ( ς ) | r ς ς + h ( ς + h s ) β 1 | A E β , β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) | r d s C 1 ς ς + h ( ς + h s ) 1 | P ( σ B 0 2 u ( ς ) ρ ) P ( σ B 0 2 u ( ς + h ) ρ ) ) | r d s
C 1 L ς ς + h ( ς + h s ) ϑ 1 d s = C 1 L h ϑ ϑ .
By combining all of the above (20), (21) and (22), we deduce that Hölder’s continuity of A η 1 ( ς ) exists. □
Theorem 4. 
Assume Theorem 3’s assumptions are satisfied. The mild solution of Equation (2) is the classical one applicable to any a ς 2 D ( A ) if ( e 1 ) is true.
Proof. 
From Lemma 2(ii), it is guaranteed that function u ( ς ) = 0 ς E β ( ( ς s ) β A ) a 2 d ς ( ς > 0 ) is a classical solution for the given problem for value a t 2 D ( A ) ,
C D ς β u = A u , ς > 0 , u ( 0 ) = a x 2 .
Step I
By verifying that
η ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) d s
is the classical solution to this following problem
C D ς β u = A u + P ( σ B 0 2 u ( ς ) ρ ) , ς > 0 , u ( 0 ) = 0 .
Thus, it follows from Theorem 3 that η C ( [ 0 , ˜ ] , J r ) . By rewriting η ( ς ) = η 1 ( ς ) + η 2 ( ς ) , then
η 1 ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) P ( σ B 0 2 u ( ς ) ρ ) d s , η 2 ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( ς ) ρ ) d s .
From Lemma 7, we can write η 1 ( ς ) D ( A ) to demonstrate the similar results for η 2 ( ς ) . By Lemma 2(iii), we find that A η 2 ( ς ) = P ( σ B 0 2 u ( ς ) ρ ) E β ( ς β A ) P ( σ B 0 2 u ( ς ) ρ ) . Since ( e 1 ) exists, | A η 2 ( ς ) | r ( 1 + C 1 ) | P ( σ B 0 2 u ( ς ) ρ ) | r ; thus,
η 2 ( ς ) D ( A ) f o r ς ( 0 , ˜ ] a n d η 2 ( ς ) C ϑ ( ( 0 , ˜ ] , J r ) .
Additionally, we verify for C D ς β η C ( ( 0 , ˜ ] , J r ) . In light of the Lemma 2(iv) with the condition η ( 0 ) = 0 , there are
C D ς β η ( ς ) = d d ς ( I ς 1 β η ( ς ) ) = d d ς E β ( ς β A ) P ( σ B 0 2 u ρ ) .
The continuous differentiability of E β ( ς β A ) P ( σ B 0 2 u ρ ) in J r remains to be proven. If we assume that 0 < h ˜ ς , we can derive:
1 h E β ( ( ς + h ) β A ) P ( σ B 0 2 u ρ ) E β ( ς β A ) P ( σ B 0 2 u ρ ) = 0 ς 1 h E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) E β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) d s + 1 h t ς + h E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) d s .
Consider
0 ς 1 h | E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) E β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) | r d s C 1 1 h 0 ς | E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) | r d s + C 1 1 h 0 ς | E β ( ( ς s ) ) β A ) P ( σ B 0 2 u ( s ) ρ ) | r d s C 1 M ( ς ) 1 h 0 ς ( ς + h s ) β s β ( 1 δ ) d s + C 1 M ( ς ) 1 h 0 ς ( ς s ) β s β ( 1 δ ) d s C 1 M ( ς ) 1 h ( ς + h ) 1 β + ς 1 β B ( 1 β , 1 β ( 1 δ ) ) .
Using Lebesgue’s LDC theorem, we arrive at the following conclusion:
lim h 0 0 ς 1 h E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) E β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) d s = 0 ς ( ς s ) β 1 A E β , β ( ( ς s ) β A ) P ( σ B 0 2 u ( s ) ρ ) d s = A η ( ς ) .
In contrast with
1 h ς ς + h E β ( ( ς + h s ) β A ) P ( σ B 0 2 u ( s ) ρ ) d s = 1 h 0 h E β ( s β A ) P ( σ B 0 2 ρ u ( ς + h s ) d s = 1 h 0 h E β ( s β A ) P ( σ B 0 2 ρ u ( ς + h s ) ) P ( σ B 0 2 ρ u ( ς s ) ) d s + 1 h 0 h E β ( s β A ) P ( σ B 0 2 ρ u ( ς s ) ) P ( σ B 0 2 ρ u ( ς ) ) d s + 1 h 0 h E β ( s β A ) P ( σ B 0 2 ρ u ( ς ) ) d s .
Lemmas 1 and 6 with property ( e 1 ) give us
| 1 h 0 ς E β ( s β A ) P ( σ B 0 2 ρ u ( ς + h s ) ) P ( σ B 0 2 ρ u ( ς s ) ) d s | r C 1 L h ϑ , | 1 h 0 t E β ( s β A ) P ( σ B 0 2 ρ u ( ς s ) ) P ( σ B 0 2 ρ u ( ς ) ) d s | r C 1 L h ϑ ϑ + 1 .
Consequently, Lemma 2(i) offers
lim h 0 1 h 0 h E β ( s ) β A P ( σ B 0 2 ρ u ( s ) ) d s = P ( σ B 0 2 ρ u ( ς ) ) .
Hence,
lim h 0 1 h ς ς + h E β ( ς + h s ) β A P ( σ B 0 2 ρ u ( s ) ) d s = P ( σ B 0 2 ρ u ( ς ) ) .
Our conclusion is that E β ( ς β A ) P ( σ B 0 2 ρ u ) is differentiable at ς + and
d d ς E β ( ς β A ) P ( σ B 0 2 ρ u ) + = A η ( ς ) + P ( σ B 0 2 ρ u ( ς ) ) .
Similarly, E β ( ς β A ) P ( σ B 0 2 ρ u ) is differential at ς and
d d ς E β ( ς β A ) P ( σ B 0 2 ρ u ) = A η ( ς ) + P ( σ B 0 2 ρ u ( ς ) ) .
We show that A η = A η 1 + A η 2 C ( ( 0 , ˜ ] , J r ) . It is obviously clear that given function
η 2 ( ς ) = P ( σ B 0 2 ρ u ( ς ) ) E β ( ς β A ) P ( σ B 0 2 ρ u ( ς ) ) .
In consideration of Lemma 1, it is continuous because of Lemma 2(iii). Additionally, Lemma 7 tells us that A η 1 ( ς ) is also continuous. Accordingly, C D ς β η C ( ( 0 , ˜ ] , J r ) .
Step II
Suppose u is a mild solution of (2). To demonstrate F ( u , u ) C ϑ ( ( 0 , ˜ ] , J r ) from (5), we must verify that A α u is Hölder-continuous in J r . Apply h > 0 in a way that 0 < ς < ς + h . Indicate ϕ ˜ ( ς ) : = E β ( ς β A ) a 2 through Lemma 2(iv) and (6), then
| A α ϕ ˜ ( ς + h ) A α ϕ ˜ ( ς ) | r = | ς ς + h s β 1 A α E β , β ( s β A ) a 2 d s | r ς ς + h s β 1 | A α δ E β , β ( s β A ) A δ a 2 | r d s C 1 ς ς + h s β 1 s β ( α δ ) | A δ a 2 | r d s = C 1 ς ς + h s β ( δ α ) + β 1 | A δ a 2 | r d s = C 1 ς ς + h s β ( 1 + δ α ) 1 d s | A δ a 2 | r = C 1 | a | H δ , r β ( 1 + δ α ) ( ς + h ) β ( 1 + δ α ) ς β ( 1 + δ α ) C 1 | a 2 | H δ , r β ( 1 + δ α ) h β ( 1 + δ α ) .
Thus, A α ϕ ˜ C ϑ ( ( 0 , ˜ ] , J r ) . Apply h in a way that ε ς < ς + h ˜ for all small ε > 0 ,
| A α η ( ς + h ) A α η ( ς ) | r | ς ς + h ( ς + h s ) β 1 A α E β , β ( ( ς + h s ) β A ) P ( σ B 0 2 ρ u ( s ) ) d s | r + | 0 ς A α ( ς + h s ) β 1 E β , β ( ( ς + h s ) β A ) ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) d s | r = η 1 ( ς ) + η 2 ( ς ) .
Lemma 6 and ( e ) are applied, and the outcome is
η 1 ( ς ) C 1 ς ς + h ( ς + h s ) β ( 1 α ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s C 1 M ( ς ) ς ς + h ( ς + h s ) β ( 1 α ) 1 s β ( 1 δ ) d s M ( ς ) C 1 β ( 1 α ) h β ( 1 α ) ς β ( 1 δ ) M ( ς ) C 1 β ( 1 α ) h β ( 1 α ) ε β ( 1 δ ) .
Finding the inequality, we may calculate η 2 ( ς ) ,
d d ς ς β 1 A α E β , β ( ς β A ) = 1 2 π i Γ μ α ς β 2 E β , β 1 ( μ ς β ) ( μ I + A ) 1 d μ = 1 2 π i Γ ψ ς β α ς β 2 E β , β 1 ( ψ ) ψ ς β I + A 1 1 ς β d ψ ,
The above equation gives the results that
| | d d ς ( ς β 1 A α E β , β ( ς β A ) ) | | C β ς β ( 1 α ) 2 .
The mean value theorem yields
| | ς β 1 A α E β , β ( ς β A ) s β 1 A α E β , β ( s β A ) | | s ς | | d d τ ( τ β 1 A α E β , β ( τ β A ) ) | | d τ C β s ς τ β ( 1 α ) 2 d τ = C β s β ( 1 α ) 1 ς β ( 1 α ) 1 .
Thus,
η 2 ( ς ) 0 ς | A α ( ς + h s ) β 1 E β , β ( ( ς + h s ) β A ) ( ς s ) β 1 E β , β ( ( ς s ) β A P ( σ B 0 2 ρ u ( s ) ) | r d s 0 ς ( ς s ) β ( 1 α ) 1 ( ς + h s ) β ( 1 α ) 1 | P ( σ B 0 2 ρ u ( s ) ) | r d s C β M ( ς ) 0 ς ( ς s ) β ( 1 α ) 1 s β ( 1 δ ) d s 0 ς + h ( ς + h s ) β ( 1 α ) 1 s β ( 1 δ ) d s + C β M ( ς ) ς ς + h ( ς + h s ) β ( 1 α ) 1 s β ( 1 δ ) d s C β M ( ς ) ς β ( δ α ) ( ς + h ) β ( δ α ) B ( β ( 1 α ) , 1 β ( 1 δ ) ) + C β M ( ς ) h β ( 1 α ) ς β ( 1 δ ) C β M ( ς ) h β ( α δ ) [ ε ( ε + h ) ] β ( δ α ) + C β M ( ς ) h β ( 1 α ) ε β ( 1 δ ) .
This guarantees that A α η C ϑ ( [ ε , ˜ ] , J r ) . Owing to random ε , A α η C ϑ ( ( 0 , ˜ ] , J r ) ,
ψ ˜ ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) F ( u ( s ) , u ( s ) ) d s .
Here, we know that | F ( u ( s ) , u ( s ) ) | r M κ ˜ 2 ( ς ) s 2 β ( α δ ) , in which the supplied function’s continuity exists but is also bounded in ( 0 , ˜ ] , and κ ˜ ( ς ) : = sup s ( 0 , ς ] s β ( α δ ) | u ( s ) | H α , r . We are able to provide the Hölder continuity of A α ψ ˜ in the same fashion in C ϑ ( ( 0 , ˜ ] , J r ) . Therefore,
A α u ( ς ) = A α φ ˜ ( ς ) + A α η ( ς ) + A α ψ ˜ ( ς ) C ϑ ( ( 0 , ˜ ] , J r ) .
Seeing as F ( u , u ) C ϑ ( ( 0 , ˜ ] , J r ) is demonstrated from the previous step II, this gives the results that C D ς β ψ C ϑ ( ( 0 , ˜ ] , J r ) , A ψ ˜ C ( ( 0 , ˜ ] , J r ) and C D ς β ψ ˜ = A ψ ˜ + F ( u , u ) . Similarly, we obtained that C D ς β u C ( ( 0 , ˜ ] , J r ) ,   A u C ( ( 0 , ˜ ] , J r ) and C D ς β u = A u + F ( u , u ) + P ( σ B 0 2 ρ u ) . Consequently, the conclusion is that u is a classic solution. □
Theorem 5. 
Suppose ( e 1 ) is true. If u is showing a classical solution of (2), then A u C ν ( ( 0 , ˜ ] , J r ) and, also, C D ς β u C ν ( ( 0 , ˜ ] , J r ) .
Proof. 
We may put u ( ς ) = φ ˜ ( ς ) + η ( ς ) + ψ ˜ ( ς ) if u exhibits the classical solution of (2). The evidence is sufficient to demonstrate A φ ˜ C β ( 1 δ ) ( ( 0 , ˜ ] , J r ) . It is necessary to demonstrate that A ϕ ˜ C β ( 1 δ ) ( [ ε , ˜ ] , J r ) is true for any ε > 0 . In fact, by choosing h such that ε ς < ς + h ˜ , using Lemma 2(iii),
| A ϕ ˜ ( ς + h ) A ϕ ˜ ( ς ) | r = | ς ς + h s β 1 A 2 E β , β ( s β A ) a 2 d s | r C 1 ς ς + h s β ( 1 δ ) 1 d s | a 2 | H δ , r = C 1 | a 2 | H δ , r β ς β ( 1 δ ) ( ς + h ) β ( 1 δ ) C 1 | a 2 | H δ , r β h β ( 1 δ ) [ ε ( ε + h ) ] β ( 1 δ ) .
As from Lemma 7, we can write η ( ς ) as
η ( ς ) = η 1 ( ς ) + η 2 ( ς ) = 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) P ( σ B 0 2 ρ u ( s ) ) P ( σ B 0 2 ρ u ( ς ) ) d s + 0 ς ( ς s ) β 1 E β , β ( ( ς s ) β A ) ( P ( σ B 0 2 ρ u ( s ) ) d s ,
in the domain of ς ( 0 , ˜ ] from Lemma 7 and (24) it follows that A η 1 ( ς ) C ν ( [ 0 , ˜ ] , J r ) and A η 2 ( ς ) C ϕ ˜ ( ( 0 , ˜ ] , J r ) , accordingly. □

8. Conclusions

The Helmholtz-Leray projection is used in this work to show the existence and uniqueness of fractional-order Navier-Stokes equations for the Cauchy problem solution. In the interim, we provide a workable local solution in S . To model phenomenon diffusion in fractal media, Navier-Stokes equations (NSEs) with time-fractional derivatives of order γ ( 0 , 1 ) are utilised. We use S to show that these equations have regular classical solutions. Additional study may build on the idea presented in this article by incorporating validity and generalising other activities. Numerous studies are being done in this interesting field, which could result in a variety of ideas and uses.

Author Contributions

Conceptualization, R.S.; Methodology, A.U.K.N.; Software, M.A.; Formal analysis, K.A.; Investigation, A.U.K.N.; Resources, M.A.; Data curation, M.A.; Writing—original draft, R.S.; Writing—review & editing, R.S.; Supervision, A.U.K.N.; Project administration, K.A.; Funding acquisition, K.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported through the Annual Funding track by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Project No. GRANT2328].

Data Availability Statement

No new data were created this study.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Gorenflo, R.; Mainardi, F. Fractional calculus. In Fractals and Fractional Calculus in Continuum Mechanics; Springer: Vienna, Austria, 1997; pp. 223–276. [Google Scholar]
  2. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  3. Machado, T.; Kiryakova, V.; Mainardi, F. A poster about the recent history of fractional calculus. Fract. Calc. Appl. Anal. 2010, 13, 329–334. [Google Scholar]
  4. El-Shahed, M.; Salem, A. On the generalized Navier-Stokes equations. Appl. Math. Comput. 2004, 156, 287–293. [Google Scholar] [CrossRef]
  5. Doering, C.R.; Gibbon, J.D. Applied Analysis of the Navier-Stokes Equations (No. 12); Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  6. Gresho, P.M.; Sani, R.L. On pressure boundary conditions for the incompressible Navier-Stokes equations. Int. J. Numer. Methods Fluids 1987, 7, 1111–1145. [Google Scholar] [CrossRef]
  7. Lemarié-Rieusset, P.G. Recent Developments in the Navier-Stokes Problem; CRC: Boca Raton, FL, USA, 2002. [Google Scholar]
  8. Von Wahl, W. The Equations of Navier-Stokes and Abstract Parabolic Equations; Vieweg+Teubner Verlag Wiesbaden: Braunschweig, Germany, 1985. [Google Scholar]
  9. Iwabuchi, T.; Takada, R. Global solutions for the Navier-Stokes equations in the rotational framework. Math. Ann. 2013, 357, 727–741. [Google Scholar] [CrossRef]
  10. Davidson, P.A. An Introduction to magnetohydrodynamics. Am. J. Phys. 2002, 70, 781. [Google Scholar] [CrossRef] [Green Version]
  11. Goedbloed, J.H.; Goedbloed, J.P.; Poedts, S. Principles of Magnetohydrodynamics: With Applications to Laboratory and Astrophysical Plasmas; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  12. Desjardins, B.; Le Bris, C. Remarks on a nonhomogeneous model of magnetohydrodynamics. Differ. Integral Equ. 1998, 11, 377–394. [Google Scholar] [CrossRef]
  13. Qayyum, M.; Abbas, T.; Afzal, S.; Saeed, S.T.; Akgül, A.; Inc, M.; Mahmoud, K.H.; Alsubaie, A.S. Heat Transfer Analysis of Unsteady MHD Carreau Fluid Flow over a Stretching/Shrinking Sheet. Coatings 2022, 12, 1661. [Google Scholar] [CrossRef]
  14. Rehman, A.U.; Riaz, M.B.; Saeed, S.T.; Jarad, F.; Jasim, H.N.; Enver, A. An exact and comparative analysis of MHD free convection flow of water-based nanoparticles via CF derivative. Math. Probl. Eng. 2022, 2022, 9977188. [Google Scholar] [CrossRef]
  15. Saeed, S.T.; Riaz, M.B.; Awrejcewicz, J.; Ahmad, H. Exact Symmetric Solutions of MHD Casson Fluid Using Chemically Reactive Flow with Generalized Boundary Conditions. Energies 2021, 14, 6243. [Google Scholar] [CrossRef]
  16. Niazi, A.U.K.; He, J.; Shafqat, R.; Ahmed, B. Existence, Uniqueness, and E q–Ulam-Type Stability of Fuzzy Fractional Differential Equation. Fractal Fract. 2021, 5, 66. [Google Scholar] [CrossRef]
  17. Shafqat, R.; Niazi, A.U.K.; Jeelani, M.B.; Alharthi, N.H. Existence and Uniqueness of Mild Solution Where α∈(1,2) for Fuzzy Fractional Evolution Equations with Uncertainty. Fractal Fract. 2022, 6, 65. [Google Scholar] [CrossRef]
  18. Alnahdi, A.S.; Shafqat, R.; Niazi, A.U.K.; Jeelani, M.B. Pattern formation induced by fuzzy fractional-order model of COVID-19. Axioms 2022, 11, 313. [Google Scholar] [CrossRef]
  19. Khan, A.; Shafqat, R.; Niazi, A.U.K. Existence Results of Fuzzy Delay Impulsive Fractional Differential Equation by Fixed Point Theory Approach. J. Funct. Spaces 2022, 2022, 4123949. [Google Scholar] [CrossRef]
  20. Abuasbeh, K.; Shafqat, R.; Niazi, A.U.K.; Awadalla, M. Local and global existence and uniqueness of solution for time-fractional fuzzy Navier-Stokes equations. Fractal Fract. 2022, 6, 330. [Google Scholar] [CrossRef]
  21. Abbas, A.; Shafqat, R.; Jeelani, M.B.; Alharthi, N.H. Significance of Chemical Reaction and Lorentz Force on Third-Grade Fluid Flow and Heat Transfer with Darcy–Forchheimer Law over an Inclined Exponentially Stretching Sheet Embedded in a Porous Medium. Symmetry 2022, 14, 779. [Google Scholar] [CrossRef]
  22. Phong, T.T.; Long, L.D. Well-posed results for nonlocal fractional parabolic equation involving Caputo-Fabrizio operator. J. Math. Comput. Sci. 2022, 26, 357–367. [Google Scholar] [CrossRef]
  23. Shafqat, R.; Niazi, A.U.K.; Yavuz, M.; Jeelani, M.B.; Saleem, K. Mild Solution for the Time-Fractional Navier-Stokes Equation Incorporating MHD Effects. Fractal Fract. 2022, 6, 580. [Google Scholar] [CrossRef]
  24. Mehmood, Y.; Shafqat, R.; Sarris, I.E.; Bilal, M.; Sajid, T.; Akhtar, T. Numerical Investigation of MWCNT and SWCNT Fluid Flow along with the Activation Energy Effects over Quartic Auto Catalytic Endothermic and Exothermic Chemical Reactions. Mathematics 2022, 10, 4636. [Google Scholar] [CrossRef]
  25. Boulares, H.; Benchaabane, A.; Pakkaranang, N.; Shafqat, R.; Panyanak, B. Qualitative Properties of Positive Solutions of a Kind for Fractional Pantograph Problems using Technique Fixed Point Theory. Fractal Fract. 2022, 6, 593. [Google Scholar] [CrossRef]
  26. Zhou, Y.; Vijayakumar, V.; Murugesu, R. Controllability for fractional evolution inclusions without compactness. Evol. Equ. Control. Theory 2015, 4, 507–524. [Google Scholar] [CrossRef]
  27. Wang, R.N.; Chen, D.H.; Xiao, T.J. Abstract fractional Cauchy problems with almost sectorial operators. J. Differ. Equ. 2012, 252, 202–235. [Google Scholar] [CrossRef] [Green Version]
  28. Podlubny, I.; Chechkin, A.; Skovranek, T.; Chen, Y.; Jara, B.M.V. Matrix approach to discrete fractional calculus II: Partial fractional differential equations. J. Comput. Phys. 2009, 228, 3137–3153. [Google Scholar] [CrossRef]
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

Abuasbeh, K.; Shafqat, R.; Niazi, A.U.K.; Awadalla, M. Mild Solutions for the Time-Fractional Navier-Stokes Equations with MHD Effects. Symmetry 2023, 15, 280. https://doi.org/10.3390/sym15020280

AMA Style

Abuasbeh K, Shafqat R, Niazi AUK, Awadalla M. Mild Solutions for the Time-Fractional Navier-Stokes Equations with MHD Effects. Symmetry. 2023; 15(2):280. https://doi.org/10.3390/sym15020280

Chicago/Turabian Style

Abuasbeh, Kinda, Ramsha Shafqat, Azmat Ullah Khan Niazi, and Muath Awadalla. 2023. "Mild Solutions for the Time-Fractional Navier-Stokes Equations with MHD Effects" Symmetry 15, no. 2: 280. https://doi.org/10.3390/sym15020280

APA Style

Abuasbeh, K., Shafqat, R., Niazi, A. U. K., & Awadalla, M. (2023). Mild Solutions for the Time-Fractional Navier-Stokes Equations with MHD Effects. Symmetry, 15(2), 280. https://doi.org/10.3390/sym15020280

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop