Analysis of a free boundary problem for vascularized tumor growth with time delays and almost periodic nutrient supply

1.   Introduction

Convex analysis is a branch of mathematics that studies the properties of convex sets and convex functions. In convex analysis, a set is considered convex if for any two points in the set, the line segment connecting those two points is also contained within the set. Convex functions are functions where the line segment between any two points on the graph of the function is above the graph. Convex analysis is used in various fields such as optimization, economics, and engineering to model and solve problems related to optimization, fairness, and risk management. The main goal of convex analysis is to provide a tractable framework for the study of optimization problems where the objective and constraint functions are convex.

Definition 1.1. ^[1] A mapping $\hbar:[\mathfrak{a}, \mathfrak{b}]\rightarrow \Re$ is called to be convex, if

where $\varpi_{1}, \varpi_{2} \in [\mathfrak{a}, \mathfrak{b}]$ and $\chi \in [0, 1].$

The Hermite-Hadamard (H-H) inequality asserts that, if a mapping $\hbar:\mathcal I\subset{\Re}\to{\Re}$ is convex in $\mathcal I$ for $\varpi_{1}, \varpi_{2}\in \mathcal I$ and $\varpi_{2} > \varpi_{1}$, then

Interested readers can refer to ^[1].

Both inequalities hold in the reversed direction if $\hbar$ is concave. The inequality (1.2) is known in the literature as the Hermite-Hadamard's inequality.

We note that the Hermite-Hadamard's inequality may be regarded as a refinement of the concept of convexity and it follows easily from Jensen's inequality. The classical Hermite-Hadamard's inequality provides estimates of the mean value of a continuous convex function $\hbar : [\varpi_{1}, \varpi_{2}]\rightarrow \Re.$ Many aesthetic inequalities on convex functions exist in the literature, among which Jensen's inequality has a special place. This inequality, which is proved under fairly simple conditions, is extensively used by researchers in fields such as information theory as well as in inequality theory. Jensen's inequality is presented as follows.

Let $0 < \zeta_{1}\leq \zeta_{2}\leq...\leq \zeta_{n}$ and let $\rho = \big(\rho_{1}, \rho_{2}, ..., \rho_{n}\big)$ be non-negative weights such that $\sum_{j = 1}^{n}\ \rho_{j} = 1$. The famous Jensen inequality (see ^[2]) in the literature states that if $\hbar$ is convex functions on the interval $\big[\, \varpi_{1}, \varpi_{2}\big]$, then

for all $\zeta_{j}\in [\, \varpi_{1}, \varpi_{2}],$ $\rho _{j}\in \left[0, 1 \right]\, and\quad\left(j = 1, 2, ..., n\right).$

It is one of the key inequality that helps to extract bounds for useful distances in information theory.

Although the studies on Jensen's inequality are a topic that many researchers have focused on, the variant put forward by Mercer is the most interesting and remarkable among them. In 2003, a new variant of Jensen's inequality was introduced by Mercer ^[3] as:

If $\hbar$ is a convex function on $\left[\, \varpi_{1}, \varpi_{2}\right]$, then

holds for all $\zeta_{j}\in \left[\, \varpi_{1}, \varpi_{2}\right],$ $\rho _{j}\in \left[0, 1 \right]\, and\, \left(j = 1, 2, ..., n\right).$

Pe$\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{{\rm{c}}}$arić et al. proposed several generalizations on Jensen-Mercer operator inequalities ^[4]. Later on, Niezgoda ^[5] have provided the several extensions to higher dimensions for Mercer's type inequalities. Recently, Jensen-Mercer's type inequality has made significant contribution on inequality theory ue to its prominent characterizations. Kian ^[6] contemplated the idea of Jensen inequality for superquadratic functions.

In ^[7], the idea of the Jensen-Mercer inequality has been used by Kian and Moslehian and the following Hermite-Hadamard-Mercer's inequality was demonstrated:

where $\hbar$ is a convex function on $[\mathfrak{a}, \mathfrak{b}].$ For more recent studies linked to the Jensen-Mercer inequality, one can refer ^[8,9,10].

