Laboratory study of polymer injection into heavy oil unconventional reservoirs to enhance oil recovery and determination of optimal injection concentration

1.   Introduction

Fractional differential equations (FDEs) are an effective mathematical tool to model and analyze many real life problems; it has been used by researchers and scientists to get better results than the integer order differential equations. Fractional order differential equations offer a superior framework for capturing the intricate dynamics of real-world phenomena compared to their integer-order counterparts. This superiority stems from the unique ability of fractional integrals and derivatives to account for the inherent hereditary and memory characteristics present in diverse processes and materials. By harnessing these fractional operators, models can more accurately depict the nuanced behaviors observed in nature, thereby enhancing our understanding and predictive capabilities across a wide range of disciplines and applications. Many fractional derivatives, including Caputo derivative, Atangana-Baleanu derivative, Coimbra derivative, and Riemann-Liouville (R-L) derivative, are frequently used to examine FDEs and fractional order stochastic differential equations (SDEs). In a similar vein, the Hilfer fractional derivative (HFD), which was just recently used to do so, was developed by Hilfer ^[21], which is a generalized version of R-L and Caputo derivatives. In actuality, fractional derivative and integrals indicate greater accuracy than integral models and also depict broader physical applications in seepage, flow in porous media, nanotechnology, fluid dynamics and traffic models ^{[6,7,12,22,24,28,31]}.

SDEs are the natural extension of deterministic systems. SDEs with impulses arise from many mathematical models of physical phenomena in different scientific fields for example, technology, physics, biology, economics, etc. They are important from the viewpoint of applications since they incorporate randomness into the mathematical description of the phenomena and provide a more accurate description of it. Certainly, in various fields like economics, bioengineering, chemistry, medicine, and biology, we often encounter situations where things change suddenly at specific points in time ^[5,23,29,32]. These abrupt changes can be explained by what we call "impulsive effects." These impulsive effects are like sudden pushes that happen at certain moments and have a big impact on the system, being studied. These pushes play a crucial role in understanding and modeling how things change in the aforementioned diverse fields. For the mathematical models of such phenomena, finding their solution is a challenging task. As a result, Boundani et al. ^[8,9] presented some specific conditions that help us to determine whether certain mathematical equations, involving randomness, can have solutions. These equations involve functional differential equations and a type of random behavior, called fractional Brownian motion.

In ^[20], Hernandez and O'Regan introduced non-instantaneous (NI) impulses. Many researchers have utilized these impulses and studied the corresponding dynamical systems ^{[3,4,27,36,39]}. To better understand NI impulses, we can think about human blood sugar levels. When they have too much or too little, glucose they get insulin medication through the bloodstream, in fact it doesn't work instantly but takes some time to be absorbed ^[37]. This gradual effect is like NI, where the impact lasts for a while in many real situations. Sudden changes don't explain things well. For instance, in the treatment of diseases with medication, we need to describe how things change over time more smoothly. That's where NI impulsive differential equations come in handy. They help us to model these gradual changes, like how drugs affect the body in pharmacotherapy.

Among the qualitative behaviors of different physical systems, different types of stabilities are the essential ones. One of these types of stabilities is Ulam-Hyers (UH) and Ulam-Hyers-Rassias (UHR) stability ^{[11,16,26,30,33]}.

In the literature, most of the results related to fractional stochastic differential equations (FSDEs) are given over infinite dimensional spaces ^{[1,2,13,18,19,25,32,35,38]}, and very few have looked at similar results in finite spaces. There is no prior work on the specific topic of non-integers order impulses in finite-dimensional equations to the problem of fractional neutral SDEs including both noises. In FSDEs incorporating both retarded and advanced arguments, a significant characteristic emerges: The rate of change of the system at the present moment is influenced not only by its past history but also by its anticipated future states. This feature underscores the intricate interplay between past, present, and future dynamics, where the system's behavior is shaped by a combination of its memory effects and anticipatory responses.

In the current manuscript, we investigate the following $\psi$-Hilfer fractional stochastic equation (HFSE):

