A finite element method for the axisymmetric three‐field Stokes system with a discontinuous pressure

A quadrilateral based velocity‐pressure‐extrastress tensor mixed finite element method for solving the three‐field Stokes system in the axisymmetric case is studied. The method derived from Fortin's Q₂ − P₁ velocity‐pressure element is to be used in connection with the standard Galerkin formulation. This makes it particularly suitable for the numerical simulation of viscoelastic flow. It is proven to be second‐order convergent in the natural weighted Sobolev norms, for the system under consideration. The crucial result that the method is uniformly stable is proven for the case of rectangular meshes. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 739–763, 1999     

Stabilized multiscale finite element method for the stationary Navier-Stokes equations

In the paper, a stabilized multiscale finite element method for the stationary incompressible Navier-Stokes equations is considered. The method is a Petrov-Galerkin approach based on the multiscale enrichment of the standard polynomial space enriched with the unusual bubble functions which no longer vanish on every element boundary for the velocity space. The stability of the P₁-P₀ triangular element (or the Q₁-P₀ quadrilateral element) is established. And the optimal error estimates of the stabilized multiscale finite element method for the stationary Navier-Stokes equations are obtained.     

Low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions

In this paper, we consider low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions. The time discretization is based on the Euler implicit scheme, and the spatial discretization is based on the low‐order element (P₁−P₁ or P₁−P₀) for the approximation of the velocity and pressure. Moreover, some error estimates for the numerical solution of fully discrete stabilized finite element scheme are obtained. Finally, numerical experiments are performed to confirm our theoretical results.     

Galerkin discontinuous approximation of the transport equation and viscoelastic fluid flow on quadrilaterals

Numerical simulation of industrial processes involving viscoelastic liquids is often based on finite element methods on quadrilateral meshes. However, numerical analysis of these methods has so far been limited to triangular meshes. In this work, we consider quadrilateral meshes. We first study the approximation of the transport equation by a Galerkin discontinuous method and prove an 𝒪(h^k+1/2) error estimates for the Q_k finite element. Then we study a differential model for viscoelastic flow with unknowns u the velocity, p the pressure, and σ the viscoelastic part of the extra-stress tensor. The approximations are ((Q₁)² transforms of) Q_k+1 continuous for u, Q_k discontinuous for σ, and P_k discontinuous for p, with k ≥ 1. Upwinding for σ is obtained by the Galerkin discontinuous method. We show that an error estimate of order 𝒪(h^k+1/2) is valid in the energy norm for the three unknowns. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 97–114, 1998     

The Crank‐Nicolson/Adams‐Bashforth scheme for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

In this article, we study the stability and convergence of the Crank‐Nicolson/Adams‐Bashforth scheme for the two‐dimensional nonstationary Navier‐Stokes equations with a nonsmooth initial data. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the implicit Crank‐Nicolson scheme for the linear terms and the explicit Adams‐Bashforth scheme for the nonlinear term. Moreover, we prove that the scheme is almost unconditionally stable for a nonsmooth initial data u₀ with div u₀ = 0, i.e., the time step τ satisfies: τ ≤ C₀ if u₀ ∈ H¹ ∩ L^∞; τ |log h| ≤ C₀ if u₀ ∈ H¹ for the mesh size h and some positive constant C₀. Finally, we obtain some error estimates for the discrete velocity and pressure under the above stability condition. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 155‐187, 2012     

Finite volume method based on stabilized finite elements for the nonstationary Navier–Stokes problem

A finite volume method based on stabilized finite element for the two‐dimensional nonstationary Navier–Stokes equations is investigated in this work. As in stabilized finite element method, macroelement condition is introduced for constructing the local stabilized formulation of the nonstationary Navier–Stokes equations. Moreover, for P₁ ? P₀ element, the H¹ error estimate of optimal order for finite volume solution (u_h,p_h) is analyzed. And, a uniform H¹ error estimate of optimal order for finite volume solution (u_h, p_h) is also obtained if the uniqueness condition is satisfied. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Mixed finite element methods for incompressible flow: Stationary Stokes equations

In this article, we develop and analyze a mixed finite element method for the Stokes equations. Our mixed method is based on the pseudostress‐velocity formulation. The pseudostress is approximated by the Raviart‐Thomas (RT) element of index k ≥ 0 and the velocity by piecewise discontinuous polynomials of degree k. It is shown that this pair of finite elements is stable and yields quasi‐optimal accuracy. The indefinite system of linear equations resulting from the discretization is decoupled by the penalty method. The penalized pseudostress system is solved by the H(div) type of multigrid method and the velocity is then calculated explicitly. Alternative preconditioning approaches that do not involve penalizing the system are also discussed. Finally, numerical experiments are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Stabilized finite element method for the viscoelastic Oldroyd fluid flows

