Browsing by Author "Yu, Huidan"
Now showing 1 - 10 of 17
Results Per Page
Sort Options
Item A new noninvasive and patient-specific hemodynamic index for the severity of renal stenosis and outcome of interventional treatment(Wiley, 2022-07) Yu, Huidan; Khan, Monsurul; Wu, Hao; Du, Xiaoping; Chen, Rou; Rollins, Dave M.; Fang, Xin; Long, Jianyun; Xu, Chenke; Sawchuk, Alan P.; Surgery, School of MedicineRenal arterial stenosis (RAS) often causes renovascular hypertension, which may result in kidney failure and life-threatening consequences. Direct assessment of the hemodynamic severity of RAS has yet to be addressed. In this work, we present a computational concept to derive a new, noninvasive, and patient-specific index to assess the hemodynamic severity of RAS and predict the potential benefit to the patient from a stenting therapy. The hemodynamic index is derived from a functional relation between the translesional pressure indicator (TPI) and lumen volume reduction (S) through a parametric deterioration of the RAS. Our in-house computational platform, InVascular, for image-based computational hemodynamics is used to compute the TPI at given S. InVascular integrates unified computational modeling for both image processing and computational hemodynamics with graphic processing unit parallel computing technology. The TPI-S curve reveals a pair of thresholds of S indicating mild or severe RAS. The TPI at S = 0 represents the pressure improvement following a successful stenting therapy. Six patient cases with a total of 6 aortic and 12 renal arteries are studied. The computed blood pressure waveforms have good agreements with the in vivo measured ones and the systolic pressure is statistical equivalence to the in-vivo measurements with p < .001. Uncertainty quantification provides the reliability of the computed pressure through the corresponding 95% confidence interval. The severity assessments of RAS in four cases are consistent with the medical practice. The preliminary results inspire a more sophisticated investigation for real medical insights of the new index. This computational concept can be applied to other arterial stenoses such as iliac stenosis. Such a noninvasive and patient-specific hemodynamic index has the potential to aid in the clinical decision-making of interventional treatment with reduced medical cost and patient risks.Item Computational Fluid Dynamics for Modeling and Simulation of Intraocular Drug Delivery and Wall Shear Stress in Pulsatile Flow(2020-08) Abootorabi, Seyedalireza; Yu, Huidan; Nematollahi, Khosrow; Yokota, HirokiThe thesis includes two application studies of computational fluid dynamics. The first is new and efficient drug delivery to the posterior part of the eye, a growing health necessity worldwide. Current treatment of eye diseases, such as age-related macular degeneration (AMD), relies on repeated intravitreal injections of drug-containing solutions. Such a drug delivery has significant cant drawbacks, including short drug life, vital medical service, and high medical costs. In this study, we explore a new approach of controlled drug delivery by introducing unique porous implants. Computational modeling contains physiological and anatomical traits. We simulate the IgG1 Fab drug delivery to the posterior eye to evaluate the effectiveness of the porous implants to control the drug delivery. The computational model was validated by established computation results from independent studies and experimental data. Overall, the results indicate that therapeutic drug levels in the posterior eye are sustained for eight weeks, similar to those performed with intravitreal injection of the same drug. We evaluate the effects of the porous implant on the time evaluation of the drug concentrations in the sclera, choroid, and retina layers of the eye. Subsequent simulations were carried out with varying porosity values of a porous episcleral implant. Our computational results reveal that the time evolution of drug concentration is distinctively correlated to drug source location and pore size. The response of this porous implant for controlled drug delivery applications was examined. A correlation between porosity and fluid properties for the porous implants was revealed in this study. The second application lays in the computational modeling of the oscillatingItem The development of polystyrene based microfluidic gas generation system(2015-05) Yuanzhi, Cao; Zhu, Likun; Yu, Huidan; El-Mounayri, Hazim A.; Anwar, SohelThe purpose of this thesis is to use experimental methods to seek deeper understanding and better performance in the self-circulating self-regulating microfluidic gas generator initially developed in Dr. Zhu’s group, by changing the major features and dimensions in the reaction channel of the device. In order to effectively conduct experiments described above, a microfabrication method that is capable of making new microfluidic devices with low cost, short time period, as well as relatively high accuracy was needed first. Developing such a fabrication method is the major part of this thesis. We initially used patterned polymer films and glass slide, and bonded them together by sequentially aligning and stacking them into a microfluidic device with patterned double-sided tapes. Later we developed a more advanced microfabrication method that used only patterned polystyrene (PS) films. The patterned PS films were obtained from a digital cutter and they were bonded into a microfluidic device by thermopress bonding method that required no heterogeneous bonding agents. This new method did not need manual assembly which greatly improved its precision (~ 100 µm), and it used only PS as device material that has favorable surface wetting property for microfluidics applications. In order to find the optimized microfluidic channel design to improve gas generating performance, we've designed and fabricated microfluidic devices with different channel dimensions using the PS fabrication method. Based on the gas generation testing results of those devices, we were able to come up with the optimal dimensions for the reaction channel that had the best gas generation performance. To obtain a more fundamental understanding about the working mechanism of our device and its bubble dynamics, we have conducted ultrafast X-ray imaging test at Advanced Photon Source (APS), Argonne National Laboratory. High speed (100 KHz) phase contrast images were captured that allowed us to observe directly inside the reaction channel on the cross section view during the self-circulating catalytic reaction. The images provided us with lots of insightful information that in turn helped the dimensional improvement for the microchannel design. The 100 KHz high speed images also gave us useful information about the dynamics of bubble development on a catalyst bed, such as growth and merging of the bubbles.Item Effects of donor-specific microvascular anatomy on hemodynamic perfusion in human choriocapillaris(Springer Nature, 2023-12-19) An, Senyou; Yu, Huidan; Islam, MD Mahfuzul; Zhang, Xiaoyu; Zhan, Yuting; Olivieri, Joseph J.; Ambati, Jayakrishna; Yao, Jun; Gelfand, Bradley D.; Mechanical and Energy Engineering, Purdue School of Engineering and TechnologyEvidence from histopathology and clinical imaging suggest that choroidal anatomy and hemodynamic perfusion are among the earliest changes in retinal diseases such as age-related macular degeneration (AMD). However, how inner choroidal anatomy affects hemodynamic perfusion is not well understood. Therefore, we sought to understand the influences of choroidal microvascular architecture on the spatial distribution of hemodynamic parameters in choriocapillaris from human donor eyes using image-based computational hemodynamic (ICH) simulations. We subjected image-based inner choroid reconstructions from eight human donor eyes to ICH simulation using a kinetic-based volumetric lattice Boltzmann method to compute hemodynamic distributions of velocity, pressure, and endothelial shear stress. Here, we demonstrate that anatomic parameters, including arteriolar and venular arrangements and intercapillary pillar density and distribution exert profound influences on inner choroidal hemodynamic characteristics. Reductions in capillary, arteriolar, and venular density not only reduce the overall blood velocity within choriocapillaris, but also substantially increase its spatial heterogeneity. These first-ever findings improve understanding of how choroidal anatomy affects hemodynamics and may contribute to pathogenesis of retinal diseases such as AMD.Item Image Based Computational Hemodynamics for Non-Invasive and Patient-Specific Assessment of Arterial Stenosis(2019-08) Khan, Md Monsurul Islam; Yu, Huidan; Wagner, Diane; Zhu, LikunWhile computed tomographic angiography (CTA) has emerged as a powerful noninvasive option that allows for direct visualization of arterial stenosis(AS), it cant assess the hemodynamic abnormality caused by an AS. Alternatively, trans-stenotic pressure gradient (TSPG) and fractional flow reserve (FFR) are well-validated hemodynamic indices to assess the ischemic severity of an AS. However, they have significant restriction in practice due to invasiveness and high cost. To fill the gap, a new computational modality, called InVascular has been developed for non-invasive quantification TSPG and/or FFR based on patient's CTA, aiming to quantify the hemodynamic abnormality of the stenosis and help to assess the therapeutic/surgical benefits of treatment for the patient. Such a new capability gives rise to a potential of computation aided diagnostics and therapeutics in a patient-specific environment for ASs, which is expected to contribute to precision planning for cardiovascular disease treatment. InVascular integrates a computational modeling of diseases arteries based on CTA and Doppler ultrasonography data, with cutting-edge Graphic Processing Unit (GPU) parallel-computing technology. Revolutionary fast computing speed enables noninvasive quantification of TSPG and/or FFR for an AS within a clinic permissible time frame. In this work, we focus on the implementation of inlet and outlet boundary condition (BC) based on physiological image date and and 3-element Windkessel model as well as lumped parameter network in volumetric lattice Boltzmann method. The application study in real human coronary and renal arterial system demonstrates the reliability of the in vivo pressure quantification through the comparisons of pressure waves between noninvasive computational and invasive measurement. In addition, parametrization of worsening renal arterial stenosis (RAS) and coronary arterial stenosis (CAS) characterized by volumetric lumen reduction (S) enables establishing the correlation between TSPG/FFR and S, from which the ischemic severity of the AS (mild, moderate, or severe) can be identified. In this study, we quantify TSPG and/or FFR for five patient cases with visualized stenosis in coronary and renal arteries and compare the non-invasive computational results with invasive measurement through catheterization. The ischemic severity of each AS is predicted. The results of this study demonstrate the reliability and clinical applicability of InVascular.Item Inlet and Outlet Boundary Conditions and Uncertainty Quantification in Volumetric Lattice Boltzmann Method for Image-Based Computational Hemodynamics(MDPI, 2022-01-10) Yu, Huidan; Khan, Monsurul; Wu, Hao; Zhang, Chunze; Du, Xiaoping; Chen, Rou; Fang, Xin; Long, Jianyun; Sawchuk, Alan P.; Surgery, School of MedicineInlet and outlet boundary conditions (BCs) play an important role in newly emerged image-based computational hemodynamics for blood flows in human arteries anatomically extracted from medical images. We developed physiological inlet and outlet BCs based on patients’ medical data and integrated them into the volumetric lattice Boltzmann method. The inlet BC is a pulsatile paraboloidal velocity profile, which fits the real arterial shape, constructed from the Doppler velocity waveform. The BC of each outlet is a pulsatile pressure calculated from the three-element Windkessel model, in which three physiological parameters are tuned by the corresponding Doppler velocity waveform. Both velocity and pressure BCs are introduced into the lattice Boltzmann equations through Guo’s non-equilibrium extrapolation scheme. Meanwhile, we performed uncertainty quantification for the impact of uncertainties on the computation results. An application study was conducted for six human aortorenal arterial systems. The computed pressure waveforms have good agreement with the medical measurement data. A systematic uncertainty quantification analysis demonstrates the reliability of the computed pressure with associated uncertainties in the Windkessel model. With the developed physiological BCs, the image-based computation hemodynamics is expected to provide a computation potential for the noninvasive evaluation of hemodynamic abnormalities in diseased human vessels.Item Inlet and Outlet Boundary Conditions and Uncertainty Quantification in Volumetric Lattice Boltzmann Method for Image-Based Computational Hemodynamics(MDPI, 2022) Yu, Huidan; Khan, Monsurul; Wu, Hao; Zhang, Chunze; Du, Xiaoping; Chen, Rou; Fang, Xin; Long, Jianyun; Sawchuk, Alan P.; Mechanical and Energy Engineering, School of Engineering and TechnologyInlet and outlet boundary conditions (BCs) play an important role in newly emerged image-based computational hemodynamics for blood flows in human arteries anatomically extracted from medical images. We developed physiological inlet and outlet BCs based on patients’ medical data and integrated them into the volumetric lattice Boltzmann method. The inlet BC is a pulsatile paraboloidal velocity profile, which fits the real arterial shape, constructed from the Doppler velocity waveform. The BC of each outlet is a pulsatile pressure calculated from the three-element Windkessel model, in which three physiological parameters are tuned by the corresponding Doppler velocity waveform. Both velocity and pressure BCs are introduced into the lattice Boltzmann equations through Guo’s non-equilibrium extrapolation scheme. Meanwhile, we performed uncertainty quantification for the impact of uncertainties on the computation results. An application study was conducted for six human aortorenal arterial systems. The computed pressure waveforms have good agreement with the medical measurement data. A systematic uncertainty quantification analysis demonstrates the reliability of the computed pressure with associated uncertainties in the Windkessel model. With the developed physiological BCs, the image-based computation hemodynamics is expected to provide a computation potential for the noninvasive evaluation of hemodynamic abnormalities in diseased human vessels.Item Laminar and Turbulent Behavior Captured by A 3-D Kinetic-Based Discrete Dynamic System(NSF-PAR, 2022-07) Zhang, Xiaoyu; McDonough, J. M.; Yu, Huidan; Surgery, School of MedicineWe have derived a 3-D kinetic-based discrete dynamic system (DDS) from the lattice Boltzmann equation (LBE) for incompressible flows through a Galerkin procedure. Expressed by a poor-man lattice Boltzmann equation (PMLBE), it involves five bifurcation parameters including relaxation time from the LBE, splitting factor of large and sub-grid motion scales, and wavevector components from the Fourier space. Numerical experiments have shown that the DDS can capture laminar behaviors of periodic, subharmonic, n-period, and quasi-periodic and turbulent behaviors of noisy periodic with harmonic, noisy subharmonic, noisy quasi-periodic, and broadband power spectra. In this work, we investigated the effects of bifurcation parameters on the capturing of the laminar and turbulent flows in terms of the convergence of time series and the pattern of power spectra. We have found that the 2nd order and 3rd order PMLBEs are both able to capture laminar and turbulent flow behaviors but the 2nd order DDS performs better with lower computation cost and more flow behaviors captured. With the specified ranges of the bifurcation parameters, we have identified two optimal bifurcation parameter sets for laminar and turbulent behaviors. Beyond this work, we are exploring the regime maps for a deeper understanding of the contributions of the bifurcation parameters to the capturing of laminar and turbulent behaviors. Surrogate models (to replace the PMLBE) are being developed using deep learning techniques to overcome the overwhelming computation cost for the regime maps. Meanwhile, the DDS is being employed in the large eddy simulation of turbulent pulsatile flows to provide dynamic sub-grid scale information.Item A new criterion of coalescence-induced microbubble detachment in three-dimensional microfluidic channel(AIP, 2021-04) Chen, Rou (陈柔 ); Zhou, Shuiyi (周书屹 ); Zhu, Likun; Zhu, Luoding; Yan, Weiwei (严微微 ); Yu, Huidan; Mechanical Engineering, School of Engineering and TechnologyThis work is motivated by an experiment of microbubble transport in a polymer microfluidic gas generation device where coalescence-induced detachment exhibits. We numerically study three-dimensional microbubble coalescence using the graphics processing unit accelerating free energy lattice Boltzmann method with cubic polynomial boundary conditions. The focus is on the coalescence-induced microbubble detachment (CIMD) in microfluidics. From the experimental observation, we identified that size inequality between two-parent bubbles and the size of the father (large) bubble are key factors to determine if a CIMD will occur. First, the analytical relationship between equilibrium contact angle and dimensionless wetting potential and experimental results of coalescence with and without CIMD are employed for the verification and validation, respectively. From eighteen experimental and computational cases, we derive a new criterion for CIMD: CIMD occurs when the two-parent bubbles are (nearly) equal with a relatively large radius. The underlying mechanism behind this criterion is explored by the time evolution of the velocity vector field, vorticity field, and kinetic energy in the entire coalescence. It is found that the symmetric capillary force drives the formation of vertical flow stream to the horizontal alignment of parent bubbles and the blockage of the downward stream due to the solid interface promotes the intensity of the upward stream. Meanwhile, large-sized parent bubbles transfer a large amount of kinetic energy from the initial free surface energy, which is essential to lead a CIMD in the post-coalescence stage. Such a new criterion is expected to impact the design and optimization of microfluidics in various applications.Item A new noninvasive and patient-specific hemodynamic index for the severity of renal stenosis and outcome of interventional treatment(Wiley, 2022) Yu, Huidan; Khan, Monsurul; Wu, Hao; Du, Xiaoping; Chen, Rou; Rollins, Dave M.; Fang, Xin; Long, Jianyun; Xu, Chenke; Sawchuk, Alan P.; Mechanical and Energy Engineering, School of Engineering and TechnologyRenal arterial stenosis (RAS) often causes renovascular hypertension, which may result in kidney failure and life-threatening consequences. Direct assessment of the hemodynamic severity of RAS has yet to be addressed. In this work, we present a computational concept to derive a new, noninvasive, and patient-specific index to assess the hemodynamic severity of RAS and predict the potential benefit to the patient from a stenting therapy. The hemodynamic index is derived from a functional relation between the translesional pressure indicator (TPI) and lumen volume reduction (S) through a parametric deterioration of the RAS. Our in-house computational platform, InVascular, for image-based computational hemodynamics is used to compute the TPI at given S. InVascular integrates unified computational modeling for both image processing and computational hemodynamics with GPU parallel computing technology. The TPI-S curve reveals a pair of thresholds of S indicating mild or severe RAS. The TPI at S=0 represents the pressure improvement following a successful stenting therapy. Six patient cases with a total of 6 aortic and 12 renal arteries are studied. The computed blood pressure waveforms have good agreements with the in-vivo measured ones and the systolic pressure is statistical equivalence to the in-vivo measurements with p<0.001. Uncertainty quantification provides the reliability of the computed pressure through the corresponding 95% confidence interval. The severity assessments of RAS in four cases are consistent with the medical practice. The preliminary results inspire a more sophisticated investigation for real medical insights of the new index. This computational concept can be applied to other arterial stenoses such as iliac stenosis. Such a noninvasive and patient-specific hemodynamic index has the potential to aid in the clinical decision-making of interventional treatment with reduced medical cost and patient risks.