FSI modeling with the DSD/SST method for the fluid and finite difference method for the structure

We present a fluid–structure interaction (FSI) modeling method based on using the deforming-spatial-domain/stabilized space–time (DSD/SST) method for the fluid mechanics part and a finite difference (FD) method for the structural mechanics part. As the structural mechanics model, we focus on the thin-shell model. The fluid mechanics equations with moving boundaries are solved with the DSD/SST method and the thin-shell structural mechanics equation is solved with a FD method, with partitioned coupling between the two parts. The coupling of the DSD/SST and FD solvers makes sure that the boundary conditions on the fluid-structure interface at the end of each time step are matched between the fluid and the structure. A hanging plate in vacuum under gravitational force is performed to validate the structure solver. In addition, a pitching plate in a uniform flow is simulated to validate the FSI solver. The present results are in reasonable agreement with data predicted by other methods.     

Solution of linear systems in arterial fluid mechanics computations with boundary layer mesh refinement

Computation of incompressible flows in arterial fluid mechanics, especially because it involves fluid–structure interaction, poses significant numerical challenges. Iterative solution of the fluid mechanics part of the equation systems involved is one of those challenges, and we address that in this paper, with the added complication of having boundary layer mesh refinement with thin layers of elements near the arterial wall. As test case, we use matrix data from stabilized finite element computation of a bifurcating middle cerebral artery segment with aneurysm. It is well known that solving linear systems that arise in incompressible flow computations consume most of the time required by such simulations. For solving these large sparse nonsymmetric systems, we present effective preconditioning techniques appropriate for different stages of the computation over a cardiac cycle.     

Wall shear stress calculations in space–time finite element computation of arterial fluid–structure interactions

The stabilized space–time fluid–structure interaction (SSTFSI) technique was applied to arterial FSI problems soon after its development by the Team for Advanced Flow Simulation and Modeling. The SSTFSI technique is based on the Deforming-Spatial-Domain/Stabilized Space–Time (DSD/SST) formulation and is supplemented with a number of special techniques developed for arterial FSI. The special techniques developed in the recent past include a recipe for pre-FSI computations that improve the convergence of the FSI computations, using an estimated zero-pressure arterial geometry, Sequentially Coupled Arterial FSI technique, using layers of refined fluid mechanics mesh near the arterial walls, and a special mapping technique for specifying the velocity profile at inflow boundaries with non-circular shape. In this paper we introduce some additional special techniques, related to the projection of fluid–structure interface stresses, calculation of the wall shear stress (WSS), and calculation of the oscillatory shear index. In the test computations reported here, we focus on WSS calculations in FSI modeling of a patient-specific middle cerebral artery segment with aneurysm. Two different structural mechanics meshes and three different fluid mechanics meshes are tested to investigate the influence of mesh refinement on the WSS calculations.     

FSI analysis of the blood flow and geometrical characteristics in the thoracic aorta

We present a fluid–structure interaction (FSI) analysis of the blood flow and geometrical characteristics in the thoracic aorta. The FSI is handled with the sequentially-coupled arterial FSI technique. The fluid mechanics equations are solved with the ST-VMS method, which is the variational multiscale version of the deforming-spatial-domain/stabilized space–time (DSD/SST) method. We focus on the relationship between the centerline geometry of the aorta and the flow field, which influences the wall shear stress distribution. The centerlines of the aorta models we use in our analysis are extracted from the CT scans, and we assume a constant diameter. Torsion-free model geometries are generated by projecting the original centerline to its averaged plane of curvature. The flow fields for the original and projected geometries are compared to examine the influence of the torsion.     

Multiscale space–time fluid–structure interaction techniques

