共查询到19条相似文献,搜索用时 305 毫秒
1.
基于近似Riemann解的有限体积ALE方法 总被引:1,自引:0,他引:1
研究二维平面坐标系和二维轴对称坐标系中四边形网格上可压缩流体力学的有限体积ALE(Arbitrary Lagrangian Eulerian)方法.数值方法采用节点中心有限体积法,数值通量采用适用于任意状态方程的HLLC(Harten-Lax-Van Leer-Collela)通量.空间二阶精度通过用WENO(weighted essentially non-oscillatory)方法对原始变量进行重构获得,时间离散采用两步显式Runge-Kutta格式.数值例子显示,方法具有良好的激波分辨能力和高精度的数值逼近能力. 相似文献
2.
3.
4.
基于欧拉框架下ADER格式,构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格式.构造r阶LADER格式包括:从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重构节点处守恒量和从1阶到r-1阶的空间导数、求拉氏框架下这些变量的Godunov值,并计算1阶到r-1阶的时间全导数,最后高精度离散积分形式的流通量函数.对光滑流场的模拟表明,LADER格式达到设计的精度;对含强间断的流场模拟表明,数值解在间断附近基本无振荡. 相似文献
5.
提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析. 相似文献
6.
给出数值求解一维双曲守恒律方程的新方法——龙格-库塔控制体积间断有限元方法(RKCVDFEM),其中空间离散基于控制体积有限元方法,时间离散基于二阶TVB Runge-Kutta技术,有限元空间选取为分段线性函数空间.理论分析表明,格式具有总变差有界(TVB)的性质,而且空间和时间离散形式上具有二阶精度.数值算例表明,数值解收敛到熵解并且对光滑解的收敛阶是最优的,优于龙格-库塔间断Galerkin方法(RKDGM)的计算结果. 相似文献
7.
8.
9.
10.
11.
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. 相似文献
12.
文章给出了一种真正多维的HLL Riemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量, 其中对流通量的计算采用迎风格式, 压力通量的计算采用HLL格式, 且将HLL格式的耗散项中的密度差用压力差代替, 从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性, 分别计算了网格界面中点和角点上的数值通量, 并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度, 计算中采用基于SDWLS梯度的线性重构获得2阶空间精度, 而时间离散采用2阶保强稳Runge-Kutta方法.数值实验表明, 相比于传统的一维HLL格式, 文章的真正多维HLL格式具有能够分辨接触间断, 以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同, 真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象. 相似文献
13.
Francesco Valentini Pierluigi Veltri Andr Mangeney 《Journal of computational physics》2005,210(2):730-751
We present a numerical algorithm for the solution of the Vlasov–Poisson system of equations, in the magnetized case. The numerical integration is performed using the well-known “splitting” method in the electrostatic approximation, coupled with a finite difference upwind scheme; finally the algorithm provides second order accuracy in space and time. The cylindrical geometry is used in the velocity space, in order to describe the rotation of the particles around the direction of the external uniform magnetic field.Using polar coordinates, the integration of the Vlasov equation is very simplified in the velocity space with respect to the cartesian geometry, because the rotation in the velocity cartesian space corresponds to a translation along the azimuthal angle in the cylindrical reference frame. The scheme is intrinsically symplectic and significatively simpler to implement, with respect to a cartesian one. The numerical integration is shown in detail and several conservation tests are presented, in order to control the numerical accuracy of the code and the time evolution of the entropy, strictly related to the filamentation problem for a kinetic model, is discussed. 相似文献
14.
Laurent Gosse 《Journal of Quantitative Spectroscopy & Radiative Transfer》2011,112(12):1995-2012
An original well-balanced (WB) Godunov scheme relying on an exact Riemann solver involving a non-conservative (NC) product is developed. It is meant to solve accurately the time-dependent one-dimensional radiative transfer equation in the discrete ordinates approximation with an arbitrary even number of velocities. The collision term is thus concentrated onto a discrete lattice by means of Dirac masses; this induces steady contact discontinuities which are integral curves of the stationary problem. One solves it by taking advantage of the method of elementary solutions mainly developed by Case, Zweifel and Cercignani. This approach produces a rather simple scheme that compares advantageously to standard existing upwind schemes, especially for the decay in time toward a Maxwellian distribution. It is possible to reformulate this scheme in order to handle properly the parabolic scaling in order to generate a so-called asymptotic-preserving (AP) discretization. Consistency with the diffusive approximation holds independently of the computational grid. Several numerical results are displayed to show the realizability and the efficiency of the method. 相似文献
15.
为实现SPH二阶粒子近似,讨论光滑函数应满足的充分条件,给出对应的光滑函数形式,分析二阶粒子近似解的唯一性,并进行数值验证.结果表明:和标准SPH方法相比,在周期边界条件下,随着粒子数加密,SPH方法二阶粒子近似的L1误差减小且L1阶趋近2.0. 相似文献
16.
17.
二维非结构自适应多重网格的Euler方程解 总被引:1,自引:0,他引:1
将原来在结构网格中提出的AUSM+(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中。利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡。为提高计算效率,采用了自适应多重网格法。在生成非结构网格时,采用了一种新颖的数据结构——数组链接表。最后给出了NACA0012翼型和带10%鼓包的管道流动的几个Euler方程的算例,并进行了分析比较。 相似文献
18.
An E-CUSP (energy-convective upwind and split pressure) scheme is developed to solve the equations of magnetohydrodynamics. Fifth order WENO reconstructions are employed to calculate the fluxes in order to achieve high order spacial accuracy. A characteristic speed of sound by averaging the fast wave speed and the acoustic speed of sound is suggested to evaluate the Mach number, which will yield robust and accurate solutions. The numerical experiments have demonstrated the accuracy and the capability of the new scheme to capture complex interactions of multiple shocks and vortices. 相似文献
19.
HE ZhiWei LI XinLiang FU DeXun & MA YanWen The Key Laboratory of High Temperature Gas Dynamics 《中国科学:物理学 力学 天文学(英文版)》2011,(3)
Based on an upwind compact difference scheme and the idea of monotonicity-preserving, a 5th order monotonicity-preserving upwind compact difference scheme (m-UCD5) is proposed. The new difference scheme not only retains the advantage of good resolution of high wave number but also avoids the Gibbs phenomenon of the original upwind compact difference scheme. Compared with the classical 5th order WENO difference scheme, the new difference scheme is simpler and small in diffusion and computation load. By emplo... 相似文献