where $D^{\gamma, \beta}_{\psi(\varepsilon)}$ is the $\psi$– Hilfer FD of order $0 < \gamma < 1$ and of type $0 < \beta\leq1$. Let $\mathscr{J} : = [0, b]$, $b > 0$. The state vector $\Gamma \in R^{n}$, $A\in R^{n\times n}$ and nonlinear functions h : $\mathscr{J} \times R^{n} \longrightarrow R^{n}$, $\Delta:\mathscr{J} \times R^{n} \longrightarrow R^{n}$, $g:\mathscr{J} \times R^{n} \longrightarrow R^{n\times n}$, $\lambda:\mathscr{J} \times R^{n} \longrightarrow R^{n\times n}$, and $h_{k}:\mathscr{J} \times R^{n} \longrightarrow R^{n}$ are measurable and bounded functions. Also, $\Gamma_{0}$ is $\mathscr{F}_{0}$ measurable $R^{n}$-valued stochastic variable and $B$ is an $n$-dimensional Wiener process.

Based on the value of $H$, the following kinds of the fractional Brownian motion (fBm) process exist:

$(1)$ if $H = \frac{1}{2}$, then the process is a Brownian motion or a Wiener process exists;

$(2)$ if $H > \frac{1}{2}$, then the increment of the process is positively correlated;

$(3)$ if $H < \frac{1}{2}$, then the increment of the process is negatively correlated.

The contributions of this paper are described as below:

● Nonlinear $\psi$-HFSE is considered in $R^{n}$.

● The existence and uniqueness results are established by using the standard Banach contraction principle.

● The weaker sufficient conditions are derived by using the generalized Schaefer FPT for the system with measure of non-compactness (MNC).

● UHR stability results are derived for $\psi$-HFSE with NI impulse.

● An example is provided for the theoretical results.

This paper is structured as follows: In Section 2, we present a number of lemmas and some fundamental definitions for fractional calculus. Section 3 derives the solution representation of $\psi$-Hilfer fractional SDEs with NI impulses. To demonstrate the key results, the generalized Schaefer's and contraction mapping principles are used in Section 4. In Section 5, UHR stability of a $\psi$-HFSE is discussed. Example is demonstrated for the validity of theoretical results in Section 6.

Notations:

● $(\Omega, \mathscr{F}, \mathscr{P})$ represents the complete probability space along a probability measure $\mathscr{P}$ on $\Omega$.

● $B(\varepsilon)$ and $B^{H}(\varepsilon)$ denote, respectively, the n-dimensional Brownian motion and fBm with Hurst index $\frac{1}{2} < H < 1$.

● $\{\mathscr{F_{\varepsilon}|\varepsilon \in \mathscr{J}}\}$ represents the filtration generated by $\{B(\varepsilon):0\leq s \leq \varepsilon \}$.

● $L_{2}(\Omega, \mathscr{F_{\varepsilon}}, \mathscr{P}, R^{n})$ : = $L_{2}(\Omega, R^{n})$ is the space of all $\mathscr{F_{\varepsilon}}$-measurable square integrable random variables with values in $R^{n}$.

Let $\mathscr{J}_{k}$ = $(\varepsilon_{k}, \varepsilon_{k+1}]$, $k = 1, 2, \dots, m$ be such that impulse times satisfy 0 = $\varepsilon_{0}$ = $s_{0} < \varepsilon_{1} < s_{1} < \varepsilon_{2} < \dots < \varepsilon_{m} \leq s_{m} < \varepsilon_{m+1}$ = b. Let $C(\mathscr{J}, L_{2}(\Omega, R^{n}))$ denote $C_{n}(\mathscr{J})$ and be the Banach space of all continuous maps from $\mathscr{J}$ into $L_{2} (\Omega, R^{n})$ of $\mathscr{F_{\varepsilon}}$-adapted square integrable functions $\Gamma(\varepsilon)$ and for its norm $\sup E||\Gamma(\varepsilon)||^{2} < \infty$.

Define the space

and there exist $\Gamma(\varepsilon^{-}_k)$ and $\Gamma(\varepsilon^{+}_k)$ with $\Gamma(\varepsilon_k) = \Gamma(\varepsilon^{-}_k)$, $k = 1, \dots, m$, endowed with the norm

Clearly, $\mathbb{Y}$ is a Banach space.

2.   Preliminaries

Lemma 2.1. ^[39] Let $p\geq 2$ and $\mathfrak{f} \in L_p (\mathscr{J}, R^{n\times n})$ such that $E |\int_0^b\mathfrak{f(s)}dB(s)|^p < \infty$, then

