Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,非线性动力学,非线性动力学,1,优选非线性动力学,优选非线性动力学,2,Beyond Perturbation,Introduction to Homotopy Analysis Method,Beyond Perturbation In,3,Outline,Concept of Homotopy in Topology,Basic ideas of Homotopy Analysis method,Examples,Applications of the theory in solving nonlinear equations,Conclusions,References,Outline Concept of Homotopy in,4,“摄动方法”的本质：,应用方程中的小（大）物理参数，将一个非线性问题转化为无穷多个线性子问题。,优点：物理意义明确；简单、易懂；,缺点:(1,)依赖小参数,当所研究问题不含小参 数时使得摄动展开法面临困难,(2)摄动展开解只在参数比较小的情况下能够给出较好的近似,随着“小参数”的增大,近似解精度下降,以致失效。,(3)无法确保解的收敛,“摄动方法”的本质：,5,怎样的近似解析方法,才是最理想的,？,不依赖小参数,确保解的收敛性，适用于强非线性问题,非线性动力学课件,6,拓扑学中的几个基本概念,拓扑,和,拓扑空间,如果对一个非空集合,X,给予适当的结构，使之能引入微积分中的极限和连续的概念，这样的结构就称为拓扑。,具有拓扑结构的空间称为拓扑空间。,引入拓扑结构的方法有多种，如邻域系、开集系、闭集系、闭包系、内部系等不同方法。,拓扑学中的几个基本概念拓扑和拓扑空间,7,同伦的基本概念,两个,拓扑空间,如果可以通过一系列,连续的形变,从一个变到另一个，那么就称这两个拓扑空间,同伦,。,同伦的基本概念,8,同伦的定义,设,X,和,Y,都是拓扑空间，,f,和,g,是,X,到,Y,的连续映 射，,即,f：XY,，,g： XY,，,如果存在连续映射,H：X,Y（,这里,=0,1），,使得对任何,xX,，,满足：,则称,f,和,g,是,同伦的,称,H,是由,f,到,g,的一个同伦,或伦移,即,H(x,0)=f(x)，H(x,1)=g(x)，,同伦是关于映射的等价关系,f（x） =,H（x,0）,g（x） =H（x,1）,H（x,q）,示意图,同伦的定义H(x,0)=f(x)，H(x,1)=g(x)，同,9,二、“同伦分析方法”简述,拓扑理论传统的同伦概念：,其中，,q,为嵌入变量,.,易知，,q=0,时,，,H,(,x;,0),=f,(,x,),;,q=1,时,，,H,(,x;,1),=,g,(,x,),.,因此，当嵌入变量,q,从0增加到1时，函数,H(x,q),从,f(x),连续变化到,g(x),.,这样，,H(x,t),建立起从,f(x),到和,g(x),之间的联系.在拓扑(,topology),理论中，这种连续的变化称为同伦(,homotopy)，,表示为,二、“同伦分析方法”简述拓扑理论传统的同伦概念：,10,Liao,提出“广义同伦”之概念：,非线性动力学课件,11,非线性动力学课件,12,Basic ideas of HAM,E,1,.,非线性代数方程,f,(,x,)=0.(,构造同伦),设 为已知的初始猜测解,嵌入变量,为一未知的嵌入变量 的函数,我们构造如下的一个单参数的非线性代数方程:,(1),当 时,上述方程为线性方程,即,Basic ideas of HAME1.非线性代数方程 f,13,当 时,方程(1)变为,则 , 就是原非线性方程,f,(,x,)=0,的解.,因此,当嵌入变量 从0变化到1时, 从初始猜测解 变化到非线性代数方程解,因此方程(1)构造了一个,的同伦.,设 存在无穷阶导数,当 时,方程(1)变为则,14,根据,Taylor,定理,有,则,如何求 ?,(2),将(1)式对,p,求一阶导数,(3),根据Taylor定理,有则如何求 ?(2)将(,15,令,得,则,将(3)式对,p,再求一次导数,(4),(5),令,得,(6),令 得则将(3)式对p再求,16,类似地，可以求得,k,阶变形导数 ，则,一阶近似公式为,（ 时为牛顿迭代公式）,类似地，可以求得k阶变形导数 ，则一阶近似公,17,E,2,.,非线性微分方程,where is a nonlinear operator, denotes independent variable, is an unknown function, respectively.,(1) Construct zero-order deformation equation,Where, ,0,1,is the embedding parameter,is a nonzero auxiliary parameter,is an auxiliary function,is an auxiliary linear operator,is an initial guess of ,is a unknown function, respectively.,(7),E2.非线性微分方程where is a non,18,Obviously, when,p,=,0 and,p,=,1, it holds,Thus as increases from 0 to 1, the solution varies from the initial guess to the solution .,Expanding in Taylor series with respect to ,one has,where,(8),Obviously, when p = 0 and p =1,19,If the auxiliary linear operator , the initial guess , the auxiliary parameter , and the auxiliary function are so properly chosen,the series (8) converges at , one has,which must be one of solutions of original nonlinear equation.,As and ,Eq(7) becomes,(9),which is used mostly in the homotopy analysis method.,(2) Construct mth-order deformation equation,If the auxiliary linear operat,20,Differentiating Eq. (7),m,times with respect to the embedding parameter,p,and then setting,p,=,0,and finally dividing them by,m,!, we have the so-called,m,th-order deformation equation,Define the vector,线性方程,(10),？,？,Differentiating Eq. (7) m tim,21,It should be emphasized that for,m,1,is governed by the linear equation (10) with the linear boundary conditions that come from original problem, which can be easily solved by symbolic computation software such as,Maple and Mathematica,.,It should be emphasized that,22,根据Taylor定理,有,的同伦.,（D）海洋工程中的应用, the solution of (11) must be expressed in the same form as (12) and the other expressions such as must be avoided.,Thus as increases from 0 to 1, the solution varies from the initial guess to the solution .,如何求 ?,同伦分析方法应用举例：非线性水波,Applications of the theory in solving nonlinear equations,(2) Construct mth-order deformation equation,三维非定常旋转黏性流动,which must be one of solutions of original nonlinear equation.,is an initial guess of ,E3.,非线性微分方程求解,According to the governing equation and the initial condition (11), the solution can be expressed by a set of base functions,(11),in the form,where is a coefficient to be determined，This provides us with the so-called,rule of solution expression, i.e., the solution of (11) must be expressed in the same form as (12) and the other expressions such as must be avoided.,(12),根据Taylor定理,有E3.非线性微分方程求解Accord,23,According to (11) and (12), we choose the linear operator,with the property,Where is constant.,From (11), we define a nonlinear operator,According to (11) and (12), we,24,According to (11) and the rule of solution expression (12), it is straightforward that the initial approximation should be in the form,(1) Construct zero-order deformation equation,Thus as increases from 0 to 1, the solution varies from the initial guess to the solution .,(2) Construct mth-order deformation equation,According to (11) and the rule,25,非线性动力学课件,26,最后得到,rule of coefficient ergodicity，H,（,）=1,最后得到 rule of coefficient ergo,27,得到一族解，通过 调节级数收敛,得到一族解，通过 调节级数收敛,28,二、“同伦分析方法”简述,“,同伦分析方法,”,特点,毋须任何小参数，可将一个非线性问题转化为无穷多个线性问题！,可自由选取辅助线性算子、初始近似:,线性子问题中的线性算子毋须与原始非线性方程中的线性算子相同或密切相关！,二、“同伦分析方法”简述“同伦分析方法”特点,29,二、“同伦分析方法”简述,初步形成一个较为完整的理论体系,（1）提出三个原则:,解表达原则（,Rule of solution expression）,解存在原则（,Rule of solution existence）,完备性原则（,Rule of coefficient ergodicity）,指导辅助线性算子、初始近似、辅助函数之选取,（2）证明了,“,收敛性定理,”,二、“同伦分析方法”简述初步形成一个较为完整的理论体系,30,同伦分析方法之优点,不同于摄动方法，“同伦分析方法”不依赖于小参数的存在，因而适用范围更广；,不同于所有其它分析方法，“同伦分析方法”本身提供了一种简单的方法调节或控制解析解级数的收敛区域；,“,同伦分析方法,”,提供选择不同基函数之自由，从而能更有效地表达非线性问题的解。,同伦分析方法之优点不同于摄动方法，“同伦分析方法”不依赖于小,31,二、“同伦分析方法”简述,广泛应用（1992年-2002年）,非线性波浪问题,边界层流动和热传导问题,非线性振动问题,极限环问题,圆球黏性阻力（,Navier-Stokes,方程）,物理、生物及宇宙学方面的非线性问题,证明“同伦分析方法”之有效性和潜力,二、“同伦分析方法”简述广泛应用（1992年-2002年）,32,( 1 )不依赖小参数,二阶近似在整个区间,内的最大误差仅为,0.48%,求解范例,( 1 )不依赖小参数求解范例,33,同伦分析方法之优点,圆球绕流问题,应用,“,同伦分析方法,”, 得到150年来与实验结果吻合得最好的圆球阻力理论公式,（2002年）,。,应用,“,同伦分析方法,”,求解一些经典非线性难题,同伦分析方法之优点圆球绕流问题,34,同伦分析方法之优点,( 2 ) 确保解的收敛性 解的收敛区域可以,调节和控制,同伦分析方法之优点( 2 ) 确保解的收敛性 解的收敛区,35,( 3 ) 有选择基函数之自由,对任何参数,我们都得到如下形式的周期解,( 3 ) 有选择基函数之自由,36,同伦分析方法之优点,Liao, S. and Tan, Y., “A general approach to obtain series solutions of nonlinear differential equations ”,Studies in Applied Mathematics, 119:1-58,2007.,同伦分析方法之优点 Liao, S. and Tan,37,非牛顿流体边界层流动,非牛顿流体边界层流动,38,非牛顿流体边界层流动,非牛顿流体边界层流动,39,三维非定常旋转黏性流动,三维非定常旋转黏性流动,40,三维非定常旋转黏性流动,Tan.Y and Liao, S. ,ASME J. Applied Mech.,74:1011-1018,2007,三维非定常旋转黏性流动 Tan.Y and Liao,41,（,B）,发现新解,( 1 ) 可渗透拉伸变形平板边界层流动：,（B）发现新解( 1 ) 可渗透拉伸变形平板边界层流动：,42,（,B）,发现新解,应用,“,同伦分析方法,”, 找到被数值方法遗漏的一个新解!,（B） 发现新解 应用 “同伦分析方法”, 找到被数值方,43,（,B）,发现新解,( 2),Cheng-Minkowycz,流动:,（B）发现新解( 2) Cheng-Minkowycz 流动,44,呈代数衰减的无穷多个解,应用 “同伦分析方法”,Liao and Magyari (2006),找到被数值方法遗漏的、,呈代数衰减的无穷多个新解！,呈代数衰减的无穷多个解 应用 “同伦分析方法”, Lia,45,（,C）,突破传统思想,（C） 突破传统思想,46,（,C）,突破传统思想,求解非线性问题时，我们所拥有的自由，远比我们过去想象的要大得多！,正面意义：提出更好的、求解非线性,问题的解析方法和数值方法,许多全新的、有趣的问题有待研究和探索,( 请见力学进展有关综述论文 ),（C） 突破传统思想 求解非线性问题时，我们所拥有的自由，,47,（,D）,海洋工程中的应用,“同伦分析方法”被成功应用于研究海洋工程中的一些基础理论问题，如：,非线性波浪；,梁的大扰度弯曲；,非线性波与非均匀流相互作用；,（D）海洋工程中的应用 “同伦分析方法”被成功应用于研究海,48,（,D）,海洋工程中的应用,（,A）,非线性,深水行进波,（D）海洋工程中的应用,49,（,D）,海洋工程中的应用,（,B）,具有间断性,的孤立波,（D）海洋工程中的应用（B）具有间断性,50,（,D）,海洋工程中的应用,（,C）,梁的大扰度弯曲,（D）海洋工程中的应用（C）梁的大扰度弯曲,51,（,D）,海洋工程中的应用,（,D）,非线性波浪与非均匀流相互作用,（D）海洋工程中的应用（D）非线性波浪与非均匀流相互作用,52,（,D）,海洋工程中的应用,非线性波浪与非均匀流相互作用,（D）海洋工程中的应用 非线性波浪与非均匀流相互作用,53,众多成功的应用实例,证实了“同伦分析方法”求解强非线性问题的有效性;,“同伦分析方法”能找到新的、甚至被数值方法忽略的解，说明了“同伦分析方法” 的巨大潜力；,众多成功的应用实例,证实了“同伦分析方法”求解强非线性问题的,54,意义,众所周知，流体力学和海洋工程中的非线性问题特别多。,“,同伦分析方法,”,的提出和完善，为流体力学和海洋工程中强非线性问题的求解提供了一个全新的、强有力的理论分析工具,。,意义 众所周知，流体力学和海洋工程中的非线性问题特别多。“,55,非线性动力学课件,56,同伦分析方法应用举例：非线性水波,同伦分析方法应用举例：非线性水波,57,同伦分析方法应用举例：,Riemman,问题,同伦分析方法应用举例：Riemman 问题,58,同伦分析方法应用举例：美式期权定价问题,同伦分析方法应用举例：美式期权定价问题,59,机遇和挑战,机遇,“,同伦分析方法,”,抛弃了小参数假设，为强非线性问题的求解提供了一个新的思路，为一些经典非线性难题的求解提供了一种新的可能性,挑战,自然界中的非线性现象千差万别，异常复杂。提出一种普遍有效的方法极为困难。必须不断地完善和改进,“,同伦分析方法,”,，使其适用于尽可能多的强非线性问题,机遇和挑战机遇挑战,60,身体健康，学习进步！,身体健康，学习进步！,