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.展开更多
In the experimental investigation of fluid-structure interactions regarding the undulatory motion like flag flapping or fish swimming,solving the force distribution on the flexible body stands as an indispensable ende...In the experimental investigation of fluid-structure interactions regarding the undulatory motion like flag flapping or fish swimming,solving the force distribution on the flexible body stands as an indispensable endeavor to gain insights into the underlying dynamic mechanisms.However,the solving process entails high-order numerical derivatives of experimental data,which poses a formidable challenge for experimental studies on fluid-structure interactions,given that the measurement noise inherent in experimental data renders the problem ill-posed.The commonly practiced regularization methods for numerical derivatives are feeble to tackle the fourth-order derivative associated with the bending force;those methods,in particular,require predetermined parameters about the unknown noise.We introduce here an empirical regularization method founded upon the kernel-term modification in the frequency domain,notably capable of determining the fourth derivative of experimental data.By leveraging the potentials of the iterative operations,our method enables the reliable estimation of an approximately optimal regularization parameter,all without reliance on any a priori knowledge about the noise characteristics.To demonstrate the reliability,robustness,and accuracy of the method,we perform rigorous numerical assessments using different data models that are infused with noise varying several orders of magnitude.Additionally,practical application of this method is achieved in the experiment on a flexible film flapping in the gusty flow,where the spatiotemporal distribution of the bending force density on the film is calculated by integrating this method with a linear reconstruction.展开更多
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.展开更多
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 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.展开更多
In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena （e.g. water waves, impacts, splashing jets, etc.）. To model these complex free surface...In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena （e.g. water waves, impacts, splashing jets, etc.）. To model these complex free surface flows is a tough and challenging task for most computational fluid dynamics （CFD） solvers which work in the Eulerian framework. As a Lagrangian and meshless method, smoothed particle hydrodynamics （SPH） offers a convenient tracking for different complex boundaries and a straightforward satisfaction for different boundary conditions. Therefore SPH is robust in modeling complex hydrodynamic problems characterized by free surface boundaries, multiphase interfaces or material discontinuities. Along with the rapid development of the SPH theory, related numerical techniques and high-performance computing technologies, SPH has not only attracted much attention in the academic community, but also gradually gained wide applications in industrial circles. This paper is dedicated to a review of the recent developments of SPH method and its typical applications in fluid-structure interactions in ocean engineering. Different numerical techniques for improving numerical accuracy, satisfying different boundary conditions, improving computational efficie- ncy, suppressing pressure fluctuations and preventing the tensile instability, etc., are introduced. In the numerical results, various typical fluid-structure interaction problems or multiphase problems in ocean engineering are described, modeled and validated. The prospective developments of SPH in ocean engineering are also discussed.展开更多
Luminescent metal-organic frameworks(MOFs)have garnered significant attention due to their structural tunability and potential applications in solid-state lighting,bioimaging,sensing,anticounterfeiting,and other field...Luminescent metal-organic frameworks(MOFs)have garnered significant attention due to their structural tunability and potential applications in solid-state lighting,bioimaging,sensing,anticounterfeiting,and other fields.Nevertheless,due to the tendency of1,4-benzenedicarboxylic acid(BDC)to rotate within the framework,MOFs composed of it exhibit significant non-radiative energy dissipation and thus impair the emissive properties.In this study,efficient luminescence of MIL-140A nanocrystals(NCs)with BDC rotors as ligands is achieved by pressure treatment strategy.Pressure treatment effectively modulates the pore structure of the framework,enhancing the interactions between the N,N-dimethylformamide vip molecules and the BDC ligands.The enhanced host-vip interaction contributes to the structural rigidity of the MOF,thereby suppressing the rotation-induced excited-state energy loss.As a result,the pressure-treated MIL-140A NCs displayed bright blue-light emission,with the photoluminescence quantum yield increasing from an initial 6.8%to 69.2%.This study developed an effective strategy to improve the luminescence performance of rotor ligand MOFs,offers a new avenue for the rational design and synthesis of MOFs with superior luminescent properties.展开更多
Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production eff...Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.展开更多
Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering s...Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.展开更多
Beryllium-containing sludge(BCS)is a typical hazardous waste from Be smelting,which can cause serious harm to ecology and human health by releasing harmful Be if it is stored long-term in environment.Nonetheless,the o...Beryllium-containing sludge(BCS)is a typical hazardous waste from Be smelting,which can cause serious harm to ecology and human health by releasing harmful Be if it is stored long-term in environment.Nonetheless,the occurrence of Be in BCS is unclear,which seriously hinders the development of pollution control technologies.In order to enhance the understanding of BCS,the occurrence of Be and the microscale interactions with coexisting phases were investigated for the first time.It was found that CaSO_(4)·2H_(2)O and amorphous SiO_(2) are the primary phases of BCS.The simulated experiments of purified materials showed that Be interacted with CaSO_(4)·2H_(2)O and amorphous SiO_(2).Be can enter into the lattice of CaSO_(4)·2H_(2)O mainly as free Be2+.Amorphous SiO_(2) can adsorb Be2+particularly at a pH range of 3–5.The dissolution behavior experiment of BCS shows that about 52%of the Be is readily extracted under acidic conditions,which refers to the Be of independent occurrence.In contrast,the remaining 48%of Be can be extracted only after the CaSO_(4)·2H_(2)O has completely dissolved.Hence,CaSO_(4)·2H_(2)O is identified as the key occurrence phase which determines the highly efficient dissolution of Be.As a result,this study enhances the understanding of BCS and lays the foundation for the development of Be separation technologies.展开更多
Investigations into the characteristics of hemodynamics will provide a better understanding of the pathology of cerebral aneurysms for clinicians.In this work,a steady state discontinuous-growth model of the cerebral ...Investigations into the characteristics of hemodynamics will provide a better understanding of the pathology of cerebral aneurysms for clinicians.In this work,a steady state discontinuous-growth model of the cerebral aneurysms was proposed.With the assumption of the fluid-structure interaction between the wall of blood vessel and blood,a fluid-structure coupling numerical simulation for this model was built using software ANSYS and CFX.The simulation results showed that as the aneurysm volume increased,a blood flow vortex came forth,the vortex region became asymptotically larger,and eddy density became gradually stronger in it.After the emergence of the vortex region,the blood flow in the vicinity of the downstream in the aneurysms volume turned into bifurcated flow,and the location of the flow bifurcated point was shifted with the aneurysm volume growing while directions of the shear stress applied to two sides of the bifurcated point were opposite.The Von Mises stress distribution along the wall of aneurysm volume decreased in the prior period and increased gradually in the later period.The maximum stress was in the neck of the volume and the minimum was on the distal end in the whole process of growth.It was shown that as the aneurysm increased the maximum deformation location of the aneurysm,vertical to the streamline,was transferred from the distal end of the aneurysm to its neck,then back to its distal end of the aneurysm again.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lat...The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lateral distance are numerically simulated using a penalty-immersed boundary method.The effects of the flapping phase and lateral distance on the propulsive performance of two fish meetings are analyzed.Results show that,when two plates meet,if their leading edges diverge laterally,the individual plate can efficiently and rapidly move apart from the other horizontally.If their leading edges converge laterally,the plate motion can be retarded,leading to high energy consumption.Moreover,an increasing lateral distance between two plates significantly weakens the fluid-structure interactions,resulting in an exponential decline in mean cruising speed.A quantitative force analysis based on vortex dynamic theory is performed to gain physics insight into the hydrodynamic interaction mechanism.It is found that lateral separation between the two leading edges enhances the vorticity generation and boundary vorticity flux on the surface of the plate,subsequently reinforcing the thrust effect and increasing horizontal velocity.This study offers insight into the hydro-dynamic mechanisms of the fluid-structure interactions among fish moving toward each other and suggests potential strategies for enhancing the maneuverability of robotic fish in complex environment.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.12127802 and 11721202)。
文摘In the experimental investigation of fluid-structure interactions regarding the undulatory motion like flag flapping or fish swimming,solving the force distribution on the flexible body stands as an indispensable endeavor to gain insights into the underlying dynamic mechanisms.However,the solving process entails high-order numerical derivatives of experimental data,which poses a formidable challenge for experimental studies on fluid-structure interactions,given that the measurement noise inherent in experimental data renders the problem ill-posed.The commonly practiced regularization methods for numerical derivatives are feeble to tackle the fourth-order derivative associated with the bending force;those methods,in particular,require predetermined parameters about the unknown noise.We introduce here an empirical regularization method founded upon the kernel-term modification in the frequency domain,notably capable of determining the fourth derivative of experimental data.By leveraging the potentials of the iterative operations,our method enables the reliable estimation of an approximately optimal regularization parameter,all without reliance on any a priori knowledge about the noise characteristics.To demonstrate the reliability,robustness,and accuracy of the method,we perform rigorous numerical assessments using different data models that are infused with noise varying several orders of magnitude.Additionally,practical application of this method is achieved in the experiment on a flexible film flapping in the gusty flow,where the spatiotemporal distribution of the bending force density on the film is calculated by integrating this method with a linear reconstruction.
基金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.
基金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.
基金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.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.U1430236,51609049)the China Postdoctoral Science Foundation(Grant No.2015M581432)the China Scholarship Council(CSC,Grant No.201506680004)
文摘In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena （e.g. water waves, impacts, splashing jets, etc.）. To model these complex free surface flows is a tough and challenging task for most computational fluid dynamics （CFD） solvers which work in the Eulerian framework. As a Lagrangian and meshless method, smoothed particle hydrodynamics （SPH） offers a convenient tracking for different complex boundaries and a straightforward satisfaction for different boundary conditions. Therefore SPH is robust in modeling complex hydrodynamic problems characterized by free surface boundaries, multiphase interfaces or material discontinuities. Along with the rapid development of the SPH theory, related numerical techniques and high-performance computing technologies, SPH has not only attracted much attention in the academic community, but also gradually gained wide applications in industrial circles. This paper is dedicated to a review of the recent developments of SPH method and its typical applications in fluid-structure interactions in ocean engineering. Different numerical techniques for improving numerical accuracy, satisfying different boundary conditions, improving computational efficie- ncy, suppressing pressure fluctuations and preventing the tensile instability, etc., are introduced. In the numerical results, various typical fluid-structure interaction problems or multiphase problems in ocean engineering are described, modeled and validated. The prospective developments of SPH in ocean engineering are also discussed.
基金supported by the National Key R&D Program of China(Grant No.2023YFA1406200)the National Natural Science Foundation of China(No.12274177 and 12304261)the China Postdoctoral Science Foundation(No.2024M751076)。
文摘Luminescent metal-organic frameworks(MOFs)have garnered significant attention due to their structural tunability and potential applications in solid-state lighting,bioimaging,sensing,anticounterfeiting,and other fields.Nevertheless,due to the tendency of1,4-benzenedicarboxylic acid(BDC)to rotate within the framework,MOFs composed of it exhibit significant non-radiative energy dissipation and thus impair the emissive properties.In this study,efficient luminescence of MIL-140A nanocrystals(NCs)with BDC rotors as ligands is achieved by pressure treatment strategy.Pressure treatment effectively modulates the pore structure of the framework,enhancing the interactions between the N,N-dimethylformamide vip molecules and the BDC ligands.The enhanced host-vip interaction contributes to the structural rigidity of the MOF,thereby suppressing the rotation-induced excited-state energy loss.As a result,the pressure-treated MIL-140A NCs displayed bright blue-light emission,with the photoluminescence quantum yield increasing from an initial 6.8%to 69.2%.This study developed an effective strategy to improve the luminescence performance of rotor ligand MOFs,offers a new avenue for the rational design and synthesis of MOFs with superior luminescent properties.
基金supported by the Chongqing Natural Science Foundation of Chongqing,China(No.CSTB2022NSCQ-MSX0333)the Science and Technology Research Program of Chongqing Municipal Education Commission(Grant No.KJZD-K202401205)+1 种基金Chongqing Three Gorges University Graduate Research and Innovation Project Funding(No.YJSKY24045)Chongqing Engineering Research Center of Disaster Prevention&Control for Banks and Structures in Three Gorges Reservoir Area(No.SXAPGC24YB14,No.SXAPGC24YB03,No.SXAPGC24YB12)。
文摘Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.
基金supported by the National Natural Science Foundation of China(Grant No.51779003)the National Key Research and Development Project of China(Grant No.2018YFB0704000)
文摘Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.
基金supported by the National Natural Science Foundation of China(No.22276219)the foundation for Innovative Research Groups of the National Natural Science Foundation of China(No.52121004)+1 种基金the major program Natural Science Foundation of Hunan Province of China(No.2021JC0001)the Fundamental Research Funds for the Central Universities of Central South University(No.2024ZZTS0063).
文摘Beryllium-containing sludge(BCS)is a typical hazardous waste from Be smelting,which can cause serious harm to ecology and human health by releasing harmful Be if it is stored long-term in environment.Nonetheless,the occurrence of Be in BCS is unclear,which seriously hinders the development of pollution control technologies.In order to enhance the understanding of BCS,the occurrence of Be and the microscale interactions with coexisting phases were investigated for the first time.It was found that CaSO_(4)·2H_(2)O and amorphous SiO_(2) are the primary phases of BCS.The simulated experiments of purified materials showed that Be interacted with CaSO_(4)·2H_(2)O and amorphous SiO_(2).Be can enter into the lattice of CaSO_(4)·2H_(2)O mainly as free Be2+.Amorphous SiO_(2) can adsorb Be2+particularly at a pH range of 3–5.The dissolution behavior experiment of BCS shows that about 52%of the Be is readily extracted under acidic conditions,which refers to the Be of independent occurrence.In contrast,the remaining 48%of Be can be extracted only after the CaSO_(4)·2H_(2)O has completely dissolved.Hence,CaSO_(4)·2H_(2)O is identified as the key occurrence phase which determines the highly efficient dissolution of Be.As a result,this study enhances the understanding of BCS and lays the foundation for the development of Be separation technologies.
基金supported by the National Natural Science Foundation of China (Grant No. 50976088)
文摘Investigations into the characteristics of hemodynamics will provide a better understanding of the pathology of cerebral aneurysms for clinicians.In this work,a steady state discontinuous-growth model of the cerebral aneurysms was proposed.With the assumption of the fluid-structure interaction between the wall of blood vessel and blood,a fluid-structure coupling numerical simulation for this model was built using software ANSYS and CFX.The simulation results showed that as the aneurysm volume increased,a blood flow vortex came forth,the vortex region became asymptotically larger,and eddy density became gradually stronger in it.After the emergence of the vortex region,the blood flow in the vicinity of the downstream in the aneurysms volume turned into bifurcated flow,and the location of the flow bifurcated point was shifted with the aneurysm volume growing while directions of the shear stress applied to two sides of the bifurcated point were opposite.The Von Mises stress distribution along the wall of aneurysm volume decreased in the prior period and increased gradually in the later period.The maximum stress was in the neck of the volume and the minimum was on the distal end in the whole process of growth.It was shown that as the aneurysm increased the maximum deformation location of the aneurysm,vertical to the streamline,was transferred from the distal end of the aneurysm to its neck,then back to its distal end of the aneurysm again.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
基金joined PI of Westlake University(Grant Nos.041030150118 and 103110556022101)Scientific Research Funding Project of Westlake University(Grant No.2021WUFP017).
文摘The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lateral distance are numerically simulated using a penalty-immersed boundary method.The effects of the flapping phase and lateral distance on the propulsive performance of two fish meetings are analyzed.Results show that,when two plates meet,if their leading edges diverge laterally,the individual plate can efficiently and rapidly move apart from the other horizontally.If their leading edges converge laterally,the plate motion can be retarded,leading to high energy consumption.Moreover,an increasing lateral distance between two plates significantly weakens the fluid-structure interactions,resulting in an exponential decline in mean cruising speed.A quantitative force analysis based on vortex dynamic theory is performed to gain physics insight into the hydrodynamic interaction mechanism.It is found that lateral separation between the two leading edges enhances the vorticity generation and boundary vorticity flux on the surface of the plate,subsequently reinforcing the thrust effect and increasing horizontal velocity.This study offers insight into the hydro-dynamic mechanisms of the fluid-structure interactions among fish moving toward each other and suggests potential strategies for enhancing the maneuverability of robotic fish in complex environment.
基金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.
基金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.
基金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.
文摘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.
基金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.
基金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.
