共查询到16条相似文献,搜索用时 46 毫秒
1.
2.
3.
4.
5.
6.
一类基于ENO插值的守恒重映算法 总被引:5,自引:3,他引:2
在大变形流体力学问题的数值模拟中,常常会涉及到计算网格的重分.基于不同网格的物理量传递便是所谓的重映技术.基于ENO插值的思想,发展了一类适用于任意网格的ENO守恒重映算法,并给出了数值结果. 相似文献
7.
任意网格重映的样条逼近算法 总被引:2,自引:1,他引:1
在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果. 相似文献
8.
在低磁场雷诺数条件下,基于电势泊松方程,发展了交错网格下可以精确计算电流和洛伦兹力(电磁力)的相容守恒格式。采用压力为变量的原始变量法求解不可压缩Navier-Stokes方程,所计算的电流满足电荷守恒定律,所计算的电磁力满足动量守恒定律。对金属流体在Hartmann数50~5000范围内验证了格式的精确性。交错网格下相容守恒格式的发展为后续MHD稳定性分析、湍流的大涡模拟及直接数值模拟提供很好的选择。 相似文献
9.
Lagrange方法中,当流场发生大变形时,跟踪流体运动的Lagrange网格发生扭曲,使计算无法进行下去,此时必须重分网格,把网格修复成较好的形状。另外,网格自适应技术中的重构、合并与加密,以及同一问题不同程序相继计算的连接,并行计算中相邻块边界区域的数据传递等,这些情况都需要利用旧网格上的物理量来确定新网格上的物理量,是一个物理量重映过程。质点重映方法是基于物理上守恒规律的一种离散的物理量守恒映射方法,既可实现分片常数分布的一阶精度重映计算,又可实现分片线性分布的二阶精度重映计算。这种方法可严格保证守恒量的守恒性,且可以实现任意多边形网格以及节点上物理量的守恒重映。但是,基于分片线性分布的二阶精度重映方法,如果新网格的守恒量没有进行保界调整,那么相应的强度量有可能在其局部的限制范围之外,破坏了原网格物理量的单调性。因而,对二阶精度的质点重映方法进行了进一步研究。在分片线性分布的基础上,将基于结构网格的保界算法扩展到非结构网格上,给出了二阶保界的质点守恒重映方法。 相似文献
10.
介绍了二维非结构网格上的守恒重映算法,重点是基于SFB/DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法,包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法,不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映,介绍了SALE和HIS算法。此外,为保证重映算法的有界性,引入了修补方法。 相似文献
11.
12.
13.
14.
15.
为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象,发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式,避免在单元内进行速度平均,从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比,网格的边界通量由节点处的多维黎曼求解器提供,克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增,给出稀疏波出现的合理判据,从而保证了热力学关系式的满足.数值实验表明:该方法能很好地消除稀疏波的过热现象,同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力. 相似文献
16.
Solution Reconstruction on Unstructured Tetrahedral Meshes Using P1 -Conservative Interpolation 下载免费PDF全文
Biao Peng Chunhua Zhou & Junqiang Ai 《advances in applied mathematics and mechanics.》2016,8(5):847-870
This paper extends an algorithm of P1-conservative interpolation on triangular
meshes to tetrahedral meshes and thus constructs an approach of solution reconstruction
for three-dimensional problems. The conservation property is achieved
by local mesh intersection and the mass of a tetrahedron of the current mesh is calculated
by the integral on its intersection with the background mesh. For each current
tetrahedron, the overlapped background tetrahedrons are detected efficiently. A mesh
intersection algorithm is proposed to construct the intersection of a current tetrahedron
with the overlapped background tetrahedron and mesh the intersection region
by tetrahedrons. A localization algorithm is employed to search the host units in background
mesh for each vertex of the current mesh. In order to enforce the maximum
principle and avoid the loss of monotonicity, correction of nodal interpolated solution
on tetrahedral meshes is given. The performance of the present solution reconstruction
method is verified by numerical experiments on several analytic functions and the
solution of the flow around a sphere. 相似文献