An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

Blade Wrap Angle Impact on Centrifugal Pump Performance:Entropy Generation and Fluid-Structure Interaction Analysis

The centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to analyze the influence of the blade wrap angle on flow characteristics and energy distribution of a centrifugal pump evaluated as a low specific speed with a value of 69.This study investigates six impellermodels that possess varying blade wrap angles(95°,105°,115°,125°,135°,and 145°)that were created while maintaining the same volute and other geometrical characteristics.The investigation of energy loss was conducted to evaluate the values of total and entropy generation rates(TEG,EGR).The fluid-structure interaction was considered numerically using the software tools ANSYS Fluent and ANSYSWorkbench.The elastic structural dynamic equation was used to estimate the structural response,while the shear stress transport k–ωturbulence model was utilized for the fluid domain modeling.The findings suggest that the blade wrap angle has a significant influence on the efficiency of the pump.The impeller featuring a blade wrap angle of 145°exhibits higher efficiency,with a notable increase of 3.76%relative to the original model.Variations in the blade wrap angle impact the energy loss,shaft power,and pump head.The model with a 145°angle exhibited a maximum equivalent stress of 14.8MPa and a total deformation of 0.084 mm.The results provide valuable insights into the intricate flow mechanism of the centrifugal pump,particularly when considering various blade wrap angles.

An extended multiple-support response spectrum method incorporating fluid-structure interaction for seismic analysis of deep-water bridges

The effects of ground motion spatial variability(GMSV)or fluid-structure interaction(FSI)on the seismic responses of deep-water bridges have been extensively examined.However,there are few studies on the seismic performance of bridges considering GMSV and FSI effects simultaneously.In this study,the original multiple-support response spectrum(MSRS)method is extended to consider FSI effect for seismic analysis of deep-water bridges.The solution of hydrodynamic pressure on a pier is obtained using the radiation wave theory,and the FSI-MSRS formulation is derived according to the random vibration theory.The influence of FSI effect on the related coefficients is analyzed.A five-span steel-concrete continuous beam bridge is adopted to conduct the numerical simulations.Different load conditions are designed to investigate the variation of the bridge responses when considering the GMSV and FSI effects.The results indicate that the incoherence effect and wave passage effect decrease the bridge responses with a maximum percentage of 86%,while the FSI effect increases the responses with a maximum percentage of 26%.The GMSV and FSI effects should be included in the seismic design of deep-water bridges.

Analysis and Optimization of Flow-Guided Structure Based on Fluid-Structure Interaction

Gases containing sulfur oxides can cause corrosion and failure of bellows used as furnace blowers in high-temperature environments.In order to mitigate this issue,the behavior of an effective blast furnace blower has been examined in detail.Firstly,the Sereda corrosion model has been introduced to simulate the corrosion rate of the related bellows taking into account the effects of temperature and SO_(2) gas;such results have been compared with effective measurements;then,the average gas velocity in the pipeline and the von Mises stress distribution of the inner draft tube have been analyzed using a Fluid-Structure Interaction model.Finally,the semi-closed internal corrosion environment caused by a 5 mm radial gap between the inner draft tube and the bellows has been considered.The gas flow rate in the residential space has been found to be low(0.5 ms–this value leads to a stable semi-closed internal corrosion environment for exhaust gas exchange);water phase in the exhaust gas is prone to accelerate the corrosion rate.On this basis,a bellows with an optimized inner draft tube has proposed,which includes corrosion-resistant honeycomb buffer rings.

Analysis of the Influence of the Blade Deformation on Wind Turbine Output Power in the Framework of a Bidirectional Fluid-Structure Interaction Model

The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.

Fluid-Structure Interaction in Problems of Patient Specific Transcatheter Aortic Valve Implantation with and Without Paravalvular Leakage Complication

