共查询到15条相似文献,搜索用时 46 毫秒
1.
给出了求解一维双曲型守恒律的一种半离散三阶中心迎风格式,并利用逐维进行计算的方法将格式推广到二维守恒律。构造格式时利用了波传播的单侧局部速度,三阶重构方法的引入保证了格式的精度。时间方向的离散采用三阶TVD Runge—Kutta方法。本文格式保持了中心差分格式简单的优点,即不需用Riemann解算器,避免了进行特征分解过程。用该格式对一维和二维守恒律进行了大量的数值试验,结果表明本文格式是高精度、高分辨率的。 相似文献
2.
3.
非线性双曲型守恒律的高精度MmB差分格式 总被引:1,自引:0,他引:1
构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的广义G odunov型差分格式。其构造思想是:首先将计算区间划分为若干个互不相交的小区间,再根据精度要求等分小区间,通过各细小区间上的单元平均状态变量,重构各等分小区间交界面上的状态变量,并加以校正;其次,利用近似R iem ann解算子求解细小区间交界面上的数值通量,并结合高阶R unge-K u tta TVD方法进行时间离散,得到了高精度的全离散方法。证明了该格式的Mm B特性。然后,将格式推广到一、二维双曲型守恒方程组情形。最后给出了一、二维Eu ler方程组的几个典型的数值算例,验证了格式的高效性。 相似文献
4.
5.
通过Mac Cormack格式和Warming-Beam的结合,构造了一种非常简单的两步二阶TVD差分格式,该差分格式更适合于使用分量形式差分计算而无须对欧拉方程组进行特征解耦。通过对流体力学方程组的大量数值试验,并与二阶ENO格式进行了比较,充分显示了该格式高精度、高分辨并且极其简单的优良特性。 相似文献
6.
将Jin's的界面方法应用到求解双曲守恒型方程的半离散中心迎风方法中,给出了一种新的求解浅水波方程的半离散中心迎风差分方法。对于源项,不是采用传统的单元均值而是采用单元界面处的值来近似,使所得格式对稳定态的求解是均衡的。且已证明所给的二阶精度的求解格式保持水深的非负性,这一特性使其能够较好的处理干河床问题。使用该方法产生的数值粘性(与O(Δ2r-1)同阶)要比交错的中心格式小(与O(Δx2r/Δt)同阶),而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小,因此适用于稳定态的求解。 相似文献
7.
8.
针对经典的五阶加权本质无振荡(WENO)格式在间断附近耗散过大以及临界点不能保精度的问题,本文提出了一种新的修正模板近似方法。改进了经典五阶WENO-JS格式中各候选子模板上数值通量的二阶多项式逼近,通过加入三次修正项使模板逼近达到四阶精度,并且通过引入可调函数φ使得新的格式具有ENO性质,理论分析新的格式具有保精度特性,通过一系列数值算例说明了新格式的高效性。 相似文献
9.
一类高精度TVD差分格式及其应用 总被引:2,自引:0,他引:2
构造了一维非线性双曲型守恒律的一个新的高精度、高分辨率的守恒型TvD差分格式。其构造思想是:首先,将计算区间划分为若干个互不相交的小区间,再根据精度要求等分小区间,通过各细小区间上的单元平均状态变量,重构各细小区间交界面上的状态变量,并加以校正;其次,利用近似Riemann解计算细小区间交界面上的数值通量,并结合高阶Runge—Kutta TVD方法进行时间离散,得到了高精度的全离散方法。证明了该格式的TVD特性。该格式适合于使用分量形式计算而无须进行局部特征分解。通过计算几个典型的问题,验证了格式具有高精度、高分辨率且计算简单的优点。 相似文献
10.
11.
直接把Nessyahu和Tadmor^[1,2]的思想推广到三维非线性双曲型守恒律情形,以交错形式Lax—Friedrichs格式为基本模块,使用二阶分片线性逼近代替一阶分片常数逼近,减少了Lax—Friedrichs格式的过多数值粘性,通过对混合导数离散形式的适当处理,构造了一类不须解Riemann问题、具有时空二阶精度高分辨率的MmB差分格式。这些差分格式很容易推广到向量系统中去。最后,一些数值模拟计算结果也证明了这些差分格式的有效性。 相似文献
12.
It has been shown that one can generate a class of nontrivial conservation laws for second-order partial differential equations using some recent results dealing with the action of any Lie–Bäcklund symmetry generator of the equivalentfirst-order system on the respective conservation law. These conservedvectors are nonlocal as they are constructed from associatednonlocal symmetries of the partial differential equation. The method canbe successfully extended to association with genuine nonlocal(potential) symmetries. However, it usually involves solving moredifficult systems of partial differential equations which may not alwaysbe easy to uncouple. 相似文献
13.
In this work we present an upwind‐based high resolution scheme using flux limiters. Based on the direction of flow we choose the smoothness parameter in such a way that it leads to a truly upwind scheme without losing total variation diminishing (TVD) property for hyperbolic linear systems where characteristic values can be of either sign. Here we present and justify the choice of smoothness parameters. The numerical flux function of a high resolution scheme is constructed using wave speed splitting so that it results into a scheme that truly respects the physical hyperbolicity property. Bounds are given for limiter functions to satisfy TVD property. The proposed scheme is extended for non‐linear problems by using the framework of relaxation system that converts a non‐linear conservation law into a system of linear convection equations with a non‐linear source term. The characteristic speed of relaxation system is chosen locally on three point stencil of grid. This obtained relaxation system is solved using composite scheme technique, i.e. using a combination of proposed scheme with the conservative non‐standard finite difference scheme. Presented numerical results show higher resolution near discontinuity without introducing spurious oscillations. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
Let v=v(x) be a non-trivial bounded steady solution of a viscous scalar conservation law u
t+f(u)
x
=u
xx on a half-line R+, with a Dirichlet boundary condition. The semi-group of this IBVP is known to be contractive for the distance d(u, u)u–u1 induced by L
1(R+). We prove here that v is asymptotically stable with respect to d: if u
0–vL
1, then u(t)–v10 as t+. When v is a constant, we show that this property holds if and only if f(v)0. These results complement our study of the Cauchy problem [2]. 相似文献
15.
Traditional mathematical models of multiphase flow in porous media use a straightforward extension of Darcys equation. The key element of these models is the appropriate formulation of the relative permeability functions. It is well known that for one-dimensional flow of three immiscible incompressible fluids, when capillarity is neglected, most relative permeability models used today give rise to regions in the saturation space with elliptic behavior (the so-called elliptic regions). We believe that this behavior is not physical, but rather the result of an incomplete mathematical model. In this paper we identify necessary conditions that must be satisfied by the relative permeability functions, so that the system of equations describing three-phase flow is strictly hyperbolic everywhere in the saturation triangle. These conditions seem to be in good agreement with pore-scale physics and experimental data. 相似文献