Analysis of a Time-Fractional Substantial Diffusion Equation of Variable Order

Abstract

A time-fractional substantial diffusion equation of variable order is investigated, in which the variable-order fractional substantial derivative accommodates the memory effects and the structure change of the surroundings of the physical processes with respect to time. The existence and uniqueness of the solutions to the proposed model are proved, based on which the weighted high-order regularity of the solutions, in which the weight function characterizes the singularity of the solutions, are analyzed.

1. Introduction

Fractional problems with substantial derivatives have attracted increasing amounts of attention due to their applications in, e.g., characterizing the functionals of anomalous diffusion and non-Brownian motions [1,2]. There are several investigations in mathematical analysis and numerical methods to such kinds of problems [3,4,5,6,7,8,9], while the corresponding studies for the variable-order fractional substantial derivative models, in which the variable fractional order accommodates, e.g., the structure change of the surroundings with respect to time, are rarely found in the literature.

In some recent works, the following Caputo variable-order time-fractional diffusion equations are investigated mathematically and numerically [10,11,12,13,14,15,16]

$$\begin{array}{c}{\displaystyle \left({\partial }_t+q{\partial }_t^{\alpha (·,t)}-\Delta \right)u(\mathit{x},t)=f(\mathit{x},t).}\end{array}$$

Here, the Caputo variable-order fractional derivative operator is defined as [17,18,19]

$$\begin{array}{c}{\displaystyle {\partial }_t^{\alpha (·,t)}g(t):={\int }_0^t\frac{{(t-s)}^{-\alpha (s,t)}}{\Gamma (1-\alpha (s,t))}{\partial }_{s}g(s)ds.}\end{array}$$

The Caputo (variable-order) fractional derivative is widely used in various applications that exhibit the power law type memory effects due to the power function integral kernel. The Caputo variable-order fractional substantial derivative operator ${}^{\sigma }{\partial }_t^{\alpha (·,t)}$, which is defined in the following for some $\sigma \ge 0$ and $0\le \alpha (s,t)<1$ [1,2,3]

$$\begin{array}{c}{\displaystyle {}^{\sigma }{\partial }_t^{\alpha (·,t)}g(t):={\int }_0^t\frac{{(t-s)}^{-\alpha (s,t)}{e}^{-\sigma (t-s)}}{\Gamma (1-\alpha (s,t))}\left({\partial }_s+\sigma \right)g(s)ds,}\end{array}$$

(1)

tempers the memory effects of the Caputo fractional derivative operator by introducing the exponential factor, such that the impacts located at s are negligible if $t-s$ is large [1,3]. Due to the complicated form of the variable-order fractional substantial derivative, existing analysis techniques for investigating the Caputo variable-order fractional derivative problems may not apply directly for problems involving (1), which motivates the current work.

In this paper, we analyze the following time-fractional diffusion equation [20,21,22,23,24,25,26,27,28] equipped with variable-order fractional substantial derivative

$$\begin{array}{c}{\displaystyle \left({\partial }_t+q{}^{\sigma }{\partial }_t^{\alpha (·,t)}-\Delta \right)u(\mathit{x},t)=f(\mathit{x},t),(\mathit{x},t)\in \Omega \times (0,T];}\\ {\displaystyle u(\mathit{x},0)=u_0(\mathit{x}),\mathit{x}\in \Omega ;u(\mathit{x},t)=0,(\mathit{x},t)\in \partial \Omega \times [0,T],}\end{array}$$

(2)

where $q\ge 0$, $T>0$, $\Omega$ is a d-dimensional domain for $1\le d\le 3$ with a smooth boundary $\partial \Omega$. In the following two sections, we rigorously prove the well-posedness and high-order regularity of the solutions to problem (2), respectively, which provide a theoretical basis for this model. In particular, a weight function is introduced in the estimate of high-order smoothing properties to characterize the initial singularity of the solutions.

In the rest of the paper, let $H^1(0,T)$ be the Hilbert space of functions with weak derivatives up to order 1 in $L^2(0,T)$, the space of square-integrable functions on $[0,T]$. The corresponding norms are defined by

