The convolution structure for Jacobi function expansions