二维自适应结构网格的变分生成方法

提出-一种二维结构网格的自适应生成方法,即采用变分方程控制网格的疏密程度以及光滑性、正交性和正规性,并且所有的控制方程具有-致的量纲.几个数值算例表明,给出的网格自适应生成方法能够较好地收敛到方程的唯-解,同时产生性质良好的二维自适应结构网格.     

弹塑性流体力学拉氏自适应网格和网格重分守恒重映数值模拟方法

拉氏自适应重分弹塑性流体力学有限元程序实现了网格完全自适应，具有良好、灵活的非结构自适应网格数据结构，实现了滑移界面两边(接触间断)网格动态调整，网格的细分和合并处理灵活，网格重分和网格自适应模块兼容、守恒重映，网格重分中采用多种方法控制新网格的质量，爆轰计算可采用Lee-Tarver的化学反应率模式。初步数值计算结果表明，弹塑性流体力学拉氏自适应重分数值模拟方法合理，计算结果正确，基本反映了流场的物理结构。     

复杂区域上的一种结构网格生成方法

讨论复杂区域上的一种结构网格生成方法,其主要思想是:以变分形式的Winslow网格生成方法为基础,通过引入网格解扭机制和网格面积均匀化技术,构造出一种新的离散泛函,进而采用一类优化算法求解这一离散泛函的极小化问题,得到所希望的网格.通过分析及大量数值实验表明,这一方法比较健壮,针对二维复杂区域通常能够生成几何品质较优的网格,它在保持Winslow方法优点的同时,克服了它的一些缺点.     

非结构多边形网格解扭和重分方法

在多维流体动力学计算中，流体运动和计算网格的关系可以分为两种情况。一是Lagrangian方法，即网格跟随流体运动；二是Eulerian方法，即流体流过固定；下动的网格。一般计算网格的运动是任意的。这就对应于任意Lagrangian—Eulerian（ALE）方法。ALE方法的核心是通过调整网格运动，使得数值模拟的精度、效率有所提高。它的主要步骤是：显式Lagrangian步；网格重分，即得到新的计算网格；物理量重映，即将Lagrangian步的计算结果变换到新网格上。在这3步中，较少研究网格重分。数值模拟和网格重分的一个基本前提是网格是合理的，或者说网格不能发生翻转，网格应当是凸的。而Lagrangian步数值模拟会造成网格扭曲，因此在网格重分前进行网格解扭是十分必要的。文中描述了通用的网格解扭、重分算法，使得解扭、重分后的网格有较好的几何品质，同时尽可能接近Lagrangian网格。     

推进方式的参考Jacobian网格优化方法

将参考Jacobian方法与前沿推进方法相结合,设计一种推进方式的参考Jacobian网格优化方法.该方法从计算区域的部分边界开始,沿着给定的方向逐步推进,每前进一步都以相邻的两行(或两列)结点为边界采用参考Jacobian方法进行优化,其中处于后方的一行(或一列)使用的是前一步优化后的结点.分析及数值实验表明,与参考Jacobian方法相比,不仅能大幅度降低执行时间,而且优化后的网格在几何品质方面相当接近,甚至更优,得到的网格与Lagrange网格更接近.     

网格重分中关于守恒重映的几个问题

首先针对离散情况,对垸恒重映中的“守恒”给出数学上严格的解释。在此基础上,对守恒重映问题中所固有的守恒量协调问题作了描述,并提出一种基于最小二乘法的处理方案,同时,还分析了容易混淆的“质量密度”和“体积密度”,网格中心量和网格节点量等基本概念,并采用两种不同的辅助网格来对网格节点量进行重映。     

三维非结构网格自动生成

采用各向异性平面非结构网格生成技术对曲面的参数平面进行三角剖分,从而得到曲面的非结构网格,作为三维非结构网格的边界网格.应用推进面法生成网格内点,增量法将生成的内点逐点插入现有网格进行网格细化,得到三维计算域的Delaunay非结构网格.讨论了非结构网格质量优化方法.给出几个算例说明方法的应用.     

大长宽比扭曲网格上辅助网格差分方法

讨论抛物型方程的离散差分格式的构造,对九点差分格式进行了适用范围的讨论,并在此基础上提出辅助网格差分方法,用于处理因网格长宽比大且扭曲较大的网格引起的计算精度与计算效率降低的问题,该方法从守恒方程出发,将九点差分格式应用于按某种合适的方式进行重分之后的网格上,减少由于网格正则性差以及网格节点上的物理量采用周围网格量的加权平均等原因所引起的计算误差,得到一个新的但其解仍然逼近原来网格上的物理量的方程组.所构造的方法便于实施,且更适合于对实际物理模型的模拟,能比较好地适应流体大变形导致的网格扭曲,数值试验表明它有较好的数值精度和稳定性.     

