一類可用Hamilton-Jacobi方法求解的非保守Hamilton系統?

王勇梅鳳翔肖靜郭永新

1)(北京理工大學宇航學院,北京 100081)

2)(廣東醫科大學信息工程學院,東莞 523808)

3)(遼寧大學物理學院,沈陽 110036)

4)(遼東學院影像物理教研室,丹東 118001)

一類可用Hamilton-Jacobi方法求解的非保守Hamilton系統?

王勇1)2)梅鳳翔1)肖靜2)郭永新3)4)?

1)(北京理工大學宇航學院,北京 100081)

2)(廣東醫科大學信息工程學院,東莞 523808)

3)(遼寧大學物理學院,沈陽 110036)

4)(遼東學院影像物理教研室,丹東 118001)

(2016年8月18日收到;2016年12月3日收到修改稿)

Hamilton-Jacobi方法通常被認為是求解完整保守Hamilton系統正則方程的重要手段,但通過現代微分幾何理論發現,這種方法的適用范圍不僅僅局限于完整保守的Hamilton系統.根據Hamilton-Jacobi理論,證明了經典Hamilton-Jacobi方法可以被推廣至一類特殊的非保守Hamilton系統,即如果非保守Hamilton系統受到非保守力,則該系統的Hamilton正則方程也可以用Hamilton-Jacobi方法求解;對于這類非保守Hamilton系統,只要能夠找到其對應的Hamilton-Jacobi方程的一個完全解,就可以得到系統正則方程的全部第一積分.經典的Hamilton-Jacobi方法則是上述方法的一個特例.

Hamilton-Jacobi理論,非保守Hamilton系統,Hamilton正則方程

1 引 言

Hamilton-Jacobi方程是一類重要的非線性偏微分方程,在經典力學、幾何光學、流體力學、粒子物理、廣義相對論、量子力學、宇宙學、最優控制、化學等諸多研究領域都有重要的應用[1?11].在經典力學中,經典Hamilton-Jacobi方法是求解完整保守Hamilton系統正則方程的重要手段[9?11],這種方法有其獨到的優點,阿諾德曾經指出:“用Jacobi解出的許多問題用別的方法是解不出來的.”根據經典Hamilton-Jacobi方法,設系統的Hamilton函數為H(t,qi,pi)(本文中所有i,j=1,2,···,n, 且采用愛因斯坦求和約定),如果可以找到Hamilton-Jacobi方程

的一個完全解S(t,qi,αi)(其中αi為任意常數),則完整保守Hamilton系統的正則方程

就能找到全部2n個第一積分,即

式中βi為n個新的任意常數.

歷史上,經典Hamilton-Jacobi方法直接源自于人們對一階偏微分方程和常微分方程組之間密切聯系的認識[11].從偏微分方程的角度看,經典Hamilton-Jacobi方法的關鍵是指出了正則方程((2)式)的積分曲線與Hamilton-Jacobi方程((1)式)的特征曲線之間的關系,因此根據偏微分方程理論,只要找到一個完全解,就可以用代數的方法解出(1)式的特征曲線,從而得到完整、保守Hamilton系統的2n個第一積分.

近些年,隨著幾何力學的發展[12?15],人們從經典Hamilton-Jacobi方法出發,逐步構建了一套有著深刻物理和數學內涵,并和Lagrange力學、Hamilton力學相對獨立的Hamilton-Jacobi理論體系[14],用現代微分幾何理論從不同角度對Hamilton-Jacobi理論進行了深入研究[14?21],深刻揭示了Hamilton-Jacobi理論的幾何內涵.其中,用微分幾何語言可以將經典Hamilton-Jacobi方法表述為如下命題[16].

命題1設S=S(t,qi)是一個定義在U×I上的作用量場(U是位形空間M的一個開集,I= [a,b]?R). 如果S是(1)式的一個解,則存在一個可微映射φ2:N→T?M×R,它定義了Hamilton形式ω=dpi∧dqi?dH∧dt的一個積分子流形.反之,若且H1(U×I)=0(其中H1表示U×I的第一上同調群),則存在一個作用量場S,它滿足Hamilton-Jacobi方程,使得成立(θ=pidqi?Hdt).

