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     

Approximation of time‐dependent,viscoelastic fluid flow: Crank‐Nicolson,finite element approximation

In this article we analyze a fully discrete approximation to the time dependent viscoelasticity equations with an Oldroyd B constitutive equation in ? , = 2, 3. We use a Crank‐Nicolson discretization for the time derivatives. At each time level a linear system of equations is solved. To resolve the nonlinearities we use a three‐step extrapolation for the prediction of the velocity and stress at the new time level. The approximation is stabilized by using a discontinuous Galerkin approximation for the constitutive equation. For the mesh parameter, h, and the temporal step size, Δt, sufficiently small and satisfying Δt ≤ Ch , existence of the approximate solution is proven. A priori error estimates for the approximation in terms of Δt and h are also derived. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 248–283, 2004     

Local well‐posedness for the ideal incompressible density‐dependent viscoelastic flow

In this paper, we study the ideal incompressible density‐dependent Oldroyd model in . We establish local in time existence and uniqueness of solutions for the ideal incompressible density‐dependent Oldroyd model. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite volume element approximation of the coupled continuum pipe‐flow/Darcy model for flows in karst aquifers

A finite volume element method is applied to approximate the continuum pipe‐flow/Darcy problem, which models the coupled conduit flow and porous media flow in Karst aquifers. A decoupled scheme is proposed for solving the coupled discretization problem. Optimal error estimates in L² norm and H¹ norm are given in this article. Some numerical examples are presented to verify the theoretical results and demonstrate the effectiveness of the decoupling approach. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 376–392, 2014     

Stability of a colocated finite volume scheme for the Navier‐Stokes equations

The aim of this article is to describe a colocated finite volume approximation of the incompressible Navier‐Stokes equation and study its stability. One of the advantages of colocated finite volume space discretizations over staggered space discretizations is that all the variables share the same location; hence, the possibility to more easily use complex geometries and hierarchical decompositions of the unknowns. The time discretization used in the scheme studied here is a projection method. First, we give the full discretization of the incompressible Navier‐Stokes equations, then, we state the stability result and prove it following the methods of Marion and Temam. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Maximum principle for the finite element solution of time‐dependent anisotropic diffusion problems

Preservation of the maximum principle is studied for the combination of the linear finite element method in space and the θ ‐method in time for solving time‐dependent anisotropic diffusion problems. It is shown that the numerical solution satisfies a discrete maximum principle when all element angles of the mesh measured in the metric specified by the inverse of the diffusion matrix are nonobtuse, and the time step size is bounded below and above by bounds proportional essentially to the square of the maximal element diameter. The lower bound requirement can be removed when a lumped mass matrix is used. In two dimensions, the mesh and time step conditions can be replaced by weaker Delaunay‐type conditions. Numerical results are presented to verify the theoretical findings. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Finite‐time stability and boundedness of switched nonlinear time‐delay systems under state‐dependent switching

This article investigates the finite‐time stability, stabilization, and boundedness problems for switched nonlinear systems with time‐delay. Unlike the existing average dwell‐time technique based on time‐dependent switching strategy, largest region function strategy, that is, state‐dependent switching control strategy is adopted to design the switching signal, which does not require the switching instants to be given in advance. Some sufficient conditions which guarantee finite‐time stable, stabilization, and boundedness of switched nonlinear systems with time‐delay are presented in terms of linear matrix inequalities. Detail proofs are given using multiple Lyapunov‐like functions. A numerical example is given to illustrate the effectiveness of the proposed methods. © 2014 Wiley Periodicals, Inc. Complexity 21: 267–275, 2015     

A space‐time mixed‐hybrid finite element method for the damped wave equation

A space‐time finite element method is introduced to solve the linear damped wave equation. The scheme is constructed in the framework of the mixed‐hybrid finite element methods, and where an original conforming approximation of H(div;Ω) is used, the latter permits us to obtain an upwind scheme in time. We establish the link between the nonstandard finite difference scheme recently introduced by Mickens and Jordan and the scheme proposed. In this regard, two approaches are considered and in particular we employ a formulation allowing the solution to be marched in time, i.e., one only needs to consider one time increment at a time. Numerical results are presented and compared with the analytical solution illustrating good performance of the present method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Nonconforming finite element methods with subgrid viscosity applied to advection‐diffusion‐reaction equations

A nonconforming (Crouzeix–Raviart) finite element method with subgrid viscosity is analyzed to approximate advection‐diffusion‐reaction equations. The error estimates are quasi‐optimal in the sense that keeping the Péclet number fixed, the estimates are suboptimal of order in the mesh size for the L²‐norm and optimal for the advective derivative on quasi‐uniform meshes. The method is also reformulated as a finite volume box scheme providing a reconstruction formula for the diffusive flux with local conservation properties. Numerical results are presented to illustrate the error analysis. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A new combined finite element‐upwind finite volume method for convection‐dominated diffusion problems

