带分数阶耗散项的Quasi-geostrophic方程周期解的大时间性态

本文主要讨论带有分数阶耗散项的quasi-geostrophic(QG)方程的周期解的大时间性态,并且对初值θ_0没有小性的要求.对于次临界(1/2α1)和临界(α=1/2)两种情况,证明了方程的周期解均具有指数级的衰减性.     

对流扩散方程的缩减基有限元法的后验误差

将缩减基(RB)方法和有限元方法相结合,在保证偏微分方程的有限元离散格式具有足够高精确度前提下,能够大幅度地降低有限元离散格式的维数,从而大大降低计算中内存容量和计算时间的消耗.针对对流扩散方程建立基于RB方法的Crank-Nicolson有限元离散格式,并给出后验误差估计结果.     

悬链线竖向集中力作用下的构形分析和张力计算

目前,悬链线在竖向集中力和均布荷载共同作用下的构形分析和受力计算的理论仍不完善。针对这一问题,通过引入悬链线的几何约束方程、力平衡方程和超越方程,建立了竖向集中力与均布荷载共同作用下的非线性方程组。采用牛顿迭代法求解方程组,得到了悬链线的构形和受力情况。为了验证理论计算的正确性,进行了算例和试验验证。结果表明,算例的计算结果与文献结论保持一致,试验测得的构形和水平张力大小与理论计算的构形和水平张力大小吻合较好。本文的理论计算可以更加简单精确地计算出悬链线在竖向集中力和均布荷载共同作用下的构形和受力情况,为实际工程提供重要的理论指导。     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

改进的PNC方法及其在非线性偏微分方程多解问题中的应用

2017年，李昭祥等提出了一种偏牛顿-校正法（Partial Newton-Correction Method，简记为PNC方法），并利用它成功地计算出了三类非线性偏微分方程的多重不稳定解.本文在PNC方法的基础上，提出并发展了一种改进的PNC方法.首先，利用Nehari流形$\mathcal{N}$与零平凡解的可分离性，建立并证明了$\mathcal{N}$的某特殊子流形$\mathcal{M}$上的全局分离定理及其推广（即局部分离定理）.全局分离定理只跟非线性偏微分算子或相应的非线性泛函本身有关，而与具体的计算方法无关.对一些典型的非线性偏微分方程多解问题（比如，Henon方程问题），该全局分离定理的分离条件，经验证是成立的.另一个方面，通过修改或补充原辅助变换的定义，去掉了原辅助变换的奇异性；接着建立并证明了某些非线性偏微分方程问题的新未知解与该非线性偏微分算子零核空间的密切关系；在证明中，去掉了在原奇异变换下所需的标准收敛（standard convergence）假设.最后，计算实例与数值结果验证了改进的PNC方法的可行性和有效性；同时表明子流形$\mathcal{M}$与已知解的可分离性是PNC方法和本文新方法能成功找到多解的关键.     

乏燃料贮存格架跌落事故冲击分析和试验验证

乏燃料贮存格架是核电厂贮存乏燃料的重要设备，在满载条件下和地震/跌落事故中，都应保持稳定和安全状态．本文基于LS-DYNA对乏燃料贮存格架进行了跌落事故冲击分析，考虑了最重重物从可能最高处意外跌落的情况．分析时考虑了碰撞、几何大变形、材料非线性等非线性因素．分析发现，浅跌落情况下贮存格架变形较大，但为局部变形，冲击载荷不影响贮存格架的安全功能．深跌落情况下，组件跌落在支座上方时支座承受的载荷最大，并在许用载荷范围内．为确保核电厂安全性和可靠性，基于分析得到的现象，设计和开展了乏燃料贮存格架浅跌落和深跌落试验．试验采用了等比例贮存腔和真实燃料组件的管座．同时对试验件进行了跌落分析，并与试验结果进行了对比，验证了分析技术的保守性和准确性．掌握的分析技术可应用于所有压水堆核电站乏燃料贮存格架的跌落事故分析．     

参外联合激励下的簇发振荡及其机理

当周期激励频率远小于系统固有频率时，会存在快慢耦合效应，与单项激励不同，参外联合激励不仅会导致快子系统平衡曲线和分岔行为的复杂化，也会产生一些特殊的非线性现象，为此，本文以两耦合Hodgkin-Huxley细胞模型为例，引入周期参外联合激励，探讨在频域不同尺度耦合时该系统的簇发振荡的特点及其分岔机制．通过建立相应的快慢子系统，得到慢变参数变化下的快子系统的各种分岔模式以及相应的分岔行为，结合转换相图，揭示耦合系统随激励幅值变化时的动力学行为及其机理．研究表明，在激励幅值较小时，系统表现为概周期振荡，两频率分别近似于快子系统平衡曲线由Hopf分岔引起的两稳定极限环的振荡频率．概周期解随激励幅值的增加进入簇发振荡，导致这些簇发振荡的主要原因是在慢变参数变化的部分区间内，存在唯一稳定的平衡曲线，使得系统的轨迹逐渐趋向该平衡曲线，产生沉寂态，并随着慢变参数的变化，由分岔进入激发态．同时，快子系统中参与簇发振荡的稳定吸引子随激励幅值的变化也会不同，导致不同形式的簇发振荡．另外，与单项激励下的情形不同，联合激励时快子系统的部分稳定吸引子掩埋在其它稳定吸引子内，从而失去对簇发振荡的影响．     

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

正负双向脉冲式爆炸及其诱导的簇发振荡

簇发振荡普遍存在.探索通向簇发振荡的可能路径是簇发研究的热点问题之一."脉冲式爆炸(pulsed-shaped explosion,PSE)"是一种最近被报道的可以诱发簇发振荡的新机制,其特征为平衡点和极限环表现出了与参数变化相关的脉冲式急剧量变.PSE会导致系统轨线急剧跃迁,从而诱发典型的簇发振荡.然而,目前报道的PSE中仅含有"单向的尖峰",未发现"双向的尖峰",且由其诱发的簇发振荡仅含单向的振荡簇.本文以多频激励Rayleigh系统为例,旨在揭示PSE的不同表现形式以及与此相关的簇发动力学.利用频率转换快慢分析法得到了Rayleigh系统的快子系统和慢变量.针对快子系统的分析表明,PSE表现出了较为复杂的动力学特性,其特征是PSE包含了正负双向两个不同的尖峰,此即所谓的正负双向PSE.其急剧量变行为,导致了系统轨线在单个振荡周期内出现正向和负向的多次跃迁,由此得到了由正负双向PSE所诱发的簇发振荡.根据吸引子类型分别揭示了点--点型和环--环型两类簇发振荡模式的产生机制.本文的研究给出了PSE的不同表现形式,丰富了多时间尺度下的簇发振荡的诱发机制.      

基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.