In this paper we solve the time-dependent incompressible Navier-Stokes equations by splitting the non-linearity and incompressibility, and using discontinuous or continuous finite element methods in space. We prove optimal error estimates for the velocity and suboptimal estimates for the pressure. We present some numerical experiments.
We present in this paper a stability study concerning finite volume schemes applied to the two-dimensional Maxwell system, using rectangular or triangular meshes. A stability condition is proved for the first-order upwind scheme on a rectangular mesh. Stability comparisons between the Yee scheme and the finite volume formulation are proposed. We also compare the stability domains obtained when considering the Maxwell system and the convection equation.
The paper is devoted to the computation of two-phase flows in a porous medium when applying the two-fluid approach. The basic formulation is presented first, together with the main properties of the model. A few basic analytic solutions are then provided, some of them corresponding to solutions of the one-dimensional Riemann problem. Three distinct Finite-Volume schemes are then introduced. The first two schemes, which rely on the Rusanov scheme, are shown to give wrong approximations in some...
We address the numerical approximation of the two-phase Stefan problem and discuss an adaptive finite element method based on rigorous a posteriori error estimation and refinement/coarsening. We also investigate how to restrict coarsening for the resulting method to be stable and convergent. We review implementation issues associated with bisection and conclude with simulations of a persistent corner singularity, for which adaptivity is an essential tool.
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
We discuss the occurrence of oscillations when using central schemes of the Lax-Friedrichs type (LFt), Rusanov’s method and the staggered and non-staggered second order Nessyahu-Tadmor (NT) schemes. Although these schemes are monotone or TVD, respectively, oscillations may be introduced at local data extrema. The dependence of oscillatory properties on the numerical viscosity coefficient is investigated rigorously for the LFt schemes, illuminating also the properties of Rusanov’s method. It turns...
We discuss the occurrence of oscillations when using central schemes of the Lax-Friedrichs type (LFt), Rusanov's method and the staggered and non-staggered second order Nessyahu-Tadmor (NT) schemes. Although these schemes are monotone or TVD, respectively, oscillations may be introduced at local data extrema. The dependence of oscillatory properties on the numerical viscosity coefficient is investigated rigorously for the LFt schemes, illuminating also the properties of Rusanov's method. It turns...
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...
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 in 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, 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...