Cytomegalovirus dynamics model with random behavior

1.   Introduction

The presence of singularities and degeneracies in elliptic equations introduces significant challenges in analyzing the behavior of solutions. These singularities, especially near the origin or boundary, can profoundly affect the properties of the operator, making solutions more sensitive to changes in the domain. For instance, when $1 < p < N$, it is known that $\tilde{u}/|y| \in L^{p}(\mathbb{R}^{N})$ if $\tilde{u} \in W^{1, p}(\mathbb{R}^{N})$, or $\tilde{u}/|y| \in L^{p}(\Omega)$ when $\tilde{u} \in W^{1, p}(\Omega)$, where $\Omega$ is a bounded domain (see Lemma 2.1 in ^[12] for further details). In this context, the solution under consideration is $\tilde{u}$, and such behavior leads to the development of Hardy-type inequalities, which are crucial for controlling the singularities of solutions near critical points, particularly when the equation includes singular potential terms (see, e.g., ^{[1,12,17,18,20]}).

Furthermore, the presence of an indefinite weight in the source term creates several challenges, mainly because it can change sign or behave irregularly. This complicates the application of standard methods for proving the existence of solutions, such as ensuring the necessary properties of the energy functional. The irregular behavior of the weight also makes it difficult to use common mathematical tools like Sobolev embeddings and variational methods. To overcome these difficulties, this manuscript employs a more flexible approach based on critical point theory ^[4], which allows establishing the existence of solutions despite the complexities introduced by the indefinite weight.

Finally, the degeneracy of differential operators, such as $p$-Laplacian or $p(x)$-Laplacian, when coupled with a weight function $\omega(x)$ inside the divergence, introduces additional complexity to the problem. The presence of $\omega(x)$, whether it is singular or merely bounded, requires a shift in the selection of appropriate functional spaces. Traditional Sobolev spaces like $W^{1, p}(\Omega)$ or $W^{1, p(x)}(\Omega)$ may no longer be adequate in such cases. To properly handle the singularities or degeneracies, it becomes necessary to consider alternative Sobolev spaces, such as $W^{1, p(x)}(\omega, \Omega)$ (see section 2 for the definition of $W^{1, p}(\omega, \Omega)$), which are specifically designed to accommodate the weight function (see ^[6] for further details). The most recent contribution to the study of the $p$ Laplacian in a bounded domain and in the whole space can be found in respectively in ^[5] and ^[3], furthermore, the degenerate $p$-Laplacian operator combined with a Hardy potential can be found in ^[16].

This paper tackles the challenges posed by degeneracy, Hardy-type singularities, and sign-changing source terms, which are common in applied mathematical models, by examining a class of weighted quasilinear elliptic Dirichlet problem involving a variable exponent $p(x)$ and an indefinite source term. The main objective is to prove the existence of three weak solutions, using a critical point theorem introduced by Bonanno and Moranno in ^[4] while accounting for the complexities introduced by the operator's degeneracy and the singularities in the equation.

This manuscript explores the multiplicity of weak solutions to a weighted elliptic equations of the form:

where $\lambda$ is a positive parameter, $1 < q < N$, and $\Omega \subset \mathbb{R}^N$ (with $N \geq 2$) is a bounded open subset with smooth boundary $\partial \Omega$. The function $u$ is a solution to a weighted quasilinear elliptic equation involving a variable exponent $p(x)\in C_+(\overline\Omega)$(see, the beginning of Section 2) and the nonlinear source term of the form $k(x) |u|^{s(x)-2} u$ which involves a weight function $k(x)$ and may exhibit singularities on $\Omega$ and can change sign, belongs to a nonstandard Lebesgue space $L^{\gamma(x)}(\Omega)$.

The operator $\Delta_{p(x), a(x, u)} u$ represents a nonlinear generalization of the classical Laplacian, defined by:

here $a(x, u)$ denotes a Carathéodory function satisfying the inequality:

with $a_1, a_2$ are two positive constants, the function $\omega(x)$ is assumed to belongs to the local Lebesgue space $L_{\text{loc}}^1(\Omega)$, and it satisfies additional growth conditions, such as $\omega^{-h(x)} \in L^1(\Omega)$, where $h(x)$ satisfies certain bounds related to the variable exponent $p(x)$. Specifically, we assume that

The nonlinearity in the equation involves the functions $k(x)$ and $s(x)$, which are assumed to satisfy the following inequality for almost every $x \in \Omega$

where $p_h(x) = \frac{h(x) p(x)}{h(x) + 1}$.

2.   Backgrounds and preliminary results

