SUMMARYA fully coupled numerical model has been developed which describes multiphase¯uid¯ow through soil: namely gas, water and a non-aqueous phase liquid (NAPL) in a deforming porous media for subsurface systems. A multiphase¯ow model, based on the two-phase¯ow model of Brooks and Corey, is presented to express the dependence of saturation and relative permeability on the capillary pressure. Non-linear saturation and relative permeability functions are incorporated into a Galerkin ®nite element model which is subsequently used to simulate multiphase immiscible¯uid¯ow under saturated and unsaturated conditions in porous media. The governing partial dierential equations, in terms of soil displacements and¯uid pressures, which are coupled and non-linear, are solved by the ®nite element method. Numerical implementation of the formulation is discussed, and example problems demonstrate the model and solution procedure. #