共查询到20条相似文献,搜索用时 788 毫秒
1.
提出基于细分和数值积分思想的一种离散的守恒重映方法——质点重映方法.密度分布可采用一阶精度的分片常数分布,或二阶精度的分片线性分布.分片线性密度分布函数采用面平均方法构造.重映过程中,借助四边形辅助网格,实现了交错网格节点量的重映.质点重映方法既适用于结构网格,也适用于非结构网格,且不要求新旧网格之间一一对应.数值结果表明,一阶精度重映算法健壮性好,但会产生较大的扩散效应;二阶精度重映算法可较好地保持密度分布的特性,但存在单调性问题.为改善二阶精度重映方法单调性,将结构网格质量守恒调整算法推广到非结构网格上,以限制新网格的质量密度.给出了一些重映的例子,并进行了误差分析. 相似文献
2.
王瑞利 《工程物理研究院科技年报》2003,(1):427-427
在计算流体力学(CFD)领域,几乎所有的方法都离不开网格,网格是各种数值方法的基础。网格质量的好坏直接影响数值结果的精度,甚至影响到数值计算的成败。为此CFD工作者发展了许多方法。如迭合网格、贴体网格和非结构网格,为了更好地数值模拟大变形问题,又进一步发展了结构/非结构混合网格的技术,尤其是发展了网格跟随流场智能化调整的网格自适应技术。这些网格技术的发展,几乎都涉及网格的变动,只要改动网格就涉及物理量的重映。重映方法一般被分为两种类型:插植重映和积分重映。所谓插植重映方法就是在计算区域D上,用已知网格上的物理量分布,通过插值理论把它插值到新网格或任意定义的规则网格上的一个过程,通过这个过程给出新网格的物理量的分布。所谓积分重映方法就是用积分的形式把旧网格上的守恒量重新映射到新网格上。如对某种体积密度分布q,简单的积分形式为 相似文献
3.
任意网格重映的样条逼近算法 总被引:2,自引:1,他引:1
在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果. 相似文献
4.
5.
6.
拉氏自适应重分弹塑性流体力学有限元程序实现了网格完全自适应,具有良好、灵活的非结构自适应网格数据结构,实现了滑移界面两边(接触间断)网格动态调整,网格的细分和合并处理灵活,网格重分和网格自适应模块兼容、守恒重映,网格重分中采用多种方法控制新网格的质量,爆轰计算可采用Lee-Tarver的化学反应率模式。初步数值计算结果表明,弹塑性流体力学拉氏自适应重分数值模拟方法合理,计算结果正确,基本反映了流场的物理结构。 相似文献
7.
介绍了二维非结构网格上的守恒重映算法,重点是基于SFB/DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法,包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法,不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映,介绍了SALE和HIS算法。此外,为保证重映算法的有界性,引入了修补方法。 相似文献
8.
一类基于ENO插值的守恒重映算法 总被引:5,自引:3,他引:2
在大变形流体力学问题的数值模拟中,常常会涉及到计算网格的重分.基于不同网格的物理量传递便是所谓的重映技术.基于ENO插值的思想,发展了一类适用于任意网格的ENO守恒重映算法,并给出了数值结果. 相似文献
9.
10.
11.
12.
二维拉氏自适应流体动力学软件LAD2D,是采用建立在拉氏自适应结构和非结构网格上的有限体积格式,可以计算平面二维和柱对称二维多物质大变形弹塑性流体动力学问题。LAD2D软件系统主要由5部分组成:主控程序、数据模块、前处理模块、主体计算模块、网格模块和后处理模块。其中主体计算采用了结构网格与非结构网格联合使用的拉氏网格体系,计算格式采用了有限体积格式。网格模块包括网格生成、自适应网格加密(AMR技术)和网格重分技术,以及网格改变后物理量守恒重映技术。LAD2D软件系统由主体程序、二维网格生成程序(GRID2D)、二维自适应网格加密程序(AMR2D)、二维自适应程序(ADAPT2D)f[I-维物理量重映程序(REMAP2D)组成。 相似文献
13.
14.
15.
ALE方法在一定程度上克服了纯Lagrange和纯Euler方法的缺点,成为计算大变形问题比较有效的方法。在模拟流固耦合、ICF和磁流体动力学等方面有广泛的应用。对于ALE方法来说,计算网格间物理量的输运(或物理量的重映)是十分关键的。在计算网格间的输运时,用二阶精度的FCT技术,使ALE方法程序计算大变形能力和计算精度都有一定程度的提高。用添加FCT技术后的程序对Taylor杆问题和炸药驱动飞片问题进行数值模拟,Taylor杆问题的计算结果与用LS-DYNA得到的计算结果一致。炸药驱动问题的计算结果与理论分析及实验结果符合较好,比修改前的程序计算变形能力有较大的提高。 相似文献
16.
基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡(WENO)格式.针对任意非结构四边形网格选取重构模板,并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格,会出现非常大的线性权和负权,使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法,对优化后得到的负线性权采用分裂方法进行处理.对于非线性权,提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权,当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证,同时一范数和无穷范数的误差绝对值不依赖于网格质量;具有强间断的数值结果证明了当前格式的有效性. 相似文献
17.
18.
提出一种在自由重映移动网格下的广义黎曼问题方法模拟反应流.该方法基于显式的自由重映移动网格广义黎曼问题的解.为保证在时间和空间上的高精度,应用广义黎曼问题方法构造数值通量.为保证反应区的高分辨率,采用变分法生成自适应移动网格.该方法不仅能够保证网格质量,而且能有效地避免任意拉格朗日—欧拉方法中由于显式重映过程而带来的数值误差.包括CJ爆轰及不稳定爆轰的数值实验说明该格式的精确性和鲁棒性,证明这种移动网格下的二阶广义黎曼问题方法可以较好地捕捉反应流的间断与光滑结构. 相似文献
19.
在流体力学数值模拟中,最基本的有Lagrange方法和Euler方法。Lagrange方法可用来计算多介质系统,能够刻划多介质界面,但网格的扭曲,翻转,长宽比失调等网格大变形是一个突出问题。在Euler方法中,计算网格是固定的,但是,当系统中包含多种介质时,一定会出现在一个Euler网格中包含多种介质的情形,网格中的物理量的处理比较困难。为提高精度.一般将Lagrange方法和Euler方法结合。这时网格最优问题是一个重要的内容。 相似文献
20.
在多维流体动力学计算中,流体运动和计算网格的关系可以分为两种情况。一是Lagrangian方法,即网格跟随流体运动;二是Eulerian方法,即流体流过固定;下动的网格。一般计算网格的运动是任意的。这就对应于任意Lagrangian—Eulerian(ALE)方法。ALE方法的核心是通过调整网格运动,使得数值模拟的精度、效率有所提高。它的主要步骤是:显式Lagrangian步;网格重分,即得到新的计算网格;物理量重映,即将Lagrangian步的计算结果变换到新网格上。在这3步中,较少研究网格重分。数值模拟和网格重分的一个基本前提是网格是合理的,或者说网格不能发生翻转,网格应当是凸的。而Lagrangian步数值模拟会造成网格扭曲,因此在网格重分前进行网格解扭是十分必要的。文中描述了通用的网格解扭、重分算法,使得解扭、重分后的网格有较好的几何品质,同时尽可能接近Lagrangian网格。 相似文献