[go: up one dir, main page]

WO2020236639A1 - Système et procédés d'estimation de flux sanguin à l'aide d'une surface de réponse et d'une modélisation d'ordre réduit - Google Patents

Système et procédés d'estimation de flux sanguin à l'aide d'une surface de réponse et d'une modélisation d'ordre réduit Download PDF

Info

Publication number
WO2020236639A1
WO2020236639A1 PCT/US2020/033239 US2020033239W WO2020236639A1 WO 2020236639 A1 WO2020236639 A1 WO 2020236639A1 US 2020033239 W US2020033239 W US 2020033239W WO 2020236639 A1 WO2020236639 A1 WO 2020236639A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood flow
vasculature
model
configurations
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/US2020/033239
Other languages
English (en)
Inventor
Sethuraman Sankaran
David Lesage
Charles Taylor
Nan XIAO
Hyun Jin Kim
David Spain
Michiel Schaap
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HeartFlow Inc
Original Assignee
HeartFlow Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by HeartFlow Inc filed Critical HeartFlow Inc
Priority to CN202080036393.5A priority Critical patent/CN113811956B/zh
Priority to AU2020278589A priority patent/AU2020278589A1/en
Priority to KR1020217040918A priority patent/KR20220009995A/ko
Priority to EP20731648.0A priority patent/EP3970162A1/fr
Priority to JP2021568016A priority patent/JP7616806B2/ja
Priority to CA3132472A priority patent/CA3132472A1/fr
Publication of WO2020236639A1 publication Critical patent/WO2020236639A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101—Computer-aided simulation of surgical operations
    • A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/107—Visualisation of planned trajectories or target regions
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Definitions

  • Various embodiments of the present disclosure relate generally to the prediction of the behavior of complex systems using a response surface and reduced order modeling, and, in particular, to efficient real-time estimation of blood flow using a response surface methodology and reduced order modeling.
  • Modeling and simulation of real-world physical phenomena may be performed to predict outcomes without invasive measurements.
  • many real-world physical phenomenma such as the flow of blood in arteries, fluid flow in porous media, and large deformation processes, may be modeled using partial differential equations. Modeling and simulation may also be used to design and optimize systems to yield a desired outcome.
  • blood flow characteristics may be relevant to assessing the health or disease of a patient.
  • hemodynamic indices may be used to assess the functional significance of lesions, blood perfusion levels, the transport of blood clots, the presence of aneurysms, and other health and disease characteristics. Hemodynamic indices may be measured invasively or assessed using blood flow simulation. While simulation techniques may be used to perform non- invasive assessments of the hemodynamics (based on available imaging data, for example), simulation techniques may also offer the potential benefit of predictive modeling of hemodynamics in response to various events (e.g., progression or regression of lesions) and predictive modeling of the outcome of planned procedures (e.g., surgical intervention). In order for predictive modeling to be realistic or clinically useful, it may be desirable or even necessary for modeling and simulation systems to be capable of computing results significantly faster than the average time needed to solve a high-fidelity model.
  • the present disclosure is, in various aspects, directed to addressing one or more of these above-referenced challenges.
  • the background description provided herein is for the purpose of generally presenting the context of the disclosure. Unless otherwise indicated herein, the materials described in this section are not prior art to the claims in this application and are not admitted to be prior art, or suggestions of the prior art, by inclusion in this section.
  • a computer-implemented method may include: performing a plurality of blood flow simulations using a first model of vascular blood flow, each of the plurality of blood flow simulations simulating blood flow in a vasculature of a patient or a geometry based on the vasculature of the patient; based on results of the plurality of blood flow simulations, generating a response surface mapping one or more first parameters of the first model to one or more second parameters of a reduced order model of vascular blood flow having lower fidelity than that of the first model;
  • a system may include a memory storing instructions; and one or more processors configured to execute the instructions to perform a method.
  • the method may include performing a plurality of blood flow simulations using a first model of vascular blood flow, each of the plurality of blood flow simulations simulating blood flow in a vasculature of a patient or a geometry based on the vasculature of the patient; based on results of the plurality of blood flow simulations, generating a response surface mapping one or more first parameters of the first model to one or more second parameters of a reduced order model of vascular blood flow having lower fidelity than that of the first model; determining values for the one or more parameters of the reduced order model mapped, by the response surface, from parameter values representing a modified state of the vasculature; and performing simulation of blood flow in the modified state of the vasculature using the reduced order model
  • a non-transitory computer-readable medium storing instructions that, when executed by one or more processors, cause the one or more processors to perform a method.
  • the method may include performing a plurality of blood flow simulations using a first model of vascular blood flow, each of the plurality of blood flow simulations simulating blood flow in a vasculature of a patient or a geometry based on the vasculature of the patient; based on results of the plurality of blood flow simulations, generating a response surface mapping one or more first parameters of the first model to one or more second parameters of a reduced order model of vascular blood flow having lower fidelity than that of the first model; determining values for the one or more parameters of the reduced order model mapped, by the response surface, from parameter values representing a modified state of the vasculature; and performing simulation of blood flow in the modified state of the vasculature using the reduced order model parameterized by the determined values for the one or more second parameters, to determine a blood flow characteristic of the modified
  • FIG. 1 depicts a flowchart of a method for estimating the behavior of a system using a response surface, according to one or more embodiments.
  • FIG. 2A illustrates a method of generating a response surface based on high-fidelity simulation, according to one or more embodiments.
  • FIG. 2B illustrates a method for predicting simulation results in real-time based on the response surface generated using the method illustrated by FIG. 2A, according to one or more embodiments.
  • FIG. 3 is a flowchart illustrating a method for modeling the effect of changing lumen geometry and boundary conditions on blood flow simulation, according to one or more embodiments.
  • FIG. 4 is a flowchart illustrating a method for modeling the effect of revascularization of coronary arteries, according to one or more embodiments.
  • FIGS. 5-6 illustrates an exemplary implementation of the method of FIG.
  • FIG. 7 illustrates an environment in which the a computer system for performing methods of the present disclosure may be implemented, according to one or more embodiments.
  • systems and methods permit a reduced order model, derived from computational fluid dynamics (CFD), to be used to simulate complex systems in real-time with arbitrary accuracy as compared to the accuracy of a high-fidelity model.
  • CFD computational fluid dynamics
  • the high-fidelity model of a physical system may be
  • the reduced order model may have a lower computational complexity than that of the high-fidelity model. Therefore, the reduced order model may be executed more quickly, so as to be more suitable for real-time simulation.
  • high-fidelity simulation using a high-fidelity model may be performed for a certain set of configurations.
  • the results of the high-fidelity simulation may then be used to parameterize a reduced order model.
  • the reduced order model may be parameterized using a response surface methodology according to the present disclosure.
  • the results of the high-fidelity simulation performed for the aforementioned set of configurations may be used to generate a response surface, which may be a mapping of parameters of the high-fidelity model to the reduced order model.
  • the response surface may then be used to parameterize the reduced order model.
  • Simulation using the parameterized reduced order model which may be real-time simulation, may be capable of predicting results significantly faster than high- fidelity simulation using the high-fidelity model, while achieving an accuracy that is arbitrarily close to that of the high-fidelity simulation.
  • the accuracy of the reduced order model and hence the accuracy of the simulation using the reduced order model, may depend on the set of configurations that was used to generate the response surface. Therefore, the accuracy of the reduced order model and the reduced order modeling may be tuned by increasing or otherwise adjusting the configurations used to generate the response surface. For example, by refining the response surface, it is possible to ensure that the simulation using the reduced order model has an accuracy within some error margin. Additionally, since the high-fidelity simulation used to generate the response surface may be computationally expensive, the high-fidelity simulation may be performed offline, prior to performing real-time simulation using the reduced order model.
  • the methods of the present disclosure may enable fast prediction of the behavior of a complex system, such as changes in hemodynamics in response to changes in the state of a patient.
  • changes in the state of a patient may be natural or planned (e.g., procedural).
  • the methods of the present disclosure may be used to produce real-time updates of FFRCT in response to a change in vessel lumen geometry.
  • This change in vessel lumen geometry may, for example, be a natural change, or a change that is expected to occur as a result of a candidate treatment.
  • the term“based on” means“based at least in part on.”
  • the singular forms“a,”“an,” and“the” include plural referents unless the context dictates otherwise.
  • the term“exemplary” is used in the sense of“example” rather than “ideal.”
  • the terms“comprises,”“comprising,”“includes,”“including,” or other variations thereof, are intended to cover a non-exclusive inclusion such that a process, method, or product that comprises a list of elements does not necessarily include only those elements, but may include other elements not expressly listed or inherent to such a process, method, article, or apparatus.
  • Relative terms, such as,“substantially” and “generally,” are used to indicate a possible variation of ⁇ 10% of a stated or understood value.
  • a reduced order model may also be referred to as a low- fidelity model or as a fast model.
  • a reduced order model that is usable for real-time simulation may also be referred to as a real-time model.
  • reduced order models and high fidelity models may be general models that can be parameterized using different parameter values, such as different values
  • parameter may refer to a parameter of any type, including boundary conditions.
  • L is an operator (e.g., a differential, integral, functional or a combination thereof)
  • u are unknowns
  • x N represents the problem dimensions
  • p represents given
  • Expressions (1 ) and (2) may represent a system, and may serve as a high-fidelity model of the system.
  • a reduced order model of the partial differential equation may approximate the operator L using a simpler operator (e.g., ordinary differential equations), reduce the dimensionality, x N , to an input space of the reduced order model, x n , in which
  • the reduced order model may be expressed as follows:
  • a goal is to have u(x n ) be a reasonable approximation to u(x w ), where x N may be a superset of x n .
  • a general approach is to perform simulations of the system, as originally formulated by expressions (1 ) and (2), for various boundary domains f B , boundary conditions &(. ), and parameter(s) such that a response surface may be used to generate an accurate approximation to the problem. Such simulations of the system may be referred to as high-fidelity simulations.
  • the domain f B may have bounds as expressed below:
  • parameter space may have bounds as expressed below:
  • M is the number of high-fidelity simulations performed.
  • M is the number of high-fidelity simulations performed.
  • Each of the M terms expressed above may correspond to as a configuration for which high fidelity simulation is to be performed. That is, the M terms may represents M configuration.
  • a“configuration” may refer to any modeling or simulation configuration, and may include any parameter (and its value).
  • a configuration may be set of value(s) of such parameter(s).
  • each of the M configurations may be represented as a set of values for the parameters of ⁇ (.), f B , and/or p.
  • the concept represented by a particular configuration may depend on the system that is being modeled. For example, if the system is blood flow through arteries of a patient, a configuration may represent a certain lumen geometry, a certain physiological state of the patient, or a combination thereof.
  • any suitable method such as a sampling method or a quadrature method, may be used in the selection of the M configurations.
  • the results of the high-fidelity simulation for the M configurations may be expressed as:
  • the response surface, R may be a mapping of the parameters of the high-fidelity model to the reduced order model:
  • R may capture the complexity of the original equations, enabling L to be a less complex operator than L.
  • the response surface R may be obtained by any suitable method. If R uses point-fitting polynomials, such as Lagrange polynomials, then the reduced order model may be constructed such that u(x n ) o u(x n ) at the M
  • the reduced order model may be built to exactly match the output of the high-fidelity model for the M configurations. This approach allows a computer to solve the faster problem of
  • a high-fidelity model may include any number of mathematical relationships. Accordingly, a high-fidelity model may include multiple different mathematical relationships of the form given by expression (1 ) described above, and may include other mathematical relationships. Similarly, a reduced order may have multiple mathematical relationships, and may have multiple different mathematical relationships of the form given by expression (3) described above.
  • a high- fidelity simulation may utilize all information available about the system in question (e.g., the full spatial and temporal representation), and the high-fidelity model used for the simulation may include any number of full-order governing equations.
  • a response surface such as response surface R, may be a mathematical relationship between a quantity or quantities of interest or parameters, and the underlying variables.
  • a response surface may be a function (e.g., a fitted function) that maps input variable(s) (e.g., parameters of a high-fidelity model) to output variables (e.g., parameters of a reduced order model).
  • the response surface may be built in a manner such that the response surface explores the parametric space using the reduced order model.
  • FIG. 1 is a flowchart illustrating a method for estimating the behavior of a system using a response surface, according to one or more embodiments.
  • Step 101 may include performing a plurality of simulations using a first model of a system.
  • the first model may be a high-fidelity model as described in this disclosure.
  • the first model may be a high-fidelity model of vascular blood flow
  • the simulations may be blood simulations that simulate blood flow in a vasculature of a patient or a vascular geometry based on the vasculature of the patient (e.g., a derived vasculature determined based on the vasculature of the patient).
  • the term“vasculature of a patient” may refer to vasculature in any portion of the body of the patient. Examples of vasculature include, but are not limited to, coronary
  • a derived vasculature may be, for example, a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having undergone a hypothetical vasculature having
  • Step 102 may include generating, based on the simulation results obtained from step 101 , a response surface mapping parameter(s) of the first model to parameter(s) of a second model having a lower fidelity than that of the first model.
  • the second model may be a model having a lower fidelity than that of the first model, such as a reduced order model as described in this disclosure. Since the first model and the second model may respectively be the high-fidelity model and the reduced order model, as described above, the response surface may be a mapping of the parameter(s) of the high-fidelity model to the parameter(s) of the reduced order model. This mapping may be a function whose output is values for the parameter(s) of the reduced order model and whose input is values of the parameter(s) of the high-fidelity model.
  • Step 103 may include determining values for the parameter(s) of the second model mapped, by the response surface, from parameter values of a
  • the parameter values of the configuration to be analyzed may be values for the aforementioned parameter(s) of the first model.
  • the first model may be a set of differential equations. Therefore, the parameter values of the configuration to be analyzed may be values of parameters (including boundary conditions) used in such differential equations.
  • the values for the second parameter(s) may be determined by the response surface as a function of the parameter values of the configuration to be analyzed.
  • Step 104 may include performing simulation using the second model parameterized by the determined values of the parameter(s) of the second model.
  • the parameter values for the configuration to be analyzed in step 103 may represent a modified state (e.g., a modified anatomical and/or physiological state) of the vasculature of the patient, in which case step 104 may determine a blood flow characteristic of the modified state of the vasculature.
  • the simulation may be performed in real-time.
  • the blood flow characteristic may be fractional flow reserve (FFR), flow magnitude, flow direction,
  • FIG. 2A illustrates a method of generating a response surface based on high-fidelity simulation.
  • the method of FIG. 2A illustrates an example implementation of the portion of the method of FIG. 1 corresponding to steps 101 and 102.
  • Step 201 may include receiving information indicative of configurations.
  • Information indicative of a configuration may include, for example, one or more geometries (e.g., a geometry in which fluid flow is to be modeled or simulated), one or more boundary conditions, and/or any other parameters that may be part of a configuration.
  • the information indicative of configurations may be indicative of a range of possible configurations, in which case the information received in step 201 may be indicative of a range of values for the aforementioned parameters
  • the information received in step 201 and may be manually input by a user or automatically determined by a process executed on the computer system.
  • Step 202 may include identifying configurations 220 for high fidelity simulation.
  • Configurations 220 may be identified based on the information received in step 201. For example, if the information received in step 201 is indicative of a range of configurations, the configurations 220 identified in step 202 may be a sample of configurations within the range of configurations. Examples of sampling and quadrature methods are discussed below in connection with the method of FIG. 3. Configurations 220 may be identified automatically, or identified based on user input.
  • Step 203 may include performing high-fidelity simulation for the identified configurations for high-fidelity simulation.
  • the configurations 220 identified in step 202 may be input into a high-fidelity model and high-fidelity simulation may be performed using the high-fidelity model parameterized in accordance with parameter values specified in the configurations 220.
  • Step 204 may include deriving parameters of a reduced order model.
  • the parameters derived in step 204 may be derived based on configurations 220 identified in step 202 and the results of the high-fidelity simulation performed using high-fidelity model.
  • Step 205 may include generating a response surface 224.
  • a response surface may be a mapping of the parameters of the high-fidelity model to the reduced order model.
  • the results of simulation using the high-fidelity model and the parameters of the reduced order model may define a correspondence between values of the parameters of the high-fidelity model and values of the parameters of the reduced order model.
  • Such correspondence may be represented as a set of points 222.
  • the response surface 224 may then be generated based on the set of points 222.
  • the response surface 224 may be a surface fitted to points 222.
  • the response surface may have an exact fit in that the surface 224 includes (intersects) all of points 222, as illustrated in FIG. 2A. Flowever, such is not a
  • Whether the surface 224 includes all of the points 222 may depend on the functional form of surface 224. As noted above, Lagrange polynomials may be used for an exact fit. In other fitting methods, it is possible for surface 224 to include only a portion of the points 223 or none of the points 223.
  • the parameters of the reduced order model derived in step 204 may be a set of parameter values that, when used in simulation using the reduced order model, yields the same results as reduced order model computes to the same results of the high-fidelity simulations. For example, if N configurations for the high fidelity simulation specified respective parameter values
  • Step 206 may include assessing an accuracy of the response surface 224.
  • Step 207 may determine, based on the accuracy assessed in step 206, whether response surface 224 is to be refined to have a higher accuracy.
  • the accuracy of a response surface may be defined by any suitable criteria.
  • accuracy may be a measure of accuracy in replicating results of the high-fidelity simulation.
  • the accuracy may be based on a closeness of results of reduced order modeling, when using the response surface 224 to parameterize the reduced order model for one or more testing configurations, to results of high-fidelity simulation for those one or more testing configurations.
  • the one or more testing configurations may include one or more configurations different from the configurations represented by the points 222 based upon which the response surface 224 was generated.
  • Step 207 may resolve in“YES” if the accuracy of the response surface 224 assessed in step 206 is insufficient (e.g., not satisfying a predefined threshold condition), and may resolve in“NO” if the accuracy of the response surface 224 is assessed as being sufficient (e.g., satisfying a predefined threshold).
  • accuracy may, for example, refer to accuracy of the reduced order model for arbitrarily defined configurations.
  • step 207 resolves in“YES” (e.g., accuracy is insufficient)
  • the method shown in FIG. 2A may proceed to step 208, which may include refining the
  • the process of refining the configurations may include adding new configurations for high fidelity simulation, removing existing configurations, and/or adjusting the values of existing configurations. For example, as shown in FIG. 2A, additional configurations may be added to the originally identified configurations 220 to improve accuracy of the response surface 224, so as to obtain a refined set of configurations 220A. Simulation using the high fidelity model (step 203) may be performed for any newly added configuration, such that the resulting response surface 224 is updated.
  • step 207 may implement a reiterated process in which the configurations for high fidelity simulation are refined (e.g., increased) in each
  • each configuration for which high-fidelity simulation is performed in step 203 may result in a corresponding point 222. Therefore, by adding additional configurations, the number of points 222 may be increased. Response surface 224 may then be fitted to the increased number of points 222 to as to potentially result in better accuracy.
  • step 207 may resolve in“NO,” and the response surface 224 may then be accepted as final response surface 224A.
  • the points 222A of final response surface 224A which may also be referred to as control points, may be more numerous than the points 222 of the initial response surface 224.
  • the final response surface 224A also serves as an example of the aforementioned response surface R, in which case the set of configurations that is used to generate the final response surface 224A would serve as an example of aforementioned set of M configurations.
  • the method illustrated in FIG. 2A may be computationally expensive, depending on the time it takes to perform high-fidelity simulation for all configurations for which high-fidelity simulation is performed. Accordingly, the method may be performed offline. For example, the final response surface 224A may be generated in advance of real-time simulations using the reduced order model.
  • FIG. 2B illustrates a method for predicting simulation results in real-time based on the response surface 224A generated using the method of FIG. 2A.
  • the method of FIG. 2B illustrates an example implementation of the portion of the method of FIG. 1 corresponding to steps 103 and 104.
  • Step 241 may including receiving a configuration to be analyzed.
  • the configuration may be defined by any suitable method.
  • the configuration may represent the settings of a certain experiment to be performed via reduced order simulation.
  • Step 242 may include probing the response surface.
  • the probing process may determine a value of a parameter (e.g., a parameter in parameter set p) of the reduced order model for the configuration to be analyzed.
  • the probing process is illustrated using point 250, which represents a value of a parameter for the reduced order model for the configuration to be analyzed.
  • point 250 may be a point that is mapped from the configuration to be analyzed.
  • the configuration to be analyzed may have values of the parameters p, b, and f B discussed above, and the response surface may determine the value of p as a function of those values of p, b, and ⁇ p B . That is, point 250 may have a position on response surface 224A
  • p, p, b, and p s may have may represents a position having the determined value of p. Since the positions on response surface 224A may be interpolated from the positions of points 224A, the position of point 250 may therefore be at an interpolated position.
  • Step 243 may include solving the reduced order model using the mapping given by the response surface 225.
  • Step 204 may include solving expression (11 ) as described above.
  • Steps 242 and 243 may be performed in real time as part of real-time simulation.
  • Step 244 may include generating and reporting the results of the simulation.
  • the results may be stored in an electronic storage device, or presented to a user (e.g., displayed on a display). Since the solving of the reduced order model may be a real-time process, the results of the reduced order model may also be presented in real-time.
  • the prediction of the behavior of a complex system may include a first process of generating the response surface, as described in relation to FIG. 2A, and a second process of fast probing of the response surface to estimate results (e.g., hemodynamic indices) for a particular configuration, as described above in relation to FIG. 2B.
  • the first process of building the response surface may be performed offline and be computationally expensive, depending on the time it takes to perform high-fidelity simulations using a high-fidelity model. The computational expense of may depend on the acceptable error for the second process.
  • FIGS. 3 and 4 illustrate further examples in which the techniques described above are applied.
  • FIG. 3 is a flowchart illustrating a method for modeling the effect of changing lumen geometry and boundary conditions on blood flow (e.g., coronary flow) simulation.
  • the method may apply various techniques described above to perform real-time estimation of blood flow in arteries (e.g., coronary arteries) under a given new configuration.
  • the given new configuration may be, for example, a lumen geometry and/or physiological state of a patient.
  • the method of FIG. 3 may be performed by any suitable computer system.
  • Step 301 may include receiving anatomical information describing a vasculature of a patient.
  • the described vasculature may include all arteries of the patient that are of interest.
  • the vasculature may be a coronary vasculature, in which case the anatomical information may describe the coronary arteries of the patient.
  • examples of other types of vasculature include, but are not limited to, peripheral vasculature, cerebral vasculature, renal vasculature, visceral vasculature, and hepatic vasculature such as portal veins.
  • the anatomical information may be received from a memory (e.g., a hard drive or other electronic storage device) of the computer system performing step 301 , or from another computer system (e.g., a computer system of a physician or third party provider) over a computer network.
  • a memory e.g., a hard drive or other electronic storage device
  • another computer system e.g., a computer system of a physician or third party provider
  • the anatomical information may include one or more images of the patient acquired using an imaging or scanning modality, and/or information extracted from (or otherwise obtained based on analysis of) such images of the patient.
  • imaging or scanning modalities include computed tomography (CT) scans, magnetic resonance (MR) imaging, micro-computed tomography (pCT) scans, micro-magnetic resonance (pMR) imaging, dual energy computed tomography scans, ultrasound imaging, single photon emission computed tomography (SPECT) scans, and positron emission tomography (PET) scans.
  • CT computed tomography
  • MR magnetic resonance
  • pCT micro-computed tomography
  • pMR micro-magnetic resonance
  • dual energy computed tomography scans dual energy computed tomography scans
  • ultrasound imaging single photon emission computed tomography (SPECT) scans
  • PET positron emission tomography
  • PET positron emission tomography
  • any model that is derived from or constructed based on such images, or other patient- specific information may be regarded as a patient-specific model. It is noted that use of the term“patient” is not intended to be limiting. A“patient” may generically be referred to as a“person.”
  • Step 302 may include generating an anatomical model of the vasculature, based on the anatomical information received in step 301.
  • the anatomical model may be of any suitable form and may model any suitable aspect of the vasculature.
  • the anatomical model may describe the patient-specific, three-dimensional geometry of the blood vessels of the vasculature as discerned from the anatomical information.
  • the anatomical model may indicate disease progression or regression, plaque rupture, thrombosis, and other characteristic of the represented vasculature(s).
  • An anatomical model of a vasculature may also be referred to as a patient-specific anatomical model or a patient-specific vascular model.
  • the anatomical model may model the characteristics of the vasculature under one or more physiological states of the patient.
  • the characteristics of the anatomical model may reflect characteristic of the vasculature when the patient is in a certain physiological state (e.g., a resting state or an exercising state).
  • a certain physiological state e.g., a resting state or an exercising state.
  • Steps 301 and 302 may be performed by the same computer system that performs the remaining steps 304 to 307 described below. Flowever, it is also possible for steps 301 and 302 to be performed by another computer system, in which case the anatomical model provided to the computer system that performs the remaining steps over a communication network. Any anatomical model received over a communication network may be stored in the memory of the computer system.
  • Step 303 may include performing high-fidelity simulation based on the anatomical model generated in step 302.
  • the simulation may be a blood flow simulation that simulates the flow of blood in the arteries as modeled by the anatomical model.
  • the high-fidelity simulation may involve detailed mathematical relationship(s) describing the system. Such mathematical relationships may include partial differential equation(s), such as the Navier-Stokes equations, in any suitable formulation.
  • the high-fidelity simulation may be performed using any suitable technique(s), such as finite element analysis, finite difference methods, lattice Boltzmann methods, etc.
  • the detailed mathematical relationship used in the high-fidelity simulation may constitute a high-fidelity model that is executed to perform the high-fidelity simulation.
  • the detailed mathematical relationships may include Navier- Stokes equations with boundary conditions and/or other parameters derived from the anatomical model.
  • the boundary conditions and/or other parameters may, for example, represent the geometry or other characteristics of the arteries as modeled by the anatomical model.
  • Step 304 may include performing high-fidelity simulation on extrema of configurations to be explored.
  • Such simulations may be a blood flow simulation that simulates the flow of blood in a structure represented by the extrema of the
  • the configurations to be explored may be any configuration that is intended to be explored (e.g., simulated or otherwise studied) using the reduced order simulation described below.
  • the extrema of the configurations to be explored may depend on the extrema of the parameter space and domain that is able to be explored using the reduced order model. Bounds may be imposed based on the limits of exploration. Such bounds may be application-specific. It is noted that the extrema of configurations to be explored, as described above, may be configurations for purposes of generating a response surface and may also be referred to as extrema
  • one or more bounds may be imposed based on anatomical limits.
  • an upper bound on the anatomical model may be imposed based on a maximally allowable dilation for a patient-specific model.
  • the patient-specific model may model the relieving of lumen narrowing at various locations, the effect of applying a higher level of nitrate, or a combination thereof.
  • a maximally allowable dilation in such treatment situations may be represented as an upper bound on the anatomical model.
  • one or more bounds imposed based on anatomical limits may represent the addition or removal of vessels.
  • an upper bound may be the maximum number of anastomoses based on available grafts.
  • one or more bounds may be imposed based on physiological limits.
  • an upper bound and/or a lower bound may be assessed based on resting-state and exercise conditions or based on other extrema of the boundary conditions.
  • an upper (or lower) bound may represent a restating state of the patient
  • an lower (or upper) bound may represent an exercising state of the patient.
  • Step 305 may include identifying one or more configurations for which high-fidelity simulation is to be performed, and performing high-fidelity simulation on the one or more identified configurations. It is noted that step 305 serves an example of steps 202 and 203 discussed above in connection with the method of FIG. 2A.
  • Any sampling or quadrature method may be used to identify the one or more configurations, including (but not restricted to): a Monte-Carlo sampling method, a Latin hypercube sampling method, a Gaussian quadrature method, a Sparse-grid quadrature method, an adaptive sparse-grid quadrature method, and combinations thereof.
  • Monte-Carlo sampling may be appropriate in sampling a large-dimensional parameter space but may converge very slowly for problems with moderate-dimensional parameter space.
  • Latin hypercube sampling may achieve separation of the parameter space and may converge better than Monte-Carlo for moderate-dimensional parameter space.
  • Gauss points may be used to generate the configurations and tensor-product interpolation may be used to scale the points to higher dimensions.
  • a sparse-grid quadrature method may be the same as the Gaussian quadrature method for one dimension but may have a sparser grid to reduce the number of simulations.
  • a adaptive sparse-grid quadrature method may be the same as the sparse-grid quadrature but may adapt to the function so that regions of shallow variations are explored less than regions of significant variations.
  • step 305 may further include performing high fidelity simulation on the one or more identified configurations.
  • Step 306 may include generating a response surface based on solutions of high-fidelity simulation.
  • the response surface may be created based on high-fidelity solutions at a plurality of configurations (e.g., M configurations) using any functional form.
  • M configurations e.g., M configurations
  • the M configurations referred to in the foregoing discussion may include any of the configurations identified in step 305, and may also include any configurations simulated in step 303 and/or step 304. Locally linear interpolation or Lagrange-polynomial interpolation may be performed to ensure that the solutions of the real-time simulation at the control points match the solution of the full simulation.
  • step 306 may utilize any of the techniques described above in connection with steps 204, 205, 206, and 208 of FIG. 2A.
  • multiple response surfaces may be generated. For example, if the high-fidelity simulation involved multiple mathematical relationship (e.g., a mathematical relationship in the form of expression (1 )) and/or the reduced order model includes multiple mathematical relationship (e.g., a mathematical relationship in the form of expression (3)), then multiple response surfaces may be generate to map between different combinations of high-fidelity and reduced order mathematical relationships. Furthermore, the response surface of step 306 may be revised by refining the configurations used to generate the response surface, as described above in connection with FIG. 1.
  • Step 307 may include performing reduced order simulation based on the response surface.
  • the reduced order simulation may be informed by interpolated values estimated by the response surface.
  • the response surface may be probed based on one or more configurations to be explored, to obtain the interpolated values.
  • the one or more configurations to be explored may depend on the application of the method.
  • the reduced order simulation may be performed using a reduced order model, which may be constructed to exactly match the output of the high-fidelity model for the M configurations where high-fidelity simulations were performed.
  • the reduced order simulation may be performed in real-time.
  • the method of FIG. 3 may include any one or more of the additional exemplary aspects described below, all of which are optional. These aspects may be implemented into one or more steps of the method described above, or implemented as additional steps of the method.
  • the method of FIG. 3 may include quantifying confidence intervals.
  • the response surface created in step 306 may be probed to run many simulations, from which confidence interval estimates for the unknown fields may be calculated.
  • the configurations to be explored may include any configuration suitable for quantifying the confidence intervals.
  • the configurations may be representative of configurations on which reduced order modeling is intended to be performed.
  • the confidence interval estimates may be, for example, used to assist a clinician performing the reduced order simulation in understanding the accuracy of the reduced order model in performing similar types of simulations.
  • the confidence interval estimates may be used to revise the response surface generated in step 306.
  • the method of FIG. 3 may include modeling disease progression and/or regression.
  • the response surface generated in step 306 may also be probed to predict the impact of lesions that might progress or regress. In turn, these may be used for patient management and monitoring.
  • the configurations to be explored may include any configurations suitable for modeling or simulating disease progression and/or regression.
  • the method of FIG. 3 may include the modeling of different physiological conditions.
  • the response surface generated in step 306 may also be probed to model different physiological conditions (e.g., resting and exercise conditions) or the effect of pharmacological agents.
  • the configurations to be explored may include any configurations suitable for modeling or simulating physical conditions.
  • FIG. 4 is a flowchart illustrating a method for modeling the effect of revascularization of coronary arteries. The method may apply various techniques described above to perform real-time computational of the effect of revascularization of coronary arteries on blood flow. An exemplary implementations of the method of FIG. 4 is illustrated by FIGS. 5-6, also discussed below. The method of FIGS. 4-6 may be performed by any suitable computer system.
  • Step 401 may include receiving anatomical information describing coronary arteries of a patient.
  • Step 401 may include any of the aspects of step 301 described above.
  • step 401 may include receiving anatomic information obtained from analysis of coronary CT scans. For example, as shown in FIG. 5, anatomical information describing anatomical characteristic of a patient, such as vessel centerlines and lumens, may be extracted from CCTA images 502 taken of the patient.
  • Step 402 may include generating patient models, which may include a base patient model and a modified patient model.
  • the base patient model may be generated based on the anatomical information received in step 401.
  • the modified patient model may be a modification of the base patient model.
  • the coronary arteries of the patient may have narrowed lumens, and the modified patient model may represent a full revascularization of the coronary arteries.
  • the base patient model may be an anatomical model that model actual anatomical characteristics (e.g., vessel centerlines and lumens) of the patient’s coronary arteries, as described by the anatomical information.
  • base patient model 503 may be generated based on vessel centerlines and lumens extracted from the CCTA images 502.
  • base patient model 503 is illustrated as having a narrowed geometry at various locations 503A, 503B and 503C of the model.
  • the narrowed geometry may model, for example, stenosis at corresponding locations of the patient’s coronary arteries.
  • the modified patient model may be the base patient model having been modified to model a change in characteristics of the patient’s coronary arteries.
  • the modified patient model may model a hypothetical condition of the patient’s coronary arteries.
  • Such condition may be, for example, an idealized condition corresponding to an extrema of configurations to be explored, in which case the modified patient model may be referred to as an idealized model.
  • idealized model 504 is an example of a modified patient model that models the coronary arteries under a condition in which the entire anatomy represented by the base patient model is revascularized.
  • the stenosis at locations 503A, 503B and 503C of the base patient model 503 are not indicated by the idealized model 504.
  • the modified patient model may be a revascularized anatomical model.
  • Step 403 may include performing high-fidelity simulation of blood flow using boundary conditions obtained based on the patient model, to simulate an effect of adenosine for hyperemia and obtain a first high-fidelity solution.
  • the boundary conditions in step 403 may be boundary conditions derived from
  • boundary conditions may include boundary conditions be obtained based (e.g., derived from) the patient model. Flowever, the present disclosure is not limited thereto, and it is also possible for some or all of the boundary conditions to be derived from other models or information.
  • the high-fidelity simulation of step 403 may be performed by constructing a computational model in the form of a high-fidelity model.
  • the computational model may include mathematical relationships, such as Navier-Stokes equations and boundary conditions derived from the patient’s anatomy, myocardium, scaling laws for resting blood flow. Such boundary conditions may simulate the effect of adenosine for hyperemia. Therefore, to perform the high-fidelity simulation, the computer system performing step 403 may solve the Navier-Stokes equations on the coronary arteries using the aforementioned boundary conditions.
  • Step 404 may include performing high-fidelity simulation of blood flow using boundary conditions obtained based on the modified patient model, to obtain a second high-fidelity solution.
  • the high-fidelity simulation of step 404 may be performed on extrema corresponding to configuration(s) in which the coronary arteries are fully revascularized. Such extrema serves as an example of the extrema of configurations to be explored described above in connection with step 304 of FIG. 3.
  • the revascularization may be one in which the entire anatomy of the patient-specific geometry is revascularized.
  • the high-fidelity simulation of step 404 may be performed using a computational model constructed as a high-fidelity model, which may include boundary conditions derived from the modified patient model described above.
  • Step 405 may include, for each of the base patient model and the modified patient model, performing an additional high-fidelity simulation of blood flow at a different flow rate, to obtain third and fourth high-fidelity solutions.
  • the high-fidelity simulation performed in step 402 may be performed for a first flow rate
  • step 405 may include a high-fidelity simulation performed in the same or substantially same manner (e.g., using boundary conditions based on the base patient model), but with a flow rate that is higher (e.g., 10%, 15%, 25%, 50%, or 75% higher) than the aforementioned first flow rate.
  • the high-fidelity simulation performed in step 403 may be performed for a second flow rate (which may be the same as the first flow rate), and step 405 may include a high-fidelity simulation that is performed in the same or substantially same manner (e.g., using boundary conditions based on the modified patient model) but with a flow rate that is higher (e.g., 10%, 15%, 25%, 50%, or 75% higher) than the aforementioned second flow rate.
  • the high-fidelity solutions obtained across steps 403-405 may be used to inform a reduced order model in which fluid resistance parameters depend on flow rate. With one additional simulation in each of the configurations associated respectively with the base patient model and the modified patient model, the fluid resistances in the reduced order model may depend linearly on flowrate.
  • Step 406 may include generating response surfaces respectively for the intercept and slope of the fluid resistance function.
  • the response surfaces may be generated based on the four high-fidelity solutions obtained across steps 403-405, and may include two response surfaces, a first response surface for the intercept of the fluid resistance function, and a second response surface for the slope of the fluid resistance function.
  • the first and second response surfaces may both be based on the one dimensional Navier-Stokes equations.
  • the first response surface may have the functional form 1/r 4 for the intercept.
  • the second response surface may have the functional form (dA/dz * 1/r 6 ) for the slope.
  • r is the local radius
  • A is the area
  • dA/dz is the gradient of area along the vessel.
  • Step 407 may include receiving a modified geometry.
  • the modified geometry may be a geometry that is to be subject to reduced order simulation, and may be a revascularized geometry including of for example, locations at which the coronary arteries are to be revascularized and final size(s) of vessel lumen(s).
  • revascularized geometry may be a simulation input that is defined by user input, or by a simulation process.
  • One or more configurations for reduced order modeling and simulation may be defined based on the revascularized geometry. For example, a value of an attribute of the vascularized geometry, such as a value of a location of revascularization and/or a value of the final size of a vessel lumen, may serve as a configuration or part of the configuration.
  • Such configurations may be used on the response surface(s) to obtain parameter of the reduced order model(s) that is used in step 408 described below.
  • Step 408 may include performing reduced order simulation based on the revascularized geometry and the two response surfaces.
  • the reduced-order simulation may be informed by interpolated values estimated using the response surfaces on the revascularized geometry.
  • the reduced order simulation may be performed in real-time, and may use one or more reduced order models constructed as described above.
  • Such reduced order model may have mathematical relationships in the form of expressions (3) and (4), and may be constructed such to yield the same results as the high fidelity model of process 510 for the four high-fidelity. It is noted that step 408 is an example of step 307 described above. Therefore, techniques described in connection with step 307 are generally applicable to step 408.
  • the output of the low fidelity simulation may be used to output the updated flowrates, blood pressures, FFR or any other quantity of interest, such as wall shear stress, for the configuration in step 403.
  • process 510 serves as an example of the high- fidelity simulations of steps 403 to 405.
  • four Navier-Stokes simulations may be performed. These simulations may include a first Navier-Stokes simulation using hyperemic boundary conditions applied based on idealized model 504, a second Navier-Stokes simulation using superemic boundary conditions applied based on idealized model 504, a third Navier-Stokes simulation using superemic boundary conditions applied based on the base patient model 503, and a fourth Navier-Stokes simulation using hyperemic boundary conditions applied based on base patient model 503. It is noted that the aforementioned boundary conditions serve as examples of simulation parameters, and that the respective simulation parameters of four
  • simulations may differ from one another in aspects other than the aforementioned boundary conditions.
  • the four sets of simulation parameters applied to the Navier-Stokes simulation may respectively result in four high-fidelity solutions, as described above in connection with step 406 of FIG. 4.
  • the four high-fidelity solutions may then be used to build response surfaces (step 520), the process of which may include deriving parameters of the reduced order model.
  • Item 504 in FIG. 4 is a visual depiction of a parameters of a reduced order model.
  • the reduced order model may be a reduced order model having mathematical relationships in the form of expressions (3) and (4), and may be constructed such that the reduced order model yields the same results as the high fidelity model of process 510 for the four high-fidelity simulations.
  • FIG. 6 illustrates the probing of the responses surfaces for reduced order modeling.
  • the response surface(s) may be probed based on configurations indicated by a modified geometry 601.
  • Modified geometry 601 may be a revascularized geometry as described above for step 407, and may be representable in a graphical form, such as a three-dimensional graphical model (e.g., a surface mesh), as shown in FIG. 6.
  • Modified geometry 601 may be an anatomical model, and may represent a particular anatomical geometry to be explored or analyzed by simulation; this geometry may, for example, be a state of the patient that is either natural or planned.
  • Modified geometry 601 may differ from the idealized 505.
  • the probing of the response surface(s) may obtain values for the parameters of the reduced order model.
  • the reduced order model may be executed to obtain a hemodynamic solution.
  • the hemodynamic solution may be graphically displayed along with the three-dimensional graphical model of the modified geometry 601.
  • the hemodynamic solution may be represented as in graphical form, and the graphics of the hemodynamic solution may be overlaid or otherwise combined with the three-dimensional graphical model of the modified geometry 610, to obtain a mapped model 602.
  • the mapped model 602 may be, for example, displayed on an electronic display. Such display may be performed in real time.
  • the methods described in this disclosure may have various clinical applications, include: planning a percutaneous coronary intervention (PCI) procedure; planning bypass graft surgery; modeling disease progression and regression of lesions; modeling positive and negative remodeling of lesions; sensitivity analysis, uncertainty quantification and/or estimation of confidence intervals for flow simulations; modeling of different physiologic conditions, such as exercise; modeling the effect of drugs, altitude or autoregulatory mechanisms.
  • PCI percutaneous coronary intervention
  • the methods described in this disclosure may be used to produce real-time updates of fractional flow reserve (FFR) (e.g., fractional flow reserve derived from computed tomography (FFRCT)) in response to a change in the vessel lumen geometry of a patient.
  • FFR fractional flow reserve
  • This change in vessel lumen geometry may be a natural change, or a change that is expected to occur as a result of a candidate treatment for a patient.
  • the lumen geometry may be represented as one or more parameters, and a user or simulation process may adjust the values of such parameters to reflect the change in vessel lumen geometry.
  • the computer system performing the simulation may identify configurations for reduced order modeling, probe response surface(s) based on the configurations to parameterize a reduced order model, and solve the reduced order model to compute value(s) of FFRCT.
  • the response surface(s) may have been generated prior to the simulation, in accordance with the methods described in this disclosure (e.g., FIGS. 3 and 4).
  • the computed value(s) of FFRCT may be output in any suitable manner (e.g., displayed on a display device, or transmitted to another computer system for display on a display device).
  • the vessel lumen geometry may be part of the coronary arteries of the patient, or part of another vasculature portion.
  • computation described in connection with expressions (1 ) to (11 ), may be performed by one or more processors of a computer system.
  • a step of a method performed by one or more processors may also be referred to as an operation.
  • FIG. 7 depicts an example of an environment in which such a computer system may be implemented as server systems 740.
  • the environment of FIG. 7 further includes a plurality of physicians 720 and third party providers 730, any of which may be connected to an electronic network 710, such as the Internet, through one or more computers, servers, and/or handheld mobile devices.
  • physicians 720 and third party providers 730 may each represent a computer system, as well as an organization that uses such a system.
  • a physician 720 may be a hospital or a computer system of a hospital.
  • Physicians 720 and/or third party providers 730 may create or otherwise obtain medical images, such as images of the cardiac, vascular, and/or organ systems, of one or more patients. Physicians 720 and/or third party providers 730 may also obtain any combination of patient-specific information, such as age, medical history, blood pressure, blood viscosity, the anatomical information described above in connection with step 301 of the method of FIG. 3, and other types of patient-specific information. Physicians 720 and/or third party providers 730 may transmit the patient- specific information to server systems 740 over the electronic network 710.
  • Server systems 740 may include one or more storage devices 760 for storing images and data received from physicians 720 and/or third party providers 730.
  • the storage devices 760 may be considered to be components of the memory of the server systems 740.
  • Server systems 740 may also include one or more processing devices 750 for processing images and data stored in the storage devices and for performing any computer-implementable process described in this disclosure.
  • Each of the processing devices 750 may be a processor or a device that include at least one processor.
  • server systems 740 may have a cloud computing platform with scalable resources for computations and/or data storage, and may run an application for performing methods described in this disclosure on the cloud computing platform.
  • any outputs may be transmitted to another computer system, such as a personal computer, for display and/or storage.
  • FIG. 1 Other examples of computer systems for performing methods of this disclosure include desktop computers, laptop computers, and mobile computing devices such as tablets and smartphones.
  • the one or more processors may be configured to perform such processes by having access to instructions (e.g., software or computer-readable code) that, when executed by the one or more processors, cause the one or more processors to perform the processes.
  • the instructions may be stored in a memory of the computer system.
  • a processor may be a central processing unit (CPU), a graphics processing unit (GPU), or another type of processing unit.
  • a computer system such as server systems 740, may include one or more computing devices. If the one or more processors of the computer system is implemented as a plurality of processors, the plurality of processors may be included in a single computing device or distribute among a plurality of computing devices. If a computer system comprises a plurality of computing devices, the memory of the computer system may include the respective memory of each computing device of the plurality of computing devices.
  • a computing device may include processor(s) (e.g., CPU, GPU, etc.).
  • processor(s) e.g., CPU, GPU, etc.
  • Memory may include volatile memory, such as RAM, and/or non-volatile memory, such as ROM and storage media.
  • volatile memory such as RAM
  • non-volatile memory such as ROM and storage media.
  • storage media include solid-state storage media (e.g., solid state drives and/or removable flash memory), optical storage media (e.g., optical discs), and/or magnetic storage media (e.g., hard disk drives).
  • the aforementioned instructions e.g., software or computer-readable code
  • the computing device may, in some embodiments, further include input device(s) (e.g., a keyboard, mouse, or touchscreen) and output device(s) (e.g., a display, printer).
  • input device(s) e.g., a keyboard, mouse, or touchscreen
  • output device(s) e.g., a display, printer
  • the aforementioned elements of the computing device may be connected to one another through a bus, which represents one or more busses.
  • the processor(s) of the computing device includes both a CPU and a GPU.
  • Instructions executable by one or more processors may be stored on a non-transitory computer-readable medium.
  • non-transitory computer-readable medium storing instructions that, when executed by one or more processors, configure or cause the one or more processors to perform the computer-implemented method.
  • Examples of non- transitory computer-readable medium include RAM, ROM, solid-state storage media (e.g., solid state drives), optical storage media (e.g., optical discs), and magnetic storage media (e.g., hard disk drives).
  • a non-transitory computer-readable medium may be part of the memory of a computer system or separate from any computer system.
  • An“electronic storage device” may include any of the non-transitory computer- readable media described above.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Mathematical Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Radiology & Medical Imaging (AREA)
  • Robotics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Physiology (AREA)
  • Operations Research (AREA)
  • Biophysics (AREA)
  • Software Systems (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)

