A partition of unity approach to fluid mechanics and fluid-structure interaction.
Ontology highlight
ABSTRACT:
SUBMITTER: Balmus M
PROVIDER: S-EPMC7610902 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Congenital bicuspid aortic valve (BAV) consists of two fused cusps and represents a major risk factor for calcific valvular stenosis. Herein, a fully coupled fluid-structure interaction (FSI) BAV model was developed from patient-specific magnetic resonance imaging (MRI) and compared against in vivo 4-dimensional flow MRI (4D Flow). FSI simulation compared well with 4D Flow, confirming direction and magnitude of the flow jet impinging onto the aortic wall as well as location and extension of secondary flows and vortices developing at systole: the systolic flow jet originating from an elliptical 1.6 cm2 orifice reached a peak velocity of 252.2 cm/s, 0.6% lower than 4D Flow, progressively impinging on the ascending aorta convexity. The FSI model predicted a peak flow rate of 22.4 L/min, 6.7% higher than 4D Flow, and provided BAV leaflets mechanical and flow-induced shear stresses, not directly attainable from MRI. At systole, the ventricular side of the non-fused leaflet revealed the highest wall shear stress (WSS) average magnitude, up to 14.6 Pa along the free margin, with WSS progressively decreasing towards the belly. During diastole, the aortic side of the fused leaflet exhibited the highest diastolic maximum principal stress, up to 322 kPa within the attachment region. Systematic comparison with ground-truth non-invasive MRI can improve the computational model ability to reproduce native BAV hemodynamics and biomechanical response more realistically, and shed light on their role in BAV patients' risk for developing complications; this approach may further contribute to the validation of advanced FSI simulations designed to assess BAV biomechanics.
| S-EPMC7990108 | biostudies-literature
Project description:PurposeTo develop a physical understanding of ureterorenoscopy irrigation, we derive mathematical models from basic physical principles and compare these predictions with the results of benchtop experiments. Mathematical modeling can be used to understand the role of inlet pressure, tip deflection, the presence of working tools, geometric properties of the instruments used, and material properties of the irrigation fluid on resulting flow rate.Materials and methodsWe develop theoretical models to describe irrigation flow in an idealized setup and compare with benchtop experiments for flow through a straight scope, a scope with a deflected tip, and a scope with a working tool inserted. The benchtop experiments were performed using Boston Scientific LithoVue ureteroscope and a variety of Boston Scientific working tools. Standard ureteroscope working channels have circular cross sections, but using theoretical models we investigate whether modifications to the cross-sectional geometry can enhance flow rates.ResultsThe theoretical flow predictions are confirmed by experimental results. Tip deflection is shown to have a negligible effect on flow rate, but the presence of working tools decreases flow significantly (for a fixed driving pressure). Flow rate is predicted to improve when tools are placed at the edge of the channel, rather than the center, and modifying the cross-sectional shape from a circle to an ellipse can further increase flow rate.ConclusionsA mathematical framework is formulated and shown to accurately predict the properties of ureteroscope irrigation flow. The theoretical approach has significant potential in quantifying irrigation flow and improving ureteroscope design.
| S-EPMC6352503 | biostudies-literature
Project description:The exact mechanism to orchestrate the action of hundreds of dynein motor proteins to generate wave-like ciliary beating remains puzzling and has fascinated many scientists. We present a 3D model of a cilium and the simulation of its beating in a fluid environment. The model cilium obeys a simple geometric constraint that arises naturally from the microscopic structure of a real cilium. This constraint allows us to determine the whole 3D structure at any instant in terms of the configuration of a single space curve. The tensions of active links, which model the dynein motor proteins, follow a postulated dynamical law, and together with the passive elasticity of microtubules, this dynamical law is responsible for the ciliary motions. In particular, our postulated tension dynamics lead to the instability of a symmetrical steady state, in which the cilium is straight and its active links are under equal tensions. The result of this instability is a stable, wave-like, limit cycle oscillation. We have also investigated the fluid-structure interaction of cilia using the immersed boundary (IB) method. In this setting, we see not only coordination within a single cilium but also, coordinated motion, in which multiple cilia in an array organize their beating to pump fluid, in particular by breaking phase synchronization.
| S-EPMC5924875 | biostudies-other
Project description:Physical properties of primary cilia membranes in living cells were examined using two independent, high-spatiotemporal-resolution approaches: fast tracking of single quantum dot-labeled G protein-coupled receptors and a novel two-photon super-resolution fluorescence recovery after photobleaching of protein ensemble. Both approaches demonstrated the cilium membrane to be partitioned into corralled domains spanning 274 ± 20 nm, within which the receptors are transiently confined for 0.71 ± 0.09 s. The mean membrane diffusion coefficient within the corrals, Dm1 = 2.9 ± 0.41 µm2/s, showed that the ciliary membranes were among the most fluid encountered. At longer times, the apparent membrane diffusion coefficient, Dm2 = 0.23 ± 0.05 µm2/s, showed that corral boundaries impeded receptor diffusion 13-fold. Mathematical simulations predict the probability of G protein-coupled receptors crossing corral boundaries to be 1 in 472. Remarkably, latrunculin A, cytochalasin D, and jasplakinolide treatments altered the corral permeability. Ciliary membranes are thus partitioned into highly fluid membrane nanodomains that are delimited by filamentous actin.
| S-EPMC6080922 | biostudies-literature
Project description:Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.
| S-EPMC7760396 | biostudies-literature
Project description:In this paper a fluid-structure interaction (FSI) experiment is presented. The aim of this experiment is to provide a challenging yet easy-to-setup FSI test case that addresses the need for rigorous testing of FSI algorithms and modeling frameworks. Steady-state and periodic steady-state test cases with constant and periodic inflow were established. Focus of the experiment is on biomedical engineering applications with flow being in the laminar regime with Reynolds numbers 1283 and 651. Flow and solid domains were defined using computer-aided design (CAD) tools. The experimental design aimed at providing a straightforward boundary condition definition. Material parameters and mechanical response of a moderately viscous Newtonian fluid and a nonlinear incompressible solid were experimentally determined. A comprehensive data set was acquired by using magnetic resonance imaging to record the interaction between the fluid and the solid, quantifying flow and solid motion.
| S-EPMC5600002 | biostudies-literature
Project description:Recirculation in venovenous extracorporeal membrane oxygenation (VV ECMO) leads to reduction in gas transfer efficiency. Studies of the factors contributing have been performed using in vivo studies and computational models. The fixed geometry of previous computational models limits the accuracy of results. We have developed a finite element computational fluid dynamics model incorporating fluid-structure interaction (FSI) that incorporates atrial deformation during atrial filling and emptying, with fluid flow solved using large eddy simulation. With this model, we have evaluated an extensive number of factors that could influence recirculation during two-site VV ECMO, and characterized their impact on recirculation, including cannula construction, insertion depth and orientation, VV ECMO configuration, circuit blood flow, and changes in volume, venous return, heart rate, and blood viscosity. Simulations revealed that extracorporeal blood flow relative to cardiac output, ratio of superior vena caval (SVC) to inferior vena caval (IVC) blood flow, position of the SVC cannula relative to the cavo-atrial junction, and orientation of the return cannula relative to the tricuspid valve had major influences (>20%) on recirculation fraction. Factors with a moderate influence on recirculation fraction (5%-20%) include heart rate, return cannula diameter, and direction of extracorporeal flow. Minimal influence on recirculation (<5%) was associated with atrial volume, position of the IVC cannula relative to the cavo-atrial junction, the number of side holes in the return cannula, and blood viscosity.
| S-EPMC8318564 | biostudies-literature
Project description:Untreated tricuspid valve regurgitation (TR) is associated with increased rates of mortality, morbidity, and hospitalization. Current pharmacological and surgical treatment options for TR are limited. MitraClip (MC), an edge-to-edge percutaneous intervention, has been reported to be effective for treatment of TR. The goal of this study was to examine the effects of MC position on TR, using a multiphysics fluid-structure-interaction (FSI) analysis. The computational set up included the tricuspid valve (TV), the chordae tendineae, the blood particles, and a tube that surrounded the leaflets and blood particles. The leaflets and chordae were modeled as hyperelastic materials, and blood was modeled using smoothed particle hydrodynamics. FSI analysis was conducted for blood flow through the closed valve for multiple simulations that account for normal, diseased, and treated conditions of the TV. To simulate the diseased TV, a group of chordae between septal and pulmonary leaflets were removed from the normal TV, which produced increased regurgitation. Four MC treated scenarios were considered: i) one MC near the annulus, ii) one MC approximately midway between the annulus and leaflet tip, iii) one MC near the leaflet tip, iv) two MCs: one approximately midway between the annulus and leaflet tip, and one close to the leaflet tip. The TR increased in diseased TV (7.5%) compared to normal TV (2.5%). All MC treated scenarios decreased TR. The MC located near the midway point between the annulus and leaflet tip led to largest decrease in TR (75.2% compared to the untreated condition). The MC located near the leaflet tip was associated with lowest reduction in TR (2.2% compared to the untreated condition). When two MCs were used, reduction in TR was relatively high (68.7%), but TR was not improved compared to the optimal single MC. MC caused high stresses in the vicinity of the clipping area in all conditions; the highest occurred when the MC was near the leaflet tips. Using a quantitative computational approach, we confirm previous clinical reports on the efficacy of MC for treatment of TR. The results of this study could lead to the design of more efficient MC interventions for TR.
| S-EPMC6706066 | biostudies-literature
Project description:A three dimensional viscous finite element model is presented in this paper for the analysis of the acoustic fluid structure interaction systems including, but not limited to, the cochlear-based transducers. The model consists of a three dimensional viscous acoustic fluid medium interacting with a two dimensional flat structure domain. The fluid field is governed by the linearized Navier-Stokes equation with the fluid displacements and the pressure chosen as independent variables. The mixed displacement/pressure based formulation is used in the fluid field in order to alleviate the locking in the nearly incompressible fluid. The structure is modeled as a Mindlin plate with or without residual stress. The Hinton-Huang's 9-noded Lagrangian plate element is chosen in order to be compatible with 27/4 u/p fluid elements. The results from the full 3d FEM model are in good agreement with experimental results and other FEM results including Beltman's thin film viscoacoustic element [2] and two and half dimensional inviscid elements [21]. Although it is computationally expensive, it provides a benchmark solution for other numerical models or approximations to compare to besides experiments and it is capable of modeling any irregular geometries and material properties while other numerical models may not be applicable.
| S-EPMC2824573 | biostudies-literature
Project description:Vascular cross-clamping is applied in many cardiovascular surgeries such as coronary bypass, aorta repair and valve procedures. Experimental studies have found that clamping of various degrees caused damage to arteries. This study examines the effects of popular clamps on vessel wall. Models of the aorta and clamp were created in Computer Assisted Design and Finite Element Analysis packages. The vessel wall was considered as a non-linear anisotropic material while the fluid was simulated as Newtonian with pulsatile flow. The clamp was applied through displacement time function. Fully coupled two-way solid-fluid interaction models were developed. It was found that the clamp design significantly affected the stresses in vessel wall. The clamp with a protrusion feature increased the overall Von Mises stress by about 60% and the compressive stress by more than 200%. Interestingly, when the protrusion clamp was applied, the Von Mises stress at the lumen (endothelium) side of artery wall was about twice that of the outer wall. This ratio was much higher than that of the plate-like clamp which was about 1.3. The flow reversal process was demonstrated during clamping. Vibrations, flow and wall shear stress oscillations were detected immediately before total vessel occlusion. The commonly used protrusion clamp increased stresses in vessel wall, especially the compressive stress. This design also significantly increased the stresses on endothelium, detrimental to vessel health. The present findings are relevant to surgical clamp design as well as the transient mechanical loading on the endothelium and potential injury. The deformation and stress analysis may provide valuable insights into the mode of tissue injury during cross-clamping.
| S-EPMC2813350 | biostudies-literature