Although fractional analysis has a history that is as old as classical analysis historically, it has attracted more attention of researchers recently. Fractional analysis, with its use in real world problems, its contribution to engineering sciences, and its structure that is open to development and different dimensions that it brings to mathematical theories, is in a continuous effort to progress. The definition of fractional order derivatives and integrals and the contribution of each new operator to the field from a different direction is one aspect that keeps the fractional analysis up to date. When the new operators are examined carefully, various features such as singularity, locality and generalization and differences in the kernel structures stand out. Although generalization and inference are the cornerstones of mathematical methods, these different features in the new fractional operators add new features to problem solutions, especially the time memory effect. Accordingly, various operators such as Riemann-Liouville, Grünwald Letnikov, Raina, Katugampola, Prabhakar, Hilfer, Caputo-Fabirizio are some of the prominent operators that reveal this potential of fractional analysis. Now we will continue by introducing the Atangana-Baleanu integral operators, which have a special place among these operators.

Some researchers have become interested in the concept of a fractional derivative in recent years. There are two types of nonlocal fractional derivatives: those with singular kernels, such as the Riemann-Liouville and Caputo derivatives, and those with nonsingular kernels, such as the Caputo-Fabrizio and Atangana-Baleanu derivatives. Fractional derivative operators with non-singular kernels, on the other hand, are very effective at resolving non-locality in real-world problems. We'll go over the Caputo-Fabrizio integral operator later.

Definition 1.2. ^[11] Let $\hbar\in H^{1}(\mathfrak{a}, \mathfrak{b})$, (where $H^{1}$ is class of first order differentiable function) $\mathfrak{b} > \mathfrak{a}$, $\alpha\in[0, 1],$ then the definition of the new Caputoifractional derivative is:

where $M(\alpha)$ is a normalizationifunction.

Moreover, theicorresponding Caputo-Fabrizioifractional integral operator is given as:

Definition 1.3. ^[12] Let $\hbar\in H^{1}(\mathfrak{a}, \mathfrak{b})$, $\mathfrak{b} > \mathfrak{a}$, i$\alpha \in [0, 1]$.

and

where $M(\alpha)$ is a normalizationifunction.

The Atangana-Baleanu derivative is a type of fractional derivative that was introduced in 2016 by Atangana and Baleanu. It is a generalization of the classical derivative, which extends the concept of differentiation to functions with non-integer order derivatives. The Atangana-Baleanu derivative is defined using the concept of the Caputo fractional derivative and provides a better approximation for real-world phenomena that exhibit memory and hereditary properties. This derivative has found applications in various fields such as physics, engineering, and biology, where it has been used to model problems involving anomalous diffusion, viscoelasticity, and electrochemistry, among others. The Atangana-Baleanu derivative has also been found to be more effective than other fractional derivatives in solving boundary value problems, particularly those with variable order derivatives. Atangana-Baleanu defined the derivative operator in both Caputo and Riemann-Liouville terms.

Definition 1.4. ^[13] Let $\nu > \mu$, $\alpha \in \lbrack 0, 1]$ and $g\in H^{1}(\mu, \nu).$ The new fractional derivative is given:

Definition 1.5. ^[13] Let $g\in H^{1}(\mathfrak{a}, \mathfrak{b})$, $\mu > \nu$, $\alpha \in [0, 1]$. The new fractionaliderivative is given as

However, Atangana-Baleanu($AB$)-fractional integral operator as follows.

Definition 1.6. ^[13] The fractional integral operator with non-local kernel of a function $g\in H^{1}(\mu, \nu)$ is defined by

where $\mathfrak{b} > \mathfrak{a}, \alpha\in[0, 1].$

The right hand side of the $AB$-fractional integral operator is

Here, $\Gamma (\alpha)$ is the Gamma function. The positivity of the normalization function $M(\alpha)$ implies that the fractional $AB$-integral of a positive function is positive. It is worth noticing that the case when the order $\alpha\rightarrow 1$ yields the classical integral and the case when $\alpha\rightarrow 0$ provides the initial function.

