共查询到18条相似文献,搜索用时 250 毫秒
1.
2.
对多车种LWR交通流模型,给出一种半离散中心迎风格式,该格式以五阶WENO-Z重构和半离散中心迎风数值通量为基础.WENO-Z重构方法的引入提高了格式的精度,并保证格式具有基本无振荡的性质.时间的离散采用保持强稳定性的Runge-Kutta方法.通过数值算例验证了格式的有效性. 相似文献
3.
提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析. 相似文献
4.
基于欧拉框架下ADER格式,构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格式.构造r阶LADER格式包括:从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重构节点处守恒量和从1阶到r-1阶的空间导数、求拉氏框架下这些变量的Godunov值,并计算1阶到r-1阶的时间全导数,最后高精度离散积分形式的流通量函数.对光滑流场的模拟表明,LADER格式达到设计的精度;对含强间断的流场模拟表明,数值解在间断附近基本无振荡. 相似文献
5.
迎风紧致格式求解Hamilton-Jacobi方程 总被引:1,自引:1,他引:0
基于Hamilton-Jacobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解. 相似文献
6.
7.
8.
数值求解二维Euler方程的有限体积法(如k-exact,WENO重构、紧致重构等),无一例外地要进行耗时的网格单元上的二维重构.然而这些二维重构最后仅用于确定网格单元边界上高斯积分点处的解值,单元上二维重构似乎并非必需的.因此,文章提出用网格边上的一维重构来取代有限体积法中网格单元上的二维重构,分别在一致矩形网格和非结构三角形网格上发展了基于网格边重构的求解二维Euler方程的新方法,称为降维重构算法.数值算例表明该算法可以计算有强激波的无黏流动问题,且有较高的计算效率. 相似文献
9.
基于近似Riemann解的有限体积ALE方法 总被引:1,自引:0,他引:1
研究二维平面坐标系和二维轴对称坐标系中四边形网格上可压缩流体力学的有限体积ALE(Arbitrary Lagrangian Eulerian)方法.数值方法采用节点中心有限体积法,数值通量采用适用于任意状态方程的HLLC(Harten-Lax-Van Leer-Collela)通量.空间二阶精度通过用WENO(weighted essentially non-oscillatory)方法对原始变量进行重构获得,时间离散采用两步显式Runge-Kutta格式.数值例子显示,方法具有良好的激波分辨能力和高精度的数值逼近能力. 相似文献
10.
龙格库塔间断有限元方法在计算爆轰问题中的应用 总被引:1,自引:1,他引:0
构造求解带源项守恒律方程组的龙格库塔间断有限元(RKDG)方法,并分别结合源项的Strang分裂法和无分裂法数值求解模型守恒律方程和反应欧拉方程.为了和有限体积型WENO方法进行比较,设计计算源项的WENO重构格式.对一维带源项守恒律的计算表明,对于非刚性问题,RKDG方法比有限体积型WENO方法的误差更小;对于刚性问题,RKDG方法对于间断面位置的捕捉更为精确.对于一二维爆轰波问题的计算结果表明,RKDG方法对爆轰波结构的分辨和爆轰波位置的捕捉能力更强. 相似文献
11.
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. 相似文献
12.
We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows. 相似文献
13.
利用函数及其高阶导数值构造五次插值函数近似网格单元内的真实解,改进数值求解双曲类偏微分方程的CIP数值算法。基于之前的一维高阶CIP数值算法思想,不同于利用时间分裂技术,发展了二维高阶CIP数值算法。改进后的算法具有五阶数值精度和显示格式的优点。 相似文献
14.
高阶CIP数值方法及其在相关物理问题中的应用 总被引:1,自引:0,他引:1
利用函数的高阶空间导数值构建其高次插值,得到高阶CIP(Constrained Interpolation Profile)数值算法,并在此基础上模拟研究等离子体物理中著名的伏拉索夫-泊松(Vlasov-Poisson)方程相关物理问题.高阶CIP数值方法具有更高数值精度,从而可以在同等精度的情况下减少计算格点数,加速数值计算速度. 相似文献
15.
16.
A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by the authors of this paper [N.K. Yamaleev, M.H. Carpenter, Third-order energy stable WENO scheme, J. Comput. Phys. 228 (2009) 3025–3047] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes including one sixth-order central scheme; ESWENO schemes up to eighth order are presented in the Appendix. We also develop new weight functions and derive constraints on their parameters, which provide consistency, much faster convergence of the high-order ESWENO schemes to their underlying linear schemes for smooth solutions with arbitrary number of vanishing derivatives, and better resolution near strong discontinuities than the conventional counterparts. 相似文献
17.
18.