Patient specific circulation model

• US 7,191,110 B1
• Filed: 09/20/1999
• Issued: 03/13/2007
• Est. Priority Date: 02/03/1998
• Status: Expired due to Term

• Current Assignee: Board of Trustees of The University of Illinois
• Original Assignee: Board of Trustees of The University of Illinois
• Inventors: Charbel, Fady T., Clark, Marlyn E., Zhao, Meide, Loth, Francis, Sadler, Lewis, Alperin, Noam
• Primary Examiner(s): JONES, HUGH M
• Application Number:
  US09/400,365
• Time in Patent Office: 2,731 Days
• Field of Search: 703/11, 703/2, 703/6, 434/262, 434/268, 434/270
• US Class Current: 703/11
• CPC Class Codes: G09B 23/28 for medicine
• CPC Class Codes: G16H 30/20 for handling medical images, e.g. DICOM, HL7 or PACS
• CPC Class Codes: G16H 50/50 for simulation or modelling of medical disorders

Abstract

A method and apparatus are provided for modeling circulation in a living subject. The method includes the steps of developing a model for living subjects in general and correcting the model to substantially conform to the overall cerebral physiology of the living subject. The method further includes the step of calculating a cerebral flow of the living subject based upon the corrected model and a selected cerebral blood flow perturbation.

Claims

1. A method of modeling circulation in a living subject, such method comprising the steps of:

simulating the fluid dynamics of an arterial network, wherein the simulation models blood flow through a plurality of arterial segments including one or more terminal efferent vessels;

adapting the simulation to substantially conform to a specific arterial anatomy of the living subject;

forcing the simulation with a forcing function made up of one or more flow-time or pressure-time signatures;

calculating a blood flows in the arterial network based upon the forced simulation;

measuring a blood flow in the living subject;

correcting the simulation based on the measured and calculated blood flows;

modifying the simulation to model a particular interventional procedure; and

,

calculating a post-procedure blood flow in a selected arterial segment using the modified simulation in order to predict an outcome of the actual interventional procedure performed in the living subject.

2. The method of modeling as in claim 1 wherein the simulated arterial network includes the Circle of Willis.

3. The method of modeling as in claim 1 wherein the step of adapting the simulation to substantially conform to the living subject'"'"'s anatomy further comprises conforming a vessel of the simulation with a corresponding vessel in an image of the living subject.

4. The method of modeling as in claim 3 wherein the step of adapting the simulation to substantially conform to the living subject'"'"'s anatomy further comprises measuring a diameter of the corresponding vessel in the image of the living subject.

5. The method of modeling as in claim 4 further comprising localizing the corresponding vessel in three-dimensional space and tracing a boundary into adjacent areas in three-dimensional space to locate respective ends of the corresponding vessel.

6. The method of modeling as in claim 1 wherein the simulation of the arterial network includes a one-dimensional, explicit, finite difference algorithm based upon a conservation of mass equation, a Navier-Stokes momentum equation, and an equation of state relating local pressure to local artery size.

7. The method of modeling as in claim 1 wherein the simulation is forced with a flow measurement obtained from the living subject.

8. The method of modeling as in claim 1 wherein the simulation is forced with a pressure-time signature obtained from a prototypical measurement.

18. The method of claim 1 further comprising the step of obtaining a flow measurement in the living subject by phase contrast magnetic resonance angiography.

19. The method of claim 1 further comprising the step of obtaining a flow measurement in the living subject by a Doppler flow measurement.

21. The method of claim 1 wherein the arterial network is divided into a plurality of sectors, wherein each sector is terminated by a terminal efferent vessel and has a primary input vessel, and further wherein the terminal efferent vessels are modeled as terminal efferent resistances, the method further comprising:

(a) measuring blood flows in the living subject corresponding to the primary input flows of the sectors;

(b) adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow for each such sector.

22. The method of claim 21 further comprising repeating the calculation of blood flows in the terminal efferent vessels and adjustment of terminal efferent resistances a selected number of times.

23. The method of claim 21 further comprising repeating the calculation of blood flows in the terminal efferent vessels and adjustment of terminal efferent resistances until the differences between the calculated terminal efferent flows and the corresponding measured primary input flows for each sector are within specified limit values.

