CN113674865B - Cardiovascular system data platform built-in processing method based on transient characteristics - Google Patents

Cardiovascular system data platform built-in processing method based on transient characteristics Download PDF

Info

Publication number
CN113674865B
CN113674865B CN202110987288.6A CN202110987288A CN113674865B CN 113674865 B CN113674865 B CN 113674865B CN 202110987288 A CN202110987288 A CN 202110987288A CN 113674865 B CN113674865 B CN 113674865B
Authority
CN
China
Prior art keywords
cardiovascular system
equation
data platform
system data
blood
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.)
Active
Application number
CN202110987288.6A
Other languages
Chinese (zh)
Other versions
CN113674865A (en
Inventor
周聪
周领
郭丽
李赟杰
吴金远
王宁
胡垠盈
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.)
FIRST AFFILIATED HOSPITAL OF JIAMUSI UNIVERSITY
Hohai University HHU
Original Assignee
FIRST AFFILIATED HOSPITAL OF JIAMUSI UNIVERSITY
Hohai University HHU
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 FIRST AFFILIATED HOSPITAL OF JIAMUSI UNIVERSITY, Hohai University HHU filed Critical FIRST AFFILIATED HOSPITAL OF JIAMUSI UNIVERSITY
Priority to CN202110987288.6A priority Critical patent/CN113674865B/en
Publication of CN113674865A publication Critical patent/CN113674865A/en
Application granted granted Critical
Publication of CN113674865B publication Critical patent/CN113674865B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Pathology (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • External Artificial Organs (AREA)

Abstract

The invention discloses a cardiovascular system data platform built-in processing method based on transient characteristics, which comprises the following steps: constructing a cardiovascular system blood flow transient equation considering dynamic friction and wall viscoelasticity; establishing a calculation grid, converting the hyperbolic partial differential equation set into a normal differential equation set according to a characteristic line method, and integrating to obtain a characteristic line equation; solving partial differentiation of the hysteresis strain of the vascular wall to time; introducing a Kagawa dynamic friction model and partial differentiation into a characteristic line equation; the calculation result is obtained and processed in consideration of the boundary conditions. The invention adopts the one-dimensional model to replace the existing three-dimensional CFD for modeling, breaks through the technical bottleneck that the existing cardiovascular system data platform is difficult to model by adopting the one-dimensional model, not only has simple modeling, strong universality and short simulation time, but also considers the dynamic friction effect and the viscoelastic effect of the vascular wall, is more in line with the actual situation, and improves the simulation precision of the platform.

Description

Cardiovascular system data platform built-in processing method based on transient characteristics
Technical Field
The invention belongs to the technical field of hydrodynamic numerical simulation calculation, relates to a cardiovascular system data platform, and in particular relates to a cardiovascular system data platform built-in processing method based on transient characteristics.
Background
The cardiovascular system is also called the "circulatory system" and is a closed system of pipes consisting of the heart and blood vessels. The heart undergoes rhythmic contraction and relaxation under the control of the nervous system, ensuring that blood circulates in a certain direction and thus can be regarded as a pump. Normally, heart pump blood is delivered through blood vessels to various parts of the whole body. The cardiovascular system data platform is used as a data integration and sharing center, not only needs to acquire various cardiovascular data, but also provides higher requirements for simulation calculation.
Modeling of the cardiovascular system of a human body is an important step for realizing simulation calculation of the cardiovascular system. The existing cardiovascular system data platform is based on three-dimensional CFD modeling to acquire the blood flow state of the whole system, and the method can realize the visualization of the blood flow of the cardiovascular system through calculation and acquire simulation calculation data. However, the three-dimensional CFD system has difficult modeling, poor universality and long calculation time, and needs to meet the calculation requirement through a large server, so that real-time simulation of large platform data is difficult to realize. The one-dimensional modeling method is a conventional calculation method in the field of engineering at present, and has the advantages of simple model, wide applicability and high calculation efficiency, and can realize real-time simulation of a system level. It should be noted that the one-dimensional modeling is a simplified version of the three-dimensional model, the problem of simplifying assumptions of the system must be fully considered, reasonable assumptions must be adopted for necessary information to make the simulation result accurate, and the existing cardiovascular system data platform cannot realize the problems.
Disclosure of Invention
The invention aims to: in order to overcome the defects in the prior art, the method for processing the cardiovascular system data platform built-in based on transient characteristics is provided, and the problems that an existing cardiovascular system data platform simulation model is complex, the data acquisition time is long and the like can be solved.
The technical scheme is as follows: in order to achieve the above object, the present invention provides a method for processing a cardiovascular system data platform based on transient characteristics, comprising the following steps:
S1: constructing a cardiovascular system blood flow transient equation considering dynamic friction and wall viscoelasticity;
s2: establishing a calculation grid, converting a hyperbolic partial differential equation set (a cardiovascular system blood flow transient equation) into a normal differential equation set according to a characteristic line method, and integrating to obtain a characteristic line equation;
S3: solving partial differentiation of the hysteresis strain of the vascular wall to time;
S4: introducing a Kagawa dynamic friction model and partial differentiation into a characteristic line equation;
s5: and (3) considering the heart boundary condition, and obtaining and processing a calculation result according to the characteristic line equation processed in the step (S4).
Further, the cardiovascular system blood flow transient equation constructed in the step S1 is:
Wherein H is the pressure of the blood vessel wall; v is the average flow rate of blood in the vessel; epsilon r is the vascular wall hysteresis strain; a is the speed of sound waves in blood; g is gravity acceleration; x is the distance along the axis of the vessel; t is time; d is the inner diameter of the vessel wall; ρ is the blood density; τ w is the vessel wall shear stress, τ w=τsus and τ u represent steady state and dynamic friction, respectively.
Further, the system of ordinary differential equations in the step S2 is specifically:
Where τ w=τuu=ρfV|V|/8+τu, f is the darcy-vessbach coefficient, when f is a fixed value, it is a steady state friction, when f=f q, it is a quasi-steady state friction, and f q is a quasi-steady state friction coefficient.
Further, in the step S2, the normal differential equation set is integrated along the positive and negative characteristic lines to obtain a characteristic line equation:
Wherein H P、HA、HB is the blood pressure of P, A, B points respectively; q P、QA、QB is the blood flow at P, A, B points, respectively; η is an integral approximation control coefficient, η=0.5 to 1; Δt is the calculation time step.
Further, the step S3 specifically includes:
Solving partial differentiation of the hysteresis strain of the vascular wall to time t, and knowing by a generalized K-V model:
Where ε r (t) is a function of ε r time; alpha is poisson's ratio; e is the thickness of the vessel wall; h (t) is a function of time of the pressure H; j (t) is the creep compliance over time t, t being a differential variable; converting the integral into a sum function, i.e.
Wherein N is the total number of integral segments; epsilon rk (t) is the kth segment value of epsilon r (t); j k and τ k are corresponding viscoelastic mechanical model parameters;
according to a convolution derivative formula, the above derivative is obtained:
I.e.
According to the element conversion method and the fractional integration, the higher-order infinitely small term is omitted, and the method can be obtainedIs a first order approximation solution of (a):
In the method, in the process of the invention,
CO_A3=CO_A1·CO_A2 (14)
CO_A4=CO_A1(1-CO_A2) (15)
Bringing equation (11) into equation (10) yields:
Further, the step S4 specifically includes:
For the Kagawa dynamic friction model:
In the method, in the process of the invention,
Wherein, the coefficient mi=(1;1.16725;2.20064;3.92861;6.78788;11.6761;20.0612;34.4541;59.1642;101.59);ni=(26.3744;72.8033;187.424;536.626;1570.6;4618.13;13601.1;40082.5;118153;348316);ν is the hemodynamic viscosity;
by taking equations (16) and (17) into the feature line equation, there is
In the method, in the process of the invention,
Let b=a/gA, r=fa Δt/2gDA 2, there is
BA=B+ηR|QA| (23)
BB=B+ηR|QB| (25)
At the time of programming, for unified writing, the following parameters can be added:
HP=CP-BPQP (27)
HP=CM+BMQP (28)
In the method, in the process of the invention,
Thus, the first and second substrates are bonded together,
Further, the step S5 specifically includes:
approximating the heart boundary as: a pump boundary and a pressurized container boundary;
For the pump boundary:
HP=HS+QP(a1+a2QP) (32)
wherein H S is a flow breaking head, a1, a2 is a characteristic curve constant;
For a pressurized container boundary, H P is a fixed value;
the boundary conditions are combined to obtain the pressure and flow rate at the boundary.
The invention improves the existing cardiovascular system data platform, combines the correlation between the heart pump blood pressure and the ejection fraction according to the one-dimensional flow assumption of blood in the cardiovascular system, considers the viscoelasticity effect and the friction effect of blood vessels, establishes the one-dimensional cardiovascular system blood flow model for calculating the data such as the heart pump blood pressure, the ejection fraction and the like when the cardiovascular system blood flows, can solve the problems of complex simulation model, long data acquisition time and the like of the existing cardiovascular system data platform, and has important theoretical significance and practical application value for realizing the simulation of the cardiovascular system large data platform.
The beneficial effects are that: compared with the prior art, the invention improves the built-in algorithm of the cardiovascular system data platform, and has the following advantages:
1. The cardiovascular system data platform provided by the invention adopts a one-dimensional model to replace the existing three-dimensional CFD for modeling, breaks through the technical bottleneck that the existing cardiovascular system data platform is difficult to model by adopting the one-dimensional model, and has the advantages of simplicity in modeling, strong universality and short simulation time.
2. The built-in algorithm of the cardiovascular system data platform considers the dynamic friction effect and the viscoelasticity effect of the vascular wall, is more in line with the actual situation, and improves the simulation precision of the platform.
Drawings
FIG. 1 is a flow chart of an implementation of the present invention;
FIG. 2 is a computational grid along the vessel axis;
FIG. 3 is a graph of aortic blood volume experiments versus simulations;
fig. 4 is a graph of ventricular pressure experiments versus simulations.
Detailed Description
The present application is further illustrated in the accompanying drawings and detailed description which are to be understood as being merely illustrative of the application and not limiting of its scope, and various modifications of the application, which are equivalent to those skilled in the art upon reading the application, will fall within the scope of the application as defined in the appended claims.
The invention provides a cardiovascular system data platform built-in processing method based on transient characteristics, as shown in figure 1, comprising the following steps:
S1: constructing a cardiovascular system blood flow transient equation considering dynamic friction and wall viscoelasticity;
Wherein H is the pressure of the blood vessel wall; v is the average flow rate of blood in the vessel; epsilon r is the vascular wall hysteresis strain; a is the speed of sound waves in blood; g is gravity acceleration; x is the distance along the axis of the vessel; t is time; d is the inner diameter of the vessel wall; ρ is the blood density; τ w is the vessel wall shear stress, τ w=τsus and τ u represent steady state and dynamic friction, respectively.
S2: establishing a computational grid along the axial direction of the blood vessel as shown in fig. 2;
converting the hyperbolic partial differential equation set into a normal differential equation set according to a characteristic line method:
Where τ w=τuu=ρfV|V|/8+τu, f is the darcy-vessbach coefficient, when f is a fixed value, it is a steady state friction, when f=f q, it is a quasi-steady state friction, and f q is a quasi-steady state friction coefficient.
And integrating the normal differential equation set along positive and negative characteristic lines to obtain a characteristic line equation:
Wherein H P、HA、HB is the blood pressure of P, A, B points respectively; q P、QA、QB is the blood flow at P, A, B points, respectively; η is an integral approximation control coefficient, η=0.5 to 1; Δt is the calculation time step.
S3: solving partial differentiation of the hysteresis strain of the vascular wall with respect to time:
Solving partial differentiation of the hysteresis strain of the vascular wall to time t, and knowing by a generalized K-V model:
Where ε r (t) is a function of ε r time; alpha is poisson's ratio; e is the thickness of the vessel wall; h (t) is a function of time of the pressure H; j (t) is the creep compliance over time t, t being a differential variable; converting the integral into a sum function, i.e.
Wherein N is the total number of integral segments; epsilon rk (t) is the kth segment value of epsilon r (t); j k and τ k are corresponding viscoelastic mechanical model parameters;
according to a convolution derivative formula, the above derivative is obtained:
I.e.
According to the element conversion method and the fractional integration, the higher-order infinitely small term is omitted, and the method can be obtainedIs a first order approximation solution of (a):
In the method, in the process of the invention,
CO_A3=CO_A1·CO_A2 (14)
CO_A4=CO_A1(1-CO_A2) (15)
Bringing equation (11) into equation (10) yields:
S4: introducing a Kagawa dynamic friction model and partial differentiation into a characteristic line equation;
For the Kagawa dynamic friction model:
In the method, in the process of the invention,
Wherein, the coefficient mi=(1;1.16725;2.20064;3.92861;6.78788;11.6761;20.0612;34.4541;59.1642;101.59);ni=(26.3744;72.8033;187.424;536.626;1570.6;4618.13;13601.1;40082.5;118153;348316);ν is the hemodynamic viscosity;
by taking equations (16) and (17) into the feature line equation, there is
In the method, in the process of the invention,
Let b=a/gA, r=fa Δt/2gDA 2, there is
BA=B+ηR|QA| (23)
BB=B+ηR|QB| (25)
At the time of programming, for unified writing, the following parameters can be added:
HP=CP-BPQP (27)
HP=CM+BMQP (28)
In the method, in the process of the invention,
Thus, the first and second substrates are bonded together,
S5: considering the heart boundary condition, according to the characteristic line equation processed in the step S4, obtaining and processing a calculation result:
approximating the heart boundary as: a pump boundary and a pressurized container boundary;
For the pump boundary:
HP=HS+QP(a1+a2QP) (32)
wherein H S is a flow breaking head, a1, a2 is a characteristic curve constant;
For a pressurized container boundary, H P is a fixed value;
the boundary conditions are combined to obtain the pressure and flow rate at the boundary.
In order to verify the use effect of the built-in model of the cardiovascular system data platform provided by the invention, the verification is performed by using the existing experimental data, and the main experimental parameters are as follows: the heart rate was 70 times/min, and the single shot blood volume was 802ml, aortic systolic pressure 1205mmHg, diastolic pressure 805mmHg, pulmonary systolic pressure 252mmHg, and diastolic pressure 81.5mmHg. The experimental and simulated comparison results are shown in fig. 3 and 4.
Fig. 3 shows experimental and simulated contrast plots of aortic blood volume over time. From fig. 3, the simulation result is well matched with the experimental result, and the model can well simulate the change of the aortic blood volume and the maximum and minimum flow values. Fig. 4 shows an experimental versus simulated comparison of ventricular pressure over time. As can be seen from fig. 4, the simulation result of the model is substantially identical to the actual ventricular pressure variation result, and the model can accurately reflect the pressure peak value in each period of the ventricle. For both fig. 3 and fig. 4, there is some degree of phase shift in the pressure and flow simulation results compared to the experiment, which may be due to inaccuracy in the simulated wave velocity but more accurate overall simulation results.

Claims (6)

1. The cardiovascular system data platform built-in processing method based on transient characteristics is characterized by comprising the following steps of:
S1: constructing a cardiovascular system blood flow transient equation considering dynamic friction and wall viscoelasticity;
S2: establishing a calculation grid, converting a cardiovascular system blood flow transient equation into a normal differential equation set according to a characteristic line method, and integrating to obtain a characteristic line equation;
S3: solving partial differentiation of the hysteresis strain of the vascular wall to time;
S4: introducing a Kagawa dynamic friction model and partial differentiation into a characteristic line equation;
s5: considering boundary conditions, and obtaining and processing a calculation result according to the characteristic line equation processed in the step S4;
the cardiovascular system blood flow transient equation constructed in the step S1 is as follows:
Wherein H is the pressure of the blood vessel wall; v is the average flow rate of blood in the vessel; epsilon r is the vascular wall hysteresis strain; a is the speed of sound waves in blood; g is gravity acceleration; x is the distance along the axis of the vessel; t is time; d is the inner diameter of the vessel wall; ρ is the blood density; τ w is the vessel wall shear stress, τ w=τsus and τ u represent steady state and dynamic friction, respectively.
2. The method for processing the cardiovascular system data platform according to claim 1, wherein the system of ordinary differential equations in step S2 is specifically:
Where τ w=τsu=ρfV|V|/8+τu, f is the darcy-vessbach coefficient, when f is a fixed value, it is a steady state friction, when f=f q, it is a quasi-steady state friction, and f q is a quasi-steady state friction coefficient.
3. The method for processing the cardiovascular system data platform according to claim 2, wherein in the step S2, the normal differential equation set is integrated along positive and negative characteristic lines to obtain characteristic line equations:
Wherein H P、HA、HB is the blood pressure of P, A, B points respectively; q P、QA、QB is the blood flow at P, A, B points, respectively; η is an integral approximation control coefficient, η=0.5 to 1; Δt is the calculation time step.
4. A method for processing a cardiovascular system data platform according to claim 3, wherein the step S3 specifically comprises:
Solving partial differentiation of the hysteresis strain of the vascular wall to time t, and knowing by a generalized K-V model:
Where ε r (t) is a function of ε r time; alpha is poisson's ratio; e is the thickness of the vessel wall; h (t) is a function of time of the pressure H; j (t) is the creep compliance over time t, t being a differential variable; converting the integral into a sum function, i.e.
Wherein N is the total number of integral segments; epsilon rk (t) is the kth segment value of epsilon r (t); j k and τ k are corresponding viscoelastic mechanical model parameters;
according to a convolution derivative formula, the above derivative is obtained:
I.e.
According to the element conversion method and the fractional integration, the higher-order infinitely small term is omitted, and the method can be obtainedIs a first order approximation solution of (a):
In the method, in the process of the invention,
CO_A3=CO_A1·CO_A2 (14)
CO_A4=CO_A1(1-CO_A2) (15)
Bringing equation (11) into equation (10) yields:
5. The method for processing the cardiovascular system data platform according to claim 4, wherein the step S4 specifically comprises:
For the Kagawa dynamic friction model:
In the method, in the process of the invention,
Wherein v is hemodynamic viscosity;
by taking equations (16) and (17) into the feature line equation, there is
In the method, in the process of the invention,
Let b=a/gA, r=fa Δt/2gDA 2, there is
BA=B+ηR|QA| (23)
BB=B+ηR|QB| (25)
At the time of programming, for unified writing, the following parameters are added:
HP=CP-BPQP (27)
HP=CM+BMQP (28)
In the method, in the process of the invention,
Thus, the first and second substrates are bonded together,
6. The method for processing the cardiovascular system data platform according to claim 5, wherein the step S5 specifically comprises:
Approximating the boundary as: a pump boundary and a pressurized container boundary;
For the pump boundary:
HP=HS+QP(a1+a2QP) (32)
wherein H S is a flow breaking head, a1, a2 is a characteristic curve constant;
For a pressurized container boundary, H P is a fixed value;
the boundary conditions are combined to obtain the pressure and flow rate at the boundary.
CN202110987288.6A 2021-08-26 2021-08-26 Cardiovascular system data platform built-in processing method based on transient characteristics Active CN113674865B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110987288.6A CN113674865B (en) 2021-08-26 2021-08-26 Cardiovascular system data platform built-in processing method based on transient characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110987288.6A CN113674865B (en) 2021-08-26 2021-08-26 Cardiovascular system data platform built-in processing method based on transient characteristics

Publications (2)

Publication Number Publication Date
CN113674865A CN113674865A (en) 2021-11-19
CN113674865B true CN113674865B (en) 2024-04-26

Family

ID=78546600

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110987288.6A Active CN113674865B (en) 2021-08-26 2021-08-26 Cardiovascular system data platform built-in processing method based on transient characteristics

Country Status (1)

Country Link
CN (1) CN113674865B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6117087A (en) * 1998-04-01 2000-09-12 Massachusetts Institute Of Technology Method and apparatus for noninvasive assessment of a subject's cardiovascular system
CN105095534A (en) * 2014-04-23 2015-11-25 北京冠生云医疗技术有限公司 Method and system for simulation of bloodstream in blood vessels
CN108765543A (en) * 2018-03-28 2018-11-06 向建平 A kind of Cardiovascular imaging reconstructing method and system based on implicit deformable model
CN110296327A (en) * 2019-06-19 2019-10-01 常州大学 A kind of pipeline leakage detection method based on transient flow frequency response analysis
CN111916183A (en) * 2019-05-10 2020-11-10 北京航空航天大学 Dynamic cardiovascular system modeling method, device, equipment and storage medium

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10980427B2 (en) * 2017-06-22 2021-04-20 Dextera Medical Method and apparatus for full-system, cardiovascular simulation and prediction

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6117087A (en) * 1998-04-01 2000-09-12 Massachusetts Institute Of Technology Method and apparatus for noninvasive assessment of a subject's cardiovascular system
CN105095534A (en) * 2014-04-23 2015-11-25 北京冠生云医疗技术有限公司 Method and system for simulation of bloodstream in blood vessels
CN108765543A (en) * 2018-03-28 2018-11-06 向建平 A kind of Cardiovascular imaging reconstructing method and system based on implicit deformable model
CN111916183A (en) * 2019-05-10 2020-11-10 北京航空航天大学 Dynamic cardiovascular system modeling method, device, equipment and storage medium
CN110296327A (en) * 2019-06-19 2019-10-01 常州大学 A kind of pipeline leakage detection method based on transient flow frequency response analysis

Non-Patent Citations (12)

* Cited by examiner, † Cited by third party
Title
180°弯曲方管牛顿流体及一种非牛顿流体湍性流动的数值模拟;冯桐, 钱伟长, 孙厚钧;应用数学和力学;19910115(第01期);6-15 *
Cardiovascular System: Modelling and Optimization;D’Arienzo Maria Pia;SEMANTIC SCHOLAR;18-117 *
Parameters affecting water-hammer wave attenuation, shape and timing;BERGANT A等;Journal of Hydraulic Research;46(3);373-381 *
The dynamic effect of pipe-wall viscoelasticity in hydraulic transients;COVAS D等;Journal of Hydraulic Research;43(1);56-70 *
Understanding aspects of cardiovascular physiology and disease via a multi-physics modelling methodology;Alberto Coccarelli;SEMANTIC SCHOLAR;23-159 *
国外有压瞬变流研究进展;***;河海科技进展;第11卷(第2期);9-13 *
工程管道与人体血管之比较分析;陈贵清;丁幼松;;唐山学院学报;20080720(第04期);45-48 *
微血管自律运动的流体力学及血液粘弹性质的影响;郭仲三;应用数学和力学;19920215(第02期);59-65 *
心血管内流场的高质量模拟实验研究(Ⅰ)──动态流体双折射测试原理与分析方法;孙永达等;实验力学;11(1);8-11 *
心血管***多元数字计算机模型及其应用前景;白净;中国生物医学工程学报(第03期);51-55 *
结合呼吸运动的改进型多分支心血管***建模与仿真;王末;王群;刘志文;;北京生物医学工程(第01期);47-52、113 *
黏弹性输水管道中含气瞬变流压力衰减分析;Alberto Coccarelli;水利学报;第49卷(第3期);303-312 *

Also Published As

Publication number Publication date
CN113674865A (en) 2021-11-19

Similar Documents

Publication Publication Date Title
CA2612317A1 (en) Mathematical circulatory system model
Yu et al. Computational fluid dynamics‐based design optimization method for archimedes screw blood pumps
Taber et al. Computational model for the transition from peristaltic to pulsatile flow in the embryonic heart tube
Darowski et al. A new hybrid (hydro-numerical) model of the circulatory system
Gregory et al. Simulation and enhancement of a cardiovascular device test rig
CN113674865B (en) Cardiovascular system data platform built-in processing method based on transient characteristics
Franzetti et al. Design of an in vitro mock circulatory loop to reproduce patient-specific vascular conditions: Toward precision medicine
Chung et al. Mechanical circulatory support devices: methods to optimize hemodynamics during use
CN109931056B (en) A kind of change control reserve dynamic analysis of gas well method
Ferrari et al. Hybrid model analysis of intra‐aortic balloon pump performance as a function of ventricular and circulatory parameters
Mamun et al. Comparative study of non-Newtonian physiological blood flow through the elastic stenotic artery with rigid wall stenotic artery
Bahraseman et al. Combining numerical and clinical methods to assess aortic valve hemodynamics during exercise
Huberts et al. A lumped model for blood flow and pressure in the systemic arteriesbased on an approximate velocity profile function
Schampaert et al. Modeling the interaction between the intra-aortic balloon pump and the cardiovascular system: the effect of timing
CN110457765B (en) Method for calculating odd vein shunt ratio in half-ventricle model based on geometric multi-scale model
CN102128163B (en) For determining the emic method of pump
Rideout et al. Computer simulation study of the pulmonary circulation
Rubtsova et al. Mathematical modeling of the Fontan blood circulation supported with pediatric ventricular assist device
CN116850445A (en) Method, device, equipment and storage medium for controlling heart pump
Chen et al. Dynamic Model and Characteristic Analysis of Viscosity‐Ultraelasticity for Bionic Vascular Network
Wang et al. Shear stress and hemolysis analysis of blood pump under constant and pulsation speed based on a multiscale coupling model
Yang et al. Numerical study of blood flow at the end-to-side anastomosis of a left ventricular assist device for adult patients
CN113128139A (en) Method and system for rapidly calculating fractional flow reserve based on simplified coronary artery zero-dimensional model and stenosis resistance prediction model
CN109994198A (en) A kind of method for reconstructing of dynamic heart threedimensional model
Wu et al. A new integrated numerical modeling approach toward piston diaphragm pump simulation: Three-dimensional model simplification and characteristic analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant