Self-pressurization of cylindrical container of cryogen is numerically analyzed. The container is axi-symmetric and heated from side wall with constant heat flux. Natural convection by external heat flux is studied numerically using finite difference method. Oxygen, nytrogen and hydrogen are working fluids in this paper. Liquid is considered incompressible fluid and vapor is assumed to behave as gas meeting with virial equation of gas. The Second virial coefficients of gas are obtained from Lennard-jones model. The important variables which have effects on self-pressurization are external heat flux, heat capacity of wall and initial ullage in container. The most important variable of them is external heat flux. The pressure rise calculated from the virial gas model is slightly different from that calculated using Ideal gas model for oxygen.