24. The method of claim 21 wherein the terminal efferent resistance of a sector is adjusted by multiplying the terminal efferent resistance by the ratio of the measured primary input flow to the calculated terminal efferent flow.

25. The method of claim 21 wherein one or more sectors has one or more secondary input or output vessels and further comprising:

calculating flows through the secondary input and output vessels; and

,

adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow adjusted for the calculated flows in the secondary input and output vessels.

9. Apparatus for modeling circulation within a living subject, such apparatus comprising:

a computerized model of an arterial network made up of a plurality of arterial segments including one or more terminal efferent vessels, wherein the apparatus includes means for calculating blood flows in the arterial network when the model is forced with a forcing function;

means for adapting the model of the arterial network to substantially conform to a specific arterial anatomy of the living subject;

means for measuring a blood flow in the living subject;

means for correcting the model based upon the calculated and measured flows;

means for modifying the model to reflect a particular interventional procedure; and

,

means for calculating a post-procedure blood flow in a selected arterial segment using the modified model in order to predict an outcome of the actual interventional procedure performed in the living subject.

10. The apparatus for modeling as in claim 9 wherein the circulation model further comprises the Circle of Willis.

11. The apparatus for modeling as in claim 9 wherein the means for measuring blood flow is a phase contrast magnetic resonance angiography flow measurement system.

12. The apparatus for modeling as in claim 11 wherein the means for adapting the model to substantially conform to the living subject'"'"'s anatomy further comprises means for selecting a vessel of the model and a corresponding vessel in an image of the living subject.

13. The apparatus for modeling as in claim 12 wherein the means for adapting the model to substantially conform to the living subject'"'"'s anatomy further comprises means for measuring a diameter of the corresponding vessel.

14. The apparatus for modeling as in claim 13 further comprising means for localizing the corresponding vessel in three-dimensional space and tracing a boundary into adjacent areas in three-dimensional space to locate respective ends of the corresponding vessel.

15. The apparatus for modeling as in claim 9 wherein the computerized simulation model includes a one-dimensional, explicit, finite difference algorithm based upon a conservation of mass equation, a Navier-Stokes momentum equation, and an equation of state relating local pressure to local artery size.

16. The apparatus for modeling as in claim 9 wherein the model is forced with a flow measurement obtained from the living subject.

17. The apparatus for modeling as in claim 9 wherein the model is forced with a pressure-time signature obtained from a prototypical measurement.

20. The apparatus for modeling as in claim 9 wherein the means for measuring blood flow is a Doppler flow measurement device.

26. A system for modeling circulation within a living subject, comprising:

a computerized model of an arterial network made up of a plurality of arterial segments including one or more terminal efferent vessels, wherein the system includes means for calculating blood flows in the arterial network when the model is forced with a forcing function;

means for adapting the model of the arterial network to substantially conform to a specific arterial anatomy of the living subject;

wherein the arterial network is divided into a plurality of sectors, wherein each sector is terminated by a terminal efferent vessel and has a primary input vessel, and further wherein the terminal efferent vessels are modeled as terminal efferent resistances;

means for measuring blood flows in the living subject corresponding to the primary input flows of the sectors; and

means for adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow for each such sector.

27. The system of claim 26 further comprising means for iteratively calculating blood flows in the terminal efferent vessel of each sector and adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow for each such sector.

28. The system of claim 26 further comprising means for iteratively calculating blood flows in the terminal efferent vessel of each sector and adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow for each such sector until the differences between the calculated terminal efferent flows and the corresponding measured primary input flows for each sector are within specified limit values.

29. The system of claim 26 wherein the terminal efferent resistance of a sector is adjusted by multiplying the terminal efferent resistance by the ratio of the measured primary input flow to the calculated terminal efferent flow.

30. The system of claim 26 wherein one or more sectors has one or more secondary input or output vessels and further comprising:

means for calculating flows through the secondary input and output vessels; and

,

means for adjusting the terminal efferent resistances for each sector in a manner which tends to make the calculated flow in the terminal efferent vessel match the measured primary input flow adjusted for the calculated flows in the secondary input and output vessels.