Stabilizing the retarded potential method for transient fluid-structure interaction problems

Abstract

The Retarded Potential (RP) method, which is a boundary element technique and non-local in both space and time, is employed to discretize the fluid domain for the analysis of transient fluid–structure interaction problems. The retarded potential analysis program RPFS is coupled to the ABAQUS non-linear finite element code to form ABAQUS/RPFS. The standard RP is inherently unstable for time steps below a critical time step that is equal to the maximum distance in the fluid divided by the wave speed. A technique referred to as the Figueiredo method is used to convert the standard RP differential-delay equations for the fluid to simply delay equations, which are more stable. The Figueiredo approach extends the stability range of the standard RP by a factor of approximately 10–20, but this time step is still not small enough to be useful for analysis. Digital signal processing methods are used to further stablize the response of the fluid by removing the oscillating high-frequency noise in the time histories of the solution without introducing phase shifting or any significant damping. Stability of the coupled system is achieved by not extrapolating the structural accelerations. ABAQUS/RPFS is applied to both a rigid and elastic sphere subjected to a plane wave, and the results using the full time histories required are completely stable and quite accurate. With this procedure, the retarded potential method may yet prove to be a valuable analysis tool for transient fluid–structure interaction problems. © 1997 John Wiley & Sons, Ltd.