A method is discussed for treating oscillations in flow through a heated boiling channel which connects two plenum regions of constant pressure difference. The approach is based on a numerical solution of pointwise difference equations representing conservation laws in fluid and metal. It is appropriate for spatially varying heat flux, flow with or without slip, and large perturbations. Details of the method, comparison with experiment, and features which distinguish it from other analytical models are discussed.