I failed to copy my reply to the list yesterday; obviously Roy has already responded to this question.
Roy: Is my memory correct in that this is the same "perfectly unnested" element (with respect to h-refinement) we looked at for shallow water equations a while back? -- John ---------- Forwarded message ---------- From: John Peterson <[email protected]> Date: Mon, Mar 22, 2010 at 6:18 PM Subject: Re: [Libmesh-users] P2-iso-P1 To: Minq Q <[email protected]> On Mon, Mar 22, 2010 at 5:29 PM, Minq Q <[email protected]> wrote: > Hi again, > >> assume you mean P2-P1 Taylor-Hood elements? Check out libMesh examples 13 > and 18. > > 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. Yeah, according to Gresho and Sani, this is like a P1 P1 (triangular) element on a 4-patch macroelement in 2D, and something equivalent in 3D. I'm not sure what the benefit is though, since they also say it's a 1st-order element. I think this is an element we looked at while doing shallow-water equations a while back, but we abandoned the idea because the function spaces don't nest under refinement... > Can I do it with libMesh? Not as-is, but I'll let Roy comment further. He has implemented what I would consider much more complicated macroelements in Libmesh, and would probably be able to give you a much better answer. ------------------------------------------------------------------------------ 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