Lemma 2.2. ^[10] Let $\varphi : \mathscr{J}\rightarrow L^0_2$ satisfy $\int_0^b ||\varphi(s)||^2_{L^0_2}ds < \infty$, then we get

Definition 2.1. ^[17] The generalized $\psi$-Hilfer FD of order $0 < \gamma < 1$ and of type $0 \leq \beta \leq 1$ is represented by

Definition 2.2. ^[10] Let $z$ be a bounded linear operator. The two parameter Mittage-Leffler (M-L) function is defined by

One of the interesting properties of the M-L function, related with their Laplace integral, is given by

That is

Lemma 2.3. ^[17] For $\gamma\in (n-1, n]$, $\beta \in [0, 1]$, the following Laplace formula for the $\psi$-Hilfer derivative is valid:

Corollary 2.1. ^[17] If $f$ is a function whose classical Laplace transform is $F(s)$, then the generalized Laplace transform of the function $f \circ \psi = f(\psi(t))$ is also $F(s)$:

Example 2.1. ^[17]

Example 2.2. Assume that $Re(\mu) > 0$ and $|\frac{A}{s^{\mu}}| < 1$. If $M^{\gamma}_{\mu, v}$ denotes the Prabhakar function, then we have

3.   Solution representation

Consider the linear deterministic system, which is represented in the following form:

By the Laplace transformation, we get

where $I$ is the identity matrix.

Substituting the $L_\psi T$ of the M-L function, one can obtain that

Therefore, the solution of (1.1) is given as follows:

Definition 3.1. ^[40] The set S is called a quasi-equicontinuous in $\mathscr{T}$ if, for $\varepsilon$ $>$ 0, there exists a $\delta$ $>$ 0, such that if $\Gamma \in S, K\in \mathscr{N}, \tau_1, \tau_2 \in \mathscr{T_{\mathscr{K}}}\bigcap\mathscr{T}$, and $|\tau_2-\tau_1| < \delta$, then $|\Gamma(\tau_2)-\Gamma(\tau_1)| < \varepsilon$.

Lemma 3.1. ^[34] The set $\text{S}\subset PC(\mathscr{J}, R^{n})$ is relatively compact if

(i) $\textrm{S}$ is uniformly bounded, i.e., $||\Gamma||_{PC} \textrm{k}$ for each $\Gamma \in \textrm{S}$ and some $\textrm{k}> 0$;

(ii) $\textrm{S}$ is quasi-equicontinous in $\mathscr{T}$.

Definition 3.2. ^[14] An operator $T$ : $\mathcal{Z} \rightarrow \mathcal{Z}$ is said to be $\chi$-condensing of any bounded set $\mathscr{B}$ of $\text{Z}$ with $\chi$-condensing($\mathscr{B}$) $>$ 0, $\chi(T(\mathscr{B})) < \chi(\mathscr{B})$, where

is Kuratowskii MNC of a bounded set $\mathscr{B}$ of Z.

Now, we state the generalized Schaefer's type FPT with $\chi$-condensing operators.

Theorem 3.1. ^[15] Let $\text{T}$ : $\text{Z} \rightarrow \text{Z}$ be an operator and $\text{Z}$ be a separable Banach space satisfying

($\textrm{A1}$) $\textrm{T}$ is $\chi$-condensing and continuous.

($\textrm{A2}$) The set $\textrm{S}$ = $\{\Gamma\in \text{Z} : \Gamma = \delta\text{T}(\Gamma)$ for some $0 < \delta < 1\}$ is bounded, then $\textrm{T}$ has a fixed point.

For convenience, define the following:

To derive the existence result, we imposed the following assumptions:

(H1) The functions $\Delta, h, g, \lambda$, and $h_{k}$ $k = 1, 2, \dots, m$ are Lipschitz continuous.

$(1)$ $E||h(\varepsilon, u_1) - h(\varepsilon, u_2)||^{2}$ $\leq M_{h}||u_1 - u_2||^{2}_\mathbb{Y}$;

$(2)$ $E||\Delta(\varepsilon, u_1, v_1) - \Delta(\varepsilon, u_2, v_2)||^{2}$ $\leq M_{f_1}||u_1-u_2||^{2}_\mathbb{Y}+M_{f_2}||v_1-v_2||^{2}_\mathbb{Y}$;

