Including Fluid Structure Interaction Research Articles

Ten Years of Passion: I.S. Gromeka’s Contribution to Science

The work and life of Ippolit Stepanovich Gromeka is reviewed. Gromeka authored a classical set of eleven papers on fluid dynamics in just ten years before a tragic illness ended his life. Sadly, he is not well known to the western scientific community because all his publications were written in Russian. He is one of the three authors who independently derived an analytical solution for accelerating laminar pipe flow. He was the first to eliminate the contradiction between the theories of Young and Laplace on capillary phenomena. He initiated the theoretical basis of helical (Beltrami) flow, and he studied the movement of cyclones and anticyclones seventeen years before Zermelo (whose work is considered as pioneering). He is also the first to analyse wave propagation in liquid-filled hoses, thereby including fluid–structure interaction.

Accurate Calculation of FFR Based on a Physics-Driven Fluid‐Structure Interaction Model

Background: The conventional FFRct numerical calculation method uses a model with a multi-scale geometry based upon CFD, and rigid walls. Therefore, important interactions between the elastic vessel wall and blood flow are not routinely considered. Changes in the resistance of coronary microcirculation during hyperaemia are likewise not typically incorporated using a fluid–structure interaction (FSI) algorithm. It is likely that both have resulted in FFRct calculation errors. Objective: In this study we incorporated both the influence of vascular elasticity and coronary microcirculatory structure on FFR, to improve the accuracy of FFRct calculation. Thus, in this study, a physics-driven 3D–0D coupled model including fluid–structure interaction was established to calculate accurate FFRct values. Methods: Based upon a novel geometric multi-scale modeling technology, a FSI simulation approach was used. A lumped parameter model (0D) was used as the outlet boundary condition for the 3D FSI coronary artery model to incorporate physiological microcirculation, with bidirectional coupling between the two models. Results: The accuracy, sensitivity, specificity, and both positive and negative predictive values of FFRDC calculated based upon the coupled 3D–0D model were 86.7, 66.7, 84.6, 66.7, and 91.7%, respectively. Compared to the calculated value using the basic CFD model (MSE = 5.9%, accuracy rate = 80%), the FFRCFD calculated based on the coupled 3D–0D model has a smaller MSE of 1.9%. Conclusion: The physics-driven coupled 3D–0D model that incorporates fluid–structure interactions not only consider the influence of the elastic vessel wall on blood flow, but also provides reliable microvascular resistance boundary conditions for the 3D FSI model. This allows for a calculation that is based upon conditions that are closer to the physiological environment, and thus improves the accuracy of FFRct calculation. It is likely that more accurate information will provide an enhanced recommendation regarding percutaneous coronary intervention (PCI) in the clinic.

A fluid–structure interaction model for assessing the safety of flood gate vibrations due to wave impacts

This paper establishes a computationally efficient model to predict flood gate vibrations due to wave impacts including fluid–structure interaction. In contrast to earlier models, composite fluid domains are included to represent the situation of a flood gate in a dewatering sluice with the presence of an overhang that causes the confined-wave impacts. The dynamic response of the gate-fluid system is derived in the frequency domain using a substructuring mode matching technique, in which the gate vibrations are first expressed in terms of in-vacuo modes while the liquid motion is described as a superposition of linear potentials. Pressure impulse theory is employed to predict the impulsive wave impact loads, which are superposed on the quasi-steady wave loads. The computational efficiency of the developed model allows for a large number of simulations. This makes it possible for the first time to perform probabilistic evaluations for this type of problems without doing concessions on the accuracy of the physical modelling of the involved fluid–structure interaction processes. This is demonstrated by application of the developed models within a probabilistic framework resulting in the explicit quantification of the failure probability of flood gates subjected to wave impacts.

A linear-elasticity-based mesh moving method with no cycle-to-cycle accumulated distortion