We present the multiscale space–time techniques we have developed for fluid–structure interaction (FSI) computations. Some of these techniques are multiscale in the way the time integration is performed (i.e. temporally multiscale), some are multiscale in the way the spatial discretization is done (i.e. spatially multiscale), and some are in the context of the sequentially-coupled FSI (SCFSI) techniques developed by the Team for Advanced Flow Simulation and Modeling (T \bigstar AFSM){({\rm T} \bigstar {\rm AFSM})}. In the multiscale SCFSI technique, the FSI computational effort is reduced at the stage we do not need it and the accuracy of the fluid mechanics (or structural mechanics) computation is increased at the stage we need accurate, detailed flow (or structure) computation. As ways of increasing the computational accuracy when or where needed, and beyond just increasing the mesh refinement or decreasing the time-step size, we propose switching to more accurate versions of the Deforming-Spatial-Domain/Stabilized Space–Time (DSD/SST) formulation, using more polynomial power for the basis functions of the spatial discretization or time integration, and using an advanced turbulence model. Specifically, for more polynomial power in time integration, we propose to use NURBS, and as an advanced turbulence model to be used with the DSD/SST formulation, we introduce a space–time version of the residual-based variational multiscale method. We present a number of test computations showing the performance of the multiscale space–time techniques we are proposing. We also present a stability and accuracy analysis for the higher-accuracy versions of the DSD/SST formulation.     

A computational approach for predicting the hydroelasticity of flexible structures based on the pressure Poisson equation

This work is concerned with the modeling of the interaction of fluid flow with flexible solid structures. The flow under consideration is governed by the Navier–Stokes equations for incompressible viscous fluids and modeled with low‐order velocity–pressure finite elements. The motion of the fluid domain is accounted for by the arbitrary Lagrangian–Eulerian formulation. The structure is represented by means of an appropriate standard finite element formulation. The spring smooth analogy is used to mesh control. The time integrating algorithm is based on the predictor–multi‐corrector algorithm. An important aspect of the present work is the introduction of a new monolithic approach based on the fluid pressure Poisson equation (PPE) to solve the hydroelasticity problem of an incompressible viscous fluid with an elastic body that is vibrating due to flow excitation. The PPE is derived to be consistent with the coupled system equation for the fluid–structure interaction (FSI). Based on this approach, an efficient monolithic method is adopted to simulate hydroelasticity between the flexible structure and the flow. The fluid pressure is implicitly derived to satisfy the incompressibility constraint, and the other unknown variables are explicitly derived. The coefficient matrix of the PPE for the FSI becomes symmetric and positive definite. To demonstrate the performance of the proposed approach, two working examples, a beam immersed in incompressible fluid and a guide vane of a Francis turbine passage, were used. The results show the validity of the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Shape optimization of pulsatile ventricular assist devices using FSI to minimize thrombotic risk

In this paper we perform shape optimization of a pediatric pulsatile ventricular assist device (PVAD). The device simulation is carried out using fluid–structure interaction (FSI) modeling techniques within a computational framework that combines FEM for fluid mechanics and isogeometric analysis for structural mechanics modeling. The PVAD FSI simulations are performed under realistic conditions (i.e., flow speeds, pressure levels, boundary conditions, etc.), and account for the interaction of air, blood, and a thin structural membrane separating the two fluid subdomains. The shape optimization study is designed to reduce thrombotic risk, a major clinical problem in PVADs. Thrombotic risk is quantified in terms of particle residence time in the device blood chamber. Methods to compute particle residence time in the context of moving spatial domains are presented in a companion paper published in the same issue (Comput Mech, doi:10.1007/s00466-013-0931-y, 2013). The surrogate management framework, a derivative-free pattern search optimization method that relies on surrogates for increased efficiency, is employed in this work. For the optimization study shown here, particle residence time is used to define a suitable cost or objective function, while four adjustable design optimization parameters are used to define the device geometry. The FSI-based optimization framework is implemented in a parallel computing environment, and deployed with minimal user intervention. Using five SEARCH/POLL steps the optimization scheme identifies a PVAD design with significantly better throughput efficiency than the original device.     

Multiscale sequentially-coupled arterial FSI technique

Multiscale versions of the Sequentially-Coupled Arterial Fluid–Structure Interaction (SCAFSI) technique are presented. The SCAFSI technique was introduced as an approximate FSI approach in arterial fluid mechanics. It is based on the assumption that the arterial deformation during a cardiac cycle is driven mostly by the blood pressure. First we compute a “reference” arterial deformation as a function of time, driven only by the blood pressure profile of the cardiac cycle. Then we compute a sequence of updates involving mesh motion, fluid dynamics calculations, and recomputing the arterial deformation. The SCAFSI technique was developed and tested in conjunction with the stabilized space–time FSI (SSTFSI) technique. Beyond providing a computationally more economical alternative to the fully coupled arterial FSI approach, the SCAFSI technique brings additional flexibility, such as being able to carry out the computations in a spatially or temporally multiscale fashion. In the test computations reported here for the spatially multiscale versions of the SCAFSI technique, we focus on a patient-specific middle cerebral artery segment with aneurysm, where the arterial geometry is based on computed tomography images. The arterial structure is modeled with the continuum element made of hyperelastic (Fung) material.     

