一维多介质可压缩流动的守恒型界面追踪方法

本文针对一维问题的ProntTracking方法,提出了一种较易实现的守恒型界面追踪方法.利用双波近似求解Riemann问题来确定界面处的数值通量,在固定的网格上采用统一的有限体积格式进行内点和交界面点的计算,通过守恒插值以及守恒量的重新分配,保证数值解在全场实现一致守恒,将该方法应用于一维多介质可压缩流动的模拟,给出了满意的数值模拟结果.     

带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解的极限

研究了带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解的极限.由于非齐次项的影响,带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解不再是自相似的.当压力和磁感强度同时消失时,它的解会收敛到零压流输运方程组的Riemann解,解中会出现δ-激波和真空现象.同时研究还得到了仅当磁感强度消失时,它的解会收敛到非齐次广义Chaplygin气体Euler方程组的Riemann解,并且解中只出现δ-激波.     

一维定常型对流占优扩散方程的一类迎风有限体积格式

本文针对一维定常型对流占优扩散方程提出了一类迎风有限体积格式 .该格式对对流项具有二阶精度 ,对扩散项保持一阶精度 ,符合对流占优扩散问题强对流、弱扩散的特点 .     

一种基于TV分裂的真正多维Riemann解法器

给出了一种真正多维的HLL Riemann解法器.采用TV(Toro-Vázquez)分裂将通量分裂成对流通量和压力通量,其中对流通量的计算采用类似于AUSM格式的迎风方法,压力通量的计算采用波速基于压力系统特征值的HLL格式,并将HLL格式耗散项中的密度差用压力差代替,来克服传统的HLL格式不能分辨接触间断的缺点.为了实现数值格式真正多维的特性,分别计算网格界面中点和角点上的数值通量,并且采用Simpson公式加权中点和角点上的数值通量来得到网格界面上的数值通量.采用基于SDWLS(solution dependent weighted least squares)梯度的线性重构来获得空间的二阶精度,时间离散采用二阶Runge-Kutta格式.数值实验表明,相比于传统的一维HLL格式,该文的真正多维HLL格式具有能够分辨接触间断,消除慢行激波波后振荡以及更大的时间步长等优点.并且,与其他能够分辨接触间断的格式（例如HLLC格式）不同的是,真正多维的HLL格式在计算二维问题时不会出现数值激波不稳定现象.     

带源项浅水波方程的高分辨率熵稳定格式

提出了一种求解带源项浅水波方程的熵稳定格式.新格式利用通量限制函数将一阶熵稳定格式和高阶熵守恒格式结合,具有熵守恒格式和熵稳定格式的优点:在解的光滑区域具有高精度,在解的间断区域避免了非物理现象的产生,同时可以准确地捕捉激波,从而达到高分辨率的效果.利用新格式计算了一维和二维的经典算例,数值结果表明,新格式是模拟带源项浅水波方程的理想方法.     

二相流计算的一种差分算法

考虑两相流的力学行为,忽略相间的耗散作用,建立了Euler型的基本控制方程．状态方程采用刚性状态方程．基于Abgrall提出的准则,在流动区域内,对可压两相流提出了一个精度较高的Euler型数值方法,数值格式是Godunov型格式,对守恒型和非守恒型方程采用HLLC型和Lax-Friedrichs型近似Riemann解算器,引入了速度驰豫和压强驰豫过程来代替两相间的相互作用．在一维情形下给出数值算例,并且和Saurel的算例进行了比较,结果表明该算法不但精确而且稳定,且在间断处没有数值振荡．     

双币种期权定价模型的一个新ADI并行差分方法

双币种期权是一种重要的金融衍生产品,其定价模型是一个含有混合导数项的二维Black-Scholes方程,研究它的数值解法有着非常重要的理论意义和实际价值.本文给出求解双币种期权定价模型的基于Craig-Sneyd分裂法的一个新ADI差分方法(C-S ADI),该方法首先将二维B1ack-Scholes方程分裂为两个一维方程和一个含有混合导数的二维方程,然后分别对一维方程构造半隐式格式,对含混合导数的二维方程构造显式格式进行计算.C-S ADI差分方法具有以下优点:并行性,无条件稳定性,收敛性及空间二阶、时间一阶的计算精度.理论分析与数值试验表明,相比于经典的Crank-Nicolson差分格式和已有的基于Douglas Rachford分裂法的ADI差分格式(D-R ADI),本文格式计算精度更高,并且由于其具有天然的并行特性,本文格式比串行的Crank-Nicolson格式节省了近1/5的计算时间,证实了该方法对求解双币种期权定价模型是有效的.     

