• Tissue-growth based Synthetic Tree Generation and Perfusion Simulation (Biomechanics and Modeling in Mechanobiology (2023) 22:1095–1112)

 Biological tissues receive oxygen and nutrients from blood vessels by developing an indispensable supply and demand relationship with the blood vessels. We implemented a synthetic tree generation algorithm by considering the interactions between the tissues and blood vessels. We first segment major arteries using medical image data and synthetic trees are generated originating from these segmented arteries. They grow into extensive  networks of small vessels to fill the supplied tissues and satisfy the metabolic demand of them. Further, the algorithm is optimized to be executed in parallel without affecting the generated tree volumes. The generated vascular trees are used to simulate blood perfusion in the tissues by performing multiscale blood flow simulations. One-dimensional blood flow equations were used to solve for blood flow and pressure in the generated vascular trees and Darcy flow equations were solved for blood perfusion in the tissues using a porous model assumption. Both equations are coupled at terminal segments explicitly. The proposed methods were applied to idealized models with different tree resolutions and metabolic demands for validation. The methods demonstrated that realistic synthetic trees were generated with significantly less computational expense compared to that of a constrained constructive optimization method. The methods were then applied to cerebrovascular arteries supplying a human brain and coronary arteries supplying the left and right ventricles to demonstrate the capabilities of the proposed methods. The proposed methods can be utilized to quantify tissue perfusion and predict areas prone to ischemia in patient-specific geometries.

• Cerebral Perfusion Simulation by Realistic Synthetic Tree Generation: Application to Healthy and Stroke Patients (submitted to Radiology)

Patient-specific cerebral vascular models with realistic morphological characteristics constructed from non-invasively acquired MRI data using growth based optimization were shown to predict perfusion territories of healthy and stroke patients realistically.

• Physiological model of the cerebrovascular system based on a supply and demand relationship between arteries and tissues (in progress)

• Unseen branch model based on physiological principles for perfusion simulation of coronary arteries (in progress)

• Convergence of phase-averaged, transitional flow in an abdominal aortic aneurysmal model (to appear in ASME Journal of Biomechanical Engineering)

Abdominal aortic aneurysm can exhibit transitional flow characteristics in laminar flow regimes. To report transitional flow characteristics, we examined the convergence of phase-averaged solutions by executing blood flow simulations of a patient-specific abdominal aortic aneurysmal model for 257 cardiac cycles with periodic, pulsatile boundary conditions. The phase-averaged solutions were computed by averaging the solutions over various numbers of cardiac cycles and compared against the ones averaged over 124 cycles. The phase-averaged solutions reported small differences when they were averaged over a large number of cardiac cycles. The instantaneous solutions, however, failed to exhibit fluctuations reported in the phase-averaged solutions. To study transitional blood flows in the aneurysmal region we need to report phase-averaged solutions as they exhibit non-periodic, disturbed flow characteristics. Additionally, when reporting phase-averaged solutions it is preferred to compute an average over a large number of cardiac cycles to be able to represent flow structures of the converged phase-averaged solutions.

• Particulate Matter Simulation (International Journal of Molecular Sciences 23 (15), 8170)

• Neutron transport using stabilized finite element methods (in progress)

• Material transport inside the cerebrospinal fluid (in progress)

• Compression simulation of lower extremities using fluid-solid interaction methods (in progress)

• Multiscale computational methods for heat transfer simulation between body tissue and blood flow (in progress)

• Vascular stent deployment and fatigue analysis (in progress) 

• Growth and remodeling of abdominal aortic aneurysmal models (in progress)

• Uncertainty quantification in coronary blood flow simulations: Impact of geometry, boundary conditions and blood viscosity (Journal of biomechanics 49 (12), 2540-2547)

Computational fluid dynamic methods are currently being used clinically to simulate blood flow and pressure and predict the functional significance of atherosclerotic lesions in patient-specific models of the coronary arteries extracted from noninvasive coronary computed tomography angiography data. One such technology, or noninvasive fractional flow reserve derived from CT data, has demonstrated high diagnostic accuracy as compared to invasively measured fractional flow reserve (FFR) obtained with a pressure wire inserted in the coronary arteries during diagnostic cardiac catheterization. However, uncertainties in modeling as well as measurement results in differences between these predicted and measured hemodynamic indices. Uncertainty in modeling can manifest in two forms – anatomic uncertainty resulting in error of the reconstructed 3D model and physiologic uncertainty resulting in errors in boundary conditions or blood viscosity. We present a data-driven framework for modeling these uncertainties and study their impact on blood flow simulations. The incompressible Navier–Stokes equations are used to model blood flow and an adaptive stochastic collocation method is used to model uncertainty propagation in the Navier–Stokes equations. We perform uncertainty quantification in two geometries, an idealized stenosis model and a patient specific model. We show that uncertainty in minimum lumen diameter (MLD) has the largest impact on hemodynamic simulations, followed by boundary resistance, viscosity and lesion length. We show that near the diagnostic cutoff, the uncertainty due to the latter three variables are lower than measurement uncertainty, while the uncertainty due to MLD is only slightly higher than measurement uncertainty. We also show that uncertainties are not additive but only slightly higher than the highest single parameter uncertainty. The method presented here can be used to output interval estimates of hemodynamic indices and visualize patient-specific maps of sensitivities.

• Reduced order model for pulsatile flow simulations (in progress)

• Data-driven modeling for blood flow simulations (in progress)

• Efficient mesh generation and adaptation for various applications (in progress)

• 개인 맞춤형 의료진단기술 기반의 심뇌혈관 M3DT 기술개발 및 검증, 한국보건산업진흥원  (연구 책임)

• 디지털 트윈 기술 기반 뇌 모니터링을 위한 멀티 스케일 뇌 역학 모델, 삼성미래기술육성사업 (연구 책임)

• 스마트 혈류역학 지표 기반 실시간 정밀의료 플랫폼 연구센터, 한국연구재단 (공동 연구)

• 사용후연료 염소화 및 용융염고속로를 이용한 보다 지속가능한 원자력 시스템 개념 연구, 한국과학기술원 (공동 연구)

• 스마트야드 전문인력양성, 중소조선연구원 (공동 연구)

• 척추온열도자의 혈액순환 효과를 예측하는 해석 모델 개발, 세라젬 (연구책임)

• 심장 역학과 결합한 생리학적인 멀티스케일 관상동맥 혈액 유동 해석 모델 개발, 한국연구재단 (연구 책임)

• 심장혈관계 질환의 진단을 위한 계산 방법 개발과 검증, 한국과학기술원 (연구 책임)

• 사지 압박 장치와 온열치료의 메커니즘을 모사하는 해석 모델 개발, 세라젬 (연구 책임)

Lab introduction slides

Lab introduction PDF

KAIST ME CVBML 연구실 소개 영상 

[메릭 웨비나] 기계 학습과 해석 방법을 이용한 심장혈관계 질환의 진단과 의료기기 해석 - 김현진 교수 (KAIST 기계공학과) 

[메릭 웨비나] 허혈 진단을 위한 개인화된 혈관 생성과 관류 해석 - 김현진 교수 (KAIST 기계공학과) 

[유체공학 부문] 심장혈관계 질환의 진단과 치료결과 예측을 위한 유체역학 - 김현진 교수 (KAIST 기계공학과)