Abstract. Many problems of fundamental and practical importance have multiple scale solutions. The direct numerical solution of multiple scale problems is difficult to obtain even with modern supercomputers. The major difficulty of direct solutions is the scale of computation. The ratio between the largest scale and the smallest scale could be as large as 10 5 in each space dimension. From an engineering perspective, it is often sufficient to predict the macroscopic properties of the multiple-scale systems, such as the effective conductivity, elastic moduli, permeability, and eddy diffusivity. Therefore, it is desirable to develop a method that captures the small scale effect on the large scales, but does not require resolving all the small scale features.The purpose of this lecture note is to review some recent advances in developing multiscale finite element (volume) methods for flow and transport in strongly heterogeneous porous media. Extra effort is made in developing a multiscale computational method that can be potentially used for practical multiscale for problems with a large range of nonseparable scales. Some recent theoretical and computational developments in designing global upscaling methods will be reviewed. The lectures can be roughly divided into 4 parts. In part 1, I will review some homogenization theory for elliptic and hyperbolic equations. This homogenization theory provides the critical guideline for designing effective multiscale methods. In part 2, I will review some recent developments of multiscale finite element (volume) methods. We also discuss the issue of upscaling one-phase, two-phase flows through heterogeneous porous media and the use of limited global information in multiscale finite element (volume) methods. In part 4, we will consider multiscale simulations of two-phase flow immiscible flows using a flow-based adaptive coordinate, and introduce a theoretical framework which enables us to perform global upscaling for heterogeneous media with long range connectivity.