$(3)$ $E||g((\varepsilon, u_1, v_1)) - g((\varepsilon, u_2, v_2))||^{2}$ $\leq M_{g_1}||u_1 - u_2||^{2}_\mathbb{Y}+ M_{g_2}||v_1 - v_2||^{2}_\mathbb{Y}$;

$(4)$ $E||\lambda(\varepsilon, u_1, v_1) - \lambda(\varepsilon, u_2, v_2)||^{2}$ $\leq M_{\lambda_1}||u_1 - u_2||^{2}_\mathbb{Y}+ M_{\lambda_2}||v_1 - v_2||^{2}_\mathbb{Y}$;

$(5)$ $E|h_{k}(\varepsilon, u_1) - h_{k}(\varepsilon, v_1)||^{2}$ $\leq M_{h{k}}||u_1- v_1||^{2}_\mathbb{Y}$,

and $h_{k} \in C((\varepsilon_{k}, s_{k}], L_{2}(\Omega, R^{n}))$, where $M_{h}, M_{f}, M_{g}, M_{\lambda}$, and $M_{h{k}}$ are positive constants, $u_1, u_2, v_1, v_2\in R^{n}$, and $\varepsilon \in \mathscr{T_{k}}$.

(H2) There exist $l_{h}$, $m_{h}, n_{h} \in \mathbb{Y}$ with $l^{*}_{h} = \sup_{\xi \in \mathcal{J}} l_{f}(\varepsilon)$, $m^*_{h} = \sup_{\xi \in \mathcal{J}} m_{h}(\varepsilon)$, and $n^*_{h} = \sup_{\xi \in \mathcal{J}}$ $n_{h}(\varepsilon)$ such that

(H3) There exist $l_{f}$, $m_{f}, n_{f} \in \mathbb{Y}$ with $l^{*}_{f} = \sup_{\xi \in \mathcal{J}} l_{f}(\varepsilon)$, $m^*_{f} = \sup_{\xi \in \mathcal{J}} m_{f}(\varepsilon)$, and $n^*_{f} = \sup_{\xi \in \mathcal{J}} n_{f}(\varepsilon)$ such that

(H4) There exist $l_{g}$, $m_{g} \in \mathbb{Y}$ with $l^{*}_{g} = \sup_{\xi \in \mathcal{J}} l_{g}(\varepsilon)$, $m^*_{g} = \sup_{\xi \in \mathcal{J}} m_{g}(\varepsilon)$, and $n^*_{g} = \sup_{\xi \in \mathcal{J}} n_{g}(\varepsilon)$ such that

(H5) There exist $l_{\lambda}$, $m_{h}, n_{h} \in \mathbb{Y}$ with $l^{*}_{\lambda} = \sup_{\xi \in \mathcal{J}} l_{\lambda}(\varepsilon)$, $m^*_{\lambda} = \sup_{\xi \in \mathcal{J}} m_{\lambda}(\varepsilon)$, and $n^*_{\lambda} = \sup_{\xi \in \mathcal{J}} n_{\lambda}(\varepsilon)$ such that

(H6) There exist $M_{h{k}} > 0$, for all $u \in R^{n}$, such that

4.   Existence and uniqueness of solution

To prove the existence and uniqueness of solution first, we need to transform the problem (1.1) into a fixed point problem and define an operator $\mathscr{H}$ : $\mathbb{Y}\rightarrow \mathbb{Y}$ by

Let $x : [0-r, b+h]\rightarrow R$ be a function defined by

For each $z \in C([0, b], R)$ with $z(0) = 0$, we denote by $u$ the function defined by

Let us set $\Gamma(\varepsilon) = z(\varepsilon)+x(\varepsilon)$ such that $y_\varepsilon = z_\varepsilon + x_\varepsilon$ for each $\varepsilon \in (0, b],$ where

Let $\Phi : \mathbb{Y} \rightarrow \mathbb{Y}$ be an operator given by

To show that the operator $\mathscr{H}$ has a fixed point, for this it is sufficient to show that the operator $\Phi$ has a fixed point and this fixed point will correspond to a solution of problem (1.1).

Theorem 4.1. Assume that the Hypothesis $(H1)$ holds, then the solution (1.1) has a unique solution of the problem (1.1), provided

where

Proof. Consider the operator $\Phi$ : $\mathbb{Y} \rightarrow \mathbb{Y}$ defined by

As $\Delta$, h, g, $\lambda$, $h_{k}$ are all continuous, we have to prove that $\Phi$ is a contraction.