Set, $\mathcal{S}(\Omega)$, the space that contains all measurable functions in $\Omega$ and

For $\tau > 0$, and $p(x)\in C_+(\overline\Omega),$ we use the following notations

In the sequel, we define the space $L^{p(x)}(\omega, \Omega)$ as follows

where $p(x)$ is a variable exponent, and $\omega(x)$ is a weight function. The space is endowed with a Luxemburg-type norm, given by:

Next, we define the corresponding variable exponent Sobolev space, which incorporates the variable exponent $p(x)$ in the functional setting.

with the norm

where $\|\nabla u\|_{p(x)} = \||\nabla u| \|_{p(x)}, |\nabla u| = (\sum\limits_{i = 1}^N \big|\frac{\partial u}{\partial x_{i}}\big|^{2})^{\frac{1}{2}}, \nabla u = \Big(\frac{\partial u}{\partial x_{1}}, \frac{\partial u}{\partial x_{2}}, ..., \frac{\partial u}{\partial x_{N}}\Big)$ is the gradient of $u$ at $(x_{1}, x_{2}, ..., x_{N}).$

Denote, by

the weighted Sobolev space and by $W^{1, p(x)}_0(\omega, \Omega)$ as the closure of $C_{0}^{\infty}(\Omega)$ in the space $W^{1, p(x)}(\omega, \Omega)$ endowed with the norm

Lemma 2.1. ^[8] If $p_1(x), p_2(x)\in C_+(\overline\Omega)$ such that $p_1(x) \leq p_2(x)$ a.e. $x\in \Omega,$ then there exists the continuous embedding $W^{1, p_{2}(x)}(\Omega)\hookrightarrow W^{1, p_{1}(x)}(\Omega)$.

Proposition 2.1 (^[9]) For $p(x)\in C_+(\overline\Omega), u, u_n \in L^{p(x)}(\Omega),$ one has

Let $0 < d(x) \in S(\Omega)$, and define the space

where $p(x)$ is a variable exponent, and $d(x)$ is a weight function. The space is equipped with a Luxemburg-type norm, defined by

Proposition 2.2 (^[10]) If $p\in C_+(\overline\Omega).$ Then

for every $u\in L^{p(x)}_{d(x)}(\Omega)$ and for a.e. $x\in \Omega$.

Combining Proposition 2.1 with Proposition 2.2, one has

Lemma 2.2. Let

For $p\in C_+(\overline\Omega), u \in W^{1, p(x) }(\omega, \Omega),$ we have

From Proposition 2.4 of ^[20], if $(\omega)$ holds, $W^{1, p(x)}(\omega, \Omega)$ is a reflexive separable Banach space.

From Theorem 2.11 of ^[15], if $(\omega)$ holds, the following embedding

is continuous, where

Combining (2.1) with Proposition 2.7 and Proposition 2.8 in ^[11], we get the following embedding

is continuous, where

Furthermore, the following embedding

is compact, when $1\leq t(x) < p_{h}^{*}(x).$

In what follows, and for any $p(x)\in C_+(\overline\Omega)$, let us denote by $p'(x): = \frac{p(x)}{p(x)-1}$, the conjugate exponent of $p(x)$.

Remark 2.1. Under Condition $(k)$, one has

● $1 < \beta(x) < p^{*}_{h}(x)$ for almost every $x\in\Omega$, where $\beta(x): = \frac{\gamma(x) s(x)}{\gamma(x)-s(x)}$, consequently

is compact.

● $1 < \alpha(x) < p^{*}_{h}(x)$ for almost every $x\in\Omega$, where $\alpha(x) = \gamma'(x)s(x)$, consequently

is compact.

Lemma 2.3 (Hölder type inequality ^[2,11]). Let $p_1, p_2, t\geq 1$ three functions that belong in $\mathcal{S}(\Omega)$ such that

If $f\in L^{p_1(x) }(\Omega)$ and $g\in L^{p_2(x) }(\Omega)$, then $fg\in L^{t(x) }(\Omega)$, moreover

Similarly, if $\frac{1}{t(x) }+\frac{1}{p_1(x) }+\frac{1}{p_2(x) } = 1$, for a.e. $x\in \Omega$, then

Lemma 2.4 (^[7]). Let $r_1(x)$ and $r_2(x)$ be measurable functions such that $r_1(x)\in L^{\infty}(\Omega)$, and $1\leq r_1(x)r_2(x)\leq\infty$, for a.e. $x\in\Omega$. Let $w\in L^{r_2(x)}(\Omega)$, $w\neq0$. Then

