Method and system for patient-specific modeling of blood flow

1. A method for non-invasive assessment of coronary artery stenosis, comprising:

extracting patient-specific anatomical measurements of coronary arteries of a patient from medical image data of the patient;

determining a model of coronary circulation comprising a reduced-order model that defines a blood pressure drop across at least one stenosis region of the coronary arteries of the patient as a function of blood flow velocity through the stenosis region;

determining first patient-specific boundary conditions of the model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient;

determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary conditions;

simulating blood flow and pressure across the at least one stenosis region of at least one coronary artery using the model of coronary circulation and the second patient-specific boundary conditions; and

determining a characteristic of blood flow through the at least one stenosis region using the simulated blood flow and pressure.

2. The method of claim 1 , wherein the reduced-order model comprises a one-dimensional computational model representing the coronary arteries and aorta of the patient.

3. The method of claim 1 , wherein the model of coronary circulation further comprises the reduced order model coupled to a one-dimensional computational model representing the coronary arteries.

4. The method of claim 3 ,

wherein the patient-specific anatomical measurements of the coronary arteries are extracted from medical image data that is acquired from the patient during a rest state;

wherein the first patient-specific boundary conditions are rest state boundary conditions;

wherein the second patient-specific boundary conditions are hyperemic boundary conditions further determined using a model of simulated hyperemia; and

wherein the simulated blood flow is hyperemic blood flow.

5. The method of claim 1 , wherein the model of coronary circulation further comprises a full-order three-dimensional computational model of the at least one stenosis region, coupled to the reduced-order model representing the at least one coronary artery.

6. The method of claim 1 , wherein the model of coronary circulation further comprises lumped models representing coronary microvascular beds, each coupled to a termination of a one-dimensional computational model representing a coronary artery branch.

7. The method of claim 1 , wherein determining first patient-specific boundary conditions of the model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient comprises:

determining a resting microvascular resistance at a termination of each of a plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient.

8. The method of claim 7 , wherein determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary comprises:

determining a hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance.

9. The method of claim 8 , wherein determining a resting microvascular resistance at the termination of each of the plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient comprises:

determining a mean arterial pressure (MAP) using a heart rate, a systolic blood pressure, and a diastolic blood pressure of the patient;

determining a resting myocardial perfusion using the heart rate and the systolic blood pressure of the patient;

determining a total resting coronary flow using the resting myocardial perfusion and a mass of the left ventricle of the patient; and

determining the resting microvascular resistance at the termination of each of the plurality of branches using the MAP and the total resting coronary flow.

10. The method of claim 9 , wherein the mass of the left ventricle is estimated from the medical image data of the patient.

11. The method of claim 8 , wherein determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance comprises:

determining a patient-specific total coronary resistance (TCR) using the non-invasive clinical measurements of the patient; and

determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined hyperemic microvascular resistance and the patient-specific TCR.

12. The method of claim 11 , wherein determining a patient-specific total coronary resistance (TCR) using the non-invasive clinical measurements of the patient comprises:

determining a resting average peak velocity using a heart rate and a systolic blood pressure of the patient;

determining a coronary flow reserve (CFR) value for each branch of the coronary arteries using a patient characteristic; and

determining the TCR for each branch using the CFR value and an estimated resting mean arterial pressure (MAP).

13. The method of claim 11 , wherein determining a patient-specific total coronary resistance (TCR) using the non-invasive clinical measurements of the patient comprises:

determining the TCR using a heart rate of the patient.

14. The method of claim 1 , wherein determining a characteristic of blood flow through the at least one stenosis region using the simulated blood flow and pressure comprises:

determining the fractional flow reserve (FFR) of the at least one stenosis region as a ratio of a mean simulated hyperemic pressure distal to the at least one stenosis and a mean simulated hyperemic aortic pressure over a cardiac cycle.

15. An apparatus comprising at least one computer system containing computer-executable programming instructions for performing a method for determining patient-specific cardiovascular information, the method comprising:

extracting patient-specific anatomical measurements of coronary arteries of a patient from medical image data of the patient;

determining a model of coronary circulation comprising a reduced-order model that defines a blood pressure drop across the at least one stenosis region as a function of blood flow velocity through the stenosis region;

determining first patient-specific boundary conditions of the model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient;

determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary conditions;

simulating blood flow and pressure across the at least one stenosis region of at least one coronary artery using the model of coronary circulation and the second patient-specific boundary conditions; and

determining a characteristic of blood flow through the at least one stenosis region using the simulated blood flow and pressure.

16. The apparatus of claim 15 , wherein the reduced-order model comprises a one-dimensional computational model representing the coronary arteries and aorta of the patient.

17. The apparatus of claim 15 , wherein determining first patient-specific boundary conditions of a model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient comprises:

determining a resting microvascular resistance at a termination of each of the plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient.

18. The apparatus of claim 17 , wherein determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary conditions comprises:

determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance.

19. The apparatus of claim 18 , wherein determining a resting microvascular resistance at the termination of each of the plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient comprises:

determining a mean arterial pressure (MAP) using a heart rate, a systolic blood pressure, and a diastolic blood pressure of the patient;

determining a resting myocardial perfusion using the heart rate and the systolic blood pressure of the patient;

determining a total resting coronary flow using the resting myocardial perfusion and a mass of the left ventricle of the patient; and

determining the resting microvascular resistance at the termination of each of the plurality of branches using the MAP and the total resting coronary flow.

20. The apparatus of claim 18 , wherein determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance comprises:

determining a patient-specific total coronary resistance (TCR) using the non-invasive clinical measurements of the patient; and

determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined hyperemic microvascular resistance and the patient-specific TCR.

21. A non-transitory computer readable medium storing computer program instructions for non-invasive assessment of coronary artery stenosis, the computer program instructions when executed by a processor cause the processor to perform operations comprising:

extracting patient-specific anatomical measurements of coronary arteries of a patient from medical image data of the patient;

determining a model of coronary circulation comprising a reduced-order model that defines a blood pressure drop across the at least one stenosis region as a function of blood flow velocity through the stenosis region;

determining first patient-specific boundary conditions of the model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient;

determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary conditions;

simulating blood flow and pressure across the at least one stenosis region of at least one coronary artery using the model of coronary circulation and the second patient-specific boundary conditions; and

determining a characteristic of blood flow through the at least one stenosis region using the simulated blood flow and pressure.

22. The non-transitory computer readable medium of claim 21 , wherein the reduced-order model comprises a one-dimensional computational model representing the coronary arteries and aorta of the patient.

23. The non-transitory computer readable medium of claim 21 , wherein determining first patient-specific boundary conditions of a model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient comprises:

determining a resting microvascular resistance at a termination of each of the plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient.

24. The non-transitory computer readable medium of claim 23 , wherein determining second patient-specific boundary conditions of the model of coronary circulation using the first patient-specific boundary comprises:

determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance.

25. The non-transitory computer readable medium of claim 24 , wherein determining a resting microvascular resistance at the termination of each of the plurality of branches using the patient-specific anatomical measurements and the non-invasive clinical measurements of the patient comprises:

determining a mean arterial pressure (MAP) using a heart rate, a systolic blood pressure, and a diastolic blood pressure of the patient;

determining a resting myocardial perfusion using the heart rate and the systolic blood pressure of the patient;

determining a total resting coronary flow using the resting myocardial perfusion and a mass of the left ventricle of the patient; and

determining the resting microvascular resistance at the termination of each of the plurality of branches using the MAP and the total resting coronary flow.

26. The non-transitory computer readable medium of claim 24 , wherein determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined resting microvascular resistance comprises:

determining a patient-specific total coronary resistance (TCR) using the non-invasive clinical measurements of the patient; and

determining the hyperemic microvascular resistance at the termination of each of the plurality of branches using the determined hyperemic microvascular resistance and the patient-specific TCR.

27. The non-transitory computer readable medium of claim 26 , wherein determining a patient-specific total coronary resistance index (TCR) using the non-invasive clinical measurements of the patient comprises:

determining a resting average peak velocity using a heart rate and a systolic blood pressure of the patient;

determining a coronary flow reserve (CFR) value for each branch of the coronary arteries using a patient characteristic; and

determining the TCR for each branch using the CFR value and an estimated resting mean arterial pressure (MAP).

28. The non-transitory computer readable medium of claim 26 , wherein determining a patient-specific total coronary resistance index (TCR) using the non-invasive clinical measurements of the patient comprises:

determining the TCR using a heart rate of the patient.

29. The non-transitory computer readable medium of claim 21 , wherein determining a characteristic of blood flow through the at least one stenosis region using the simulated blood flow and pressure comprises:

determining the fractional flow reserve (FFR) of the at least one stenosis region as a ratio of a mean simulated hyperemic pressure distal to the at least one stenosis and a mean simulated hyperemic aortic pressure over a cardiac cycle.

30. A method for non-invasive assessment of coronary artery stenosis, comprising:

extracting patient-specific anatomical measurements of the coronary arteries from medical image data of a patient acquired during a hyperemia state;

determining a model of coronary circulation comprises a reduced-order model that defines a blood pressure drop across the at least one stenosis region as a function of blood flow velocity through the stenosis region;

determining patient-specific hyperemic boundary conditions of the model of coronary circulation representing the coronary arteries using the patient-specific anatomical measurements and non-invasive clinical measurements of the patient at hyperemia;

simulating hyperemic blood flow and pressure across the at least one stenosis region of at least one coronary artery using the model of coronary circulation and the patient-specific hyperemic boundary conditions; and

determining a characteristic of blood flow through the at least one stenosis region using the simulated hyperemic blood flow and pressure.