?? ?? = Δ?? の解は、熱核
② H 型群上の半線形熱方程式が実際にどのような
自然現象(社会現象)を記述しているのかを解の数値実験
などを行い探求する。(今後)
釧路高専 創造工学科 一般教育部門
2017年11月24日(金)
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 1 / 35
1 Y. Oka,Local well-posedness for semilinear heat equations on H type group, submitted to Taiwanese Journal of Mathematics.
2 A. Kaplan,Fundamental solutions for a class of hypoelliptic PDE generated by composition of quadratic forms, Trans. Amer. Math.
Soc., 258, (1980), 147-153.
3 H. Brezis and T. Cazenave,A nonlinear heat equation with singular initial date, J. Anal. Math. 68 (1996),277-304.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 2 / 35
初期値問題
(∂t−∆Hd p
)
u(g,t) =|u|r−1u, g∈Hdp,t >0,
u(g,0) =u0(g)∈Lq(Hdp). (1) の解の一意存在を考察する。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 3 / 35
初期値問題(1)をデュアメルの原理を用いて、積分方程式 u(t) =et∆Hdpu0+
∫ t
0
e(t−σ)∆Hdp|u(σ)|r−1u(σ)dσ の問題に書き換えて考察する。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 4 / 35
(z,s)◦(z′,s′) = (
z+z′,s+s′+1 2[z,z′]
)
, (2)
z = (x,y), z′ = (x′,y′)∈R2d,s,s′∈Rp を持つR2d+pである。 ここで,[z,z′]j=⟨
z,Ujz′⟩
,j =1,2,· · ·,pである. ただし,Ujは次の性質を満たすものである.
1 Uj,j =1,· · ·,pは2d×2dの交代かつ直交行列である.
2 すべてのi,j ∈ {1,2,· · ·,p},i ,jに対し,UiUj+UjUi =0が成り立つ. このとき,集合GをHdp と表すこととする。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 5 / 35
注意
H型の群に関して,中心の次元とその直交補空間の次元との関係は, p+1≤2d (A. Kaplan and F. Ricci, 1983)
行列Uiの条件が、線形独立な交代行列の場合、カルノー(Carnot)群と 呼ばれる.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 6 / 35
Xj = ∂
∂xj
+1 2
∑p
k=1
∑2d l=1
zlUlk,j
∂
∂sk
Yj = ∂
∂yj
+1 2
∑p
k=1
∑2d l=1
zlUlk,j+d
∂
∂sk,
ここで,zl=xl,zl+d =yl(l =1,2,· · ·,d)であり,Uik,j,Uki,j+dは行列Ukの (i,j)成分や(i,j+d)成分を表す.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 7 / 35
U1 = 0 −1 1 0 とすると,
[z,z′] =⟨
z,U1z′⟩
=
⟨(x y )
,
(0 −1
1 0 ) (x′
y′ )⟩
=
⟨(x y )
, (−y′
x′ )⟩
=−xy′+yx′ となるので, (2)より群Gの積の法は,
(z,s)◦(z′,s′) = (
z+z′,s+s′+1
2(yx′−xy′) )
, ここで,z= (x,y),z′ = (x′,y′)∈R2,s∈R1である.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 8 / 35
また,左不変ベクトル場は, X = ∂
∂x +1 2y ∂
∂s, Y = ∂
∂y −1 2x ∂
∂s, S = ∂
∂s. となる.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 9 / 35
∂x 2 ∂s ∂y 2 ∂s ∂y 2 ∂s ∂x 2 ∂s
= ∂2
∂x∂yu−1 2
∂
∂s
∂
∂x(xu) + 1 2y ∂2
∂y∂su− 1 4xy ∂2
∂s2u
− ∂2
∂x∂yu− 1 2
∂
∂s
∂
∂y(yu) +1 2x ∂2
∂x∂su+1 4xy ∂2
∂s2u
=−1 2
∂
∂s (
u+x ∂
∂xu )
−1 2
∂
∂s (
u+y ∂
∂yu )
+1 2y ∂2
∂y∂su+1 2x ∂2
∂x∂su
=− ∂
∂su
=−Su
∴ [X,Y] =−S, ([X,X] = [Y,Y] = [S,S] = [X,S] = [Y,S] =0)
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 10 / 35
x
y
P A
B
s方向の測地線はらせん
xy平面上の測地線は直線 ρ(x,y,s) ={
(x2+y2)2+s2}14 ハイゼンベルグ群上の距離関数 O
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 11 / 35
(4×4)の交代行列U1を,
U1 =
0 −1 0 0
1 0 0 0
0 0 0 −1
0 0 1 0
, U2 =
0 0 1 0
0 0 0 −1
−1 0 0 0
0 1 0 0
とすると, (2)より,群Gの積の法は,
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 12 / 35
z+z′ s1+s1′ +1
2
(−x1x2′+x2x1′−y1y2′+y2y1′)
s2+s2′ +1 2
(x1y1′−x2y2′ −x1′y1+y2x2′)
,
ここで,z= (x1,x2,y1,y2),z′= (x1′,x2′,y1′,y2′)∈z⊥,
s= (s1,s2),s′ = (s1′,s2′)∈zである. また,左不変ベクトル場は,
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 13 / 35
X2 = ∂
∂x2
+1 2 (
−x1 ∂
∂s1
+y2 ∂
∂s2
) ,
S1 = ∂
∂y1 − 1 2 (
y2 ∂
∂s1
+x1 ∂
∂s2
) ,
S2 = ∂
∂y2 − 1 2 (
−y1 ∂
∂s1 −x2 ∂
∂s2
)
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 14 / 35
● ∆Hd
p =
∑2d i=1
Xi2をサブラプラシアンという。
例:H11の場合は,∆H1
1 =X2+Y2となる。
● N =2d+2pを等質次元という。
例:H11の場合は,N=2×2+2×1=4
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 15 / 35
u(g,0) =u0(g)∈L (Hp).
(定理1 時間局所解の存在性) Y. Oka
q>N(r−1)/2かつq≥1,N≥4とする。このとき、すべてのu0 ∈Lq(Hdp) に対し、正定数T と(1)の解u∈C([0,T];Lq(Hdp))が存在する。
(定理2 時間局所解の存在性) Y. Oka
q=N(r −1)/2かつq>1,N≥4とする。このとき、すべてのu0 ∈Lq(Hdp) に対し、正定数T と(1)の解u∈C([0,T];Lq(Hdp))が存在する。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 16 / 35
q>N(r−1)/2(resp.q=N(r−1)/2)かつq≥r(resp.q>r),N ≥4を仮 定する。このとき、解
u(t) =et∆u0+
∫ t
0
e(t−σ)∆|u(σ)|r−1u(σ)dσ の一意性が空間C([0,T];Lq(Hdp))で成り立つ。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 17 / 35
に対し次の評価式が成り立つ。
∫
Hdp
|f(g)h(g)|dg≤ ∥f∥Lq(Hdp)∥g∥Lq′(Hdp).
(補題2)ヤングの不等式
1≤q≤ ∞とする。このとき、f ∈L1(Hdp)とh∈Lq(Hdp)に対し次の評価式 が成り立つ。
∥f∗h∥Lq(Hdp) ≤ ∥f∥L1(Hdp)∥h∥Lq(Hdp).
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 18 / 35
の解u(x,t)は、熱核htとの畳込みで次の様に表せる:
u(x,t) =et∆Hdpu0=
∫
Hdp
u0(g′)ht(g′−1g)dg′
(補題3)熱核の評価(D. S. Jerison and A. S«anchez-Calle 1986)
∆Hd
pに付随する熱核をht(g)とする。このとき、∆Hd
p に依存した正定数C1 とCI,lが存在して次が成り立つ。
|∂ltXIht(g)| ≤CI,lt−l−|2I|−N2e−c1ρ
(g)2 t ,
ここで、I= (i1,· · ·,im),|I|=mそしてXI =Xi1Xi2· · ·Xim. MoreoverρはH 型群のCaratheodory distanceである。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 19 / 35
∥e φ∥Lβ(Hdp) ≤Ct ∥φ∥Lα(Hdp), t >0, がすべてのφ∈Lα(Hdp)に対して成り立つ。
●Lq−Lm 評価として補題4を使うと、
∥e(t−σ)∆Hdp(|u(σ)|r−1u(σ))∥Lm(Hdp)≤C(t−σ)−N2(1q−m1)∥|u(σ)|r−1u(σ)∥Lq(Hdp)
=C(t−σ)−N2(1q−m1)
∫
Hdp
|u(g, σ)|qrdg
1 qr
r
=C(t−σ)−N2(1q−m1)∥u(σ)∥rLqr(Hdp).
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 20 / 35
Hp Hp
が成り立つ。一方で、半群et∆Hdp のL2有界性より
∥et∆Hdpφ∥L2(Hdp) ≤C∥φ∥L2(Hdp). が成り立つ。Riesz-Thorin補完定理によって、
∥et∆Hdpφ∥Lβ(Hdp) ≤Ct−N2(1α−1β)∥φ∥Lα(Hdp).
が成り立つ。空間S(Hdp)はLα(Hdp)で稠密なので、φ∈Lα(Hdp)においてこ
の評価式が成り立つ。 □
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 21 / 35
T :X →X が,
∃k ∈(0,1)s.t.d(Tx,Ty)≤kd(x,y), x,y ∈X (3) を満たすとき、T はXにおいてただ1つの不動点x0、すなわちTx0 =x0と なる点をもつ。
※ 条件(3)を満たす写像T を縮小写像という。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 22 / 35
定理1の証明の流れ
1 関数空間YT とBM+1の設定する。
2 次で定義される写像Φ:BM+1→YT
Φ[u](t) =et∆Hdpu0+
∫ t 0
e(t−σ)∆Hdp(|u(σ)|r−1u(σ))dσ, がBM+1からBM+1への写像であることを示す。
3 写像ΦがBM+1からBM+1への縮小写像であることを示す。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 23 / 35
∥u∥YT =max sup
0<t<T∥u(t)∥Lq(Hdp), sup
0<t<T
tλ∥u(t)∥Lqr(Hdp) , λ= N(r−1) 2qr <1 で定義する。次に、YT の部分集合として空間BM+1を
BM+1 ={u| ∥u∥YT ≤M+1} とする。ここで、M=max{M1,M2} s.t.
∥u0∥YT ≤M1and∥et∆Hdpu0∥YT(≤C∥u0∥YT)≤M2. とする。Mは∥u0∥YT に依存し、t には独立である。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 24 / 35
最初に写像ΦがBM+1からBM+1への中への写像であることを示す。補題 4(Lq−Lm 評価)とq>N(r −1)/2によって,すべてのu∈BM+1に対し、
∫ t
0
e(t−σ)∆Hdp(|u(σ)|r−1u(σ))dσ Lm(Hdp)
≤C
∫ t
0
(t−σ)−N2(1q−m1)∥u(σ)∥rLqr(Hdp)dσ
≤C(M+1)r
∫ t
0
(t−σ)−N2(1q−m1)σ−rλ dσ
=C(M+1)rt1−rλ−N2(1q−m1)
∫ 1
0
(1−σ)−N2(q1−m1)σ−rλdσ がm≥qにおいて成り立つ。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 25 / 35
tN2(q1−m1) ∫ t
0
e(t−σ)∆Hdp(|u(σ)|r−1u(σ))dσ Lm(Hdp)
≤C(M+1)rT1−rλ. (4) となる。いま、(4)内でm=qもしくはm=qrとすると,次が成り立つ。
∫ t
0
e(t−σ)∆Hdp(|u(σ)|r−1u(σ))dσ Lq(Hdp)
≤C1(M+1)rT1−rλ または
tλ ∫ t
0
e(t−σ)∆Hdp(|u(σ)|r−1u(σ))dσ Lqr(Hdp)
≤C2(M+1)rT1−rλ.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 26 / 35
となる。十分小さいT >0に対して、
max{C1,C2}(M+1)rT1−rλ ≤1.
が成り立つ。よって、写像ΦはBM+1からBM+1への写像である。
次に写像ΦがBM+1からYT への縮小写像であることを示す。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 27 / 35
と補題4(L −L 評価)によって、u1,u2∈BM+1に対し、次が成り立つ。
∥Φ[u1](t)−Φ[u2](t)∥Lm(Hdp)≤C3(M+1)r−1t1−rλ−N2(q1−m1)∥u1−u2∥YT (6) ここで、定数C3 >0.いま(6)においてm=qもしくはm=qrとする。
このとき、次が成り立つ。
∥Φ[u1](t)−Φ[u2](t)∥Lq(Hdp)≤C3(M+1)r−1T1−rλ∥u1−u2∥YT
または
tλ∥Φ[u1](t)−Φ[u2](t)∥Lqr(Hdp) ≤C3(M+1)r−1T1−rλ∥u1−u2∥YT
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 28 / 35
∥Φ[u1](t)−Φ[u2](t)∥YT ≤C4(M+1)r−1T1−rλ∥u1−u2∥YT
ここで、正定数C4 >0. 1−rλ >0なので十分小さいT >0に対し、
C4(M+1)r−1T1−rλ ≤ 1 2.
よって、写像Φは十分小さいT >0に対し縮小写像になる。 従って、縮小 写像の原理(バナッハの不動点定理)より、 写像Φの不動点uがBM+1内 に 存在することが示された。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 29 / 35
最後に、u∈C([0,T];Lq(Hdp))を示す。実際、u∈BM+1かつrλ <1より
|u|r−1u∈L1((0,T);Lq(Hdp)).がいえる。これから、u∈C([0,T];Lq(Hdp))が
導かれる。これにて定理1の証明は終了。 □
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 30 / 35
Gronwallの補題を用いて、定理1で得られた解の初期値連続性
∥u(t)−v(t)∥Lq(Hdp)≤C∥u0−v0∥Lq(Hdp), a.a.t ∈[0,T]. も示すことができる。
定理2に関しては、技術的な難しさはあるが、定理2とほぼ同様に証 明できる。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 31 / 35
u(t) =et∆u0+
∫ t
0
e(t−σ)∆|u(σ)|r−1u(σ)dσ の一意性が空間C([0,T];Lq(Hdp))で成り立つ。
(Caseq>N(r−1)/2): Letu,v∈C([0,T];Lq(Hdp))be two solutions. Then we have
u(t)−v(t) =
∫ t 0
e(t−σ)∆(|u(σ)|r−1u(σ)− |v(σ)|r−1v(σ))dσ.
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 32 / 35
L (Hp) L (Hp) L (Hp)
によって、次が成り立つ。
∥u(t)−v(t)∥Lq(Hdp)≤C
∫ t 0
(t −σ)−θ∥|u|r−1u− |v|r−1v∥Lqr(Hd p)dσ
≤C′
∫ t 0
(t −σ)−θ
(∥u∥rL−q1(Hdp)+∥v∥rL−q1(Hdp)
)∥u−v∥Lq(Hdp),
(7) ここで、θ=N(r −1)/2q<1かつ正定数C, C′はt に独立である。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 33 / 35
ψ(t) = sup
0≤σ≤t∥u(t)−v(t)∥Lq(Hdp)とする。このとき(7)によって、次が成り 立つ。
ψ(t)≤CMr−1T1−θ
1−θψ(t). (8)
いま、T′は0<T′<T となる十分小さいとし、t ∈[0,T′]とする。 このと き、(8)によって、ψ(t) =0が導ける。 同様なことを有限回繰り返し, t ∈[0,T]においてψ(t) =0. □
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 34 / 35
(∂t−∆H1 1
)
u(g,t) =|u|r−1u, g∈H11,t >0, u(g,0) =u0(g)
の解uは局所可解である。ここで、H11(=R3)は1次元のハイゼンベルグ群 を表す。つまり,初期値問題
(
∂t−∂2x−∂2y− 1
4(x2+y2)∂2s−(
y∂x−x∂y
)∂s
)
u(x,y,s,t) =|u|r−1u, u(x,y,s,0) =u0(x,y,s)
の解uは(x,y,s)∈H11かつt >0において局所可解である。
岡(釧路高専) 第2章 H型群に付随する半線形熱方程式 2017年11月24日(金) 35 / 35