Paravalvular Leakage(PVL)has been recognized as one of the most dangerous complications in relation to Transcathether Aortic Valve Implantation(TAVI)activities.However,data available in the literature about Fluid Structure Interaction(FSI)for this specific problem are relatively limited.In the present study,the fluid and structure responses of the hemodynamics along the patient aorta model and the aortic wall deformation are studied with the aid of numerical simulation taking into account PVL and 100%TAVI valve opening.In particular,the aorta without valve(AWoV)is assumed as the normal condition,whereas an aorta with TAVI 26 mm for 100%Geometrical Orifice Area(GOA)is considered as the patient aorta with PVL complication.A 3D patient-specific aorta model is elaborated using the MIMICS software.Implantation of the identical TAVI valve of Edward SAPIEN XT 26(Edwards Lifes ciences,Irvine,California)is considered.An undersized 26 mm TAVI valve with 100%valve opening is selected to mimic the presence of PVL at the aortic annulus.The present research indicates that the existence of PVL can increase the blood velocity,pressure drop and WSS in comparison to normal conditions,thereby paving the way to the development of recirculation flow,thrombus formation,aorta wall collapse,aortic rupture and damage of endothelium.

Improved frequency modeling and solution for parallel liquid-filled pipes considering both fluid-structure interaction and structural coupling

The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.

Fluid-structure interaction simulation of three-dimensional flexible hydrofoil in water tunnel

The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.

Effects of renal artery stenosis on realistic model of abdominalaorta and renal arteries incorporating fluid-structureinteraction and pulsatile non-Newtonian blood flow

The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.

An improved algorithm for fluid-structure interaction of high-speed trains under crosswind

Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.

Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.

Fluid-Structure Interaction Simulation of Aqueous Outflow System in Response to Juxtacanalicular Meshwork Permeability Changes with a Two-Way Coupled Method

Elevated intraocular pressure appears to have a broader impact on increased resistance to aqueous humor outflow through the conventional aqueous outflow system(AOS).However,there is still no consensus about exact location of the increased outflow resistance of aqueous humor,and the mechanism is not perfect.In addition,it is difficult to accurately obtain hydrodynamic parameters of aqueous humor within the trabecular meshwork outflow pathways based on the current technology.In this paper,a two-way fluid-structure interaction simulation was performed to study the pressure difference and velocity in the superficial trabecular meshwork,juxtacanalicular meshwork(JCM)and Schlemm’s canal in response to JCM permeability changes.We obtained the JCM permeability of normal intraocular pressure varied between 1×10?15 m2 and 10×10?15 m2 while permeability of the JCM ranged from 2×10?16 m2 and 3×10?16 m2 under conditions of high intraocular pressure.The study indicated that the fluid dynamics parameters in trabecular meshwork and Schlemm’s canal are most significantly affected by the changes of JCM permeability.Moreover,the study demonstrates that the finite element modeling of AOS provides a practical means for studying the outflow dynamics and the biomechanical environment of the AOS.

Dynamic Adaptive Finite Element Analysis of Acoustic Wave Propagation Due to Underwater Explosion for Fluid-structure Interaction Problems

In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.

OpenIFEM:A High Performance Modular Open-Source Software of the Immersed Finite Element Method for Fluid-Structure Interactions

We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.

3D numerical simulation on fluid-structure interaction of structure subjected to underwater explosion with cavitation

In the underwater-shock environment, cavitation occurs near the structural surface. The dynamic response of fluid-structure interactions is influenced seriously by the cavitation effects. It is also the difficulty in the field of underwater explosion. With the traditional boundary element method and the finite element method （FEM）, it is difficult to solve the nonlinear problem with cavitation effects subjected to the underwater explosion. To solve this problem, under the consideration of the cavitation effects and fluid compressibility, with fluid viscidity being neglected, a 3D numerical model of transient nonlinear fluid-structure interaction subjected to the underwater explosion is built. The fluid spectral element method （SEM） and the FEM are adopted to solve this model. After comparison with the FEM, it is shown that the SEM is more precise than the FEM, and the SEM results are in good coincidence with benchmark results and experiment results. Based on this, combined with ABAQUS, the transient fluid-structure interaction mechanism of the 3D submerged spherical shell and ship stiffened plates subjected to the underwater explosion is discussed, and the cavitation region and its influence on the structural dynamic responses are presented. The paper aims at providing references for relevant research on transient fluid-structure interaction of ship structures subjected to the underwater explosion.

