In this talk, I will present our recent work on mathematical models, asymptotic analysis and numerical simulation for degenerate dipolar quantum gas. As preparatory steps, I begin with the three-dimensional Gross-Pitaevskii equation with a long-range dipolar interaction potential which is used to model the degenerate dipolar quantum gas and reformulate it as a Gross-Pitaevskii-Poisson type system by decoupling the two-body dipolar interaction potential which is highly singular into short-range (or local) and long-range interactions (or repulsive and attractive interactions). Based on this new mathematical formulation, we prove rigorously existence and uniqueness as well as nonexistence of the ground states, and discuss the existence of global weak solution and finite time blowup of the dynamics in different parameter regimes of dipolar quantumgas. In addition, a backward Euler sine pseudospectral method ispresented for computing the ground states and a time-splitting sine pseudospectral method is proposed for computing the dynamics of dipolar BECs. Our new numerical methods avoid evaluating integrals with high singularity and thus they are more efficient and accurate than those numerical methods currently used. In addition, new mathematical formulations intwo-dimensions and one dimension for dipolar quantum gas are obtained when the external trapping potential is highly confined in one or two directions.