On Mon, 22 Mar 2010, Minq Q wrote: > I am sorry, The elements I am asking is not P1-iso-P2. The correct is > P1-iso-P2/P1 element. Here, the pressure and velocity variables are discrete > by the linear elements. But the velocity elements are "twice finer" than the > pressure element. That mean, 4 velocity elements = 1 pressure element. > > Can I do it with libMesh?
Only with some difficulty - we don't have those elements right now, so you'd have to create a new FEType, code the shape functions, and add a new macroelement quadrature rule for the integration. None of that is easy, but it's definitely doable. Coincidentally, that all describes my own first contributions to libMesh, in the CLOUGH finite elements. You'd have it a bit easier than I did: the Clough-Tocher shape functions are hairier, you've now got a macroelement quadrature rule (QClough) to use as an example, and now adaptivity works with arbitrary C0 elements without any special effort on your part. That may be more work than you wanted, but if you do give it a try, let us know and we'd be happy to add it in to the main library. --- Roy ------------------------------------------------------------------------------ Download Intel® Parallel Studio Eval Try the new software tools for yourself. Speed compiling, find bugs proactively, and fine-tune applications for parallel performance. See why Intel Parallel Studio got high marks during beta. http://p.sf.net/sfu/intel-sw-dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
