Flying fish is a family of unique aerial-aquatic animals, which can both swim in the water and glide over the sea surface. Most previous studies on their aerodynamic characteristics were based on field observations or measurements of their morphometric parameters. In the present study, we consider three different flying fish models, of which the preliminary one mimics the Cypselurus hiraii in the pectoral fin morphology, following a previous wind tunnel experiment (Park and Choi 2010 J. Exp. Biol. 213 3269-79). Their aerodynamic performances are numerically studied by the computational fluid dynamics (CFD) method. The maximum lift force coefficient of 1.03 is reached at the angle of attack [Formula: see text], and the maximum lift-to-drag ratio of 4.7 is achieved at [Formula: see text]. By choosing appropriately the center of gravity, the flying fish model is proved to be longitudinally stable, according to the negative slope of pitching moment profile. Furthermore, we build a three-degrees-of-freedom (3-DOF) dynamic model in the longitudinal plane based on the aerodynamic coefficients obtained in our simulations, to predict its gliding performance. The results show that the flying fish can achieve a distance up to 45.4 m, and reach a height of 13.2 m, indicating an extraordinary gliding performance. Our numerical simulations are consistent with previous experimental results and theoretical prediction, which can be taken as the basis of further research on robotic flying fish.