A cell edge‐based linearity‐preserving scheme for diffusion problems on two‐dimensional unstructured grids

We propose a new finite volume scheme for 2D anisotropic diffusion problems on general unstructured meshes. The main feature lies in the introduction of two auxiliary unknowns on each cell edge, and then the scheme has both cell‐centered primary unknowns and cell edge‐based auxiliary unknowns. The auxiliary unknowns are interpolated by the multipoint flux approximation technique, which reduces the scheme to a completely cell‐centered one. The derivation of the scheme satisfies the linearity‐preserving criterion that requires that a discretization scheme should be exact on linear solutions. The resulting new scheme is then called as a cell edge‐based linearity‐preserving scheme. The optimal convergence rates are numerically obtained on unstructured grids in case that the diffusion tensor is taken to be anisotropic and/or discontinuous. Copyright © 2017 John Wiley & Sons, Ltd.     

A positivity‐preserving nonstandard finite difference scheme for the damped wave equation

A positivity‐preserving nonstandard finite difference scheme is constructed to solve an initial‐boundary value problem involving heat transfer described by the Maxwell‐Cattaneo thermal conduction law, i.e., a modified form of the classical Fourier flux relation. The resulting heat transport equation is the damped wave equation, a PDE of hyperbolic type. In addition, exact analytical solutions are given, special cases are mentioned, and it is noted that the positivity condition is equivalent to the usual linear stability criteria. Finally, solution profiles are plotted and possible extensions to a delayed diffusion equation and nonlinear reaction‐diffusion systems are discussed. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

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     

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.     

Monotone combined edge finite volume–finite element scheme for Anisotropic Keller–Segel model

In this article, a new numerical scheme for a degenerate Keller–Segel model with heterogeneous anisotropic tensors is treated. It is well‐known that standard finite volume scheme not permit to handle anisotropic diffusion without any restrictions on meshes. Therefore, a combined finite volume‐nonconforming finite element scheme is introduced, developed, and studied. The unknowns of this scheme are the values at the center of cell edges. Convergence of the approximate solution to the continuous solution is proved only supposing the shape regularity condition for the primal mesh. This scheme ensures the validity of the discrete maximum principle under the classical condition that all transmissibilities coefficients are positive. Therefore, a nonlinear technique is presented, as a correction of the diffusive flux, to provide a monotone scheme for general tensors. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1030–1065, 2014     

A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Post processing of solution and flux for the nodal mimetic finite difference method

We develop and analyze a post processing technique for the family of low‐order mimetic discretizations based on vertex unknowns for the numerical treatment of diffusion problems on unstructured polygonal and polyhedral meshes. The post processing works in two steps. First, from the nodal degrees of freedom, we reconstruct an elemental‐based vector field that approximates the gradient of the exact solution. Second, we solve a local problem for each mesh vertex associated with a scheme degree of freedom to determine a post processed normal flux that is conservative and divergence preserving. Theoretical results and numerical experiments for two‐dimensional (2D) and 3D benchmark problems show optimal convergence rates. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 336–363, 2015     

A Vertex-Centered and Positivity-Preserving Finite Volume Scheme for Two-Dimensional Three-Temperature Radiation Diffusion Equations on General Polygonal Meshes

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multimaterial system and explain the exchange of energy among electrons, ions and photons. In this paper, we suggest a new positivity-preserving finite volume scheme for 2-D 3-T radiation diffusion equations on general polygonal meshes. The vertex unknowns are treated as primary ones for which the finite volume equations are constructed. The edge-midpoint and cell-centered unknowns are used as auxiliary ones and interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. By comparison, most existing positivity-preserving finite volume schemes are cell-centered and based on the convex decomposition of the co-normal. Here, the co-normal decomposition is not convex in general, leading to a fixed stencil of the flux approximation and avoiding a certain search algorithm on complex grids. Moreover, the new scheme effectively alleviates the numerical heat-barrier issue suffered by most existing cell-centered or hybrid schemes in solving strongly nonlinear radiation diffusion equations. Numerical experiments demonstrate the second-order accuracy and the positivity of the solution on various distorted grids. For the problem without analytic solution, the contours of the numerical solutions obtained by our scheme on distorted meshes accord with those on smooth quadrilateral meshes.     

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

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

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

