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 s...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.展开更多
We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive ...We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.展开更多
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 phenome...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.展开更多
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.展开更多
This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier...This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.展开更多
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 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 perfo...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.展开更多
Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The b...Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
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.展开更多
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.展开更多
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.展开更多
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 ...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.展开更多
The fluid-structure interaction(FSI)in aircraft hydraulic pipeline systems is of great concern because of the damage it causes.To accurately predict the vibration characteristic of long hydraulic pipelines with curved...The fluid-structure interaction(FSI)in aircraft hydraulic pipeline systems is of great concern because of the damage it causes.To accurately predict the vibration characteristic of long hydraulic pipelines with curved segments,we studied the frequency-domain modeling and solution method for FSI in these pipeline systems.Fourteen partial differential equations(PDEs)are utilized to model the pipeline FSI,considering both frequency-dependent friction and bending-flexibility modification.To address the numerical instability encountered by the traditional transfer matrix method(TMM)in solving relatively complex pipelines,an improved TMM is proposed for solving the PDEs in the frequency domain,based on the matrix-stacking strategy and matrix representation of boundary conditions.The proposed FSI model and improved solution method are validated by numerical cases and experiments.An experimental rig of a practical hydraulic system,consisting of an aircraft engine-driven pump,a Z-shaped aero-hydraulic pipeline,and a throttle valve,was constructed for testing.The magnitude ratio of acceleration to pressure is introduced to evaluate the theoretical and experimental results,which indicate that the proposed model and solution method are effective in practical applications.The methodology presented in this paper can be used as an efficient approach for the vibrational design of aircraft hydraulic pipeline systems.展开更多
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 conditi...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.展开更多
In this work,we employ fluid-structure interaction(FSI)systems with immersed flexible structures with or without free surfaces to explore both Singular Value Decomposition(SVD)-based model reduction methods and mode s...In this work,we employ fluid-structure interaction(FSI)systems with immersed flexible structures with or without free surfaces to explore both Singular Value Decomposition(SVD)-based model reduction methods and mode superposition methods.For acoustoelastic FSI systems,we adopt a three-field mixed finite element formulation with displacement,pressure,and vorticity moment unknowns to effectively enforce the irrotationality constraint.We also propose in this paper a new Inf-Sup test based on the lowest non-zero singular value of the coupling matrix for the selection of reliable sets of finite element discretizations for displacement and pressure as well as vorticity moment.Our numerical examples demonstrate that mixed finite element formulations can be effectively used to predict resonance frequencies of fully coupled FSI systems within different ranges of respective physical motions,namely,acoustic,structural,and slosh motions,without the contamination of spurious(non-physical)modes with nonzero frequencies.Our numerical results also confirm that SVD-based model reduction methods can be effectively used to reconstruct from a few snapshots of transient solutions the dominant principal components with moderate level of signal to noise ratio,which may eventually open doors for simulation of long-term behaviors of both linear and nonlinear FSI systems.展开更多
During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing str...During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing structural fatigue.In this study,a forced vibration analysis of tubing under CO_(2) fracturing conditions is carried out by taking into account the fluid-structure coupling and related interaction forces by means of the method of characteristics(MOC).The results show that for every 1 m^(3)/min increase in pumping displacement,the fluid flow rate increases up to 3.67 m/s.The flow pressure in the pipe tends to be consistent with the pumping pressure at the initial stage and then decreases with an increase in the pump starting time.When the pumping pressure increases by 10 MPa,the additional stress in the tubing increases by 11.8%,and the peak value of the additional stress at the bottom of the well is the largest.The temperature in the tubing grows with well depth,which causes a phase change in CO_(2) due to heat absorption.At this time the pressure in the tubing decreases,the fluid flow rate increases by about 1.12 m/s,and the additional stress grows by about 1.5 MPa.展开更多
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.展开更多
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.展开更多
文摘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.
基金This work was supported by the National Natural Science Foundation of China(Grant No.12172160)Shenzhen Science and Technology Program(Grant No.JCYJ20220818100600002)+1 种基金South-ern University of Science and Technology(Grant No.Y01326127)the Department of Science and Technology of Guangdong Province(Grant Nos.2020B1212030001 and 2021QN020642).
文摘We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.
文摘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.
文摘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.
文摘This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.
文摘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.
基金National Natural Science Foundation of China under Grant Nos.51427901 and 51678407Tianjin Municipal Education Commission under Grant No.2021KJ055Fundamental Research Funds for the Central Universities of China under Grant No.2000560616。
文摘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.
文摘Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
文摘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.
文摘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.
文摘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.
基金Project supported by the Program for New Century Excellent Talents in University (No. NCET-10-0054)the Fok Ying-Tong Education Foundation,China (No. 121073)+1 种基金the National Natural Science Foundation of China (No. 10976008)the State Key Program of National Natural Science of China (No. 50939002)
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.51975025 and 51890822)the Young Elite Scientists Sponsorship Program by China Association for Science and Technology(No.2016QNRC001)the National Key Research and Development Program of China(No.2019YFB2004500)。
文摘The fluid-structure interaction(FSI)in aircraft hydraulic pipeline systems is of great concern because of the damage it causes.To accurately predict the vibration characteristic of long hydraulic pipelines with curved segments,we studied the frequency-domain modeling and solution method for FSI in these pipeline systems.Fourteen partial differential equations(PDEs)are utilized to model the pipeline FSI,considering both frequency-dependent friction and bending-flexibility modification.To address the numerical instability encountered by the traditional transfer matrix method(TMM)in solving relatively complex pipelines,an improved TMM is proposed for solving the PDEs in the frequency domain,based on the matrix-stacking strategy and matrix representation of boundary conditions.The proposed FSI model and improved solution method are validated by numerical cases and experiments.An experimental rig of a practical hydraulic system,consisting of an aircraft engine-driven pump,a Z-shaped aero-hydraulic pipeline,and a throttle valve,was constructed for testing.The magnitude ratio of acceleration to pressure is introduced to evaluate the theoretical and experimental results,which indicate that the proposed model and solution method are effective in practical applications.The methodology presented in this paper can be used as an efficient approach for the vibrational design of aircraft hydraulic pipeline systems.
基金Supported by the National Natural Science Foundation of China (50774064)
文摘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.
文摘In this work,we employ fluid-structure interaction(FSI)systems with immersed flexible structures with or without free surfaces to explore both Singular Value Decomposition(SVD)-based model reduction methods and mode superposition methods.For acoustoelastic FSI systems,we adopt a three-field mixed finite element formulation with displacement,pressure,and vorticity moment unknowns to effectively enforce the irrotationality constraint.We also propose in this paper a new Inf-Sup test based on the lowest non-zero singular value of the coupling matrix for the selection of reliable sets of finite element discretizations for displacement and pressure as well as vorticity moment.Our numerical examples demonstrate that mixed finite element formulations can be effectively used to predict resonance frequencies of fully coupled FSI systems within different ranges of respective physical motions,namely,acoustic,structural,and slosh motions,without the contamination of spurious(non-physical)modes with nonzero frequencies.Our numerical results also confirm that SVD-based model reduction methods can be effectively used to reconstruct from a few snapshots of transient solutions the dominant principal components with moderate level of signal to noise ratio,which may eventually open doors for simulation of long-term behaviors of both linear and nonlinear FSI systems.
基金funded by National Natural Science Foundation of China(Grant No.52105209)funded by the National Engineering Research Centre for Oil and Gas Drilling and Completion Technology(Grant No.F2023199).
文摘During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing structural fatigue.In this study,a forced vibration analysis of tubing under CO_(2) fracturing conditions is carried out by taking into account the fluid-structure coupling and related interaction forces by means of the method of characteristics(MOC).The results show that for every 1 m^(3)/min increase in pumping displacement,the fluid flow rate increases up to 3.67 m/s.The flow pressure in the pipe tends to be consistent with the pumping pressure at the initial stage and then decreases with an increase in the pump starting time.When the pumping pressure increases by 10 MPa,the additional stress in the tubing increases by 11.8%,and the peak value of the additional stress at the bottom of the well is the largest.The temperature in the tubing grows with well depth,which causes a phase change in CO_(2) due to heat absorption.At this time the pressure in the tubing decreases,the fluid flow rate increases by about 1.12 m/s,and the additional stress grows by about 1.5 MPa.
文摘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.
文摘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.
