• 検索結果がありません。

On the Well Posedness and Refined Estimates for the Global Attractor of the TYC Model

N/A
N/A
Protected

Academic year: 2022

シェア "On the Well Posedness and Refined Estimates for the Global Attractor of the TYC Model"

Copied!
29
0
0

読み込み中.... (全文を見る)

全文

(1)

Volume 2010, Article ID 405816,29pages doi:10.1155/2010/405816

Research Article

On the Well Posedness and Refined Estimates for the Global Attractor of the TYC Model

Rana D. Parshad

1

and Juan B. Gutierrez

2

1Department of Mathematics and Computer Science, Clarkson University, Potsdam, NY 13676, USA

2Mathematical Biosciences Institute, Ohio State University, Columbus, OH 43210, USA

Correspondence should be addressed to Rana D. Parshad,rparshad@clarkson.edu Received 14 July 2010; Accepted 2 November 2010

Academic Editor: Sandro Salsa

Copyrightq2010 R. D. Parshad and J. B. Gutierrez. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The Trojan Y Chromosome strategyTYC is a theoretical method for eradication of invasive species. It requires constant introduction of artificial individuals into a target population, causing a shift in the sex ratio that ultimately leads to local extinction. In this work we demonstrate the existence of a unique weak solution to the infinite dimensional TYC system. Furthermore, we obtain improved estimates on the upper bounds for the Hausdorffand fractal dimensions of the global attractor of the TYC system, via the use of weighted Sobolev spaces. These results confirm that the TYC eradication strategy is a sound theoretical method of eradication of invasive species in a spatial setting. It also provides a solid ground for experiments in silico and validates the use of the TYC strategy in vivo.

1. Introduction

An exotic species is a species that resides outside its native habitat. When it causes some sort of measurable damage, it is often referred to as invasive. The recent globalization process has expedited the pace at which exotic species are introduced into new environments. Once established, these species can be extremely difficult to manage and almost impossible to eradicate1,2. Studies have indicated that the losses caused by invasive species could be as much as $120 billion/year by 20043. The effect of these invaders is thus devastating4.

Current approaches for controlling exotic fish species are limited to general chemical control methods applied to small water bodies and/or small isolated populations that kill native fish in addition to the target fish5. For example, the piscicide Rotenone has been used to eradicate exotic fish, but at the expense of killing all the endogenous fish, making it necessary to restock native fish from other sources1,2.

A genetic strategy to cause extinction of invasive species was proposed by Gutierrez and Teem6. This strategy is relevant to species amenable to sex reversal and with an XY sex- determination system, in which males are the heterogametic sexcarrying one X chromosome

(2)

and one Y chromosome, XY and females are the homogametic sex carrying two X chromosomes, XX. The strategy relies on the fact that variations in the sex chromosome number can be produced through genetic manipulation, for example, a normal and fertile male bearing two Y chromosomessupermale, YY 7–10. Also hormone treatments can be used to reverse the sex, resulting in a feminized YY supermale5,11,12.

The eradication strategy requires adding a sex-reversed “Trojan” female individual bearing two Y chromosomes, that is, feminized supermales r, at a constant rate μ to a target population of an invasive species, containing normal females and males denoted as f and m, respectively. Matings involving the introduced r generate a disproportionate number of males over time. The higher incidence of males decrease the female to male ratio.

Ultimately, the number offdecline to zero, causing local extinction. This theoretical method of eradication is known as Trojan Y ChromosomeTYCstrategy.

The original model considered by Gutierrez and Teem was an ODE model. Spatial spread is ubiquitous in aquatic settings and was thus considered by Gutierrez et al. 13, resulting in a PDE model. In14, we considered the PDE model and showed the existence of a global attractor for the system, which isH2Ωregular, attracting orbits uniformly in the L2Ωmetric. We showed that this attractor supports a state, in which the female population is driven to zero, thus resulting in local extinction. Recall the TYC model with spatial spread takes the following form14:

∂f

∂t DΔf1

2fmβLδf, f

∂Ω 0, 1.1

∂m

∂t DΔm

1 2fm1

2rmfs

βLδm, m|∂Ω 0, 1.2

∂s

∂t DΔs 1

2rmrs

βLδs, s|∂Ω 0, 1.3

∂r

