三维双曲型守恒律非结构网格最小二乘有限体积方法研究

来源 :国防科技大学 | 被引量 : 0次 | 上传用户:a1a1a1a1a1a1a1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
双曲型守恒律方程组在流体力学、空气动力学、航空航天和生物学等领域都有着广泛的应用,一般情况下其解析解不存在,这就推动了此类方程数值方法的深入研究与发展,例如有限差分方法、有限元方法、有限体积方法等。有限体积方法离散方式灵活多样,适用于对复杂区域进行复杂问题的计算,目前已成功运用于计算流体力学、气象学、计算电磁学、半导体计算以及计算生物学等工程领域。考虑到双曲型守恒律方程组的解中包含激波、接触间断等特征,即在局部区域物理量发生突变,因此网格和算法的局部自适应技术在此类方程中的应用也得到了广泛关注。鉴于以上讨论,本文针对三维Euler方程,在非结构四面体网格及移动网格上发展了一类非振荡有限体积方法。本文主要工作如下:1、对于三维双曲型守恒律,介绍了有限体积方法的基本原理、网格和控制体、TVD型Runge-Kutta时间离散方法、数值流通量等内容。2、构造了一类在非结构四面体网格上求解三维Euler方程的非振荡有限体积方法。采用时空分离方法,空间上以每个四面体单元作为控制体,采用Roe通量(或HLLC通量)进行离散,同时为提高数值精度,运用最小二乘方法通过求解超定线性方程组来重构每个四面体网格上物理量的线性插值多项式,并添加限制器以使其满足半离散形式的局部极值原理;时间上采用两步Runge-Kutta TVD方法进行离散。在保证高分辨率的前提下,在空间计算时该方法比照ENO格式避免了模板的选取,极大减少了计算量。3、讨论了一种基于弹簧振子的三维移动网格方法,并将上述非振荡有限体积方法成功应用于该移动网格上,同时找到了一种守恒的加权重构方法对每个时间步网格上的守恒物理量进行更新。
其他文献
我国复杂的地质条件致使地质灾害频发,地质灾害监测技术及监测装置的研究为防灾减灾提供理论和数据支持。岩土体发生滑坡、坍塌等地质灾害时其内部必定会发生相对位移产生裂缝。岩土体深部位移监测作为地质灾害监测最为重要的一部分,为地质灾害监测及预警预报提供准确、有效的数据支持。针对目前现有的岩土体深部位移监测技术的不足,如适用场合、测量范围、安装困难以及耐用性等,本文基于互感效应、霍尔效应和重力测斜等理论原理
四维变分同化技术是提高数值天气预报时效的重要方法,切线性模式和伴随模式是同化系统的关键技术。本文针对我国自主研发的多尺度通用数值预报模式系统GRAPES,为建立GRAPES全球四维变分同化系统的切线和伴随模式,开展了以下工作:基于GRAPES全球模式绝热动力学框架和四维伴随变分资料同化理论,建立了GRAPES全球模式的切线性模式(GRAPES TL)和伴随模式(GRAPES AD),切线性模式精度
石灰行业是一个古老的行业,在中国已经有几千年的历史。远在周朝时期,古人就用蛤的外壳烧制石灰,蛤壳的主要成分为碳酸钙,煅烧时,二氧化碳全部逸出后留下的蜃灰就是石灰。在崇尚厚葬的古代,将蜃灰作为胶凝材料来修建陵墓非常普遍。而中国古代的巨型工程——万里长城,它的许多地段被后人发现是用石灰砌筑而成的。明代的《天工开物》和清代的《营造法原》就详细地记载了
期刊
分子的高次谐波是分子与强激光场相互作用产生的非线性效应之一,氢分子离子作为最简单的分子成为这一课题中研究的热点对象。本文通过数值求解含时Schr?o?d inger方程对氢分子离子的高次谐波谱进行了数值模拟,研究了氢分子离子的高次谐波产生过程中的干涉效应,高次谐波的辐射强度对原子核间距的依赖关系及利用高次谐波谱探测氢分子离子的振动频率。在Born-Oppenheimer近似下,对双原子分子高次谐波
由多种轨道类型构成的混合星座卫星导航系统能以较少的卫星数目实现对区域的良好覆盖,导航星座中包含多种轨道类型卫星,但目前卫星导航算法多针对等高近圆单一轨道类型的Walker星座,在混合星座中并不适用。本文基于混合星座卫星导航系统,研究在差分定位、与全球卫星导航系统组合导航等方面的应用问题。混合星座中不同轨道类型卫星的观测值精度不一致,本文利用用户测距误差进行加权单点定位可以获得更高的精度。同样地,伪
在高功率微波和脉冲功率技术工业应用的推动下,重频强流脉冲加速器逐步成为脉冲功率技术研究的热点和发展趋势之一。初级能源回路的功率供给与能源补充是实现加速器重频运行的一项关键技术。论文研究了利用中间储能与直接供电两种方式实现重频强流脉冲加速器的初级能源系统。首先对利用储能电容实现加速器重频运行的初级能源进行了理论、模拟和实验研究;然后对利用直接供电方式实现加速器重频运行的初级能源进行了电路分析;最后对
在现在已经发现的六味夸克中,当取重夸克Q = c ,b质量极限( mQ→∞)时,介子具有明显的自旋-味对称性,这个时候的有效理论称为重夸克有效理论(HQET) ;而当取轻夸克( q = u , d ,s)质量极限( mq→0)时,QCD拉格朗日量则会表现出手征对称性,这个时候的有效理论称为手征微扰理论(ChPT)。在QCD的近似对称下我们可以同时采取这两种质量极限建立一个有效理论——重手征理论,用
本文在理论分析的基础上,对炸药水下爆炸气泡脉动进行了数值模拟研究,模拟了水下爆炸冲击波的传播及气泡脉动过程。应用AUTODYN有限元程序建立一维计算模型,对1 kg球形TNT装药水下10米爆炸气泡脉动进行了数值模拟,冲击波峰值压力、气泡第一次脉动周期的数值模拟值均与经验公式符合的比较好。应用AUTODYN有限元程序建立二维轴对称计算模型,对12 g球形TNT装药水下10米爆炸气泡脉动进行了数值模拟
在卫星导航定位系统中,导航卫星向用户发播的导航信息包括测距码和数据码,数据码包括卫星广播星历参数、卫星钟差参数和电离层延迟改正参数等,这些参数是由主控站经数据处理后注入到导航卫星的,它们对用户的导航定位精度起着至关重要的作用。因此,本文对导航定位系统中与主控站相关的一些数据处理问题进行了研究。在距离观测量的数据处理中,针对目前广泛使用的Hatch滤波平滑伪距方法,本文进行了仿真分析,并与其它伪距平
2014铝合金是航天材料中应用最广泛的一种铝合金。本文主要采用实验的方法对2014铝合金的层裂特性进行研究,侧重于考虑不同的热处理过程对材料层裂强度的影响,希望对材料的制造工艺和使用过程提供有益的指导。论文通过实验测量了材料的层裂强度,并对实验回收试样进行了显微观察,对层裂的损伤演化过程进行了分析,结合铝合金热处理理论对铝合金材料层裂特性影响因素进行了初步讨论。通过上述工作,本文取得的主要成果有: