SUMMARYA numerical simulation of reaction injection molding (RIM) of polymeric foam is developed, using a ÿnite volume method (FVM). In this study we predict mold ÿlling with a variable-density uid that ÿlls a mold by self-expansion. We deal with two-dimensional, isothermal cases. With the assumptions of ideal mixing and rapid bubble nucleation, the foam is modelled as a continuum with a time-dependent density. The continuum is assumed to be a Newtonian uid. We develop a pressure-based FVM for unstructured meshes that includes the SIMPLE algorithm with treatment of uid compressibility. Cellbased, co-located storage is used for all physical variables. To treat the moving interface, an explicit high-resolution interface capturing method is used. Foam ow in a slit is investigated, and the numerical calculations are in good agreement with an approximate analytic solution. For fountain ow in a rectangular cavity, the shape of the ow front is atter and the traces of the particles are more complicated for an expanding foam than for a constant-density uid. An example of mold ÿlling by an expanding foam demonstrates the geometric exibility of the method.