Analytical and numerical analyses of a viscous strain gradient problem involving type Ⅱ thermoelasticity

1.   Introduction

The theme of fractional calculus (FC) has appeared as a broad and interesting research point due to its broad applications in science and engineering. FC is now greatly evolved and embraces a wide scope of interesting findings. To obtain detailed information on applications and recent results about this topic, we refer to ^[1,2,3,4,5] and the references therein.

Some researchers in the field of FC have realized that innovation for new FDs with many non-singular or singular kernels is very necessary to address the need for more realistic modeling problems in different fields of engineering and science. For instance, we refer to works of Caputo and Fabrizio ^[6], Losada and Nieto ^[7] and Atangana-Baleanu ^[8]. The class of FDs and fractional integrals (FIs) concerning functions is a considerable branch of FC. This class of operators with analytic kernels is a new evolution proposed in ^[1,9,10]. Every one of these operators is appropriate broader to cover various kinds of FC and catch diversified behaviors in fractional models. Joining the previous ideas yields another, significantly wide, which is a class of function-dependent-kernel fractional derivatives. This covers both of the two preceding aforesaid classes see ^[11,12].

In another context, the Langevin equations (LEs) were formulated by Paul Langevin in 1908 to describe the development of physical phenomena in fluctuating environments ^[13]. After that, diverse generalizations of the Langevin equation have been deliberated by many scholars we mention here some works ^[14,15,16]. Recently, many researchers have investigated sufficient conditions of the qualitative properties of solutions for the nonlinear fractional LEs involving various types of fractional derivatives (FDs) and by using different types of methods such as standard fixed point theorems (FPTs), Leray-Schauder theory, variational methods, etc., e.g. ^{[17,18,19,20,21,22,23,24]}. Some recent results on the qualitative properties of solutions for fractional LEs with the generalized Caputo FDs can be found in ^{[25,26,27,28,29,30]}, e.g., Ahmad et al. ^[25] established the existence results for a nonlinear LE involving a generalized Liouville-Caputo-type

Seemab et al. ^[26] investigated the existence and UHR stability results for a nonlinear implicit LE involving a $\vartheta$ -Caputo FD

The Hyers-Ulam (U-H) stability and existence for various types of generalized FDEs are established in the papers ^{[31,32,33,34,35,36,37,38,39,40,41,42,43]}.

Motivated by the above works and inspired by novel developments in $\vartheta$ -FC, in the reported research, we investigate the existence, uniqueness, and U-H-type stability of the solutions for the nonlinear fractional Langevin differential equation (for short FLDE) described by

where $\mathcal{D}_{a^{+}}^{\varepsilon; \vartheta }$ denote the $\vartheta$-Caputo FD of order $\varepsilon \in \{\alpha_{1}, \alpha_{2} \}$ such that $\alpha_{1} \in (2, 3]$, $\alpha_{2} \in (1, 2]$, $\xi \in (a, b)$, $\delta > 0$, $\mathcal{I}_{a^{+}}^{\alpha_{1};\vartheta }$ is the $\vartheta$-fractional integral of the Riemann-Liouville (RL) type, $\mathcal{F}:\mathrm{J}\times \mathbb{R}\rightarrow \mathbb{R}$ is a continuous function, and $\lambda, \kappa$ $\in \mathbb{R},$ $\xi \in (a, b).$

The considered problem in this work is more general, in other words, when we take certain values of function $\vartheta$, the problem (1.3) is reduced to many problems in the frame of classical fractional operators. Also, the gained results here are novel contributes and an extension of the evolution of FDEs that involving a generalized Caputo operator, especially, the study of stability analysis of Ulam-Hyers type of fractional Langevin equations is a qualitative addition to this work. Besides, analysis of the results was restricted to a minimum of assumptions.

Here is a brief outline of the paper. Section 2 provides the definitions and preliminary results required to prove our main findings. In Section 3, we establish the existence, uniqueness, and stability in the sense of Ulam for the system (1.3). In Section 5, we give some related examples to light the gained results.