A parallel finite volume scheme preserving positivity for diffusion equation on distorted meshes

Parallel domain decomposition methods are natural and efficient for solving the implicity schemes of diffusion equations on massive parallel computer systems. A finite volume scheme preserving positivity is essential for getting accurate numerical solutions of diffusion equations and ensuring the numerical solutions with physical meaning. We call their combination as a parallel finite volume scheme preserving positivity, and construct such a scheme for diffusion equation on distorted meshes. The basic procedure of constructing the parallel finite volume scheme is based on the domain decomposition method with the prediction‐correction technique at the interface of subdomains: First, we predict the values on each inner interface of subdomains partitioned by the domain decomposition. Second, we compute the values in each subdomain using a finite volume scheme preserving positivity. Third, we correct the values on each inner interface using the finite volume scheme preserving positivity. The resulting scheme has intrinsic parallelism, and needs only local communication among neighboring processors. Numerical results are presented to show the performance of our schemes, such as accuracy, stability, positivity, and parallel speedup.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2159–2178, 2017     

High‐order asymptotic‐preserving schemes for linear systems: Application to the Goldstein–Taylor equations

In this work, we address the numerical approximation of linear systems with possibly stiff source terms which induce an asymptotic diffusion limit. More precisely, we are interested in the design of high‐order asymptotic‐preserving schemes. Our approach is based on a very simple modification of the numerical flux associated with the usual HLL scheme. This alteration can be understood as a numerical diffusion reduction technique and allows to capture the correct asymptotic behavior in the diffusion limit and to consider uniformly high‐order extensions. We more specifically consider the case of the Goldstein–Taylor model but the overall approach is shown to be easily adapted to more general systems.     

Convergence of an implicit Voronoi finite volume method for reaction–diffusion problems

We investigate the convergence of an implicit Voronoi finite volume method for reaction–diffusion problems including nonlinear diffusion in two space dimensions. The model allows to handle heterogeneous materials and uses the chemical activities of the involved species as primary variables. The numerical scheme works with boundary conforming Delaunay meshes and preserves positivity and the dissipative property of the continuous system. Starting from a result on the global stability of the scheme (uniform, mesh‐independent global upper, and lower bounds), we prove strong convergence of the chemical activities and their gradients to a weak solution of the continuous problem. To illustrate the preservation of qualitative properties by the numerical scheme, we present a long‐term simulation of the Michaelis–Menten–Henri system. Especially, we investigate the decay properties of the relative free energy over several magnitudes of time, and obtain experimental orders of convergence for this quantity. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 141–174, 2016     

Analysis of the nonlinear scheme preserving the maximum principle for the anisotropic diffusion equation on distorted meshes

Science China Mathematics - In this paper, a nonlinear finite volume scheme preserving the discrete maximum principle for the anisotropic diffusion equation on distorted meshes is described. We...     

Efficient energy‐preserving scheme of the three‐coupled nonlinear Schrödinger equation

An energy‐preserving scheme is proposed for the three‐coupled nonlinear Schrödinger (T‐CNLS) equation. The T‐CNLS equation is rewritten into the classical Hamiltonian form. Then the spatial variable is discretized by using high‐order compact method to convert it into a finite‐dimensional Hamiltonian system. Next, a second‐order averaged vector field (AVF) method is employed in time which results in an energy‐preserving scheme. Some theoretical results such as convergence are investigated. In addition, it provides some numerical examples to illustrate the robustness and reliability of the theoretical results. It also explores the role of the parameters in the model and initial condition on the wave propagation.     

A finite volume method on general meshes for a degenerate parabolic convection–reaction–diffusion equation