Abstract

L'invention concerne des systèmes et des procédés de simulation de flux sanguin. Par exemple, un procédé peut comprendre la réalisation d'une pluralité de simulations de flux sanguin à l'aide d'un premier modèle de flux sanguin vasculaire, chacune de la pluralité de simulations de flux sanguin simulant un flux sanguin dans un système vasculaire d'un patient ou une géométrie basée sur le système vasculaire du patient ; sur la base des résultats de la pluralité de simulations de flux sanguin, la génération d'une surface de réponse mettant en correspondance un ou plusieurs premiers paramètres du premier modèle avec un ou plusieurs seconds paramètres d'un modèle d'ordre réduit de sang vasculaire ; la détermination des valeurs pour le ou les paramètres du modèle d'ordre réduit mis en correspondance, par la surface de réponse, à partir de valeurs de paramètre représentant un état modifié du système vasculaire ; et la réalisation d'une simulation à l'aide du modèle d'ordre réduit paramétré par les valeurs déterminées, pour déterminer une caractéristique de flux sanguin de l'état modifié du système vasculaire.
PCT/US2020/033239 2019-05-17 2020-05-15 Système et procédés d'estimation de flux sanguin à l'aide d'une surface de réponse et d'une modélisation d'ordre réduit Ceased WO2020236639A1 (fr)