$$\begin{array}{c}{\displaystyle {\parallel g\parallel }_{L^2(0,T)}:={\left({\int }_0^T{g}^2(t)dt\right)}^{1/2},}\\ {\displaystyle {\parallel g\parallel }_{H^1(0,T)}:={\left({\parallel g\parallel }_{L^2(0,T)}^2+{\parallel {\partial }_{t}g\parallel }_{L^2(0,T)}^2\right)}^{1/2}.}\end{array}$$

$H^1(0,T;\mathcal{D})$ for some Banach space $\mathcal{D}$ stands for the space of the functions in $H^1(0,T)$ with respect to the norm of $\mathcal{D}$. Q denotes a generic constant that may assume different values at different cases.

2. Existence and Uniqueness

We prove the existence and uniqueness of the solutions to problem (2) in this section. Throughout the paper, we assume that there exists a constant $0<{\alpha }^{\ast }<1$, such that

$$0\le \alpha (s,t)\le {\alpha }^{\ast }\mathrm{for}0\le s\le t\le T.$$

Theorem 1.

If $f\in L^2(0,T;L^2(\Omega ))$ and $u_0\in {\check{H}}^1(\Omega )$, problem (2) has a unique solution $u\in H^1(0,T;L^2(\Omega ))\cap L^2(0,T;{\check{H}}^2(\Omega ))$ and

$${\parallel u\parallel }_{H^1(0,T;L^2(\Omega ))}+{\parallel u\parallel }_{L^2(0,T;{\check{H}}^2(\Omega ))}\le Q\left({\parallel f\parallel }_{L^2(0,T;L^2(\Omega ))}+{\parallel u_0\parallel }_{{\check{H}}^1(\Omega )}\right)$$

where Q is independent from the data and the solutions and the norm ${\parallel ·\parallel }_{{\check{H}}^p(\Omega )}$ for $p\ge 0$ is defined in terms of the eigenpairs ${\{{\lambda }_i,{\varphi }_i(\mathit{x})\}}_{i=1}^{\infty }$ of $-\Delta$ by

$${\parallel v\parallel }_{{\check{H}}^p(\Omega )}^2:=\sum _{i=1}^{\infty }{\lambda }_i^p{(v,{\varphi }_i)}^2.$$

Proof. 

If we expand u in terms of ${\left\{{\varphi }_i(\mathit{x})\right\}}_{i=1}^{\infty }$ as

$$u(\mathit{x},t)=\sum _{i=1}^{\infty }{u}_i(t){\varphi }_i(\mathit{x}),u_i(t)=(u,{\varphi }_i),$$

then ${\left\{u_i(t)\right\}}_{i=1}^{\infty }$ satisfy the following ordinary differential equation

$${\partial }_{t}w(t)+q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w(t)+\lambda w(t)=h(t),t\in (0,T];w(0)=w_0$$

(3)

where w, $\lambda$, h and $w_0$ refer to $u_i$, ${\lambda }_i$, $(f,{\varphi }_i)$ and $(u_0,{\varphi }_i)$, respectively, for $i\ge 1$. We consider the fractional term as part of the right-hand side term in (3) to obtain

$${\partial }_{t}w(t)+\lambda w(t)=-q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w(t)+h(t).$$

An application of the variation-of-constants formula yields

$$w=e^{-\lambda t}{w}_0+e^{-\lambda t}\ast \left(-q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w+h\right),$$

(4)

where * stands for the convolution on $[0,t]$ as follows

$$g_1(t)\ast g_2(t):={\int }_0^t{g}_1(s)g_2(t-s)ds.$$

Then, we differentiate (4) with respect to t, multiply $e^{-\kappa t}$ for some $\kappa \ge 0$ on the resulting equation and reformulate $e^{-\kappa t}·{}^{\sigma }{\partial }_t^{\alpha (·,t)}$ as a sum of a linear operator $\mathcal{S}$ and a function $\mathcal{H}(t)$