We propose a finite volume method on general meshes for the discretization of a degenerate parabolic convection–reaction–diffusion equation. Equations of this type arise in many contexts, such as for example the modeling of contaminant transport in porous media. The diffusion term, which can be anisotropic and heterogeneous, is discretized using a recently developed hybrid mimetic mixed framework. We construct a family of discretizations for the convection term, which uses the hybrid interface unknowns. We consider a wide range of unstructured possibly nonmatching polyhedral meshes in arbitrary space dimension. The scheme is fully implicit in time, it is locally conservative and robust with respect to the Péclet number. We obtain a convergence result based upon a priori estimates and the Fréchet–Kolmogorov compactness theorem. We implement the scheme both in two and three space dimensions and compare the numerical results obtained with the upwind and the centered discretizations of the convection term numerically.     

Block‐centered upwind multistep difference method and convergence analysis for numerical simulation of oil reservoir

The three‐dimensional displacement of two‐phase flow in porous media is a preliminary problem of numerical simulation of energy science and mathematics. The mathematical model is formulated by a nonlinear system of partial differential equations to describe incompressible miscible case. The pressure is defined by an elliptic equation, and the concentration is defined by a convection‐dominated diffusion equation. The pressure generates Darcy velocity and controls the dynamic change of concentration. We adopt a conservative block‐centered scheme to approximate the pressure and Darcy velocity, and the accuracy of Darcy velocity is improved one order. We use a block‐centered upwind multistep method to solve the concentration, where the time derivative is approximated by multistep method, and the diffusion term and convection term are treated by a block‐centered scheme and an upwind scheme, respectively. The composite algorithm is effective to solve such a convection‐dominated problem, since numerical oscillation and dispersion are avoided and computational accuracy is improved. Block‐centered method is conservative, and the concentration and the adjoint function are computed simultaneously. This physical nature is important in numerical simulation of seepage fluid. Using the convergence theory and techniques of priori estimates, we derive optimal estimate error. Numerical experiments and data show the support and consistency of theoretical result. The argument in the present paper shows a powerful tool to solve the well‐known model problem.     

Two positivity preserving flux limited,second‐order numerical methods for a haptotaxis model

Two numerical methods for a one‐dimensional haptotaxis model, which exploit the use of van Leer flux limiter, are developed and analyzed. Sufficient conditions time step size and flux limiting are given for such formulation to ensure the non‐negativity of the discrete solution and second‐order accuracy in space. Another advantage is that we avoid solving large nonlinear systems of algebraic equations. The discrete preservation of total conservation of cell density, concentration, and logarithmic density is also verified for the numerical solution. Numerical results concerning accuracy, convergence rate, positivity, and conservation properties are presented and discussed. Similar approach could be applied efficiently in the corresponding two‐ and three‐dimensional problems. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Virtual element methods for parabolic problems on polygonal meshes

The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time‐dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2110–2134, 2015     

A virtual element method-based positivity-preserving conservative scheme for convection–diffusion problems on polygonal meshes

In this paper, we propose a positivity-preserving conservative scheme based on the virtual element method (VEM) to solve convection–diffusion problems on general meshes. As an extension of finite element methods to general polygonal elements, the VEM has many advantages such as substantial mathematical foundations, simplicity in implementation. However, it is neither positivity-preserving nor locally conservative. The purpose of this article is to develop a new scheme, which has the same accuracy as the VEM and preserves the positivity of the numerical solution and local conservation on primary grids. The first step is to calculate the cell-vertex values by the lowest-order VEM. Then, the nonlinear two-point flux approximations are utilized to obtain the nonnegativity of cell-centered values and the local conservation property. The new scheme inherits both advantages of the VEM and the nonlinear two-point flux approximations. Numerical results show that the new scheme can reach the optimal convergence order of the virtual element theory, that is, the second-order accuracy for the solution and the first-order accuracy for its gradient. Moreover, the obtained cell-centered values are nonnegative, which demonstrates the positivity-preserving property of our new scheme.