非结构网格上一类满足局部极值原理的三阶精度有限体积方法

对二维标量双曲型守恒律方程,发展了一类满足局部极值原理的非结构网格有限体积格式.其构造思想是,以单调数值通量为基础,通过应用基于最小二乘法的二次重构和极值限制器,使数值解满足局部极值原理.为保证数值解在光滑区域达到三阶精度,该格式可结合局部光滑探测器使用.本文从理论上分析了格式的稳定性条件,数值实验验证了格式的精度和对间断的分辨能力.     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

Time‐splitting methods with charge conservation for the nonlinear Dirac equation

In this work, four numerical time‐splitting methods are proposed for the (1 + 1)‐dimensional nonlinear Dirac equation. All of these methods (or schemes) are proved to satisfy the charge conservation in the discrete level. To enhance the computation efficiency, the block Thomas algorithm is adopted. Numerical experiments are given to test the accuracy order for these schemes, to simulate numerically the binary collision including two standing waves and two moving solitons, meanwhile, the dynamic properties for the nonlinear Dirac equation are discussed. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1582–1602, 2017     

Existence of solution of a finite volume scheme preserving maximum principle for diffusion equations

In this article, a cell‐centered finite volume scheme preserving maximum principle for diffusion equations with scalar coefficients is developed. The construction of the scheme consists of three steps: at first the discrete normal flux is obtained by a linear combination of two single‐sided fluxes, then the tangential term of the normal flux is modified by using a nonlinear combination of two single‐sided tangential fluxes, finally the auxiliary unknowns in the tangential fluxes are calculated by the convex combinations of the cell‐centered unknowns. It is proved that this nonlinear scheme satisfies the discrete maximum principle (DMP). Moreover, the existence of a solution of the nonlinear scheme is proved by using the Brouwer's fixed point theorem and the bounded estimates. Numerical experiments are presented to show that the scheme not only satisfies DMP, but also obtains the second‐order accuracy and conservation.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 80–96, 2018     

非粘性土壤中溶质运移问题的守恒混合有限元法及其数值模拟

奉文建立了非粘性土壤水中溶质运移问题的守恒混合元格式,讨论了广义解和混合元解的存在唯一性,并给出了误差估计.数值模拟结果表叫,用该方法模拟溶质运移问题是合理有效的,不仅提高了通量的模拟精度,而且使计算稳定.     

非饱和水流问题的混合元法及其数值模拟

1.引 言 均质土壤中的地下水流动可归结为非饱和土壤水的流动,是土壤水未完全充满孔隙时的流动,是多孔介质流体运动的一种重要形式.非饱和流动的预报在大气科学、土壤学、农业     

Solving hyperbolic conservation laws using multiquadric quasi‐interpolation

In this article, we apply the univariate multiquadric (MQ) quasi‐interpolation to solve the hyperbolic conservation laws. At first we construct the MQ quasi‐interpolation corresponding to periodic and inflow‐outflow boundary conditions respectively. Next we obtain the numerical schemes to solve the partial differential equations, by using the derivative of the quasi‐interpolation to approximate the spatial derivative of the differential equation and a low‐order explicit difference to approximate the temporal derivative of the differential equation. Then we verify our scheme for the one‐dimensional Burgers' equation (without viscosity). We can see that the numerical results are very close to the exact solution and the computational accuracy of the scheme is ??(τ), where τ is the temporal step. We can improve the accuracy by using the high‐order quasi‐interpolation. Moreover the methods can be generalized to the other equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A linearized high‐order difference scheme for the fractional Ginzburg–Landau equation

The numerical solution for the one‐dimensional complex fractional Ginzburg–Landau equation is considered and a linearized high‐order accurate difference scheme is derived. The fractional centered difference formula, combining the compact technique, is applied to discretize fractional Laplacian, while Crank–Nicolson/leap‐frog scheme is used to deal with the temporal discretization. A rigorous analysis of the difference scheme is carried out by the discrete energy method. It is proved that the difference scheme is uniquely solvable and unconditionally convergent, in discrete maximum norm, with the convergence order of two in time and four in space, respectively. Numerical simulations are given to show the efficiency and accuracy of the scheme. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 105–124, 2017     

Antidiffusive Lagrangian‐remap schemes for models of polydisperse sedimentation