Furthermore, in recent years, several significant inequalities for interval-valued functions (H-H, Ostrowski, Simpson, and others) have been studied. In ^[14,15], to obtain Ostrowski type inequalities, Chalco-Cano et al. used the Hukuhara derivative for interval-valued functions. For more results, we can see the literature in ^{[16,17,18,19]}.

2.   Calculus on interval and inequalities

This section contains notation for interva analysis as well as the background information. The space of all closed intervals of $\Re$ is denoted by $\mathcal I_{c}$ and $\Delta$ is a bounded element on $\mathcal I_{c}$. As we have

where $\underline{\alpha}, \bar{\alpha} \in \Re$ and $\underline{\alpha} \leq\bar{\alpha}.$ $L(\Delta) = \bar{\alpha}-\underline{\alpha}$ can be used to describe the length of the interval $\Delta = [\underline{\alpha}, \bar{\alpha}].$ The endpoint of left and right of interval $\Delta$ are denoted by the numbers $\underline{\alpha}$ and $\bar{\alpha}$, respectively. The interval $\Delta$ is said to be degenerate where $\underline{\alpha} = \bar{\alpha}$ and the form $\Delta = \alpha = [\underline{\alpha}, \bar{\alpha}]$ is used. Also, if $\underline{\alpha} > 0,$ we can say $\Delta$ is positive and if $\bar{\alpha} < 0,$ we can say that $\Delta$ is negative. Since $\mathcal I^+_{c}$ and $\mathcal I^-_{c}$ denote the sets of all closed positive intervals and closed negative intervals of $\Re,$ respectively, the Pompeiu Hausdroff distance between the $\Delta$ and $\Lambda$ is defined by

The set $(\mathcal I_{c}, \mathcal P_{G})$ is a complete metric space. As far as we know (see ^[20]), the absolute vale of $\Delta$ is $|\Delta|,$ which is the maximum values of the end points:

The following are the concepts for fundamental interval arithmetic operations for the intervals $\Delta$ and $\Lambda$.

The scalar multiplication of the interval $\Delta$ is defined by

where $\gamma \in \Re.$ The interval $\Delta$ in an opposite form is

where $\gamma = -1.$ In general, $- \Delta$ is not additive inverse for $\Delta,$ i.e., $\Delta-\Delta \neq 0.$

Definition 2.1. ^[21] For some kind of intervals $\Delta, \Lambda \in \mathcal I_{c},$ we denote the $G$-difference of $\Delta, \Lambda$ as the $\Omega\in\mathcal I_{c},$ we have

It appears beyond controversy that

where $L(\Delta) = \bar{\alpha}-\underline{\alpha}$ and $L(\Lambda) = \bar{\beta}-\underline{\beta}$. A large number of algebraic properties are produced by the operation definitions, allowing $\mathcal I_{c}$ to be a space of quasi-linear (see ^[22]).

The following are some of these characteristics (see ^[23,24,25]).

$(1)$ (Law of associative under +) $(\Delta+\Lambda)+\Theta = \Delta+(\Lambda+\Theta);$ for all $\Delta, \Lambda, \Theta \in \mathcal I_{c}.$

$(2)$ (Additive element) $\Delta+0 = 0+\Delta = \Delta;$ for all $\Delta \in \mathcal I_{c}.$

$(3)$ (Law of commutative under +) $\Delta+\Lambda = \Lambda+\Delta;$ for all $\Delta, \Lambda \in \mathcal I_{c}.$

$(4)$ (Law of cancelation under +) $\Delta+\Theta = \Lambda+\Theta \Rightarrow \Delta = \Lambda;$ for all $\Delta, \Lambda, \Theta \in \mathcal I_{c}.$

$(5)$ (Law of associative under $\times$) $(\Delta. \Lambda).\Theta = \Delta.(\Lambda.\Theta);$ for all $\Delta, \Lambda, \Theta \in \mathcal I_{c}.$