Good mesh moving methods are always part of what makes moving-mesh methods good in computation of flow problems with moving boundaries and interfaces, including fluid–structure interaction. Moving-mesh methods, such as the space–time (ST) and arbitrary Lagrangian–Eulerian (ALE) methods, enable mesh-resolution control near solid surfaces and thus high-resolution representation of the boundary layers. Mesh moving based on linear elasticity and mesh-Jacobian-based stiffening (MJBS) has been in use with the ST and ALE methods since 1992. In the MJBS, the objective is to stiffen the smaller elements, which are typically placed near solid surfaces, more than the larger ones, and this is accomplished by altering the way we account for the Jacobian of the transformation from the element domain to the physical domain. In computing the mesh motion between time levels t_n and t_{n+1} with the linear-elasticity equations, the most common option is to compute the displacement from the configuration at t_n. While this option works well for most problems, because the method is path-dependent, it involves cycle-to-cycle accumulated mesh distortion. The back-cycle-based mesh moving (BCBMM) method, introduced recently with two versions, can remedy that. In the BCBMM, there is no cycle-to-cycle accumulated distortion. In this article, for the first time, we present mesh moving test computations with the BCBMM. We also introduce a version we call “half-cycle-based mesh moving” (HCBMM) method, and that is for computations where the boundary or interface motion in the second half of the cycle consists of just reversing the steps in the first half and we want the mesh to behave the same way. We present detailed 2D and 3D test computations with finite element meshes, using as the test case the mesh motion associated with wing pitching. The computations show that all versions of the BCBMM perform well, with no cycle-to-cycle accumulated distortion, and with the HCBMM, as the wing in the second half of the cycle just reverses its motion steps in the first half, the mesh behaves the same way.

Blast Response of Single-Degree-of-Freedom Systems Including Fluid-Structure Interaction

AbstractThe state of the practice for the blast-resistant design of structures uses a simplified single-degree-of-freedom (SDOF) analysis and estimates the blast loads assuming reflection from an u...

Application of a Krylov subspace method for an efficient solution of acoustic transfer functions

Solving acoustic radiation problems, arising from systems including fluid–structure interaction, is of interest in many engineering applications. Computing frequency response functions over a large frequency range is a concern in such applications. A method which solves the Helmholtz equation for multiple frequencies in one step is the matrix-Padé-via-Lanczos connection for unsymmetric systems, as presented by Wagner et al. [1]. The present work is based on Ref. [1] and presents a method for efficiently computing frequency responses over a frequency range for coupled structural-acoustic problems, where the structure and the acoustic near field are discretized with finite elements and an analytical Dirichlet-to-Neumann map approximates the far field. The method is based on a Krylov-subspace projection technique which derives a matrix-valued Padé approximation for a restricted area in the near field and the pressure field on a spherical boundary. On the spherical boundary, where the finite domain is truncated, the non-local modified Dirichlet-to-Neumann operator is applied as a low-rank update matrix. The present contribution extends this method and incorporates new techniques for a more stable model reduction through the Lanczos algorithm and a novel weighted adaptive windowing technique. Further, structural damping is incorporated, for computing the acoustic radiation of a harmonically excited plate. These computed results are compared with acoustic measurements in an anechoic chamber and verified with computational results obtained with a commercial code that uses the perfectly matched layer method.

Sloshing Impact Response Analysis for Insulation System of LNG CCS Considering Elastic Support Effects of Hull Structures

The sloshing pressure acting on a membrane-type LNG CCS is a typical irregular impact load, and the structural response of a tank system induced by sloshing also shows very complex behavior, including fluid structure interaction. Therefore, it is not easy to accurately estimate the sloshing impact pressures and resulting structural response. Moreover, a huge time consuming process to deal with the enormous pressure data obtained during a model tank test and the following structural analysis would be inevitable. To reduce the computation time for structural analysis, in this study, a rational structural modeling strategy was considered, and a simplified scheme to analyze the dynamic structural responses of an LNG CCS was introduced, which was based on the concept of the linear combination of the triangular response functions obtained by a transient response analysis of structures under unit triangular impact pressure. A structural analysis of a real Mark III membrane type insulation system under the sloshing impact pressure time histories obtained by model tests was performed using the various proposed structural models and simplified analysis scheme. The results were investigated in detail, including the elastic support effects of the hull structure.

Assessment of URANS surface effect ship models for calm water and head waves

Surface effect ship (SES) air cushion and seal models are implemented in an URANS hydrodynamics solver. The air cushion is modeled either as a prescribed pressure patch, or as a compressible isothermal/adiabatic ideal stagnant air with fan and leakage flows. The seals are either discretized as hinged bodies or modeled as 2D planing surfaces with hydrodynamic interaction. Verification and validation studies are performed using T-Craft experimental data for calm water resistance, sinkage and trim at Froude number (Fr)=0.1–0.6; impulsive heave and pitch decay at Fr=0; and wave-induced resistance and motion predictions in head waves at Fr=0 and 0.6. The compressible air cushion model with fan and leakage flows perform better than those without the fan and leakage flows and the prescribed pressure patch model. The hinged seal model performs better than the 2D planing surface model, but is computationally expensive for time accurate simulations. Therefore, the 2D planing surface model is used for the validation studies. SES simulations on grids with 5.3M cells show grid verification intervals of 6%, which are comparable to those reported for displacement and semi-planing hull studies on similar grid sizes. On an average calm water and impulsive motion predictions compare within 8.5% of the experimental data, and wave-induced motion predictions show somewhat larger error of 13.5%. The errors levels are mostly comparable to those for displacement and semi-planing and planing hulls. The study identifies that most critical advancement needed for SES simulations is the seal modeling including fluid structure interaction.