在最近的研究中發現,命題1不僅適用于形如(1)式的經典Hamilton-Jacobi方程,而且可以完全類似地推廣到形如

的更一般的偏微分方程的討論中,這說明有可能可以把求解完整保守Hamilton系統的經典Hamilton-Jacobi方法推廣,應用于求解更一般的Hamilton系統的運動問題中.本文將證明可以用Hamilton-Jacobi方法求解受到一個非保守力Fi=μ(t)pi(其中μ(t)是一個時間因子)的非保守Hamilton系統的運動問題,并給出兩個算例.

2 受到非保守力Fi=μ(t)pi的非保守Hamilton系統的Hamilton-Jacobi方法

若一個Hamilton系統受到一個非保守力Fi=μ(t)pi,則系統的Lagrange方程為

對應的Hamilton方程為

此類問題的求解也可以化為求解一個一階偏微分方程完全解的問題,即有如下命題.

命題2設力學系統的Hamilton函數為H(t,qi,pi),構造如下一階偏微分方程:

式中bi為n個任意常數,即可得到(6)式所對應的包含2n個任意常數ai和bi的2n個第一積分.

要證明上述命題,首先可以直接驗證(6)式確實是(7)式的特征線方程的一部分,因此只需驗證由(7)式和(8)式可以推出(6)式即可.

由(8)式中第一個等式可得

將(9)式對時間求導可得

將(7)式分別關于ai和an+1求偏導數并考慮(8)式中第二個等式后可得

聯立(8)式、(10)式和(11)式,消去bi和μ(t)后得

只要行列式

即可由(12)式得到(6)式的第一個等式.

再由(8)式中第二個等式對時間求偏導可得

將(7)式對qi求偏導數并考慮(8)式中第二個等式后可得

聯立(14)式和(15)式后即得(6)式的第二個等式.證畢.

到此已經證明了(8)式確實是(6)式所對應的2n個第一積分.需要說明的是,從物理意義的角度考慮,上述2n個第一積分中只能包含2n個任意常數,即n個bi和(ai,an+1)中的n個任意常數,說明(3)式所求出的完全解S(t,qi,ai,an+1)至少關于其中一個常數是線性的,在本文中將這個常數記作an+1.

可以看出,經典Hamilton-Jacobi方法實際上是命題2在μ(t)=0時的特例,此時Hamilton-Jacobi方程的解具有的形式,(8)式就退化為經典Hamilton-Jacobi方法中的(3)式.

此外,還可以證明如下命題.

命題3可以用形如(4)式的一階非線性Hamilton-Jacobi方程求解非保守Hamilton系統只有命題2中的情況.

要證明命題3,假設非保守Hamilton系統所受的非保守力為Fi,記T?M和T?×R為系統的余切叢和擴展后的余切叢,對應的局部坐標系分別為和則系統在其余切叢T?M上的矢量場為

將形如(4)式的一階非線性Hamilton-Jacobi方程看作流形T?×R中的約束子流形

則該約束子流形在流形T??M×R中定義了一個由

所張成的可積一維余分布??,與(4)式的特征線所對應的特征矢量場

屬于與一維余分布??所對偶的可積分布?.由一階偏微分方程的幾何理論可知,若切叢T?M上的矢量場(16)式可用形如(4)式的一階非線性Hamilton-Jacobi方程求解,則必須存在一個從T?M到T?×R上的映射φ,使得矢量場(16)式的推前

(其中f4和f5為與映射φ有關的兩個函數)等于Hamilton-Jacobi方程所對應的特征矢量場(19)式,即必須滿足

3 算 例

算例1沿直線做一維運動的質量為1的粒子,所受外力為F=μp,其中μ為不為零的常數,求粒子的運動方程.

對應的偏微分方程為

設S(t,q)=W(t)+Q(q),則可解得(24)式的一個完全解為

式中a1和a2是兩個任意常數.令