Priority Applications (6)

Application Number Priority Date Filing Date Title
CN202080036393.5A CN113811956B (zh) 2019-05-17 2020-05-15 用于使用响应表面和降阶建模来估计血液流动的系统和方法
AU2020278589A AU2020278589A1 (en) 2019-05-17 2020-05-15 System and methods for estimation of blood flow using response surface and reduced order modeling
KR1020217040918A KR20220009995A (ko) 2019-05-17 2020-05-15 반응 표면 및 차수 감소 모델링을 사용한 혈류 추정 시스템 및 방법
EP20731648.0A EP3970162A1 (fr) 2019-05-17 2020-05-15 Système et procédés d'estimation de flux sanguin à l'aide d'une surface de réponse et d'une modélisation d'ordre réduit
JP2021568016A JP7616806B2 (ja) 2019-05-17 2020-05-15 応答曲面および減次モデル化を使用した血流の推定のためのシステムおよび方法
CA3132472A CA3132472A1 (fr) 2019-05-17 2020-05-15 Systeme et procedes d'estimation de flux sanguin a l'aide d'une surface de reponse et d'une modelisation d'ordre reduit

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201962849489P 2019-05-17 2019-05-17
US62/849,489 2019-05-17

Publications (1)

Publication Number Publication Date
WO2020236639A1 true WO2020236639A1 (fr) 2020-11-26

