This paper deals with the problem of numerical approximation in the Cauchy-Dirichlet problem for a scalar conservation law with a flux function having finitely many discontinuities. The well-posedness of this problem was proved by Carrillo [J. Evol. Eq. 3 (2003) 687–705]. Classical numerical methods do not allow us to compute a numerical solution (due to the lack of regularity of the flux). Therefore, we propose an implicit Finite Volume method based on an equivalent formulation of the initial...

We present the convergence analysis of locally divergence-free discontinuous Galerkin methods for the induction equations which appear in the ideal magnetohydrodynamic system. When we use a second order Runge Kutta time discretization, under the CFL condition $\Delta t\sim h^{4/3}$, we obtain error estimates in $L^2$ of order $𝒪(\Delta t^2+h^{m+1/2})$ where $m$ is the degree of the local polynomials.

We present the convergence analysis of locally divergence-free discontinuous Galerkin methods for the induction equations which appear in the ideal magnetohydrodynamic system. When we use a second order Runge Kutta time discretization, under the CFL condition $\Delta t\sim h^{4/3}$, we obtain error estimates in L2 of order $𝒪(\Delta t^2+h^{m+1/2})$ where m is the degree of the local polynomials.

The convergence of Rothe’s method in Hölder spaces is discussed. The obtained results are based on uniform boundedness of Rothe’s approximate solutions in Hölder spaces recently achieved by the first author. The convergence and its rate are derived inside a parabolic cylinder assuming an additional compatibility conditions.

We formulate a finite element method for the computation of solutions to an anisotropic phase-field model for a binary alloy. Convergence is proved in the $H^1$-norm. The convergence result holds for anisotropy below a certain threshold value. We present some numerical experiments verifying the theoretical results. For anisotropy below the threshold value we observe optimal order convergence, whereas in the case where the anisotropy is strong the numerical solution to the phase-field equation does not...

Numerical solution of fractional order diffusion problems with homogeneous Dirichlet boundary conditions is investigated on a square domain. An appropriate extension is applied to have a well-posed problem on ${ℝ}^2$ and the solution on the square is regarded as a localization. For the numerical approximation a finite difference method is applied combined with the matrix transformation method. Here the discrete fractional Laplacian is approximated with a matrix power instead of computing the complicated...

Iterative approximation algorithms are successfully applied in parametric approximation tasks. In particular, reduced basis methods make use of the so-called Greedy algorithm for approximating solution sets of parametrized partial differential equations. Recently, a priori convergence rate statements for this algorithm have been given (Buffa et al. 2009, Binev et al. 2010). The goal of the current study is the extension to time-dependent problems, which are typically approximated using the POD–Greedy...

The problem of modeling acoustic waves scattered by an object with Neumann boundary condition is considered. The boundary condition is taken into account by means of the fictitious domain method, yielding a first order in time mixed variational formulation for the problem. The resulting system is discretized with two families of mixed finite elements that are compatible with mass lumping. We present numerical results illustrating that the Neumann boundary condition on the object is not always correctly...

The problem of modeling acoustic waves scattered by an object with Neumann boundary condition is considered. The boundary condition is taken into account by means of the fictitious domain method, yielding a first order in time mixed variational formulation for the problem. The resulting system is discretized with two families of mixed finite elements that are compatible with mass lumping. We present numerical results illustrating that the Neumann boundary condition on the object is not always...

We propose and analyse two convergent fully discrete schemes to solve the incompressible Navier-Stokes-Nernst-Planck-Poisson system. The first scheme converges to weak solutions satisfying an energy and an entropy dissipation law. The second scheme uses Chorin's projection method to obtain an efficient approximation that converges to strong solutions at optimal rates.

The incompressible MHD equations couple Navier-Stokes equations with Maxwell's equations to describe the flow of a viscous, incompressible, and electrically conducting fluid in a Lipschitz domain $\Omega \subset {ℝ}^3$. We verify convergence of iterates of different coupling and decoupling fully discrete schemes towards weak solutions for vanishing discretization parameters. Optimal first order of convergence is shown in the presence of strong solutions for a splitting scheme which decouples the computation of velocity...

A semidiscretization in time of a fourth order nonlinear parabolic system in several space dimensions arising in quantum semiconductor modelling is studied. The system is numerically treated by introducing an additional nonlinear potential. Exploiting the stability of the discretization, convergence is shown in the multi-dimensional case. Under some assumptions on the regularity of the solution, the rate of convergence proves to be optimal.

A semidiscretization in time of a fourth order nonlinear parabolic system in several space dimensions arising in quantum semiconductor modelling is studied. The system is numerically treated by introducing an additional nonlinear potential. Exploiting the stability of the discretization, convergence is shown in the multi-dimensional case. Under some assumptions on the regularity of the solution, the rate of convergence proves to be optimal.