• 沒有找到結果。

非線性擾動下三維有界域半線性波方式爆炸解之穩定性研究

N/A
N/A
Protected

Academic year: 2021

Share "非線性擾動下三維有界域半線性波方式爆炸解之穩定性研究"

Copied!
12
0
0

加載中.... (立即查看全文)

全文

(1)

行政院國家科學委員會專題研究計畫 成果報告

非線性擾動下三維有界域半線性波方式爆炸解之穩定性研

研究成果報告(精簡版)

計 畫 類 別 : 個別型

計 畫 編 號 : NSC 99-2115-M-004-001-

執 行 期 間 : 99 年 08 月 01 日至 100 年 07 月 31 日

執 行 單 位 : 國立政治大學應用數學學系

計 畫 主 持 人 : 李明融

共 同 主 持 人 : 謝宗翰、白仁德

報 告 附 件 : 國外研究心得報告

處 理 方 式 : 本計畫涉及專利或其他智慧財產權,2 年後可公開查詢

中 華 民 國 100 年 09 月 23 日

(2)

Stability of positive solutions for some semilinear wave

equations under nonlinear perturbation near blow-up

solutions in 3-space dimension

u

u

p+

u

q =

0

Meng-Rong Li

Department of Mathematical Sciences National Chengchi University

Abstract In this research we treat the stability of positive solutions of some particular semilinear wave equations under nonlinear perturbation in bounded do-main near blow-up solutions in 3-space dimension.

1.

Introduction

Consider the initial value problem for the semilinear wave equation of the type u + g (u) = 0 in [0; T ) R3;

(0.1)

u (0; ) = u0; _u (0; ) = u1;

(0.2)

where g : R ! R is a real valued function, the initial data are given su¢ ciently smooth functions and u := utt 4u; 4 is the Laplace operator. The linear

case g(u) = mu, where m is a constant, corresponds to the classical Klein Gordon equation in relativistic particle physics; the constant m is interpreted as the mass and is assumed to be nonnegative generally. To model also nonlinear phenomena like quantization, in the 1950s equations of (0.1) type with nonlinearities like g(u) = mu + u3; m 0; were proposed as models in relativistic quantum mechanics with

local interaction. Solutions could be considered as real or complex valued functions. In the latter case it was assumed that the nonlinearity commutes with the phase; that is, g ei'u = "ei'g (u) for ' 2 R and that g(0) = 0. In this case, g may be expressed g(u) = uf juj2 , which gives the study of equation (0:1) [J ]. In the

(3)

noncoercive case it is easy to construct solutions of (0:1) with smooth initial data that blow up in …nite time; for instance, for any > 0 the function u(t; x) = (1 t) 1=m solves the equation u + (1 + )u juj2m

= 0; m 2 N and blows up at t = 1. Modifying the initial data o¤ fx : jxj 2g, say, we even posses a singular solution with C1-data having compact support.

In this study we want to deal with the stability of positive solutions for the semilinear wave equation

(1:1) u = up+ uq in [0; T ) ; R3

with boundary value null and initial values u (0; ) = u0( ) 2 H2( ) \ H01( ) and

_u (0; ) = u1( ) 2 H01( ) ; where p; q 2 (1; 1) and is a bounded domain in R3:

We will use the following notations: := @ @t; Du := ( _u; ru) ; 4u := @2u @2x 1 + @ 2u @2x 2 + @ 2u @2x 3 ; a (t) := Z u2(t; x) dx; E (t) := Z jDuj2 2 p + 1u p+1 2 q + 1u q+1 (t; x) dx:

For a Banach space X and 0 < T 1 we set

Ck(0; T; X) = Space of Ck functions : [0; T ) ! X;

H1 := C1 0; T; H01( ) \ C2 0; T; L2( ) :

The existence result to the equation (1:1) is proved [Li 3] and the positive solu-tion blows-up in …nite time if 0 [Li 2] ; this means that the positive solutions for the semilinear wave equation

u = up in [0; T ) ; (1:2)