Case 1. For $z$, $y \in \mathbb{Y}$ and for $\varepsilon \in [0, \varepsilon_{1}]$, we have

This implies

Thus we take

So, we get

Case 2. For $\varepsilon \in (\varepsilon_{k}, s_{k}]$, $k = 1, 2, \dots, m$, we have

Take $L_{k} : = M_{h{k}}$ and therefore,

Case 3. For $z, y \in \mathbb{Y}$ and for $\varepsilon \in (s_{k}, \varepsilon_{k+1}]$, we have

Thus,

From the above three cases, we obtain that $||(\Phi z)(\varepsilon)-(\Phi y)(\varepsilon)||^2$ $\leq$ $L_{0} ||z(\varepsilon)- y(\varepsilon)||^{2}_{\mathbb{Y}}$ as per (4.1), $\Phi$ is contraction, and, thus, (1.1) has a unique fixed point $z$, which is a solution to problem (1.1).

Remark 4.1. Banach contraction principle provides not only the existence results, but also uniqueness is assured. Also, the nonlinear functions could satisfy only Lipschitz conditions to prove existence and uniqueness results, even though conditions are stronger.

Remark 4.2. It is to be noted that the assumption $L_0 < 1$ in Theorem 4.1 shows a restrictive smallness on Lipschitz constants for the nonlinear functions h, g, $\Delta,$ and $\lambda$ when compared with the periods of time, while the impulses are active or vice-versa. In order to relax such a kind of smallness, the generalized Schaefer's FPT is introduced.

Theorem 4.2. Assume that $(H2)-(H5)$ hold, then the nonlinear operator $\Phi:\mathbb{Y}\rightarrow \mathbb{Y}$ has a fixed point, which is a solution of problem (1.1)

Proof. Since $\Phi$ is well defined, we will present the proof by the following four steps.

First, we prove that $\Phi$ is completely continuous (CC), that is, $\mathcal{T}$ is continuous, maps bounded sets into bounded sets, and maps bounded sets into quasi-equicontinuous sets.

Step 1. To prove $\Phi$ is continuous.

Since h, g, $\Delta, \lambda$, and $h_{\mathcal{K}}$ are all continuous, then it is clear that the operator $\Phi$ is continuous on $\mathcal{J}$.

Step 2. To prove $\Phi$ maps bounded sets in $\mathbb{Y}$.

In fact, it is sufficient to show that for any $r > 0$, there exists an $\eta > 0$ such that for each $z \in \mathcal{B}_{r} = \{z \in \mathbb{Y}: ||z||^2_{\mathbb{Y}}\leq z\}.$

Case 1. For $\varepsilon \in [0, \varepsilon_1], z \in \mathbb{Y}$

Using $(H2)-(H4)$ and Lemmas 2.1 and 2.2, one can enumerate

Case 2. For $\varepsilon \in (\varepsilon_{k}, s_{k}]$, k = 1, 2, $\dots$, m,

Case 3. For $\varepsilon \in (s_{k}, \varepsilon_{k+1}]$, $k = 1, 2, \dots, m$, $\Gamma \in \mathbb{Y}$, we have

Thus

Let $\eta$ = $\{\eta_{0}, \eta_{k}, \eta_{b}\}$, $k = 1, 2, \dots, m$, then $\{\psi^{1-v}(\Phi z)(\varepsilon):z \in \mathscr{B} \}$ is a bounded set in $\mathbb{Y}$, i.e., $||\Phi z||^{2}_{\mathbb{Y}}$ $\leq \eta.$

Step 3. $\Phi$ is a quasi-equicontinous set in $\mathbb{Y}$.

Case 1. Let $\tau_{1}$, $\tau_{2}$ $\in [0, \varepsilon_{1}]$ with $0\leq\tau_{1}\leq\tau_{2}\leq\varepsilon_{1}$. One could establish the following estimate:

We conclude that as $\tau_2 - \tau_1 \rightarrow 0$ with $\varepsilon$ sufficiently small, the right hand side of the above inequality tends to zero independently of $z \in \mathscr{B_{r}}$. Furthermore, the similar results are true for $\varepsilon_{k} < \tau_1 < \tau_2 \leq s_{k}$ and $s_{k} < \tau_1 < \tau_2 \leq \varepsilon_{k+1}$ for $k = 1, 2, \dots, m$. It proves the equicontinuous of $\Phi$ on $\mathbb{Y}$. Thus, for $\tau_1, \tau_2 \in [0, b] \bigcap (\varepsilon_{k}, \varepsilon_{k+1}]$, $k = 1, 2, \dots, m$, whenever $\mathscr{B}$ is a bounded set of $\mathbb{Y}$ as in Step 2, let $z \in \mathscr{B}_r$, then