∂t DΔrμδr, r|∂Ω 0. 1.4

Here,Ω⊂R3is a bounded domain. Also

L 1−

fmrs K

, 1.5

whereKis the carrying capacity of the ecosystem,Dis a diffusivity coefficient,δis a birth coefficienti.e., what proportion of encounters between males and females result in progeny, and δ is a death coefficienti.e., what proportion of the population is dying at any given moment. We assume initial data is positive and inL2Ω. At the outset we would like to point out that the difficulty in analyzing1.1–1.4lies in the nonlinear termsLfm,L1/2fm 1/2rmfsandL1/2rmrs. See15for a PDE dealing with similar nonlinearities, albeit in the setting of a fluid-saturated porous medium. We will also assume positivity of solutions as negativef, m, r, sdo not make sense in the biological context. We also provide a rigoros proof to this end.

(3)

In the current paper we will show that the TYC model,1.1–1.4, possesses a unique weak solutionf, m, r, s. By this we mean that there existf, m, r, ssuch that the following is satisfied in the distributional sense:

d dt

f, v D

∇f,∇v δ

f, v 1

2fmβL, v

, d

dtm, vD∇m,∇vδm, v 1 2fm 1

2rmfs

βL, v

, d

dtr, vD∇s,∇vδs, v 1

2rmrs

βL, v

, d

dts, vD∇s,∇vδs, v μ, v

.

1.6

Here, · is the standard inner product in L2Ω. Furthermore the above hold for all vH01Ω. Our main result is summarized in the following theorem.

Theorem 1.1. Consider the Trojan Y Chromosome model,1.1–1.4. There exists a unique weak solutionf, m, r, sto the system for positive initial data inL2Ω, such that

f, m, r, s

C

0, T;L2Ω

L

0, T;L2Ω

L2

0, T;H01Ω , ∂f

∂t,∂m

∂t ,∂r

∂t,∂s

∂t

L2

0, T;H−1Ω ,

1.7

for allT >0. Furthermore,f, m, r, sare continuous with respect to initial data.

Our strategy to prove the above is as follows: we first derive a priori estimates for thef, m, r, svariables. We then show existence of a solution to1.1. Note, showing existence of a solution to1.1requires a priori estimates onm, r, salso. The key here isLemma 4.1 which enables convergence of the nonlinear term Lfm. Next we show uniqueness of the solution to 1.1. The procedure to show existence and uniqueness of solutions to 1.2–

1.4follow similarly. We then consider the question of sharpening the upper bounds on the Hausdorffand fractal dimension of the global attractor for the system, derived in14. This constitutes our second main result,Theorem 7.2. Lastly, we offer some concluding remarks.

In all estimates made hence, forth,Cis a generic constant that can change in its value from line to line and sometimes within the same line if so required.

2. A Bound in L

Ω

The biology of the system dictates that the solutions are bounded in the supremum norm by the carrying capacity. We now provide a proof via a maximum principle argument.

(4)

Lemma 2.1. Consider the Trojan Y Chromosome model,1.1–1.4. The solutionsf, m, r, sof the system are bounded as follows:

f

K,

|m|K,

|s|K,

|r|K.

2.1

Proof. The proof relies heavily on the form of the nonlinearity in the system. We concentrate on the nonlinear term in1.1,

F

f, m, r, s fm

1−fmrs K

. 2.2

The analysis for the other terms is similar. As is biologically viable, we assumesf,m,r, and sare always positive, thus, we have

f >0, m >0, r >0, s >0. 2.3

Assuming positive initial data,f0 >0,m0 >0,r0 >0, ands0 > 0, the solution at later times remains positive. In order to prove this let us assume the contrary, that isf0>0,m0 >0,r0>0, ands0>0, but sayfcan become negative at a later time. Consider an interior minimum point in the parabolic cyliderΩ×0, T, that is somex, t, such thatf attains a minimum there, and thatfx, t < 0,mx, t < 0,rx, t < 0, andsx, t < 0. Under this setting, from standard calculus, we have

∂f

∂tx, t 0, Δfx, t≥0, 2.4

furthermore,

−δfx, t>0,

βfx, tmx, t

1−fx, t mx, t rx, t sx, t K

>0.

2.5

(5)

Thus from1.1, we have

∂f

∂tx, t 0