Seismic Design of a Double Deck Floating Roof Type Used for Liquid Storage Tanks

Sloshing response of a cylindrical liquid storage tank with the double deck type floating roof (DDFR) subjected to seismic excitation is considered in this paper. The aim of the paper is to clarify the significant parameters that should be considered in the seismic design of a DDFR and proposing a practical seismic design procedure for evaluating the dynamic stresses inside a DDFR. A numerical method including fluid–structure interaction and the geometry details of a DDFR tank are established. The geometric nonlinear effects on the seismic behavior of the DDFR as well as the accuracy of common analytical solution suggested in the literature are examined by the numerical model. The numerical results show that the geometric nonlinear effects can considerably reduce the seismic stress in DDFR, but have no significant effect on the liquid hydrodynamic pressure exerted on the DDFR and the roof's vertical displacement. It is also revealed that not only the general displacement of DDFR but also the local effects of liquid hydrodynamic pressure on the bottom plate should be considered for seismic design of a DDFR. Finally, a design procedure for the evaluation of dynamic stress in the DDFR due to the seismic loads is proposed and discussed.

Study of composite plate in water with transient and steady state motions

Hydrodynamic loading and dynamic response of a composite plate were measured inside a tow tank while the speed of the plate increased from zero to various steady-state values and remained at the steady state values. Strain responses of the plate were recorded using strain gages. To determine the effect of the composite plate with fluid–structure interaction, a stiffer and heavier aluminum plate was compared to the composite plate. The composite plate was tested under various setting angles relative to the towing direction. Some experimental data were compared to numerical results at steady state motions, and they agreed well each other. The study showed that the hydrodynamic loading was greater on the polymer composite plate than on the aluminum plate of the same size. Furthermore, the hydrodynamic force was much greater during the transient period of the increasing speed than at the final steady state speed. The ratio of the peak strain in the composite plate during the transient period to the strain at the steady state motion depended on the plate speed and the setting angles significantly. Therefore, design and analysis of composite structures moving in water should consider the hydrodynamic loading during the transient motion including fluid–structure interaction.

ST and ALE-VMS methods for patient-specific cardiovascular fluid mechanics modeling

This paper provides a review of the space–time (ST) and Arbitrary Lagrangian–Eulerian (ALE) techniques developed by the first three authors' research teams for patient-specific cardiovascular fluid mechanics modeling, including fluid–structure interaction (FSI). The core methods are the ALE-based variational multiscale (ALE-VMS) method, the Deforming-Spatial-Domain/Stabilized ST formulation, and the stabilized ST FSI technique. A good number of special techniques targeting cardiovascular fluid mechanics have been developed to be used with the core methods. These include: (i) arterial-surface extraction and boundary condition techniques, (ii) techniques for using variable arterial wall thickness, (iii) methods for calculating an estimated zero-pressure arterial geometry, (iv) techniques for prestressing of the blood vessel wall, (v) mesh generation techniques for building layers of refined fluid mechanics mesh near the arterial walls, (vi) a special mapping technique for specifying the velocity profile at an inflow boundary with non-circular shape, (vii) a scaling technique for specifying a more realistic volumetric flow rate, (viii) techniques for the projection of fluid–structure interface stresses, (ix) a recipe for pre-FSI computations that improve the convergence of the FSI computations, (x) the Sequentially-Coupled Arterial FSI technique and its multiscale versions, (xi) techniques for calculation of the wall shear stress (WSS) and oscillatory shear index (OSI), (xii) methods for stent modeling and mesh generation, (xiii) methods for calculation of the particle residence time, and (xiv) methods for an estimated element-based zero-stress state for the artery. Here we provide an overview of the special techniques for WSS and OSI calculations, stent modeling and mesh generation, and calculation of the residence time with application to pulsatile ventricular assist device (PVAD). We provide references for some of the other special techniques. With results from earlier computations, we show how these core and special techniques work.

