Screenshot instructions:
Windows
Mac
Red Hat Linux
Ubuntu
Click URL instructions:
Rightclick on ad, choose "Copy Link", then paste here →
(This may not be possible with some types of ads)
From: Vasilis Vavourakis <vasvav@gm...>  20110110 17:50:57

hi all and happy new year :) i'd like to ask if anyone has implemented libmesh in Boundary Element Methods. i'm currently trying to built a simple 2D BEM solver for the Laplace PDE "Grad^2(u) = 0". since i'm not that a libmesh expert i face some difficulties in the implementation: (1) in BEM mostly boundary integrations occur, which means that in 2D problems only 1D discretization is needed. however, since only Edge elements are "loaded" to the MeshBase class, when i try to evaluate the normal vectors to each quadrature point for each boundary (edge) element through the "get_normals" function then i get an error. one option would be to load a 2D mesh (of triangles and quadrilaterals) and perform only boundary integration to all nonneighboring side elements, but this way i would load into the System the nodes located inside the domain. any ideas on that? (2) the system matrix "A" of the final linear system "A*x=b" in not banded as in FEM, but full populated... this means that parallelization is not that easy to implement, right? i have more questions but i think that the above are the most "difficult" ones i have so far... in case anyone in the past has worked in that area, their ideas and help would be much appreciated! many thanks, Vas 