Δfx, tδfx, t βfx, tmx, t

1−fx, t mx, t rx, t sx, t K

>000 0.

2.6

This is clearly a contradiction. Thus even at an interior minimumf >0, hencef >0 everywhere else. The same argument can be applied on the equations describing them,r, andsvariables. Actually the equation forris exactly solvable and is seen to be positive. Thus our assumption via2.3is feasible. Thus we proceed with our proof via maximum principle.

Despite not biologically viable, assume for purposes of analysis that

fK≥1, mK≥1. 2.7

We now define the positive and negative parts off−Kas fK

x

⎧⎨

fK, f > K, 0, otherwise, fK

x

⎧⎨

fK, f < K, 0, otherwise.

2.8

We now multiply1.1byf−Kxand integrate by parts to yield d

dtfK

2

2fK

2

2δfK

2

2

ΩF

f, m, r, s fK

xdx. 2.9 Whenf < Kthe right-hand side is zero. Whenf > K, assumingfKwhere >0, and m > kvia2.3, we have

ΩF

f, m, r, s fK

xdx

Ωfm

1−fmrs K

fK

xdx

ΩKK

1−fmrs K

dx

ΩKK

1−2K2δ K

dx

≤ |Ω|K2

−1−2δ K

≤ |Ω|K2

−2δ K

≤0.

2.10

(6)

Hence, via Poincar´e’s Inequality, we obtain d

dtfK

2

2 CδfK

2

2≤0. 2.11

Application of Gronwall’s Lemma now yields fK

2

2e−Cδtf0K

2

2. 2.12

We can now considert → ∞to yield

fK

0. 2.13

The same argument on the negative part offyields, fK

0. 2.14

Since the positive and negative parts offcan be no more thanK, we obtain f

K. 2.15

The same technique works onm, s, andrand is trivially seen to be bounded from the form of1.4.

3. A Priori Estimates

3.1. A Priori Estimates forfn

In order to prove the well posedness we follow the standard approach of projecting onto a finite dimensional subspace. This reduces the PDE to a finite dimensional system of ODE’s.

It is on this truncated system that we make a priori estimates. Essentially The truncation for ftakes the form

fnt n

j 1

fnjtwj. 3.1

Here wj are the eigenfunctions of the negative Laplacian, so −Δwi λiwi. A similar truncation can be performed for m, r and s. Thus, essentially the following holds for all 1≤jn,

∂fn

∂t DΔfnPn

F

fn, mn, rn, sn

δfn, 3.2

fn0 Pn

f0

. 3.3

(7)

HerePnis the projection onto the space of the firstneigenvectors. Note in general fn, Pn

F

fn

Pn fn

, F

fn fn, F

fn

. 3.4

We multiply3.2byfnand integrate by parts overΩ. We thus obtain 1

2 dfn2

2

dt −D∇fn2

2β

2 Ωmnfn2dx

Ωmnfn2fnmnrnsn

K dx

δfn2

2. 3.5

Via the positivity offn,mn,rn,sn, andKit follows that

Ωmnfn2fnmnrnsn

K dx

Ωmnfn2fn

Kdx. 3.6

This estimate is used in3.5to yield 1

2 dfn2

2

dt D∇fn2

2δfn2

2 β

2K

Ωmnfn3dxβ 2

Ωmnfn2dx. 3.7 We now use Young’s Inequality to obtain

1 2

dfn2

2

dt D∇fn2

2δfn2

2 β

2K

Ωmnfn3dxβ 2K

Ωmnfn3dxβK2 2

Ωmndx. 3.8 Using

|mn|≤ |m|K, 3.9

we obtain the following 1 2

dfn2

2

dt D∇fn2

2δfn2

2βK3

2 |Ω|. 3.10

The use of Poincar´e’s Inequality yields dfn2

2

dt CDδfn2

2βK3|Ω|. 3.11

Now, we can apply Gronwall’s Lemma to yield fnt2

2e−CDδtf02

2 βK3|Ω|

CDδC, ∀t≥0. 3.12

(8)

On the other hand we can integrate3.10from 0 toTto obtain 1

2fnT2

2D T

0

∇fn2

2dtδ T

0

fn2

2dtT

0

βK3|Ω|dtfn02

2. 3.13

This immediately yields T

