3 results
Search Results
Now showing 1 - 3 of 3
- ItemMAC schemes on triangular Delaunay meshes(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2011) Eymard, Robert; Fuhrmann, Jürgen; Linke, AlexanderWe study two classical generalized MAC schemes on unstructured triangular Delaunay meshes for the incompressible Stokes and Navier-Stokes equations and prove their convergence for the first time. These generalizations use the duality between Voronoi and triangles of Delaunay meshes, in order to construct two staggered discretization schemes. Both schemes are especially interesting, since compatible finite volume discretizations for coupled convection-diffusion equations can be constructed which preserve discrete maximum principles. In the first scheme, called tangential velocity scheme, the pressures are defined at the vertices of the mesh, and the discrete velocities are tangential to the edges of the triangles. In the second scheme, called normal velocity scheme, the pressures are defined in the triangles, and the discrete velocities are normal to the edges of the triangles. For both schemes, we prove the convergence in $L^2$ for the velocities and the discrete rotations of the velocities for the Stokes and the Navier-Stokes problem. Further, for the normal velocity scheme, we also prove the strong convergence of the pressure in $L^2$. Linear and nonlinear numerical examples illustrate the theoretical predictions.
- ItemConvergence of a finite volume scheme to the eigenvalues of a Schrödinger operator(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2007) Koprucki, Thomas; Eymard, Robert; Fuhrmann, JürgenWe consider the approximation of a Schrödinger eigenvalue problem arising from the modeling of semiconductor nanostructures by a finite volume method in a bounded domain $OmegasubsetR^d$. In order to prove its convergence, a framework for finite dimensional approximations to inner products in the Sobolev space $H^1_0(Omega)$ is introduced which allows to apply well known results from spectral approximation theory. This approach is used to obtain convergence results for a classical finite volume scheme for isotropic problems based on two point fluxes, and for a finite volume scheme for anisotropic problems based on the consistent reconstruction of nodal fluxes. In both cases, for two- and three-dimensional domains we are able to prove first order convergence of the eigenvalues if the corresponding eigenfunctions belong to $H^2(Omega)$. The construction of admissible meshes for finite volume schemes using the Delaunay-Voronoï method is discussed. As numerical examples, a number of one-, two- and three-dimensional problems relevant to the modeling of semiconductor nanostructures is presented. In order to obtain analytical eigenvalues for these problems, a matching approach is used. To these eigenvalues, and to recently published highly accurate eigenvalues for the Laplacian in the L-shape domain, the results of the implemented numerical method are compared. In general, for piecewise $H^2$ regular eigenfunctions, second order convergence is observed experimentally.
- ItemConvergence of a finite volume scheme for the biharmonic problem(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Eymard, Robert; Gallou¨et, T.; Herbin, R.; Linke, A.We propose a finite volume scheme for the approximation of a biharmonic problem with Dirichlet boundary conditions. We prove that the piece-wise constant approximate solution converges to the exact solution, as well as the discrete approximate of the gradient and the discrete approximate of the Laplacian of the exact solution. These results are confirmed by numerical results.