共查询到18条相似文献,搜索用时 171 毫秒
1.
提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点. 相似文献
2.
3.
提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析. 相似文献
4.
粒子输运方程的线性间断有限元方法 总被引:1,自引:0,他引:1
将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象. 相似文献
5.
6.
精细地捕捉激波及分离流对航空叶型/翼型设计具有重要意义。WENO(Weighted Essentially Non-Oscillatory)格式以其高间断分辨率和低耗散的特性得到广泛青睐。格式的构造都是基于均匀网格,而在实际应用中考虑计算量的问题,通常需要在梯度较大的地方进行局部加密,生成非均匀网格。坐标变换法是一种有效地将格式运用于非均匀网格的方法,而坐标变换的精度对格式特性的影响尚未得知,存在一定的盲目性。针对该问题,基于非均匀网格,开展不同精度坐标变换方法对有限差分WENO格式的精度、耗散影响研究;同时结合工程实际问题,研究不同精度非均匀网格处理方法对激波分辨率和大尺度涡捕捉能力的影响,并给出工程应用建议。 相似文献
7.
高精度、高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义。为了提高三阶WENO-Z格式在极值点处的计算精度,首先通过理论推导给出三阶WENO格式满足收敛精度的充分条件。采用泰勒级数展开的方式,推导确定所构造格式的参数。通过精度测试证明改进格式在光滑流场区域能收敛到三阶精度。选用Sod激波管、Rayleigh-Taylor不稳定性等经典算例证实了提出的改进格式WENO-NN3相较其他格式(WENO-SJ3、WENO-Z3和WENO-N3)具有精度高、耗散低、对流场结构分辨率高的特性。 相似文献
8.
9.
10.
WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。 相似文献
11.
G. Capdeville 《Journal of computational physics》2008,227(4):2430-2454
This paper proposes a new WENO procedure to compute problems containing both discontinuities and a large disparity of characteristic scales.In a one-dimensional context, the WENO procedure is defined on a three-points stencil and designed to be sixth-order in regions of smoothness. We define a finite-volume discretization in which we consider the cell averages of the variable and its first derivative as discrete unknowns. The reconstruction of their point-values is then ensured by a unique sixth-order Hermite polynomial. This polynomial is considered as a symmetric and convex combination, by ideal weights, of three fourth-order polynomials: a central polynomial, defined on the three-points stencil, is combined with two polynomials based on the left and the right two-points stencils.The symmetric nature of such an interpolation has an important consequence: the choice of ideal weights has no influence on the properties of the discretization. This advantage enables to formulate the Hermite interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, nonlinear weights are then defined.To deal with the peculiarities of the Hermite interpolation near discontinuities, we define a new procedure in order for the nonlinear weights to smoothly evolve between the ideal weights, in regions of smoothness, and one-sided weights, otherwise.The resulting scheme is a sixth-order WENO method based on central Hermite interpolation and TVD Runge–Kutta time-integration. We call this scheme the HCWENO6 scheme.Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In these experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems. 相似文献
12.
13.
G. Capdeville 《Journal of computational physics》2008,227(5):2977-3014
This paper proposes a new WENO procedure to compute multi-scale problems with embedded discontinuities, on non-uniform meshes.In a one-dimensional context, the WENO procedure is first defined on a five-points stencil and designed to be fifth-order accurate in regions of smoothness. To this end, we define a finite-volume discretization in which we consider the cell averages of the variable as the discrete unknowns. The reconstruction of their point-values is then ensured by a unique fifth-order polynomial. This optimum polynomial is considered as a symmetric and convex combination, by ideal weights, of four quadratic polynomials.The symmetric nature of the resulting interpolation has an important consequence: the choice of ideal weights has no influence on the accuracy of the discretization. This advantage enables to formulate the interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, non-oscillatory weights are then calculated from the ideal weights.We adapt this procedure for the non-linear weights to maintain the theoretical convergence properties of the optimum reconstruction, whatever the problem considered.The resulting scheme is a fifth-order WENO method based on central interpolation and TVD Runge–Kutta time-integration. We call this scheme the CWENO5 scheme.Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In those experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems. Finally, the new algorithm is directly extended to bi-dimensional problems. 相似文献
14.
15.
The paper extends weighted essentially non-oscillatory (WENO) methods to three dimensional mixed-element unstructured meshes, comprising tetrahedral, hexahedral, prismatic and pyramidal elements. Numerical results illustrate the convergence rates and non-oscillatory properties of the schemes for various smooth and discontinuous solutions test cases and the compressible Euler equations on various types of grids. Schemes of up to fifth order of spatial accuracy are considered. 相似文献
16.
In this work, an adaptive central-upwind 6th-order weighted essentially non-oscillatory (WENO) scheme is developed. The scheme adapts between central and upwind schemes smoothly by a new weighting relation based on blending the smoothness indicators of the optimal higher order stencil and the lower order upwind stencils. The scheme achieves 6th-order accuracy in smooth regions of the solution by introducing a new reference smoothness indicator. A number of numerical examples suggest that the present scheme, while preserving the good shock-capturing properties of the classical WENO schemes, achieves very small numerical dissipation. 相似文献
17.
This paper presents efficient second-order kinetic schemes on unstructured meshes for both compressible unsteady and incompressible steady flows. For compressible unsteady flows, a time-dependent gas distribution function with a discontinuous particle velocity space at a cell interface is constructed and used for the evaluations of both numerical fluxes and conservative flow variables. As a result, a compact scheme on the unstructured meshes is developed. For incompressible steady flows, a continuous second-order gas-kinetic BGK type scheme is presented, for which the time-dependent gas distribution function with a continuous particle velocity is used on unstructured meshes. The efficiency of the schemes lies in the fact that the slopes of the flow variables inside each cell can be constructed using values of the flow variables within that cell only without involving neighboring cells. Therefore, even with the stencil of a first-order scheme, a high resolution method is constructed. Numerical examples are presented which are compared with the benchmark solutions and the experimental measurements. 相似文献
18.
Jiming Wu Zihuan Dai Zhiming Gao Guangwei Yuan 《Journal of computational physics》2010,229(9):3382-3401
In this paper, we employ the so-called linearity preserving method, which requires that a difference scheme should be exact on linear solutions, to derive a nine-point difference scheme for the numerical solution of diffusion equation on the structured quadrilateral meshes. This scheme uses firstly both cell-centered unknowns and vertex unknowns, and then the vertex unknowns are treated as a linear combination of the surrounding cell-centered unknowns, which reduces the scheme to a cell-centered one. The weights in the linear combination are derived through the linearity preserving approach and can be obtained by solving a local linear system whose solvability is rigorously discussed. Moreover, the relations between our linearity preserving scheme and some existing schemes are also discussed, by which a generalized multipoint flux approximation scheme based on the linearity preserving criterion is suggested. Numerical experiments show that the linearity preserving schemes in this paper have nearly second order accuracy on many highly skewed and highly distorted structured quadrilateral meshes. 相似文献