求解二维Euler方程的旋转通量混合格式

为提高求解二维Euler方程数值结果的分辨率,提出了一种旋转通量混合格式.该算法采用旋转通量法的类一维处理思想,通量函数选用满足热力学第二定律的熵稳定数值通量和具有良好鲁棒性的HLL数值通量耦合的混合格式,时间方向采用三阶强稳定Runge-Kutta方法进行推进.该旋转通量混合格式具有结构简单、分辨率高的优点,数值结果表明了该算法的良好特性.     

一维Euler方程的特征有限体积格式

提出了一种用于求解一维标量方程和无粘Euler方程组的高阶有限体积格式.其中时间离散采用Sjanpson数值积分公式从而实现时间上的高阶.利用特征线理论得到网格节点在各个时间层沿着特征线的位置,而积分公式中的节点值通过三阶和五阶的中心加权本质无震荡重构得到.最后,给出了几个数值算例验证此方法的高精度和收敛性以及捕获激波的能力.     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

A well-balanced approximate Riemann solver for compressible flows in variable cross-section ducts

A well-balanced approximate Riemann solver is introduced in this paper in order to compute approximations of one-dimensional Euler equations in variable cross-section ducts. The interface Riemann solver is grounded on the VFRoe-ncv scheme, and it enforces the preservation of Riemann invariants of the steady wave. The main properties of the scheme are detailed. We provide numerical results to assess the validity of the scheme, even when the cross-section is discontinuous. A first series is devoted to analytical test cases, and the last results correspond to the simulation of a bubble collapse.     

A Finite Volume Scheme for Savage-Hutter Equations on Unstructured Grids

A Godunov-type finite volume scheme on unstructured grids is proposed to numerically solve the Savage-Hutter equations in curvilinear coordinate. We show the direct observation that the model isn't a Galilean invariant system. At the cell boundary, the modified Harten-Lax-van Leer (HLL) approximate Riemann solver is adopted to calculate the numerical flux. The modified HLL flux is not troubled by the lack of Galilean invariance of the model and it is helpful to handle discontinuities at free interface. Rigidly the system is not always a hyperbolic system due to the dependence of flux on the velocity gradient. Even so, our numerical results still show quite good agreements with reference solutions. The simulations for granular avalanche flows with shock waves indicate that the scheme is applicable.     

Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems

We propose a new well-balanced unstaggered central finite volume scheme for hyperbolic balance laws with geometrical source terms. In particular we construct a new one and two-dimensional finite volume method for the numerical solution of shallow water equations on flat/variable bottom topographies. The proposed scheme evolves a non-oscillatory numerical solution on a single grid, avoids the time consuming process of solving Riemann problems arising at the cell interfaces, and is second-order accurate both in space and time. Furthermore, the numerical scheme follows a well-balanced discretization that first discretizes the geometrical source term according to the discretization of the flux terms, and then mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The resulting scheme exactly satisfies the C-property at the discrete level. The proposed scheme is then applied and classical one and two-dimensional shallow water equation problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

Well-balanced central finite volume methods for the Ripa system

We propose a new well-balanced central finite volume scheme for the Ripa system both in one and two space dimensions. The Ripa system is a nonhomogeneous hyperbolic system with a non-zero source term that is obtained from the shallow water equations system by incorporating horizontal temperature gradients. The proposed numerical scheme is a second-order accurate finite volume method that evolves a non-oscillatory numerical solution on a single grid, avoids the process of solving Riemann problems arising at the cell interfaces, and follows a well-balanced discretization that ensures the steady state requirement by discretizing the geometrical source term according to the discretization of the flux terms. Furthermore the proposed scheme mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The proposed scheme is then applied and classical one and two-dimensional Ripa problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

A flux-limiter method for dam-break flows over erodible sediment beds

