Study on fluid-structure interaction in liquid oxygen feeding pipe systems using finite volume method

The fluid-structure interaction may occur in space launch vehicles,which would lead to bad performance of vehicles,damage equipments on vehicles,or even affect astronauts＇ health.In this paper,analysis on dynamic behavior of liquid oxygen （LOX） feeding pipe system in a large scale launch vehicle is performed,with the effect of fluid-structure interaction （FSI） taken into consideration.The pipe system is simplified as a planar FSI model with Poisson coupling and junction coupling.Numerical tests on pipes between the tank and the pump are solved by the finite volume method.Results show that restrictions weaken the interaction between axial and lateral vibrations.The reasonable results regarding frequencies and modes indicate that the FSI affects substantially the dynamic analysis,and thus highlight the usefulness of the proposed model.This study would provide a reference to the pipe test,as well as facilitate further studies on oscillation suppression.

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.

ALE Fractional Step Finite Element Method for Fluid-Structure Nonlinear Interaction Problem

A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.

Methodology for Comparing Coupling Algorithms for Fluid-Structure Interaction Problems

The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.

Three-Dimensional Fluid-Structure Interaction Case Study on Cubical Fluid Cavity with Flexible Bottom

In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lombardi(2012), and Trimarchi(2012). We focus on a 3D test case in which we evaluated the sensitivity of several input parameters on the fluid and structural results. In particular, this analysis provides a starting point from which we can look deeper into specific aspects of these simulations and analyze more realistic cases, e.g., in sails design. In this study, using the commercial software ADINATM, we addressed a well-known unsteadiness problem comprising a square box representing the fluid domain with a flexible bottom modeled with structural shell elements. We compared data from previously published work whose authors used the same numerical approach, i.e., a partitioned approach coupling a finite volume solver(for the fluid domain) and a finite element solver(for the solid domain). Specifically, we established several benchmarks and made comparisons with respect to fluid and solid meshes, structural element types, and structural damping, as well as solution algorithms. Moreover, we compared our method with a monolithic finite element solution method. Our comparisons of new and old results provide an outline of best practices for such simulations.

DYNAMIC CHARACTERISTIC ANALYSIS OF A 3-D SEMI-SUBMERGED BODY AS A FLUID-STRUCTURE INTERACTION SYSTEM

An Arnoldi's method with new iteration pattern,which was designed for solving a large unsymmetric eigenvalue problem introduced by displacement-pressure FE (Finite Element) pattern of a fluid-structure interaction system,was adopted here to get the dynamic characteristics of the semi-submerged body. The new iteration pattern could be used efficiently to obtain the Arnoldi's vectors in the shift-frequency technique,which was used for the zero-frequency problem. Numerical example showed that the fluid-structure interaction is one of the important factors to the dynamic characteristics of large semi-submerged thin-walled structures.

Modeling and Simulation of Valve Cycle in Vein Using an Immersed Finite Element Method

A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.

Conjugate Gradient Method to Solve Fluid Structure Interaction Problem

In this paper, we propose a method to solve coupled problem. Our computational method is mainly based on conjugate gradient algorithm. We use finite difference method for the structure and finite element method for the fluid. Conjugate gradient method gives suitable numerical results according to some papers.

Numerical Simulation of the Solitary Wave Interacting with an Elastic Structure Using MPS-FEM Coupled Method

Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.

Multi-resolution technique integrated with smoothed particle element method (SPEM) for modeling fluid-structure interaction problems with free surfaces

Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.

An Efficient Parallel/Unstructured-Multigrid Implicit Method for Simulating 3D Fluid-Structure Interaction

A finite volume(FV)method for simulating 3D Fluid-Structure Interaction(FSI)is presented in this paper.The fluid flow is simulated using a parallel unstructured multigrid preconditioned implicit compressible solver,whist a 3D matrix-free implicit unstructured multigrid finite volume solver is employed for the structural dynamics.The two modules are then coupled using a so-called immersed membrane method(IMM).Large-Eddy Simulation(LES)is employed to predict turbulence.Results from several moving boundary and FSI problems are presented to validate proposed methods and demonstrate their efficiency。

液态金属堆内小长径比同轴双层薄壁结构的流固耦合试验研究

