一维守恒的Lagrangian ADER格式 A 1D Conservative Lagrangian ADER Scheme.pdfVIP

第30卷第4期 计 算 物 理 V01．30．NO．4 CHINESE OFCOMPUTATIONALPHYSICS 2013年7月 JOURNAL July，2013 文章编号：1001—246X(2013104-0501-08 ADER格式 一维守恒的Lagrangian 程军波1， EleuterioF．Toro of Trento Mathematics (1．北京应用物理与计算数学研究所，北京100094；2．LaboratoryApplied University，Trento，haly) 摘 要：基于欧拉框架下ADER格式，构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格 式．构造r阶LADER格式包括：从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重 构节点处守恒量和从1阶到r—l阶的空间导数、求拉氏框架下这些变量的Godunov值，并计算1阶到r一1阶的时间 全导数，最后高精度离散积分形式的流通量函数．对光滑流场的模拟表明，LADER格式达到设计的精度；对含强间 断的流场模拟表明，数值解在间断附近基本无振荡． 关键词：中心型拉格朗日格式；任意高精度格式；时空高精度 中图分类号：0351 文献标识码：A O 引言 2001年，Toro等[1。引入具有任意精度的ADER(arbitrarilyhigh—order 有同等的高精度． 现有ADER格式建立在欧拉框架下，网格没有移动．我们希望将其推广到拉格朗日框架．与欧拉框架相 比，拉格朗日框架下网格随流体一起运动．根据物理量所定义位置的不同，拉格朗日方法大致分为两类：①为 运动量(如位移、速度和加速度)定义在网格节点、状态量(如密度、压力、内能等)定义在网格中心的交错型 拉格朗日方法口3；②为所有物理量(包括运动量和状态量)均定义在网格中心的中心型拉格朗日方法(有时 也称为非交错型拉格朗日方法)．中心型拉格朗日方法将所有物理量都定义在单元中心，由此可从守恒方程 出发，借用守恒型格式的一些理论结果，如采用Godunov格式来计算所有物理量．关于中心型拉格朗日方法 有大量文献Hq】．从这些文献中我们发现，现有中心型拉格朗日方法最高精度只有三阶．这远远比不上欧拉 型高精度Godunov方法． 我们希望提高中心型拉格朗日方法的离散精度。为此，我们借鉴欧拉框架下ADER格式的思想，构造拉 格朗日框架下的ADER格式，从而使格式在时空两个方向上都具有较高精度．本课题研究高精度的中心型 拉格朗日格式，最终建立守恒、基本不振荡的高精度中心型拉格朗日格式． 1 一维拉格朗日ADER格式 1．1基本推导过程 考虑可压缩欧拉方程 堕+曼!!塑：o． at d戈 其中， U=(p，pu，pE)7， (2) 收稿日期：2012—09—11；修回日期：2013—03—13 基金项目：国家自然科学基金及中国工程物理研究院科学发展基金重点基金(2009A0202014)资助项目 万方数据 502 计 算 物 理 第30卷 F=(pu，pu2+P，(pE+P)u)T． (3) 其中，p、p,u．E分别表示压力、密度、速度和单位质量总能量，E