2.   Preliminaries and lemmas

We start this part by giving some basic definitions and results required for fractional analysis.

Consider the space of real and continuous functions $\mathcal{U} = C(\mathrm{J}, \mathbb{R})$ space with the norm

Let $\vartheta \in \mathcal{C}^{1} = C^{1}(\mathrm{J}, \mathbb{R})$ be an increasing differentiable function such that $\vartheta ^{\prime }(\varsigma)\neq 0$, for all $\varsigma \in \mathrm{J.}$

Now, we start by defining $\vartheta$-fractionals operators as follows:

Definition 2.1. ^[1] The $\vartheta$-RL fractional integral of order $\alpha_{1} > 0$ for an integrable function $\omega \colon \mathrm{J} \longrightarrow \mathbb{R}$ is given by

Definition 2.2. ^[1] Let $\alpha_{1} \in(n-1, n)$, $n\in \mathbb{N}$, $\omega: \mathrm{J} \rightarrow \mathbb{R}$ is an integrable function, and $\vartheta \in C^{n}(\mathrm{J}, \mathbb{R})$, the $\vartheta$-RL FD of a function $\omega$ of order $\alpha_{1}$ is given by

where $n = [\alpha_{1}]+1$ and $D_{\varsigma} = \frac{d}{dt}$.

Definition 2.3. ^[9] For $\alpha_{1} \in(n-1, n)$, and $\omega, \vartheta \in C^{n}(\mathrm{J}, \mathbb{R})$, the $\vartheta$-Caputo FD of a function $\omega$ of order $\alpha_{1}$ is given by

where $n = [\alpha_{1}]+1$ for $\alpha_{1} \notin \mathbb{N}$, $n = \alpha_{1}$ for $\alpha_{1} \in \mathbb{N}$, and $\omega _{\vartheta }^{[n]}(\varsigma) = \left(\frac{ D_{\varsigma}}{\vartheta ^{\prime }(\varsigma)}\right) ^{n}\omega (\varsigma).$

From the above definition, we can express $\vartheta$-Caputo FD by formula

Also, the $\vartheta$-Caputo FD of order $\alpha_{1}$ of $\omega$ is defined as

For more details see [9,Theorem 3].

Lemma 2.4. ^[1] For $\alpha_{1}, \alpha_{2} > 0,$ and $\omega \in C(\mathrm{J}, \mathbb{R})$, we have

Lemma 2.5. ^[44] Let $\alpha_{1} > 0.$

If $\omega \in C(\mathrm{J}, \mathbb{R}),$ then

and if $\omega \in C^{n-1}(\mathrm{J}, \mathbb{R}),$ then

for all $\varsigma \in\mathrm{J}.$ Moreover, if $m \in \mathrm{N}$ be an integer and $\omega \in \mathrm{C}^{n+m}(\mathrm{J}, \mathbb{R})$ a function. Then, the following holds:

Observe that from Eq (2.3), if $\omega_{\vartheta}^{[k]}(a) = 0,$ for all $k = n, n+1, \ldots, n+m-1$ we can get the following relation

Lemma 2.6. ^[1,9] For $\varsigma > a, \alpha_{1} \geq 0,$ $\alpha_{2} > 0,$ we have

● $\mathcal{I}_{a^{+}}^{\alpha_{1};\vartheta }(\vartheta (\varsigma)-\vartheta (a))^{\alpha_{2} -1} = \frac{\Gamma (\alpha_{2})}{\Gamma (\alpha_{2} +\alpha_{1})}(\vartheta (\varsigma)-\vartheta (a))^{\alpha_{2} +\alpha_{1} -1},$

● $\mathit{\mbox{}}^{c}\mathcal{D}_{a^{+}}^{\alpha_{1};\vartheta }(\vartheta (\varsigma)-\vartheta (a))^{\alpha_{2} -1} = \frac{\Gamma (\alpha_{2})}{\Gamma (\alpha_{2} -\alpha_{1})}(\vartheta (\varsigma)-\vartheta (a))^{\alpha_{2} -\alpha_{1} -1},$