Fluid–structure interaction analysis of bioprosthetic heart valves: significance of arterial wall deformation

We propose a framework that combines variational immersed-boundary and arbitrary Lagrangian–Eulerian methods for fluid–structure interaction (FSI) simulation of a bioprosthetic heart valve implanted in an artery that is allowed to deform in the model. We find that the variational immersed-boundary method for FSI remains robust and effective for heart valve analysis when the background fluid mesh undergoes deformations corresponding to the expansion and contraction of the elastic artery. Furthermore, the computations presented in this work show that the arterial wall deformation contributes significantly to the realism of the simulation results, leading to flow rates and valve motions that more closely resemble those observed in practice.     

Fluid–structure interaction modeling of ringsail parachutes with disreefing and modified geometric porosity

Fluid?Cstructure interaction (FSI) modeling of parachutes poses a number of computational challenges. These include the lightness of the parachute canopy compared to the air masses involved in the parachute dynamics, in the case of ringsail parachutes the geometric porosity created by the construction of the canopy from ??rings?? and ??sails?? with hundreds of ??ring gaps?? and ??sail slits,?? in the case of parachute clusters the contact between the parachutes, and ??disreefing?? from one stage to another when the parachute is used in multiple stages. The Team for Advanced Flow Simulation and Modeling (T?AFSM) has been successfully addressing these computational challenges with the Stabilized Space?CTime FSI (SSTFSI) technique, which was developed and improved over the years by the T?AFSM and serves as the core numerical technology, and a number of special techniques developed in conjunction with the SSTFSI technique. The quasi-direct and direct coupling techniques developed by the T?AFSM, which are applicable to cases with nonmatching fluid and structure meshes at the interface, yield more robust algorithms for FSI computations where the structure is light. The special technique used in dealing with the geometric complexities of the rings and sails is the homogenized modeling of geometric porosity (HMGP), which was developed and improved in recent years by the T?AFSM. The surface-edge-node contact tracking (SENCT) technique was introduced by the T?AFSM as a contact algorithm where the objective is to prevent the structural surfaces from coming closer than a minimum distance in an FSI computation. The recently-introduced conservative version of the SENCT technique is more robust and is now an essential technology in the parachute cluster computations carried out by the T?AFSM. As an additional computational challenge, the parachute canopy might, by design, have some of its panels and sails removed. In FSI computation of parachutes with such ??modified geometric porosity,?? the flow through the ??windows?? created by the removal of the panels and the wider gaps created by the removal of the sails cannot be accurately modeled with the HMGP and needs to be actually resolved during the FSI computation. In this paper we focus on parachute disreefing, including the disreefing of parachute clusters, and parachutes with modified geometric porosity, including the reefed stages of such parachutes. We describe the additional special techniques we have developed to address the challenges involved and report FSI computations for parachutes and parachute clusters with disreefing and modified geometric porosity.     

Computation of residence time in the simulation of pulsatile ventricular assist devices

A continuum-based model of particle residence time for moving-domain fluid mechanics and fluid–structure interaction (FSI) computations is proposed, analyzed, and applied to the simulation of an adult pulsatile ventricular assist device (PVAD). Residence time is a quantity of clinical interest for blood pumps because it correlates with thrombotic risk. The proposed technique may be easily implemented in any flow or FSI solver. In the context of PVADs the results of the model may be used to assess how efficiently the pump moves the blood through its interior. Three scalar measures of particle residence time are also proposed. These scalar quantities may be used in the PVAD design with the goal of reducing thrombotic risk.     

A combined fluid–structure interaction and multi–field scalar transport model for simulating mass transport in biomechanics

Mass transport processes are known to play an important role in many fields of biomechanics such as respiratory, cardiovascular, and biofilm mechanics. In this paper, we present a novel computational model considering the effect of local solid deformation and fluid flow on mass transport. As the transport processes are assumed to influence neither structure deformation nor fluid flow, a sequential one‐way coupling of a fluid–structure interaction (FSI) and a multi‐field scalar transport model is realized. In each time step, first the non‐linear monolithic FSI problem is solved to determine current local deformations and velocities. Using this information, the mass transport equations can then be formulated on the deformed fluid and solid domains. At the interface, concentrations are related depending on the interfacial permeability. First numerical examples demonstrate that the proposed approach is suitable for simulating convective and diffusive scalar transport on coupled, deformable fluid and solid domains. Copyright © 2014 John Wiley & Sons, Ltd.     