$(6)$ (Law of commutative under $\times$) $\Delta.\Lambda = \Lambda.\Delta;$ for all $\Delta, \Lambda \in \mathcal I_{c}.$

$(7)$ (Multiplicative element) $\Delta.1 = 1.\Delta = \Delta;$ for all $\Delta \in \mathcal I_{c}.$

$(8)$ (The first law of distributivity) $\lambda(\Delta+\Lambda) = \lambda\Delta+ \lambda\Lambda;$ for all $\Delta, \Lambda \in \mathcal I_{c}, \lambda \in \Re.$

$(9)$ (The second law of distributivity) $(\lambda+\gamma)\Delta = \lambda\Delta+ \gamma\Delta;$ for all $\Delta \in \mathcal I_{c}, for \ all \ \lambda, \gamma \in \Re.$

Apart from these features, the distributive law always applies to intervals. For example,

$\Delta = [1, 2],$ $\Lambda = [2, 3]$ and $\Theta = [-2, -1].$ We have

where

The inclusion of $\subseteq,$ is another distance feature, which is desired by

Moore defined the Riemann integral for functions with interval values in ^[24]. $IR_{([\varpi_{1}, \varpi_{2}])}$ denotes the sets of Riemann-Liouville integrable interval-valued functions and $R_{([\varpi_{1}, \varpi_{2}])}$ denotes the real-valued functions on $[\varpi_{1}, \varpi_{2}]$, respectively. The following theorem establishes a connection between Riemann integrable (R)functions and (R)-integrable functions (see ^[26], pp. 131).