u (0; ) = u0( ) 2 H2( ) \ H01( ) ;

_u (0; ) = u1( ) 2 H01( ) ;

is stable under nonlinear perturbation uq providing p > 1; q > 1; > 0; but it is

not clearly whether it is also true for any p > 1; q > 1; < 0 ? If so, we would want to estimate the blow-up time and the blow-up rate under such a situation.

It is also important to study the asymptotic behavior of the solution u ; the velocity and the rate of the approximation for approaches to zero.

Such questions are also not easy to answer even under the case for the ordinary di¤erential equation

u00= up(c + u0(t)q) ; (1.3)

u (0) = u0; u0(0) = u1;

where p > 1; q > 1; c > 0; > 0: We have studied the blow-up behavior of the solution for problem (1:3) and got some estimates on blow-up time and blow-up rate [Li4] but it is di¢ cult to …nd the real blow-up time (life-span). Further literature could be fund in [S], [R], [W1] and [W2].

(4)

3

In this study we hope that our ideals used in [Li 2]; [Li 4]; [Li 5],[Li 7] ; [Li8]; [LiLinShieh]; [ShiehLi] and [SLLLW ] can do help us dealing such problem (1:1) on our topics.

2. De…nition and Fundamental Lemma

There are many de…nitions of the weak solutions of the initial-boundary problems of the wave equation, we use here as following.

De…nition 2.1: For p > 1; u 2 H1 is called a positive weakly solution of equation (1:1), if Z t 0 Z ( _u (r; x) _' (r; x) + (up+ uq) (r; x) ' (r; x)) dxdr = Z t 0 Z ru (r; x) r' (r; x) dxdr 8' 2 H1 and Z t 0 Z u (r; x) (r; x) dxdr 0 for each positive 2 C1

0 ([0; T ) ).

Remark 2.2:

1) The de…nition 1.1 is resulted from the multiplying with ' to the equation (1:1) and integrating in from 0 to t.

2) From the local Lipschitz functions up+ uq; p > 1; q > 1 the initial-boundary value problem (1:1) possesses a unique solution in H1 [Li1]. Hereafter we use the notations: 1 C := 1= sup kuk2= @u @x 2: u 2 H 1 0( ) ; q = sup kukq= @u @x 2: u 2 H 1 0( ) \ Lq( ) ; q > 1:

In this study we need the following lemmas

Lemma 2.3:

Suppose that u 2 H1 is a weakly positive solution of (1:1) with

E (0) = 0 for p > 1; q > 1; then for a (0) > 0 we have: (i) a 2 C2(R+) and E (t) = E (0) 8t 2 [0; T ).

(ii) a0(t) > 0 8t 2 [0; T ), provided a0(0) > 0.

(iii) a0(t) > 0 8t 2 (0; T ), if a0(0) = 0.

(iv) For a0(0) < 0; there exists a constant t0> 0 with

a0(t) > 0 8t > t0

and a0(t 0) = 0:

(5)

Lemma 2.4

:

Suppose that u is a positive weakly solution in H1 of equation (1:1) with u (0; ) = 0 = _u (0; ) in L2( ). For p > 1; q > 1; > 0; we have u 0

in H1.

According to Lemma 2.4, we discuss the following theme (3) E (0) = 0; a (0) > 0 and a0(0) 0 or a0(0) < 0. (4) E (0) < 0; a (0) > 0 and a0(0) 0 or a0(0) < 0.

3. Estimates for the Life-Span

3.1. Estimates for the Life-Span of the Solutions of (1.1) under Null-Energy. We study the case that E (0) = 0, p > 1; q > 1; > 0 and divide it into two parts

(i) a (0) > 0; a0(0) 0 and (ii) a (0) > 0; a0(0) < 0:

Remark 3.. 1) The local existence and uniqueness of solutions of equation (1:1) in H1 are known [Li 2].

2) For = 0, x p > 1 and E (0) = 0; the life-span of the positive solution u 2 H1 of equation (1.1) is bounded by T 1:= k21sin 1 k2 k1a p 1 4 (0) ! with k1:= p 1 4 a p 1 4 (0) q a0(0) a 2(0) + 4C2; k 2:= p 1 2 C ; 1