$$\begin{array}{c}{\displaystyle e^{-\kappa t}·{}^{\sigma }{\partial }_t^{\alpha (·,t)}w:={\int }_0^t\frac{{(t-s)}^{-\alpha (s,t)}{e}^{-(\sigma +\kappa )(t-s)}}{\Gamma (1-\alpha (s,t))}}\hfill \\ {\displaystyle \times \left(e^{-\kappa s}{\partial }_{s}w(s)+\sigma {\int }_0^s{e}^{-\kappa (s-y)}{e}^{-\kappa y}{\partial }_{y}w(y)dy+\sigma e^{-\kappa s}{w}_0\right)ds}\hfill \\ {\displaystyle =:\mathcal{S}(e^{-\kappa t}{\partial }_{t}w)+\mathcal{H}(t),}\hfill \\ {\displaystyle \mathcal{H}(t):=\sigma w_0{e}^{-\kappa t}{\int }_0^t\frac{{(t-s)}^{-\alpha (s,t)}{e}^{-\sigma (t-s)}}{\Gamma (1-\alpha (s,t))}ds}\hfill \end{array}$$

to obtain an integral equation in terms of $\tilde{w}:=e^{-\kappa t}{\partial }_{t}w$ as follows

$$\begin{array}{c}{\displaystyle \tilde{w}=-\lambda e^{-(\lambda +\kappa )t}{w}_0+\left(-q\mathcal{S}(\tilde{w})-q\mathcal{H}(t)+e^{-\kappa t}h\right)}\hfill \\ {\displaystyle -\lambda e^{-(\lambda +\kappa )t}\ast \left(-q\mathcal{S}(\tilde{w})-q\mathcal{H}(t)+e^{-\kappa t}h\right).}\hfill \end{array}$$

(5)

In order to prove the well-posedness of this integral equation in $L^2(0,T)$, it suffices to show that the linear operator $-qS+\lambda q{e}^{-(\lambda +\kappa )t}\ast \mathcal{S}$ is a contraction in $L^2(0,T)$. For $v\in L^2(0,T)$, we apply

$${|(t-s)}^{-\alpha (s,t)}{|=(t-s)}^{-{\alpha }^{\ast }}{(t-s)}^{{\alpha }^{\ast }-\alpha (s,t)}\le max\{1,T\}{(t-s)}^{-{\alpha }^{\ast }}$$

and Young’s convolution inequality to bound $\mathcal{S}v$ by

$$\begin{array}{cc}\hfill & {\displaystyle {\parallel \mathcal{S}v\parallel }_{L^2(0,T)}}\hfill \\ \hfill & {\displaystyle \le Q\parallel {\int }_0^t\frac{e^{-(\sigma +\kappa )(t-s)}}{{(t-s)}^{{\alpha }^{\ast }}}\left(\left|v(s)\right|+\sigma {\int }_0^s{e}^{-\kappa (s-y)}\left|v(y)\right|dy\right)ds{\parallel }_{L^2(0,T)}}\hfill \\ \hfill & {\displaystyle \le Q\parallel e^{-(\sigma +\kappa )t}{t}^{-{\alpha }^{\ast }}{\parallel }_{L^1(0,T)}{\parallel v\parallel }_{L^2(0,T)}}\hfill \\ \hfill & {\displaystyle \le Q{(\sigma +\kappa )}^{-(1-{\alpha }^{\ast })}{\parallel v\parallel }_{L^2(0,T)}.}\hfill \end{array}$$

The operator $\lambda e^{-(\lambda +\kappa )t}\ast \mathcal{S}$ is then accordingly bounded by

$$\parallel \lambda e^{-(\lambda +\kappa )t}{\ast \mathcal{S}v\parallel }_{L^2(0,T)}\le {\parallel \mathcal{S}v\parallel }_{L^2(0,T)}\le Q{(\sigma +\kappa )}^{-(1-{\alpha }^{\ast })}{\parallel v\parallel }_{L^2(0,T)}.$$

We incorporate the above two equations to conclude that for $\kappa$ large enough, $-qS+\lambda q{e}^{-(\lambda +\kappa )t}\ast \mathcal{S}$ is a contraction in $L^2(0,T)$ and, thus, the integral Equation (5) has a unique solution $\tilde{w}\in L^2(0,T)$. We further apply this contractability on (5) to obtain a stability estimate

$$\begin{array}{cc}\hfill {\displaystyle \parallel \tilde{w}{\parallel }_{L^2(0,T)}}& {\displaystyle \le Q\left({\lambda }^{1/2}|w_0{|+\parallel \mathcal{H}\parallel }_{L^2(0,T)}+{\parallel h\parallel }_{L^2(0,T)}\right)}\hfill \\ \hfill & {\displaystyle \le Q\left({\lambda }^{1/2}|w_0{|+\parallel h\parallel }_{L^2(0,T)}\right).}\hfill \end{array}$$

