A numerical code, Gemini, based on the implicit multifield method (IMF) of Harlow and Amsden for Eulerian two-fluid modelling, is used to simulate the fluid dynamics of bubbling fluidized beds, assuming no turbulence in the gas or solid phase. The paper gives a formulation of the equations of motion and empirical closure laws in general curvilinear coordinates for calculation of the fluid dynamics in beds with complex internal geometries. A special discretization method for general curvilinear structured grids with multiblock connectivity is implemented, and two-dimensional non-stationary calculations are performed for a bed with a cross-sectional width of 0.3 m, containing two horizontal heat exchanger tubes. The local visible bubble flow and the gas and particle motion around the tubes are briefly discussed and compared with experimental fluid dynamic results at different pressure levels.