High‐resolution monotone schemes based on quasi‐characteristics technique

In this article, we consider a new technique that allows us to overcome the well‐known restriction of Godunov's theorem. According to Godunov's theorem, a second‐order explicit monotone scheme does not exist. The techniques in the construction of high‐resolution schemes with monotone properties near the discontinuities of the solution lie in choosing of one of two high‐resolution numerical solutions computed on different stencils. The criterion for choosing the final solution is proposed. Results of numerical tests that compare with the exact solution and with the numerical solution obtained by the first‐order monotone scheme are presented. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 262–276, 2001     

On splitting‐based mass and total energy conserving arbitrary order shallow‐water schemes

A new method for numerical solution to the shallow‐water equations is suggested. The method allows constructing a family of finite difference schemes of different approximation order that conserve the mass and the total energy. Our approach is based on the method of splitting, and unlike others it permits to derive conservative numerical schemes after discretizing all the partial derivatives, both spatial and temporal. The schemes thus appear to be fully discrete, both in time and in space. Besides, due to a simple structure of the matrices appeared therewith, the method provides essential benefits in the computational cost of solution and is easy‐to‐implement in the Cartesian and spherical geometries. Numerical results confirm functionality and efficiency of the developed method. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Pressure‐based adaption indicator for compressible euler equations

We consider the Euler equations of gas dynamics and develop a new adaption indicator, which is based on the weak local residual measured for the nonconservative pressure variable. We demonstrate that the proposed indicator is capable of automatically detecting discontinuities and distinguishing between the shock and contact waves when they are isolated from each other. We then use the developed indicator to design a scheme adaption algorithm, according to which nonlinear limiters are used only in the vicinity of shocks. The new adaption algorithm is realized using a second‐order limited and a high‐order nonlimited central‐upwind scheme. We demonstrate robustness and high resolution of the designed method on a number of one‐ and two‐dimensional numerical examples. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1844–1874, 2015     

Asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation terms

We devise a new class of asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation source terms governed by a relaxation time ε. As an alternative to classical operator‐splitting techniques, the objectives of these schemes are twofold: first, to give accurate numerical solutions for large, small, and in‐between values of ε and second, to make optional the choice of the numerical scheme in the asymptotic regime ε tends to zero. The latter property may be of particular interest to make easier and more efficient the coupling at a fixed spatial interface of two models involving very different values of ε. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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.     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

An upwind finite volume scheme and its maximum‐principle‐preserving ADI splitting for unsteady‐state advection‐diffusion equations

We develop an upwind finite volume (UFV) scheme for unsteady‐state advection‐diffusion partial differential equations (PDEs) in multiple space dimensions. We apply an alternating direction implicit (ADI) splitting technique to accelerate the solution process of the numerical scheme. We investigate and analyze the reason why the conventional ADI splitting does not satisfy maximum principle in the context of advection‐diffusion PDEs. Based on the analysis, we propose a new ADI splitting of the upwind finite volume scheme, the alternating‐direction implicit, upwind finite volume (ADFV) scheme. We prove that both UFV and ADFV schemes satisfy maximum principle and are unconditionally stable. We also derive their error estimates. Numerical results are presented to observe the performance of these schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 211–226, 2003     

Indecomposable quasi‐characteristics scheme on pyramidal stencil and its application for numerical simulation of two‐phase flows through heterogeneous porous medium

A new high‐resolution indecomposable quasi‐characteristics scheme with monotone properties based on pyramidal stencil is considered. This scheme is based on consideration of two high‐resolution numerical schemes approximated governing equations on the pyramidal stencil with different kinds of dispersion terms approximation. Two numerical solutions obtained by these schemes are analyzed, and the final solution is chosen according to the special criterion to provide the monotone properties in regions where discontinuities of solutions could arise. This technique allows to construct the high‐order monotone solutions and keeps both the monotone properties and the high‐order approximation in regions with discontinuities of solutions. The selection criterion has a local character suitable for parallel computation. Application of the proposed technique to the solution of the time‐dependent 2D two‐phase flows through the porous media with the essentially heterogeneous properties is considered, and some numerical results are presented. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 44–55, 2002     

