On Wed, 16 Dec 2015, Harshad Sahasrabudhe wrote: > Sorry I meant to say that the control volume will go beyond the > element if the element has obtuse angles. I want to solve an > elliptic second order PDE with vertex centered finite volume method.
Ah! Vertex centered makes your question make sense, but also changes things. libMesh doesn't have any built in features for dealing with dual meshes, so you'd have to handle that all on your own. Are you sure you're not creating an XY problem here? Your original concern was the non-diagonal mass matrix that FEM gives you, right? Wouldn't simply mass lumping with a Lagrange basis be an easy fix? --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
