A Numerical Simulation of Peristaltic motion in the Ureter Using Fluid Structure Interactions

An axisymmetric model with fluid-structure interactions (FSI) is introduced and solved to perform ureter flow and stress analysis. The Navier-Stokes equations are solved for the fluid and a linear elastic model for ureter is used. The finite element equations for both the structure and the fluid were solved by the Newton-Raphson iterative method. Our results indicated that shear stresses were high around the throat of moving contracted wall. The pressure gradient magnitude along the ureter wall and the symmetry line had the maximum value around the throat of moving contracted wall which decreased as the peristalsis propagates toward the bladder. The flow rate at the ureter outlet at the end of the peristaltic motion was about 650 mm³/s. During propagation of the peristalsis toward the bladder, the inlet backward flow region was limited to the areas near symmetry line but the inner ureter backward flow regions extended to the whole ureter contraction part. The backward flow was vanished after 1.5 seconds of peristalsis propagation start up and after that time the urine flow was forward in the whole ureter length, so reflux is more probable to be present at the beginning of the wall peristaltic motion.