内部带特征约束的任意平面域的三角形网格生成方法

基于Ｄｅｌａｕｎａｙ三角化方法 ,给出了内部带有任意复杂特征约束的任意平面域的三角形网格生成方法。区域内部特征约束可以是点、直线、圆、圆弧、平面样条曲线以及由它们组合构成的封闭环或者不封闭环。     

理想不可压缩MHD的广义变分原理

本文应用钱伟长方法,在前工作的基础上,在一定的条件下,探讨了理想不可压缩的磁流体力学(MHD)的广义变分原理。     

热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.     

基于变分方法的混沌系统参数估计

提出一种基于变分原理的估计混沌系统未知参数的方法,对以x= F(x,θ) 为控制方程的所有混沌系统具有普适性.首先将混沌系统方程引入到目标泛函中;接着利用变分原理导出了混沌系统的伴随方程和待辨识参数泛函梯度的通用公式;然后设计了估计混沌系统未知参数的算法;最后对典型的Lorenz混沌系统和超混沌Chen系统的未知参数进行了估计.数值仿真结果表明该方法是一种非常有效的估计混沌系统未知参数的方法. 关键词： 混沌系统 参数估计 变分方法 伴随方程     

用分层求解的方法快速形成二维数值网格

讨论了用分层求解的思想来快速构造二维数值网格。把以往文献中针对于外场区域的解抛物型方程的方法推广到封闭的二维区域,并把分层求解的思想应用于更为复杂的用椭圆型方程构造网格的数值求解中,给出了一些闭域的网格构造图。     

非线性映射参数辨识的离散变分方法

提出一种辨识非线性映射系统中未知参数的离散变分方法, 对以x_k+1 = F(x_k,θ)为状态控制方程的所有映射混沌系统具有通用性. 对典型的Logistic映射和Henón映射中的未知参数进行了估计, 仿真结果表明了该方法的有效性.     

基于变分原理的三维扩散方程离散化格式

研究三维扩散方程的数值模拟.在非正规六面体网格上,用变分原理建立扩散流形式的差分格式.将扩散流通量作为未知函数,求其泛函极值,并与原函数温度联立求解.详细推导了计算公式,并给出了数值试验结果.     

贴体坐标网格生成技术的研究

本文对数值求解微分方程法生成贴体坐标网格的技术进行了探讨和尝试。文中不仅对数值求解双调和方程的两种不同方法作了介绍,还利用该方程生成了典型的Ｈ型、Ｃ型、Ｏ型网格。结果表明,双调和方程法在控制网格线间距及与边界的正交性方面具有很强的能力,所生成的网格质量令人满意,且可用于复杂的物理区域。因而,采用双调和方程生成贴体坐标网格是流场数值模拟中一种有效的、高质量、具有广泛通用性的网格生成技术。     

广义经典力学系统的第一积分与变分方程特解的联系

给出了广义经典力学系统的变分方程,研究了变分方程的解与系统的第一积分之间的联系,并证明可由系统的第一积分求得变分方程的特解. 关键词： 广义经典力学系统 第一积分 变分方程     

基于二阶离散变分方法的非线性映射参数估计

提出一种估计非线性映射未知参数的二阶离散变分方法.首先针对非线性离散混沌系统, 利用变分方法导出了伴随方程和目标泛函梯度, 以此为基础利用二阶离散变分方法给出了二阶伴随方程和精确计算Hessian矩阵-向量乘积的显式表达式; 其次设计了估计非线性映射未知参数的新算法, 并以此对Hyperhenón映射和二维抛物映射中的未知参数进行了精确的估计. 数值仿真结果表明了该方法的有效性和优点. 关键词： 非线性映射 参数估计 二阶离散变分方法 伴随方程     

强磁场中氢原子能级的另一种绝热变分计算

采用绝热近似和变分原理，引入描述Coulomb作用的等效模型势，给出绝热变分方程，计算了均匀强磁场(1≤β≤2000)中氢原子的基态和几个较低激发态的能级，并与绝热近似方法和公认为最精确的类多组态Hartree-Fock方法的结果进行了比较，给出的基矢比Landau基矢具有更好的收敛性. 关键词：