Then, the w defined by

$$w:={\int }_0^t{e}^{\kappa s}\tilde{w}(s)ds+w_0\in H^1(0,T)$$

with the estimate

$${\parallel w\parallel }_{H^1(0,T)}\le Q{\parallel \tilde{w}\parallel }_{L^2(0,T)}\le Q\left({\lambda }^{1/2}|w_0{|+\parallel h\parallel }_{L^2(0,T)}\right)$$

satisfies the differential Equation (3). The unique $H^1$ solution to (3) follows from that of (5).

Based on the above derivations, we finally estimate the solution u to problem (2) by

$$\begin{array}{c}{\displaystyle {\parallel u\parallel }_{H^1(0,T;L^2(\Omega ))}^2\le Q\sum _{i=1}^{\infty }\parallel u_i{\parallel }_{H^1(0,T)}^2\le Q\sum _{i=1}^{\infty }({\lambda }_i{\left|(u_0,{\varphi }_i)\right|}^2}\hfill \\ {\displaystyle +\parallel (f,{\varphi }_i){\parallel }_{L^2(0,T)}^2)=Q({\parallel f\parallel }_{L^2(0,T;L^2(\Omega ))}^2+{\parallel u_0\parallel }_{{\check{H}}^1(\Omega )}^2).}\hfill \end{array}$$

The estimate of the solutions in ${\parallel ·\parallel }_{L^2(0,T;{\check{H}}^2(\Omega ))}$ follows from

$$\begin{array}{cc}\hfill {\displaystyle {\parallel u\parallel }_{L^2(0,T;{\check{H}}^2(\Omega ))}}& {\displaystyle ={\parallel \Delta u\parallel }_{L^2(0,T;L^2(\Omega ))}}\hfill \\ & {\displaystyle =\parallel {\partial }_{t}u+q{}^{\sigma }{\partial }_t^{\alpha (·,t)}u-f{\parallel }_{L^2(0,T;L^2(\Omega ))}}\hfill \\ & {\displaystyle \le Q({\parallel f\parallel }_{L^2(0,T;L^2(\Omega ))}+{\parallel u_0\parallel }_{{\check{H}}^1(\Omega )}).}\hfill \end{array}$$

The uniqueness of the solutions to problem (2) follows from that of the differential Equation (3). □

3. Weighted Regularity

We prove the weighted regularity of the solutions to problem (2) in the following theorem.

Theorem 2.

Suppose α has bounded first-order derivatives for each coordinate, $u_0\in {\check{H}}^{s+3}(\Omega )$, $f\in L^2(0,T;{\check{H}}^{s+2}(\Omega ))$ and ${\partial }_{t}f\in L^2(0,T;{\check{H}}^s(\Omega ))$ for some $s\ge 0$. Then, the following estimate holds

$$\begin{array}{c}{\displaystyle \parallel t^{\delta }{\partial }_t^2{u\parallel }_{L^2(0,T;{\check{H}}^s(\Omega ))}}\hfill \\ {\displaystyle \le Q(\parallel u_0{\parallel }_{{\check{H}}^{s+3}(\Omega )}+\parallel {\partial }_t{f\parallel }_{L^2(0,T;{\check{H}}^s(\Omega ))}+{\parallel f\parallel }_{L^2(0,T;{\check{H}}^{s+2}(\Omega ))})}\hfill \end{array}$$

where

$$\delta :=max\{\alpha (0,0)-1/2+\theta ,0\}with0<\theta \ll 1$$

and Q is independent from the data and the solutions.

Proof. 

We differentiate (4) twice in time to obtain

$$\begin{array}{c}{\displaystyle {\partial }_t^{2}w={\lambda }^2{e}^{-\lambda t}{w}_0+{\partial }_t\left(-q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w+h\right)}\hfill \\ {\displaystyle -\lambda \left(-q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w+h\right)+{\lambda }^2{e}^{-\lambda t}\ast \left(-q{}^{\sigma }{\partial }_t^{\alpha (·,t)}w+h\right).}\hfill \end{array}$$

(6)

The ${\partial }_t{}^{\sigma }{\partial }_t^{\alpha (·,t)}w$ could be simplified as

