The disclosed embodiments include a method, apparatus, and computer program product for generating hybrid computational meshes around complex and discrete fractures for the purpose of reservoir simulation. For example, one embodiment is a computer-implemented method for modeling three-dimensional (3D) geological fractures. The method includes the steps of receiving a set of fractures with geometry that has been discretized in a two-dimensional plane by a collection of line segments. The method generates closed loops around all the line segments associated with a fracture using a set of stadia and further generates shape elements within the closed loops. A constrained mesh around the closed loops of the set of fracture is generated to fill in a remainder space of the two-dimensional plane. The cell elements in the two-dimensional plane can then be extruded in a third dimension to form a three-dimensional model usable for performing reservoir simulations.