This article is concerned with a mathematical model of tumor growth governed by 2 nd order diffusion equation
. The source of mitotic inhibitor is almost periodic and time-dependent within the tissue. The system is set up with the initial condition
C(
r, 0) =
C
0(
r) and Robin type inhomogeneous boundary condition
. Under certain conditions we show that there exists a unique solution for this model which is almost periodic.
Since 1980, researchers in mathematics and biology have proposed and studied several deterministic mathematical models for tumor growth by diffusion equation. The main assumption is that these models serve as simplified but complementary description for one aspect of a complex biological phenomenon: the growth and stability of tissue. Given some simplified conditions, the study of these models is focused on describing qualitatively the early stages of growth of tissue. The diffusion equation is the type of a linear partial differential equation.
For example, Shymko, Glass [
∂ C ∂ t = D ∂ 2 C ∂ x 2 − λ C + P S ( x ) (1)
with the source term S ( x ) = 1 − 2 L | x | for | x | ≤ L 2 and 0 otherwise. Here L is the initial length of the chalone-producing tissue, being confined to the domain | x | ≤ L / 2 . The production rate of the chalone is P > 0 per unit length, the diffusion coefficient is D and the decay rate is λ which is proportional to its concentration C ( x , t ) . Since S ( x ) is not constant, the source of mitotic inhibitor is not uniformly distributed within the tissue (in contrast to many earlier results). It was found that stable and unstable regimes of growth become significantly modified from the uniform-source case. Consequently, this model is very sensitive to the type of source term assumed. In fact, many of the existing deterministic models of tumor cell growth are proposed by an ordinary differential equation (ODE) coupled to one or more equations of reaction and diffusion type. The ODE derives from mass conservation applied to the tumor and describes the evolution of the tumor boundary and the reaction-diffusion equations describe the distribution of nutrients (oxygen and glucose) and growth inhibitory factors (chalones) [
Clearly, a more realistic model requires a higher dimension because systems governing tumor growth are best served in a three dimensional domain. Consequently, Britton and Chaplain studied a more generalized system below [
∂ C ∂ t = f ( C ) + D ∇ 2 C + λ S ( r ) in Ω × ( 0, ∞ ) , (2)
D ∂ C ∂ n + P C = 0 in ∂ Ω × ( 0, ∞ ) , P ≥ 0 (3)
C ( r ,0 ) = C 0 ( r ) ≥ 0 for r ∈ Ω (4)
where C = C ( r , t ) is the concentration of some chemical inhibitors in a bounded n-dimensional region Ω ( n = 1 , 2 , 3 ). We note that D , λ > 0 , S ≥ 0 , P is the permeability of the tissue surface. Using maximum principles for parabolic and elliptic operators, the authors examined the effect of growth inhibitory factor. It was shown that if f , C 0 and Ω satisfy the conditions of the parabolic comparison theorem, then C is always non-negative and unique. Also, the concentration decreases monotonically in the open interval ( 0, R ) provided that S ′ ( r ) ≤ 0 and f differentiable with f ′ < 0 . This model is certainly a big improvement compared to one-dimensional system.
The following model adds a more general, time-dependent source function S ( r , t )
∂ C ∂ t = f ( C ) + D ∇ 2 C + λ S ( r , t ) in Ω × ( 0, ∞ ) , (5)
D ∂ C ∂ n + P C = 0 in ∂ Ω × ( 0 , ∞ ) , P ≥ 0 , (6)
C ( r ,0 ) = C 0 ( r ) ∈ L 2 ( Ω ) . (7)
Here the source function is time-dependent. The boundary condition is of homogeneous Robin type. Existence and uniqueness of an almost periodic solution for this model were studied in [
In this paper, we study a similar system similar to (5)-(7) under non-homogeneous Robin boundary condition D ∂ C ∂ n + P C = K ( r , t ) in ∂ Ω × ( 0, ∞ ) and initial condition C ( r , 0 ) = C 0 ( r ) ∈ L 2 ( Ω ) . We prove that, under certain conditions of initial and boundary data, there exists a unique solution which is almost periodic. The existence is obtained via continuous contraction semigroup and fixed point theorem and uniqueness is obtained via integral estimates on L2 norm of C.
In this section, we discuss the existence and uniqueness of almost periodic solution when the source function is almost periodic and the dynamics of the system. We consider the following system with a time-dependent source function S ( r , t ) and Robin inhomogeneous boundary data:
∂ C ∂ t = f ( C ) + D ∇ 2 C + λ S ( r , t ) in Ω × ( 0 , ∞ ) , (8)
D ∂ C ∂ n + P C = K ( r , t ) in ∂ Ω × ( 0, ∞ ) , (9)
C ( r ,0 ) = C 0 ( r ) ∈ L 2 ( Ω ) . (10)
Here D , P > 0 , K ( r , t ) ∈ C 2 ( Ω × [ 0, T ] ) for any T > 0 . First, we prove the following existence theorem for (8)-(10).
Theorem 2.1. Assume that C 0 ∈ L 2 ( Ω ) and f ′ ( C ) ≤ − α < 0 , P > 0 . Let f ( ⋅ ) : H 1 ( Ω ) → L 2 ( Ω ) be local Lipschitz continuous, L 2 ( Ω ) -norm of K is uniformly bounded in time, and S is temporally almost periodic with its L 2 ( Ω ) -norm uniformly bounded in time. Then there exists a unique almost periodic solution for (8)-(10) such that C ∈ L ∞ ( 0, T ; L 2 ) ∩ L 2 ( 0, T ; H 1 ( Ω ) ) for any T > 0 .
PROOF. We use a transformation C = G + H where H is a smooth function in Ω satisfying the boundary condition D ∂ H ∂ n + P H = K ( r , t ) on ∂ Ω . Then (8) is converted to
∂ G ∂ t = f ( G + H ) + D ∇ 2 G + λ S ( r , t ) + D ∇ 2 H − ∂ H ∂ t (11)
where G satisfies the homogeneous boundary condition D ∂ G ∂ n + P G = 0 . Similar to the system studied in [
Let ‖ S ‖ ∞ be the L ∞ norm of S and ‖ C ‖ , ‖ ∇ C ‖ be the L 2 norms of C and ∇ C . We differentiate ‖ C ‖ 2 and substitute (10) to get
1 2 d d t ‖ C ‖ 2 = ∫ Ω C ⋅ C t d x = ∫ Ω f ( C ) ⋅ C d x + D ∫ Ω ∇ C ⋅ C d x + λ ∫ Ω S C d x = ∫ Ω f ( C ) C d x − D ‖ ∇ C ‖ 2 + ∫ ∂ Ω K ( r , t ) C d S − P ∫ ∂ Ω | C | 2 d S + λ ∫ Ω S ( r , t ) C d x (12)
Since L 2 ( Ω ) -norm of K is uniformly bounded in time, there exists a number c 0 such that
∫ ∂ Ω K ( r , t ) C d S − P ∫ ∂ Ω | C | 2 d S ≤ 1 2 P ∫ ∂ Ω | K | 2 d S + P 2 ∫ ∂ Ω | C | 2 d S − P ∫ ∂ Ω | C | 2 d S d S ≤ c 0 (13)
Substitute (13) in (12) we have
1 2 d d t ‖ C ‖ 2 = ∫ Ω f ( C ) ⋅ C d x + D ∫ Ω ∇ C ⋅ C d x + ∫ ∂ Ω K ( r , t ) C d S − P ∫ ∂ Ω | C | 2 d S + λ ∫ Ω S C d x ≤ − α ‖ C ‖ 2 − D ‖ ∇ C ‖ 2 + c 0 + λ 2 2 α ‖ S ‖ 2 + α 2 ‖ C ‖ 2 ≤ − α 2 ‖ C ‖ 2 − D ‖ ∇ C ‖ 2 + c 0 + λ 2 2 α ‖ S ‖ 2 . (14)
This implies that
‖ C ‖ 2 ≤ e − α t ‖ C 0 ‖ 2 + ( 1 − e − α t ) ( 2 c 0 α + λ 2 α 2 ‖ S ‖ 2 ) ≤ ‖ C 0 ‖ 2 + ( 2 c 0 + λ 2 α 2 ‖ S ‖ 2 ) = M (15)
when t ≥ T for some sufficiently large T. Therefore, all solutions C enter the following bounded set in L 2 ( Ω )
B = { C : ‖ C ‖ ≤ M } . (16)
Suppose that C 1 and C 2 are two solutions with same initial value C 0 ( r ) and boundary value K ( r , t ) . Then C = C 1 − C 2 satisfies the following system
∂ C ∂ t = f ( C 1 ) − f ( C 2 ) + D ∇ 2 C (17)
D ∂ C ∂ n + P C = 0 in ∂ Ω × ( 0 , ∞ ) , P > 0 (18)
with C ( r ,0 ) = 0 .
A quick calculation shows
1 2 d d t ‖ C ‖ 2 = ∫ Ω f ′ ( η ) | C | 2 d x − D ‖ ∇ C ‖ 2 − P ∫ ∂ Ω | C | 2 d S ≤ − α ‖ C ‖ 2 − D ‖ ∇ C ‖ 2 (19)
This implies that
1 2 d d t ‖ C ‖ 2 + α ‖ C ‖ 2 ≤ 0 (20)
therefore we know C ≡ 0 and the solution is unique since C ( r ,0 ) = 0 .
On the other hand, if we assume that C 1 and C 2 are two solutions with C 01 , C 02 as initial values respectively and same boundary value. Then C = C 1 − C 2 satisfies (17) and (18) with initial value C ( r ,0 ) = C 01 ( r ) − C 02 ( r ) . It is easy to check that C satisfies (19) and (20). Therefore,
‖ C ‖ 2 ≤ e − 2 α t ‖ C 01 − C 02 ‖ 2 . (21)
Define the solution operator S t ,0 : L 2 → L 2 by S t , 0 C 0 = C ( t ) for t ≥ 0 , where C ( t ) is the solution of (8)-(10). By (21), S t ,0 possesses strong contraction property with absorbing sets (16).
Recall that a function φ : I R → X where ( X , d X ) is a metric space, is called almost periodic [
d X ( φ ( t + τ ) , φ ( t ) ) ≤ ε (22)
for all t ∈ I R and τ ∈ M ε . Almost periodic functions play an important role in the theory of nonautonomous dynamical systems.
By similar arguments in [
We study a mathematical model of tumor growth presented by a diffusion equation with appropriate initial and boundary conditions. The boundary value is of Robin type and is inhomogeneous. We show that there exists a unique solution that is almost periodic. The main method is to show that a pullback attractor defines a unique almost periodic solution for the system.
This research was supported by William R. Kenan Jr. Professorship, Staley Small Grant and Wellesley College Faculty Award.
The author declares no conflicts of interest regarding the publication of this paper.
Bu, C. (2022) Existence and Uniqueness of Almost Periodic Solution for a Mathematical Model of Tumor Growth. Journal of Applied Mathematics and Physics, 10, 1013-1018. https://doi.org/10.4236/jamp.2022.104068