Numerical simulation for the propagation of nonlinear pulsatile waves in arteries.

The behavior of nonlinear pulsatile flow of incompressible blood contained in an elastic tube is examined. The theory takes into account the nonlinear convective terms of the Navier-Stokes equations. The motion of the arterial wall is characterized by a set of linearized differential equations. The region bounded by the flexible arterial wall is mapped into a fixed area in which numerical discretization takes place. The finite element method (Galerkin weighted residual approach) is used for the solution of this nonlinear system. The results obtained are pressure distribution, velocity profile, flow rate and wall displacements along the elastic tube (20 cm long).