池式液态金属堆的主容器与热屏之间具有狭窄的流体间隙,该结构属于含窄缝间隙的小长径比同轴双层柔性壳体。窄缝间隙中流体与结构之间相互作用产生的流固耦合效应在抗震分析时必须加以考虑。现有研究主要针对长径比较大和流体间隙相对较大的圆柱壳体,对于小长径比和极小流体间隙的同轴柔性壳体研究较为缺乏。本文参考液态金属反应堆的结构设计了三种不同尺寸的窄缝间隙,开展窄缝间隙条件小长径比同轴双层壳体模型的振动试验,测量结构在不同间隙尺寸下的加速度,得到模态频率和主要振型。然后,使用有限元法进行模态分析,通过与试验结果的对比,验证了有限元法的准确性。最后,计算出不同间隙尺寸下模型的附加质量;且随着间隙尺寸的减小,模型附加质量随之增大。本研究可为类似的含窄缝小长径比同轴双层柔性壳体结构的抗震设计提供数据支持,对于液态金属堆的抗震分析具有重要意义。

VIBRATION CHARACTERISTICS OF FLUID-STRUCTURE INTERACTION OF CONICAL SPIRAL TUBE BUNDLE

Heat transfer enhancement is achieved by flow-induced vibration in elastic tube bundles heat exchangers. For a further understanding of heat transfer enhancement mechanism and tube structure optimization, it is of importance to study the vibration characteristics of fluid-structure interaction of tube bundles. The finite element method is applied in the study of fluid-structure interaction of a new type elastic heat transfer element, i.e., the dimensional conical spiral tube bundle. The vibration equation and element matrix for the tube are set up by the regulation of different helical angles and coordinate transformation, together with the simplification of the joint body of the two pipes. The vibration characteristics of conical spiral tube bundle are analyzed at different velocities of the tube-side flow, and the critical velocity of vibration buckling is obtained. The results show that the natural frequency of the tube bundle decreases as the flow speed increases, especially for the first order frequency, and the critical velocity of vibration buckling is between 1.2665 m/s-1.2669 m/s. The vibration mode of conical spiral tube bundle is mainly z-axial, which is feasible to be induced and controlled.

Fluid-structure interaction in Z-shaped pipe with different supports

Fluid-structure interaction(FSI)has a strong relation with layout of fluid delivery system.FSI is liable to cause local damage.Thus,FSI analysis is necessary in many cases,especially for flexible pipe systems.FSI modeling consists of eight governing equations and then completely solved via the finite volume method(FVM).Friction,Poisson and joint couplings were discussed in detail to reveal the influence of a Z-shaped pipe with different supports and elbows on FSI.After the feasibility of solving FSI by FVM was verified,the different effects of free,fixed and elastic supports on FSI in the commonly used and simplified Z-shaped pipe were further analyzed.Results indicated that different support stiffness lead to various FSI responses.If coupling occurs at the elbow and less support is considered,then the pipe has a relatively large amplitude and complex pressure fluctuation.

DYNAMIC ANALYSIS OF FLUID-STRUCTURE INTERACTION OF ENDOLYMPH AND CUPULA IN THE LATERAL SEMICIRCULAR CANAL OF INNER EAR

The semicircular canals, composed of lateral, anterior and posterior canals in the inner ear, are the sensors of equilibrium during head rotation movements in the three-dimensional space. Semicircular canals are filled with endolymph confined by the cupula. The study of the relationship between endolymph flow and cupular deformation is important in revealing the semicircular canals biomechanical behavior. To date, there are few studies focusing on the transient endolymph flow and cupular deformation in response to a head rotation motion. The lateral semicircular canal is mainly responsible for the sense of the horizontal rotation movement. In order to figure out the intricate dynamics in the lateral semicircular canal during the head rotation motion, the time evolutions of both endolymph flow and cupular deformation are analyzed in this article by using a fully coupled fluid-structure interaction model. It is shown that the cupular deformation provides cues for understanding the physiology of sensing the head rotation.

FLUID-STRUCTURE INTERACTION OF HYDRODYNAMIC DAMPER DURING THE RUSH INTO THE WATER CHANNEL

The hydrodynamic damper is a device to decrease the motion of armament carrier by use of the water resistance. When hydrodynamic damper rushes into the water channel with high velocity, it is a complicated flow phenomenon with fluid-structure interaction, free surface and moving interface. Numerical simulation using the Smoothed Particle Hydrodynamics （SPH） method coupled with the Finite Element （FE） method was successfully conducted to predict the dynamic characteristics of hydrodynamic damper. The water resistance, the pressure in the interface and the stress of structure were investigated, and the relationship among the peak of water resistance, initial velocity and actual draught was also discussed. The empirical formula was put forward to predict the water resistance. And it is found that the resistance coefficient is commonly in the range of 0.3 ≤ C ≤ 0.5, when the initial velocity is larger than 50 m/s. It can be seen that the SPH method coupled with the FE method has many obvious advantages over other numerical methods for this complicated flow problem with fluid-structure interaction.

