Trivariate spline approximation of divergence-free vector fields