Family

ID=71069955

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2020/033239 Ceased WO2020236639A1 (fr) 2019-05-17 2020-05-15 Système et procédés d'estimation de flux sanguin à l'aide d'une surface de réponse et d'une modélisation d'ordre réduit

Country Status (8)

Country Link
US (3) US11707325B2 (fr)
EP (1) EP3970162A1 (fr)
JP (1) JP7616806B2 (fr)
KR (1) KR20220009995A (fr)
CN (1) CN113811956B (fr)
AU (1) AU2020278589A1 (fr)
CA (1) CA3132472A1 (fr)
WO (1) WO2020236639A1 (fr)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018515167A (ja) * 2015-04-02 2018-06-14 ハートフロー， インコーポレイテッド 血管網と潅流組織の機能的関連性を特定し可視化するシステム及び方法
WO2021141921A1 (fr) 2020-01-07 2021-07-15 Cleerly, Inc. Systèmes, procédés et dispositifs d'analyse d'images médicales, de diagnostic, de stratification de risque, de prise de décision et/ou de suivi de maladie
US11969280B2 (en) 2020-01-07 2024-04-30 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US20220392065A1 (en) 2020-01-07 2022-12-08 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US20250143657A1 (en) 2022-03-10 2025-05-08 Cleerly, Inc. Systems, devices, and methods for non-invasive image-based plaque analysis and risk determination
US12406365B2 (en) 2022-03-10 2025-09-02 Cleerly, Inc. Systems, devices, and methods for non-invasive image-based plaque analysis and risk determination
US12440180B2 (en) 2022-03-10 2025-10-14 Cleerly, Inc. Systems, devices, and methods for non-invasive image-based plaque analysis and risk determination
US20250217981A1 (en) 2022-03-10 2025-07-03 Cleerly, Inc. Systems, methods, and devices for image-based plaque analysis and risk determination
CN115221804B (zh) * 2022-06-08 2025-07-11 中国科学院深圳先进技术研究院 人工肺血流的流场计算方法、装置、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120041739A1 (en) 2010-08-12 2012-02-16 Heartflow, Inc. Method and System for Patient-Specific Modeling of Blood Flow
US20160125154A1 (en) * 2014-11-04 2016-05-05 Heartflow, Inc. Systems and methods for simulation of occluded arteries and optimization of occlusion-based treatments
US20170018081A1 (en) * 2015-07-17 2017-01-19 Heartflow, Inc. Systems and methods for assessing the severity of plaque and/or stenotic lesions using contrast distribution predictions and measurements

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19814371A1 (de) * 1998-03-31 1999-10-14 Pulsion Verwaltungs Gmbh & Co Verfahren zur in-vivo Bestimmung der Compliance-Funktion und des systemischen Blutflusses eines Lebewesens und Vorrichtung zur Durchführung der Verfahren
US8200466B2 (en) * 2008-07-21 2012-06-12 The Board Of Trustees Of The Leland Stanford Junior University Method for tuning patient-specific cardiovascular simulations
US10328266B2 (en) * 2012-05-31 2019-06-25 Zoll Medical Corporation External pacing device with discomfort management
US10398386B2 (en) * 2012-09-12 2019-09-03 Heartflow, Inc. Systems and methods for estimating blood flow characteristics from vessel geometry and physiology
US9424395B2 (en) * 2013-03-04 2016-08-23 Heartflow, Inc. Method and system for sensitivity analysis in modeling blood flow characteristics
US9805463B2 (en) * 2013-08-27 2017-10-31 Heartflow, Inc. Systems and methods for predicting location, onset, and/or change of coronary lesions
US20150164451A1 (en) * 2013-12-18 2015-06-18 Heartflow, Inc. Systems and methods for predicting coronary plaque vulnerability from patient-specific anatomic image data
US9501622B2 (en) * 2014-03-05 2016-11-22 Heartflow, Inc. Methods and systems for predicting sensitivity of blood flow calculations to changes in anatomical geometry
US9390232B2 (en) * 2014-03-24 2016-07-12 Heartflow, Inc. Systems and methods for modeling changes in patient-specific blood vessel geometry and boundary conditions
US20160196384A1 (en) * 2015-01-06 2016-07-07 Siemens Aktiengesellschaft Personalized whole-body circulation in medical imaging
US10517678B2 (en) * 2015-10-02 2019-12-31 Heartflow, Inc. System and method for diagnosis and assessment of cardiovascular disease by comparing arterial supply capacity to end-organ demand
CN105326486B (zh) * 2015-12-08 2017-08-25 博动医学影像科技（上海）有限公司 血管压力差与血流储备分数的计算方法及系统
WO2017147500A1 (fr) * 2016-02-26 2017-08-31 Heartflow, Inc. Systèmes et procédés d'identification et de modélisation de vaisseaux non résolus dans des modèles hémodynamiques spécifiques au patient basés sur des images
US10945606B2 (en) * 2016-09-20 2021-03-16 Heartflow, Inc. System and methods for estimation of blood flow characteristics using reduced order model and machine learning
US20220107256A1 (en) * 2019-02-06 2022-04-07 Universität Zürich Method and apparatus for predicting fluid flow through a subject conduit

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120041739A1 (en) 2010-08-12 2012-02-16 Heartflow, Inc. Method and System for Patient-Specific Modeling of Blood Flow
US20160125154A1 (en) * 2014-11-04 2016-05-05 Heartflow, Inc. Systems and methods for simulation of occluded arteries and optimization of occlusion-based treatments
US20170018081A1 (en) * 2015-07-17 2017-01-19 Heartflow, Inc. Systems and methods for assessing the severity of plaque and/or stenotic lesions using contrast distribution predictions and measurements

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GIANLUCA GERACI ET AL: "Towards a Multi-fidelity Hemodynamic Model Pipeline for the Analysis of Cardiovascular Flow Under Uncertainty", CONFERENCE PROCEEDINGS ARTICLE, 1 July 2017 (2017-07-01), United States, XP055720899 *
PARIS PERDIKARIS ET AL: "Model inversion via multi-fidelity Bayesian optimization: a new paradigm for parameter estimation in haemodynamics, and beyond", JOURNAL OF THE ROYAL SOCIETY. INTERFACE, vol. 13, no. 118, 1 May 2016 (2016-05-01), GB, pages 20151107, XP055720882, ISSN: 1742-5689, DOI: 10.1098/rsif.2015.1107 *
PETER BENNER ET AL: "A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems", SIAM REVIEW., vol. 57, no. 4, 1 January 2015 (2015-01-01), pages 483 - 531, XP055720890, ISSN: 0036-1445, DOI: 10.1137/130932715 *

