A new scheme for discretizing the P1 model on unstructured polygonal meshes is proposed. This scheme is designed such that its limit in the diffusion regime is the MPFA-O scheme which is proved to be a consistent variant of the Breil-Maire diffusion scheme. Numerical tests compare this scheme with a derived GLACE scheme for the P1 system.

We present an asynchronous multi-domain time integration algorithm with a dual domain decomposition method for the initial boundary-value problems for a parabolic equation. For efficient parallel computing, we apply the three-field domain decomposition method with local Lagrange multipliers to ensure the continuity of the primary unknowns at the interface between subdomains. The implicit method for time discretization and the multi-domain spatial decomposition enable us to use different time steps...

We present the current Reduced Basis framework for the efficient numerical approximation of parametrized steady Navier–Stokes equations. We have extended the existing setting developed in the last decade (see e.g. [S. Deparis, SIAM J. Numer. Anal. 46 (2008) 2039–2067; A. Quarteroni and G. Rozza, Numer. Methods Partial Differ. Equ. 23 (2007) 923–948; K. Veroy and A.T. Patera, Int. J. Numer. Methods Fluids 47 (2005) 773–788]) to more general affine and nonaffine parametrizations (such as volume-based...

This paper deals with nonlinear diffusion problems which include the Stefan problem, the porous medium equation and cross-diffusion systems. We provide a linear scheme for these nonlinear diffusion problems. The proposed numerical scheme has many advantages. Namely, the implementation is very easy and the ensuing linear algebraic systems are symmetric, which show low computational cost. Moreover, this scheme has the accuracy comparable to that of the wellstudied nonlinear schemes and make it possible...

An energy analysis is carried out for the usual semidiscrete Galerkin method for the semilinear equation in the region$\Omega$ (E) ${\left(t{u}_t\right)}_t={\sum }_{i,j=1}{\left(a_{ij}\left(x\right)u_{x_i}\right)}_{x_j}-a_0\left(x\right)u+f\left(u\right)$, subject to the initial and boundary conditions, $u=0$ on $\partial \Omega$ and $u(x,0)=u_0$. (E) is degenerate at $t=0$ and thus, even in the case $f\equiv 0$, time derivatives of $u$ will blow up as $t\to 0$. Also, in the case where $f$ is locally Lipschitz, solutions of (E) can blow up for $t>0$ in finite time. Stability and convergence of the scheme in $W^{2,1}$ is shown in the linear case without assuming $u_{tt}$ (which can blow up as $t\to 0$ is...

We consider the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water. The difficulty in this system comes from the coupling terms involving some derivatives of the unknowns that make the system nonconservative, and eventually nonhyperbolic. Due to these terms, a numerical scheme obtained by performing an arbitrary scheme to each layer, and using time-splitting or other similar techniques leads to instabilities in...

We study error estimates and their convergence rates for approximate solutions of spectral Galerkin type for the equations for the motion of a viscous chemical active fluid in a bounded domain. We find error estimates that are uniform in time and also optimal in the L2-norm and H1-norm. New estimates in the H(-1)-norm are given.

We describe a numerical method for the equation $u_t+p{u}_x-\epsilon u_{xx}=f$ in $(0,1)\times (0,T)$ with Dirichlet boundary and initial conditions which is a combination of the method of characteristics and the finite-difference method. We prove both an a priori local error-estimate of a high order and stability. Example 3.3 indicates that our approximate solutions are disturbed only by a minimal amount of the artificial diffusion.

In this paper Rothe’s classical method is extended so that it can be used to solve some linear parabolic boundary value problems in non-cylindrical domains. The corresponding existence and uniqueness theorems are proved and some further results and generalizations are discussed and applied.

In this paper, we discuss an hp-discontinuous Galerkin finite element method (hp-DGFEM) for the laser surface hardening of steel, which is a constrained optimal control problem governed by a system of differential equations, consisting of an ordinary differential equation for austenite formation and a semi-linear parabolic differential equation for temperature evolution. The space discretization of the state variable is done using an hp-DGFEM, time and control discretizations are based on a discontinuous Galerkin...

In this paper, we discuss an hp-discontinuous Galerkin finite element method (hp-DGFEM) for the laser surface hardening of steel, which is a constrained optimal control problem governed by a system of differential equations, consisting of an ordinary differential equation for austenite formation and a semi-linear parabolic differential equation for temperature evolution. The space discretization of the state variable is done using an hp-DGFEM, time and control discretizations are based on a discontinuous Galerkin...

An implicit-explicit (IMEX) method is developed for the numerical solution of reaction-diffusion equations with pure Neumann boundary conditions. The corresponding method of lines scheme with finite differences is analyzed: explicit conditions are given for its convergence in the ‖·‖∞ norm. The results are applied to a model for determining the overpotential in a proton exchange membrane (PEM) fuel cell.

In this paper we study the finite element approximations to the parabolic and hyperbolic integrodifferential equations and present an immediate analysis for global superconvergence for these problems, without using the Ritz projection or its modified forms.