Effects of fluid structure interaction in a three dimensional model of the spinal subarachnoid space

It is unknown whether spinal cord motion has a significant effect on cerebrospinal fluid (CSF) pressure and therefore the importance of including fluid structure interaction (FSI) in computational fluid dynamics models (CFD) of the spinal subarachnoid space (SAS) is unclear. This study aims to determine the effects of FSI on CSF pressure and spinal cord motion in a normal and in a stenosis model of the SAS. A three-dimensional patient specific model of the SAS and spinal cord were constructed from MR anatomical images and CSF flow rate measurements obtained from a healthy human being. The area of SAS at spinal level T4 was constricted by 20% to represent the stenosis model. FSI simulations in both models were performed by running ANSYS CFX and ANSYS Mechanical in tandem. Results from this study show that the effect of FSI on CSF pressure is only about 1% in both the normal and stenosis models and therefore show that FSI has a negligible effect on CSF pressure.

Computational fluid dynamics as a modern tool for engineering characterization of bioreactors

Since design, construction and evaluation of bioreactors for large-scale production is costly and time consuming, computational methods may give some insights into the fluid mechanics within bioreactors. Thus, critical limiting factors, such as insufficient mixing as well as inhomogeneous nutrient and oxygen mass transfer, may be identified early in the design process. Although advanced experimental techniques such as laser Doppler anemometry and particle image velocimetry are also reliable, they are too time consuming to characterize the complete flow pattern in industrial scales and rely on optical accessibility. Therefore, the knowledge of flow characteristics provided by computational fluid dynamics (CFD) is indispensable for the rational design of bioreactors. Based on previously published reviews, the present work summarizes the latest publications on the usage of CFD to characterize and scale-up bioreactors used in biotechnological processes. Selected models that are used to predict the fluid flow pattern and key engineering parameters of commonly used bioreactors are described. Related issues, such as grid dependency of CFD results and the requirement for experimental verification are also addressed. Finally, an overview of proposed but not yet feasible CFD applications is presented, including fluid–structure interaction, the use of direct numerical simulation and the coupling fluid flow and chemical reactions.

Seismic analysis of rectangular water-containing structures with floating ice blocks

This paper presents a new formulation to investigate the effects of floating ice blocks on seismically-excited rectangular water-containing structures. The proposed method is based on a sub-structuring approach, where the flexible containing structure and ice-added mass are modeled using finite elements, while hydrodynamic effects are modeled analytically through interaction forces at the water–structure and water–ice interfaces, thus eliminating the need for reservoir finite element discretization. In addition to accounting for the influence of floating ice blocks and container walls' flexibility, the developed frequency- and time-domain techniques also include the effects of container geometrical or material asymmetry as well as the coupling between convective and impulsive components of hydrodynamic pressure. The proposed formulation is illustrated through a numerical example illustrating the dynamic response of symmetric and asymmetric water-containing structures covered with floating ice blocks. Obtained time- and frequency-domain responses are successfully validated against advanced finite element analyses including fluid–structure interaction capabilities. For the water-containing structures studied, the results show that the presence of floating ice blocks affects the frequency content and amplitudes of the dynamic responses corresponding to convective and impulsive modes. ► New formulation for seismically-excited ice-covered water-containing structures ► Reservoir modeled analytically thus no need for its finite element discretization ► Effects of floating ice blocks, asymmetry, convective and impulsive modes ► Obtained dynamic responses successfully validated against finite element solutions ► Example highlighting the effects of floating ice blocks on dynamic response

Three-dimensional modal analysis of an idealized human head including fluid–structure interaction effects

A three-dimensional analytical and numerical method is presented in this article for the analysis of the acoustic fluid–structure interaction systems including, but not limited to, the brain, cerebro-spinal fluid (CSF), and skull. The model considers a three-dimensional acoustic fluid medium interacting with two solid domains. This article deals with the analytical and numerical computation of eigenproperties for an idealized human head model including fluid–structure interaction phenomena. We determine in the present work the natural frequencies and the modes shapes of the system of the brain, cerebro-spinal fluid (CSF), and skull. Two models are presented in this study: an elastic skull model and a rigid model. In the analysis, a potential technique is used to obtain in three-dimensional cylindrical coordinates a general solution for a solid problem. A finite element method analysis is also used to check the validity of the present method. The results from the proposed method are in good agreement with numerical solutions. The effects of the fluid thickness and compressibility on the natural frequencies are also investigated.

