共查询到18条相似文献,搜索用时 125 毫秒
1.
提出基于细分和数值积分思想的一种离散的守恒重映方法——质点重映方法.密度分布可采用一阶精度的分片常数分布,或二阶精度的分片线性分布.分片线性密度分布函数采用面平均方法构造.重映过程中,借助四边形辅助网格,实现了交错网格节点量的重映.质点重映方法既适用于结构网格,也适用于非结构网格,且不要求新旧网格之间一一对应.数值结果表明,一阶精度重映算法健壮性好,但会产生较大的扩散效应;二阶精度重映算法可较好地保持密度分布的特性,但存在单调性问题.为改善二阶精度重映方法单调性,将结构网格质量守恒调整算法推广到非结构网格上,以限制新网格的质量密度.给出了一些重映的例子,并进行了误差分析. 相似文献
2.
介绍了二维非结构网格上的守恒重映算法,重点是基于SFB/DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法,包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法,不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映,介绍了SALE和HIS算法。此外,为保证重映算法的有界性,引入了修补方法。 相似文献
3.
4.
5.
随着流体问题的复杂和数值模拟高精度的要求,常常需要不同方法的耦合或不同网格体系的转化、叠加.只要改动网格或不同体系耦合(网格重构、细化与粗化、结构网格与非结构网格的互换)就涉及两种体系网格物理量的耦合或重映.将曲面拟合或插值的方法引入到散乱物理量重映计算中,给出几种物理量重映的算法.这些算法不受网格体系的限制,简单易用. 相似文献
6.
7.
一类基于ENO插值的守恒重映算法 总被引:5,自引:3,他引:2
在大变形流体力学问题的数值模拟中,常常会涉及到计算网格的重分.基于不同网格的物理量传递便是所谓的重映技术.基于ENO插值的思想,发展了一类适用于任意网格的ENO守恒重映算法,并给出了数值结果. 相似文献
8.
Lagrange方法中,当流场发生大变形时,跟踪流体运动的Lagrange网格发生扭曲,使计算无法进行下去,此时必须重分网格,把网格修复成较好的形状。另外,网格自适应技术中的重构、合并与加密,以及同一问题不同程序相继计算的连接,并行计算中相邻块边界区域的数据传递等,这些情况都需要利用旧网格上的物理量来确定新网格上的物理量,是一个物理量重映过程。质点重映方法是基于物理上守恒规律的一种离散的物理量守恒映射方法,既可实现分片常数分布的一阶精度重映计算,又可实现分片线性分布的二阶精度重映计算。这种方法可严格保证守恒量的守恒性,且可以实现任意多边形网格以及节点上物理量的守恒重映。但是,基于分片线性分布的二阶精度重映方法,如果新网格的守恒量没有进行保界调整,那么相应的强度量有可能在其局部的限制范围之外,破坏了原网格物理量的单调性。因而,对二阶精度的质点重映方法进行了进一步研究。在分片线性分布的基础上,将基于结构网格的保界算法扩展到非结构网格上,给出了二阶保界的质点守恒重映方法。 相似文献
9.
10.
任意网格重映的样条逼近算法 总被引:2,自引:1,他引:1
在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果. 相似文献
11.
12.
13.
14.
Stéphane Del Pino 《Journal of computational physics》2011,230(5):1793-1821
In this paper, we present a method to compute compressible flows in 2D. It uses two steps: a Lagrangian step and a metric-based triangular mesh adaptation step. Computational mesh is locally adapted according to some metric field that depends on physical or geometrical data. This mesh adaptation step embeds a conservative remapping procedure to satisfy consistency with Euler equations. The whole method is no more Lagrangian.After describing mesh adaptation patterns, we recall the metric formalism. Then, we detail an appropriate remapping procedure which is first-order and relies on exact intersections.We give some hints about the parallel implementation. Finally, we present various numerical experiments which demonstrate the good properties of the algorithm. 相似文献
15.
提出一种在自由重映移动网格下的广义黎曼问题方法模拟反应流.该方法基于显式的自由重映移动网格广义黎曼问题的解.为保证在时间和空间上的高精度,应用广义黎曼问题方法构造数值通量.为保证反应区的高分辨率,采用变分法生成自适应移动网格.该方法不仅能够保证网格质量,而且能有效地避免任意拉格朗日—欧拉方法中由于显式重映过程而带来的数值误差.包括CJ爆轰及不稳定爆轰的数值实验说明该格式的精确性和鲁棒性,证明这种移动网格下的二阶广义黎曼问题方法可以较好地捕捉反应流的间断与光滑结构. 相似文献
16.
《Journal of computational physics》2003,184(1):266-298
An accurate conservative interpolation (remapping) algorithm is an essential component of most arbitrary Lagrangian–Eulerian (ALE) methods. In this paper we describe a local remapping algorithm for a positive scalar function. This algorithm is second-order accurate, conservative, and sign preserving. The algorithm is based on estimating the mass exchanged between cells at their common interface, and so is equally applicable to structured and unstructured grids. We construct the algorithm in a series of steps, clearly delineating the assumptions and errors made at each step. We validate our theory with a suite of numerical examples, analyzing the results from the viewpoint of accuracy and order of convergence. 相似文献
17.
A novel subgridding algorithm for the Method of Lines is presented. It enables to refine a mesh in subgridded regions by any positive integer number. High-accuracy of the proposed algorithm is achieved by second-order finite difference operators. The interpolation of the missing field components at the interface coarse/fine grid is performed with second-order accuracy as well. Numerical results are provided to check the accuracy of the proposed formulation. 相似文献
18.
In this paper, a new flux-based one-step hybrid remapping method for multi-material arbitrary Lagrangian–Eulerian (ALE) approach is introduced. In the vicinity of material interfaces, the swept region is intersected with pure material polygons in the Lagrangian mesh to construct the material fluxes. Far from interfaces, the fluxes are constructed in a standard swept-region manner without intersections. This method is conservative, second-order accurate and linearity-preserving (in case of straight material interfaces), and faster than method based on intersections, as shown on selected numerical examples. 相似文献