Iterative operator‐splitting methods for nonlinear differential equations and applications

In this article, we consider iterative operator‐splitting methods for nonlinear differential equations with bounded and unbounded operators. The main feature of the proposed idea is the embedding of Newton's method for solving the split parts of the nonlinear equation at each step. The convergence properties of such a mixed method are studied and demonstrated. We confirm with numerical applications the effectiveness of the proposed scheme in comparison with the standard operator‐splitting methods by providing improved results and convergence rates. We apply our results to deposition processes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1026–1054, 2011     

Conservative semi‐Lagrangian finite volume schemes

Semi‐Lagrangian finite volume schemes for the numerical approximation of linear advection equations are presented. These schemes are constructed so that the conservation properties are preserved by the numerical approximation. This is achieved using an interpolation procedure based on area‐weighting. Numerical results are presented illustrating some of the features of these schemes. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:403–425, 2001     

Monotone first‐order weighted schemes for scalar conservation laws

In this work, we present a monotone first‐order weighted (FORWE) method for scalar conservation laws using a variational formulation. We prove theoretical properties as consistency, monotonicity, and convergence of the proposed scheme for the one‐dimensional (1D) Cauchy problem. These convergence results are extended to multidimensional scalar conservation laws by a dimensional splitting technique. For the validation of the FORWE method, we consider some standard bench‐mark tests of bidimensional and 1D conservation law equations. Finally, we analyze the accuracy of the method with L¹ and L^∞ error estimates. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Splitting methods for Hamilton‐Jacobi equations

We explain how the exploitation of several kinds of operator splitting methods, both local and global in time, lead to simple numerical schemes approximating the solution of nonlinear Hamilton‐Jacobi equations. We review the existing local methods which have been used since the early 80's and we introduce a new method which is global in time. We show some numerical experiments. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

SERK2v2: A new second‐order stabilized explicit Runge‐Kutta method for stiff problems

Traditionally, explicit numerical algorithms have not been used with stiff ordinary differential equations (ODEs) due to their stability. Implicit schemes are usually very expensive when used to solve systems of ODEs with very large dimension. Stabilized Runge‐Kutta methods (also called Runge–Kutta–Chebyshev methods) were proposed to try to avoid these difficulties. The Runge–Kutta methods are explicit methods with extended stability domains, usually along the negative real axis. They can easily be applied to large problem classes with low memory demand, they do not require algebra routines or the solution of large and complicated systems of nonlinear equations, and they are especially suited for discretizations using the method of lines of two and three dimensional parabolic partial differential equations. In Martín‐Vaquero and Janssen [Comput Phys Commun 180 (2009), 1802–1810], we showed that previous codes based on stabilized Runge–Kutta algorithms have some difficulties in solving problems with very large eigenvalues and we derived a new code, SERK2, based on sixth‐order polynomials. Here, we develop a new method based on second‐order polynomials with up to 250 stages and good stability properties. These methods are efficient numerical integrators of very stiff ODEs. Numerical experiments with both smooth and nonsmooth data support the efficiency and accuracy of the new algorithms when compared to other well‐known second‐order methods such as RKC and ROCK2. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Understanding Saul'yev‐type unconditionally stable schemes from exponential splitting

Saul'yev‐type asymmetric schemes have been widely used in solving diffusion and advection equations. In this work, we show that Saul'yev‐type schemes can be derived from the exponential splitting of the semidiscretized equation which fundamentally explains their unconditional stability. Furthermore, we show that optimal schemes are obtained by forcing each scheme's amplification factor to match that of the exact amplification factor. A new second‐order explicit scheme is found for solving the advection equation with the identical amplification factor as the implicit Crank–Nicolson algorithm. Other new schemes for solving the advection–diffusion equation are also derived.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1961–1983, 2014     

SWAMP: A two‐dimensional hydrodynamic and quality modeling platform for shallow waters

