# An Existence Theorem for a Nonlocal Global Pandemic Model for Insect-Borne Diseases

Authors: John R. Cannon; Daniel J. Galiffa
Journal: International Journal of Differential Equations
 (2014)

---

## Abstract

We construct and analyze a nonlocal global pandemic model that comprises a system of two nonlocal integrodifferential equations (functional differential equations) and
an ordinary differential equation. This model was constructed by considering a spherical coordinate transformation of a previously established epidemiology model that can be applied to insect-borne
diseases, like yellow fever. This transformation amounts to a nonlocal boundary value problem on the unit sphere and therefore can be interpreted as a global pandemic model for insect-borne diseases. We ultimately show that a weak solution to the weak formulation of this model exists using a fixed point argument, which calls upon the construction of a weak formulation and the existence
and uniqueness of an auxiliary problem.


---

## Body

## 1. Introduction

In [1], Cannon and Galiffa analyzed the nonlocal model:
(1)
-
α
(
∫
0
1
u
(
t
)
d
t
)
u
′′
=
f
(
x
)
,
x
∈
(
0,1
)
u
(
0
)
=
a
,
u
(
1
)
=
b
,
with α
=
α
(
q
) being a positive function of q defined over -
∞
<
q
<
∞, f
(
x
) defined over 0
≤
x
≤
1, and a and b real constants. The analysis contained in [1] included both existence and uniqueness theorems for the analytic solution, u, to (1) and a discretization, as well as the construction and implementation of a finite-difference procedure. This model was a special case of [2], which has applications to thermodynamics.The following homogeneous extension of model (1) was considered in [3]:
(2)
-
α
(
∫
0
1
u
(
t
)
d
t
)
u
′′
(
x
)
+
[
u
(
x
)
]
2
n
+
1
=
0
,
x
∈
(
0,1
)
,
u
(
0
)
=
a
,
u
(
1
)
=
b
.
As in (1), the existence and uniqueness of the analytic solution, u, to (2), as well as a discretization, were established and another finite-difference procedure was developed and executed; however, the analysis was significantly different than in [1].Upon the completion of [3], the authors realized the potential for additional applications to physical problems via mathematical models consisting of variations of the nonlocal term α
(
∫
0
1
‍
u
(
t
)
d
t
). In particular, a model was constructed and analyzed in [4, 5] that can be applied to insect-borne diseases, for example, yellow fever, which is as follows:
(3)
u
t
=
u
x
x
-
α
(
∫
a
b
w
(
x
,
t
)
d
x
)
u
+
γ
w
,
hhhh
-
∞
<
x
<
∞
,
t
≤
T
,
v
t
=
v
x
x
+
α
(
∫
a
b
w
(
x
,
t
)
d
x
)
u
-
β
v
,
hhhh
-
∞
<
x
<
∞
,
t
≤
T
,
w
t
=
β
v
-
γ
w
,
-
∞
<
x
<
∞
,
t
≤
T
,
u
(
x
,
0
)
=
f
(
x
)
,
-
∞
<
x
<
∞
,
v
(
x
,
0
)
=
0
,
-
∞
<
x
<
∞
,
w
(
x
,
0
)
=
0
,
-
∞
<
x
<
∞
with β and γ arbitrary positive constants.We next give an overview of this model, which is based on [5]. Biologically, β represents the rate at which mobile infecteds become immobile and γ represents the rate at which immobile infecteds recover. The unknowns u, v, and w, respectively, represent a susceptible class, a mobile infected class, and an immobile infected class. In addition, neither infected class is infective; that is, susceptibles become infected via contact with vectors (disease-spreading agents, e.g.,Aedes aegypti mosquitoes), which was not modeled directly. The nonnegative function α
(
q
), 0
≤
q
<
∞, represents the per capita rate at which susceptibles contract the disease in a disease-free population, which is decreasing with q. The argument ∫
a
b
‍
w
(
x
,
t
)
d
x, representing the number of the population immobilized from the disease, when substituted into α brings about a reduction in the rate of infection via actions taken at large and at the individual level, to reduce the risk of infections. Even more specifically, ∫
a
b
‍
w
(
x
,
t
)
d
x counts the number of immobile infecteds within a single fixed interval [
a
,
b
] and influences the rate α
(
q
) from a biological standpoint. For example, if all news stories about immobile infecteds concern individuals living only in region [
a
,
b
], then susceptibles will take precautions based only on the reported number ∫
a
b
‍
w
(
x
,
t
)
d
x.In this paper, we consider an extension of model (3) via a spherical coordinate transformation that results in the following model on the unit sphere S:
(4)
u
t
=
1
sin
ϕ
∂
∂
ϕ
(
sin
ϕ
∂
u
∂
ϕ
)
+
csc
2
ϕ
∂
2
u
∂
θ
2
-
α
(
∫
Ω
‍
w
(
ϕ
,
θ
,
t
)
sin
ϕ
d
ϕ
d
θ
)
u
+
γ
w
,
hhhhhhhhhhhhhhhhh
on
⁡
S
×
(
0
,
T
]
,
v
t
=
1
sin
ϕ
∂
∂
ϕ
(
sin
ϕ
∂
w
∂
ϕ
)
+
csc
2
ϕ
∂
2
v
∂
θ
2
+
α
(
∫
Ω
‍
w
(
ϕ
,
θ
,
t
)
sin
ϕ
d
ϕ
d
θ
)
u
-
β
v
,
hhhhhhhhhhhhhhhhh
on
⁡
S
×
(
0
,
T
]
,
w
t
=
β
v
(
ϕ
,
θ
,
t
)
-
γ
w
(
ϕ
,
θ
,
t
)
,
on
⁡
S
×
(
0
,
T
]
,
u
(
ϕ
,
θ
,
0
)
=
f
(
ϕ
,
θ
)
on
⁡
S
,
v
(
ϕ
,
θ
,
0
)
=
0
on
⁡
S
,
w
(
ϕ
,
θ
,
0
)
=
0
on
⁡
Ω
⊂
S
.
Physically, the system comprising (4) can be interpreted as a preliminary model for a global pandemic of an insect-borne disease. First, we note that the system of spherical coordinates is the natural and historical coordinate system for the surface of the Earth. Second, we note that diffusion has been used to describe epidemics in the past; compare [6–9] for a few of several such examples. Third, the system proposed herein is preliminary, since we take no account of the oceans, as to do so will require boundary conditions to be assigned to the coasts of the Earth's land masses. Fourth, spherical coordinates come into play for relatively short distances on the Earth's surface. For example, a Cartesian coordinate system (
x
,
y
,
z
) with origin located at the John F. Kennedy Space Center will find the plane (
x
,
y
,
0
) approximately 20,000 feet over Tampa, FL. Fifth, we have made no attempt to model air travel, which can be handled by the introduction of spatially located sources and sinks. Sixth, we note that mosquitoes have been found alive in shipping containers filled with used tires shipped to the United States from Asia, as well as other locations, which provides a vector for tropical Asian diseases on the American continent (cf. http://www.epa.state.il.us/land/tires/mosquito-borne-illnesses.html, for one such example). Moreover, we also note that passengers entering Australia are sprayed with insecticide disinfectant before they can be allowed to deplane (http://abcnews.go.com/).Finally, our interest here is to show that a solution exists to (4) in order to enable the study of the possible effect of the death toll of a disease on the dynamics of the epidemic itself. Thus, the system (4), which incorporates a coefficient containing a functional of part of the solution of the system, requires a nonlinear analysis, which we will outline later in this introduction.Several papers have recently been published regarding the mathematical analysis of pandemics, as well as the construction of mathematical models comprising a system of differential equations applicable to such biological phenomena. Several of these papers influenced our current work and it is therefore worth briefly describing them before we further discuss (4).For example, in [10] the authors explain how the utilization of antiviral drugs was intended to moderate the severity of a new strain of an influenza pandemic and that the success of these drugs was dependent upon the timely onset of therapy (within 48 hours) after the appearance of clinical symptoms. They discuss further that this requirement can be understood by a compartmental model that examines the density of infected individuals in terms of the time elapsed since the onset of symptoms. Therefore, based on this compartmental model, a system of delay-differential equations, with both discrete and disturbed delays, was constructed and analyzed. This system is as follows:
(5)
d
d
t
S
(
t
)
=
-
β
S
(
t
)
Q
(
t
)
,
d
d
t
E
(
t
)
=
β
S
(
t
)
Q
(
t
)
-
μ
E
E
(
t
)
,
for t
≥
0, where susceptible and exposed classes are, respectively, denoted by S and E, with β representing the baseline transmission rate, 1
/
μ
E representing the incubation period, and Q
(
t
) representing the force of infection yet to be formulated.As a precaution to a potential H5N1 pandemic, the authors in [11] simulate a possible outbreak in Italy by analyzing the structure and the effect of including more specific details like heterogeneities and stochasticity. More specifically, the heterogeneity was with respect to age and spatial distribution, and the stochasticity regulated the epidemic dynamics when the number of the infected population is low. This resulted in the SEIR model:
(6)
d
d
t
S
i
p
=
-
S
i
p
∑
j
,
q
β
i
j
p
q
I
j
q
N
j
q
,
d
d
t
I
i
p
=
ν
E
i
p
-
γ
I
i
p
,
d
d
t
E
i
p
=
S
i
p
∑
j
,
q
β
i
j
p
q
I
j
q
N
j
q
-
ν
E
i
p
,
d
d
t
R
i
p
=
I
i
p
,
with i
,
j
=
1
,
…
,
12 and p
,
q
=
1
,
…
,
20, where 1
/
ν and 1
/
γ, respectively, represent the mean length of the latency and infectivity periods.Other fairly recent papers related to the modeling of pandemics are worth mentioning here as well. For example, in [12] genetic algorithms were used in order to determine optimal vaccination strategies for pandemic influenza. Also, in [13] in order to analyze recurrent outbreaks of the avian H5N1 influenza virus in Asia, which constitute a global pandemic threat, the authors analyzed the great influenza pandemic of 1918 in Geneva, Switzerland, which is also referred to as the “Spanish flu.” Lastly, [14] shows how previous constraints regarding continuous space contact models describing a class of multitype epidemics can be removed in order to construct a general proof of the pandemic bonds for the multitype epidemic on the lattice Z
N. For further interesting results related to recent studies of pandemics consider [15–17].Our current model is quite different than the aforementioned ones, since it comprises two nonlocal integrodifferential equations, also often referred to as functional differential equations, and an ordinary differential equation. Moreover, in order to establish the existence of the weak solution to the weak formulation of (4), an interesting Hilbert space is constructed, as established in Section 3. In addition, several facets of the analysis needed in the study of Laplace's equation on the unit sphere are used throughout. We also emphasize that, to the very best of our knowledge, there are currently no published results related to nonlocal global pandemic models.This paper is organized as follows. In Section2, we derive a weak formulation of (4). Based on this analysis, an auxiliary problem is constructed in Section 3, from which we show that its solution exists and is unique. In Section 4, we prove the existence of a solution to the weak formulation of model (4), using the analysis of Section 3. We conclude the paper with Section 5, wherein we briefly discuss future considerations.


## 2. A Weak Formulation

We begin by multiplying each equation in (4) by sin
ϕ, resulting in
(7)
u
t
sin
ϕ
=
∂
∂
ϕ
(
sin
ϕ
∂
u
∂
ϕ
)
+
csc
ϕ
∂
2
u
∂
θ
2
-
α
(
∫
Ω
‍
w
(
ϕ
,
θ
,
t
)
sin
ϕ
d
ϕ
d
θ
)
u
sin
ϕ
+
γ
w
sin
ϕ
,
on
⁡
S
×
(
0
,
T
]
,
(8)
v
t
sin
ϕ
=
∂
∂
ϕ
(
sin
ϕ
∂
w
∂
ϕ
)
+
csc
ϕ
∂
2
v
∂
θ
2
+
α
(
∫
Ω
‍
w
(
ϕ
,
θ
,
t
)
sin
ϕ
d
ϕ
d
θ
)
u
sin
ϕ
-
β
v
sin
ϕ
,
on
⁡
S
×
(
0
,
T
]
,
(9)
w
t
=
β
v
(
ϕ
,
θ
,
t
)
-
γ
w
(
ϕ
,
θ
,
t
)
,
on
⁡
S
×
(
0
,
T
]
,
(10)
u
(
ϕ
,
θ
,
0
)
=
f
(
ϕ
,
θ
)
on
⁡
S
,
v
(
ϕ
,
θ
,
0
)
=
0
on
⁡
S
,
h
h
h
h
w
(
ϕ
,
θ
,
0
)
=
0
on
⁡
Ω
⊂
S
.
Moreover, u, v, and w, as well as ∂
u
/
∂
θ, ∂
v
/
∂
θ, and ∂
w
/
∂
θ, are all continuous across the data line, and all of the derivatives with respect to ϕ of u, v, and w exist for every direction θ when ϕ
=
0
,
π.In considering a weak solution to (7)–(9), we need to establish appropriate test functions. Consider g
(
ϕ
,
θ
,
t
)
∈
C
1 in ϕ, θ, and t and let g and g
θ be continuous across the data line. Now multiply (7) by g
(
ϕ
,
θ
,
t
) and let g
(
ϕ
,
θ
,
T
)
=
0; then
(11)
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
g
u
t
sin
ϕ
d
ϕ
d
θ
d
t
=
-
∫
-
π
π
‍
∫
0
π
‍
g
(
ϕ
,
θ
,
0
)
f
(
ϕ
,
θ
)
sin
ϕ
d
ϕ
d
θ
-
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
‍
u
g
t
sin
ϕ
d
ϕ
d
θ
d
t
.
The next term in (7) yields
(12)
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
g
∂
∂
ϕ
(
sin
ϕ
u
ϕ
)
d
ϕ
d
θ
d
t
=
∫
0
T
‍
∫
-
π
π
g
u
ϕ
sin
ϕ
∣
ϕ
=
0
ϕ
=
π
d
θ
d
t
-
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
g
u
ϕ
sin
ϕ
d
ϕ
d
θ
d
t
.
We also have
(13)
∫
0
T
‍
∫
0
π
‍
∫
-
π
π
csc
ϕ
g
u
θ
θ
d
θ
ϕ
d
ϕ
d
t
=
∫
0
T
‍
∫
0
π
csc
ϕ
g
u
θ
∣
θ
=
-
π
θ
=
π
d
ϕ
d
t
-
∫
0
T
‍
∫
0
π
‍
∫
-
π
π
csc
ϕ
g
θ
u
θ
d
θ
d
ϕ
d
t
.
Rearranging terms we obtain
(14)
-
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
u
g
t
sin
ϕ
d
ϕ
d
θ
d
t
-
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
g
ϕ
u
ϕ
sin
ϕ
d
ϕ
d
θ
d
t
+
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
csc
ϕ
g
θ
u
θ
d
ϕ
d
θ
d
t
,
=
∫
-
π
π
‍
∫
0
π
g
(
ϕ
,
θ
,
0
)
f
(
ϕ
,
θ
)
sin
ϕ
d
ϕ
d
θ
-
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
α
(
∫
Ω
‍
w
sin
ϕ
d
ϕ
d
θ
)
g
ϕ
u
sin
ϕ
d
ϕ
d
θ
d
t
+
γ
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
g
ϕ
w
sin
ϕ
d
ϕ
d
θ
,
for all g
∈
C
1
(
Q
T
) and Q
T
=
{
(
ϕ
,
θ
,
t
)
∣
0
≤
ϕ
≤
π
,
-
π
≤
θ
≤
π
,
0
<
t
≤
T
}.Now letD
t
:
=
{
(
ϕ
,
θ
,
t
)
∣
0
≤
ϕ
≤
π
,
-
π
≤
θ
≤
π
}. For ease of notation we assign the following:
(15)
∫
Q
T
u
g
t
sin
ϕ
d
Q
T
:
=
∫
0
T
‍
∫
-
π
π
‍
∫
0
π
u
g
t
sin
ϕ
d
ϕ
d
θ
d
t
,
∫
D
t
u
g
t
sin
ϕ
d
D
t
:
=
∫
-
π
π
‍
∫
0
π
u
g
t
sin
ϕ
d
ϕ
d
θ
,
∫
Ω
w
d
Ω
:
=
∫
Ω
‍
∫
w
(
ϕ
,
θ
,
t
)
sin
ϕ
d
ϕ
d
θ
,
where
Ω
⊆
S
.
We can then write the weak form of (7) as follows:
(16)
-
∫
Q
T
u
g
t
sin
ϕ
d
Q
T
-
∫
Q
T
u
ϕ
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
u
θ
g
θ
csc
ϕ
d
Q
T
=
∫
D
0
g
f
sin
ϕ
d
D
0
-
∫
Q
T
α
(
∫
Ω
‍
w
d
Ω
)
u
g
sin
ϕ
d
Q
T
+
γ
∫
Q
T
g
w
sin
ϕ
d
Q
T
,
for all g
∈
C
1
(
Q
T
) and g
(
ϕ
,
θ
,
T
)
=
0, or equivalently for all g
∈
W
1,1, which is the closure of C
1
(
Q
T
) with g
(
ϕ
,
θ
,
T
)
=
0 in the norm
(17)
∥
g
∥
W
1,1
2
:
=
∫
Q
T
g
2
sin
ϕ
d
Q
T
+
∫
Q
T
‍
g
θ
2
sin
ϕ
d
Q
T
+
∫
Q
T
g
θ
2
csc
ϕ
d
Q
T
+
∫
Q
T
g
t
2
sin
ϕ
d
Q
T
.
For the weak form of (8) we have
(18)
-
∫
Q
T
‍
v
g
t
sin
ϕ
d
Q
T
+
∫
Q
T
‍
v
ϕ
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
‍
v
θ
g
θ
csc
ϕ
d
Q
T
=
∫
Q
T
‍
α
(
∫
Ω
w
d
Ω
)
u
g
sin
ϕ
d
Q
T
-
β
∫
Q
T
‍
v
g
sin
ϕ
d
Q
T
,
∀
g
∈
W
1,1
.
On the other hand, the weak form of (9) is the pointwise integral representation for almost all (
ϕ
,
θ
,
t
)
∈
Q
T, which is
(19)
w
=
β
∫
0
t
v
e
-
γ
(
t
-
τ
)
d
τ
.
The integral representation (19) of w is a nonlinear functional equation for w as ∫
Ω
‍
w
d
Ω occurs in the function α in both (16) and (18) for u and v.


## 3. The Auxiliary Problem

We replace∫
Ω
‍
w
d
Ω in the function α in (16) and (18) with the continuous function q
(
t
) to yield the auxiliary problem of finding u
(
ϕ
,
θ
,
t
;
q
), v
(
ϕ
,
θ
,
t
;
q
), and w
(
ϕ
,
θ
,
t
;
q
) satisfying
(20)
-
∫
Q
T
‍
u
g
t
sin
ϕ
d
Q
T
+
∫
Q
T
‍
u
ϕ
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
‍
u
θ
g
θ
csc
ϕ
d
Q
T
=
∫
D
0
g
f
sin
ϕ
d
D
0
-
∫
Q
T
α
(
q
(
t
)
)
u
g
sin
ϕ
d
Q
T
+
γ
∫
Q
T
‍
g
w
sin
ϕ
d
Q
T
,
∀
g
∈
W
1,1
,
-
∫
Q
T
v
g
t
sin
ϕ
d
Q
T
+
∫
Q
T
v
ϕ
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
v
θ
g
θ
csc
ϕ
d
Q
T
=
∫
Q
T
α
(
q
(
t
)
)
u
g
sin
ϕ
d
Q
T
-
β
∫
Q
T
v
g
sin
ϕ
d
Q
T
,
∀
g
∈
W
1,1
,
w
(
ϕ
,
θ
,
t
;
q
)
=
β
∫
0
t
v
(
ϕ
,
θ
,
τ
;
q
)
e
-
γ
(
t
-
τ
)
d
τ
for almost all (
ϕ
,
θ
,
t
)
∈
Q
T. It should also be noted that we have already assumed that g is periodic in θ. Also, we can extend g as an odd function in ϕ with -
π
≤
ϕ
≤
0 and regard it as periodic in ϕ. Therefore, we can utilize the functions sin
n
ϕ
sin
m
θ and sin
n
ϕ
cos
⁡
m
θ in order to generate an indexed set of functions h
k
(
ϕ
,
θ
) by first counting the functions sin
n
ϕ
cos
⁡
m
θ according to an orderly method of counting the ordered pairs (
n
,
m
) for n
=
1,2
,
… and m
=
0,1
,
2
,
… in the closure of the x-y plane via the diagonal pattern
(21)
(
1,0
)
⟶
(
1,1
)
⟶
(
2,0
)
⟶
(
3,0
)
⟶
(
2,1
)
⟶
(
1,2
)
⟶
(
1,3
)
⟶
(
2,2
)
⟶
(
3,1
)
⟶
(
4,0
)
⟶
(
5,0
)
⟶
(
4,1
)
⟶
(
3,2
)
⟶
(
2,3
)
⟶
(
1,4
)
⟶
(
1,5
)
⟶
(
2,4
)
⋯
.
Next, take the number index assigned to each pair (
n
,
m
) and double it, yielding the h
k-terms for even integers. Then, for the odd integers, as we did above, count the functions sin
n
ϕ
cos
⁡
m
θ according to an orderly method of counting the ordered pairs (
n
,
m
) for n
=
1,2
,
… and m
=
1,2
,
… in the closure of the x-y plane via the diagonal pattern
(22)
(
1,1
)
⟶
(
1,2
)
⟶
(
2,1
)
⟶
(
3,1
)
⟶
(
2,2
)
⟶
(
1,3
)
⟶
(
1,4
)
⟶
(
2,3
)
⟶
(
3,2
)
⟶
(
4,1
)
⟶
(
5,1
)
⟶
(
4,2
)
⟶
(
3,3
)
⟶
(
2,4
)
⟶
(
1,5
)
⟶
(
1,6
)
⟶
(
2,5
)
⋯
.
So, we have a list of h
k
(
ϕ
,
θ
) for k
=
1,2
,
…. The first six of such are as follows:
(23)
h
1
=
sin
(
ϕ
)
sin
(
θ
)
,
h
2
=
sin
(
ϕ
)
,
h
3
=
sin
(
ϕ
)
sin
(
2
θ
)
,
h
4
=
sin
(
ϕ
)
cos
⁡
(
θ
)
,
h
5
=
sin
(
2
ϕ
)
sin
(
θ
)
,
h
6
=
sin
(
2
ϕ
)
,
…
.Since the integrals involving theh
k
2-term, given by
(24)
∫
-
π
π
‍
∫
0
π
‍
sin
2
n
ϕ
cos
⁡
2
m
θ
sin
ϕ
d
ϕ
d
θ
,
∫
-
π
π
‍
∫
0
π
‍
sin
2
n
ϕ
sin
2
m
θ
sin
ϕ
d
ϕ
d
θ
,
converge, and the integrals involving the ∂
h
k
2
/
∂
ϕ-term, which are
(25)
∫
-
π
π
‍
∫
0
π
‍
n
2
cos
⁡
2
n
ϕ
cos
⁡
2
m
θ
sin
ϕ
d
ϕ
d
θ
,
∫
-
π
π
‍
∫
0
π
‍
n
2
cos
⁡
2
n
ϕ
sin
2
m
θ
sin
ϕ
d
ϕ
d
θ
,
also converge, we then consider the integrals involving the ∂
h
k
2
/
∂
θ-term below:
(26)
∫
-
π
π
‍
∫
0
π
‍
m
2
csc
ϕ
sin
2
n
ϕ
sin
2
m
θ
d
ϕ
d
θ
,
∫
-
π
π
‍
∫
0
π
‍
m
2
csc
ϕ
sin
2
n
ϕ
cos
⁡
2
m
θ
d
ϕ
d
θ
and remark that
(27)
∫
0
π
‍
m
2
csc
ϕ
sin
2
n
ϕ
d
ϕ
<
n
2
m
2
π
4
8
.
Hence,
(28)
∫
-
π
π
‍
∫
0
π
h
k
2
sin
ϕ
d
ϕ
d
θ
+
∫
-
π
π
‍
∫
0
π
(
∂
h
k
∂
ϕ
)
2
sin
ϕ
d
ϕ
d
θ
+
∫
-
π
π
‍
∫
0
π
(
∂
h
k
∂
θ
)
2
csc
ϕ
d
ϕ
d
θ
<
∞
,
for k
=
1,2
,
…. Thus, via the Gram-Schmidt process, we can generate a countable orthonormal basis for the Hilbert space L
2
(
D
0
) with the inner product
(29)
〈
f
(
ϕ
,
θ
)
,
g
(
ϕ
,
θ
)
〉
D
0
=
∫
D
0
f
g
sin
ϕ
d
D
0
+
∫
D
0
∂
f
∂
ϕ
∂
g
∂
ϕ
sin
ϕ
d
D
0
+
∫
D
0
∂
f
∂
θ
∂
g
∂
θ
csc
ϕ
d
D
0
,
where D
0
=
{
(
ϕ
,
θ
)
∣
0
≤
ϕ
≤
π
,
-
π
≤
θ
≤
π
}. Henceforth, we denote the orthonormal basis formed from h
k
(
ϕ
,
θ
) via the Gram-Schmidt process by {
h
k
(
ϕ
,
θ
)
}.Now, define(30)
u
N
(
ϕ
,
θ
,
t
)
=
∑
k
=
1
N
‍
a
k
N
h
k
(
ϕ
,
θ
)
,
v
N
(
ϕ
,
θ
,
t
)
=
∑
k
=
1
N
‍
b
k
N
h
k
(
ϕ
,
θ
)
,
w
N
(
ϕ
,
θ
,
t
)
=
∑
k
=
1
N
‍
c
k
N
h
k
(
ϕ
,
θ
)
.
From there, we construct a system of ODEs as follows:
(31)
∫
D
0
u
n
h
i
sin
ϕ
d
D
0
+
∫
D
0
∂
u
N
∂
ϕ
∂
h
i
∂
ϕ
sin
ϕ
d
D
0
+
∫
D
0
∂
u
N
∂
θ
∂
h
i
∂
θ
csc
ϕ
d
D
0
=
-
∫
D
0
‍
α
(
q
(
t
)
)
u
N
h
i
sin
ϕ
d
D
0
+
γ
∫
D
0
w
N
h
i
sin
ϕ
d
D
0
,
∫
D
0
v
n
h
i
sin
ϕ
d
D
0
+
∫
D
0
∂
v
N
∂
ϕ
∂
h
i
∂
ϕ
sin
ϕ
d
D
0
+
∫
D
0
∂
v
N
∂
θ
∂
h
i
∂
θ
csc
ϕ
d
D
0
=
∫
D
0
α
(
q
(
t
)
)
u
N
h
i
sin
ϕ
d
D
0
-
β
∫
D
0
v
N
h
i
sin
ϕ
d
D
0
,
∫
D
0
w
n
h
i
sin
ϕ
d
D
0
=
β
∫
D
0
v
N
h
i
sin
ϕ
d
D
0
hhhh
-
γ
∫
D
0
w
N
h
i
sin
ϕ
d
D
0
,
a
i
N
(
0
)
=
〈
f
,
h
i
〉
D
0
,
b
i
N
(
0
)
=
c
i
N
(
0
)
=
0
,
i
=
1,2
,
…
,
N
.
The equations above become
(32)
d
a
i
N
(
t
)
d
t
+
∑
l
=
1
N
‍
A
i
,
l
N
a
l
N
(
t
)
=
-
α
(
q
(
t
)
)
a
i
N
(
t
)
+
γ
c
i
N
(
t
)
,
d
b
i
N
(
t
)
d
t
+
∑
l
=
1
N
‍
B
i
,
l
N
b
l
N
(
t
)
=
α
(
q
(
t
)
)
a
i
N
(
t
)
-
β
b
i
N
(
t
)
,
d
c
i
N
(
t
)
d
t
=
β
b
i
N
(
t
)
-
γ
c
i
N
(
t
)
,
a
i
N
(
0
)
=
〈
f
,
h
i
〉
D
0
,
b
i
N
(
0
)
=
c
i
N
(
0
)
=
0
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
i
=
1,2
,
…
,
N
.
As the above equations are linear, for each positive integer N there exists a solution on [
0
,
T
] since α
(
q
(
t
)
) is bounded for q
(
t
)
≥
0. Hence, for each N, u
N, v
N, and w
N are determined uniquely. We note that
(33)
c
i
N
(
t
)
=
β
∫
0
t
b
i
N
(
τ
)
e
γ
(
t
-
τ
)
d
τ
.
Upon multiplying each of the equations in the system (32) by the coefficient corresponding to the t-derivative in the equation and summing each equation over i
=
1
,
…
,
N, we achieve
(34)
〈
d
u
N
d
t
,
u
N
〉
D
0
+
∫
D
0
‍
(
∂
u
N
∂
ϕ
)
2
sin
ϕ
d
D
0
+
∫
D
0
‍
(
∂
u
N
∂
θ
)
2
csc
ϕ
d
D
0
=
-
∫
D
0
‍
α
(
q
(
t
)
)
u
N
2
sin
ϕ
d
D
0
+
γ
∫
D
0
‍
w
N
u
N
sin
ϕ
d
D
0
,
〈
d
v
N
d
t
,
v
N
〉
D
0
+
∫
D
0
‍
(
∂
v
N
∂
ϕ
)
2
sin
ϕ
d
D
0
+
∫
D
0
‍
(
∂
v
N
∂
θ
)
2
csc
ϕ
d
D
0
=
∫
D
0
‍
α
(
q
(
t
)
)
u
N
v
N
sin
ϕ
d
D
0
+
β
∫
D
0
‍
v
N
2
sin
ϕ
d
D
0
,
〈
d
w
N
d
t
,
w
N
〉
D
0
=
β
∫
D
0
‍
v
N
w
N
sin
ϕ
d
D
0
hhhhhhhhhhhhhh
-
γ
∫
D
0
‍
w
N
2
sin
ϕ
d
D
0
and integrating each of the above equations from 0 to t we obtain
(35)
1
2
∥
u
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
+
∫
Q
t
(
∂
u
N
∂
ϕ
)
2
sin
ϕ
d
Q
t
+
∫
Q
t
‍
(
∂
u
N
∂
θ
)
2
csc
ϕ
d
Q
t
=
1
2
∥
u
N
(
ϕ
,
θ
,
0
)
∥
D
0
2
-
∫
Q
t
α
(
q
(
t
)
)
u
N
2
sin
ϕ
d
Q
t
+
γ
∫
Q
t
‍
w
N
u
N
sin
ϕ
d
Q
t
,
1
2
∥
v
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
+
∫
Q
t
(
∂
v
N
∂
ϕ
)
2
sin
ϕ
d
Q
t
+
∫
Q
t
‍
(
∂
v
N
∂
θ
)
2
csc
ϕ
d
Q
t
=
∫
Q
t
α
(
q
(
t
)
)
u
N
v
N
sin
ϕ
d
Q
t
-
β
∫
Q
t
v
N
2
sin
ϕ
d
Q
t
,
1
2
∥
w
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
=
β
∫
Q
t
v
N
w
N
sin
ϕ
d
Q
t
hhhh
-
γ
∫
Q
t
w
N
2
sin
ϕ
d
Q
t
.
Multiplying the above equations by 2, utilizing the inequality 2
a
b
≤
a
2
+
b
2 on the right-hand sides of each equation, adding the inequalities, and dropping the Q
t integrals on the left-hand side of the resulting inequalities, we have
(36)
F
(
t
)
≤
∥
f
∥
D
0
2
+
C
1
∫
0
t
F
(
τ
)
d
τ
,
where
(37)
F
(
t
)
:
=
∥
u
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
+
∥
v
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
+
∥
w
N
(
ϕ
,
θ
,
t
)
∥
D
0
2
,
C
1
=
3
∥
α
∥
∞
+
3
β
+
2
γ
+
1
and ∥
u
N
(
ϕ
,
θ
,
0
)
∥
D
0
2
≤
∥
f
∥
D
0
2.Then, by Gronwall’s lemma we see that(38)
y
′
≤
∥
f
∥
D
0
2
+
C
1
y
d
d
t
(
y
e
-
C
1
t
)
≤
∥
f
∥
D
0
2
e
-
C
1
t
y
(
t
)
≤
∥
f
∥
D
0
2
C
1
[
e
C
1
t
-
1
]
.
Therefore,
(39)
F
(
t
)
≤
∥
f
∥
D
0
2
e
C
1
t
.
In addition, we have
(40)
F
(
t
)
+
2
∫
Q
t
(
∂
u
N
∂
ϕ
)
2
sin
ϕ
d
Q
t
+
2
∫
Q
t
(
∂
u
N
∂
θ
)
2
csc
ϕ
d
Q
t
+
2
∫
Q
t
(
∂
v
N
∂
ϕ
)
2
sin
ϕ
d
Q
t
+
2
∫
Q
t
(
∂
v
N
∂
θ
)
2
csc
ϕ
d
Q
t
≤
∥
f
∥
D
0
2
+
C
1
∫
0
t
F
(
τ
)
d
τ
and we observe that the entire left-hand side of the above inequality is also bounded by ∥
f
∥
D
0
2
e
C
1
t for each t
∈
[
0
,
T
] and so is each individual term on the left-hand side.From there, we see that(41)
∫
Q
T
u
N
2
sin
ϕ
d
Q
T
,
∫
Q
T
v
N
2
sin
ϕ
d
Q
T
,
∫
Q
T
w
N
2
sin
ϕ
d
Q
T
are all bounded uniformly with respect to N. Likewise, u
N and v
N are also uniformly bounded in the Hilbert space V with inner product
(42)
〈
h
,
g
〉
V
=
∫
Q
T
h
g
sin
ϕ
d
Q
T
+
∫
Q
T
h
ϕ
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
h
θ
g
θ
csc
ϕ
d
Q
T
,
where h
ϕ, g
ϕ, h
θ, and g
θ are the first partial derivatives of h and g, with respect to the subscripted variables.Next, we take (33) and multiply it by h
i
(
ϕ
,
θ
) and sum with i from 1 to N and obtain
(43)
w
N
(
ϕ
,
θ
,
t
)
=
β
∫
0
t
v
N
(
ϕ
,
θ
,
τ
)
e
-
γ
(
t
-
τ
)
d
τ
,
hhhhhhhhhh
∀
(
ϕ
,
θ
,
t
)
∈
Q
T
.
Thus, by a Cantor diagonalization process, there exist sequences {
u
N
}, {
v
N
}, and {
w
N
}, which converge weakly in L
2
(
Q
T
) with {
u
N
} and {
v
N
} also converging weakly in V. Let u, v, and w denote the weak limits and let
(44)
g
m
(
ϕ
,
θ
,
t
)
=
∑
i
=
1
m
e
i
(
t
)
h
i
(
ϕ
,
θ
)
denote an approximation in W
2
1,1
(
Q
T
) to g
(
ϕ
,
θ
,
t
)
∈
C
1
(
Q
T
) with e
i
(
T
)
=
0 for all i
=
1,2
,
…; that is, g
m
→
g strongly in W
2
1,1
(
Q
T
) and g
m
(
ϕ
,
θ
,
T
)
=
g
(
ϕ
,
θ
,
T
)
=
0. Taking the system of ODEs and multiplying each by e
i
(
t
) and summing each block of equations with respect to i transform the system as follows:
(45)
〈
∂
u
N
∂
t
,
g
m
〉
D
0
+
∫
D
0
∂
u
N
∂
ϕ
∂
u
N
∂
ϕ
g
m
sin
ϕ
d
D
0
+
∫
D
0
∂
u
N
∂
θ
∂
u
N
∂
θ
g
m
csc
ϕ
d
D
0
=
∫
D
0
‍
α
(
q
(
t
)
)
u
N
g
m
sin
ϕ
d
D
0
+
γ
∫
D
0
‍
w
N
g
m
sin
ϕ
d
D
0
,
〈
∂
v
N
∂
t
,
g
m
〉
D
0
+
∫
D
0
‍
∂
v
N
∂
ϕ
∂
v
N
∂
ϕ
g
m
sin
ϕ
d
D
0
+
∫
D
0
‍
∂
v
N
∂
θ
∂
v
N
∂
θ
g
m
csc
ϕ
d
D
0
=
∫
D
0
‍
α
(
q
(
t
)
)
u
N
g
m
sin
ϕ
d
D
0
-
β
∫
D
0
‍
v
N
g
m
sin
ϕ
d
D
0
,
∫
Q
T
‍
(
w
N
-
β
∫
0
t
‍
v
N
e
-
γ
(
t
-
τ
)
d
τ
)
g
m
d
Q
T
=
0
.
Integrating with respect to t from 0 to T and performing integration by parts on the first term in each equation, we observe that u
N, v
N, and w
N satisfy the weak formulation for g
m. Then holding g
m fixed and letting N
→
∞ we see that the weak limits u, v, and w satisfy the original system (20) for each m. As g
m
→
g
∈
C
1
(
Q
t
) in the norm of W
2
1,1
(
Q
t
), then u, v, and w satisfy the weak formulation for all g
∈
C
1
(
Q
t
)
∩
W
2
1,1 and g
(
ϕ
,
θ
,
T
)
=
0. Thus, as g is dense in W
2
1,1 the weak limits u, v, and w together form a weak solution to the auxiliary problem (20).Upon considering the energy inequality for parabolic equations (cf. [18]), it can be shown that any weak solution for u, v, and w satisfies an energy inequality similar to the one for {
u
N
,
v
N
,
w
N
}. Hence, u, v, and w satisfy the same bounds as u
N, v
N, and w
N. Consequently, by the linearity of the system f
(
ϕ
,
θ
)
≡
0, we have u
=
v
=
w
=
0 and thus, the weak solution is unique.


## 4. Existence of a Solution to the Auxiliary Problem

As the auxiliary problem possesses a unique weak solution for each continuousq
(
t
), we can now define a mapping F
:
C
(
[
0
,
T
]
)
→
C
(
[
0
,
T
]
) via
(46)
(
F
q
)
(
t
)
=
∫
Ω
w
(
ϕ
,
θ
,
t
;
q
)
d
Ω
.
First, we show that (
F
q
)
(
t
) is bounded and continuous on [
0
,
T
]. From integrating the pointwise representation for w
(
ϕ
,
θ
,
t
;
q
), we see that
(47)
∫
Ω
w
d
Ω
=
β
∫
Ω
∫
0
t
v
e
-
γ
(
t
-
τ
)
d
τ
d
Ω
and, taking into account the fact that 0
<
e
-
γ
(
t
-
τ
)
≤
1, we have
(48)
∫
Ω
‍
w
d
Ω
≤
β
∫
0
t
‍
∫
Ω
‍
|
v
|
d
Ω
d
τ
=
β
|
∫
0
t
‍
∫
Ω
‍
|
v
|
sin
ϕ
d
θ
d
ϕ
d
τ
|
≤
β
(
∫
0
t
‍
∫
Ω
‍
sin
ϕ
d
θ
d
ϕ
d
τ
)
1
/
2
×
(
∫
0
t
‍
∫
Ω
‍
v
2
sin
ϕ
d
θ
d
ϕ
d
τ
)
1
/
2
≤
β
T
1
/
2
[
mes
(
Ω
)
]
1
/
2
(
∫
Q
T
‍
v
2
sin
ϕ
d
Q
T
)
1
/
2
≤
C
2
,
via the Schwartz inequality, where C
2 is a positive constant that depends on C
1 of Section 3, T, and mes
(
Ω
).For the continuity of(
F
q
)
(
t
), we consider 0
≤
t
1
<
t
2
≤
T, and thus
(49)
|
∫
Ω
‍
w
(
ϕ
,
θ
,
t
1
;
q
)
d
Ω
-
∫
Ω
‍
w
(
ϕ
,
θ
,
t
2
;
q
)
d
Ω
|
≤
β
∫
Ω
‍
|
∫
0
t
1
‍
v
(
ϕ
,
θ
,
τ
;
q
)
e
-
γ
(
t
1
-
τ
)
d
τ
h
h
h
h
h
h
h
h
h
-
β
∫
0
t
2
‍
v
(
ϕ
,
θ
,
τ
;
q
)
e
-
γ
(
t
2
-
τ
)
d
τ
|
d
Ω
≤
β
∫
Ω
‍
|
∫
t
1
t
2
‍
v
e
-
γ
(
t
2
-
τ
)
d
τ
|
d
Ω
+
β
∫
Ω
‍
|
∫
0
t
1
‍
v
[
e
-
γ
(
t
1
-
τ
)
-
e
-
γ
(
t
2
-
τ
)
]
d
τ
|
d
Ω
≤
β
∫
Ω
‍
∫
t
1
t
2
‍
|
v
|
d
τ
d
Ω
+
β
∫
Ω
‍
∫
0
t
1
‍
|
v
|
e
-
γ
(
t
1
-
τ
)
[
1
-
e
γ
(
t
1
-
τ
)
-
γ
(
t
2
-
τ
)
]
d
τ
d
Ω
≤
β
∫
Ω
‍
∫
t
1
t
2
‍
|
v
|
d
τ
d
Ω
+
β
∫
Ω
‍
∫
0
t
1
‍
|
v
|
[
1
-
e
-
γ
(
t
2
-
t
1
)
]
d
τ
d
Ω
≤
β
∫
Ω
‍
∫
t
1
t
2
‍
|
v
|
d
τ
d
Ω
+
β
∫
Ω
‍
∫
0
t
1
‍
|
v
|
γ
(
t
2
-
t
1
)
d
τ
d
Ω
≤
C
3
(
t
2
-
t
1
)
1
/
2
+
C
4
|
t
2
-
t
1
|
,
where we used the fact that 0
<
e
-
γ
(
t
-
τ
)
≤
1, estimated 1
-
e
-
γ
(
t
2
-
t
1
) via the mean value theorem, and estimated the last two integrals exactly as the integral in (48). We note that C
3 and C
4 have the same dependence as C
2, which is independent of q since C
1 is independent of q.It therefore follows from (48) and (49) that F maps C
(
[
0
,
T
]
) into a convex subset of itself. All that remains for an application of Schauder's fixed point theorem is to show that F is continuous in the norm ∥
·
∥
∞ of the separable Banach space B
(
C
(
[
0
,
T
]
)
,
∥
·
∥
∞
).To show that(
F
q
)
(
t
) is continuous in ∥
·
∥
∞, we consider
(50)
(
F
q
1
)
(
t
)
-
(
F
q
2
)
(
t
)
=
β
∫
Ω
‍
∫
0
t
[
v
(
ϕ
,
θ
,
τ
;
q
1
)
-
v
(
ϕ
,
θ
,
τ
;
q
2
)
]
e
-
γ
(
t
-
τ
)
d
τ
d
Ω
,
where v
(
ϕ
,
θ
,
τ
;
q
i
) is the weak solution of (20) for continuous q
i
(
t
), 0
≤
t
≤
T, for i
=
1,2. Forming the difference of (20) for q
i
(
t
), for i
=
1,2, we have
(51)
-
∫
Q
T
(
u
1
-
u
2
)
g
t
sin
ϕ
d
Q
T
+
∫
Q
T
(
u
1
ϕ
-
u
2
ϕ
)
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
(
u
1
θ
-
u
2
θ
)
g
θ
csc
ϕ
d
Q
T
=
-
∫
Q
T
[
α
(
q
1
)
u
1
-
α
(
q
2
)
u
2
]
g
sin
ϕ
d
Q
T
+
γ
∫
Q
T
(
w
1
-
w
2
)
g
sin
ϕ
d
Q
T
,
∀
g
∈
W
1,1
,
(52)
-
∫
Q
T
(
v
1
-
v
2
)
g
t
sin
ϕ
d
Q
T
+
∫
Q
T
(
v
1
ϕ
-
v
2
ϕ
)
g
ϕ
sin
ϕ
d
Q
T
+
∫
Q
T
(
v
1
θ
-
v
2
θ
)
g
θ
csc
ϕ
d
Q
T
=
∫
Q
T
[
α
(
q
1
)
u
1
-
α
(
q
2
)
u
2
]
g
sin
ϕ
d
Q
T
-
β
∫
Q
T
(
v
1
-
v
2
)
g
sin
ϕ
d
Q
T
,
∀
g
∈
W
1,1
,
(53)
w
1
(
ϕ
,
θ
,
t
;
q
1
)
-
w
2
(
ϕ
,
θ
,
t
;
q
2
)
=
β
∫
0
t
‍
[
v
1
(
ϕ
,
θ
,
t
;
q
1
)
-
v
2
(
ϕ
,
θ
,
t
;
q
2
)
]
e
-
γ
(
t
-
τ
)
d
τ
,
where the u
i, v
i, and w
i depend upon the q
i, for i
=
1,2, as indicated by the w
i and v
i expressed in (52). Also, in (51) we can write
(54)
∫
Q
T
[
α
(
q
1
)
u
1
-
α
(
q
2
)
u
2
]
g
sin
ϕ
d
Q
T
=
∫
Q
T
[
α
(
q
1
)
-
α
(
q
2
)
]
u
1
g
sin
ϕ
d
Q
T
+
∫
Q
T
α
(
q
2
)
[
u
1
-
u
2
]
g
sin
ϕ
d
Q
T
.
Hence, if the u, v, w, and α
(
q
) as in (20) were, respectively, replaced by u
1
-
u
2, v
1
-
v
2, w
1
-
w
2, and α
(
q
2
) above, and if the term ∫
D
0
‍
g
f
sin
ϕ
d
D
0 was deleted, then the right-hand side of (52) would only contain the term
(55)
∫
Q
T
‍
[
α
(
q
1
)
-
α
(
q
2
)
]
u
1
g
sin
ϕ
d
Q
T
to be estimated, instead of the initial data term ∫
D
0
‍
g
f
sin
ϕ
d
D
0. Now the term (55) can be estimated as follows:
(56)
|
∫
Q
T
‍
[
α
(
q
1
)
-
α
(
q
2
)
]
u
1
g
sin
ϕ
d
Q
T
|
≤
[
∫
Q
T
‍
[
α
(
q
1
(
τ
)
)
-
α
(
q
2
(
τ
)
)
]
2
u
1
2
sin
ϕ
d
Q
T
]
1
/
2
×
[
∫
Q
T
‍
g
2
sin
ϕ
d
Q
T
]
1
/
2
≤
1
2
∫
Q
T
‍
[
α
(
q
1
(
τ
)
)
-
α
(
q
2
(
τ
)
)
]
2
u
1
2
sin
ϕ
d
Q
T
+
1
2
∫
Q
T
‍
g
2
sin
ϕ
d
Q
T
≤
1
2
(
∥
α
∥
∞
+
∥
α
∥
∞
)
∥
d
α
d
q
∥
∞
∥
q
1
-
q
2
∥
∞
×
∫
Q
T
‍
u
1
2
sin
ϕ
d
Q
T
+
1
2
∫
Q
T
‍
g
2
sin
ϕ
d
Q
T
≤
C
5
∥
q
1
-
q
2
∥
∞
+
1
2
∫
Q
T
‍
g
2
sin
ϕ
d
Q
T
,
where C
5 depends on the bound for ∥
u
1
∥
Q
T which is given by (39) for the weak solution of (20) and is independent of q
1. The second term on the right-hand side of (54) and (56) involving the g becomes u
1
-
u
2 in the energy inequality associated with (52). We now have a replay of (37), (38), and (39), with ∥
f
∥
D
0
2 replaced with C
5
∥
q
1
-
q
2
∥
∞. From the estimate of ∥
v
1
-
v
2
∥
Q
T
2 resulting from the Galerkin estimate for ∥
u
1
-
u
2
∥
Q
T
2
+
∥
v
1
-
v
2
∥
Q
T
2
+
∥
w
1
-
w
2
∥
Q
T
2 in terms of C
6
∥
q
1
-
q
2
∥
∞, we obtain the estimation of (50) below:
(57)
∥
F
(
q
1
)
-
F
(
q
2
)
∥
∞
≤
C
7
∥
q
1
-
q
2
∥
∞
.
This leads to the following statement.Theorem 1.
From (57), one infers the continuity of the mapping F and the validity of the application of the Schauder fixed point theorem, yielding the existence of a weak solution of the weak formulation of (7), (8), and (9), given by (16), (18), and (19).Proof.
See the preceding analysis.


## 5. Future Considerations

We leave open for consideration the establishment of a unique solution to the weak formulation developed in this paper. This will entail a numerous amount of estimates, which are too tedious to be included in the present work. We refer the interested reader to Chapter 3, Section 2, of [18], as we anticipate that analogues of the energy inequalities derived therein can be generalized to ones that apply to our current problem. We also anticipate that these inequalities will ultimately produce a contraction mapping via a bootstrapping procedure. The constants in these inequalities will have to be estimated closely enough to show the dependence on (small) t, in order to permit the construction of the aforementioned contraction mapping.



---

*Source: 187685-2014-07-24.xml*