One‐dimensional models of gravity‐driven sedimentation of polydisperse suspensions with particles that belong to N size classes give rise to systems of N strongly coupled, nonlinear first‐order conservation laws for the local solids volume fractions. As the eigenvalues and eigenvectors of the flux Jacobian have no closed algebraic form, characteristic‐wise numerical schemes for these models become involved. Alternative simple schemes for this model directly utilize the velocity functions and are based on splitting the system of conservation laws into two different first‐order quasi‐linear systems, which are solved successively for each time iteration, namely, the Lagrangian and remap steps (so‐called Lagrangian‐remap [LR] schemes). This approach was advanced in (Bürger, Chalons, and Villada, SIAM J Sci Comput 35 (2013), B1341–B1368) for a multiclass Lighthill–Whitham‐Richards traffic model with nonnegative velocities. By incorporating recent antidiffusive techniques for transport equations a new version of these Lagrangian‐antidiffusive remap (L‐AR) schemes for the polydisperse sedimentation model is constructed. These L‐AR schemes are supported by a partial analysis for N = 1. They are total variation diminishing under a suitable CFL condition and therefore converge to a weak solution. Numerical examples illustrate that these schemes, including a more accurate version based on MUSCL extrapolation, are competitive in accuracy and efficiency with several existing schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1109–1136, 2016     

多边形网格上扩散方程新的单调格式

 本文在星形多边形网格上, 构造了扩散方程新的单调有限体积格式.该格式与现有的基于非线性两点流的单调格式的主要区别是, 在网格边的法向流离散模板中包含当前边上的点, 在推导离散法向流的表达式时采用了定义于当前边上的辅助未知量, 这样既可适应网格几何大变形, 同时又兼顾了当前网格边上物理量的变化. 在光滑解情形证明了离散法向流的相容性.对于具有强各向异性、非均匀张量扩散系数的扩散方程, 证明了新格式是单调的, 即格式可以保持解析解的正性. 数值结果表明在扭曲网格上, 所构造的格式是局部守恒和保正的, 对光滑解有高于一阶的精度, 并且, 针对非平衡辐射限流扩散问题, 数值结果验证了新格式在计算效率和守恒精度上优于九点格式.     

Some temporal second order difference schemes for fractional wave equations

In this article, motivated by Alikhanov's new work (Alikhanov, J Comput Phys 280 (2015), 424–438), some difference schemes are proposed for both one‐dimensional and two‐dimensional time‐fractional wave equations. The obtained schemes can achieve second‐order numerical accuracy both in time and in space. The unconditional convergence and stability of these schemes in the discrete H¹‐norm are proved by the discrete energy method. The spatial compact difference schemes with the results on the convergence and stability are also presented. In addition, the three‐dimensional problem is briefly mentioned. Numerical examples illustrate the efficiency of the proposed schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 970–1001, 2016     

An enhanced‐physics‐based scheme for the NS‐α turbulence model

We study a new enhanced‐physics‐based numerical scheme for the NS‐alpha turbulence model that conserves both energy and helicity. Although most turbulence models (in the continuous case) conserve only energy, NS‐alpha is one of only a very few that also conserve helicity. This is one reason why it is becoming accepted as the most physically accurate turbulence model. However, no numerical scheme for NS‐alpha, until now, conserved both energy and helicity, and thus the advantage gained in physical accuracy by modeling with NS‐alpha could be lost in a computation. This report presents a finite element numerical scheme, and gives a rigorous analysis of its conservation properties, stability, solution existence, and convergence. A key feature of the analysis is the identification of the discrete energy and energy dissipation norms, and proofs that these norms are equivalent (provided a careful choice of filtering radius) in the discrete space to the usual energy and energy dissipation norms. Numerical experiments are given to demonstrate the effectiveness of the scheme over usual (helicity‐ignoring) schemes. A generalization of this scheme to a family of high‐order NS‐alpha‐deconvolution models, which combine the attractive physical properties of NS‐alpha with the high accuracy gained by combining α‐filtering with van Cittert approximate deconvolution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Analysis of new conservative difference scheme for two-dimensional Rosenau-RLW equation

In this article, numerical solution for the Rosenau-RLW equation in 2D is considered and a conservative Crank–Nicolson finite difference scheme is proposed. Existence of the numerical solutions for the difference scheme has been shown by Browder fixed point theorem. A priori bound and uniqueness as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability and a second-order accuracy on both space and time of the difference scheme are proved. Numerical experiments are given to support our theoretical results.     

Locally divergence-preserving schemes of higher order

Violation of the divergence constraint on the magnetic flux density in magnetohydrodynamical (MHD) simulations leads to stability problems. It is therefore of great importance to numerically respect this intrinsic constraint. Since the divergence preservation is a local phenomenon inherent in the MHD-system it is appealing to mimic this property numerically by a locally divergence-preserving scheme. A common numerical technique for simulation of the MHD-system of conservation laws is the finite volume (FV) method. In [SISC 26 2005 pp. 1166] a local procedure to redistribute the numerical fluxes in a FV-scheme so that a discrete divergence operator vanishes was presented. This procedure stabilizes the base scheme and respects the accuracy to the second order level. The present note describes a development of the above procedure that complies with the finite volume framework, preserves a fourth order discrete divergence operator locally and retains the accuracy of a generic semi-discrete finite volume scheme up to fourth order. The redistribution of the numerical magnetic field fluxes is formulated in a standard conservative setting, making it trivial to implement the divergence-preserving modification in an existing FV-scheme; see [JCP 227 2008 pp.3405] for the details. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A LEGENDRE PSEUDOSPECTRAL METHOD FOR SOLVINGNONLINEAR KLEIN-GORDON EQUATION

