Solving Hamilton-Jacobi-Bellman (HJB) equations is essential in feedback optimal control. Using the solution of HJB equations, feedback optimal control laws can be implemented in real-time with minimum computational load. However, except for systems with two or three state variables, numerically solving HJB equations for general nonlinear systems is unfeasible due to the curse of dimensionality. In this paper, we develop a new computational method of solving HJB equations. The method is causality free, which enjoys the advantage of perfect parallelism on a sparse grid. Compared with dense grids, a sparse grid has a significantly reduced size which is feasible for systems with relatively high dimensions, such as 6-D HJB equations for the attitude control of rigid bodies. The method is applied to the optimal attitude control of a satellite system using momentum wheels. The accuracy of the numerical solution is verified at a set of randomly selected sample points.