Also Published As

Publication number Publication date
AU2020278589A1 (en) 2022-01-06
CN113811956A (zh) 2021-12-17
EP3970162A1 (fr) 2022-03-23
US20240407849A1 (en) 2024-12-12
CA3132472A1 (fr) 2020-11-26
US20230310085A1 (en) 2023-10-05
US11707325B2 (en) 2023-07-25
JP7616806B2 (ja) 2025-01-17
US12096990B2 (en) 2024-09-24
CN113811956B (zh) 2025-08-19
US20200360088A1 (en) 2020-11-19
JP2022533345A (ja) 2022-07-22
KR20220009995A (ko) 2022-01-25

Similar Documents

Publication Publication Date Title
US12096990B2 (en) Systems and methods for estimation of blood flow using response surface and reduced order modeling
AU2020244586B2 (en) Method and system for sensitivity analysis in modeling blood flow characteristics
JP7241790B2 (ja) 血管網と潅流組織の機能的関連性を特定し可視化するシステム及び方法
US20210161384A1 (en) System and methods for estimation of blood flow characteristics using reduced order model and machine learning
EP3117771B1 (fr) Calcul direct de biomarqueurs dérivés d'images
US11127138B2 (en) Automatic detection and quantification of the aorta from medical images
US20100130878A1 (en) Systems, apparatus and processes for automated blood flow assessment of vasculature
Colebank et al. Influence of image segmentation on one-dimensional fluid dynamics predictions in the mouse pulmonary arteries
Maul et al. Transient hemodynamics prediction using an efficient octree-based deep learning model
US20170004278A1 (en) System and method for sparse pressure/flowrate reduced modeling of hemodynamics
Tanade et al. HarVI: Real-time intervention planning for coronary artery disease using machine learning
JP2022519876A (ja) 対象となる導管を通る流体の流れを決定するシステムおよび方法
Eulzer et al. Instantaneous Visual Analysis of Blood Flow in Stenoses Using Morphological Similarity
Geddes et al. Digital twins for noninvasively measuring predictive markers of right heart failure
CN120015335A (zh) 一种体位性脑血流变化的评估方法、系统、终端及存储介质
JP2024155812A (ja) 血管ツインを用いた血管挙動の予測
Sainz De Mena et al. New tools and methods to inform patient-specific cancer models

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20731648

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3132472

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2021568016

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 20217040918

Country of ref document: KR

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2020731648

Country of ref document: EP

Effective date: 20211217

ENP Entry into the national phase

Ref document number: 2020278589

Country of ref document: AU

Date of ref document: 20200515

Kind code of ref document: A

WWG Wipo information: grant in national office

Ref document number: 202080036393.5

Country of ref document: CN