Thus, $\Phi$ is quasi-equicontinuous, then $\Phi(z)$ is relatively compact by Lemma 3.1, which implies that $\Phi(z)$ is CC.

Step 4. $\Phi$ has a Prior bound.

It remains to estimate that the set $\prod(\Phi) = \{z \in \mathbb{Y}:z = \Theta\Phi z, 0 < \Theta < 1\}$ is bounded.

Let $z \in \Pi(\Phi)$, then $z = \Theta\Phi(z)$ for some $\Theta\in (0, 1)$, by following the proof of Step 2 that $||z||_{\mathbb{Y}} \leq \eta.$ This proves that the set $\Pi(\Phi)$ is bounded. Hence, by Theorem 3.1, $\Phi$ has a fixed point, which is the required solution on $\mathscr{J}$.

5.   UHR stability results

Here, we derive the UHR stability for (1.1). Let $\omega > 0$, $\phi\geq 0$, and $\zeta \in PC(\mathscr{J}, R^{n})$ be nondecreasing. Consider the following inequalities.

Define a vector space as $\chi$:

Definition 5.1. System (1.1) is UHR stable with respect to $(\zeta, \phi)$ if there exists $C_{(M, L, P, \zeta)} > 0$ such that for each solution $\tilde{\Gamma} \in \chi$ of the inequality (5.1), there exists solution $\Gamma \in PC(\mathscr{J}, R^{n})$ of (1.1) with

Remark 5.1. A function $\tilde{\Gamma} \in \chi$ is a solution of (5.1) $\Leftrightarrow$ there is $\mathscr{Q} \in \bigcap_{i = 0}^{\mathscr{K}}$ $((s_{k}, \varepsilon_{k+1}], R^{n}]$ and $q \in \bigcap_{i = 0}^{\mathscr{K}} ((\varepsilon_{k}, s_{k}], R^{n})$ such that:

(1) $E||\mathscr{Q}(\varepsilon)||^2 \leq \omega\zeta(\varepsilon), \varepsilon \in \bigcap_{i = 0}^{\mathscr{K}}(s_{k}, \varepsilon_{k+1}]; E||q(\varepsilon)||^2 \leq \omega \phi, \varepsilon \in \bigcap_{i = 0}^{\mathscr{K}}(t_{k}, s_{k}];$

(2) $D^{\gamma, \beta}_{\psi(\varepsilon)}[\tilde{\Gamma}(\varepsilon)-h(\varepsilon, \tilde{\Gamma}(\varepsilon))] = A\tilde{\Gamma}(\varepsilon)+ \Delta(\varepsilon, \tilde{\Gamma}(\varepsilon), D^{\gamma, \beta}_{\psi(\varepsilon)}\tilde{\Gamma}(s)))+\int_{0}^{\varepsilon}g(s, \tilde{\Gamma}(s), D^{\gamma, \beta}_{\psi(\varepsilon)}\tilde{\Gamma}(s)))dB(s)$ $+\int_{0}^{\varepsilon}\lambda(s, \tilde{\Gamma}(s), D^{\gamma, \beta}_{\psi(\varepsilon)}\tilde{\Gamma}(s)))dB^{H}(s)+\mathscr{Q}(\varepsilon), \quad \varepsilon\ \in (s_{k}, \varepsilon_{k+1}], \quad k = 0, 1, 2, \dots, m,$

(3) $\tilde{\Gamma}(\varepsilon) = h_{k}(\varepsilon, \tilde{\Gamma}(\varepsilon))+q(\varepsilon), \varepsilon \in (\varepsilon_{k}, s_{k}], k = 0, 1, 2, \dots, m.$

By Remark 5.1, we have

Lemma 5.1. Let $\beta \in [0, 1]$, $\gamma \in (0, 1)$. If a function $\tilde{\Gamma} \in \chi$ is a solution of (5.1) then we have:

By Remark 5.1, we have:

Case 1. For $\varepsilon \in [0, \varepsilon_1]$, we have

