This paper proposes 2.5-dimensional polymer micromachined insect-mimetic wings based on a fluid-structure interaction(FSI)design concept that enables natural deformations like cambering and pitching under fluid forces...This paper proposes 2.5-dimensional polymer micromachined insect-mimetic wings based on a fluid-structure interaction(FSI)design concept that enables natural deformations like cambering and pitching under fluid forces.Instead of directly employing an analysis for the FSI,an iterative structural Design Window(DW)search is used to reduce the computational cost significantly.A DW search using the iterative method refines the initial design by addressing fabrication challenges and tuning it to meet manufacturability constraints.The successful fabrication and demonstration of the final design solution for a wing demonstrates the effectiveness of the iterative DW search based on the FSI design concept.Furthermore,a pixel model is introduced to convert an unstructured to a structured mesh for the FSI analysis to further reduce the computational cost.The camber and pitching error between the unstructured and structured meshes is minimized to achieve insect-like aerodynamic performance by adjusting the elastic moduli of center and root veins.Finally,an analysis for the FSI is conducted,based on the parameters obtained from the pixel model to evaluate the flight performance on the basis of the lift,camber,and pitching required by an actual insect to maneuver and hover.展开更多
Submarine pipelines are critical infrastructures for offshore energy transport and communications. Understanding their structural response to near-field explosions is crucial for enhancing their blast resistance and o...Submarine pipelines are critical infrastructures for offshore energy transport and communications. Understanding their structural response to near-field explosions is crucial for enhancing their blast resistance and operational safety. This study presents a computational study on the interaction between explosion-induced bubbles and a seabed-mounted pipeline. A recently developed computational framework is employed, which couples a compressible fluid solver with a finite element structural solver via a partitioned procedure. An embedded boundary method and a level-set method are employed to handle the fluid-structure and gas-liquid interfaces. Using this framework, we analyze the flow field evolution, bubble dynamics, and transient pipe deformation. Two distinct response modes are identified: periodic oscillation under low-pressure loading and downward collapse triggered by high-pressure loading and bubble jet impact. Specifically, under high-pressure conditions, the pipe initially deforms inward, generating a localized high-pressure zone within the concave region. During structural rebound, the trapped fluid is expelled upward, giving rise to a bubble jet. Further parametric studies on the pipe's internal pressure, wall thickness, and support angle reveal several key insights. A higher internal pressure delays structural collapse, and a greater pipe thickness results in more uniform implosion morphologies. The support angle strongly influences the collapse dynamics, with the shortest collapse time occurring at 60 °. These findings offer new insights for the protective design of submarine pipelines.展开更多
The flow control at low Reynolds numbers is one of the most promising technologies in the field of aerodynamics,and it is also an important source of the innovation for novel aircraft.In this study,a new way of nonlin...The flow control at low Reynolds numbers is one of the most promising technologies in the field of aerodynamics,and it is also an important source of the innovation for novel aircraft.In this study,a new way of nonlinear flow control by interaction between two flexible flaps is proposed,and their flow control mechanism is studied employing the self-constructed immersed boundary-lattice Boltzmann-finite element method(IB-LB-FEM).The effects of the difference in material properties and flap length between the two flexible flaps on the nonlinear flow control of the airfoil are discussed.It is suggested that the relationship between the deformation of the two flexible flaps and the evolution of the vortex under the fluid-structure interaction(FSI).It is shown that the upstream flexible flap plays a key role in the flow control of the two flexible flaps.The FSI effect of the upstream flexible flap will change the unsteady flow behind it and affect the deformation of the downstream flexible flap.Two flexible flaps with different material properties and different lengths will change their own FSI characteristics by the induced vortex,effectively suppressing the flow separation on the airfoil’s upper surface.The interaction of two flexible flaps plays an extremely important role in improving the autonomy and adjustability of flow control.The numerical results will provide a theoretical basis and technical guidance for the development and application of a new flap passive control technology.展开更多
The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on mo...The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.展开更多
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.展开更多
This paper is concerned with the computational results of two-dimensional axisymmetric rigid and elastic wall formulation. In this paper, steady flow in a stenotic vessel is simulated and compared to available numeric...This paper is concerned with the computational results of two-dimensional axisymmetric rigid and elastic wall formulation. In this paper, steady flow in a stenotic vessel is simulated and compared to available numerical data with COMSOL Multiphysics software. Numerical results for a 2D axisymmetric vessel of 45% area reduction indicate that as the area is reduced with the decreasing of cross-section, the maximum axial velocity at post stenotic decreases until the end of the artery but the radial velocity increases upto 4 mm from the stenosis throat and then decreases. Overall, comparison is carried out on hemodynamics for elastic and rigid wall of steady flow. Our investigated findings may enable risk factor for patients with attacked cardiovascular diseases and can play an important role to detect a solution to such kinds of diseases.展开更多
The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional can...The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.展开更多
The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstr...The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.展开更多
The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more ...The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.展开更多
The inflation of a five-ring cone parachute with the airflow velocity of 18 m/s is studied based on the simplified arbitrary Lagrange Euler （SALE）/fluid-structure interaction （FSI） method. The numerical results of...The inflation of a five-ring cone parachute with the airflow velocity of 18 m/s is studied based on the simplified arbitrary Lagrange Euler （SALE）/fluid-structure interaction （FSI） method. The numerical results of the canopy shape, stability, opening load, and drag area are obtained, and they are well consistent with the experimental data gained from wind tunnel tests. The method is then used to simulate the opening process under different velocities. It is found that the first load shock affected by the velocity often occurs at the end of the initial inflation stage. For the first time, the phenomena that the inflation distance proportion coefficient increases and the dynamic load coefficient decreases, respectively, with the increase in the velocity are revealed. The above proposed method is competent to solve the large deformation problem without empirial coefficients, and can collect more space-time details of fluid-structure-motion information when it is compared with the traditional method.展开更多
A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic a...A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic and structural equations are marched simultaneously in time. The three-dimensional （3D） unsteady Reynolds average Navier-Stokes （RANS） equations are solved with a multiblock finite volume scheme on dynamic deforming grids to evaluate the aerodynamic force. Dual time-stepping technique and an efficient implicit scheme with multigrid are employed to march the solution in time. The blade vibration is modeled with an aeroelasticity model in which blade motion is computed by linear combination of responses of each mode under unsteady loads. The code is validated in prediction of the unsteady flow flutter behavior of an oscillating cascade and is applied to flutter analysis of a transonic fan at the design speed.展开更多
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.展开更多
Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication ...Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.展开更多
Paravalvular Leakage(PVL)has been recognized as one of the most dangerous complications in relation to Transcathether Aortic Valve Implantation(TAVI)activities.However,data available in the literature about Fluid Stru...Paravalvular Leakage(PVL)has been recognized as one of the most dangerous complications in relation to Transcathether Aortic Valve Implantation(TAVI)activities.However,data available in the literature about Fluid Structure Interaction(FSI)for this specific problem are relatively limited.In the present study,the fluid and structure responses of the hemodynamics along the patient aorta model and the aortic wall deformation are studied with the aid of numerical simulation taking into account PVL and 100%TAVI valve opening.In particular,the aorta without valve(AWoV)is assumed as the normal condition,whereas an aorta with TAVI 26 mm for 100%Geometrical Orifice Area(GOA)is considered as the patient aorta with PVL complication.A 3D patient-specific aorta model is elaborated using the MIMICS software.Implantation of the identical TAVI valve of Edward SAPIEN XT 26(Edwards Lifes ciences,Irvine,California)is considered.An undersized 26 mm TAVI valve with 100%valve opening is selected to mimic the presence of PVL at the aortic annulus.The present research indicates that the existence of PVL can increase the blood velocity,pressure drop and WSS in comparison to normal conditions,thereby paving the way to the development of recirculation flow,thrombus formation,aorta wall collapse,aortic rupture and damage of endothelium.展开更多
The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid ...The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.展开更多
This paper presents a monolithic approach to the thermal fluidstructure interaction （FSI） with nonconforming interfaces. The thermal viscous flow is governed by the Boussinesq approximation and the incompressible Na...This paper presents a monolithic approach to the thermal fluidstructure interaction （FSI） with nonconforming interfaces. The thermal viscous flow is governed by the Boussinesq approximation and the incompressible NavierStokes equations. The motion of the fluid domain is accounted for by an arbitrary LagrangianEulerian （ALE） strategy. A pseudosolid formulation is used to manage the deformation of the fluid do main. The structure is described by the geometrically nonlinear thermoelastic dynamics. An efficient data transfer strategy based on the Gauss points is proposed to guarantee the equilibrium of the stresses and heat along the interface. The resulting strongly coupled set of nonlinear equations for the fluid, solution procedure. A numerical example efficiency of the methodology. structure, and heat is solved by a monolithic is presented to demonstrate the robustness and展开更多
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.展开更多
A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in t...A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.展开更多
In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images ...In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.展开更多
Elevated intraocular pressure appears to have a broader impact on increased resistance to aqueous humor outflow through the conventional aqueous outflow system(AOS).However,there is still no consensus about exact loca...Elevated intraocular pressure appears to have a broader impact on increased resistance to aqueous humor outflow through the conventional aqueous outflow system(AOS).However,there is still no consensus about exact location of the increased outflow resistance of aqueous humor,and the mechanism is not perfect.In addition,it is difficult to accurately obtain hydrodynamic parameters of aqueous humor within the trabecular meshwork outflow pathways based on the current technology.In this paper,a two-way fluid-structure interaction simulation was performed to study the pressure difference and velocity in the superficial trabecular meshwork,juxtacanalicular meshwork(JCM)and Schlemm’s canal in response to JCM permeability changes.We obtained the JCM permeability of normal intraocular pressure varied between 1×10?15 m2 and 10×10?15 m2 while permeability of the JCM ranged from 2×10?16 m2 and 3×10?16 m2 under conditions of high intraocular pressure.The study indicated that the fluid dynamics parameters in trabecular meshwork and Schlemm’s canal are most significantly affected by the changes of JCM permeability.Moreover,the study demonstrates that the finite element modeling of AOS provides a practical means for studying the outflow dynamics and the biomechanical environment of the AOS.展开更多
基金supported by the Japan Society for the Promotion of Science KAKENHI under grant number 23H00475.
文摘This paper proposes 2.5-dimensional polymer micromachined insect-mimetic wings based on a fluid-structure interaction(FSI)design concept that enables natural deformations like cambering and pitching under fluid forces.Instead of directly employing an analysis for the FSI,an iterative structural Design Window(DW)search is used to reduce the computational cost significantly.A DW search using the iterative method refines the initial design by addressing fabrication challenges and tuning it to meet manufacturability constraints.The successful fabrication and demonstration of the final design solution for a wing demonstrates the effectiveness of the iterative DW search based on the FSI design concept.Furthermore,a pixel model is introduced to convert an unstructured to a structured mesh for the FSI analysis to further reduce the computational cost.The camber and pitching error between the unstructured and structured meshes is minimized to achieve insect-like aerodynamic performance by adjusting the elastic moduli of center and root veins.Finally,an analysis for the FSI is conducted,based on the parameters obtained from the pixel model to evaluate the flight performance on the basis of the lift,camber,and pitching required by an actual insect to maneuver and hover.
基金supported by the National Key R&D Program of China(Grant No.2024YFC3013200)the Shenzhen Peacock Plan(Grant No.QD2023006C).
文摘Submarine pipelines are critical infrastructures for offshore energy transport and communications. Understanding their structural response to near-field explosions is crucial for enhancing their blast resistance and operational safety. This study presents a computational study on the interaction between explosion-induced bubbles and a seabed-mounted pipeline. A recently developed computational framework is employed, which couples a compressible fluid solver with a finite element structural solver via a partitioned procedure. An embedded boundary method and a level-set method are employed to handle the fluid-structure and gas-liquid interfaces. Using this framework, we analyze the flow field evolution, bubble dynamics, and transient pipe deformation. Two distinct response modes are identified: periodic oscillation under low-pressure loading and downward collapse triggered by high-pressure loading and bubble jet impact. Specifically, under high-pressure conditions, the pipe initially deforms inward, generating a localized high-pressure zone within the concave region. During structural rebound, the trapped fluid is expelled upward, giving rise to a bubble jet. Further parametric studies on the pipe's internal pressure, wall thickness, and support angle reveal several key insights. A higher internal pressure delays structural collapse, and a greater pipe thickness results in more uniform implosion morphologies. The support angle strongly influences the collapse dynamics, with the shortest collapse time occurring at 60 °. These findings offer new insights for the protective design of submarine pipelines.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.92371201,52192633,11872293,and 92152301)the Natural Science Basic Research Program of Shaanxi(Grant Nos.2024JC-YBQN-0008,and 2022JC-03)+1 种基金Shaanxi Key Research and Development Program(Grant No.2022ZDLGY02-07)the Joint Natural Science Foundation of China with Guangdong Province for TianHe-II Supercomputer Resources,and the Research Start-up Foundation of Xi’an University of Science and Technology for the High-Level Talent.
文摘The flow control at low Reynolds numbers is one of the most promising technologies in the field of aerodynamics,and it is also an important source of the innovation for novel aircraft.In this study,a new way of nonlinear flow control by interaction between two flexible flaps is proposed,and their flow control mechanism is studied employing the self-constructed immersed boundary-lattice Boltzmann-finite element method(IB-LB-FEM).The effects of the difference in material properties and flap length between the two flexible flaps on the nonlinear flow control of the airfoil are discussed.It is suggested that the relationship between the deformation of the two flexible flaps and the evolution of the vortex under the fluid-structure interaction(FSI).It is shown that the upstream flexible flap plays a key role in the flow control of the two flexible flaps.The FSI effect of the upstream flexible flap will change the unsteady flow behind it and affect the deformation of the downstream flexible flap.Two flexible flaps with different material properties and different lengths will change their own FSI characteristics by the induced vortex,effectively suppressing the flow separation on the airfoil’s upper surface.The interaction of two flexible flaps plays an extremely important role in improving the autonomy and adjustability of flow control.The numerical results will provide a theoretical basis and technical guidance for the development and application of a new flap passive control technology.
基金support of the National Natural Science Foundation of China(No.51675406)the Basic Research Project Group,China(No.514010106-205)。
文摘The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.
基金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.
文摘This paper is concerned with the computational results of two-dimensional axisymmetric rigid and elastic wall formulation. In this paper, steady flow in a stenotic vessel is simulated and compared to available numerical data with COMSOL Multiphysics software. Numerical results for a 2D axisymmetric vessel of 45% area reduction indicate that as the area is reduced with the decreasing of cross-section, the maximum axial velocity at post stenotic decreases until the end of the artery but the radial velocity increases upto 4 mm from the stenosis throat and then decreases. Overall, comparison is carried out on hemodynamics for elastic and rigid wall of steady flow. Our investigated findings may enable risk factor for patients with attacked cardiovascular diseases and can play an important role to detect a solution to such kinds of diseases.
基金Project supported by the National Natural Science Foundation of China(No.10832007)the Shanghai Leading Academic Discipline Project(No.B206)
文摘The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.
文摘The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.
基金Project supported by the National Natural Science Foundation of China(No.11972112)the Fundamental Research Funds for the Central Universities of China(Nos.N2103024 and N2103002)the Major Projects of Aero-Engines and Gasturbines(No.J2019-I-0008-0008)。
文摘The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.
基金supported by the National Natural Science Foundation of China(No.11172137)the Aeronautical Science Foundation of China(No.20122910001)
文摘The inflation of a five-ring cone parachute with the airflow velocity of 18 m/s is studied based on the simplified arbitrary Lagrange Euler （SALE）/fluid-structure interaction （FSI） method. The numerical results of the canopy shape, stability, opening load, and drag area are obtained, and they are well consistent with the experimental data gained from wind tunnel tests. The method is then used to simulate the opening process under different velocities. It is found that the first load shock affected by the velocity often occurs at the end of the initial inflation stage. For the first time, the phenomena that the inflation distance proportion coefficient increases and the dynamic load coefficient decreases, respectively, with the increase in the velocity are revealed. The above proposed method is competent to solve the large deformation problem without empirial coefficients, and can collect more space-time details of fluid-structure-motion information when it is compared with the traditional method.
文摘A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic and structural equations are marched simultaneously in time. The three-dimensional （3D） unsteady Reynolds average Navier-Stokes （RANS） equations are solved with a multiblock finite volume scheme on dynamic deforming grids to evaluate the aerodynamic force. Dual time-stepping technique and an efficient implicit scheme with multigrid are employed to march the solution in time. The blade vibration is modeled with an aeroelasticity model in which blade motion is computed by linear combination of responses of each mode under unsteady loads. The code is validated in prediction of the unsteady flow flutter behavior of an oscillating cascade and is applied to flutter analysis of a transonic fan at the design speed.
基金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 Foundations of China(Nos.50821063 and 50823004)973 Program(No.2007CB714701)the Fundamental Research Funds for the Central Universities(No.2010XS34)
文摘Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.
基金Universiti Putra Malaysia,for providing funds for this project through Grant UPM GP-IPM/2019/9675000.
文摘Paravalvular Leakage(PVL)has been recognized as one of the most dangerous complications in relation to Transcathether Aortic Valve Implantation(TAVI)activities.However,data available in the literature about Fluid Structure Interaction(FSI)for this specific problem are relatively limited.In the present study,the fluid and structure responses of the hemodynamics along the patient aorta model and the aortic wall deformation are studied with the aid of numerical simulation taking into account PVL and 100%TAVI valve opening.In particular,the aorta without valve(AWoV)is assumed as the normal condition,whereas an aorta with TAVI 26 mm for 100%Geometrical Orifice Area(GOA)is considered as the patient aorta with PVL complication.A 3D patient-specific aorta model is elaborated using the MIMICS software.Implantation of the identical TAVI valve of Edward SAPIEN XT 26(Edwards Lifes ciences,Irvine,California)is considered.An undersized 26 mm TAVI valve with 100%valve opening is selected to mimic the presence of PVL at the aortic annulus.The present research indicates that the existence of PVL can increase the blood velocity,pressure drop and WSS in comparison to normal conditions,thereby paving the way to the development of recirculation flow,thrombus formation,aorta wall collapse,aortic rupture and damage of endothelium.
基金the financial support provided by the Office of Naval Research(ONR) through grant number N00014-09-1-1204 (managed by Dr. Ki-Han Kim)supported in part by the National Research Foundation of Korea (NRF)grant funded by the Korea government (MEST) through the GCRC-SOP Grant No. 2012-0004783
文摘The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.
文摘This paper presents a monolithic approach to the thermal fluidstructure interaction （FSI） with nonconforming interfaces. The thermal viscous flow is governed by the Boussinesq approximation and the incompressible NavierStokes equations. The motion of the fluid domain is accounted for by an arbitrary LagrangianEulerian （ALE） strategy. A pseudosolid formulation is used to manage the deformation of the fluid do main. The structure is described by the geometrically nonlinear thermoelastic dynamics. An efficient data transfer strategy based on the Gauss points is proposed to guarantee the equilibrium of the stresses and heat along the interface. The resulting strongly coupled set of nonlinear equations for the fluid, solution procedure. A numerical example efficiency of the methodology. structure, and heat is solved by a monolithic is presented to demonstrate the robustness and
文摘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.
文摘A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.
基金The project supported by the National Natural Science Foundation of China(10672036,10472025 and 10421002)the Natural Science Foundation of Liaoning Province(20032109)
文摘In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.
基金This work study is financially supported by the National Natural Science Foundation of China(Nos.31570952,10802053,81471702)the Natural Science Foundation of Beijing(No.3122010)。
文摘Elevated intraocular pressure appears to have a broader impact on increased resistance to aqueous humor outflow through the conventional aqueous outflow system(AOS).However,there is still no consensus about exact location of the increased outflow resistance of aqueous humor,and the mechanism is not perfect.In addition,it is difficult to accurately obtain hydrodynamic parameters of aqueous humor within the trabecular meshwork outflow pathways based on the current technology.In this paper,a two-way fluid-structure interaction simulation was performed to study the pressure difference and velocity in the superficial trabecular meshwork,juxtacanalicular meshwork(JCM)and Schlemm’s canal in response to JCM permeability changes.We obtained the JCM permeability of normal intraocular pressure varied between 1×10?15 m2 and 10×10?15 m2 while permeability of the JCM ranged from 2×10?16 m2 and 3×10?16 m2 under conditions of high intraocular pressure.The study indicated that the fluid dynamics parameters in trabecular meshwork and Schlemm’s canal are most significantly affected by the changes of JCM permeability.Moreover,the study demonstrates that the finite element modeling of AOS provides a practical means for studying the outflow dynamics and the biomechanical environment of the AOS.