可得兩個第一積分

式中b是任意常數.

本例也可以對(23)式直接求解,求得的結果為

考慮到式中c1和c2為任意常數,(28)式和(27)式等價,說明本例用Hamilton-Jacobi方法求解的結果是正確的.

算例2系統的Hamilton函數為sin2q·p2,若初始時刻q(0)∈(0,π),求:1)系統不受外力時的自由運動;2)系統所受外力為F=μ(t)p時的運動.

1)系統不受外力時

對應的偏微分方程為

設S(t,q)=W(t)+Q(q),則可解得(30)式的一個完全解為

式中a1和a2是兩個任意常數.令

式中b是任意常數,可得兩個第一積分:

2)當系統受到外力F=μ(t)p時

系統的Hamilton方程為

對應的偏微分方程為

設S(t,q)=W(t)+Q(q),則可解得(35)式的一個完全解為

式中a1和a2是兩個任意常數.令

式中b是任意常數,可得兩個第一積分

4 結 論

本文給出了求解受到一個非保守力Fi=μ(t)pi的非保守Hamilton系統運動方程的Hamilton-Jacobi方法,這種方法是Hamilton-Jacobi理論在非保守Hamilton系統中的一個應用.從幾何角度看,Hamilton-Jacobi方法本質上是通過Hamilton-Jacobi方程把系統Hamilton方程所描述的余切叢T?M上的矢量場提升到流形T?M×R中的一個約束子流形上,只要設法求出提升后的矢量場的積分曲線,則該積分曲線在余切叢T?M上的投影就是系統Hamilton方程的解.證明上述方法可以應用于求解受到一個非保守力Fi=μ(t)pi的Hamilton系統的運動問題.經典的Hamilton-Jacobi方法是這一方法的一個特例.

[1]Benamou J 1996J.Comput.Phys.128 463

[2]Fleming W H,Rishel R 1975Deterministic and Stochastic Optimal Control(Berlin:Spinger)pp80–105

[3]Feng C J,Wang P,Wang X M 2015Acta Phys.Sin.64 030502(in Chinese)[封晨潔,王鵬,王旭明2015物理學報64 030502]

[4]Fedkiw R P,Aslam T,Merrima B,Osher S 1999J.Comput.Phys.152 457

[5]Yang S Z,Lin K 2010Sci.China40 507(in Chinese)[楊樹政,林愷2010中國科學40 507]

[6]Yang S Z,Lin K 2013Acta Phys.Sin.62 060401(in Chinese)[楊樹政,林愷 2013物理學報62 060401]

[7]Kim J H,Lee H W 2000Can.J.Phys.77 411

[8]Joulin G,Mitani T 1981Comb.Flame.40 235

[9]Arnold V I.1978Mathematical Methods of Classical Mechanics(New York:Spriner-Verlag)pp161–271

[10]Mei F X 2013Analytical Mechanics(Vol.1)(Beijing:Beijing Institute of Technology Press)pp272–287(in Chinese)[梅鳳翔 2013分析力學 (上冊)(北京:北京理工大學出版社)第272—287頁]

[11]Courant R,Hilbert D 1989Methods of Mathematical Physics(Vol.2)(New York:John Wiley&Sons)pp62–153

[12]Guo Y X,Luo S K,Mei F X 2004Adv.Mech.34 477(in Chinese)[郭永新,羅紹凱,梅鳳翔 2004力學進展 34 477]

[13]Guo Y X,Liu S X,Liu C,Luo S K,Wang Y 2007J.Math.Phys.48 082901

[14]Marmo G,Morandi G,Mukunda N 1990La Rivista del Nuovo Cimento13 1

[15]Wang H 2013 arXiv:1305.3457v2[math.SG]

[16]Westenholtz C N 1981Differential Forms in Mathematical Physics(Amsterdam:North-Horland Publishing Company)pp389–439

[17]Barbero-Lin?án M,de León M,Martin de Diego D 2012Monatsh.Math.171 269

[18]Marmo G,Morandi G,Mukunda N 2009J.Geom.Mech.1 317