In this article, we develop a combined finite element‐weighted upwind finite volume method for convection‐dominated diffusion problems in two dimensions, which discretizes the diffusion term with the standard finite element scheme, and the convection and source terms with the weighted upwind finite volume scheme. The developed method leads to a totally new scheme for convection‐dominated problems, which overcomes numerical oscillation, avoids numerical dispersion, and has high‐order accuracy. Stability analyses of the scheme are given for the problems with constant coefficients. Numerical experiments are presented to illustrate the stability and optimal convergence of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 799–818, 2016     

Shape reconstruction for the time‐dependent Navier‐Stokes flow

This article deals with the shape reconstruction of a bounded domain with a viscous incompressible fluid driven by the time‐dependent Navier‐Stokes equations. For the approximate solution of the ill‐posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of the domain derivative. Numerical examples indicate the feasibility of our method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

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     

An efficient splitting technique for two‐layer shallow‐water model

We consider the numerical approximation of the weak solutions of the two‐layer shallow‐water equations. The model under consideration is made of two usual one‐layer shallow‐water model coupled by nonconservative products. Because of the nonconservative products of the system, which couple both one‐layer shallow‐water subsystems, the usual numerical methods have to consider the full model. Of course, uncoupled numerical techniques, just involving finite volume schemes for the basic shallow‐water equations, are very attractive since they are very easy to implement and they are costless. Recently, a stable layer splitting technique was introduced [Bouchut and Morales de Luna, M2AN Math Model Numer Anal 42 (2008), 683–698]. In the same spirit, we exhibit new splitting technique, which is proved to be well balanced and non‐negative preserving. The main benefit issuing from the here derived uncoupled method is the ability to correctly approximate the solution of very severe benchmarks. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1396–1423, 2015     

Ideal density‐dependent incompressible viscoelastic flow in the critical Besov spaces

In this paper, we consider the well‐posedness issue for the density‐dependent incompressible viscoelastic fluids of the Oldroyd model for the ideal case in space dimension greater than 2. We obtain the local well‐posedness of this model under the assumption that the initial density is bounded away from zero in the critical Besov spaces by means of the Littlewood‐Paley theory and Bony's paradifferential calculus. In particular, we obtain a Beale‐Kato‐Majida–type regularity criterion.     

Minimum time‐step criteria for the Galerkin finite element methods applied to one‐dimensional parabolic partial differential equations

The finite element method has been well established for numerically solving parabolic partial differential equations (PDEs). Also it is well known that a too large time step should not be chosen in order to obtain a stable and accurate numerical solution. In this article, accuracy analysis shows that a too small time step should not be chosen either for some time‐stepping schemes. Otherwise, the accuracy of the numerical solution cannot be improved or can even be worsened in some cases. Furthermore, the so‐called minimum time step criteria are established for the Crank‐Nicolson scheme, the Galerkin‐time scheme, and the backward‐difference scheme used in the temporal discretization. For the forward‐difference scheme, no minimum time step exists as far as the accuracy is concerned. In the accuracy analysis, no specific initial and boundary conditions are invoked so that such established criteria can be applied to the parabolic PDEs subject to any initial and boundary conditions. These minimum time step criteria are verified in a series of numerical experiments for a one‐dimensional transient field problem with a known analytical solution. The minimum time step criteria developed in this study are useful for choosing appropriate time steps in numerical simulations of practical engineering problems. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Stability of a combined finite element ‐ finite volume discretization of convection‐diffusion equations

We consider a time‐dependent and a stationary convection‐diffusion equation. These equations are approximated by a combined finite element – finite volume method: the diffusion term is discretized by Crouzeix‐Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the nonstationary case, we use an implicit Euler approach for time discretization. This scheme is shown to be L²‐stable uniformly with respect to the diffusion coefficient. In addition, it turns out that stability is unconditional in the time‐dependent case. These results hold if the underlying grid satisfies a condition that is fulfilled, for example, by some structured meshes. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 402–424, 2012     

Adequacy of finite difference schemes for convection‐diffusion equations

Finite difference schemes for the numerical solution of singularly perturbed convection problems on uniform grids are studied in the limit case where the viscosity and the meshsize approach zero at the same time. The present error estimates are given in terms of order of magnitude in the above limit process and are useful in a priori choosing adequate schemes and meshsizes for boundary‐layer problems and problems with closed characteristics. Published 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 280–295, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10007     

Finite element approximations with quadrature for second‐order hyperbolic equations

In this article, the effect of numerical quadrature on the finite element Galerkin approximations to the solution of hyperbolic equations has been studied. Both semidiscrete and fully discrete schemes are analyzed and optimal estimates are derived in the L^∞(H¹), L^∞(L²) norms, whereas quasi‐optimal estimate is derived in the L^∞(L^∞) norm using energy methods. The analysis in the present paper improves upon the earlier results of Baker and Dougalis [SIAM J Numer Anal 13 (1976), pp 577–598] under the minimum smoothness assumptions of Rauch [SIAM J Numer Anal 22 (1985), pp 245–249] for a purely second‐order hyperbolic equation with quadrature. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 537–559, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10022     

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