Theorem 2.1. A mapping of interval-valued $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I_{c}$ with $\hbar(\chi) = [\underline{\hbar}(\chi), \bar{\hbar}(\chi)].$ The mapping $\hbar \in IR_{([\varpi_{1}, \varpi_{2}]}$ $\Leftrightarrow$ $\underline{\hbar}(\chi), \bar{\hbar}(\chi) \in R_{([\varpi_{1}, \varpi_{2}]}$ and

In ^[19,27], Zhao et al. defined the following convex interval-valued function.

Definition 2.2. For all $x, \upsilon \in [\varpi_{1}, \varpi_{2}]$ and $v\in (0, 1),$ a mapping $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I^+_{c}$ is $h$-convex stated as

where $h: [c, d]\rightarrow \Re$ is a nonnegative mapping with $h\neq0$ and $(0, 1)\subseteq [c, d].$ We shall show that the set of all $h$-convex interval-valued functions with $SX(h, [\varpi_{1}, \varpi_{2}], \mathcal I^+_{c})$.

The standard definition of a convex interval-valued function is (2.4) with $h(\nu) = \nu$ (see^[28]). If $h(\nu) = \nu^s$ in (2.4), then we get the definition of an $s$-convex interval-valued function (see ^[29]).

Zhao et al. used the h-convexity of interval-valued functions in ^[19] and obtained the following H-H inclusion:

Theorem 2.2. If $SX(h, [\varpi_{1}, \varpi_{2}], \mathcal I^+_{c})$ and $h(\frac{1}{2})\neq0,$ then we have

Remark 2.1. The inclusion 2.5 becomes the following if $h(\vartheta) = \vartheta$

which was discovered by Sadowska in ^[28].

The inclusion 2.5 becomes the following if $h(\vartheta) = \vartheta^s$

which was discovered by Osuna-Gómez et al. in ^[30].

Definition 2.3. ^[28] A mapping $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I_{c}$ is said to be convex interval-valued, if for all $x, v \in [\varpi_{1}, \varpi_{2}], \nu\in(0, 1),$ we have

Theorem 2.3. A mapping $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I_{c}$ is said to be convex interval-valued, if and only if $\underline{\hbar}$ is a convex function on $[\varpi_{1}, \varpi_{2}]$ and $\bar \hbar$ is a concave function on $[\varpi_{1}, \varpi_{2}].$

Theorem 2.4. (Jensen's Inclusion Interval Valued Function) ^[31] Let $0 < \zeta_{1}\leq \zeta_{2}\leq...\leq \zeta_{n}$ and $\hbar$ be a convex function on the interval-valued mapping on an interval containing $\rho_{k}$, then the following inclusion holds:

where $\sum_{j = 1}^{n}\ \rho_{j} = 1,$ $\rho _{j}\in \left[0, 1\right].$

The inequality (2.9) is extended to convex interval-valued functions by the authors in ^[31] as follows.

Theorem 2.5. (Jensen-Mercer Inclusion Interval Valued Function) ^[31] Let $\hbar$ be a convex function on the interval-valued mapping on $[\varpi_{1}, \varpi_{2}]$ such that $L(\varpi_{2})\geq L(\varpi_{o})$ for all $\varpi_{o}\in [\varpi_{1}, \varpi_{2}],$ then the following inclusion

is true.

3.   Main results

In this section, we prove the Mercer's Hermite-Hadamard type inclusion for convex interval-valued functions via Atangana-Baleanu fractional integrals.

Theorem 3.1. Let $\alpha \in (0, 1).$ Suppose that if $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I^+_{c}$ is an interval valued convex function such that $\hbar(\chi) = [\underline{\hbar}(\chi), \bar{\hbar}(\chi)]$ and $\mathcal L(\varpi_{2})\geq \mathcal L(\varpi_{0}),$ $\forall \varpi_{0} \in [\varpi_{1}, \varpi_{2}],$ then the following inclusions for $AB$(Atangana-Baleanu) fractiona integral operator hold:

where $\mathcal M\left(\alpha \right) > 0$ is a normalization function.

Proof. Since $\hbar$ is an interval valued convex function, then for all $u, v \in [\varpi_{1}, \varpi_{2}],$ we have

By Using

and

for all $\mathcal X, \mathcal V \in [\varpi _{1}, \varpi _{2}]$ and $\chi \in [0, 1]$, we get

Now, multiplying both sides of (3.4) by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

and the proof of the first inclusion in (3.1) is completed. To prove that the second inclusion in (3.1), first we note that $\hbar$ is an interval-valued convex function, so we have

and

Adding (3.7) and (3.8), we obtain the following from Jensen-Mercer inclusion.

Now, multiplying both sides of (3.9) by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

Concatenating the Eqs (3.6) and (3.10), we can get (3.1).

Corollary 3.1. By the same assumptions as defined in Theorem 3.1 with $\underline{\hbar}(\chi) = \bar{\hbar}(\chi)$, we have

Theorem 3.2. Let $\alpha \in (0, 1).$ Suppose that if $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I^+_{c}$ is an interval valued convex function such that $\hbar(\chi) = [\underline{\hbar}(\chi), \bar{\hbar}(\chi)]$ and $\mathcal L(\varpi_{2})\geq \mathcal L(\varpi_{0}),$ $\forall \varpi_{0} \in [\varpi_{1}, \varpi_{2}],$ then the following inclusions for the AB (Atangana-Baleanu) fractional integral operator hold:

where $\mathcal M\left(\alpha \right) > 0$ is a normalization function.

Proof. Since $\hbar$ is an interval valued convex function, we have for all $u, v \in [\varpi_{1}, \varpi_{2}],$

By Using

and

for all $\mathcal X, \mathcal V \in [\varpi _{1}, \varpi _{2}]$ and $\chi \in [0, 1]$, we get

Now, multiplying both sides of above by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

and the proof of the first inclusion (3.12) is completed. To prove the second inclusion in (3.12), first we note that $\hbar$ is an interval-valued convex function, we have

and

Adding above equations, we get

Multiplying both sides of above by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

Concatenating the Eqs (3.15) and (3.20), we can get (3.12). This ends the proof.

Corollary 3.2. By the same assumptions as defined in Theorem 3.2 with $\underline{\hbar}(\chi) = \bar{\hbar}(\chi)$, we have

Remark 3.1. If $\alpha = 1$ in Corollary 3.2, then it reduces to Kian and Moslehian in Theorem 2.1 ^[6].

Corollary 3.3. By the same assumptions as defined in Theorem 3.2 with $\mathcal X = \varpi _{1}$ and $\mathcal V = \varpi _{2}$, we have

Theorem 3.3. Suppose that if $\hbar:[\varpi_{1}, \varpi_{2}]\rightarrow \mathcal I^+_{c}$ is an interval valued convex function suchithat $\hbar(\chi) = [\underline{\hbar}(\chi), \bar{\hbar}(\chi)]$ and $\mathcal L(\varpi_{2})\geq \mathcal L(\varpi_{0}),$ $\forall \varpi_{0} \in [\varpi_{1}, \varpi_{2}],$ then the following inclusions for AB(Atangana-Baleanu) fractional integral operator hold:

Proof. Since $\hbar$ is an interval valued convex function, we have for all $u, v \in [\varpi_{1}, \varpi_{2}],$

and

for all $\mathcal X, \mathcal V \in [\varpi _{1}, \varpi _{2}]$ and $\chi \in [0, 1]$, we get

Now, multiplying both sides of above by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

The proof of the first inclusion in (3.28) is completed. To prove that the second inclusion (3.28) from Jensen-Mercer inclusion, we have

and

Adding above equations, we get

Multiplying both sides of (3.31) by $\frac{\alpha}{\mathcal M\left(\alpha \right)\Gamma\left(\alpha \right)}\chi^{\alpha-1}$ and integrating by inclusion with respect to $\chi$ over $[0, 1],$ we obtain

Concatenating the Eqs (3.28) and (3.32), we can get (3.23). This ends the proof.

Corollary 3.4. By the same assumptions as defined in Theorem 3.3 with $\underline{\hbar}(\chi) = \bar{\hbar}(\chi)$, we have

Corollary 3.5. By the same assumptions as defined in Theorem 3.3 with $\mathcal X = \varpi _{1}$ and $\mathcal V = \varpi _{2}$, we have

4.   Applications to matrix

We denote by $C^n$ the set of $n\times n$ complex matrices, $M_n$ the algebra of $n\times n$ complex matrices, and by $M^+_n$ the strictly positive matrices in $M_n.$ That is, $\mathcal A\in M^+_n$ if $\langle \mathcal Ax, x\rangle > 0$ for all nonzero $x\in C^n.$ In ^[32], Sababheh proved that the function $\psi \left(\theta \right) = \left\vert \left\vert \mathcal A^{\theta }\mathcal Y\mathcal B^{1-\theta }+\mathcal A^{1-\theta }\mathcal Y\mathcal B^{\theta }\right\vert \right\vert, \mathcal A, \mathcal B\ \in M_{n}^{+}\, \ \mathcal Y\in M_{n},$ is convex for all $\theta \in [0, 1].$

Example 4.1. By using Theorem 3.1 with satisfies the conditions in Theorem 3.1, we have

Example 4.2. Under the same assumptions of above example, if we consider Theorem 3.2, we have

Example 4.3. Under the same assumptions of above example, if we consider Theorem 3.3, we have

5.   Conclusions

The use of Atangana-Baleanu fractional integrals allowed us to demonstrate Hermite-Hadamard-Mercer inclusions for convex interval-valued functions. It is a fascinating and original problem where prospective scholars can discover related inequality for different convexities and fractional integrals. It would be intriguing to apply these findings to other convexities. In this exciting area of disparities and analysis, we anticipate that much study will centre on our recently revealed approach. It is possible to broaden and develop the amazing techniques and wonderful ideas in this text to coordinates and fractional integrals. This is how we intend to carry out our research in the long run.

Data availability

The data used to support the findings of this study are available from the corresponding author upon request.

Authors' Contributions

All authors contribute equally in this paper.

Conflict of interest

The authors declare that they have no conflict of interest.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work grant code: 22UQU4331214DSR01.

Funding

This paper does not receive any external funding.