0

∇fn2

2dtT

0

βK3|Ω|dtfn02

2

T

0

βK3|Ω|dtf02

2

C.

3.14

Thus, via3.12and3.14, we obtain fnL

0, T;L2Ω

, 3.15

fnL2

0, T;H01Ω

. 3.16

3.2. Estimate for the Time Derivative offn

We multiply3.2by awH01Ωto yield ∂fn

∂t , w

−D

∇fn,∇w

F

fn, mn, rn, sn

, Pnw

δ fn, w

. 3.17

We estimate the nonlinear term as follows:

F fn

, Pnw

Ωmnfn

1−fnmnrnsn K

Pnwdx

ΩmnfnPnwdx

≤ |mn|

ΩfnPnwdx

Kfn

4|Pnw|4/3

Cfn

4|w|H1

0.

3.18

This follows via the compact embedding ofH01ΩL4/3Ω. Thus, we have ∂fn

∂t 2

H−1Ωfn2

4. 3.19

(9)

Integrating both sides of the above in the time interval0, Tyields T

0

∂fn

∂t 2

H−1ΩdtT

0

fn2

4dtC T

0

∇fn2

2dtC. 3.20

This follows from the derived estimate via3.16and the compact embedding ofH01Ω L4Ω. Thus, we obtain

∂fn

∂tL2

0, T;H−1Ω

. 3.21

We can now via3.15and3.16extract a subsequencefnjsuch that

fnj f inL

0, T;L2Ω , fnj f inL2

0, T;H01Ω , fnj −→f inL2

0, T;L2Ω .

3.22

The convergence in the last equation follows via the compact embedding ofH01ΩL2Ω.

3.3. A Priori Estimates form,r, and s

The a priori estimates form,randsare very similar to the estimates forf. We omit the details here and present the results.

The truncation formsatisfies the following a priori estimates:

mnL

0, T;L2Ω , mnL2

0, T;H01Ω ,

∂sn

∂tL2

0, T;H−1Ω .

3.23

We can now extract a subsequencemnj such that

mnj s inL

0, T;L2Ω , mnj s inL2

0, T;H01Ω , mnj −→s inL2

0, T;L2Ω .

3.24

(10)

The last inequality follows via the compact embedding of

H01ΩL2Ω. 3.25

The truncation forssatisfies the following a priori estimates:

snL

0, T;L2Ω , snL2

0, T;H01Ω ,

∂sn

∂tL2

0, T;H−1Ω .

3.26

We can now extract a subsequencesnjsuch that

snj

s inL

0, T;L2Ω , snj s inL2

0, T;H01Ω , snj −→s inL2

0, T;L2Ω .

3.27

The last inequality follows via the compact embedding of

H01ΩL2Ω. 3.28

The truncation forrsatisfies the following a priori estimates:

rnL

0, T;L2Ω , rnL2

0, T;H01Ω ,

∂rn

∂tL2

0, T;H−1Ω .

3.29

We can now extract a subsequencernjsuch that

rnj r inL

0, T;L2Ω , rnj r inL2

0, T;H01Ω , rnj−→r inL2

0, T;L2Ω .

3.30

(11)

The last inequality follows via the compact embedding of

H01ΩL2Ω. 3.31

4. Existence of Solution

4.1. Preliminaries

We recast1.1in the following form:

∂f

∂t DΔfF

f, r, m, s

δf, f

∂Ω 0. 4.1

Here,

F

f, m, r, s β 2

1−fmrs K

fm. 4.2

Note that the key element in proving the existence will be to show convergence of the nonlinear termFfnj, mnj, rnj, snjtoFf, m, r, s. To this end we state the following lemma, Lemma 4.1. Consider the non linear termsFf1, m1, r1, s1andFf2, m2, r2, s2as defined via4.2.

The following estimate for their difference holds F

f1, m1, r1, s1

F

f2, m2, r2, s2

2Cf1f2

2|m1m2|2|s1s2|2|r1r2|2 . 4.3

Proof. Via4.2, we have that F

f1, m2, r2, s2

F

f2, m2, r2, s2 f1m1f2m2

f12m1f22m2

m21f1m22f2 f1m1r1f2m2r2f1m1s1f2m2s2

f1m1m2 m2 f1f2

f12m1m2 m2

