非定常自由面流激波解的二阶守恒算法

将计算双曲型守恒律弱解的Lax-Wendroff型TVD格式推广到断面形状沿程任意变化的一般浅水方程组，构造了二阶精度的差分格式．新格式适用于模拟天然河道中溃坝洪水波的传播．提供了表明方法性能的算例，实际天然梯级水库溃坝问题的数值实验表明格式稳定，适应性强．     

求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

长水波近似方程组的动力学行为和辛几何算法

长水波近似方程组作为一类重要的非线性方程有着许多广泛的应用前景,特别是在浅水非线性色散波的研究中具有重要意义.给出了长水波近似方程组的动力学行为,并基于Hamilton空间体系的多辛理论研究了长水波近似方程组的数值解法,讨论了利用Preissmann方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,格式满足多辛守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

浅水长波近似方程组的拟小波解

以浅水长波近似方程组为例,提出了拟小波方法求解(1 1)维非线性偏微分方程组数值解,该方程用拟小波离散格式离散空间导数,得到关于时间的常微分方程组,用四阶Runge-K utta方法离散时间导数,并将其拟小波解与解析解进行比较和验证.     

一类交错网格的Gauss型格式

本文在交错网格的情况下 ,利用 Gauss型求积公式构造了一类不需解 Riemann问题的求解一维单个双曲守恒律的二阶显式 Gauss型差分格式 ,证明了该格式在CFL条件限制下为 TVD格式 ,并证明了这类格式的收敛性 ,然后将格式推广到方程组的情形 .由于在交错网格的情况下构造的这类差分格式 ,不需要求解 Riemann问题 ,因此这类格式与诸如 Harten等的 TVD格式相比具有如下优点 :由于不需要完整的特征向量系 ,因此可用于求解弱双曲方程组 ,计算更快、编程更加简便等 .     

求解Hamilton-Jacobi方程的有限元方法

本文将Galerkin二次有限元应于Hamilton-Jacobi方程，得到了求解Hamilton-Jacobi方程的数值格式。这些格式是TVD型的，在更强的条件下，基半离散格式的数值解收敛于Hamilton-Jacobi方程的粘性解。数值结果表明这类格式具有较高分辨导数间断的能力。     

浅水流动与污染物扩散的高分辨率计算模型

将组合型TVD格式应用于守恒型的浅水方程和污染物扩散方程，建立了二者耦合求解的高分辨率有限体积计算模型。给出了溃坝水流、明渠突扩流和污染物输运计算的典型算例，并与实验数据或其它数值结果进行了比较，证实了该模型的有效性，表明它不但能处理有激波的非恒定流问题，也能较好地计算具有任意边界的一般的浅水流动和污染物扩散问题，为浅水流动和水环境模拟提供了精度高、稳定性好、普适性强的数值方法。     

基于拟Shannon小波浅水长波近似方程组的数值解

本文研究了浅水长波近似方程组初边值问题的数值解.利用小波多尺度分析和区间拟Shannon小波,对浅水长波近似方程组空间导数实施空间离散,用时间步长自适应精细积分法对其变换所的非线性常微分方程组进行求解,得到了浅水长波近似方程组的数值解,并将此方法计算的结果与其解析解进行比较和验证.     

拟线性抛物方程组具有界面外推的并行本性差分方法

构造了拟线性抛物型方程组初边值问题的一类具有界面外推的并行本性差分格式. 为给出子区域间界面上的值或者与界面相邻点处的值,给出了两类时间外推的方式, 得到了二阶精度无条件稳定的并行差分格式. 并且不作启示性假定,证明了所构造的并行差分格式的离散向量解的存在性和 唯一性. 而且在格式的离散向量解对原始问题的已知离散数据连续依赖的意义下, 证明了并行差分格式的解按离散W^(2,1)₂(Q_Δ)范数是无条件稳定的.最后证明了具有界面外推的并行本性差分格式的离散向量解收敛到原始拟线性抛物问题的唯一广义解. 给出了数值例子,数值结果表明所构造的格式是无条件稳定的, 具有二阶精度,且具有高度并行性.     

多步Runge-Kutta型TVB时间离散

近年来TVD,TVB和ENO方法出现并得到广泛应用,见[1]—[8].特别,在[6]—[8]中利用线方法和时间离散的结合构造了TVD,TVB和 ENO差分格式.整个构造过程较Harten的工作简化得多,从而开辟了一条构造高精度无振荡差分格式的新途径.他在[6],[8]中讨论了线性多步TVB时间离散,在[7]中又讨论了Runge-Kutta型TVD时间离散,并得到了时间离散在TVD,TVB意义下所应满足的条件.本     

Implicit numerical method for computing three-dimensional rarefied gas flows on unstructured meshes

The method based on the numerical solution of a model kinetic equation is proposed for analyzing three-dimensional rarefied gas flows. The basic idea behind the method is the use of a second-order accurate TVD scheme on hybrid unstructured meshes in physical space and a fast implicit time discretization method without iterations at the upper level. The performance of the method is illustrated by computing test examples of three-dimensional rarefied gas flows in variously shaped channels in a wide range of Knudsen numbers.     

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.     

