共查询到20条相似文献,搜索用时 140 毫秒
1.
论及高分辨分子动力学通向量分裂(KFVS)有限体积方法的推广。在方法中提出了适当修改Maxwell平衡分布用以修复Euler方程。基于熟知的Euler方程与Boltzm方程的关系,提出了一类求解多分量Euler方程的高分辨分子动力学通向量分裂(KFVS)有限体积方法。应用该方法不需要求解任何Riemann问题或求解附加的非守恒压力方程也不需要任何非守恒修正。数值计算表明,数值解在物质界面附近无振荡,激波速度也正确,显示出方法的高精度及其稳健性。 相似文献
2.
3.
4.
5.
龙格库塔间断有限元方法在计算爆轰问题中的应用 总被引:1,自引:1,他引:0
构造求解带源项守恒律方程组的龙格库塔间断有限元(RKDG)方法,并分别结合源项的Strang分裂法和无分裂法数值求解模型守恒律方程和反应欧拉方程.为了和有限体积型WENO方法进行比较,设计计算源项的WENO重构格式.对一维带源项守恒律的计算表明,对于非刚性问题,RKDG方法比有限体积型WENO方法的误差更小;对于刚性问题,RKDG方法对于间断面位置的捕捉更为精确.对于一二维爆轰波问题的计算结果表明,RKDG方法对爆轰波结构的分辨和爆轰波位置的捕捉能力更强. 相似文献
6.
对于一个原子平均体积为V,温度为T的热力学系统,体系的Helmboltz自由能可以写为F(V,T)=Ec(V)+Fion(V,T)+Fel(V,T)+Fman(V,T)其中Ec为0-K冷能。对于其中的电子热激发贡献Fel,目前流行的有三种计算方案,即:Moruzzi的Debye-Grueneisen方案、Moriarty的MGPT方案和Wasserman的CELL模型方案。Debye-Gruene 相似文献
7.
8.
9.
用有限元方法求解双曲守恒律 总被引:1,自引:1,他引:0
分片线性插值有限元给出了求解双曲守恒律的计算方法。有别于不连续有限元方法求解双曲守恒律在相邻单元边界上求Riemann解,利用双曲守恒律的Hamilton-Jacobi方程形式,直接应用有限元求解,在CFL下,证明了计算格式满足极大值原理,并且是TVD格式。数值例子在文后给出。此外,方法推广到流体力学方程组和高维问题,将在另文中邓以讨论。 相似文献
10.
11.
The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate
下载免费PDF全文
![点击此处可从《中国物理 B》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume (RKCV) discontinuous finite element method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discontinuous Galerkin (RKDG) method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm. 相似文献
12.
In [16], [17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions), it is more difficult to design high order schemes which do not produce negative density or pressure. In this paper, we first show that our framework to construct positivity-preserving high order schemes in [16], [17] can also be applied to Euler equations with a general equation of state. Then we discuss an extension to Euler equations with source terms. Numerical tests of the third order Runge–Kutta DG (RKDG) method for Euler equations with different types of source terms are reported. 相似文献
13.
Euler方程某些问题的解具有自相似特点,可以使用更准确的方法求解.提出了两种数值方法,分别称为自相似和准自相似方法,新方法可以使用现有守恒律方程的数值格式,无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程,一维计算结果显示数值解基本等同于精确解,二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程,新方法可以求解不具有自相似性但接近自相似的问题,并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究,高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示. 相似文献
14.
We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or discontinuous Galerkin (DG) method with first order Euler forward time discretization solving one-dimensional scalar conservation laws. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. It is straightforward to extend the method to two and higher dimensions on rectangular meshes. We also show that the same limiter can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. Numerical tests for both the WENO finite volume scheme and the DG method are reported. 相似文献
15.
A Switch Function-Based Gas-Kinetic Scheme for Simulation of Inviscid and Viscous Compressible Flows
下载免费PDF全文
![点击此处可从《advances in applied mathematics and mechanics.》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Yu Sun Shu Chang Liming Yang & C. J. Teo 《advances in applied mathematics and mechanics.》2016,8(5):703-721
In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented
for the simulation of inviscid and viscous compressible flows. With the finite
volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS
is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by
the conventional smooth function approximation. Unlike the traditional gas-kinetic
scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS),
the numerical dissipation is controlled with a switch function in the present scheme.
That is, the numerical dissipation is only introduced in the region around strong shock
waves. As a consequence, the present SF-GKS can well capture strong shock waves
and thin boundary layers simultaneously. The present SF-GKS is firstly validated by
its application to the inviscid flow problems, including 1-D Euler shock tube, regular
shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous
transonic and hypersonic flow problems. Good agreement between the present
results and those in the literature verifies the accuracy and robustness of SF-GKS. 相似文献
16.
17.
用改进的耦合型Level Set方法计算一维双介质可压缩流动 总被引:2,自引:1,他引:1
用带有虚拟流体(Ghost Fluid)修正的Level Set方法计算了一维可压缩双介质流动,把描述流动的Euler方程和描述流体界面运动的Level Set方程耦合起来,得到一个整体的守恒律系统,应用高分辨率差分格式求解;为了解决流体界面附近的数值跳动问题,在界面附近引入了虚拟流体方法的Isobaric修正,并给出了算例. 相似文献
18.
《Journal of computational physics》2006,211(1):154-178
In this paper, the spectral volume method is extended to the two-dimensional Euler equations with curved boundaries. It is well-known that high-order methods can achieve higher accuracy on coarser meshes than low-order methods. In order to realize the advantage of the high-order spectral volume method over the low order finite volume method, it is critical that solid wall boundaries be represented with high-order polynomials compatible with the order of the interpolation for the state variables. Otherwise, numerical errors generated by the low-order boundary representation may overwhelm any potential accuracy gains offered by high-order methods. Therefore, more general types of spectral volumes (or elements) with curved edges are used near solid walls to approximate the boundaries with high fidelity. The importance of this high-order boundary representation is demonstrated with several well-know inviscid flow test cases, and through comparisons with a second-order finite volume method. 相似文献
19.
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. 相似文献