f12f22

m21

f1f2 f2

m21m22

m1r1 f1f2

f2m2r1r2 f2r1m1m2 m1s1

f1f2

f2m2s1s2 f2s1m1m2.

4.4

We supress the dependence of the right-hand side on the constantβ/2 for convenience.

(12)

This follows from standard algebraic manipulation. Application of Holder’s and Minkowski’s inequalities yield

F

f1, m2, r2, s2

F

f2, m2, r2, s2

2

f1

|m1m2|2|m2|f1f2

2f12

|m1m2|2 |m2|f1f2

f1f2

2|m1|2f1f2

2

f2

|m1m2||m1m2|2 |m1||r1|f1f2

2f2

|m2||r1r2|2|r1|f2

|m1m2|2 |s1||m1|f1f2

2f2

|m2||s1s2|2|s1|f2

|m1m2|2

Cf1f2

2|m1m2|2|s1s2|2|r1r2|2 .

4.5

4.2. Passage to Weak Limit

As, we have made the a priori estimates on the truncations, we will attempt to pass to the weak limit, as is the standard practice. We will focus on1.1. Recall via Galerkin truncation we are seeking an approximate solution of the form

fnt n

j 1

fnjtwj, 4.6

such that, for each 1≤jn, and for allφC0 0, T, the following holds:

dfn dt, φwj

D

∇fnj,∇wjφt δ

fnj, φtwj

F

fnj , φwj

, 4.7

fn0 Pn f0

. 4.8

Here and henceforth we assume Ffnj PnFfnj, wherePn is the projection operator onto the firstneigenvectors. Upon passage to the weak limit of4.7, we will have obtained

df dt, wj

D

∇f,∇wj

δ

f, wj F

f , wj

. 4.9

This will imply the existence of a weak solutionfto1.1. We proceed as follows. Consider a φC0 0, T. We multiply4.7byφtand integrate by parts in time to yield

T

0

fnj, φtwj

dt −D T

0

∇fnj,∇wjφt dt

T

0

F

fnj , φtwj

dt

δ T

0

fnj, φtwj

dt.

4.10

We will first show convergence of the nonlinear term. This is stated via the following lemma.

(13)

Lemma 4.2. Consider the nonlinear termFf, m, r, sas defined via4.2. The following convergence result holds:

jlim→ ∞

T

0

ΩF

fnj, mnj, snj, rnj

φtwjdxdt T

0

ΩF

f, m, s, r

φtwjdxdt, 4.11

for allφC0 0, T. Proof. Consider

lim

j→ ∞

T

0

ΩF fnj

φtwjdxdtT

0

ΩF f

φtwjdxdt

C T

0

Ω

F

fnj

, φwj

F

f

, φwj2dxdt

wj

T

0

Ω

F f

F

fnj2dxdt

C T

0

ffnj2

2m−mnj2

2r−rnj2

2s−snj2

2

dt

C ffnj

L20,T;L2mmnj

L20,T;L2ssnj

L20,T;L2

C rrnj

L20,T;L2

C0000 0.

4.12

This follows viaLemma 4.1and because, we have demonstrated fnj −→f inL2

0, T;L2Ω , mnj −→m inL2

0, T;L2Ω , snj −→s inL2

0, T;L2Ω , rnj−→r inL2

0, T;L2Ω .

4.13

(14)

Thus the convergence of the nonlinear term has been established. Now, taking the limit asj → ∞in4.10, we obtain

jlim→ ∞

T

0

fnj, φtwj

dtD T

0

∇fnj,∇wjφt dt

δ T

0

fnj, φwj dt

T

0

F

fnj , φwj

dt T

0

f, φtwj

dtD T

0

∇f,∇wjφ dt

δ T

0

f, φwj dt

T

0

F f

, φwj dt 0.

4.14

The last term on the right-hand side can be bounded as follows T

0

F f

, φwj

dt

T

0

Ω

f2wjdt

wj

2f

L20,T;L4Ω

wj

H01Ωf

L20,T;H01Ω

Cwj

H10Ω.

4.15

This follows by the compact embedding ofH01Ω L4Ω L2Ω. This implies that, we have continuity with respect towj. Thus, we obtain that for anyvH01Ω the following holds

T

0

f, φtv dtD

T

0

