Create 3D triangulation bounded by volume

I have two sets of data (3 really) and I want to create a delaunay triangulation that is constrained by a defined volume for the purposes of interpolation.

1. *.stl representation of a solid (watertight surface triangulation)
2. cloud data for points that lie on the surface or inside of the the solid defined by the *.stl file
3. nodes of a finite element model (point set in 3D space) to interpolate the data from 2 on to

For the sake of simplicity lets say that I have a cylinder with a hole through the middle (i.e. section of pipe). I have a surface boundary triangulation of the cylinder (*.stl), as well as data for points inside of the solid. If I take the points inside of the solid and use the scatteredInterpolant class a triangulation of the points is created. However, this triangulation will not actually represent the solid as there are tetrahedrons that will pass through the hole in the center of the cylinder based on the connectivity that matlab determines. What I want to know is.....

Is there a way to tell matlab that the tetrahedrons in the triangulation cannot have edges that penetrate the volume defined by the surface boundary (*.stl file)?

I have looked at alpha Volumes which actually works well for simple shapes, but when things get more complex using a single radius value to create the volume does not work so well.

It boils down to me wanting to create a delaunay triangulation bounded by an existing delaunay triangulation.