Let's define the functional $\mathcal{I}_{\lambda}\colon W^{1, p(x)}_0(\omega, \Omega)\to \mathbb{R}$ as

where

and

It is noted that, based on Remark 2.1 and Lemma 2.4, the aforementioned functionals are both well-defined and continuously Gâteaux differentiable (see ^[14] for further details). The Gâteaux derivatives are as follows

and

Furthermore, $\mathcal{M}'(u)$ is compact in the dual space $(W^{1, p(x)}_0(\omega, \Omega))^*$ (see ^[14]).

$u\in W^{1, p(x)}_0(\omega, \Omega)$ is said to be a weak solution of the problem $(1.1)$ if, the following holds for every $v\in W^{1, p(x)}_0(\omega, \Omega)$.

Lemma 2.5. $\mathcal{L}'$ is a strictly monotone coercive functional that belongs in $(W^{1, p(x)}_{0}(\omega, \Omega))^*.$

Proof. For any $u \in W^{1, p(x)}_{0}(\omega, \Omega)\setminus {\{0\}}$, by Lemma 2.2, one has

thus

then $\mathcal{L}'$ is coercive in view of $p(x)\in C_+(\overline\Omega)$.

According to (2.2) of ^[13], for all $x, y \in \mathbb{R}^{N}$, there is a positive constant $C_{p}$ such that

and

where $\langle., .\rangle$ is the usual inner product in $\mathbb{R}^{N}.$ Thus, for any $u, v\in X$ satisfying $u\neq v,$ by standard arguments we can obtain

hence, one has $\mathcal{L}'$ is strictly monotone in $W^{1, p(x)}_0(\omega, \Omega)$. □

Lemma 2.6. The functional $\mathcal{L}'$ is a mapping of $(S_{+})$-type, i.e. if $u_{n}\rightharpoonup u$ in $W^{1, p(x)}_{0}(\omega, \Omega),$ and $\overline{\lim}_{n\rightarrow \infty}\langle \mathcal{L}'(u_{n})-\mathcal{L}'(u), u_{n}-u)\rangle\leq 0,$ then $u_{n}\rightarrow u$ in $W^{1, p(x)}_{0}(\omega, \Omega).$

Proof. Let $u_{n}\rightharpoonup u$ in $W^{1, p(x)}_{0}(\omega, \Omega),$ and $\overline{\lim}_{n\rightarrow \infty}\langle \mathcal{L}'(u_{n})-\mathcal{L}'(u), u_{n}-u\rangle\leq 0.$

Noting that $\mathcal{L}'$ is strictly monotone in $W^{1, p(x)}_{0}(\omega, \Omega),$ we have

while

thus we get

Further, by (1.2) one has

then $u_{n}\rightarrow u$ in $W^{1, p(x)}_{0}(\omega, \Omega)$ via Lemma 3.2 in ^[19]. □

Lemma 2.7. $\mathcal{L}'$ is an homeomorphism.

Proof. The strict monotonicity of $\mathcal{L}'$ implies that it is injective. Since $\mathcal{L}'$ is coercive, it is also surjective, and hence $\mathcal{L}'$ has an inverse mapping.

