This paper is concerned with optimal power flow (OPF), which is the problem of optimizing the transmission of electricity in power systems. Our main contributions are as follows: (i) we propose a novel parabolic relaxation, which transforms non-convex OPF problems into convex quadratically-constrained quadratic programs (QCQPs) and can serve as an alternative to the common practice semidefinite programming (SDP) and second-order cone programming (SOCP) relaxations, (ii) we propose a penalization technique which is compatible with the SDP, SOCP, and parabolic relaxations and guarantees the recovery of feasible solutions for OPF, under certain assumptions. The proposed penalized convex relaxation can be used sequentially to find feasible and near-globally optimal solutions for challenging instances of OPF. Extensive numerical experiments on small and large-scale benchmark systems corroborate the efficacy of the proposed approach. By solving a few rounds of penalized convex relaxation, fully feasible solutions are obtained for benchmark test cases from [1]-[3] with as many as 13659 buses. In all cases, the solutions obtained are not more than 0.32% worse than the best-known solutions.