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 beh...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.展开更多
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 s...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.展开更多
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 th...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.展开更多
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 ...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.展开更多
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), Lomb...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.展开更多
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...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.展开更多
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 techn...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.展开更多
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 th...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.展开更多
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 e...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.展开更多
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 meth...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.展开更多
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,wh...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。展开更多
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 imp...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(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 mode...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.展开更多
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 fille...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.展开更多
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 phenom...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.展开更多
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 cond...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.展开更多
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 ma...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.展开更多
文摘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.
文摘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.
文摘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.
文摘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.
文摘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.
文摘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.
基金by Key Aviation Scientific and Technological Laboratory of High-speed Hydrodynamic under grant MJ-2015-F-028.
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China(51379125,51490675,11432009,51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of the Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for Doctoral Students
文摘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.
基金supported by the National Numerical Wind Tunnel Project (Grant No. NNW2019ZT2-B02)the National Natural Science Foundation of China (Grant Nos. 12032002,51779003,and 11902005)the SinoGerman Mobility Programme (Grant No. M-0210)。
文摘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.
文摘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。
基金supported by the National Basic Research Program of China (973 Program,Grant No.2007CB206900)
文摘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(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.
基金supported by the National Natural Science Foundation of China (Grant No. 30971528)the Shanghai Committee of Science and Technology of China (Grant No.10ZR1403500)supported by the Graduate Innovation Fund of Fudan University (Grant No.EYH2126023)
文摘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.
文摘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.
文摘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.
文摘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.