C := 1= sup kuk2= kDuk2: u 2 H

1 0( ) :

3.2. Estimates for the Life-Span of the Solutions of equation (1.1) under Negativ-Energy. We use the following result and those argumentations of proof are not true for positive energy, so under positive energy we need another method to show the similar results.

Lemma 3

: Suppose that u 2 H1 is a positive weakly solution of equation (1.1)

with a (0) > 0 and E (0) < 0 for = 0. Then (i) for a0(0) 0, we have a0(t) > 0 8t > 0:

(ii) for a0(0) < 0; there exists a constant t5> 0 with

a0(t) > 0 8t > t5; a0(t5) = 0

and

t5 t6:=

a0(0)

(p 1) ( 2 E (0));

where is the positive root of the equation 2

p + 1

p+1

(6)

5

4. Stability of positive solutions of equation (1.1) near blow-up solutions under Negativ-Energy

In this study we use our ideals used in [Li 2]; [Li 4]; [Li 5],[Li 7] ; [Li8]; [LiLinShieh]; [ShiehLi] and [SLLLW ] to deal such problem (1:1) on our topics under negative energy and

obtain the following results:

Theorem 4.1:

Suppose that u 2 H1 is a weakly positive solution of (1:1) with

E (0) 0 for p > 1; q > 1; then for a (0) > 0 we have:

The equation (1:1) is stable for ! 0+; this means that weakly positive solution u of (SL) blows up in …nite time for ! 0+:

Theorem 4.2:

Suppose that u 2 H1 is a weakly positive solution of (1:1) with

E (0) 0 for p > 1; q > 1; then for a (0) > 0 we have:

The equation (1:1) is stable for p > q; ! 0 ; this means that weakly positive solution u of (1:1) blows up in …nite time for p > q; ! 0 :

Theorem 4.3:

Suppose that u 2 H1 is a weakly positive solution of (1:1) with

E (0) 0 for p > 1; q > 1; then for a (0) > 0 we have:

The equation (1:1) is unstable for p < q; ! 0 ; this means that some weakly positive solution u of (1:1) blow up in …nite time for p < q; ! 0 ; but also there were some global weakly positive solution u of (1:1) for p < q; ! 0 :

Remark:

The decade rate of the di¤erence of life-spans T of u and T of u; can not be estimated very well for ! 0; thus it will be a good topic on asymptotic behavior near the blow-up solutions.

Reference:

[J] Jörgens,K.: Das Anfangswertproblem im Größen fur eine Klasse nichtlinearer Wellengleichungen. M. Z.77. pp.295-307 (1961)

[Li1] Meng-Rong Li: On the Semi-Linear Wave Equations (I) : Taiwanese Journal of Math. Vol. 2, No. 3, pp. 329-345, Sept. 1998

[Li2] Meng-Rong Li: Estimates for the Life-Span of the Solutions of some Semi-linear Wave Equations. Communications on Pure and Applied Analysis vol.7, no. 2, pp.417-432. (2008).

[Li3] Meng-Rong Li: Nichtlineare Wellengleichungen 2. Ordnung auf beschränk-ten Gebiebeschränk-ten. PhD-Dissertation Tübingen 1994.

[Li 4] Meng-Rong Li: Blow-up solutions to the nonlinear second order di¤erential equation u00= up(c1+ c2u0(t)q). Taiwanese Journal of Mathematics, vol.12, no.3,

pp.599-622, June 2008.

[Li5] Renjun Duan, Meng-Rong Li; Tong Yang: Propagation of Singularities in the Solutions to the Boltzmann Equation near Equilibrium, Mathematical Models and Methods in Applied Sciences (M3AS), vol.18, no.7, pp.1093-1114.(2008)

(7)

[Li6] Meng-Rong Li; Brain Pai: Quenching problem in some semilinear wave equations. Acta math. scientia vol.28, no.3, pp.523-529, July 2008.

