Numerical schemes for unsteady fluid flow through collapsible tubes

The study of fluid flow through compliant tubes is a fluid-structure type problem, in which a dynamic equilibrium is maintained between the fluid and the tube wall. The analogy between this flow and gas dynamics initiated the use of a number of numerical methods which were originally developed to solve compressible flow in rigid ducts. In this study we investigate the solutions obtained by applying the Lax-Wendroff and MacCormack schemes to one-dimensional incompressible flow through a straight collapsible tube. The time-evolving numerical results were compared with exact steady-state solutions. For boundary conditions which were held fixed after a prescribed rise time, the unsteady numerical solution converges to the exact steady-state solution with very good accuracy. The stability and accuracy of all the methods depend on the amount of viscous pressure loss dictated by wall friction. Flows with undamped oscillations cannot, however, be solved with these techniques.