The equilibrium selective adsorption and fluxes of oxygen/nitrogen binary gas mixtures through carbon membranes are investigated at 303 K, respectively, using a grand canonical Monte Carlo simulation and a dual control volume grand canonical molecular dynamics method. The carbon membrane pores are modeled as slit-like pores with a two-dimensional structure where carbon atoms are placed according to the structure of graphite layers. The effect of the membrane thickness, bulk pressure, and pore width on the equilibrium selective adsorption and dynamic separation factor is discussed. Meanwhile a new iteration approach to calculate the flux and dynamic separation factor of binary gas mixtures through membranes is proposed, by which we can simulate the permeation and fluxes of gases through the membranes in the presence of pressure gradient and consider the effect of pressure and composition of low-pressure side in the meantime. The simulated results show that bulk pressure and membrane thickness have no effect on the equilibrium selectivity, but they have a great effect on the fluxes and dynamic separation factors of gases. The pore width impacts the equilibrium selectivity and dynamic separation factors strongly, especially when the pore width is very small. Molecular sieving dominates the separation of oxygen/nitrogen in non-equilibrium simulations. But due to the comparable molecular size of oxygen and nitrogen, we have to modify the carbon membranes in order to improve dynamic separation of atmosphere.