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 non-neighboring 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 ------------------------------------------------------------------------------ Gaining the trust of online customers is vital for the success of any company that requires sensitive data to be transmitted over the Web. Learn how to best implement a security strategy that keeps consumers' information secure and instills the confidence they need to proceed with transactions. http://p.sf.net/sfu/oracle-sfdevnl _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