Partitioned Method of Insect Flapping Flight for Maneuvering Analysis

This study proposed a partitioned method to analyze maneuvering of insects during flapping flight.This method decomposed the insect flapping flight into wing and body subsystems and then coupled them via boundary conditions imposed on the wing’s base using one-way coupling.In the wing subsystem,the strong coupling of the flexible wings and surrounding fluid was accurately analyzed using the finite element method to obtain the thrust forces acting on the insect’s body.The resulting thrust forces were passed from the wing subsystem to the body subsystem,and then rigid body motion was analyzed in the body subsystem.The rolling,yawing,and pitching motions were simulated using the proposed method as follows:In the rolling simulation,the difference of the stroke angle between the right and left wings caused a roll torque.In the yawing simulation,the initial feathering angle in the right wing only caused a yaw torque.In the pitching simulation,the difference between the front-and back-stroke angles in both the right and left wings caused a pitch torque.All three torques generated maneuvering motion comparable with that obtained in actual observations of insect flight.These results demonstrate that the proposed method can adequately simulate the fundamental maneuvers of insect flapping flight.In the present simulations,the maneuvering mechanisms were investigated at the governing equation level,which might be difficult using other approaches.Therefore,the proposed method will contribute to revealing the underlying insect flight mechanisms.

A BEM/FEM Coupling Approach for Fluid-Structure Interaction Simulation of Cell Motion

In this paper, accurate and efficient simulation of cell motion in a biological fluid flow is investigated. The membrane of a moving cell is represented by athin shell composed of incompressible neo-Hookean elastic materials and the liquidsaround the membrane are approximated as incompressible Newtonian flows with lowReynolds numbers. The biofluid mechanics is approximated by the Stokes flow equations. A low-order BEM model is developed for the two biological fluids coupled atthe membrane surface. The moving boundary problem in fluid mechanics can be effectively solved using the BEM with a GMRES solver. The FEM model based on a flatthin shell element is further developed to predict the membrane load due to the largedeformation of a moving cell. Computational efficiency is greatly improved due tothe one-dimensional reduction in the present BEM and FEM models. The BEM solverfor the biological fluids is coupled with the FEM solver for the cell membrane at themembrane surface. The position of the membrane surface nodes is advanced in time byusing the classical fourth-order Runge-Kutta method. Numerical instability is avoidedby using a relatively small time step. Further numerical instabilities in the FEM solveris alleviated by using various techniques. The present method is applied to the FSIproblems of cell motion in a cylindrical flow. Numerical examples can illustrate thedistinct accuracy, efficiency and robustness of the present method. Furthermore, theimportance of bending stiffness of a cell membrane for stable cell motion simulation isemphasized. It is suggested that the present approach be an appealing alternative forsimulating the fluid-structure interaction of moving cells.

储液容器跌落事故的有限元分析

考虑了流固耦合效应在储液容器跌落分析中的作用,采用耦合欧拉-拉格朗日(CEL)方法模拟液体在储液容器跌落过程中的惯性效应以及液体对容器所产生的侧向液动压力,综合考虑了流固耦合效应对储液容器跌落过程中容器本身的变形和动态响应的影响。结果表明,该方法可以较为准确地评价容器的安全性和结构设计的合理性,也为同类结构的设计评价提供参考。

充液系统液体-多体耦合动力响应分析

提出了充液系统的液体-多体耦合力学模型,基于ALE有限元法和多体系统动力学理论,发展了液体-多体耦合动力响应分析的一种有效方法.对于液体子系统,将其运动分解为随同贮箱的大位移运动和相对贮箱的大幅晃动,引入贮箱固连参考系中的任意拉格朗日-欧拉(ALE)运动学描述,建立了贮箱固连非惯性参考系中液体的ALE有限元方程,对液体有限元方程的缩聚大大减少了液体子系统的计算规模.为了计及液体阻尼的影响,引入了液体修正的Rayleigh阻尼,避免了伪阻尼力的出现.对于多体子系统,应用多体系统动力学理论建立动力学方程.在此基础上详细导出了液体-多体耦合动力学方程,并采用预估-多重校正算法(PMA)和时间步长控制算法进行迭代求解,既保证了迭代收敛,又提高了计算效率.所给算例成功求解了液体运输车辆系统的液体-多体耦合动力响应,深入分析了有关参数对系统动力响应的影响,获得了一些结论.