In this paper, a new stabilized finite element method based on two local Gauss integrations is considered for the two-dimensional viscoelastic fluid motion equations, arising from the Oldroyd model for the non-Newtonian fluid flows. This new stabilized method presents attractive features such as being parameter-free, or being defined for non-edge-based data structures. It confirms that the lowest equal-order P ₁???P ₁ triangle element and Q ₁???Q ₁ quadrilateral element are compatible. Moreover, the long time stabilities and error estimates for the velocity in H ¹-norm and for the pressure in L ²-norm are obtained. Finally, some numerical experiments are performed, which show that the new method is applied to this model successfully and can save lots of computational cost compared with the standard ones.     

Numerical and theoretical study of weak Galerkin finite element solutions of Turing patterns in reaction–diffusion systems

In this paper, we introduce numerical schemes and their analysis based on weak Galerkin finite element framework for solving 2‐D reaction–diffusion systems. Weak Galerkin finite element method (WGFEM) for partial differential equations relies on the concept of weak functions and weak gradients, in which differential operators are approximated by weak forms through the Green's theorem. This method allows the use of totally discontinuous functions in the approximation space. In the current work, the WGFEM solves reaction–diffusion systems to find unknown concentrations (u, v) in element interiors and boundaries in the weak Galerkin finite element space WG(P₀, P₀, RT₀) . The WGFEM is used to approximate the spatial variables and the time discretization is made by the backward Euler method. For reaction–diffusion systems, stability analysis and error bounds for semi‐discrete and fully discrete schemes are proved. Accuracy and efficiency of the proposed method successfully tested on several numerical examples and obtained results satisfy the well‐known result that for small values of diffusion coefficient, the steady state solution converges to equilibrium point. Acquired numerical results asserted the efficiency of the proposed scheme.     

Global analysis of virus dynamics model with logistic mitosis,cure rate and delay in virus production

In this paper, we study a virus dynamics model with logistic mitosis, cure rate, and intracellular delay. By means of construction of a suitable Lyapunov functionals, obtained by linear combinations of Volterra—type functions, composite quadratic functions and Volterra—type functionals, we provide the global stability for this model. If R₀, the basic reproductive number, satisfies R₀ ≤ 1, then the infection‐free equilibrium state is globally asymptotically stable. Our system is persistent if R₀ > 1. On the other hand, if R₀ > 1, then infection‐free equilibrium becomes unstable and a unique infected equilibrium exists. The local stability analysis is carried out for the infected equilibrium, and it is shown that, if the parameters satisfy a condition, the infected equilibrium can be unstable and a Hopf bifurcation can occur. We also have that if R₀ > 1, then the infected equilibrium state is globally asymptotically stable if a sufficient condition is satisfied. We illustrate our findings with some numerical simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

On stability of the P n mod /P n element for incompressible flow problems

It is well known that finite element spaces used for approximating the velocity and the pressure in an incompressible flow problem have to be stable in the sense of the inf-sup condition of Babuška and Brezzi if a stabilization of the incompressibility constraint is not applied. In this paper we consider a recently introduced class of triangular nonconforming finite elements of nth order accuracy in the energy norm called P _n ^mod elements. For n ≤ 3 we show that the stability condition holds if the velocity space is constructed using the P _n ^mod elements and the pressure space consists of continuous piecewise polynomial functions of degree n. This research has been supported by the Grant Agency of the Czech Republic under the grant No. 201/05/0005 and by the grant MSM 0021620839.     

解Stokes问题的P1非协调四边形元的稳定化方法

本文研究了用(~P)_1-Q_0元(其中(~P)_1表示P_1非协调四边形元)解Stokes问题的非协调混合有限元稳定化逼近方法.(~P)_1-Q_0元不满足LBB条件(见[7,14] ),因而其不能直接用来求解Stokes问题.受[3] 的启发,我们提出了一种用(~P)_1-Q_0元解Stokes问题的稳定化方法,证明了这种方法的稳定性和离散问题解的存在唯一性,得到了最优误差估计.文章最后给出的数值算例验证了我们的理论结果.     

Modeling diseases with latency and nonlinear incidence rates: global dynamics of a multi‐group model

In this paper, we perform global stability analysis of a multi‐group SEIR epidemic model in which we can consider the heterogeneity of host population and the effects of latency and nonlinear incidence rates. For a simpler version that assumes an identical natural death rate for all groups, and with a gamma distribution for the latency, the basic reproduction number is defined by the theory of the next generation operator and proved to be a sharp threshold determining whether or not disease spread. Under certain assumptions, the disease‐free equilibrium is globally asymptotically stable if R₀≤1 and there exists a unique endemic equilibrium which is globally asymptotically stable if R₀>1. The proofs of global stability of equilibria exploit a matrix‐theoretic method using Perron eigenvetor, a graph‐theoretic method based on Kirchhoff's matrix tree theorem and Lyapunov functionals. Copyright © 2015 John Wiley & Sons, Ltd.     

