In this paper, a class of cell centered finite volume schemes, on general unstructured meshes, for a linear convection-diffusion problem, is studied. The convection and the diffusion are respectively approximated by means of an upwind scheme and the so called diamond cell method [4]. Our main result is an error estimate of order h, assuming only the W2,p (for p>2) regularity of the continuous solution, on a mesh of quadrangles. The proof is based on an extension of the ideas developed in...

We study a finite volume method, used to approximate the solution of the linear two dimensional convection diffusion equation, with mixed Dirichlet and Neumann boundary conditions, on Cartesian meshes refined by an automatic technique (which leads to meshes with hanging nodes). We propose an analysis through a discrete variational approach, in a discrete H1 finite volume space. We actually prove the convergence of the scheme in a discrete H1 norm, with an error estimate of order O(h) (on meshes...

We extend the analysis of the recently proposed nonlinear EIS scheme applied to the partial eigenvalue problem. We address the case where the Rayleigh quotient iteration is used as the smoother on the fine-level. Unlike in our previous theoretical results, where the smoother given by the linear inverse power method is assumed, we prove nonlinear speed-up when the approximation becomes close to the exact solution. The speed-up is cubic. Unlike existent convergence estimates for the Rayleigh quotient...

In this paper, two algorithms are proposed to solve systems of algebraic equations generated by a discretization procedure of the weak formulation of boundary value problems for systems of nonlinear elliptic equations. The first algorithm, Newton-CG-MG, is suitable for systems with gradient mappings, while the second, Newton-CE-MG, can be applied to more general systems. Convergence theorems are proved and application to the semiconductor device modelling is described.

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...

The numerical minimization of the functional $\mathcal{F}\left(u\right)={\int }_{\mathrm{\Omega }}\varphi \left(x,{\nu }_u\right)\left|Du\right|+{\int }_{\partial \mathrm{\Omega }}\mu ud{\mathcal{H}}^{n-1}-{\int }_{\mathrm{\Omega }}\kappa udx$, $u\in BV\left(\mathrm{\Omega };\left\{-1,1\right\}\right)$ is addressed. The function $\varphi$ is continuous, has linear growth, and is convex and positively homogeneous of degree one in the second variable. We prove that $\mathcal{F}$ can be equivalently minimized on the convex set $BV\left(\mathrm{\Omega };\left[-1,1\right]\right)$ and then regularized with a sequence ${\left\{{\mathcal{F}}_{ϵ}\left(u\right)\right\}}_{ϵ}$, of stricdy convex functionals defined on $BV\left(\mathrm{\Omega };\left[-1,1\right]\right)$. Then both $\mathcal{F}$ and ${\mathcal{F}}_{ϵ}$, can be discretized by continuous linear finite elements. The convexity property of the functionals on $BV\left(\mathrm{\Omega };\left[-1,1\right]\right)$ is useful in the numerical minimization...

The least Steklov eigenvalue d1 for the biharmonic operator in bounded domains gives a bound for the positivity preserving property for the hinged plate problem, appears as a norm of a suitable trace operator, and gives the optimal constant to estimate the L2-norm of harmonic functions. These applications suggest to address the problem of minimizing d1 in suitable classes of domains. We survey the existing results and conjectures about this topic; in particular, the existence of a convex domain...

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multi-scale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a second-order...

We analyze the charge and spin distributions induced in an interacting electron system confined inside a semiconductor quantum wire with spin orbit interaction in the presence of an external magnetic field. The wire, assumed to be infinitely long, is obtained through lateral confinement in three different materials: GaAs, InAs, and InSb. The spin-orbit coupling, linear in the electron momentum is of both Rashba and Dresselhaus type. Within the Hartree-Fock approximation the many-body Hamiltonian...

In order to handle the flow of a viscous incompressible fluid in a porous medium with cracks, the thickness of which cannot be neglected, we consider a model which couples the Darcy equations in the medium with the Stokes equations in the cracks by a new boundary condition at the interface, namely the continuity of the pressure. We prove that this model admits a unique solution and propose a mixed formulation of it. Relying on this formulation, we describe a finite element discretization and derive...

In order to handle the flow of a viscous incompressible fluid in a porous medium with cracks, the thickness of which cannot be neglected, we consider a model which couples the Darcy equations in the medium with the Stokes equations in the cracks by a new boundary condition at the interface, namely the continuity of the pressure. We prove that this model admits a unique solution and propose a mixed formulation of it. Relying on this formulation, we describe a finite element discretization and derive...

This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is much cheaper...

This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is...