A New Well-Balanced Finite Volume CWENO Scheme for Shallow Water Equations over Bottom Topography

In this article,we develop a new well-balanced finite volume central weighted essentially non-oscillatory(CWENO)scheme for one-and two-dimensional shallow water equations over uneven bottom.The well-balanced property is of paramount importance in practical applications,where many studied phenomena can be regarded as small perturbations to the steady state.To achieve the well-balanced property,we construct numerical fluxes by means of a decomposition algorithm based on a novel equilibrium preserving reconstruction procedure and we avoid applying the traditional hydrostatic reconstruction technique accordingly.This decomposition algorithm also helps us realize a simple source term discretization.Both rigorous theoretical analysis and extensive numerical examples all verify that the proposed scheme maintains the well-balanced property exactly.Furthermore,extensive numerical results strongly suggest that the resulting scheme can accurately capture small perturbations to the steady state and keep the genuine high-order accuracy for smooth solutions at the same time.

求解双曲守恒律方程的五阶紧凑CWENO格式

提出一种新的求解双曲守恒律方程的五阶紧凑CWENO格式 ,基于Godunov方法的思想 ,该格式中每个模板上的重构多项式是单元均值意义下的插值多项式 .对空间方向上的重构 ,引入了五阶紧凑CWENO格式 ,重构多项式是基于不同模板上的插值多项式的凸组合 .该方法的核心是 :首先构造一个最优 4次多项式 ,通过凸组合的形式使解在光滑区域达到五阶精度 ,在间断区域 ,凸组合的权值会自适应地选择单个模板上的三阶插值多项式 ,从而避免了伪震荡的产生 (WENO思想 ) .这种新的五阶重构格式是基于非常紧凑的五点模板构造的 .最后此格式的精确性、稳定性及高分辨率性通过一维算例给以验证 .

基于五阶CWENO重构的中心迎风格式

提出一种五阶CWENO重构,将其与中心迎风数值通量相结合,得到一种求解双曲型守恒律方程的五阶中心迎风格式。该格式构造简单,无需Riemann求解器,易于编程实现。运用该格式求解标量守恒律方程和Euler方程组,数值结果表明,该格式具有五阶精度,分辨率高,能准确计算出解的各种复杂结构。

Application of a fourth-order relaxation scheme to hyperbolic systems of conservation laws

A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourthorder central weighted essentially nonoscillatory （CWENO） reconstruction for one-dimensional cases, which is generalized to two-dimensional cases by the dimension-by-dimension approach. The large stability domain Runge-Kutta-type solver ROCK4 is used for time integration. The resulting method requires neither the use of Riemann solvers nor the computation of Jacobians and therefore it enjoys the main advantage of the relaxation schemes. The high accuracy and high-resolution properties of the present method are demonstrated in one- and two-dimensional numerical experiments.

求解一维理想磁流体方程的5阶紧凑CWENO中心迎风格式

通过结合求解双曲型守恒律的5阶紧凑CWENO格式和半离散中心迎风格式,推广应用于求解一维理想磁流体力学方程,得到计算一维理想磁流体力学方程的5阶紧凑CWENO中心迎风格式.

高阶保号熵稳定格式

高阶熵稳定格式构造的一个核心任务是如何保证熵变量在进行高阶重构前后的符号不变.该文构造了高阶保号熵稳定格式(熵守恒通量采用Fjordholm方案,耗散部分的熵变量采用三阶紧致CWENO重构),证明了基于该重构的熵变量在跳跃间断处满足保号性.数值结果表明,该格式达到三阶精度、分辨率高、鲁棒性强且无振荡产生.

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。

一维浅水波方程的高分辨率熵相容算法