1.IntroductionAsweknow,theKlein-Gordonequationisanimportantmathematicalmodelinquantummechanics.Itisoftheformwherefi=(--1,1)",x=(xl,x25'jx.),bisarealnumber.AssumethatU000=UI(x)=0onOnandAsin[1],itcanbeshownthatifUOEHI(n)nH(n),UIEL'(fl)andfEL'(o,T;L'(n)),then(1.1)hasuniquesolutionUELoo(o,T;Hi(n)nH(fl)).oftheotherhand,somefinitedifferenceschemeswereproposedin12,31withstrictproo]ofgeneralizedstabilityandconvergence.TheirnumericalsolutionskeepthediscretEconservations.Oneofspecialcases(…     

Multisymplectic Fourier pseudo-spectral integrators for Klein-Gordon-Schrödinger equations

A multisymplectic Fourier pseudo-spectral scheme,which exactly preserves the discrete multisymplectic conservation law,is presented to solve the Klein-Gordon-Schrdinger equations.The scheme is of spectral accuracy in space and of second order in time.The scheme preserves the discrete multisymplectic conservation law and the charge conservation law.Moreover,the residuals of some other conservation laws are derived for the geometric numerical integrator.Extensive numerical simulations illustrate the numerical behavior of the multisymplectic scheme,and demonstrate the correctness of the theoretical analysis.     

Entropy production in second-order three-point schemes

Summary We discuss semi-discrete three-point finite difference methods for the numerical solution of system of conservation laws which are second order accurate in space in the sense of truncation error. Particular discretizations of the numerical entropy flux associated with such schemes are studied clarifying the importance of this discretization with regard to the production of numerical entropy. Using a numerical entropy flux constructed in a canonical way we prove that a wide class of finite difference methods cannot satisfy a discrete entropy inequality. Together with a well known result of Schonbek concerning Lax-Wendroff type schemes our result indicates a strong relationship between entropy production and oscillations in numerical solutions.The research reported here was supported by a grant from the Stiftung Volkswagenwerk, Federal Republic of Germany. It is a part of the doctoral thesis of the above author, Universität Stuttgart, 1991.     

High Resolution Schemes and Discrete Entropy Conditions for 2-D Linear Conservation Laws

In this paper, fully discrete entropy conditions of a class of high resolution schemes with the MmB property are discussed by using the theory of proper discrete entropy flux for the linear scalar conservation laws in two dimensions. The theoretical resluts show that the high resolution schemes satisfying fully discrete entropy conditions with proper discrete entropy flux cannot preserve second order accuracy in the case of two dimensions.     

High-resolution schemes for bubbling flow computations

This paper is concerned with the effects of numerical schemes on the simulation of dense gas-particle two-phase flows. The first-order upwind difference, the central difference, the second-order upwind difference, the central difference plus artificial dissipation, the deferred correction, the quadratic upstream interpolation (for convective kinematics) and the monotone upstream-centered schemes for conservation laws with different flux limiters were all accomplished to simulate the multi-phase flows. It was found that numerical schemes may significantly affect the solution accuracy and numerical convergence. The monotone upstream-centered schemes for conservation laws are the best choice of all, because they can effectively suppress the non-physical oscillations with the introduction of adaptive numerical dissipation into numerical solutions.     

A Large Time Step explicit scheme (CFL>1) on unstructured grids for 2D conservation laws: Application to the homogeneous shallow water equations

In this work, a Large Time Step (LTS) explicit finite volume scheme designed to allow CFL > 1 is applied to the numerical resolution of 2D scalar and systems of conservation laws on triangular grids. Based on the flux difference splitting formulation, a special concern is put on finding the way of packing the information to compute the numerical solution when working on unstructured grids. Not only the cell areas but also the length of the interfaces and their orientation are questions of interest to send the information from each edge or interface. The information to update the cell variables is computed according to the local average discrete velocity and the orientation of the edges of the cells involved. The performance of these ideas is tested and compared with the conventional explicit first order and second order schemes in academic configurations for the 2D linear scalar equation and for 2D systems of conservation laws (in particular the shallow water equations) without source terms. The LTS scheme is demonstrated to preserve or even gain accuracy and save computational time with respect to the first order scheme.