Cardiovascular disease (CVD) is the leading cause of mortality worldwide and its incidence is rising due to an aging population. The development and progression of CVD is directly linked to adverse vascular hemodynami...
详细信息
Cardiovascular disease (CVD) is the leading cause of mortality worldwide and its incidence is rising due to an aging population. The development and progression of CVD is directly linked to adverse vascular hemodynamics and biomechanics, whose in-vivo measurement remains challenging but can be simulated numerically and experimentally. The ability to evaluate these parameters in patient-specific CVD cases is crucial to better predict future disease progression, risk of adverse events, and treatment efficacy. While significant progress has been made toward patient-specific hemodynamic simulations, blood vessels are often assumed to be rigid, which does not consider the compliant mechanical properties of vessels whose malfunction is implicated in disease. In an effort to simulate the biomechanics of flexible vessels, fluid-structure interaction (FSI) simulations have emerged as promising tools for the characterization of hemodynamics within patient-specific cardiovascular anatomies. Since FSI simulations combine the blood's fluid domain with the arterial structural domain, they pose novel challenges for their experimental validation. This paper reviews the scientific work related to FSI simulations for patient-specific arterial geometries and the current standard of FSI model validation including the use of compliant arterial phantoms, which offer novel potential for the experimental validation of FSI results.
A series of fluidstructureinteraction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structu...
详细信息
A series of fluidstructureinteraction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluidstructureinteraction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure.
Non-Darcian effects on natural convective flow and heat transfer in a square enclosure filled with a porous medium was analyzed numerically using fluid-structure interaction (FSI) model. The transport equations were s...
详细信息
Non-Darcian effects on natural convective flow and heat transfer in a square enclosure filled with a porous medium was analyzed numerically using fluid-structure interaction (FSI) model. The transport equations were solved for various pertinent parameters using a finite element formulation based on the Galerkin method of weighted residuals. Such parameters included Rayleigh number, porosity, elasticity of the flexible wall, and the effective thermal conductivity of the porous medium. Further, the fluid domain was described by an Arbitrary-Lagrangian-Eulerian (ALE) formulation that is fully coupled to the structure domain. Different flow models for porous media such as Darcy's law model and Darcy-Forchheimer model were considered in this investigation. Comparisons of isotherms, streamlines, and average Nusselt number are made between rigid and FSI models. The results of this investigation showed that Rayleigh number and the elasticity of the flexible wall had a profound effect on the shape of the flexible wall and consequently on the heat transfer enhancement within the enclosure. (C) 2012 Elsevier Ltd. All rights reserved.
In this article, we present a partitioned procedure for fluid-structure interaction problems in which contacts among different deformable bodies can occur. A typical situation is the movement of a thin valve (e.g. the...
详细信息
In this article, we present a partitioned procedure for fluid-structure interaction problems in which contacts among different deformable bodies can occur. A typical situation is the movement of a thin valve (e.g. the aortic valve) immersed in an incompressible viscous fluid (e.g. the blood). In the proposed strategy the fluid and structure solvers are considered as independent "black-boxes" that exchange forces and displacements;the structure solvers are moreover not supposed to manage contact by themselves. The hypothesis of non-penetration among solid objects defines a non-convex optimization problem. To solve the latter, we use an internal approximation algorithm that is able to directly handle the cases of thin structures and self-contacts. A numerical simulation on an idealized aortic valve is finally realized with the aim of illustrating the proposed scheme. (C) 2008 Elsevier B.V. All rights reserved.
In this paper a problem belonging to the moving boundary class is tackled with a 2-D application of computational fluid dynamics techniques. The motion of an isolated rigid particle freely suspended in an incompressib...
详细信息
In this paper a problem belonging to the moving boundary class is tackled with a 2-D application of computational fluid dynamics techniques. The motion of an isolated rigid particle freely suspended in an incompressible Newtonian fluid in a narrow channel is studied numerically at a low Reynolds number, yet different from zero. The actual problem consists of two coupled problems: the motion of the viscous fluid and that of the rigid particle suspended and convected with the fluid. The full Navier-Stokes equations (i.e. both transient and convective terms are included) are solved in the fluid domain by means of the finite element method while the motion of the particle is determined on the basis of a rigid act of motion. Results from simulations corresponding to differential initial positions of the particle are shown in this paper: they allow one to study the rotational motions of the particle as well as its displacements. The goal of the paper is to analyse the lateral displacement behaviour of the particle, already observed in experimental studies in microcirculation. In particular, lateral migrations are supposed to be due to inertial forces acting in the fluid around the moving particle combined with the proximity of the resting wall (wall effect). Preliminary results are in fairly good agreement with those available in the literature.
In this article we discuss the application of a Lagrange multiplier based fictitious domain method for the simulation of the motion of two rigid flaps in an unsteady flow generated by pressure gradients. The distribut...
详细信息
In this article we discuss the application of a Lagrange multiplier based fictitious domain method for the simulation of the motion of two rigid flaps in an unsteady flow generated by pressure gradients. The distributed Lagrange multiplier technique can be an important numerical tool to design a mechanical heart valve and investigate the flow around rigid flaps without assuming the motion of the flaps in advance. Here, we derive a mathematical formulation of a fluid-structure interaction model that includes the generalized Neumann boundary conditions on the upstream and downstream boundaries along with rigid flaps rotating around the fixed points. The solution method includes the finite element approximation for space and the Marchuk-Yanenko operator splitting scheme for time discretization. This study presents the numerical results obtained for flap motion for a simple sinusoidal wave. Furthermore, these simulations are extended to apply to a more complex biological system involving the systolic phase of the pulse pressure. (c) 2007 Elsevier Ltd. All rights reserved.
The ascending branch of the aorta is one of the most stressed organ of the arterial system. We aim to design a biomechanical model for analysing the aorta dynamics under a shock. The model includes the aorta layers an...
详细信息
The ascending branch of the aorta is one of the most stressed organ of the arterial system. We aim to design a biomechanical model for analysing the aorta dynamics under a shock. The model includes the aorta layers and the influence of the blood pressure. We undertake a three-dimensional modal analysis of the coupled aorta-blood system. We determine in the present work the coupled natural frequencies and the modes shapes of the system of the aorta and blood. Three models are presented in this study: three-layers model, two-layers model and one layer model. For the analytical solving a potential technique is used to obtain a general solution for an aorta domain. The finite element model is then validated by these original analytical solutions. The results from the proposed method are in good agreement with numerical solutions. (C) 2014 Elsevier Ltd. All rights reserved.
Calcific aortic valve disease (CAVD) is characterized by calcification accumulation and thickening of the aortic valve cusps, leading to stenosis. The importance of fluid flow shear stress in the initiation and regula...
详细信息
Calcific aortic valve disease (CAVD) is characterized by calcification accumulation and thickening of the aortic valve cusps, leading to stenosis. The importance of fluid flow shear stress in the initiation and regulation of CAVD progression is well known and has been studied recently using fluid-structure interaction (FSI) models. While cusp calcifications are three-dimensional (3D) masses, previously published FSI models have represented them as either stiffened or thickened two-dimensional (2D) cusps. This study investigates the hemodynamic effect of these calcifications employing FSI models using 3D patient-specific calcification masses. A new reverse calcification technique (RCT) is used for modeling different stages of calcification growth based on the spatial distribution of calcification density. The RCT is applied to generate the 3D calcification deposits reconstructed from a patient-specific CT scans. Our results showed that consideration of 3D calcification deposits led to both higher fluid shear stresses and unique fluid shear stress distribution on the aortic side of the cusps that may have an impact on the calcification growth rate. However, the flow did not seem to affect the geometry of the calcification during the growth phase.
Purpose The purpose of this study is to model the dynamic characteristics of an opened supersonic disk-gap-band parachute. Design/methodology/approach A fluid-structure interaction (FSI) method with body-fitted mesh i...
详细信息
Purpose The purpose of this study is to model the dynamic characteristics of an opened supersonic disk-gap-band parachute. Design/methodology/approach A fluid-structure interaction (FSI) method with body-fitted mesh is used to simulate the supersonic parachute. The compressible flow is modeled using large-eddy simulation (LES). A contact algorithm based on the penalty function with a virtual contact domain is proposed to solve the negative volume problem of the body-fitted mesh. Automatic unstructured mesh generation and automatic mesh moving schemes are used to handle complex deformations of the canopy. Findings The opened disk-gap-band parachute is simulated using Mach 2.0, and the simulation results fit well with the wind tunnel test data. It is found that the LES model can successfully predict large-scale turbulent vortex in the flow. This study also demonstrates the capability of the present FSI method as a tool to predict shock oscillation and breathing phenomenon of the canopy. Originality/value The contact algorithm based on the penalty function with a virtual contact domain is proposed for the first time. This methodology can be used to solve the negative volume problem of the dynamic mesh in the flow field.
In this paper we present our aerodynamics and fluid-structure interaction (FSI) computational techniques that enable dynamic, fully coupled, 3D FSI simulation of wind turbines at full scale, and in the presence of the...
详细信息
In this paper we present our aerodynamics and fluid-structure interaction (FSI) computational techniques that enable dynamic, fully coupled, 3D FSI simulation of wind turbines at full scale, and in the presence of the nacelle and tower (i.e., simulation of the "full machine"). For the interaction of wind and flexible blades we employ a nonmatching interface discretization approach, where the aerodynamics is computed using a low-order finite-element-based ALE-VMS technique, while the rotor blades are modeled as thin composite shells discretized using NURBS-based isogeometric analysis (IGA). We find that coupling FEM and IGA in this manner gives a good combination of efficiency, accuracy, and flexibility of the computational procedures for wind turbine FSI. The interaction between the rotor and tower is handled using a non-overlapping sliding-interface approach, where both moving-and stationary-domain formulations of aerodynamics are employed. At the fluid-structure and sliding interfaces, the kinematic and traction continuity is enforced weakly, which is a key ingredient of the proposed numerical methodology. We present several simulations of a three-blade 5 MW wind turbine, with and without the tower. We find that, in the case of no tower, the presence of the sliding interface has no effect on the prediction of aerodynamic loads on the rotor. From this we conclude that weak enforcement of the kinematics gives just as accurate results as the strong enforcement, and thus enables the simulation of rotor-tower interaction (as well as other applications involving mechanical components in relative motion). We also find that the blade passing the tower produces a 10-12% drop (per blade) in the aerodynamic torque. We feel this finding may be important when it comes to the fatigue-life analysis and prediction for wind turbine blades.
暂无评论