Finite volume methods for dam-break flows over erodible sediment beds require a monotone numerical flux. In the present study we present a new flux-limiter scheme based on the Lax–Wendroff method coupled with a non-homogeneous Riemann solver and a flux limiter function. The non-homogeneous Riemann solver consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The proposed method satisfy the conservation property such that the discretization of the flux gradients and the source terms are well-balanced in the numerical solution of suspended sediment models. The flux-limiter method provides accurate results avoiding numerical oscillations and numerical dissipation in the approximated solutions. Several standard test examples are considered to verify the performance and the accuracy of the proposed method.     

Asymptotic preserving HLL schemes

This work concerns the derivation of HLL schemes to approximate the solutions of systems of conservation laws supplemented by source terms. Such a system contains many models such as the Euler equations with high friction or the M1 model for radiative transfer. The main difficulty arising from these models comes from a particular asymptotic behavior. Indeed, in the limit of some suitable parameter, the system tends to a diffusion equation. This article is devoted to derive HLL methods able to approximate the associated transport regime but also to restore the suitable asymptotic diffusive regime. To access such an issue, a free parameter is introduced into the source term. This free parameter will be a useful correction to satisfy the expected diffusion equation at the discrete level. The derivation of the HLL scheme for hyperbolic systems with source terms comes from a modification of the HLL scheme for the associated homogeneous hyperbolic system. The resulting numerical procedure is robust as the source term discretization preserves the physical admissible states. The scheme is applied to several models of physical interest. The numerical asymptotic behavior is analyzed and an asymptotic preserving property is systematically exhibited. The scheme is illustrated with numerical experiments. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1396–1422, 2011     

THE EARLY INFLUENCE OF PETER LAX ON COMPUTATIONAL HYDRODYNAMICS AND AN APPLICATION OF LAX-FRIEDRICHS AND LAX-WENDROFF ON TRIANGULAR GRIDS IN LAGRANGIAN COORDINATES

We give a brief discussion of some of the contributions of Peter Lax to Computational Fluid Dynamics. These include the Lax-Friedrichs and Lax-Wendroff numerical schemes. We also mention his collaboration in the 1983 HLL Riemann solver. We develop two-dimensional Lax-Friedrichs and Lax-Wendroff schemes for the Lagrangian form of the Euler equations on triangular grids. We apply a composite scheme that uses a Lax-Friedrichs time step as a dissipative filter after several Lax-Wendroff time steps. Numerical results for Noh’s infinite strength shock problem, the Sedov blast wave problem, and the Saltzman piston problem are presented.     

MUSTA: A multi-stage numerical flux

Numerical methods for conservation laws constructed in the framework of finite volume and discontinuous Galerkin finite elements require, as the building block, a monotone numerical flux. In this paper we present some preliminary results on the MUSTA approach [E.F. Toro, Multi-stage predictor–corrector fluxes for hyperbolic equations, Technical Report NI03037-NPA, Isaac Newton Institute for Mathematical Sciences, University of Cambridge, UK, 17th June, 2003] for constructing upwind numerical fluxes. The scheme may be interpreted as an un-conventional approximate Riemann solver that has simplicity and generality as its main features. When used in its first-order mode we observe that the scheme achieves the accuracy of the Godunov method used in conjunction with the exact Riemann solver, which is the reference first-order method for hyperbolic systems. At least for the scalar model hyperbolic equation, the Godunov scheme is the best of all first-order monote schemes, it has the smallest truncation error. Extensions of the scheme of this paper are realized in the framework of existing approaches. Here we present a second-order TVD (TVD for the scalar case) extension and show numerical results for the two-dimensional Euler equations on non-Cartesian geometries. The schemes find their best justification when solving very complex systems for which the solution of the Riemann problem, in the classical sense, is too complex, too costly or is simply unavailable.     

An Efficient Method for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations

We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography,and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results in a well-balanced scheme that can capture the steady state solution with a remarkable accuracy. This method approximates the source terms via the numerical fluxes produced by an (approximate) Riemann solver for the homogeneous hyperbolic systems with slight additional computation complexity using Newton‘s iterations and numerical integrations. This method solves well the subor super-critical flows, and with a transonic fix, also handles well the transonic flows over the concentration. Numerical examples provide strong evidence on the effectiveness of this new method for both unsteady and steady state calculations.