1. Introduction
Consider the following initial-boundary value problem for the Navier-Stokes equations in two dimensions with non-local viscosity. It means, find a vector function
\[ {u}: \Omega \times [0,T] \rightarrow \mathbb{R}^2,\]
and a scalar function
\[p: \Omega \times [0,T] \rightarrow \mathbb{R},\]
satisfying
\begin{align}
&\dfrac{d {u}(x,t)}{dt} – c(l(u_{1}(x,t)),l(u_{2}(x,t))) \Delta {u}(x,t)+( {u}(x,t) \cdot \nabla ) {u}(x,t) + \nabla p(x) = {f}(x,t) & \mbox{ in } \Omega \times (0,T), \label{eqi1} \\
\end{align}
(1)
\begin{align}
&div( {u}(x))=0 & \mbox{ on } \Omega, \label{eqi2}\\
\end{align}
(2)
\begin{align}
& {u}(x,t)= {g} & \mbox{ on } \partial \Omega,\label{eqi3} \\
\end{align}
(3)
\begin{align}
& {u}(x,0)= {u}_0(x) &\mbox{ in } \Omega, \label{eqi4}
\end{align}
(4)
where \(\Omega\) is a domain sufficiently regular, \(\partial \Omega \) its boundary well regular, and we have that \(c(l(u_{1}(x,t)),l(u_{2}(x,t)))\) satisfies these hypotheses:
Given \(x=(x_1,x_2), y=(y_1,y_2) \in \mathbb{R}^2,\)
- (A1) \(0< c_- \leq c(x_1,x_2) \leq c_+,\)
- (A2) \(|c(x) – c(y)| \leq A_1 |x_1-y_1| + A_2 | x_2 – y_2|, \mbox{ for some } A_1,A_2>0,\)
and \(l:L^2(\Omega) \rightarrow \mathbb{R}\) is a continuous linear functional, defined by \(u \mapsto \int_\Omega u d\Omega.\)
We mentioned that the existence, uniqueness and exponential decay of the solution to the problem (1)-(4) were studied by Ferreira, Shahrouzi, Andrade and Panni in [1].
The motivation to study this kind of problem is we can describe motion of fluids which viscosity depends of time and satisfies the hypotheses \((A1)-(A2)\), and, when \(c(l(u_{1}(x,t)),l(u_{2}(x,t)))=\mu\), constant, we obtain the regular Navier-Stokes equations in two dimensions. This non-local term was introduced by Chipot [2], and it arrives naturally when we study the growth of a bacteria population, one kind of this problem was suggested by Ladyzhenskaya [3] where \(c(t):= \mu_0 + \mu_1 \| {u}(t)\|^2\) when \(\mu_0\) and \(\mu_1\) are positive constants.
The rest of the paper is organized as follows. In §2, we recall some notations, the weak formulation, lemmas and theorems. In §3, we study the existence of strong solutions to the problem (1)-(4). In §4, we introduce the existence of periodic solutions using Brouwer’s fixed point theorem. The conclusions of the paper are presented in §5.
2. Preliminaries
In this section we introduce some notations, achieve weak formulation and enunciate some important results.
2.1. Notations
Let \(\Omega \subset \mathbb{R}^2\) be a regular domain with \(\partial \Omega\) a well regular boundary. We denote the inner product in \(H^1_0(\Omega)\) by \(((\cdot,\cdot))\) and \((\cdot,\cdot)\) in \(L^2(\Omega)\), and norms respectively by \(\|\cdot\|\) and \(|\cdot|\). By \( {H}^1_0(\Omega)\) we denote \((H_0^1(\Omega))^2\) and, \( {L}^2(\Omega)\) by \((L(\Omega))^2\). The set \( {\mathfrak{V}}\) is the set of all distributions \( {u}: (\mathcal{D}(\Omega))^2 \mapsto \mathbb{R}^2\), which its divergent is null, in other words, \( {\mathfrak{V}}(\Omega):= \{ {u} \in (\mathcal{D}(\Omega))^2 ; div( {u}) =0\}\). Also, we denote the closure of \( {\mathfrak{V}}(\Omega)\) in \( {H}^1_0(\Omega)\) by \( {V}\) and the closure of \( {\mathfrak{V}}(\Omega)\) in \((L(\Omega))^2\) by \(H\).
A well known propriety of non-local term, see [4,5,6,7,8], is that this term commutes with spatial integral sign
\[\int_\Omega c(l(u_1),l(u_2)) {u} d\Omega = c(l(u_1),l(u_2)) \int_\Omega {u}d\Omega.\]
2.2. Weak formulation
Consider \( {v} \in {V}\). Doing inner-product in \(L^2(\Omega)\) with Eq. (1) we get,
\begin{align*}
&\left(\dfrac{d {u}(x,t)}{dt}, {v}(x)\right)- c(l(u_1),l(u_2))(\nabla {u}(x,t), {v}(x) ) + (( {u}(x,t) \cdot \nabla ) {u}(x,t), {v}(x))+ (\nabla p(x), {v}(x)), = ( {f}(x,t), {v}).
\end{align*}
By green first identity and integration by parts,
\begin{align*}
&\dfrac{d}{dt}\left( {u}(x,t), {v}(x)\right) + c(l(u_1),l(u_2))(( {u}(x,t), {v}(x) ))+ (( {u}(x,t) \cdot \nabla ) {u}(x,t), {v}(x)) = ( {f}(x,t), {v}(x)).
\end{align*}
Now we define a bilinear form \(a( {u}, {v}):=(( {u}, {v} ))\) and a trilinear form \(b( {u}, {v}, {w}):= (( {u} \cdot \nabla ) {v}, {w})\), and then we obtain the weak form of Eq. (1),
\begin{equation*}
\dfrac{d}{dt}\left( {u}, {v}\right) + c(l(u_1),l(u_2))a ( {u}, {v} ) + b( {u}, {u}, {v}) = ( {f}, {v}).
\end{equation*}
2.3. Some results
Lemma 1.[9] Let \( {u} \in {L}^2(0,T; {V})\), then the function \(B {u}\) defined by,
\[ \langle B {u}(t), {v} \rangle := b( {u}, {v}, {w}), \forall {v} \in {V}, \text{for a.e. } t \in [0,T],\]
belongs to \( {L}^1(0,T; {V}’).\)
Lemma 2. Let \(\Omega \subset \mathbb{R}^n\) be an bounded Lipschitz open set in \(\mathbb{R}^2\).
- If a distribution \(p\) has all its first-order derivatives \(D_i p\) in \(L^2(\Omega),\) then \(p \in L^2(\Omega)\) and
\[\|p\|_{L^2(\Omega) \backslash \mathbb{R} } \leq c(\Omega) |\nabla p|_{L^2(\Omega)}.\]
- If a distribution \(p\) has all its first-order derivatives in \(H^{-1}(\Omega).\) Then \(p \in L^2(\Omega)\) and
\[\|p\|_{L^2(\Omega) \backslash \mathbb{R}} \leq c \|\nabla p\|_{H^{-1}(\Omega)},\]
where \(L^2(\Omega)\backslash \mathbb{R}:= \left\{ p \in L^2(\Omega) \big\rvert \int_\Omega p(x) dx =0 \right\}.\)
Problem 1. For \( {f}\) and \( {u}_0\) given, with
\begin{align}
& {f} \in L^2(0,T; {V}’), \label{prob1}\\
\end{align}
(5)
\begin{align}
& {u}_0 \in H,\label{prob2}
\end{align}
(6)
to find \( {u}\) satisfying,
\begin{align*}
{u} \in L^2(0,T; {V}), {u}’ \in L^1(0,T; {V}’),
\end{align*}
\begin{align*}
&\dfrac{d}{dt}\left( {u}(x,t), {v}(x)\right) + c(l(u_1),l(u_2))(( {u}(x,t), {v}(x) ))
+ (( {u}(x,t) \cdot \nabla ) {u}(x,t), {v}(x)) = ( {f}(x,t), {v}(x)),
\end{align*}
for any \( {v} \in {V}\).
3. Existence of strong solutions
Suppose the existence of weak solutions to the problem (1)-(4). Our goal in this section is recover the pressure and prove the existence of strong solutions.
Theorem 2.
Given \( {f}\) and \( {u}_0\) satisfying (5) and (6). Suppose that \( {u}\) is a solution of the Problem 1 and \[ {f}-c(l(u_1),(u_2))A {u} -B {u} – {u}’ \in L^2(0,T;V’),\] then the solution \( {u}\) is also strong.
Proof.
Let,
\[ {U}(t):=\int_0^t {u}(s)ds, {F}(t):=\int_0^t {f}(s) ds \hspace{0.2cm} \mbox{ and } \beta(t):=\int_0^t B( {u}(s), {u}(s))ds \in V’.\]
Since \( {u}, {f},B {u} \in L^2(0,T;V’)\) then,
\begin{equation}
{U}, {F} \hspace{0.2cm} \mbox{ and } \beta \in C^0(0,T;V’) \mbox{ these are absolute continuous}.\label{106}
\end{equation}
(7)
Integrating \(c(l(u_1),(u_2))A {u} +B {u}+ {u}’ = {f}\), and, by (7), we get
\[ {u}(t) – {u}(0) + c(l(u_1),l(u_2))\int_0^t A {u}(s)ds + \int_0^t B {u}(s)ds = \int_0^t {f}(s) ds \mbox{ in V’.}\]
Then,
\[ {u}(t)- {u}_0 + c(l(u_1),l(u_2))A {U}(t) + \beta(t) = {F}(t) \mbox{ in }V’, \forall t \in [0,T].\]
So, for each \(\phi \in {\mathfrak{V}}\),
\begin{equation}
\langle {u}(t) – {u}(0) + c(l(u_1),l(u_2))A {U}(t)+\beta(t) – {F}(t),\phi \rangle =0. \label{107}
\end{equation}
(8)
Define,
\begin{equation}
{S}(t) := {u}(t) – {u}_0 + c(l(u_1),l(u_2))A {U}(t) + \beta(t) – {F}(t) \in V’. \label{108}
\end{equation}
(9)
For each \(t \in [0,T]\) it is possible to extend \( {S}(t)\) on a functional \( {T}(t)\in H^{-1}(\Omega)\) such as,
\begin{equation}
\langle {T}(t), {v} \rangle = \langle {S}(t), {v} \rangle, \forall {v} \in V.\label{109}
\end{equation}
(10)
But, from (8) and (10) we can conclude that,
\[\langle {T}(t), {\phi} \rangle =0, \forall {\phi} \in \mathfrak{V}.\]
From Lemma 2 results that \(\exists P(t) \in L^2(\Omega)\) satisfying,
\begin{equation}
{T}(t)= \nabla P(t) \mbox{ in }H^{-1}(\Omega). \label{110}
\end{equation}
(11)
So, from (10) and (11) we get,
\begin{equation}
\nabla P(t)\big\rvert_V \equiv {S}(t) \mbox{ in } V’, \forall t \in [0,T]. \label{111}
\end{equation}
(12)
Replacing (12) in (9),
\[ {u}(t) – {u}_0 + c(l(u_1),l(u_2))A {U}(t) + \beta(t) – {F}(t) = \nabla P(t) \mbox{ in }V’, \forall v \in [0,T].\]
As the expression on the left belongs to the space \( C ^ 0 (0, T; V ‘) \) we have \( \nabla P \in C ^ 0 (0, T; V’) \), and hence we can derive the above equation in the sense of distributions, with this:
\[ {u}’ + c(l(u_1),l(u_2))A {u} – {f} + B {u}= \nabla \dfrac{\partial P}{\partial t} \mbox{ in } L^2(0,T;V’).\]
Therefore is possible to say that equality above is given a.e. in \( (0, T) \). Setting \(p(x,t) = – \dfrac{\partial P}{\partial t},\) results in,
\[ {u}’ + c(l(u_1),l(u_2)) A {u} +B {u} = {f} -\nabla p \in L^2(0,T;V’).\]
4. Existence of periodic solutions
The purpose of this section is to prove the existence of periodic solutions to the Navier-Stokes equations.
Theorem 3.
Let \(\Omega \subset \mathbb{R}^2\) a bounded open set with boundary \(\partial \Omega\) well regular and \(Q:=[0,T] \times \Omega\). Consider the following problem,
\begin{align}
\begin{cases}
\dfrac{\partial {u}}{\partial t} – c(l(u_1),l(u_2)) \Delta {u} + ( {u} \cdot \Delta) {u} +\nabla p = {f} & \mbox{ in } Q, \\
div( {u}) =0 &\mbox{ in } Q, \\
{u}=0 & \mbox{ on } \partial \Omega,\\
{u}(x,0)= {u}(x,T), \forall x \in \Omega,
\end{cases} \label{perisollabel}
\end{align}
(13)
where \( {f} \in L^2(0,T;V’)\). This problem admits weak solution in \( {u}:Q \rightarrow \mathbb{R}^2,\) \( {u} \in L^2(0,T;H)\cap L^\infty(0,T;H)\) and \( {u}’ \in L^2(0,T;V’).\)
Proof.
The weak formulation of (13) is given by,
\begin{align}
\begin{cases}&\langle {u}'(t), {v} \rangle + c(l(u_1),l(u_2))(( {u}(t), {v})) + b( {u}(t), {u}(t), {v}) = \langle {f}(t), {v} \rangle \mbox{ in } \mathcal{D}'(0,T), \forall {v} \in V, \\
& {u}(0) = {u}(T).\label{P03}
\end{cases}\end{align}
(14)
Consider \(\{ {w}_1,\cdots, {w}_n,\cdots\}\) a base of \(V\). We truncate the series in \( m \)-th term, which leads to the approximate solution space \( V_m \). Setting \( {u}_m(t):=g_{im}(t) {w}_i,\)
\begin{align}
\begin{cases}
&( {u}’_m(t), {w}_j) + c(l(u_1),l(u_2))(( {u}_m(t), {w}_j)) + b( {u}_m(t), {u}_m(t), {w}_j) = \langle {f}(t), {w}_j \rangle, \\
& {u}_m(0)= {v} \in V_m,\label{p09}
\end{cases}\end{align}
(15)
where \(j=1,\cdots, m\).
The approximate system above has a global solution, since by similar procedure to the case of the existence of solutions [1], we obtain the following inequality,
\[| {u}(t)|^2 + \int_0^t \| {u}_m(s)\|^2ds \leq | {v}| + \frac{1}{c_-}\| {f}\|_{L^2(0,T;V’)}\leq c(m),\]
as \(m\) is fixed, we can extend \( {u}(t)\) in \([0,T].\) Our goal is to show that, among all solutions of the approximate equation, there is at least one \( {u}_m \) solution that satisfies periodicity,
\[ {u}_m(0) = {u}_m(T).\]
To do this, just prove that for every \( m \in \mathbb{N}, \) the application,
\begin{align}
& {\tau}_m: V_m \rightarrow V_m \nonumber \\
& {v} \mapsto {\tau}_m( {v}) = {u}_m(T), \nonumber
\end{align}
has a single fixed point, because in this case there will be a single function \( {v} \in V_m \) such that
\begin{equation}
{u}_m(T)= {\tau}_m( {v}) = {v} = {u}_m(0), \hspace{0.2cm} \forall m \in \mathbb{N}. \label{p10}
\end{equation}
(16)
Thus (16) we have a \( ( {u}_m) \) sequence of approximate solutions such that they all satisfy the periodicity condition.
Lemma 3. Exists \(\rho_0 >0\) such as \( {\tau}_m \overline{(B_{\rho_0}(0))} \subset \overline{B_{\rho_0}(0)}.\)
Proof. Using the \(H\) induced topology in \(V_m\), it suffices to prove that
\begin{align*}
\exists\;\;\;\; \rho_0 > 0 \mbox{ such that } | {\tau}_m( {v})|_H \leq \rho_0; \forall v \in V_m, \mbox{ where } | {v}|_H \leq \rho_0.
\end{align*}
Applying the energy method,
\begin{align*}
\frac{1}{2} \dfrac{d}{dt}| {u}_m(t)|^2 + c_- \| {u}_m(t)\|^2 &+ b( {u}_m(t), {u}_m(t), {u}_m(t))\\
&\leq \frac{1}{2} \dfrac{d}{dt}| {u}_m(t)|^2 + c(l(u_1),l(u_2)) \| {u}_m(t)\|^2 + b( {u}_m(t), {u}_m(t), {u}_m(t))\\
&=\langle {f}(t), {u}_m(t) \rangle\\
&\leq\| {f}(t)\|_{V’} | {u}_m(t)|,
\end{align*}
implies that,
\begin{align*}
\frac{1}{2} \dfrac{d}{dt}| {u}_m(t)|^2 + c_- \| {u}_m(t)\|^2 \leq \frac{1}{2c_-}\| {f}(t)\|^2_{V’} + \frac{c_-}{2}\| {u}_m(t)\|^2,
\end{align*}
then,
\begin{equation}
\frac{1}{2} \dfrac{d}{dt}| {u}_m(t)|^2 + \frac{c_-}{2} \| {u}_m(t)\|^2 \leq \frac{1}{c_-} \| {f}(t)\|^2_{V’}. \label{periodica2}
\end{equation}
(17)
As \(V \hookrightarrow H,\) exists \(c_0 >0\) such as,
\begin{equation}
c_0^2| {u}_m(t)|^2 \leq \| {u}_m(t)\|^2.\label{periodica3}
\end{equation}
(18)
Thus from (17) and (18) we get,
\[\dfrac{d}{dt}| {u}_m(t)|^2 + c_0^2 c_- | {u}_m(t)|^2 \leq \frac{1}{c_-} \| {f}(t)\|^2_{V’}.\]
Multiplying both sides by \(e^{c_0^2c_-t}\):
\begin{align*}
\dfrac{d}{dt}(| {u}_m(t)|^2e^{c_0^2c_-t}) \leq \frac{1}{c_-}\| {f}(t)\|_{V’}^2 e^{c_0^2c_-t}.
\end{align*}
Integrating from \(0\) to \(T\) we get,
\begin{align*}
| {u}_m(t)|^2e^{c_0^2c_-t} \leq | {u}_m(0)|^2 + \dfrac{1}{c_-}\int_0^T \| {f}(t)\|^2_{V’}e^{c_0^2c_-t}dt,
\end{align*}
which means,
\begin{equation*}
| {u}_m(t)|^2 \leq e^{-c_0^2c_-T}| {u}_m(0)|^2 + \frac{1}{c_-}\int_0^T \| {f}(t)\|^2_{V’}dt,
\end{equation*}
then,
\begin{equation*}
| {u}_m(t)|^2\leq e^{-c_0^2 c_- T}| {u}_m(0)| + \frac{1}{c_-}\| {f}\|^2_{L^2(0,T;V’)}.
\end{equation*}
Denoting \(\theta = e^{-c_0^2 c_- T}\) and \(c= \frac{1}{c_-}\| {f}\|^2_{L^2(0,T;V’)},\) we can write
\[| {u}_m(t)|^2 \leq \theta| {u}_m(0)|^2 +c,\]
so,
\[| {\tau}_m( {v})|^2 \leq \theta | {v}|^2 + c, \hspace{0.2cm} \forall v \in V_m.\]
Now, how \(0 < \theta < 1\) then \(0 < 1-\theta0\), big enough that \(c< (1-\theta)\rho^2_0.\) So if \(| {v}| < \rho_0\) then,
\[\theta | {v}|^2 + c \leq \theta \rho_0^2 + (1-\theta)\rho^2_0 = \rho^2_0\,,\]
where,
\[| {\tau}_m( {v})|^2 \leq \rho_0^2, \hspace{0.2cm}\forall m \in \mathbb{N},\]
which proves this lemma.
Lemma 4. The application \( {\tau}_m: V_m \mapsto V_m\) defined in (15) is continuous.
Proof. Let \( {v}_1, {v}_2 \in V_m\) and \( {u}_m, {z}_m\) solutions of the approximate problem with initial data \( {v}_1\) and \( {v}_2\), respectively. Our goal is to show that the solutions are Lipschitz-continuous, \(| {\tau}_m( {v}_1) – {\tau}_m( {v}_2)|\leq c_m | {v}_1- {v}_2|\) for some \(c_m>0.\)
\begin{align*}
&( {u}’_m(t), {w}_j) + c(l(u_1),l(u_2))(( {u}_m(t), {w}_j)) + b( {u}_m(t), {u}_m(t), {w}_j) = \langle {f}(t), {w}_j \rangle,\\
&( {z’}_m(t), {w}_j) + c(l(z_1),l(z_2))(( {z}_m(t), {w}_j)) + b( {z}_m(t), {z}_m(t), {w}_j) = \langle {f}(t), {w}_j \rangle.
\end{align*}
Doing the difference between these equations and defining \( {\eta}_m = {z}_m – {u}_m,\)
\begin{align*}
( {\eta}_m, {w}_j) &+ c(l(u_1),l(u_2))(( {u}_m(t), {w}_j)) -a(l(z_1),l(z_2))(( {z}_m(t), {w}_j))\\
&+b( {u}_m(t), {u}_m(t), {w}_j)-b( {z}_m(t), {z}_m(t), {w}_j) =0,
\end{align*}
we proceed as in [1]
\begin{equation*}
\frac{d | {\eta}_m|^2}{dt} – | {\eta}_m|^2\left(\frac{2}{c_-}\|u_{2m}(t)\|^2 + \frac{K^2}{c_-}\| {z}_m \|^2\right)\leq 0.
\end{equation*}
Defining \(\theta_m(t) = \left(\frac{2}{c_-}\|u_{2m}(t)\|^2 + \frac{K^2}{c_-}\| {z}_m \|^2\right)\), we get,
\begin{equation*}
\frac{d | {\eta}_m|^2}{dt} – | {\eta}_m|^2\theta_m(t)\leq 0.
\end{equation*}
Multiplying both sides of inequality \(e^{-\int_0^t \theta_m(s)ds},\)
\begin{equation*}
\dfrac{d}{dt}\left( | {\eta}_m(t)|^2e^{-\int_0^t \theta_m(s)ds} \right) \leq 0.
\end{equation*}
Integrating the inequality from \(0\) to \(T\),
\[| {\eta}_m(T)|^2e^{-\int_0^t \theta_m(s)ds} – | {\eta}_m(0)|^2\leq 0.\]
Defining \(c_m = e^{-\int_0^t \theta_m(s)ds},\)
\[| {\eta}_m(T)|^2 \leq c_m| {\eta}_m(0)|^2.\]
By other hand,
\[ {\eta}_m(s) = {u}_m(s) – {z}_m(s),\]
so,
\[| {u}_m(T)- {z}_m(T)|^2 \leq c_m | {u}_m(0) – {z}_m(0)|^2.\]
Then,
\[| {\tau}_m( {v}_1) – {\tau}_m( {v}_2)| \leq c_m| {v}_1 – {v}_2|,\]
which is what we want to prove.
The hypotheses of Brouwer’s fixed point theorem are satisfied by virtue of Lemmas 3 and 4, so we have
\[ {\tau}_m: \overline{B_{\rho_0}(0)} \rightarrow \overline{B_{\rho_0}(0)},\]
admits a fixed point, which means, there is a \( {v} \in \overline{B_{\rho_0}(0)}\) such as \( {\tau}_m( {v}) = {v},\) so, \( {u}_m(0) = {u}_m(T).\)
Then, for each \(m\in \mathbb{N},\) there is a least one \( {u}_m(t)\) such as \( {u}_m(0) \in \overline{B_{\rho_0}(0)}\) and, \(\forall j=1,\cdots,m\),
\begin{align*}
\begin{cases}
( {u}’_m(t), {w}_j)+ c(l(u_1),l(u_2)) (( {u}_m(t), {w}_j)) + b( {u}_m(t), {u}_m(t), {w}_j) = \langle {f}(t), {w}_j \rangle, \\
{u}_m(0) = {u}_m(t).
\end{cases}
\end{align*}
From the fact that \( {u}_m(0) \in \overline{B_{\rho_0}(0)}\) we can repeat the estimates getting a subsequence \(( {u}_{\nu})\) of \(( {u}_m)\) such as
\begin{align}
& {u}_\nu \stackrel{\star}\rightarrow {u} \mbox{ in } L^\infty(0,T;H), \label{convu1}\\
\end{align}
(19)
\begin{align}
& {u}_\nu \rightarrow {u} \mbox{ in } L^2(0,T;V),\\
\end{align}
(20)
\begin{align}
& {u}’_\nu \rightarrow {u}’ \mbox{ in } L^2(0,T;V’) \label{convu2}.
\end{align}
(21)
From the convergence results (19) – (21), by passing the limit in the approximate equation desired in (14). Similarly to the proof of the initial condition in the previous case, we prove that \( {u}(0) = {u}(T), \) which concludes the statement.
5. Conclusions
We studied the Navier-Stokes equations with non-local viscosity, considering a bounded domain \(\Omega \subset \mathbb{R}^2\) with smooth boundary \(\partial \Omega\). Using Faedo-Galerkin’s method and Brouwer’s fixed point theorem, we proved the strong solutions and periodic solutions.
Author Contributions:
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Conflicts of Interest:
”The authors declare no conflict of interest.”
Data Availability:
All data required for this research is included within this paper.
Funding Information:
The third author was supported by FCT – Fundação para a Ciência e a Tecnologia, through Centro de Matemática e Aplicações – Universidade da Beira Interior, under the Grant Number UI/BD/150794/2020, and also supported by MCTES, FSE and UE.