Outputting auxiliary cell-level data along with solution data

[mrjonm...@gmail.com]

I have a model where I'm using the boundary_id to specify boundary regions where certain boundary conditions apply, and I'd like to add the boundary_id numbers to the output, so I can see them spatially in Paraview along with my solution output. 

I'm currently trying to use a DataOutFaces object along with a class derived from DataPostprocessor to get these values and output them to .vtu files. I'm having trouble debugging it, and I was curious whether there is a simpler way to output data that are defined for an entire cell or face. 

Along with that, it would also be helpful to be able to output some global scalar data (e.g. simulation time) into the output files. If there's an easy way to to that, I'd appreciate it as well. In my current approach, I'm just going to set a variable named "time" to the appropriate value for all DoFs, but it seems like an inefficient way to do it. 

Thanks,

Jonathan

[Praveen C]

You can use this

https://dealii.org/developer/doxygen/deal.II/structGridOutFlags_1_1Vtu.html

https://dealii.org/developer/doxygen/deal.II/structDataOutBase_1_1VtkFlags.html

to set solution time.

Here is an example from my code

https://github.com/cpraveen/dflo/blob/master/src_mpi/output.cc#L66

To visualize boundary faces, VisIt has an option to make use of material id of cells and faces, which allows you to view subset of the data, e.g., plot only faces with material id = 100. You must save material id into vtk/vtu file and then use that in VisIt. I have never done this through deal.II but maybe there is a better way to do it in deal.II which I am not aware of.

Best

praveen

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dealii+un...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Corbin Foucart]

Hi Praveen,

Thanks for the helpful answer. I'm new to Visit; would you mind specifying how to view the subset of data (e.g., faces with boundary ID 0)? 

Best,

Corbin