∇f,∇vφt dtδ

T

0

f, φtv dt

T

0

F f

, φtv

dt. 4.16

This yields the existence of anfsuch that the following is true in a distributional sense d

dt f, v

D

∇f,∇v δ

f, v F

f , v

, ∀v∈H01Ω. 4.17

In other words there exists a weak solutionfto1.1. Since fL

0, T;L2Ω

L2

0, T;H01Ω ,

∂f

∂tL2

0, T;H−1Ω ,

4.18

(15)

it follows via standard PDE theory, see16,17, that fC

0, T;L2Ω

. 4.19

This establishes that the solution belongs to the requisite functional spaces.

4.3. Continuity with Respect to Initial Data and Uniqueness of Solutions We now show continuity with respect to initial data of the solution via the following lemma.

Lemma 4.3. Consider the Trojan Y Chromosome model. For positive initial data inL2Ω, any weak solutionf, m, s, rof the Trojan Y Chromosome model is continuous with respect to initial data, that is,

f0 f0, m0 m0, s0 s0, r0 r0. 4.20

Proof. We will show the details forf, and the other variables follow suit accordingly. We take a test functionφC10, Tsuch that

φ0 1, φT 0. 4.21

With this choice ofφtin4.17, we integrate the first term twice by parts to yield

T

0

f, φtv dtD

T

0

∇f,∇vφt dtδ

T

0

f, φtv dt f0, v

D T

0

∇f,∇vφt dtδ

T

0

f, φtv dt.

4.22

Note that the truncation satisfies T

0

fnj, φtv dtD

T

0

∇fnj,∇vφt dtδ

T

0

fnj, φtv dt

fnj0, v D

T

0

∇fnj,∇vφt dtδ

T

0

fnj, φtv dt.

4.23

Thus, taking the limit asj → ∞in4.28just as done earlier yields

T

0

f, φtv dtD

T

0

∇f,∇vφt dtδ

T

0

f, φtv dt f0, v

D T

0

∇f,∇vφt dtδ

T

0

f, φtv dt.

4.24

(16)

Thus, we obtain

f0, v f0, v

, ∀v∈H01Ω. 4.25

This yields

f0 f0, 4.26

as is required.

We now state the uniqueness result via the following lemma.

Lemma 4.4. Consider the Trojan Y Chromosome model. For positive initial data inL2Ωany weak solutionf, m, s, rof the Trojan Y Chromosome model is unique.

Proof. We work out the case for thefvariable, uniqueness for the others follow similarly. We consider the difference of two solutionsf1andf2to1.1. We denote

w f1f2, 4.27

andwsatisfies the following equation:

dw

dtDΔwδw F f1

F f2

, 4.28

w0 f10−f20 0. 4.29

We can multiply4.28bywand integrate by parts overΩto yield d|w|22

dt D|∇w|22δ|w|22

Ω

F f1

F f2

wdx. 4.30

Via the uniformL2estimates onm, r, s, see14, andLemma 4.1, we obtain d|w|22

dt D|∇w|22δ|w|22Cf1f2

2|w|2CK|w|22. 4.31 This yields

d|w|22

dt D|∇w|22δ|w|22−CK|w|22≤0. 4.32 Now using Poincar´e’s Inequality, we obtain,

d|w|22

dt C|w|22≤0. 4.33

(17)

The use of Gronwall’s Lemma yields that for anyt >0 the following estimate holds:

|wt|22e−Dδ−Ct|w0|22≤0. 4.34

Equation4.17in conjunction withLemma 4.4yieldsTheorem 1.1.

5. Weighted Sobolev Spaces

The purpose of this section is to introduce weighted Sobolev spaces into the framework of our present problem. We will show thatr given by1.4, remains bounded in the norms of these spaces. This will enable us to state a theorem about the existence of weak solution in the weighted spaces. This in turn will entail making refined estimates on the dimension of the global attractor for TYC system, when the phase space is a weighted Sobolev space. This will be achieved via the elegant technique of projecting the trace operator onto a weighted Sobolev space. We first make certain requisite definitions.

Definition 5.1. The weighted Sobolev spaceWωxk,p , with weight functionωx, is defined to be the space consisting of all functionsusuch that

|α|≤k

Ω|Dαu|pωxdx

1/p

<∞. 5.1