$$\begin{array}{cc}\hfill {\displaystyle \left|{\partial }_t{}^{\sigma }{\partial }_t^{\alpha (·,t)}w(t)\right|}& {\displaystyle =\left|{\partial }_t{\int }_0^t\frac{s^{-\alpha (t-s,t)}{e}^{-\sigma s}}{\Gamma (1-\alpha (t-s,t))}\left({\partial }_y+\sigma \right)w(y){|}_{y=t-s}ds\right|}\hfill \\ \hfill & {\displaystyle =|\frac{t^{-\alpha (0,t)}{e}^{-\sigma t}}{\Gamma (1-\alpha (0,t))}\left({\partial }_{t}w(0)+\sigma w(0)\right)}\hfill \\ \hfill & {\displaystyle +{\int }_0^t{\partial }_t\left(\frac{s^{-\alpha (t-s,t)}{e}^{-\sigma s}}{\Gamma (1-\alpha (t-s,t))}\right)\left({\partial }_y+\sigma \right)w(y){|}_{y=t-s}ds}\hfill \\ \hfill & {\displaystyle +{\int }_0^t\frac{s^{-\alpha (t-s,t)}{e}^{-\sigma s}}{\Gamma (1-\alpha (t-s,t))}\left({\partial }_y^2+\sigma {\partial }_y\right)w(y){|}_{y=t-s}ds|}\hfill \\ \hfill & {\displaystyle \le Q{t}^{-\alpha (0,t)}\left(|{\partial }_{t}w(0)|+|w(0)|\right)+Q{\int }_0^t\frac{|{\partial }_s^{2}u(s)|}{{(t-s)}^{\alpha (s,t)}}ds.}\hfill \end{array}$$

(7)

As

$$t^{-\alpha (0,t)}=t^{-\alpha (0,0)}{t}^{\alpha (0,0)-\alpha (0,t)}\le Q{t}^{-\alpha (0,0)},$$

we multiply $t^{\delta }$ on both sides of (7) and apply

$${(t-s)}^{\alpha (s,t)}\ge Q{(t-s)}^{\alpha (t,t)}$$

to find

$$\begin{array}{cc}\hfill {\displaystyle \left|t^{\delta }{\partial }_t{}^{\sigma }{\partial }_t^{\alpha (·,t)}w(t)\right|}& {\displaystyle \le Q{t}^{\delta -\alpha (0,0)}\left(|{\partial }_{t}w(0)|+|w(0)|\right)}\hfill \\ \hfill & {\displaystyle +Q{\int }_0^t\frac{s^{\delta }\left|{\partial }_s^{2}u(s)\right|}{{(t-s)}^{{\alpha }^{\ast }}}ds+Q{\int }_0^t\frac{{(t-s)}^{\delta }\left|{\partial }_s^{2}u(s)\right|}{{(t-s)}^{\alpha (t,t)}}ds.}\hfill \end{array}$$

(8)

Since

$$\begin{array}{c}{\displaystyle {\int }_0^t\frac{{(t-s)}^{\delta }\left|{\partial }_s^{2}u(s)\right|}{{(t-s)}^{\alpha (t,t)}}ds}\hfill \\ {\displaystyle ={\int }_0^t\frac{{(t-s)}^{\delta -\alpha (t,t)/2}}{s^{\delta }}\frac{s^{\delta }\left|{\partial }_s^{2}u(s)\right|}{{(t-s)}^{\alpha (t,t)/2}}ds}\hfill \\ {\displaystyle \le {\left({\int }_0^t\frac{{(t-s)}^{2\delta -\alpha (t,t)}}{s^{2\delta }}ds\right)}^{1/2}{\left({\int }_0^t\frac{s^{2\delta }{\left|{\partial }_s^{2}u(s)\right|}^2}{{(t-s)}^{\alpha (t,t)}}ds\right)}^{1/2}}\hfill \\ {\displaystyle \le Q{\left(t^{1-\alpha (t,t)}\right)}^{1/2}{\left({\int }_0^t\frac{s^{2\delta }{\left|{\partial }_s^{2}u(s)\right|}^2}{{(t-s)}^{\alpha (t,t)}}ds\right)}^{1/2}}\hfill \\ {\displaystyle \le Q{\left({\int }_0^t\frac{s^{2\delta }{\left|{\partial }_s^{2}u(s)\right|}^2}{{(t-s)}^{{\alpha }^{\ast }}}ds\right)}^{1/2},}\hfill \end{array}$$