Water quality two‐dimensional models are often partitioned into separate modules with separate hydraulic and biological units. In most cases this approach results in poor flexibility whenever the biological dynamics has to be adapted to a specific situation. Conversely, an integrated approach is pursued in this article, producing a two‐dimensional hydraulic‐water quality model, named Shallow Water Analysis and Modeling Program (SWAMP) designed for shallow water bodies. The major objective of the work is to create a comprehensive two‐dimensional water quality assessment tool, based on an open framework and combining easy programming of additional procedures with a user‐friendly interface. The model is based on the numerical solution of the partial differential equations describing advection‐diffusion and biological processes on a two‐dimensional rectangular finite elements mesh. The hydraulics and advection‐diffusion modules model were validated both with experimental tracer data collected at a constructed wetland site and a comparison with a commercial hydrodynamic software, showing good agreement in both cases. Moreover, the model was tested in critical conditions for mass conservation, such as time‐varying wet boundary, showing a considerable numerical robustness. In the last part of the article water quality simulations are presented, though validation data are not yet available. Nevertheless, the observed model response demonstrates general consistency with expected results and the advantages of integrating the hydraulic and quality modules. The interactive graphical user interface (GUI) is also shown to represent a simple and effective connective tool to the integrated package. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 663–687, 2002; DOI 10.1002/num.10014     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

A nonconventional Eulerian‐Lagrangian single‐node collocation method with Hermite polynomials for unsteady‐state advection‐diffusion equations

We developed a nonconventional Eulerian‐Lagrangian single‐node collocation method (ELSCM) with piecewise‐cubic Hermite polynomials as basis functions for the numerical simulation to unsteady‐state advection‐diffusion transport partial differential equations. This method greatly reduces the number of unknowns in the conventional collocation method, and generates accurate numerical solutions even if very large time steps are taken. The method is relatively easy to formulate. Numerical experiments are presented to show the strong potential of this method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 271–283, 2003.     

A finite‐volume scheme for the multidimensional quantum drift‐diffusion model for semiconductors

A finite‐volume scheme for the stationary unipolar quantum drift‐diffusion equations for semiconductors in several space dimensions is analyzed. The model consists of a fourth‐order elliptic equation for the electron density, coupled to the Poisson equation for the electrostatic potential, with mixed Dirichlet‐Neumann boundary conditions. The numerical scheme is based on a Scharfetter‐Gummel type reformulation of the equations. The existence of a sequence of solutions to the discrete problem and its numerical convergence to a solution to the continuous model are shown. Moreover, some numerical examples in two space dimensions are presented. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1483–1510, 2011     

A finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors

An implicit Euler finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors is analyzed. The model consists of strongly coupled parabolic equations for the electron density matrix or, alternatively, of weakly coupled equations for the charge and spin‐vector densities, coupled to the Poisson equation for the electric potential. The equations are solved in a bounded domain with mixed Dirichlet–Neumann boundary conditions. The charge and spin‐vector fluxes are approximated by a Scharfetter–Gummel discretization. The main features of the numerical scheme are the preservation of nonnegativity and bounds of the densities and the dissipation of the discrete free energy. The existence of a bounded discrete solution and the monotonicity of the discrete free energy are proved. For undoped semiconductor materials, the numerical scheme is unconditionally stable. The fundamental ideas are reformulations using spin‐up and spin‐down densities and certain projections of the spin‐vector density, free energy estimates, and a discrete Moser iteration. Furthermore, numerical simulations of a simple ferromagnetic‐layer field‐effect transistor in two space dimensions are presented. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 819–846, 2016     

Numerical solution of Fokker‐Planck equation using the cubic B‐spline scaling functions

In this article a numerical technique is presented for the solution of Fokker‐Planck equation. This method uses the cubic B‐spline scaling functions. The method consists of expanding the required approximate solution as the elements of cubic B‐spline scaling function. Using the operational matrix of derivative, the problem will be reduced to a set of algebraic equations. Some numerical examples are included to demonstrate the validity and applicability of the technique. The method is easy to implement and produces very accurate results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009