● $\mathit{\mbox{}}^{c}\mathcal{D}_{a^{+}}^{\alpha_{1};\vartheta }(\vartheta (\varsigma)-\vartheta (a))^{k} = 0, \forall k\in \{0, \dots, n-1\}, n\in \mathbb{N}.$

Theorem 2.7. (Banach's FPT ^[45]). Let $(R, d)$ be a nonempty complete metric space with a contraction mapping $\mathcal{G}: R \rightarrow R$ i.e., $\mathrm{d}(\mathcal{G} z, \mathcal{G} \varkappa) \leq \mathrm{L} \mathrm{\; d}(z, \varkappa)$ for all $z, \varkappa \in R,$ where $\mathrm{L} \in(0, 1)$ is a constant. Then $\mathcal{G}$ possesses a unique fixed point.

Theorem 2.8. (Kransnoselskii's FPT ^[46]). Let $\mathrm{E}$ be a Banach space. Let $S$ is a nonempty convex, closed and bounded subset of $\mathrm{E}$ and let $A_{1}, A_{2}$ be mapping from $S$ to $\mathrm{E}$ such that:

(i) $A_{1} z+A_{2} \varkappa \in S$ whenever $z, \varkappa \in S$

(ii) $A_{1}$ is continuous and compact;

(iii) $A_{2}$ is a contraction.

Then there exists $z \in S$ such that $z = A_{1} z+A_{2} z$.

3.   Main results

This portion interests in the existence, uniqueness, and Ulam stability of solutions to the suggested problem (1.3).

The next auxiliary lemma, which attentions the linear term of a problem (1.3), plays a central role in the afterward findings.

Lemma 3.1. Let $\alpha_{1}\in(2, 3]$, $\alpha_{2}\in(1, 2]$. Then the linear BVP

has a unique solution defined by

where

and

with

Proof. Applying the RL operator $\mathcal{I}_{a^{+}}^{\alpha_{2};\vartheta }$ to (3.1) it follows from Lemma 2.5 that

Again, we apply the RL operator $\mathcal{I}_{a^{+}}^{\alpha_{1};\vartheta }$ and use the results of Lemma 2.5 to get

where $c_{0}, c_{1}, c_{2}, c_{3}, c_{4}\in \mathbb{R}.$ By utilizing the boundary conditions in (3.1) and (3.7), we obtain $c_{2} = 0, c_{3} = 0, c_{4} = 0$.

Hence,

Now, by using the conditions $\mathcal{D}_{a^{+}}^{\alpha_{1};\vartheta }z(a) = \mathcal{I}_{a^{+}}^{\delta; \vartheta }z(\gamma)$ and $\mathcal{D} _{a^{+}}^{\alpha_{1};\vartheta }z(b)+\kappa z(b) = 0$, we get

Substituting $c_{0}$ and $c_{1}$ in (3.8), we finish with (3.2).

The reverse direction can be shown easily with the help of results in Lemmas 2.5 and 2.6, i.e. Eq (3.2) solves problem (3.1). This ends the proof.

Now, we shall need to the following lemma:

Lemma 3.2. The functions $\mu$ and $\nu$ are continuous functions on $\mathrm{J}$ and satisfy the following properties:

(1) $\mu^{\ast} = \max _{0 \leq \varsigma \leq b}|\mu(\varsigma)|$,

(2) $\nu^{\ast} = \max _{0 < \varsigma < b}|\nu(\varsigma)|$,

where $\mu$ and $\nu$ are defined by Lemma 3.1.

Here, we give the following hypotheses:

(H1) The function $\mathcal{F}:\mathrm{J} \times \mathbb{R} \rightarrow \mathbb{R}$ is continuous.

(H2) There exists a constant $\mathcal{L} > 0$ such that

(H3) There exist positive functions $h(\varsigma)\in C\left(\mathrm{J}, \mathbb{R }^{+}\right)$ with bounds $\Vert h \Vert$ such that

For simplicity, we denote

As a result of Lemma 3.1, we have the subsequent lemma:

Lemma 3.3. Suppose that $\mathcal{F}: \mathrm{J} \times \mathbb{R} \rightarrow \mathbb{R}$ is continuous. A function $z(\varsigma)$ solves (1.3) if and only if it is a fixed-point of the operator $\mathcal{G}: \mathcal{U} \rightarrow \mathcal{U}$ defined by

Now, we are willing to give our first result which based on Theorem 2.7.

Theorem 3.4. Suppose that (H1) and (H2) hold. If $\Delta < 1$, where $\Delta$ is given by (3.11), then there exists a unique solution for (1.3) on the interval $\mathrm{J}$.

Proof. Thanks to Lemma 3.1, we consider the operator $\mathcal{G}: \mathcal{U} \rightarrow \mathcal{U}$ defined by (3.13). Thus, $\mathcal{G}$ is well defined as $\mathcal{F}$ is a continuous. Then the fixed point of $\mathcal{G}$ coincides with the solution of FLDE $(1.3).$ Next, the Theorem 2.7 will be used to prove that $\mathcal{G}$ has a fixed point. For this end, we show that $\mathcal{G}$ is a contraction.

Let $\mathcal{B}_{\mathcal{R}} = \{z \in \mathcal{U}:\|z \| \leq \mathcal{R}\}$, where $\mathcal{R} > \frac{\mathcal{M} \Lambda_{1}}{1-\mathcal{L}\Lambda_{1}- \Lambda_{2}}$. Since

we obtain

which implies that $\|\mathcal{G} z\| \leq \mathcal{R}$, i.e.,

Now, let $z, \varkappa \in \mathcal{U}$. Then, for every $\varsigma \in \mathrm{J},$ using (H2), we can get

Also note that

Using the above arguments, we get

As $\Delta < 1,$ we derive that $\mathcal{G}$ is a contraction. Hence, by Theorem 2.7, $\mathcal{G}$ has a unique fixed point which is a unique solution of FLDE (1.3). This ends the proof.

Now, we apply the Theorem 2.8 to obtain the existence result.

Theorem 3.5. Let us assume (H1)–(H3) hold. Then FLDE (1.3) has at least one solution on $\mathrm{J}$ if $\Lambda_{3} < 1$, where it is supposed that

Proof. By the assumption (H3), we can fix

where $B_{\rho } = \{z\in \mathcal{U}:\Vert z\Vert\leq \rho \}.$ Let us split the operator $\mathcal{G}:\mathcal{U}\rightarrow \mathcal{U}$ defined by (3.13) as $\mathcal{G} = \mathcal{G}_{1}+\mathcal{G}_{2},$ where $\mathcal{G}_{1}$ and $\mathcal{G}_{2}$ are given by

and

The proof will be split into numerous steps:

Step 1: $\mathcal{G}_{1}(z)+\mathcal{G}_{2}(z)\in B_{\rho }$.

Hence

which shows that $\mathcal{G}_{1}z+\mathcal{G}_{2}z_{1}\in B_{\rho }.$

Step 2: $\mathcal{G}_{2}$ is a contraction map on $B_{\rho }$.

Due to the contractility of $\mathcal{G}$ as in Theorem 3.4, then $\mathcal{G}_{2}$ is a contraction map too.

Step 3: $\mathcal{G}_{1}$ is completely continuous on $B_{\rho }$.

From the continuity of $\mathcal{F}\left(\cdot, z(\cdot)\right)$, it follows that $\mathcal{G}_{1}$ is continuous.

Since

we get $\left\Vert \mathcal{G}_{1}z\right\Vert\leq \mathcal{N}$ which emphasize that $\mathcal{G}_{1}$ uniformly bounded on $B_{\rho }$.

Finally, we prove the compactness of $\mathcal{G}_{1}$.

For $z\in B_{\rho }\ $and $\varsigma \in \mathrm{J}$, we can estimate the operator derivative as follows:

where we used the fact

Hence, for each $\varsigma _{1}, \varsigma _{2}\in \mathrm{J}$ with $a < \varsigma _{1} < \varsigma _{2} < b$ and for $z\in B_{\rho }$, we get

which as $(\varsigma _{2}-\varsigma _{1})$ tends to zero independent of $z$. So, $\mathcal{G}_{1}$ is equicontinuous. In light of the foregoing arguments along with Arzela–Ascoli theorem, we derive that $\mathcal{G}_{1}$ is compact on $B_{\rho }$. Thus, the hypotheses of Theorem 2.8 holds, so there exists at least one solution of (1.3) on $\mathrm{J}$.

4.   Ulam-Hyers stability analysis for the $\vartheta$-Caputo FLDE (1.3)

In the current section, we are interested in studying Ulam-Hyers (U-H) and the generalized Ulam-Hyers stability types of the problem (1.3).

Let $\varepsilon > 0$. We consider the next inequality:

Definition 4.1. FLDE (1.3) is stable in the frame of U-H type if there exists $c_{\mathcal{F}}\in \mathbb{R}^{+}$ such that for every $\varepsilon > 0$ and for each solution $\tilde{z} \in \mathcal{U}$ of the inequality (4.1) there exists a solution $z \in \mathcal{U}$ of (1.3) with

Definition 4.2. FLDE (1.3) has the generalized U-H stability if there exists $C_{\mathcal{F}}:$ $\mathrm{C}\left(\mathbb{R}_{+}, \mathbb{R}_{+}\right)$ along with $C_{\mathcal{F}}(0) = 0$ such that for every $\varepsilon > 0$ and for each solution $\tilde{z} \in \mathcal{U}$ of the inequality (4.1), a solution $z \in C(\mathrm{J}, \mathbb{R})$ of (1.3) exists uniquely for which

Remark 4.3. A function $\tilde{z} \in \mathcal{U}$ is a solution of the inequality (4.1) if and only if there exists a function $\varrho \in \mathcal{U}$ (which depends on solution $\tilde{z}$) such that

1.$|\varrho(\varsigma)| \leq \varepsilon, \varsigma \in\mathrm{J}$.

2.$\mathcal{D}_{a^{+}}^{\alpha_{2}; \vartheta}\left(\mathcal{D}_{a^{+}}^{a, \vartheta}-\lambda\right) \tilde{z}(\varsigma) = \mathcal{F}(\varsigma, \tilde{z} (\varsigma))+\varrho(\varsigma), \quad \varsigma \in\mathrm{J}.$

Theorem 4.4. Let $\Delta < 1$, (H1) and (H2) hold. Then the FLDE (1.3) is U-H stable on $\mathrm{J}$ and consequently generalized U-H stable.

Proof. For $\varepsilon > 0$ and $\tilde{z} \in C(\mathrm{J}, \mathbb{R})$ be a function which fulfills the inequality (4.1). Let $z \in \mathcal{U}$ the unique solution of

By Lemma 3.1, we have

Since we have assumed that $\widetilde{z}$ is a solution of (4.1), hence by Remark 4.3

Again by Lemma 3.1, we have

On the other hand, for any $\varsigma\in\mathrm{J}$

Using part (i) of Remark 4.3 and (H2), we get

where $\Delta$ is defined by $(3.11).$ In consequence, it follows that

If we let $c_{\mathcal{F}} = \left(\frac{(\vartheta(b)-\vartheta(a))^{\alpha_{1}+\alpha_{2}}}{ \Gamma(\alpha_{1}+\alpha_{2}+1)}\frac{\left[\nu^{\ast}(\kappa-\lambda)+1\right]}{ (1-\Delta)}+\frac{\nu^{\ast}}{(1-\Delta)}\frac{(\vartheta(b)-\vartheta(a))^{ \alpha_{2}}}{\Gamma(\alpha_{2}+1)}\right),$ then, the U-H stability condition is satisfied. More generally, for $C_{\mathcal{F}}(\varepsilon) = \left(\frac{ (\vartheta(b)-\vartheta(a))^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\frac{\left[ \nu^{\ast}|\kappa-\lambda|+1\right]}{(1-\Delta)}+\frac{\nu^{\ast}}{ (1-\Delta)}\frac{(\vartheta(b)-\vartheta(a))^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)}\right) \varepsilon$; $C_{\mathcal{F}}(0) = 0$ the generalized U-H stability condition is also fulfilled.

5.   Examples

This section is intended to illustrate the reported results with relevant examples.

Example 5.1. We formulate the system of FLDE in the frame of Caputo type:

In this case we take

Obviously, the hypothesis (H1) of the Theorem 3.4 is fulfilled. On the opposite hand, for each $\varsigma \in[0, 1], z, \varkappa \in \mathbb{R}$ we have

Hence, $\left(H_{2}\right)$ holds with $\mathcal{L} = 0.1.$ Thus, we find that $\Delta = 0.7635 < 1.$ Since all the assumptions in Theorem 3.4 hold, the FLDE (5.1) has a unique solution on $[0, 1].$ Moreover, Theorem 4.4 ensures that the FLDE (1.3) is U-H stable and generalized U-H stable.

Example 5.2. We formulate the system of FLDE in the frame of Hadamard type:

Here

Obviously, the assumption (H1) of the Theorem 3.4 holds. On the other hand, for any $\varsigma \in[1, e], z, \varkappa \in \mathbb{R}$ we get

Consequently, $\left(H_{2}\right)$ holds with $\mathcal{L} = 0.1.$ Besides, by computation directly we find that $\Delta = 0.8731 < 1$.

To illustrate Theorem 3.5, it is clear that the function $f$ satisfies (H1) and (H3) with $\Vert h\Vert = 0.1$. In addition, $\Lambda_3 \approx 0.7823 < 1$. It follows from theorem 3.5 that the FLDE (5.2) has a unique solution on $[1, e]$.

Example 5.3. We formulate the system of FLDE in the frame of $\vartheta$-Caputo type:

Take

For any $\varsigma \in[0, 1], z, \varkappa \in \mathbb{R}$ we obtain

Hence, (H2) holds with $\mathcal{L} = 0.125.$ Moreover, by computation directly we find that $\Delta = 0.7133 < 1$. It follows from Theorem 3.4 that FLDE (5.4) has a unique solution on $[0, 1].$

To illustrate Theorem 3.5, it is clear that the function $\mathcal{F}(\varsigma, z)$ given by (5.5) satisfies the hypotheses (H1)–(H3) with $\Vert h\Vert = 0.125$. and $\Lambda_3 \approx 0.3247 < 1$. It follows from Theorem 3.5 that the FLDE (5.2) has a unique solution on $[0, 1]$.

6.   Conclusions

In this reported article, we have considered a class of nonlinear Langevin equations involving two different fractional orders in the frame of Caputo function-dependent-kernel fractional derivatives with antiperiodic boundary conditions. The existence and uniqueness results are established for the suggested problem. Our perspective is based on properties of $\vartheta$-Caputo's derivatives and applying of Krasnoselskii's and Banach's fixed point theorems. Moreover, we discuss the Ulam-Hyers stability criteria for the at-hand problem. Some related examples illustrating the effectiveness of the theoretical results are presented. The results obtained are recent and provide extensions to some known results in the literature. Furthermore, they cover many fractional Langevin equations that contain classical fractional operators.

Acknowledgments

This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, Saudi Arabia through the Fast-track Research Funding Program.

Conflict of interest

The authors declare that there is no conflict of interests regarding the publication of this paper.