Thus

From above, we obtain

Case 2. For $\varepsilon \in (\varepsilon_{k}, s_{k}]$, we have

Case 3. For $\varepsilon \in (s_{k}, \varepsilon_{k+1}]$, we get

Thus

In order to prove stability, we assume the following assumptions:

(H7) There exist positive constants $M_{k} (k = 1, 2, \dots, m)$ such that

(H8) Let $\zeta \in C(\mathscr{J}, R^{n})$ be a nondecreasing function if there exists $c_\zeta > 0$ such that $\int_{0}^{\varepsilon}\zeta(s)ds < c_{\zeta}\zeta(\varepsilon)$, $\forall$ $\varepsilon \in \mathscr{J}.$

Lemma 5.2. Let $P_0 = P \cup {0}$ where $p = {1, 2, \dots, m}$, and the following inequality holds

where $\Gamma, a, b \in PC(R^+, R^+)$, a is nondecreasing and $b(\varepsilon) > 0$, $\alpha_{k} > 0$, $k \in P$. For $\varepsilon \in R^{+}$,

where $\alpha = \sup_{K \in P}\{\alpha_K\}$ and $\varepsilon_0 = 0.$

Theorem 5.1. If the assumptions $(H1), (H7)$, and $(H8)$ are satisfied, then (1.1) is UHR stable with respect to $(\zeta, \phi)$.

Proof. Let $\tilde{\Gamma} \in \chi$ be a solution of inequality (5.1) and $\Gamma$ be the unique solution of (1.1).

Case 1. For $\varepsilon \in [0, \varepsilon_1]$, we have:

where $c_{p} = \frac{b^{2\gamma}}{\gamma^2} M_2.$

Case 2. For $\varepsilon \in (\varepsilon_{k}, s_{k}]$, we have

Case 3. For $\varepsilon \in (s_{k}, \varepsilon_{k+1}]$, we have

Hence, for $\varepsilon \in [0, \varepsilon_1]$, we have:

For $\varepsilon \in (\varepsilon_k, s_k]$, we get

For $\varepsilon \in (s_k, \varepsilon_{k+1}]$, k = 1, 2, $\dots$, m,

Combining (5.2)–(5.4), one can get an inequality of the form given in Lemma 5.2. For $\varepsilon \in \mathscr{J}$, since $\varepsilon \in (\varepsilon_k, \varepsilon_{k+1}]$, $K \in P_0$, we have

where

Thus,

where $C_{(M, L, p, {\zeta})}$ is a constant depending on $M, L, p, \zeta$. Hence, (1.1) is UHR stable w.r.t $(\zeta, \phi)$.

6.   Example

Consider the following nonlinear $\psi$-HFSE with NI impulses driven by both noises. This type of fractional SDE can be applied in pharmacotherapy.

where $\Gamma \in R^2$, $\Gamma \in \frac{1}{2}$, $\beta = \frac{3}{4}$, $0 < \varepsilon_0 < s_0 < \varepsilon_1 < \dots < s_m = 1$ are prefixed numbers, $\mathscr{J} = [0, 1]$. Here

Calculate the M-L function by using ^[12]. We have

where

Also calculate

where

Therefore, we need to check the hypotheses of nonlinear functions:

We get $M_1 = 0.0862$, $M_2 = 0.3824$, $M_{hk} = 0.065$, and $\gamma = 0.5$. Hence we have $L = max\{L_0, L_k, L_k*\} = 0.52 < 1.$ Also,

In this example, all the conditions stated in Theorems 4.1 and 5.1 are satisfied, so the example has a unique solution and is also UHR stable.

7.   Conclusions

With the help of Schaefer's FPT and the Banach contraction principles, we obtained the existence and uniqueness results for HSFEs with retarded and advanced arguments, selected the non-instantaneous impulses with both multiplicative and fractional noises, and obtained the UHR stability for HSFEs. UHR stability gives bounds between the exact and approximation solution, which is why this theory is very important in the numerical analysis as well as in approximation theory. We are hoping that our findings will have a great importance in the mentioned theories. Finally, a case study is provided to demonstrate the efficacy of the suggested outcomes. In the future, we can use the findings to investigate the controllability of HSFEs.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was supported by Anhui Province Natural Science Research Foundation (2023AH051810).

Conflicts of interest

The authors declare no conflicts of interest.