[Li7] Meng-Rong Li: Estimates for the life-span of the solutions of some semilin-ear wave equation u up= 0 in one space dimension. Communications on Pure

and Applied Analysis, 2012 to appear.

[Li8 ] Meng-Rong Li; Chuen-Hsin Chang: Analysis on the Duration of Deep Earthquakes, Taiwanese Journal of Mathematics, Vol.12, No.7, Oct. 2008, pp.1667-1680.

[LiShiehYuLiLi] M.R. Li; T.H. Shieh; etc : Parabola Method in Ordinary Dif-ferential Equation, Taiwanese Journal of Mathematics, 2012 to appear.

[PaiLiChangChiu] J.T. Pai; Meng-Rong Li; Y.L. Chang; S.M. Chiu: A mathe-matical model of enterprise competitive abilitry and performance through a partic-ular Emden-Fowler equation ( to appear in Acta Mathematica scientia 2011,31B(5), SCIE.)

[LiLinShieh] Meng-Rong Li; Y.J. Lin; T.H. Shieh: The ‡ux model of the move-ment of tumor cells and health cells using a system of nonlinear heat equations. Journal of Computational Biology, January 2011, ahead of print.

[ShiehLi ] T.H. Shieh; Meng-Rong Li: Numeric treatment of contact disconti-nuity with multi-gases Journal of computational and applied mathematics (Journal of Computational and Applied Mathematics, 2009 August pp.656-673.

[SLLLW] T.H. Shieh; T.M. Liou; Meng-Rong Li; C.H. Liu; W.J. Wu: Analysis on numerical results for stage separation with di¤erent exhaust holes, International communications in heat and mass transfer, 2009 April pp.342-345.

[R] Racke R.: Lectures on nonlinear Evolution Equations: Initial value problems. Aspects of Math. Braunschweig Wiesbaden Vieweg(1992).

[S] Strauss W.A.: Nonlinear Wave Equations. A.M.S. Providence(1989). Di-mensions. J. Di¤erential Equations52.pp.378-406(1984).

[W1] von Wahl W.: Klassische Lösungen nichtlinearer Wellengleichungen im Großen.M.Z.112.pp.241-279(1969).

[W2] von Wahl W.: Klassische Lösungen nichtlinearer gedämpfter Wellengle-ichungen im Großen. Manuscripta Math.3. pp7-33(1970).

(8)

一月十七日(一) 與尤釋賢教授 討論 波方程式解的一些性質 (I)

一月十八日(二) 與尤釋賢教授 討論 波方程式解的一些性質 (II)

一月十九日(三) 與尤釋賢教授 討論 波方程式解的奇異性 (I)

一月二十日(四) 與尤釋賢教授 討論 波方程式解的奇異性 (II)

一月二十一日(五) 與尤釋賢教授 討論 波方程式解的奇異集 (I)

一月二十二日(六) 與尤釋賢教授 討論 波方程式解的奇異集 (II)

一月二十三日(日) 與尤釋賢教授 共同討論行進波與震波

之交互作用(I)

一月二十四日(一) 與尤釋賢教授 共同討論行進波與震波

之交互作用(II)

並贈語曰

星州晚暮

侵晨蘇醒夜朦朧 美人依舊臨虛空

埜花不知寄何處 星州大學綠萬叢

細索問題本原地 菩提知慧無底洞

醉月早忘千年愁 與君共享夕陽紅

神奇無雨

今日最神奇 此間無煙雨 禿筆忘畊雲 宿墨離塵汙

星州新夜月 行者問所須 勸君壹壺酒 千憂自茲去

庚寅冬禿筆風語于星州怡閣

(9)

國科會補助計畫衍生研發成果推廣資料表

日期:2011/09/20

國科會補助計畫

計畫名稱: 非線性擾動下三維有界域半線性波方式爆炸解之穩定性研究 計畫主持人: 李明融 計畫編號: 99-2115-M-004-001- 學門領域: 微分方程

無研發成果推廣資料

(10)

99 年度專題研究計畫研究成果彙整表

計畫主持人:

李明融

計畫編號:

99-2115-M-004-001-計畫名稱:

非線性擾動下三維有界域半線性波方式爆炸解之穩定性研究

量化

成果項目

實際已達成

數(被接受

或已發表)

預期總達成

數(含實際已

達成數)

本計畫實

際貢獻百

分比

單位

備 註

質 化 說

明:如 數 個 計 畫

共 同 成 果、成 果

列 為 該 期 刊 之

封 面 故 事 ...

期刊論文

1

1

100%

研究報告/技術報告

1

1

100%

研討會論文

0

0

100%

論文著作

專書

0

0

100%

申請中件數

0

0

100%

專利

已獲得件數

0

0

100%

件數

0

0

100%

技術移轉

權利金

0

0

100%

千元

碩士生

3

0

100%

博士生

2

0

100%

博士後研究員

0

0

100%

國內

參與計畫人力

(本國籍)

專任助理

0

0

100%

人次

期刊論文

2

2

100%

研究報告/技術報告

0

0

100%

研討會論文

0

0

100%

論文著作

專書

0

0

100%

章/本

申請中件數

0

0

100%

專利

已獲得件數

0

0

100%

件數

0

0

100%

技術移轉

權利金

0

0

100%

千元

碩士生

0

0

100%

博士生

0

0

100%

博士後研究員

0

0

100%

國外

參與計畫人力

(外國籍)

專任助理

0

0

100%

人次

(11)

其他成果

(

無法以量化表達之成

果如辦理學術活動、獲

得獎項、重要國際合

作、研究成果國際影響

力及其他協助產業技

術發展之具體效益事

項等,請以文字敘述填

列。)

成果項目

量化

名稱或內容性質簡述

測驗工具(含質性與量性)

0

課程/模組

0

電腦及網路系統或工具

0

教材

0

舉辦之活動/競賽

0

研討會/工作坊

0

電子報、網站

0

目 計畫成果推廣之參與(閱聽)人數

0

(12)

國科會補助專題研究計畫成果報告自評表

請就研究內容與原計畫相符程度、達成預期目標情況、研究成果之學術或應用價

值(簡要敘述成果所代表之意義、價值、影響或進一步發展之可能性)

、是否適

合在學術期刊發表或申請專利、主要發現或其他有關價值等,作一綜合評估。

1. 請就研究內容與原計畫相符程度、達成預期目標情況作一綜合評估

■達成目標

□未達成目標(請說明,以 100 字為限)

□實驗失敗

□因故實驗中斷

□其他原因

說明:

2. 研究成果在學術期刊發表或申請專利等情形:

論文:□已發表 □未發表之文稿 ■撰寫中 □無

專利:□已獲得 □申請中 ■無

技轉:□已技轉 □洽談中 ■無

其他:(以 100 字為限)

3. 請依學術成就、技術創新、社會影響等方面,評估研究成果之學術或應用價

值(簡要敘述成果所代表之意義、價值、影響或進一步發展之可能性)(以

500 字為限)

參考文獻

相關文件

The first row shows the eyespot with white inner ring, black middle ring, and yellow outer ring in Bicyclus anynana.. The second row provides the eyespot with black inner ring

When the spatial dimension is N = 2, we establish the De Giorgi type conjecture for the blow-up nonlinear elliptic system under suitable conditions at infinity on bound

Quadratically convergent sequences generally converge much more quickly thank those that converge only linearly.

“Find sufficiently accurate starting approximate solution by using Steepest Descent method” + ”Compute convergent solution by using Newton-based methods”. The method of

denote the successive intervals produced by the bisection algorithm... denote the successive intervals produced by the

Abstract In this paper, we consider the smoothing Newton method for solving a type of absolute value equations associated with second order cone (SOCAVE for short), which.. 1

In this paper we establish, by using the obtained second-order calculations and the recent results of [23], complete characterizations of full and tilt stability for locally

In this paper we establish, by using the obtained second-order calculations and the recent results of [25], complete characterizations of full and tilt stability for locally