An efficient parallel computation of unsteady incompressible viscous flow with elastic moving and compliant boundaries on unstructured grids

This paper presents the development and validation of a parallel unstructured‐grid fluid–structure interaction (FSI) solver for the simulation of unsteady incompressible viscous flow with long elastic moving and compliant boundaries. The Navier–Stokes solver on unstructured moving grid using the arbitrary Lagrangian Eulerian formulation is based on the artificial compressibility approach and a high‐order characteristics‐based finite‐volume scheme. Both unsteady flow and FSI are calculated with a matrix‐free implicit dual time‐stepping scheme. A membrane model has been formulated to study fluid flow in a channel with an elastic membrane wall and their interactions. This model can be employed to calculate arbitrary wall movement and variable tension along the membrane, together with a dynamic mesh method for large deformation of the flow field. The parallelization of the fluid–structure solver is achieved using the single program multiple data programming paradigm and message passing interface for communication of data. The parallel solver is used to simulate fluid flow in a two‐dimensional channel with and without moving membrane for validation and performance evaluation purposes. The speedups and parallel efficiencies obtained by this method are excellent, using up to 16 processors on a SGI Origin 2000 parallel computer. A maximum speedup of 23.14 could be achieved on 16 processors taking advantage of an improved handling of the membrane solver. The parallel results obtained are compared with those using serial code and they are found to be identical. Copyright © 2005 John Wiley & Sons, Ltd.     

FSI modeling of the reefed stages and disreefing of the Orion spacecraft parachutes

Orion spacecraft main and drogue parachutes are used in multiple stages, starting with a “reefed” stage where a cable along the parachute skirt constrains the diameter to be less than the diameter in the subsequent stage. After a period of time during the descent, the cable is cut and the parachute “disreefs” (i.e. expands) to the next stage. Fluid–structure interaction (FSI) modeling of the reefed stages and disreefing involve computational challenges beyond those in FSI modeling of fully-open spacecraft parachutes. These additional challenges are created by the increased geometric complexities and by the rapid changes in the parachute geometry during disreefing. The computational challenges are further increased because of the added geometric porosity of the latest design of the Orion spacecraft main parachutes. The “windows” created by the removal of panels compound the geometric and flow complexity. That is because the Homogenized Modeling of Geometric Porosity, introduced to deal with the flow through the hundreds of gaps and slits involved in the construction of spacecraft parachutes, cannot accurately model the flow through the windows, which needs to be actually resolved during the FSI computation. In parachute FSI computations, the resolved geometric porosity is significantly more challenging than the modeled geometric porosity, especially in computing the reefed stages and disreefing. Orion spacecraft main and drogue parachutes will both have three stages, with computation of the Stage 1 shape and disreefing from Stage 1 to Stage 2 for the main parachute being the most challenging because of the lowest “reefing ratio” (the ratio of the reefed skirt diameter to the nominal diameter). We present the special modeling techniques and strategies we devised to address the computational challenges encountered in FSI modeling of the reefed stages and disreefing of the main and drogue parachutes. We report, for a single parachute, FSI computation of both reefed stages and both disreefing events for both the main and drogue parachutes. In the case of the main parachute, we also report, for a 2-parachute cluster, FSI computation of the disreefing from Stage 2 to Stage 3. With results from these computations, we demonstrate that we have to a great extent overcome one of the most formidable challenges in FSI modeling of spacecraft parachutes.     

A parallel iterative partitioned coupling analysis system for large-scale acoustic fluid–structure interactions