Remark 5.2. Here,Dα is theαth weak derivative ofu. In particular, we are interested in the following spaces for our application:

L2ωΩ

u:

Ωωx|u|2dx 1/2

<

,

H0,ω1 Ω

u:|u|2,ω|∇u|2,ω<.

5.2

Also, we denote

Ωωx|u|2dx1/2 |u|2,ω. We defineHω−1Ωto be the dual ofH0,ω1 Ω.

5.1. Estimates forr in Weighted Sobolev Spaces Recall the equation forr

∂r

∂t DΔrδrμ, r|∂Ω 0. 5.3

(18)

We chooseωx eμx,μ >0, multiply5.3byreμx, and integrate by parts overΩto yield

1 2

d dt

Ω|r|2eμxdx −D

Ω|∇r|2eμxdxD

Ω∇r· ∇reμxdx−δ

Ω|r|2eμxdx μ

Ωreμxdx

≤ −D

Ω|∇r|2eμxdx D 2

Ω|∇r|2eμxdxμ2 2

Ω|r|2eμxdx

δ

Ω|r|2eμxdxμ

Ωreμxdx

≤ −D 2

Ω|∇r|2eμxdxδ

Ω|r|2eμxdxC μ2K2

2 μK |Ω|.

5.4

These follow via integration by parts, the estimate |r|K, and the Cauchy-Schwartz inequality. Thus, we obtain

1 2

d

dt|r|22,ωD

2|∇r|22,ωδ|r|22,ωC μ2K2

2 μK |Ω|. 5.5

The use of Poincaire’s Inequality gives us 1

2 d|r|22,ω

dt

D 2 δ

|∇r|22,ωC μ2K2

2 μK |Ω|. 5.6

Now, we can apply the Gronwall Lemma to yield

|rt|22,ωe−CDδt|r0|22,ωμ2K2/2μK

CDδ , ∀t≥0. 5.7

On the other hand we can integrate5.5from 0 toT to obtain 1

2|rT|22,ωD 2

T

0

|∇r|22,ωdtδ T

0

|r|22,ωdtT

0

μ2K2

2 μK |Ω|dt. 5.8

This immediately yields T

0

|∇r|22,ωdtT

0

μ2K2

2 μK |Ω|dt. 5.9

(19)

Thus, via5.7and5.9, we have that

|r|L0,T;L2ωΩC <∞, 5.10

|r|L20,T;H10,ωΩC <∞. 5.11

5.2. Estimate for the Time Derivative ofr in Weighted Sobolev Space We multiply5.3by awH0,ω1 Ωto yield

∂r

∂t, w

2,ω

−D∇r,∇w2,ωδr, w2,ω w, μ

2,ω, ∂r

∂t

Hω−1Ωμ|w|2,ω.

5.12

Integrating both sides in time from 0 toTyields T

0

∂r

∂t 2

Hω−1Ωdtμ T

0

|w|22,ω

dt. 5.13

Because of the estimate via5.11and the embedding of

H0,ω1 ΩL2ωΩ, 5.14

we have

∂r

∂tL2

0, T;Hω−1Ω

< C <∞. 5.15

Thus it follows via the standard functional analysis theory, see16, that rC

0, T;L2ωΩ

. 5.16

These estimates show that r remains bounded in the appropriate weighted spaces introduced earlier and thus enables us to state the following theorem.

Theorem 5.3. Consider1.4in the TYC system. For positiver0L2ωΩ, there exists a unique weak solutionrto the system with

rC

0, T;L2ωΩ

L

0, T;L2ωΩ

L2

0, T;H0,ω1 Ω ,

∂r

∂tL2

0, T;Hω−1Ω .

5.17

Furthermore, the solutions are continuous with respect to initial data.

(20)

The uniqueness and convergence result by mimicking the method of proof for Theorem 1.1.

6. Existence of Global Attractor in Weighted Sobolev Space

We recall the following spaces from14, as the natural phase space for our problem:

H L2Ω×L2Ω×L2Ω×L2Ω, Y H01Ω×H01Ω×H01Ω×H01Ω, X H2Ω×H2Ω×H2Ω×H2Ω.

6.1

We next state the following definition.

Definition 6.1. Consider a semigroupStacting on a phase spaceM, then the global attractor A ⊂Mfor this semigroup is an object that satisfies