Numerical modelling of free-surface shallow flows over irregular topography with complex geometry

A well-balanced Godunov-type finite volume algorithm is developed for modelling free-surface shallow flows over irregular topography with complex geometry. The algorithm is based on a new formulation of the classical shallow water equations in hyperbolic conservation form. Unstructured triangular grids are used to achieve the adaptability of the grid to the geometry of the problem and to facilitate localised refinement. The numerical fluxes are calculated using HLLC approximate Riemann solver, and the MUSCL-Hancock predictor–corrector scheme is adopted to achieve the second-order accuracy both in space and in time where the solutions are continuous, and to achieve high-resolution results where the solutions are discontinuous. The novelties of the algorithm include preserving well-balanced property without any additional correction terms and the wet/dry front treatments. The good performance of the algorithm is demonstrated by comparing numerical and theoretical results of several benchmark problems, including the preservation of still water over a two-dimensional hump, the idealised dam-break flow over a frictionless flat rectangular channel, the circular dam-break, and the shock wave from oblique wall. Besides, two laboratory dam-break cases are used for model validation. Furthermore, a practical application related to dam-break flood wave propagation over highly irregular topography with complex geometry is presented. The results show that the algorithm can correctly account for free-surface shallow flows with respect to its effectiveness and robustness thus has bright application prospects.     

Numerical solution of laminar incompressible generalized Newtonian fluids flow

This paper deals with the numerical solution of laminar viscous incompressible flows for generalized Newtonian fluids in the branching channel. The generalized Newtonian fluids contain Newtonian fluids, shear thickening and shear thinning non-Newtonian fluids. The mathematical model is the generalized system of Navier-Stokes equations. The finite volume method combined with an artificial compressibility method is used for spatial discretization. For time discretization the explicit multistage Runge-Kutta numerical scheme is considered. Steady state solution is achieved for t → ∞ using steady boundary conditions and followed by steady residual behavior. For unsteady solution a dual-time stepping method is considered. Numerical results for flows in two dimensional and three dimensional branching channel are presented.     

TVD scheme for computing open channel wave flows

For the shallow water equations in the first approximation (Saint-Venant equations), a TVD scheme is developed for shock-capturing computations of open channel flows with discontinuous waves. The scheme is based on a special nondivergence approximation of the total momentum equation that does not involve integrals related to the cross-section pressure force and the channel wall reaction. In standard divergence difference schemes, most of the CPU time is spent on the computation of these integrals. Test computations demonstrate that the discontinuity relations reproduced by the scheme are accurate enough for actual discontinuous wave propagation to be numerically simulated. All the qualitatively distinct solutions for a dam collapsing in a trapezoidal channel with a contraction in the tailwater area are constructed as an example.     

Implicit finite element schemes for stationary compressible particle-laden gas flows

The derivation of macroscopic models for particle-laden gas flows is reviewed. Semi-implicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.     

A conservative extension of the method of characteristics for 1-D shallow flows

The method of characteristics (MOC) has been used for a long time in open channels and pipes flows. It is based on non-conservative equations, and hence it cannot be used directly for solving discontinuous shallow flows. In this paper we develop a conservative version of the MOC scheme for 1-D shallow flows by imposing the conservation law at the interpolation step. The conservation property of the scheme ensures the production of an accurate shock modeling and enables the MOC scheme to simulate dam-break type flows. By using a proper interpolation function, the proposed method can also produce quite accurate low-oscillatory results. A number of challenging test cases show considerable improvement compared to the traditional non-conservative MOC scheme in the case of dam-break type and trans-critical flow simulations.     

Second-order schemes for conservation laws with discontinuous flux modelling clarifier–thickener units

Continuously operated clarifier–thickener (CT) units can be modeled by a non-linear, scalar conservation law with a flux that involves two parameters that depend discontinuously on the space variable. This paper presents two numerical schemes for the solution of this equation that have formal second-order accuracy in both the time and space variable. One of the schemes is based on standard total variation diminishing (TVD) methods, and is addressed as a simple TVD (STVD) scheme, while the other scheme, the so-called flux-TVD (FTVD) scheme, is based on the property that due to the presence of the discontinuous parameters, the flux of the solution (rather than the solution itself) has the TVD property. The FTVD property is enforced by a new nonlocal limiter algorithm. We prove that the FTVD scheme converges to a BV _t solution of the conservation law with discontinuous flux. Numerical examples for both resulting schemes are presented. They produce comparable numerical errors, while the FTVD scheme is supported by convergence analysis. The accuracy of both schemes is superior to that of the monotone first-order scheme based on the adaptation of the Engquist–Osher scheme to the discontinuous flux setting of the CT model (Bürger, Karlsen and Towers in SIAM J Appl Math 65:882–940, 2005). In the CT application there is interest in modelling sediment compressibility by an additional strongly degenerate diffusion term. Second-order schemes for this extended equation are obtained by combining either the STVD or the FTVD scheme with a Crank–Nicolson discretization of the degenerate diffusion term in a Strang-type operator splitting procedure. Numerical examples illustrate the resulting schemes.