In many engineering fields, dynamic response in fluid–structure interaction (FSI) is important, and some of the FSI phenomena are treated as acoustic FSI (AFSI) problems. Dynamic interactions between fluids and structures may change dynamic characteristics of the structure and its response to external excitation parameters such as seismic loading. This paper describes a parallel coupling analysis system for large-scale AFSI problems using iterative partitioned coupling techniques. We employ an open source parallel finite element analysis system called ADVENTURE, which adopts an efficient preconditioned iterative linear algebraic solver. In addition, we have recently developed a parallel coupling tool called ADVENTURE_Coupler to efficiently handle interface variables in various parallel computing environments. We also employ the Broyden method for updating interface variables to attain robust and fast convergence of fixed-point iterations. This paper describes key features of the coupling analysis system developed, and we perform tests to validate its performance for several AFSI problems. The system runs efficiently in a parallel environment, and it is capable of analyzing three-dimensional-complex-shaped structures with more than 20 million degrees-of-freedom (DOFs). Its numerical results also show good agreement with experimental results.     

一种基于压力泊松方程的流体结构紧耦合算法

从流固耦合系统的整体控制方程出发,推导出与流体控制方程一致的耦合等价方程,并得到基于耦合方程的压力泊松方程,通过求解耦合系统压力泊松方程和一致的等价方程就能获得耦合系统的解,而不需要直接求解整体耦合系统的控制方程,有利于降低求解自由度。预估-校正多步迭代格式用于耦合系统的时间推进,克服了传统迭代耦合方法由于时间不同步而产生较大数值误差的不足。应用该方法对附带局部突起的主动脉弓动脉瘤进行流固耦合分析,验证数值方法的可行性。     

Parallel implementation of spectral element method for Lamb wave propagation modeling

The proposed spectral element method implementation is based on sparse matrix storage of local shape function derivatives calculated at Gauss–Lobatto–Legendre points. The algorithm utilizes two basic operations: multiplication of sparse matrix by vector and element‐by‐element vectors multiplication. Compute‐intensive operations are performed for a part of equation of motion derived at the degree of freedom level of 3D isoparametric spectral elements. The assembly is performed at the force vector in such a way that atomic operations are minimized. This is achieved by a new mesh coloring technique The proposed parallel implementation of spectral element method on GPU is applied for the first time for Lamb wave simulations. It has been found that computation on multicore GPU is up to 14 times faster than on single CPU. Copyright © 2015 John Wiley & Sons, Ltd.     

Fluid–structure interaction simulation of pulsatile ventricular assist devices

In this paper we present a collection of fluid–structure interaction (FSI) computational techniques that enable realistic simulation of pulsatile Ventricular Assist Devices (VADs). The simulations involve dynamic interaction of air, blood, and a thin membrane separating the two fluids. The computational challenges addressed in this work include large, buckling motions of the membrane, the need for periodic remeshing of the fluid mechanics domain, and the necessity to employ tightly coupled FSI solution strategies due to the very strong added mass effect present in the problem. FSI simulation of a pulsatile VAD at realistic operating conditions is presented for the first time. The FSI methods prove to be robust, and may be employed in the assessment of current, and the development of future, pulsatile VAD designs.     

薄膜结构气弹动力稳定性研究

将扁壳的无矩理论和流体的理想势流理论结合起来对薄膜结构的气弹动力稳定性进行了研究,提出了结构失稳的判别准则,确定了结构失稳临界风速。首先应用扁壳的无矩理论建立了薄膜结构的动力平衡方程。然后假设来流为均匀的理想势流,考虑流固耦合作用,对风向沿结构拱向和垂向时分别采用不同的气弹模型确定了作用于薄膜表面的气动力,得到了两种情况下薄膜结构的气弹动力耦合作用方程。利用Bubnov-Galerkin方法将此耦合作用方程转化为一常系数二阶微分方程,并根据Routh-Hurwitz稳定性准则确定了薄膜的失稳临界风速。最后通过对临界风速的影响因素进行分析,得到了一些重要结论,并提出了防止薄膜结构气弹失稳的一些基本措施。     

A fully-coupled fluid-structure interaction simulation of cerebral aneurysms

This paper presents a computational vascular fluid-structure interaction (FSI) methodology and its application to patient-specific aneurysm models of the middle cerebral artery bifurcation. A fully coupled fluid-structural simulation approach is reviewed, and main aspects of mesh generation in support of patient-specific vascular FSI analyses are presented. Quantities of hemodynamic interest such as wall shear stress and wall tension are studied to examine the relevance of FSI modeling as compared to the rigid arterial wall assumption. We demonstrate the importance of including the flexible wall modeling in vascular blood flow simulations by performing a comparison study that involves four patient-specific models of cerebral aneurysms varying in shape and size.