3D Results in TOUGH2

[Ethan]

Hello,

How are the continuous contour images in 3-D Results in PetraSim constructed from discrete grid block values? Are they interpolated linearly between grid block centers? Is there any type of weighting scheme?

Thanks,

Ethan

[Charlie Thornton]

Ethan,

The simulation output gives us data for the cell centers. We use that data to interpolate values for the cell corners. Based on the interpolated corner data, we create the isosurfaces using VTK's ContourFilter class which I believe uses an implementation of the marching cubes algorithm for 3D volumes. The interpolation is linear - based on the distance between the corner and neighboring cell centers.

- Charlie