针对一维浅水波方程的数值求解问题,构造了三阶、四阶高分辨率熵相容算法:数值通量函数是满足总熵衰变特点的熵相容通量,时间方向是具有强稳定特点的优化三阶Runge-Kutta算法,空间方向分别采用三阶、四阶CWENO(Central Weighted Essentially Non-oscillatory)重构法进行了离散,通过若干算例研究了新算法的性能.结果表明:新算法能够准确捕捉激波、稀疏波,计算结果与准确解符合很好,且在强稀疏波算例中新算法能够有效避免膨胀激波现象的产生.新算法是熵相容的,易于编程实现,且计算结果可靠,便于向高维推广.

Hamilton-Jacobi方程的高分辨率解法

Ham ilton-Jacob i方程经常应用于控制论和微分对策等方面.它与双曲型守恒律有非常紧密的联系,在一维的情形下,两者几乎完全相同.在多维的空间中,类似的情形依然存在.因此,可以通过这种联系来设计不同的近似算法来求解Ham ilton-Jacob i方程.本文利用CWENO算法成功地求解了包括控制最优化问题在内的许多标量问题,结果表明,这种算法在解的光滑区域具有很高的精度,并能很好地解决具有不连续偏导数的计算问题,数值算例结果也表明这种算法是收敛的.

求解复杂交通流模型的低耗散中心迎风格式

针对非均匀道路上的多车种LWR交通流模型,提出一种低耗散中心迎风格式。以4阶中心加权基本无震荡重构和低耗散中心迎风数值通量为基础,通过构造不同形式的全局光滑因子及增大非光滑模板对应的非线性权重优化数值格式的耗散特性,并采用Runge-Kutta方法对半离散数值格式在时间方向上进行离散使其保持4阶精度。对非均匀道路上多车种LWR交通流模型的车道数变化和交通信号灯控制问题进行数值模拟,结果表明该格式具有4阶求解精度,且分辨率高。

A CHARACTERISTIC-BASED FINITE VOLUME SCHEME FOR SHALLOW WATER EQUATIONS

We propose a new characteristic-based finite volume scheme combined with the method of Central Weighted Essentially Non-Oscillatory （CWENO） reconstruction and characteristics, to solve shallow water equations. We apply the scheme to simulate dam-break problems. A number of challenging test cases are considered, such as large depth differences even wet/dry bed. The numerical solutions well agree with the analytical solutions. The results demonstrate the desired accuracy, high-resolution and robustness of the presented scheme.

SOLUTION OF 2D SHALLOW WATER EQUATIONS BY GENUINELY MULTIDIMENSIONAL SEMI-DISCRETE CENTRAL SCHEME

A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization and the optimal third-order Strong Stability Preserving （SSP） Runge-Kutta method for time integration. The third-order compact Central Weighted Essentially Non-Oscillatory （CWENO） reconstruction was adopted to guarantee the non-oscillatory behavior of the presented scheme and improve the resolution. Two kinds of source terms were considered in this work. They were evaluated using different approaches. The resulting scheme does not require Riemann solvers or characteristic decomposition, hence it retains all the attractive features of central schemes such as simplicity and high resolution. To evaluate the performance of the presented scheme, several numerical examples were tested. The results demonstrate that our method is efficient, stable and robust.

交通流LWR模型的高分辨率熵相容算法

针对交通流LwR模型的数值求解，基于Ismail熵相容理论，推导出相应的熵相容格式；并在空间方向采用CWENO重构，时间方向采用优化三阶Runge—Kutta方法改进该算法，进而得到高精度、高分辨率及数值稳定的熵相容方法。最终将新构造的数值求解格式应用于多个实际交通现象的数值求解当中。数值结果显示：熵相容算法及高分辨率熵相容算法能够避免稀疏一激波现象地产生，在间断区域没有非物理振荡；且高分辨率熵相容算法数值稳定性强，比已有的几种算法分辨率高（激波仅需2～3个过渡点，稀疏波、接触间断与参考值吻合很好），新算法具有TVB（Total Variations Bounded）性，是一种具有实际应用价值的算法。