On Mar 31, 2013, at 11:35 PM, "Manav Bhatia" <[email protected]> wrote:
> Thanks, Derek.
>
> A related question about libMesh:
>
> If I need to output the solution values on the boundary nodes, can I use
> BoundaryMesh along with the IO classes for that?
>
> Meaning, if I do the following
>
> BoundaryMesh b_mesh(mesh.mesh_dimension()-1);
> mesh.boundary_info->sync(b_mesh);
> VTKIO(b_mesh).write_equation_systems("boundary_output.pvtu",
> equation_systems);
>
>
> Would this write the solution values on the boundary nodes?
Almost. Unfortunately it is a little more manual than this, but that's the
right idea. In your example the equation systems object can only be associated
with one mesh, so it won't work. The fix is to create a second equation systems
object that lives on the boundary and then extract the trace values.
This is not automated because in general the user will want physics-dependent
things like shear stress, heat flux, whatever...
I suppose we could automate the values-only case though.
-Ben
------------------------------------------------------------------------------
Own the Future-Intel® Level Up Game Demo Contest 2013
Rise to greatness in Intel's independent game demo contest.
Compete for recognition, cash, and the chance to get your game
on Steam. $5K grand prize plus 10 genre and skill prizes.
Submit your demo by 6/6/13. http://p.sf.net/sfu/intel_levelupd2d
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
