Theoretical Convergence of Numerical Schemes for Couped 3D-1D Partial Differential Equations