Current Search: Fluid-Structure interaction (x)
View All Items
- Title
- Investigation of Multiscale Fluid Structure Interaction Modeling of Flow in Arterial Systems.
- Creator
-
Sotelo, Sebastian, Kassab, Alain, Ilie, Marcel, Divo, Eduardo, University of Central Florida
- Abstract / Description
-
The study of hemodynamic patterns in large blood vessels, such as the ascending aortic artery, brachiocephalic trunk, right carotid artery and right subclavian artery presents the challenging complexity of vessel wall compliance induced by the high levels of shear stress gradients and blood flow pulsatility. Accurate prediction of hemodynamics in such conditions requires a complete Fluid Structure Interaction (FSI) analysis that couples the fluid flow behavior throughout the cardiac cycle...
Show moreThe study of hemodynamic patterns in large blood vessels, such as the ascending aortic artery, brachiocephalic trunk, right carotid artery and right subclavian artery presents the challenging complexity of vessel wall compliance induced by the high levels of shear stress gradients and blood flow pulsatility. Accurate prediction of hemodynamics in such conditions requires a complete Fluid Structure Interaction (FSI) analysis that couples the fluid flow behavior throughout the cardiac cycle with the structural response of the vessel walls. This research focuses on the computational study of a Multiscale Fluid-Structure Interaction on the arterial wall by coupling Finite Volumes Method (FVM) predictions of the Fluid Dynamics within the artery with Finite Elements Method (FEM) predictions of the Elasto-Dynamics response of the arterial walls and 1-D closed loop electrical circuit system to generate the dynamic pressure pulse. To this end, a commercial FVM Computational Fluid Dynamics (CFD) code (STAR-CCM+ 7.09.012) will be coupled through an external interface with a commercial FEM Elasto-Dynamics code (ABAQUS V6.12). The coupling interface is written in such a way that the wall shear stresses and pressures predicted by the CFD analysis will be passed as boundary conditions to the FEM structural solver. The deformations predicted by the FEM structural solver will be passed to the CFD solver to update the geometry in an implicit manner before the following iteration step. The coupling between the FSI and the 1-D closed loop lump parameter circuit updated the pressure pulse and mass flow rates generated by the circuit in an explicit manner after the periodic solution in the FSI analysis had settled. The methodology resulting from this study will be incorporated in a larger collaborative research program between UCF and ORHS that entails optimization of surgical implantation of Left Ventricular Assist Devices (LVAD) cannulae and bypass grafts with the aim to minimize thrombo-embolic events. Moreover, the work proposed will also be applied to another such collaborative project focused on the computational fluid dynamics modeling of the circulation of congenitally affected cardiovascular systems of neonates, specifically the Norwood and Hybrid Norwood circulation of children affected by the hypoplastic left heart syndrome.
Show less - Date Issued
- 2013
- Identifier
- CFE0004753, ucf:49794
- Format
- Document (PDF)
- PURL
- http://purl.flvc.org/ucf/fd/CFE0004753
- Title
- Transient CFD analysis of autorotation using hybrid LES and adaptive mesh morphing techniques.
- Creator
-
Coronado Domenge, Patricia, Das, Tuhin, Kassab, Alain, Kumar, Ranganathan, Leishman, J., Bhattacharya, Aniket, University of Central Florida
- Abstract / Description
-
Large-Eddy Simulation (LES) based turbulence modeling is a developing area of research in Fluid-Structure Interaction (FSI). There is considerable scope for further scientific research in this field and this dissertation aims to extend it to the study of flow-induced motion. The emphasis of this work is on autorotation, an important category of flow-induced motion that is commonly seen in energy applications such as wind turbines and in aviation applications such as the autogyro. In contrast...
Show moreLarge-Eddy Simulation (LES) based turbulence modeling is a developing area of research in Fluid-Structure Interaction (FSI). There is considerable scope for further scientific research in this field and this dissertation aims to extend it to the study of flow-induced motion. The emphasis of this work is on autorotation, an important category of flow-induced motion that is commonly seen in energy applications such as wind turbines and in aviation applications such as the autogyro. In contrast to existing works on FSI that typically assume prescribed motion of structures in a flow field, this research develops LES based FSI studies for large-scale flow-induced motions as seen in autorotation. The uniqueness of the formulation and modeling approach lies in the development of a numerically stable computational scheme that incorporates a moving and morphing mesh structure. The method is first demonstrated for the autorotation of a square flat plate and then extended to a rotor structure similar to that of a helicopter.In order to simulate an autorotating square flat plate, a coupled Computational Fluid Dynamics (CFD) - Rigid Body Dynamics (RBD) model is proposed, employing the delayed-detached-eddy simulation (DDES) and the Smagorinsky turbulence models to resolve subgrid-scale stresses (SGS). The plate is allowed to spin freely about its center of mass. Computational results are compared to experimental measurements and Reynolds Average Navier-Stokes (RANS) simulations found in the literature. When compared to RANS, the results from the LES models provide better predictions of the pressure coefficient. Moreover, LES accurately captures the transient behavior of the plate, and close correspondence is found between the predicted and measured moment coefficients. The qualitative prediction of vortex structures and the quantitative computation of pressure coefficients are in good agreement with experimental results. Hybrid models, such as improved Delayed-Detached-Eddy Simulation (iDDES), are shown to provide very similar results to those of pure LES. Therefore hybrid models are found to be a good alternative to use for the simulation of FSI in autorotation, saving valuable computational time . The iDDES method combines both RANS and LES, dividing the flow domain into LES far away from a solid wall and RANS near a solid wall, overcoming the computational costs of pure LES.Encouraging results from this effort prompted the extension to a realistic scenario, namely the autorotation of a flapping-blade rotor in a prevailing wind field. A coupled CFD - Multi Body Dynamics (MBD) model is developed to study the complex FSI of an autorotating 3-blade rotor, similar to that of a helicopter, employing the iDDES turbulence model. In addition to the rotor being allowed to spin freely about its axis, each of the individual blades is free to rotate about hinges at the root. This adds degrees of freedom to the kinematics of the rotor and necessitates localized mesh morphing around the blades to capture the FSI with accuracy. The model is validated against experimental data and shows excellent agreement. The experimental apparatus consists of a flapping blade rotor and a fixture used to mount it at different angles of incidence with respect to the wind field. The rotor is instrumented with a DC motor that is operated in generator mode. The setup is dual-purpose, providing speed measurement using the motor's back-emf and regenerative braking by varying the current draw. Overall, the presented research can help obtain accurate values of aerodynamic parameters at a high spatial resolution that would be otherwise difficult to acquire in experiments. Ultimately this approach can be a cost effective means of aerodynamic modeling in applications involving large scale FSI.
Show less - Date Issued
- 2016
- Identifier
- CFE0006088, ucf:50952
- Format
- Document (PDF)
- PURL
- http://purl.flvc.org/ucf/fd/CFE0006088
- Title
- A FLUID STRUCTURE INTERACTION MODEL OF INTRACORONARY ATHEROSCLEROTIC PLAQUE RUPTURE.
- Creator
-
TEUMA-MELAGO, Eric, Ilegbusi, Olusegun, University of Central Florida
- Abstract / Description
-
Plaque rupture with superimposed thrombosis is the primary cause of acute coronary syndromes of unstable angina, myocardial infarction and sudden death. Although intensive studies in the past decade have shed light on the mechanism that causes unstable atheroma, none has directly addressed the clinical observation that most myocardial infarction (MI) patients have moderate stenoses (less than 50%). Considering the important role the arterial wall compliance and pulsitile blood flow play in...
Show morePlaque rupture with superimposed thrombosis is the primary cause of acute coronary syndromes of unstable angina, myocardial infarction and sudden death. Although intensive studies in the past decade have shed light on the mechanism that causes unstable atheroma, none has directly addressed the clinical observation that most myocardial infarction (MI) patients have moderate stenoses (less than 50%). Considering the important role the arterial wall compliance and pulsitile blood flow play in atheroma rupture, fluid-structure interaction (FSI) phenomenon has been of interest in recent studies. In this thesis, the impact is investigated numerically of coupled blood flow and structural dynamics on coronary plaque rupture. The objective is to determine a unique index that can be used to characterize plaque rupture potential. The FSI index, developed in this study for the first time derives from the theory of buckling of thin-walled cylinder subjected to radial pressure. Several FSI indices are first defined by normalizing the predicted hemodynamic endothelial shear stress by the structural stresses, specifically, by the maximum principal stress (giving the ratio ), and the Von Mises stress (giving the ratio ). The predicted at the location of maximum (i.e { }) denoted , is then chosen to characterize plaque rupture through systematic investigation of a variety of plaque characteristics and simulated patient conditions. The conditions investigated include varying stenosis levels ranging from 20% to 70%, blood pressure drop ranging from 3125 Pa/m to 9375 Pa/m, fibrous cap thickness ranging from to , lipid pool location ranging from the leading to the trailing edge of plaque, lipid pool volume relative to stenosis volume ranging from 24% to 80%, Calcium volume relative to stenosis volume ranging from 24% to 80% and arterial remodeling. The predicted varies with the stenosis severity and indicates that the plaques investigated are prone to rupture at approximately 40-45% stenosis levels. It predicts that high pressure significantly lowers the threshold stenosis rate for plaque rupture. In addition, the plaque potential to rupture increases for relatively thin fibrous cap, lipid core located near the leading plaque shoulder, and dramatically for relative lipid pool volume above 60%. However, calcium deposit has marginal effect on plaque rupture. Overall, the predicted results are consistent with clinical observations, indicating that the has the potential to characterize plaque rupture when properly established. In the appendix, the unsteady flow in a collapsible tube model of a diseased artery is solved analytically. The novelty of our approach is that the set of governing equations is reduced to a single integro-differential equation in the transient state. The equation was solved using the finite difference method to obtain the pressure and compliant wall behavior. The analytical approach is less computer-intensive than solving the full set of governing equations. The predicted membrane deflection is quite large at low inlet velocity, suggesting possible approach to breakdown in equilibrium. As the transmural pressure increases with wall deflection, bulges appear at the ends of the membrane indicating critical stage of stability, consistent with previous studies. An increase in wall thickness reduces the wall deflection and ultimately results in its collapse. The collapse is due to breakdown in the balance of wall governing equation. An increase in internal pressure is required to maintain membrane stability.
Show less - Date Issued
- 2006
- Identifier
- CFE0001471, ucf:47084
- Format
- Document (PDF)
- PURL
- http://purl.flvc.org/ucf/fd/CFE0001471
- Title
- Multi-scale fluid-structure interaction model analysis of patient-specific geometry for optimization of lvad outflow graft implantation: an investigation aimed at reducing stroke risk.
- Creator
-
Prather, Ray, Kassab, Alain, Mansy, Hansen, Bai, Yuanli, Divo, Eduardo, DeCampli, William, University of Central Florida
- Abstract / Description
-
A Left Ventricular Assist Device (LVAD), is a mechanical pump capable of(&)nbsp;providing circulatory myocardium relief when used as bridge-to-transplantation by reducing the workload of a failing heart, with the additional bonus of allowing for cardiac recovery when used as destination therapy. The newer generations of continuous flow VADs are essentially axial or radial flow pumps, and while these devices are capable their efficiency depends upon fluid composition and flow field patterns....
Show moreA Left Ventricular Assist Device (LVAD), is a mechanical pump capable of(&)nbsp;providing circulatory myocardium relief when used as bridge-to-transplantation by reducing the workload of a failing heart, with the additional bonus of allowing for cardiac recovery when used as destination therapy. The newer generations of continuous flow VADs are essentially axial or radial flow pumps, and while these devices are capable their efficiency depends upon fluid composition and flow field patterns. The most devastating complication of VAD therapy is caused by embolization of thrombi formed within the LVAD or inside the heart into the brain leading to stroke. Anticoagulation management and improved LVADs design has reduced stroke incidence, however, investigators have recently reported the incidence of thromboembolic cerebral events is still significant and ranges from 14% to 47% over a period of 6-12 months. Blood clots may cause obstruction of critical vessels, such as cerebral arteries, reducing brain oxygenation and resulting in devastating consequences like major neurocognitive malfunction and complications which can be fatal.The hypothesis that incidence of stroke can be significantly reduced by adjusting the VAD outflow cannula implantation to direct dislodged thrombi away from the cerebral vessels has been recently supported by a series of steady flow computations assuming rigid vessel walls for the vasculature. Such studies have shown as much as a 50% reduction in embolization rates depending on outflow cannula implantation. In this study, a pulsatile fully compliant vessel wall model is developed to further establish this hypothesis. A time-dependent multi-scale Eulerian Computational Fluid Dynamics (CFD) analysis of patient-specific geometry models of the VAD-bed vasculature is coupled with a 3D Finite Element Analysis (FEA) of the mechanical response of the vascular walls to establish the VAD assisted hemodynamics. A Lagrangian particle tracking algorithm is used to determine the embolization rates of thrombi emanating from the cannula or other possible thrombogenic locations such as the aortic root. This multiscale Eulerian-Lagrangian pulsatile fluid-structure coupled paradigm allows for a fully realistic model of the hemodynamics of interest. The patient-specific geometries obtained from CT scan are implemented into the numerical domain in two modes. In the 3D CFD portion of the problem, the geometry accounts solely for the flow volume where the fluid is modelled as constant density and non-Newtonian under laminar pulsatile flow conditions. The blood-thrombus ensemble in treated as a two-phase flow, handled by an Eulerian-Lagrangian coupled scheme to solve the flow field and track particle transport. Thrombi are modelled as constant density spherical particles. Particle interactions are limited to particle-to-wall and particle-to-fluid, while particle-to-particle interaction are neglected for statistical purposes. On the other hand, with the help of Computer Aided Design (CAD) software a patient-specific aortic wall geometry with variable wall thickness is brought into the numerical domain. FEA is applied to determine the aortic wall cyclic displacement under hydrodynamic loads. To properly account for wall deformation, the arterial wall tissue incorporates a hyperelastic material model based on the anisotropic Holzapfel model for arteries. This paradigm is referred to as Fluid Structure Interaction (FSI) and allows structural analysis in conjunction with flow investigation to further monitor pathological flow patterns. The FSI model is driven by time dependent flow and pressure boundary conditions imposed at the boundaries of the 3D computational domain through a 50 degree of freedom 0D lumped parameter model (LPM) electric circuit analog of the peripheral VAD-assisted circulation.Results are presented for a simple vessel model of the ascending aorta to validate the anisotropic fiber orientation implementation. Arterial wall dilation is measured between 5-20% in the range reported in literature. Hemodynamics of the VAD assisted flow in a patient-derived geometry computed using rigid vessels walls are compared to those for a linearly elastic vessel wall model and a hyperelastic anisotropic vessel wall model. Moreover, the thromboembolization rates are presented and compared for pulsatile hemodynamics in rigid and compliant wall models. Pulsatile flow solutions for embolization probabilities corroborate the hypothesis that tailoring the LVAD cannula implantation configuration can significantly reduce thromboembolization rates, and this is consistent with indications from previous steady-flow calculations.
Show less - Date Issued
- 2018
- Identifier
- CFE0007077, ucf:52017
- Format
- Document (PDF)
- PURL
- http://purl.flvc.org/ucf/fd/CFE0007077