Pressure wave propagation in liquids inside containers of arbitrary three-dimensional shape can be solved by the computer code BEWAVE based on the Boundary Element Retarded Potential Method. This code may be applied to acoustic wave propagation in homogenous, inviscid but compressible fluids. BEWAVE computes, using a semi-implicit time-stepping technique, the transient solution of the velocity potential, velocity components tangential or normal to the boudary and pressure at the whole boundary. The dynamic loading is due to a pressure source representing an explosion or a dynamic boundary condition as defined by the user. Several new developments will be discussed. The implementation of a fast linear solver optimises the use of available core and disc space of the computer system. Furthermore, a cost-efffective procedure is discussed for checking on possible intersection of the line-of-sight between any pair of nodal points with the structures to see whether the influence between these points has to be excluded. The boundary conditions currently allow for rigid walls, given pressure, rupture discs, non-reflecting openings and prescribed flow. Pre- and post-processing has been optimised through a coupling with the interactive colour graphics package PATRAN. TM The boundary of the fluid domain, which may consist of both flat and curve parts, must be divided in quadratic quadrilateral or triangular elements. Recently the coupling of BEWAVE to a finite-difference program for the computation of one-dimensional pressure waves in complex piping systems has been investigated. The coupling method will be discussed and an industrial example for a nuclear steam generator in which a violent reaction is assumed to occur and which is coupled to a piping system will be given. The occurence and importance of the three-dimensional effects in this example will be discussed.