Next, we show that the inverse mapping $(\mathcal{L}')^{-1}$ is continuous.

Let $\tilde{f}_n, \tilde{f} \in (W^{1, p(x)}_0(\omega, \Omega))^*$ such that $\tilde{f}_n \to \tilde{f}$. We aim to prove that $(\mathcal{L}')^{-1}(\tilde{f}_n) \to (\mathcal{L}')^{-1}(\tilde{f})$.

Indeed, let $(\mathcal{L}')^{-1}(\tilde{f}_n) = u_n$ and $(\mathcal{L}')^{-1}(\tilde{f}) = u$, so that $\mathcal{L}'(u_n) = \tilde{f}_n$ and $\mathcal{L}'(u) = \tilde{f}$. By the coercivity of $\mathcal{L}'$, the sequence $u_n$ is bounded. Without loss of generality, assume $u_n \rightharpoonup u_0$, which implies

Thus, $u_n \to u_0$ because $\mathcal{L}'$ is of $(S_+)$-type, which ensures that $\mathcal{L}'(u_n) \to \mathcal{L}'(u_0)$. Combining this with $\mathcal{L}'(u_n) \to \mathcal{L}'(u)$, we deduce that $\mathcal{L}'(u) = \mathcal{L}'(u_0)$. Since $\mathcal{L}'$ is injective, it follows that $u = u_0$, and hence $u_n \to u$. Therefore, we have $(\mathcal{L}')^{-1}(\tilde{f}_n) \to (\mathcal{L}')^{-1}(\tilde{f})$, proving that $(\mathcal{L}')^{-1}$ is continuous.

□

The following critical point theorems constitute the principal tools used to obtain our result.

Theorem 2.1. ([4, Theorem 3.6]). Let $X$ be a reflexive real Banach space and assume the following

● $\mathcal{L}: X \to \mathbb{R}$ be a coercive functional that is continuously Gateaux differentiable and weakly lower semicontinuous in the sequential sense

● The Gateaux derivative of $\mathcal{L}$ has a continuous inverse on the dual space $X^*$.

● $\mathcal{M}: X \to \mathbb{R}$ is a continuously Gateaux differentiable functional whith a compact Gateaux derivative.

Furthermore, suppose that

There exist a positive constant $d$ and a point $\overline{v} \in X$ such that $d 06 \mathcal{L}(\overline{v})$, and the following conditions are satisfied:

Then, for any $\lambda \in \Lambda_d$, $\mathcal{L} - \lambda \mathcal{M}$ has at least three distinct critical points in $X$.

3.   Main results

In this section, a theorem about the existence of at least three weak solutions to the problem (1.1) is obtained.

Recall the Hardy inequality (refer to Lemma 2.1 in ^[12] for more details), which asserts that for $1 < t < N$, the following inequality holds:

where the optimal constant $\mathcal{H}$ is given by:

By combining this with Lemma 2.1 and using the fact that $1 < q < p_h(x) < N$, we deduce the continuous embeddings

which leads to the inequality

where $\mathcal{H} = \left(\frac{N-q}{q} \right)^q$.

We are now ready to present our primary result. To this end, we define

for each $x \in \Omega$, here $B(x, \tilde{\mathfrak{D}})$ denotes a ball centered at $x$ with radius $\tilde{\mathfrak{D}}$. It is clear that there exists a point $x^0 \in \Omega$ such that $B(x^0, R) \subseteq \Omega$, where

In the remainder, assume that $k(x)$, fulfill this requirement

where $k_0$ is a positive constant, the symbol $\tilde{m}$ will represent the constant

with $\Gamma$ denoting the Gamma function.

Theorem 3.1. Assume that $p^- > s^+$, and, there exist two positive constants $d$ and $\delta > 0$, such that

and

then for any $\lambda \in] A_{\delta}, B_{d}[$, problem (1.1) has at least three weak solutions.

Proof. It is worth noting that the functional $\mathcal{L}$ and $\mathcal{M}$ associated with problem (1.1) and defined in (2.2) and (2.3), satisfy the regularity assumptions outlined in Theorem 2.1. We will now establish the fulfillment of conditions $(a_1)$ and $(a_2)$. To this end, let's consider

and consider $v_d \in X$ such that

Then, by the definition of $\mathcal{L}$, we have

Therefore, $\mathcal{L}(v_{\delta}) > d$. However, it is important to consider the following

In addition, for each $u\in\mathcal{L}^{-1}(]-\infty, d])$, we have

therefore,

Furthermore, we can deduce using Lemmas 2.3, 2.4 and Remark 2.1 the following

where $c_{\gamma' s}$ is the constant from the continuous embedding of $W^{1, p(x)}_0(\omega, \Omega)$ into $W^{1, \gamma'(x) s(x)}(\Omega)$.

This leads to the following result

and

Furthermore, we can establish the coerciveness of $\mathcal{I}_{\lambda}$ for any positive value of $\lambda$ by employing inequality $(3.1)$ once more. This yields the following result

When $\|u\|$ is great enough, the following can be inferred

By considering the fact that $p^- > s^+$, we can reach the desired conclusion. In conclusion, considering the aforementioned fact that

since all assumptions of Theorem 2.1 are fulfilled, it can be deduced that for any $\lambda \in \bar{\Lambda}_d$, the function $\mathcal{L}-\lambda \mathcal{M}$ possesses at least three critical points that belong in $X: = W^{1, p}_0(\omega, \Omega)$. Consequently these critical points are exactly weak solutions of problem $(1.1)$. □

Author contributions

Khaled Kefi: Conceptualization, Methodology, Writing–original draft, Supervision; Nasser S. Albalawi: Conceptualization, Methodology, Writing–original draft, Supervision. All authors have read and agreed to the published version of the manuscript.

Use of Generative-AI tools declaration

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

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at Northern Border University, Arar, KSA for funding this research work through the project number NBU-FPEJ-2025-1706-01.

Conflict of interest

The authors declare that they have no conflicts of interest.