we apply $\parallel e^{-\kappa t}{·\parallel }_{L^2(0,T)}$ for some $\kappa >0$ on (8) to find

$$\begin{array}{cc}\hfill {\displaystyle \parallel e^{-\kappa t}{t}^{\delta }{\partial }_t{}^{\sigma }{\partial }_t^{\alpha (·,t)}w(t){\parallel }_{L^2(0,T)}}& {\displaystyle \le Q\left(|{\partial }_{t}w(0)|+|w(0)|\right)}\hfill \\ \hfill & {\displaystyle +Q\parallel t^{-{\alpha }^{\ast }}{e}^{-\kappa t}{\parallel }_{L^1(0,T)}\parallel e^{-\kappa t}{t}^{\delta }{\partial }_t^{2}u{\parallel }_{L^2(0,T)}}\hfill \\ \hfill & {\displaystyle +Q{\left({\int }_0^T{\int }_0^t\frac{e^{-2\kappa (t-s)}}{{(t-s)}^{{\alpha }^{\ast }}}{\left(e^{-\kappa s}{s}^{\delta }{\partial }_s^{2}u(s)\right)}^{2}dsdt\right)}^{1/2}}\hfill \\ \hfill & {\displaystyle \le Q\left(|{\partial }_{t}w(0)|+|w(0)|\right)+Q{\kappa }^{-(1-{\alpha }^{\ast })/2}\parallel e^{-\kappa t}{t}^{\delta }{\partial }_t^{2}u{\parallel }_{L^2(0,T)}.}\hfill \end{array}$$

(9)

We, therefore, multiply $t^{\delta }$ on both sides of (6), apply $\parallel e^{-\kappa t}{·\parallel }_{L^2(0,T)}$ on the resulting equation, choosing $\kappa$ large enough and employ (9) and ${(t\lambda )}^{\delta }{e}^{-\lambda t/2}\le Q$, as well as similar (and simpler) proofs as above to obtain

$$\begin{array}{c}{\displaystyle \parallel e^{-\kappa t}{t}^{\delta }{\partial }_t^2{w\parallel }_{L^2(0,T)}\le Q\left({\lambda }^{3/2}|w_0|+\parallel {\partial }_t{h\parallel }_{L^2(0,T)}+\lambda {\parallel h\parallel }_{L^2(0,T)}\right).}\hfill \end{array}$$

We employ this result to obtain the estimate of ${\partial }_t^{2}u$ as

$$\begin{array}{c}{\displaystyle \parallel t^{\delta }{\partial }_t^2{u\parallel }_{L^2(0,T;{\check{H}}^s(\Omega ))}^2=\sum _{i=1}^{\infty }{\lambda }_i^s{\parallel t^{\delta }{\partial }_t^2{u}_i\parallel }_{L^2(0,T)}^2}\hfill \\ {\displaystyle \le Q\sum _{i=1}^{\infty }\left({\lambda }_i^{s+3}|(u_0,{\varphi }_i){|}^2+\parallel ({\partial }_{t}f,{\varphi }_i){\parallel }_{L^2(0,T)}^2+{\lambda }_i^2{\parallel (f,{\varphi }_i)\parallel }_{L^2(0,T)}^2\right)}\hfill \\ {\displaystyle =Q\left(\parallel u_0{\parallel }_{{\check{H}}^{s+3}(\Omega )}^2+\parallel {\partial }_t{f\parallel }_{L^2(0,T;L^2(\Omega ))}^2+{\parallel f\parallel }_{L^2(0,T;{\check{H}}^2(\Omega ))}^2\right),}\hfill \end{array}$$

which completes the proof. □

4. Conclusions

In this paper we rigorously prove the well-posedness and high-order regularity of the solutions to the time-fractional diffusion equation equipped with variable-order fractional substantial derivative. In particular, a weight function is introduced in the estimate of high-order smoothing properties to characterize the initial singularity of the solutions. Based on these theoretical estimates, we will carry out numerical analysis for the proposed model in the near future in order to provide supports for practical implementations.