Numerical Simulation of Flow Through Biofluid Devices

Abstract

The results of a numerical simulation of flow through an artificial heart and through an artificial tilting-disk heart valve are presented. The simulation involves solving the incompressible Navier-Stokes equations; the solution process is described. The details and difficulties of mod eling these particular geometries are discussed. The arti ficial heart geometry uses a single moving grid, and the valve computation uses an overlaid-grid approach with one moving grid and one stationary grid. The equations must be solved iteratively for each discrete time step of the computations, requiring a significant amount of computing time. It is particularly difficult to analyze and present the fluid physics represented by these calcula tions because of the time-varying nature of the flow, and because the flows are internal. The use of three- dimensional graphics and scientific visualization tech niques have become instrumental in solving these problems.

Get full access to this article

View all access options for this article.

References

Beneck, J.A. , Buning, P.G. , and Steger, J.L. 1985. A 3-D chimera grid embedding technique. AIAA Paper 85-1523-CP.

Bjorstad, P.E. 1980. Numerical solution of the bi-harmonic equation. Ph.D. dissertation, Stanford University.

Dougherty, F.C. , Benek, J.A. , and Steger, J.L. 1985. On applications of chimera grid schemes to store separation . NASA TM 88193. Moffett Field, Calif.: NASA Ames Research Center.

Figliola, R.S. , and Mueller, T.J. 1981. On the hemolytic and thrombogenic potential occluder prosthetic heart valves from in-vitro measurements. J. Biomech. Engrg. 103:83-90.

Idelsohn, S.R. , Costa, L.E. , and Ponso, R. 1985. A comparative computational study of blood flow through prosthetic heart valves using the finite element method. J. Fluid Mech. 18:97-115.

Liepsch, D. , Moravec, S. , Rastogi, A.K. , and Vlachos, N.S. 1983. Measurement and calculations of laminar flow in a ninety-degree bifurcation. J. Biomech. 15:753-766.

Liepsch, D. , Steiger, H.J. , Poll, A. , and Reulen, H.-J. 1987. Hemodynamics stress in lateral sacular aneurysms. J. Biorheology 24:689-710.

Peskin, C.S. 1977. Numerical analysis of blood flow in the heart. J. Comput. Phys. 25:220-252.

Peskin, C.S. 1982. The fluid dynamics of heart valves: experimental, theoretical and computational methods. Ann. Rev. Fluid Mech. 14:235-259.

10.

Roe, P.L. 1981. Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys. 43:357-372.

11.

Rogers, S.E. , and Kwak, D. 1989. Numerical solution of the incompressible Navier-Stokes equations for steady-state and time-dependent problems. AIAA Paper 89-0463.

12.

Rogers, S.E. , and Kwak, D. 1990. An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations. AIAA J. 28(2):253-262.

13.

Tarbell, J.M. , Gunshinan, J.P. , Geselowitz, D.B. , Rosenburg, G. , Shung, K.K. , and Pierce, W.S. 1986. Pulse ultrasonic Doppler velocity measurements inside a left ventricular assist device. J. Biomech. Engrg., Trans. ASME 108:232-238.

14.

Tiderman, W.G. , Steinle, M.J. , and Phillips, W.M. 1986. Two-component laser velocimeter measurements downstream of heart valve prostheses in pulsatile flow. J. Biomech. Engrg., Trans. ASME 108:59-64.

15.

Underwood, F.N. , and Mueller, T.J. 1977. Numerical study of the steady axisymmetric flow through a disk-type prosthetic heart valve in a constant diameter chamber . J. Biomech. Engrg. 99:91-97.

16.

Underwood, F.N. , and Mueller, T.J. 1979. Numerical study of the steady axisymmetric flow through a disk-type prosthetic heart valve in a aortic-shaped chamber. J. Biomech. Engrg. 101:198-204.

17.

Yoganathan, A.P. , Concoran, W.H. , and Harrison, E.E. 1979a. In vitro velocity measurements in the vicinity of aortic prostheses. J. Bio-mechanics 12:135-152.

18.

Yoganathan, A.P. , Concoran, W.H. , and Harrison, E.E. 1979b. Pressure drops across prosthetic aortic heart valves under steady and pulsatile flow. J. Bio-mechanics 12:153-164.