Hamilton力学的辛算法
- 格式:ppt
- 大小:2.26 MB
- 文档页数:44
哈密尔顿系统的辛几何算法哈密尔顿系统是一类具有特殊的物理意义的动力系统,其在物理学、力学、动力学和计算力学等领域有着广泛应用。
哈密尔顿系统通常具有一组关于位置和动量的相变量,其演化满足哈密尔顿方程。
由于哈密尔顿系统具有良好的保持量和结构稳定性,因此在数值模拟中的算法设计尤其重要。
辛几何算法是一类特殊的数值演化方法,其以保持哈密尔顿系统相变量守恒和辛结构稳定性为目标,常常用于哈密尔顿系统的数值积分。
辛几何算法最早由李约瑟于 1988 年提出,其不仅能够在数值计算中保持相变量的守恒,还能够在哈密尔顿系统的长期演化中保持辛结构稳定性。
辛几何算法主要由两个部分组成,即辛映射和辛算子。
辛映射指的是从一个相变量向下一个相变量的映射,它通常满足“保相量”和辛结构不变性的特点。
保相量指的是相变量在变化过程中的守恒,而辛结构不变性则指的是哈密尔顿系统在演化过程中的辛不变性。
而辛算子则是这个辛映射的数值逼近,常常采用辛波发方法、显式和隐式辛算法等方法。
在演化哈密尔顿系统时,辛几何算法通常采用显式辛算法进行数值模拟。
显式辛算法的主要思路是采用辛映射和辛算子的组合,来实现对哈密尔顿系统的数值模拟。
在模拟过程中,辛几何算法需要保证每一步的演化都是辛的,这样系统才能保持哈密尔顿量以及其他相变量不变。
因此,辛几何算法在数值模拟中的应用非常广泛。
然而,辛几何算法的实现却比较困难。
在数值模拟时,辛几何算法需要考虑一系列问题,如相变量的数值守恒、哈密尔顿量的捕获和重构、快速演化、长时间演化、难以计算的高维效应等等。
这些问题都需要采用一些特殊的技巧和策略来解决。
总的来说,辛几何算法是一种特殊的数值演化方法,其以保持哈密尔顿系统相变量守恒和辛结构稳定性为目标,在计算力学、物理学、动力学等领域有着广泛应用。
辛算法研究冯康教授于1984年提出Hamilton系统的辛几何算法,首次将保持Hamilton 系统几何结构的思想引入数值分析,随后引来了国内外在这方面的极大兴趣。
冯康教授开辟了一个新的研究领域,并在近十年的时间里带领他的研究小组在辛几何算法等保结构算法及其在数值分析中的应用等方面进行了广泛深入的研究,取得了丰硕的成果。
现在人们已越来越意识到保结构算法的重要性,事实上保结构算法已在很多领域包括天体力学,量子化学,非线性波,不可压流体,大气物理和地物勘探数据处理等,找到很好的应用。
另一方面,在这些领域中的应用反过来又必将促进辛几何算法等保结构算法本身的不断完善何不断发展。
大量数值实验结果已证明,较之传统算法,保结构算法对保结构动力系统的计算具有令人信服的优势。
在辛几何算法理论方面,仍有一些悬而未决的难点问题和最新提出的一些重点问题有待解决;另一方面,在辛算法的基础上,现在国际上又兴起了专门解保守型偏微分方程的多辛算法。
孤立子波动方程是广泛应用于物理领域的非常重要的守恒型偏微分方程,对其复杂孤立子波的数值模拟二十多年来一直是国际上一个热门课题,数值结果已经显示辛算法,多辛算法的优势和可行性,但至尽仍然存在许多有待解决的问题。
我们主要研究的是多辛算法,现阶段的研究水平已完全与国际科研水平持平,并保持同步发展。
我们把多辛几何算法应用到孤子方程中去,数值研究表明多辛算法和常微分方程情形一样具有长期跟踪能力, 不会带来人为污染, 能正确反应孤子碰撞问题.我们发表在Phys.A Mathematical gen(2000), 33:18, 3613--3626上的论文给出了一个三层12点格式. 此文发表后不久, 两美国学者Schultz M ,Trimper S 在此刊同卷41期上发表论文:“动力运动产生孤子”, 文中称我们的方法是著名的方法。
李群算法是最近才发展起来的一种很有发展前途的数值算法。
它用来求解齐次流形上的常微分方程,使得所求的数值解仍在同一流形上李群算法是冯康先生辛几何算法的拓扩,把保几何结构的思想推广到一般李群上。
辛算法在电磁计算中的应用摘要近几年,随着计算机性能的飞速发展和计算物理中各种新型算法的出现,各种电磁场数值方法层出不穷,但很多算法面临着计算时间长、储存空间不足及计算精度低等方面的困难。
Hamilton系统理论是当代数学物理中的一个重要的工具。
一切守恒的物理过程,总能表示成适当的Hamilton系统。
辛算法正是保持Hamilton系统内在性质的一种新型数值方法,该算法在长时间的数值计算中,具有一般数值方法无可比拟的计算优势。
本文首先介绍了电磁学的基本背景和电磁计算的研究,然后介绍了辛算法。
接着,介绍了辛算法在Maxwell方程中的应用,然后在无耗煤质和散射存在时的情况下分析了辛时域有限差分法的计算式。
最后,以真空中一维的高斯脉冲电磁波为例用辛算法进行了数值运算。
关键词:电磁计算;辛算法;Hamilton系统;Maxwell方程一.引言电磁场理论的应用遍及地理学、生命科学、医学、材料科学和信息科学等几乎所有技术学科领域。
计算电磁学是以电磁场理论为基础,以高性能的计算技术为手段,运用计算数学提供的各种方法,解决复杂电磁场理论和工程问题的应用科学。
因此,开展计算电磁学的研究不仅可以产生国际水平的基础研究成果,更重要的是可以促进我国民用和军用电磁学相关领域的发展。
早在1864年,Maxwell在前人理论和实验的基础上建立了统一的电磁场理论,并用数学模型揭示了自然界一切宏观电磁现象所遵循的普遍规律,这就是Maxwell方程组,它包括微分形式和积分形式。
简单地说,所有的宏观电磁问题都可以归结为Maxwell方程组在各种边界条件下的求解问题。
计算电磁学自20世纪60年代兴起,至今40余年。
纵观整个电磁理论发展的过程,电磁学的发展可以分为两个阶段。
以20世纪60年代为分界点,之前可以称为经典电磁学阶段,在这个时期,电磁场理论和工程中的许多问题大多采用解析或渐进的方法进行处理,即在几种可分离变量的坐标系中求解Maxwell方程组或其退化形式,最后得到解析解。
哈密尔顿系统有限元的守恒性和辛性质哈密尔顿系统是最重要的动力系统。
冯康院士曾指出,一切真实的无耗散的物理过程都可表示为这样或那样的Hamilton形式,它们都是常微分或偏微分方程组。
Hamilton系统有两个最重要的特性:守恒性和辛结构。
在数值计算中能否保持这些特性具有重要意义。
1983年冯康研究此问题,他惊讶地发现,此前这里竟是一片空白,许多经典的算法都不适应,计算几万步后有时已面目全非。
他1984年首创性提出辛差分算法,并作了深入系统的研究,开辟了一大片研究新领域。
以后国内外许多学者作了多方面推广和广泛应用。
冯康的这项首创工作得到了国际一致公认。
但是任何离散算法,一般不可能同时保辛又保能量(Ge-Masden定理)。
辛差分算法很好地保辛,但只在格式精度意义下保能量。
而许多学者认为,有时保能量更重要。
因此我们转向有限元法,却发现至今有关研究极少。
而我们的研究表明,有限元总是保能量的,对线性系统也是辛的,对非线性系统是高精度保辛的,而且长时间计算稳定且精度高,效果非常好。
这些结论已刻划了有限元的基本特征。
因此有限元法是与辛差分算法完全不同的另一种算法,从另一方面弥补了辛差分算法的不足。
本研究是对辛算法的一次重要推进。
本文主要创新点如下:(1).首次系统深入研究任意m次有限元解非线性Hamilton系统,证明了在任何节点上能量总是守恒的,因此在相平面上轨道总是稳定的。
并首次提出用超收敛分析方法研究有限元的辛性质;(2).对线性Hamilton系统的任意m次有限元,得到了一个深刻的高阶超收敛O(h<sup>2m+1</sup>)新估计,首次证明m次有限元的节点值是2m阶对角Páde 逼近,因而是辛格式。
此结果与冯康等研究的辛差分格式结论一致;(3).对非线性Hamilton系统的任意m次有限元法,构造了新的辅助问题,并得到误差估计和负范数估计,首次证明m次有限元对每一次步进是高精度O(h<sup>2m+1</sup>)意义下近似保辛的。
哈密尔顿偏微分方程的多辛算法
本文主要研究了哈密尔顿(Hamilton)偏微分方程的多辛算法。
文中首先介绍了辛算法在求解Hamilton系统中的重要性,随后详细阐叙了Hamilton系统的发展历史和辛算法的发展历史与现状,并简单的介绍了本文所进行的工作。
其次,介绍了Hamilton系统的一些基本概念与保持Hamilton系统辛结构的辛算法即辛Runge-Kutta(RK)方法以及相关的辛方法,并给出了稳定性分析常用到的方法—变量分离法以及一些判断稳定性的常用定理。
而后,在前面介绍的知识基础上对膜自由振动方程应用多辛的Runge-Kutta-Nystr?m(RKN)方法。
首先提出了膜自由振动方程的一个多辛形式,进而构造多辛的RKN格式并证明了该格式满足离散的多辛守恒律;为了便于数值实验随后构造了一个显式辛格式并且给出这个显式格式是稳定的一个充分条件,通过数值实验说明多辛的RKN 方法不仅对解有长时间的模拟而且能够保持一些重要的物理守恒量。
最后,讨论了非线性Boussinesq方程忽略它的非线性部分的多辛RKN方法。
先给出了方程的一个多辛形式,接着构造了多辛的RKN格式与相应的离散多辛守恒律,为了数值模拟给出了一个显式的辛RKN格式进而给出格式为稳定的一个充分条件,最后用数值实验说明了辛算法离散方程的优越性。
mindlin板动力学问题的hamilton体系及其辛解法Mindlin板动力学问题的Hamilton体系及其辛解法一、简介Mindlin板动力学问题是由Richard Mindlin在1943年提出的,它是一个复杂的动力学问题,是一种多物理量耦合的动力学系统,主要涉及到板的弯曲和剪切,以及应力、应变和位移。
Mindlin板动力学问题中,Hamilton体系是一种可以描述Mindlin板动力学问题的数学模型,它通过求解相关的动力学方程将板的应力,应变和位移联系起来,从而帮助我们更好的理解板的动力学行为。
二、Mindlin板动力学问题的Hamilton体系1、定义 Mindlin板动力学问题的Hamilton体系是一种用于描述Mindlin板动力学问题的数学模型,它由一系列相关的动力学方程组组成。
其中,包括有:板上每一点的位移方程、板上每一点的力对位移的响应方程、板上每一点的应力和应变之间的关系方程以及板上每一点的力和力矩之间的关系方程等。
2、位移方程位移方程是描述板动力学问题的基本方程,它表示板上每一点的位移u,在时间t上的变化状态,即∂u/∂t=v,其中v是板上每一点的速度,由于板的位移受到板的自重和外界力的影响,因此可以得到:∂v/∂t=-g-f (1)其中,g表示板的自重,f表示外界力。
3、力对位移的响应方程力对位移的响应方程描述的是板上每一点的力对位移的反作用,即力f对位移u的响应,从而使位移随时间发生变化。
对于Mindlin板,由于板的弯曲和剪切都会对位移产生影响,因此可以得到:f=Kuu+Kuv (2)其中,Kuu表示板的弯曲,Kuv表示板的剪切,Kuu和Kuv分别是板的弯曲和剪切系数。
4、应力和应变之间的关系方程应力和应变之间的关系方程描述的是板上每一点的应力σ和应变ε之间的关系,即σ=Eε,其中E是板的杨氏模量。
5、力和力矩之间的关系方程力和力矩之间的关系方程描述的是板上每一点的力f和力矩τ之间的关系,即τ=Gf,其中G是板的刚度矩阵。
薄板问题的Hamilton体系和辛几何方法
邹贵平
【期刊名称】《应用基础与工程科学学报》
【年(卷),期】1996(0)4
【摘要】本文通过薄板问题混合能变分原理,选用状态变量及其对偶变量,导出了一般的Hamilton型广义变分原理和Hamilton正则方程,这样就突破了欧几里德空间的限制,在Hamilton力学的数学框架辛几何空间中,对全状态相变量进行分离变量,并采用共轭辛正交归一关系,给出任意支承条件下薄板问题的辛精确解.
【总页数】9页(P335-343)
【关键词】Hamilton体系;辛几何;正则方程;弹性薄板
【作者】邹贵平
【作者单位】同济大学工程力学与技术科学系
【正文语种】中文
【中图分类】O343
【相关文献】
1.复合材料力学的Hamilton体系和辛几何方法(Ⅰ)——一般原理 [J], 钟万勰;欧阳华江
2.复合材料力学的Hamilton体系和辛几何方法(Ⅱ)——平面问题 [J], 钟万勰;欧阳华江
3.反对称铺设层合板动力问题的Hamilton体系及辛几何解法 [J], 邹贵平
4.复合材料力学的Hamilton体系和辛几何方法(Ⅲ)——弯曲问题和板的振动 [J], 欧阳华江;钟万勰
5.考虑剪切效应层合板的Hamilton体系及辛几何方法 [J], 邹贵平
因版权原因,仅展示原文概要,查看原文内容请购买。
哈密尔顿矩阵特征谱问题的辛算法
丁克伟
【期刊名称】《合肥工业大学学报(自然科学版)》
【年(卷),期】2000(023)003
【摘要】文章基于前人的工作,在哈密尔顿矩阵约化过程中,采用了辛相似变换,使得哈密尔顿矩阵在辛相似变换下仍保持Hamilton结构,这样从根本上确保了特征值的正确性和稳定性,也能保证特征值成对出现且在每个半平面上都只求得n个特征值,不至于出现特征值在小扰动下跨过虚轴的混乱局面.
【总页数】5页(P336-340)
【作者】丁克伟
【作者单位】大连理工大学,工程力学系,辽宁,大连,116023
【正文语种】中文
【中图分类】O151.21
【相关文献】
1.测试Hamiltonian矩阵结构问题的辛算法 [J], 丁克伟
panion矩阵的伪谱问题 [J], 张敏王正盛徐贵力;
3.Hamiltonian矩阵平方约化求解特征问题的辛算法 [J], 丁克伟
4.哈密尔顿偏微分方程多辛算法 [J], 王雨顺;洪佳林
5.Hamiltonian矩阵特征谱问题的辛算法 [J], 丁克伟
因版权原因,仅展示原文概要,查看原文内容请购买。
《辛对称Hamilton算子的可逆性与Fredholm性》篇一一、引言在数学物理和偏微分方程的研究中,Hamilton算子扮演着重要的角色。
其不仅在经典力学中有着广泛的应用,也在量子力学、光学和电磁学等领域中发挥着关键作用。
辛对称Hamilton算子作为Hamilton算子的一种特殊形式,其可逆性和Fredholm性更是研究的热点。
本文旨在探讨辛对称Hamilton算子的可逆性与Fredholm性,为相关领域的研究提供理论支持。
二、辛对称Hamilton算子的定义与性质辛对称Hamilton算子是一种特殊的偏微分算子,具有辛对称性质。
其定义涉及复数域上的函数空间及相应的偏导数运算。
该算子具有自伴性、正定性等基本性质,这些性质使得辛对称Hamilton算子在物理和工程领域具有广泛的应用。
三、可逆性的研究可逆性是算子理论中的一个基本概念,对于辛对称Hamilton 算子而言,其可逆性与其定义域、值域及算子的其他性质密切相关。
本文将从定义域和值域的角度出发,探讨辛对称Hamilton算子的可逆性条件。
通过严密的数学推导,得出辛对称Hamilton算子可逆的充分必要条件。
四、Fredholm性的研究Fredholm性是算子理论中的另一个重要概念,与算子的谱结构、本征值及本征函数等密切相关。
本文将研究辛对称Hamilton 算子的Fredholm性质,包括其Fredholm指标的计算及性质。
通过分析算子的谱结构,得出辛对称Hamilton算子为Fredholm算子的条件。
五、数值分析与实例验证为了验证理论的正确性,本文将通过数值分析的方法,对辛对称Hamilton算子的可逆性和Fredholm性进行实例验证。
通过对比理论计算结果与实际数值结果,验证了本文所提出理论的正确性和有效性。
六、结论本文通过对辛对称Hamilton算子的可逆性和Fredholm性的研究,得出以下结论:1. 辛对称Hamilton算子的可逆性与其定义域、值域及算子的其他性质密切相关。