Study of velocity effects on parachute inflation performance based on fluid-structure interaction method

The inflation of a five-ring cone parachute with the airflow velocity of 18 m/s is studied based on the simplified arbitrary Lagrange Euler （SALE）/fluid-structure interaction （FSI） method. The numerical results of the canopy shape, stability, opening load, and drag area are obtained, and they are well consistent with the experimental data gained from wind tunnel tests. The method is then used to simulate the opening process under different velocities. It is found that the first load shock affected by the velocity often occurs at the end of the initial inflation stage. For the first time, the phenomena that the inflation distance proportion coefficient increases and the dynamic load coefficient decreases, respectively, with the increase in the velocity are revealed. The above proposed method is competent to solve the large deformation problem without empirial coefficients, and can collect more space-time details of fluid-structure-motion information when it is compared with the traditional method.

Patient-Specific Echo-Based Fluid-Structure Interaction Modeling Study of Blood Flow in the Left Ventricle with Infarction and Hypertension

Understanding cardiac blood flow behaviors is of importance for cardiovascular research and clinical assessment of ventricle functions.Patient-specific Echo-based left ventricle(LV)fluid-structure interaction(FSI)models were introduced to perform ventricle mechanical analysis,investigate flow behaviors,and evaluate the impact of myocardial infarction(MI)and hypertension on blood flow in the LV.Echo image data were acquired from 3 patients with consent obtained:one healthy volunteer(P1),one hypertension patient(P2),and one patient who had an inferior and posterior myocardial infarction(P3).The nonlinear Mooney-Rivlin model was used for ventricle tissue with material parameter values chosen to match echo-measure LV volume data.Using the healthy case as baseline,LV with MI had lower peak flow velocity(30%lower at beginejection)and hypertension LV had higher peak flow velocity(16%higher at begin-filling).The vortex area(defined as the area with vorticity>0)for P3 was 19%smaller than that of P1.The vortex area for P2 was 12%smaller than that of P1.At peak of filling,the maximum flow shear stress(FSS)for P2 and P3 were 390%higher and 63%lower than that of P1,respectively.Meanwhile,LV stress and strain of P2 were 41%and 15%higher than those of P1,respectively.LV stress and strain of P3 were 36%and 42%lower than those of P1,respectively.In conclusion,FSI models could provide both flow and structural stress/strain information which would serve as the base for further cardiovascular investigations related to disease initiation,progression,and treatment strategy selections.Large-scale studies are needed to validate our findings.

Failure pressure calculation of fracturing well based on fluid-structure interaction

Failure pressure is a key parameter in reservoir hydrofracturing operation. Existing analytical methods for calculating the failure pressure are based on the assumption that borehole fluid is under two extreme conditions: non-infiltration or complete infiltration. The assumption is not suitable for the actual infiltration process, and this will cause a great error in practical calculation. It shows that during the injection process, the dynamic variation in effective stress-dependent permeability has an influence on the infiltration, and the influence also brings about calculation errors. Based on the fluid-structure interaction and finite element method (FEM), considering partial infiltration during injection process, a numerical model for calculating rock failure pressure is established. According to the analysis of permeability test results and response-surface method, a new variation rule of rock permeability with the change of effective stress is presented, and the relationships among the permeability, confining pressure and pore pressure are proposed. There are some differences between the dynamic value of permeability-effective-stress coefficient observed herein and the one obtained by the classical theory. Combining with the numerical model and the dynamic permeability, a coupling method for calculating failure pressure is developed. Comparison of field data and calculated values obtained by various methods shows that accurate values can be obtained by the coupling method. The coupling method can be widely applied to the calculation of failure pressure of reservoirs and complex wells to achieve effective fracturing operation.