[19]Vitagliano L 2012Int.J.Geom.Methods Mod.Phys.9 1260008

[20]de León M,Vilari?o S 2014Int.J.Geom.Methods Mod.Phys.11 1450007

[21]Ohsawa T,Bloch A M 2009J.Geom.Mech.1 461

PACS:45.20.Jj,02.40.Yy DOI:10.7498/aps.66.054501

A kind of non-conservative Hamilton system solved by the Hamilton-Jacobi method?

Wang Yong1)2)Mei Feng-Xiang1)Xiao Jing2)Guo Yong-Xin3)4)?

1)(School of Aerospace Engineering,Beijing Institute of Technology,Beijing 100081,China)
2)(School of Information Engineering,Guangdong Medical University,Dongguan 523808,China)
3)(College of Physics,Liaoning University,Shenyang 110036,China)
4)(Department of Medical Imaging Physics,Eastern Liaoning University,Dandong 118001,China)

18 August 2016;revised manuscript

3 December 2016)

The Hamilton-Jacobi equation is an important nonlinear partial differential equation.In particular,the classical Hamilton-Jacobi method is generally considered to be an important means to solve the holonomic conservative dynamics problems in classical dynamics.According to the classical Hamilton-Jacobi theory,the classical Hamilton-Jacobi equation corresponds to the canonical Hamilton equations of the holonomic conservative dynamics system.If the complete solution of the classical Hamilton-Jacobi equation can be found,the solution of the canonical Hamilton equations can be found by the algebraic method.From the point of geometry view,the essential of the Hamilton-Jacobi method is that the Hamilton-Jacobi equation promotes the vector field on the cotangent bundleT?Mto a constraint submanifold of the manifoldT?M×R,and if the integral curve of the promoted vector field can be found,the projection of the integral curve in the cotangent bundleT?Mis the solution of the Hamilton equations.According to the geometric theory of thefirst order partial differential equations,the Hamilton-Jacobi method may be regarded as the study of the characteristic curves which generate the integral manifolds of the Hamilton 2-formω.This means that there is a duality relationship between the Hamilton-Jacobi equation and the canonical Hamilton equations.So if an action field,defined onU×I(Uis an open set of the configuration manifoldM,I?R),is a solution of the Hamilton-Jacobi equation,then there will exist a differentiable mapφfromM×RtoT?M×Rwhich defines an integral submanifold for the Hamilton 2-formω.Conversely,ifφ?ω=0 andH1(U×I)=0(H1(U×I)is the first de Rham group ofU×I),there will exist an action fieldSsatisfying the Hamilton-Jacobi equation.Obviously,the above mentioned geometric theory can not only be applicable to the classical Hamilton-Jacobi equation,but also to the general Hamilton-Jacobi equation,in which some first order partial differential equations correspond to the non-conservative Hamiltonian systems.The geometry theory of the Hamilton-Jacobi method is applied to some special non-conservative Hamiltonian systems,and a new Hamilton-Jacobi method is established.The Hamilton canonical equations of the non-conservative Hamiltonian systems which are applied with non-conservative forceFi=μ(t)pican be solved with the new method.If a complete solution of the corresponding Hamilton-Jacobi equation can be found,all the first integrals of the non-conservative Hamiltonian system will be found.The classical Hamilton-Jacobi method is a special case of the new Hamilton-Jacobi method.Some examples are constructed to illustrate the proposed method.

Hamilton-Jacobi theory,non-conservative Hamilton system,Hamilton canonical equation

PACS:45.20.Jj,02.40.Yy

10.7498/aps.66.054501

?國家自然科學基金(批準號:11572145,11272050,11572034)和廣東省自然科學基金(批準號:2015A030310127)資助的課題.

?通信作者.E-mail:yxguo@lnu.edu.cn

*Project supported by the National Natural Science Foundation of China(Grant Nos.11572145,11272050,11572034)and the Natural Science Foundation of Guangdong Province,China(Grant No.2015A030310127).

?Corresponding author.E-mail:yxguo@lnu.edu.cn