Based on computational fluid dynamics Euler multiphase model, gas-liquid two-phase co-current upflow through packed bed was investigated. Tung & Siefken model, Schulenberg model and Attou k-fluid model were applied, respectively. The gas-particle, liquid-particle and gas-liquid drag force programs were implemented to CFD three-dimensional numerical simulation with user defined functions. A comparison of the experiment data has shown that Tung & Siefken's model was more suitable for the numerical simulation of gas-liquid two-phase upflow through porous media packed with 2~8 mm sphere under normal temperature and pressure. The CFD numerical simulated values of pressure drop and liquid hold-up by Tung & Siefken's model were consistent with the experiment results.