Stability and convergence of the mark and cell finite difference scheme for Darcy‐Stokes‐Brinkman equations on non‐uniform grids

In this paper, we consider the mark and cell (MAC) method for Darcy‐Stokes‐Brinkman equations and analyze the stability and convergence of the method on nonuniform grids. Firstly, to obtain the stability for both velocity and pressure, we establish the discrete inf‐sup condition. Then we introduce an auxiliary function depending on the velocity and discretizing parameters to analyze the super‐convergence. Finally, we obtain the second‐order convergence in L2 norm for both velocity and pressure for the MAC scheme, when the perturbation parameter ? is not approaching 0. We also obtain the second‐order convergence for some terms of ∥·∥_? norm of the velocity, and the other terms of ∥·∥_? norm are second‐order convergence on uniform grid. Numerical experiments are carried out to verify the theoretical results.     

A new stabilized finite volume method for the stationary Stokes equations

In this paper we develop and study a new stabilized finite volume method for the two-dimensional Stokes equations. This method is based on a local Gauss integration technique and the conforming elements of the lowest-equal order pair (i.e., the P ₁–P ₁ pair). After a relationship between this method and a stabilized finite element method is established, an error estimate of optimal order in the H ¹-norm for velocity and an estimate in the L ²-norm for pressure are obtained. An optimal error estimate in the L ²-norm for the velocity is derived under an additional assumption on the body force. This work is supported in part by the NSF of China 10701001 and by the US National Science Foundation grant DMS-0609995 and CMG Chair Funds in Reservoir Simulation.     

A coupling of nonconforming and mixed finite element methods for Biot's consolidation model

In this article, we develop a nonconforming mixed finite element method to solve Biot's consolidation model. In particular, this work has been motivated to overcome nonphysical oscillations in the pressure variable, which is known as locking in poroelasticity. The method is based on a coupling of a nonconforming finite element method for the displacement of the solid phase with a standard mixed finite element method for the pressure and velocity of the fluid phase. The discrete Korn's inequality has been achieved by adding a jump term to the discrete variational formulation. We prove a rigorous proof of a‐priori error estimates for both semidiscrete and fully‐discrete schemes. Optimal error estimates have been derived. In particular, optimality in the pressure, measured in different norms, has been proved for both cases when the constrained specific storage coefficient c₀ is strictly positive and when c₀ is nonnegative. Numerical results illustrate the accuracy of the method and also show the effectiveness of the method to overcome the nonphysical pressure oscillations. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

P4‐decompositions of regular graphs

In this article, we show that every simple r‐regular graph G admits a balanced P₄‐decomposition if r ≡ 0(mod 3) and G has no cut‐edge when r is odd. We also show that a connected 4‐regular graph G admits a P₄‐decomposition if and only if |E(G)| ≡ 0(mod 3) by characterizing graphs of maximum degree 4 that admit a triangle‐free Eulerian tour. © 1999 John Wiley & Sons, Inc. J Graph Theory 31: 135–143, 1999     

On the nested refinement of quadrilateral and hexahedral finite elements and the affine approximation

Summary. It is shown in this paper that the optimal approximation property of bilinear or trilinear finite elements would be retained when the affine mapping is used to replace the Q₁ mapping on each element, if the grids are refined nestedly. The new method truncates the quadratic and cubic terms in reference mappings and produces constant Jacobians and Jacobian matrices. This would avoid a shortcoming of the quadrilateral and hexahedral elements where the integrals of rational functions have to be computed or approximated. Numerical tests verify the analysis.Mathematics Subject Classification (2000): 65N30, 65N50, 65N55Revised version received January 29, 2004     

Minimal extensions of Π01 classes

A minimal extension of a Π⁰₁ class P is a Π⁰₁ class Q such that P ? Q, Q – P is infinite, and for any Π⁰₁ class R, if P ? R ? Q, then either R – P is finite or Q – R is finite; Q is a nontrivial minimal extension of P if in addition P and Q′ have the same Cantor‐Bendixson derivative. We show that for any class P which has a single limit point A, and that point of degree ≤ 0 , P admits a nontrivial minimal extension. We also show that as long as P is infinite, then P does not admit any decidable nontrivial minimal extension Q. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixed finite element methods for general quadrilateral grids

We study a new mixed finite element of lowest order for general quadrilateral grids which gives optimal order error in the H(div)-norm. This new element is designed so that the H(div)-projection Π_h satisfies ∇ · Π_h = P_hdiv. A rigorous optimal order error estimate is carried out by proving a modified version of the Bramble-Hilbert lemma for vector variables. We show that a local H(div)-projection reproducing certain polynomials suffices to yield an optimal L²-error estimate for the velocity and hence our approach also provides an improved error estimate for original Raviart-Thomas element of lowest order. Numerical experiments are presented to verify our theory.