iAis compact inM.

iiAis invariant, that is,StA A, t≥0.

iiiIfBis bounded inM, then

distMStB,A−→0, t−→ ∞. 6.2

We showed in 14 that there exists a H, X global attractor for the TYC system.

That is an attractor that is compactX, and attracts bounded subsets inHin theXtopology.

Furthermore we showed this attractor had finite fractal and Hausdorffdimension. Our goal now is to improve these estimates, on a somewhat different attractor, via the technique of weighted Sobolev spaces. To this end we define

H! L2Ω×L2Ω×L2Ω×L2ωΩ,

"

Y H01Ω×H01Ω×H01Ω×H0,ω1 Ω. 6.3

Hereωis the weight as introduced earlier. We will first demonstrate the existence of aH,! H! attractor for the TYC system. We will then provide estimates for its Hausdorffand fractal dimensions. The following proposition is stated next.

Proposition 6.2. Consider the TYC system,1.1–1.4. There exists aH,!H! global attractorA"

for the this system which is compact and invariant inH and attracts bounded subsets of" H!in theH! metric.

The proof follows readily by applying the techniques of14to the weighted spaces in question. Recall that there are two essential ingredients to show the existence of a global attractor. The existence of a bounded absorbing set and the asymptotic compactness of the semigroup, see18. Thus we will just focus onr, as the proof for the other variables is the

(21)

same as in14. We will prove the above proposition via two lemmas. The first of these is stated next.

Lemma 6.3. Consider the equation forr,1.4, in the TYC system. Forr0L2ωΩthere exists a bounded absorbing set forrinL2ωΩ.

Proof. Recall, via5.7, we have

|rt|22,ωe−CDδt|r0|22,ωμ2K2/2μK

CDδ , ∀t≥0, 6.4

Now consider a timet1such that

t1 max

⎜⎝0,ln

|r0|22,ω CDδ

⎟⎠. 6.5

It follows that for any timet > t1the following uniform estimate holds

|rt|22,ω≤1μ2K2/2μK

CDδC. 6.6

This gives us a bounded absorbing set forrinL2ωΩ.

We next state the following lemma.

Lemma 6.4. The semigroupStfor the TYC system,1.1–1.4, is asymptotically compact inH.! Proof. We again demonstrate the proof forr. Multiply5.3by−Δreμxand integrate by parts overΩto yield

1 2

d dt

Ω|∇r|2eμxdx≤ −D

Ω|Δr|2eμxdxδ

Ω∇r· ∇reμxdx μ

Ω

∇r∂r

∂teμx dx.

6.7

Now Poincaire’s Inequality along with Cauchy-Schwartz imply that 1

2 d dt

Ω|∇r|2eμxdxCDδ

Ω|∇r|2eμxdxC

|∇r|22 ∂r

∂t 2

2

. 6.8

However directly from1.4and the compact Sobolev embedding of

H2ΩH01Ω. 6.9

参照

関連したドキュメント

Many interesting graphs are obtained from combining pairs (or more) of graphs or operating on a single graph in some way. We now discuss a number of operations which are used

Key words: Benjamin-Ono equation, time local well-posedness, smoothing effect.. ∗ Faculty of Education and Culture, Miyazaki University, Nishi 1-1, Gakuen kiharudai, Miyazaki

It is suggested by our method that most of the quadratic algebras for all St¨ ackel equivalence classes of 3D second order quantum superintegrable systems on conformally flat

We present a Sobolev gradient type preconditioning for iterative methods used in solving second order semilinear elliptic systems; the n-tuple of independent Laplacians acts as

This paper is devoted to the investigation of the global asymptotic stability properties of switched systems subject to internal constant point delays, while the matrices defining

Next, we prove bounds for the dimensions of p-adic MLV-spaces in Section 3, assuming results in Section 4, and make a conjecture about a special element in the motivic Galois group

Transirico, “Second order elliptic equations in weighted Sobolev spaces on unbounded domains,” Rendiconti della Accademia Nazionale delle Scienze detta dei XL.. Memorie di

Then it follows immediately from a suitable version of “Hensel’s Lemma” [cf., e.g., the argument of [4], Lemma 2.1] that S may be obtained, as the notation suggests, as the m A