Here is a suggested answer from Wolfram's team
Needs["NDSolve`FEM`"]
Needs["OpenCascadeLink`"]
cyl = OpenCascadeShape[Cylinder[{{-9, 0, 0}, {9, 0, 0}}, 7]];
dg = DiscretizeGraphics[ExampleData[{"Geometry3D", "SpaceShuttle"}]];
orbiter = OpenCascadeShape[ToBoundaryMesh[dg]];
shape = OpenCascadeShapeDifference[cyl, orbiter];
bmesh = OpenCascadeShapeSurfaceMeshToBoundaryMesh[shape];
bmesh["Wireframe"]