The Donalda gold deposit in the southern part of the Archean Abitibi greenstone belt consists mainly of a subhonzontal gold-quartz vein perpendicular to subvertical shear zones.The 0.3—0.5 m thick vein is characteriz...The Donalda gold deposit in the southern part of the Archean Abitibi greenstone belt consists mainly of a subhonzontal gold-quartz vein perpendicular to subvertical shear zones.The 0.3—0.5 m thick vein is characterized by vein-parallel banding structures indicating multiple episodes of fracture opening and mineral precipitation.Measurement of the c-axis of primary growth quartz indicates that quartz preferentially grew perpendicular to the fracture,suggesting open space filling and/or extensional nature of the fracture.Measurement of the orientations of microfractures,veinlets and fluid—inclusion planes（FIPs） crosscutting primary growth quartz indicates that the vein minerals were subject to a vertical maximum principal stress（σ_1）,which is inconsistent with the subhorizontalσ_1 inferred from the regional stress field with N—S shortening.This apparent discrepancy is explained by invoking episodic fluid pressure fluctuation between supralithostatic and hydrostatic regimes accompanied by episodic opening and closing of the sub-horizontal fracture.When fluid pressure was higher than the lithostatic value,the fracture was opened and primary growth minerals were precipitated,whereas when fluid pressure decreased toward the hydrostatic value,the hanging wall of the fracture collapsed,causing collision of protruding primary growth minerals from both sides of the fracture and resulting in formation of vein-parallel deformation bands.The columns where the two facing sides of the fracture collided were subject to higher-than-lithostatic stress due to the bridging effect and reduced support surface area,explaining the development of verticalσ_1.This hypothesis is consistent the fault-valve model,and explains the flipping ofσ_1 without having to change the regional stress field.展开更多
In order to enable wing morphing(e.g.change in camber or folds)without incurring additional weight to the aircraft,lightweight flexible materials such as membrane are needed.However,the research on fluid-structure cou...In order to enable wing morphing(e.g.change in camber or folds)without incurring additional weight to the aircraft,lightweight flexible materials such as membrane are needed.However,the research on fluid-structure coupling of membranes has mainly focused on parachutes in low-speed conditions,while that in supersonic flow conditions is lacking.Here,the degraded shell method is proposed to study membrane deformation by using shell element,which is more effective than using membrane elements directly.A fluid-structure interaction computational framework is proposed,whereby the aerodynamic module is composed of either the piston theory or computational fluid dynamics.A rectangular membrane of length 0.4 m and width 0.6 m is investigated in supersonic conditions.The characteristics of the limit cycle and steady deformation are analyzed,considering the effects of angle of attack and dynamic pressure.It is found that the structural response exhibits significant differences under various angles of attack.Furthermore,initial relaxation of membrane has significant influence on the structural deformation.Finally,the aeroelastic scaling method for membrane structures is derived,providing guidance for the design of wind tunnel models.This study provides a theoretical foundation for the analysis and application of membrane structures under supersonic conditions in future research.展开更多
Shock tunnels are indispensable facilities for hypersonic aerodynamic experimentation.Within these systems,the diaphragm plays a pivotal role,as its rupture process critically influences shock wave generation quality,...Shock tunnels are indispensable facilities for hypersonic aerodynamic experimentation.Within these systems,the diaphragm plays a pivotal role,as its rupture process critically influences shock wave generation quality,experimental repeatability,and facility reliability.A thorough understanding of diaphragm rupture dynamics is therefore essential for optimizing shock tunnel design,improving experimental accuracy,and ensuring operational safety.To address the complex challenge of fully coupled multiphysics analysis in high-pressure-ratio shock tunnels,this study introduces a high-fidelity,three-dimensional,fully coupled Fluid-Structure Interaction(FSI)simulation framework.This framework seamlessly integrates the Dual Conservation Element and Solution Element(Dual-CESE)method,the Immersed Boundary Method(IBM),and the JohnsonCook(J-C)material constitutive and failure model.The combined approach enables synchronized simulation and analysis of the entire diaphragm rupture sequence—including pre-deformation,crack initiation and propagation,and fully developed petaling deformation—alongside the formation and evolution of the associated supersonic flow field.The simulation results show strong agreement with experimental observations,with the post-rupture geometric morphology accurately replicated and a shock wave velocity deviation of only 2.55%from experimental measurements.The study uncovers the dynamic failure mechanisms,revealing that nonlinear pressure loading initiates cracking within the diaphragm.It further elucidates how the nonlinearly coupled interactions between petaling dynamics and fracture morphology directly impact shock wave formation and evolution.This computational framework provides a novel and robust methodology for advancing shock tunnel design and conducting comprehensive reliability assessments.展开更多
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.展开更多
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.展开更多
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.展开更多
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%.展开更多
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.展开更多
In the study of thin-walled wedge shells subjected to lateral water entry impact loads,similarity theory is employed to predict the behavior of the prototype by analyzing the dynamic response of the model.However,the ...In the study of thin-walled wedge shells subjected to lateral water entry impact loads,similarity theory is employed to predict the behavior of the prototype by analyzing the dynamic response of the model.However,the traditional similarity law for structural impact fails to accurately describe the relationship between the dynamic responses of the model and the prototype,mainly due to the ignoring of the fluid-structure coupling effect and geometric distortion in the thickness direction.To address this limitation,this study employs dimensional analysis to derive the scaling factors for achieving structural similarity during water entry while considering the effects of fluid-structure coupling.Based on the analysis of the governing equations for thin-walled wedge shells under water entry impact loads,the scaling factor accounting for geometric distortion for water entry of wedge shells is determined.Numerical models are established to verify the applicability of the geometric distortion similarity model with different materials.The results demonstrate the effectiveness of the proposed similarity model,which significantly reduces the disparities in displacement peaks and energy between the scaled model and the prototype.In addition,experimental platforms are constructed to further verify the proposed similarity model by performing vertical water entry tests on thickness-distorted and material-distorted specimens.展开更多
This study addresses the optimization of automated yarn handling in textile manufacturing by examining the related suction process through a combined numerical and experimental approach.In particular,a three-dimension...This study addresses the optimization of automated yarn handling in textile manufacturing by examining the related suction process through a combined numerical and experimental approach.In particular,a three-dimensional model of the suction nozzle was coupled with an equivalent linear-elastic beam representation of the yarn,and a Fluent-IDW-Abaqus weakly coupled fluid-structure interaction(FSI)framework was employed to capture the yarn’s release and dynamic response under negative-pressure suction.High-speed imaging experiments validated the simulations,demonstrating excellent agreement in displacements and velocities.According to the results,increasing the initial suction pressure from -0.04 MPa to -0.06 MPa reduces adsorption time by approximately 62% and markedly dampens yarn-end vibrations,enhancing suction performance.Pressures beyond -0.06 MPa,however,induce overshoot and nozzle collisions,increasing the risk of entanglement and mechanical damage.The outcomes of a statistical analysis are also presented to further quantify the interplay among energy consumption,suction efficiency,and operational success under varying pressures,thereby providing a rigorous foundation for the optimal selection of pressure parameters in automated yarn-handling systems.展开更多
The intrinsic interaction mechanism of flutter between the flow and structure of a rectangular plate remains a mystery from the viewpoint of unsteady flow.The present study provides a novel insight into this interacti...The intrinsic interaction mechanism of flutter between the flow and structure of a rectangular plate remains a mystery from the viewpoint of unsteady flow.The present study provides a novel insight into this interaction mechanism based on an adequate understanding of the formation and evolution of the flapping leading-edge vortex(LEV).A series of wind tunnel tests was conducted to investigate the nonlinear flutter instability of an 8∶1 rectangular plate.The complete flow fields around the model throughout the flutter process were obtained by a particle image velocimetry(PIV)technique using two synchronous cameras with an interpolation and resampling method.To acquire the flow structures corresponding to the characteristic frequency of flutter,the spectral proper orthogonal decomposition(SPOD)method was extended to a noninertial frame to reconstruct the low-rank flow field during flutter and extract the characteristic flow pattern coupled with oscillations.It was found that when the 8∶1 rectangular plate undergoes flutter,the LEVs exhibit a periodic flapping phenomenon induced by the structure oscillations.A two-dimensional correlation analysis of the flapping LEVs was conducted for different inflow velocities.The results demonstrate that there is a substantial phase lead phenomenon in the LEV evolution downstream for a higher inflow velocity.This phenomenon may be related to a phase offset of aerodynamic forces,and finally,it gives rise to flutter.展开更多
This paper presents a simplified design tool based on semi-analytical formulations to investigate the dynamic response of an immersed composite cylinder subjected to a far-field underwater explosion.The cylinder is si...This paper presents a simplified design tool based on semi-analytical formulations to investigate the dynamic response of an immersed composite cylinder subjected to a far-field underwater explosion.The cylinder is simply supported,fully submerged and filled with air inside.A classical shell theory using a Double Fourier series solution combined with the first-order Doubly Asymptotic Approximation(DAA1)formulation is adapted to model the fluid-structure interaction.An explicit non-standard finite difference scheme is applied to solve the coupled differential equations in time domain.The validity of DAA1 model is established by comparing the LS-DYNA/USA finite element results with existing experimental data from the literature.Then the proposed semi-analytical solutions are compared to the LS-DYNA/USA results,showing good correlation with a discrepancy of 7%for peak deflections and±9%for maximum stresses at the stand-off point for cylinders with relatively small length over radius ratios.Parametric studies examining the effect of different loading conditions,areal masses,and material configurations reveal that a large charge mass located far from the composite panel turns out to be more damaging than a small mass located nearby due to a broader pressure-time profile.Finally,the proposed model demonstrates a significant reduction in computation time,being approximately 30 times faster than its numerical counterpart,LS-DYNA/USA,making it a valuable tool for the preliminary design stages.展开更多
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.展开更多
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 of...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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles a...This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles are analyzed synchronously using FLUENT software, and the effects of fluid-structure interaction on the performance of the labyrinth seal are revealed. The results indicate that with the growth of tooth profile angle, the leakage rate of labyrinth seal tends to decrease first, and then increase. The results of fluid-structure interaction analysis are close to those of actual engineering. The effect of fluid-structure interaction makes tiny deformation in calculation mesh of piston and cylinder structure, and the coupling interaction affects the performance of the labyrinth seal.展开更多
文摘The Donalda gold deposit in the southern part of the Archean Abitibi greenstone belt consists mainly of a subhonzontal gold-quartz vein perpendicular to subvertical shear zones.The 0.3—0.5 m thick vein is characterized by vein-parallel banding structures indicating multiple episodes of fracture opening and mineral precipitation.Measurement of the c-axis of primary growth quartz indicates that quartz preferentially grew perpendicular to the fracture,suggesting open space filling and/or extensional nature of the fracture.Measurement of the orientations of microfractures,veinlets and fluid—inclusion planes（FIPs） crosscutting primary growth quartz indicates that the vein minerals were subject to a vertical maximum principal stress（σ_1）,which is inconsistent with the subhorizontalσ_1 inferred from the regional stress field with N—S shortening.This apparent discrepancy is explained by invoking episodic fluid pressure fluctuation between supralithostatic and hydrostatic regimes accompanied by episodic opening and closing of the sub-horizontal fracture.When fluid pressure was higher than the lithostatic value,the fracture was opened and primary growth minerals were precipitated,whereas when fluid pressure decreased toward the hydrostatic value,the hanging wall of the fracture collapsed,causing collision of protruding primary growth minerals from both sides of the fracture and resulting in formation of vein-parallel deformation bands.The columns where the two facing sides of the fracture collided were subject to higher-than-lithostatic stress due to the bridging effect and reduced support surface area,explaining the development of verticalσ_1.This hypothesis is consistent the fault-valve model,and explains the flipping ofσ_1 without having to change the regional stress field.
文摘In order to enable wing morphing(e.g.change in camber or folds)without incurring additional weight to the aircraft,lightweight flexible materials such as membrane are needed.However,the research on fluid-structure coupling of membranes has mainly focused on parachutes in low-speed conditions,while that in supersonic flow conditions is lacking.Here,the degraded shell method is proposed to study membrane deformation by using shell element,which is more effective than using membrane elements directly.A fluid-structure interaction computational framework is proposed,whereby the aerodynamic module is composed of either the piston theory or computational fluid dynamics.A rectangular membrane of length 0.4 m and width 0.6 m is investigated in supersonic conditions.The characteristics of the limit cycle and steady deformation are analyzed,considering the effects of angle of attack and dynamic pressure.It is found that the structural response exhibits significant differences under various angles of attack.Furthermore,initial relaxation of membrane has significant influence on the structural deformation.Finally,the aeroelastic scaling method for membrane structures is derived,providing guidance for the design of wind tunnel models.This study provides a theoretical foundation for the analysis and application of membrane structures under supersonic conditions in future research.
基金supported by the National Key R&D Program of China(No.2021YFC3100700)。
文摘Shock tunnels are indispensable facilities for hypersonic aerodynamic experimentation.Within these systems,the diaphragm plays a pivotal role,as its rupture process critically influences shock wave generation quality,experimental repeatability,and facility reliability.A thorough understanding of diaphragm rupture dynamics is therefore essential for optimizing shock tunnel design,improving experimental accuracy,and ensuring operational safety.To address the complex challenge of fully coupled multiphysics analysis in high-pressure-ratio shock tunnels,this study introduces a high-fidelity,three-dimensional,fully coupled Fluid-Structure Interaction(FSI)simulation framework.This framework seamlessly integrates the Dual Conservation Element and Solution Element(Dual-CESE)method,the Immersed Boundary Method(IBM),and the JohnsonCook(J-C)material constitutive and failure model.The combined approach enables synchronized simulation and analysis of the entire diaphragm rupture sequence—including pre-deformation,crack initiation and propagation,and fully developed petaling deformation—alongside the formation and evolution of the associated supersonic flow field.The simulation results show strong agreement with experimental observations,with the post-rupture geometric morphology accurately replicated and a shock wave velocity deviation of only 2.55%from experimental measurements.The study uncovers the dynamic failure mechanisms,revealing that nonlinear pressure loading initiates cracking within the diaphragm.It further elucidates how the nonlinearly coupled interactions between petaling dynamics and fracture morphology directly impact shock wave formation and evolution.This computational framework provides a novel and robust methodology for advancing shock tunnel design and conducting comprehensive reliability assessments.
基金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.
基金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 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.
基金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%.
基金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 Natural Science Foundation of China(Grant Nos.11972309,12272320 and 12272309)the Overseas Expertise Introduction Project for Discipline Innovation(the 111 Project)(Grant No.BP0719007).
文摘In the study of thin-walled wedge shells subjected to lateral water entry impact loads,similarity theory is employed to predict the behavior of the prototype by analyzing the dynamic response of the model.However,the traditional similarity law for structural impact fails to accurately describe the relationship between the dynamic responses of the model and the prototype,mainly due to the ignoring of the fluid-structure coupling effect and geometric distortion in the thickness direction.To address this limitation,this study employs dimensional analysis to derive the scaling factors for achieving structural similarity during water entry while considering the effects of fluid-structure coupling.Based on the analysis of the governing equations for thin-walled wedge shells under water entry impact loads,the scaling factor accounting for geometric distortion for water entry of wedge shells is determined.Numerical models are established to verify the applicability of the geometric distortion similarity model with different materials.The results demonstrate the effectiveness of the proposed similarity model,which significantly reduces the disparities in displacement peaks and energy between the scaled model and the prototype.In addition,experimental platforms are constructed to further verify the proposed similarity model by performing vertical water entry tests on thickness-distorted and material-distorted specimens.
文摘This study addresses the optimization of automated yarn handling in textile manufacturing by examining the related suction process through a combined numerical and experimental approach.In particular,a three-dimensional model of the suction nozzle was coupled with an equivalent linear-elastic beam representation of the yarn,and a Fluent-IDW-Abaqus weakly coupled fluid-structure interaction(FSI)framework was employed to capture the yarn’s release and dynamic response under negative-pressure suction.High-speed imaging experiments validated the simulations,demonstrating excellent agreement in displacements and velocities.According to the results,increasing the initial suction pressure from -0.04 MPa to -0.06 MPa reduces adsorption time by approximately 62% and markedly dampens yarn-end vibrations,enhancing suction performance.Pressures beyond -0.06 MPa,however,induce overshoot and nozzle collisions,increasing the risk of entanglement and mechanical damage.The outcomes of a statistical analysis are also presented to further quantify the interplay among energy consumption,suction efficiency,and operational success under varying pressures,thereby providing a rigorous foundation for the optimal selection of pressure parameters in automated yarn-handling systems.
基金The National Natural Science Foundation of China(No.52178470)the National Key Research and Development Program of China(No.2022YFC3005303).
文摘The intrinsic interaction mechanism of flutter between the flow and structure of a rectangular plate remains a mystery from the viewpoint of unsteady flow.The present study provides a novel insight into this interaction mechanism based on an adequate understanding of the formation and evolution of the flapping leading-edge vortex(LEV).A series of wind tunnel tests was conducted to investigate the nonlinear flutter instability of an 8∶1 rectangular plate.The complete flow fields around the model throughout the flutter process were obtained by a particle image velocimetry(PIV)technique using two synchronous cameras with an interpolation and resampling method.To acquire the flow structures corresponding to the characteristic frequency of flutter,the spectral proper orthogonal decomposition(SPOD)method was extended to a noninertial frame to reconstruct the low-rank flow field during flutter and extract the characteristic flow pattern coupled with oscillations.It was found that when the 8∶1 rectangular plate undergoes flutter,the LEVs exhibit a periodic flapping phenomenon induced by the structure oscillations.A two-dimensional correlation analysis of the flapping LEVs was conducted for different inflow velocities.The results demonstrate that there is a substantial phase lead phenomenon in the LEV evolution downstream for a higher inflow velocity.This phenomenon may be related to a phase offset of aerodynamic forces,and finally,it gives rise to flutter.
基金supported by French Defense Innovation Agency(AID-DGA)(Grant No.ANR-21-ASM2-0002-02)in the framework of the Astrid Maturation SUCCESS+project,a collaborative French research project.
文摘This paper presents a simplified design tool based on semi-analytical formulations to investigate the dynamic response of an immersed composite cylinder subjected to a far-field underwater explosion.The cylinder is simply supported,fully submerged and filled with air inside.A classical shell theory using a Double Fourier series solution combined with the first-order Doubly Asymptotic Approximation(DAA1)formulation is adapted to model the fluid-structure interaction.An explicit non-standard finite difference scheme is applied to solve the coupled differential equations in time domain.The validity of DAA1 model is established by comparing the LS-DYNA/USA finite element results with existing experimental data from the literature.Then the proposed semi-analytical solutions are compared to the LS-DYNA/USA results,showing good correlation with a discrepancy of 7%for peak deflections and±9%for maximum stresses at the stand-off point for cylinders with relatively small length over radius ratios.Parametric studies examining the effect of different loading conditions,areal masses,and material configurations reveal that a large charge mass located far from the composite panel turns out to be more damaging than a small mass located nearby due to a broader pressure-time profile.Finally,the proposed model demonstrates a significant reduction in computation time,being approximately 30 times faster than its numerical counterpart,LS-DYNA/USA,making it a valuable tool for the preliminary design stages.
文摘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 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.
文摘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 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.
文摘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.
基金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 Science and Technology Projects of Liaoning Province(No.2012219020)the China Postdoctoral Science Foundation(No.2013M541249)
文摘This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles are analyzed synchronously using FLUENT software, and the effects of fluid-structure interaction on the performance of the labyrinth seal are revealed. The results indicate that with the growth of tooth profile angle, the leakage rate of labyrinth seal tends to decrease first, and then increase. The results of fluid-structure interaction analysis are close to those of actual engineering. The effect of fluid-structure interaction makes tiny deformation in calculation mesh of piston and cylinder structure, and the coupling interaction affects the performance of the labyrinth seal.