Practical dynamic analysis of structures laterally vibrating in contact with water

This paper proposes a practical formulation to investigate the dynamic response of structures laterally vibrating in contact with water on one or both sides. The proposed technique accounts for structure’s flexibility, soil flexibility, varying water levels and various boundary conditions. Simplified procedures are developed for practical assessment of the vibration periods, hydrodynamic loads and seismic response of structure-water systems including higher mode effects. The efficiency of the proposed techniques is validated through examples of structure–water systems with different configurations. We show that the proposed methods give excellent results when compared to more advanced finite element solutions including fluid–structure interaction capabilities.

Optimal transportation meshfree approximation schemes for fluid and plastic flows

AbstractWe develop an optimal transportation meshfree (OTM) method for simulating general solid and fluid flows, including fluid–structure interaction. The method combines concepts from optimal transportation theory with material‐point sampling and max‐ent meshfree interpolation. The proposed OTM method generalizes the Benamou–Brenier differential formulation of optimal mass transportation problems to problems including arbitrary geometries and constitutive behavior. The OTM method enforces mass transport and essential boundary conditions exactly and is free from tension instabilities. The OTM method exactly conserves linear and angular momentum and its convergence characteristics are verified in standard benchmark problems. We illustrate the range and scope of the method by means of two examples of application: the bouncing of a gas‐filled balloon off a rigid wall; and the classical Taylor‐anvil benchmark test extended to the hypervelocity range. Copyright © 2010 John Wiley & Sons, Ltd.

Shape optimal design of arch dams using an adaptive neuro-fuzzy inference system and improved particle swarm optimization

An efficient methodology is proposed to find the optimal shape of arch dams including fluid–structure interaction subject to earthquake ground motion. In order to reduce the computational cost of optimization process, an adaptive neuro-fuzzy inference system (ANFIS) is built to predict the dam effective response instead of directly evaluating it by a time-consuming finite element analysis (FEA). The presented ANFIS is compared with a widespread neural network termed back propagation neural network (BPNN) and it appears a better performance generality for estimating the dam response. The optimization task is implemented using an improved version of particle swarm optimization (PSO) named here as IPSO. In order to assess the effectiveness of the proposed methodology, the optimization of a real world arch dam is performed via both IPSO–ANFIS and PSO–BPNN approaches. The numerical results demonstrate the computational advantages of the proposed IPSO–ANFIS for optimal design of arch dams when compared with the PSO–BPNN approach.

An overview of models for viscothermal wave propagation, including fluid structure interaction

In acoustics, the standard wave propagation models neglect the effects of viscosity and thermal conductivity. When waves propagate in narrow tubes or thin layers these simplifications might not be accurate. This paper presents an overview of models that take the effects of inertia, viscosity, thermal conductivity and compressibility into account. Based on the use of dimensionless parameters, three classes of models are outlined. The most important dimensionless parameter is the shear wave number, an unsteady Reynolds number that indicates the ratio between inertial and viscous effects. These viscothermal wave propagation models can be coupled to structural models to capture the fluid structure interaction. Analytical solutions can be found for these coupled analysis cases only for simple geometries and boundary conditions. For more complex geometries, numerical models were developed. Examples of applications of these models are also presented.

Elasto-plastic earthquake response of arch dams including fluid–structure interaction by the Lagrangian approach

Elasto-plastic earthquake response of arch dams including fluid–structure interaction by the Lagrangian approach is mainly investigated in this study. To this aim, three-dimensional eight-noded version of Lagrangian fluid finite element including the effects of compressible wave propagation and surface sloshing motion, and three-dimensional version of Drucker–Prager model based on associated flow rule assumption were programmed in FORTRAN language by authors and incorporated into the program NONSAP. Two new components added into the NONSAP were tested on a simple fluid tank and a simple fluid–structure system and obtained very reasonable results. Elasto-plastic earthquake analysis of an arch dam including fluid–structure interaction by the Lagrangian approach was performed using Wilson- θ method. The El Centro N–S component of Imperial Valley earthquake in 1940 was chosen as a ground motion. The component considered was applied in upstream–downstream direction. Displacements at various nodes on dam body, hydrodynamic pressures at three fluid elements near the dam–reservoir interface and yield functions at two elements predicted cracks in dam body were obtained. The elasto-plastic response of the arch dam was also compared to its linear response. It is seen that Drucker–Prager model in the representing elasto-plastic behaviour of arch dams and in estimating location of cracks in dam body is a fast tool.