average gradient of FVM on unstructured mesh
 

When I compute viscous term, I need to use gradient on each edge of control volume. If I use cell-centered discretization, I can use average of gradients at cell centers on the two sides of an edge. How can I compute the gradient of speeds on a boundary edge. In ghost boundary cell, there is no gradient. Thanks for your help in advance.