As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises s...As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.展开更多
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.展开更多
Understanding the deterioration behaviors and mechanisms of rocks under thermo-hydromechanical(THM)interactions is crucial for mitigating slope instability.In this study,the physicomechanical properties of silty mudst...Understanding the deterioration behaviors and mechanisms of rocks under thermo-hydromechanical(THM)interactions is crucial for mitigating slope instability.In this study,the physicomechanical properties of silty mudstone subjected to THM interactions were investigated by triaxial tests.The underlying micro-mechanisms were revealed using microscopic tests.The triaxial test results indicate that the strength parameters of silty mudstone decrease by 89.50%(deformation modulus),78.15%(peak strength),70.58%(cohesion),and 48.65%(friction angle)under 16 THM cycles,a load of 300 kPa,and alternating between 0℃water immersion and 60℃drying.The SEM test results indicate that the deterioration of silty mudstone strength primarily results from hydrothermal-expansion softening and cracking driven by the TLHM interactions.The specimens manifest shear failure under confining pressure exceeding 140 kPa.Furthermore,a new constitutive model considering hydrothermalexpansion strain and non-linear deformation characteristics was developed.The discrepancy between the experimentally measured peak strength and the damage constitutive model prediction remains below 5%.The proposed model is verified to be in satisfactory agreement with the experimental results.The self-designed THM apparatus overcomes the limitations of traditional investigations,enabling simultaneous consideration of thermal,hydraulic,and mechanical interactions.展开更多
Recycling of waste rubber(WR)is crucial for the sustainable development of the rubber industry.The enhancement of interfacial interactions is the main strategy for waste polymer recycling.However,there is a lack of me...Recycling of waste rubber(WR)is crucial for the sustainable development of the rubber industry.The enhancement of interfacial interactions is the main strategy for waste polymer recycling.However,there is a lack of methods for enhancing the interfacial interactions for WR recycling because WR contains abundant inert C―H bonds.Herein,we designed thioctic acid inverse vulcanization copolymers to endow recycled WR with dynamic disulfide interfacial interactions,significantly improving the mechanical properties of recycled WR.These disulfide interfacial interactions among the recycled WR tend to exchange,which dramatically increases the fractocohesive length and prevents stress concentration near the crack tips.When recycled WR is subjected to external stress,the loads are redistributed across a broad region of adjacent regions instead of being concentrated on a limited length scale,which resists crack propagation.This work effectively recycled WR,providing a strategy for solvent-free reaction-derived inverse vulcanization copolymers to improve the toughness of WR recycling.展开更多
For hypersonic air-breathing vehicles,the V-shaped leading edges(VSLEs)of supersonic combustion ramjet(scramjet)inlets experience complex shock interactions and intense aerodynamic loads.This paper provides a comprehe...For hypersonic air-breathing vehicles,the V-shaped leading edges(VSLEs)of supersonic combustion ramjet(scramjet)inlets experience complex shock interactions and intense aerodynamic loads.This paper provides a comprehensive review of flow characteristics at the crotch of VSLEs,with particular focus on the transition of shock interaction types and the variation of wall heat flux under different freestream Mach numbers and geometric configurations.The mechanisms governing shock transition,unsteady oscillations,hysteresis,and three-dimensional effects in VSLE flows are first examined.Subsequently,thermal protection strategies aimed at mitigating extreme heating loads are reviewed,emphasizing their relevance to practical engineering applications.Special attention is given to recent studies addressing thermochemical nonequilibrium effects on VSLE shock interactions,and the limitations of current research are critically assessed.Finally,perspectives for future investigations into hypersonic VSLE shock interactions are outlined,highlighting opportunities for advancing design and thermal management strategies.展开更多
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.展开更多
The coupled chemo-mechanical impact of supercritical CO_(2)-H_(2)O(ScCO_(2)-H_(2)O)reactions on fracture geometry and nonlinear flow regimes in deep shale under confining pressures remains inadequately quantified.This...The coupled chemo-mechanical impact of supercritical CO_(2)-H_(2)O(ScCO_(2)-H_(2)O)reactions on fracture geometry and nonlinear flow regimes in deep shale under confining pressures remains inadequately quantified.This study systematically investigates the effects of ScCO_(2)-H_(2)O-shale interactions on fracture morphology and flow properties under confining pressures from 15 MPa to 40 MPa by integrating XRD(X-ray diffraction),micro-CT,3D surface profilometry,and multistage steady-state flow experiments.The results demonstrate that ScCO_(2)-H_(2)O exposure drives pyrite/feldspar dissolution and localized clay precipitation,resulting in fracture branching and macroscopic aperture regularization.Critically,confining pressure dictates the net hydraulic response:under low confining pressure(15-25 MPa),dissolution dominates,enhancing permeability,flow efficiency(Q/VP),and pre-linear flow behavior(n<1).At high confining pressures(30-40 MPa)mechanical compaction and mineral precipitation amplify flow resistance,shifting the flow regime toward quasi-linear behavior,as inertial effects become negligible compared to dominant viscous forces and increased flow resistance.Confining pressure thus critically mediates the dissolution-precipitation balance during ScCO_(2)-H_(2)O treatment,with an optimal window of 15-25 MPa identified for enhancing conductivity while minimizing clogging risk.These findings provide a quantitative framework for predicting stress-dependent flow evolution in chemically altered shale fractures.展开更多
Effectively controlling the deformation and temperature of heated structures is crucial for achieving highperformance active cooling through fluid flow.In this study,the topology optimization design of structures cons...Effectively controlling the deformation and temperature of heated structures is crucial for achieving highperformance active cooling through fluid flow.In this study,the topology optimization design of structures considering fluid–structure interactions and heat transfer performance was investigated,and then optimized designs of two-dimensional/three-dimensional cooling impingement systems obtained using the proposed method were obtained.In the optimization model,the objective function was constructed as a weighted combination of the mechanical deformations at specific locations and the average temperature within the designated solid channel structures.Additionally,explicit functional interpolation models were introduced to establish connections between the thermal,fluid,and solid properties,along with the element densities.In the analysis model,the strongly coupled structural mechanical deformation and fluid velocity field were analyzed via a dynamicgrid-based finite element model with a Winslow elliptic smoother to automatically track the fluid–structure interface during the process of optimization.To solve the optimization problems,the globally convergent moving asymptotic optimizer method was used to adjust the design variables on the basis of the sensitivity analysis.A demonstration of the efficacy of the proposed algorithm is provided through the presentation of several optimization examples.Furthermore,two-and three-dimensional cooling impingement systems were designed with the proposed method.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China(No.12372233)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University,China(No.25GH01020005)the“111 Project”of China(No.B17037)。
文摘As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.
文摘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 Natural Science Foundation of China(52378440,52078067,52078066,42477143,52408458)the Key Science and Technology Program in the Transportation Industry(2022-MS1-032,2022-MS5-125)+4 种基金the Natural Science Foundation of Hunan Province(2023JJ10045)the Outstanding Innovative Youth Training Program of Changsha City(kq2305023)Scientific Research Foundation of Hunan Provincial Education Department(24B0292)Water Resources Science and Technology Project of Hunan Province(XSKJ2023059-41)the Guangxi Key Research and Development Program(AB23075184)。
文摘Understanding the deterioration behaviors and mechanisms of rocks under thermo-hydromechanical(THM)interactions is crucial for mitigating slope instability.In this study,the physicomechanical properties of silty mudstone subjected to THM interactions were investigated by triaxial tests.The underlying micro-mechanisms were revealed using microscopic tests.The triaxial test results indicate that the strength parameters of silty mudstone decrease by 89.50%(deformation modulus),78.15%(peak strength),70.58%(cohesion),and 48.65%(friction angle)under 16 THM cycles,a load of 300 kPa,and alternating between 0℃water immersion and 60℃drying.The SEM test results indicate that the deterioration of silty mudstone strength primarily results from hydrothermal-expansion softening and cracking driven by the TLHM interactions.The specimens manifest shear failure under confining pressure exceeding 140 kPa.Furthermore,a new constitutive model considering hydrothermalexpansion strain and non-linear deformation characteristics was developed.The discrepancy between the experimentally measured peak strength and the damage constitutive model prediction remains below 5%.The proposed model is verified to be in satisfactory agreement with the experimental results.The self-designed THM apparatus overcomes the limitations of traditional investigations,enabling simultaneous consideration of thermal,hydraulic,and mechanical interactions.
基金financially supported by the National Natural Science Foundation of China(No.52363007)。
文摘Recycling of waste rubber(WR)is crucial for the sustainable development of the rubber industry.The enhancement of interfacial interactions is the main strategy for waste polymer recycling.However,there is a lack of methods for enhancing the interfacial interactions for WR recycling because WR contains abundant inert C―H bonds.Herein,we designed thioctic acid inverse vulcanization copolymers to endow recycled WR with dynamic disulfide interfacial interactions,significantly improving the mechanical properties of recycled WR.These disulfide interfacial interactions among the recycled WR tend to exchange,which dramatically increases the fractocohesive length and prevents stress concentration near the crack tips.When recycled WR is subjected to external stress,the loads are redistributed across a broad region of adjacent regions instead of being concentrated on a limited length scale,which resists crack propagation.This work effectively recycled WR,providing a strategy for solvent-free reaction-derived inverse vulcanization copolymers to improve the toughness of WR recycling.
基金funded by the Research Fund of National Key Laboratory of Aerospace Physics in Fluids,grant number 2024-APF-KFZD-01Guangdong Basic and Applied Basic Research Foundation,grant number 2025A1515012081+1 种基金National Natural Science Foundation of China,grant number 12002193Shandong Provincial Natural Science Foundation,China,grant number ZR2019QA018.
文摘For hypersonic air-breathing vehicles,the V-shaped leading edges(VSLEs)of supersonic combustion ramjet(scramjet)inlets experience complex shock interactions and intense aerodynamic loads.This paper provides a comprehensive review of flow characteristics at the crotch of VSLEs,with particular focus on the transition of shock interaction types and the variation of wall heat flux under different freestream Mach numbers and geometric configurations.The mechanisms governing shock transition,unsteady oscillations,hysteresis,and three-dimensional effects in VSLE flows are first examined.Subsequently,thermal protection strategies aimed at mitigating extreme heating loads are reviewed,emphasizing their relevance to practical engineering applications.Special attention is given to recent studies addressing thermochemical nonequilibrium effects on VSLE shock interactions,and the limitations of current research are critically assessed.Finally,perspectives for future investigations into hypersonic VSLE shock interactions are outlined,highlighting opportunities for advancing design and thermal management strategies.
基金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.
基金support from the Science and Technology Innovation Program of Hunan Province(Grant No.2023RC1021)the Natural Science Foundation of Sichuan Province(Grant No.2025YFHZ0323).-。
文摘The coupled chemo-mechanical impact of supercritical CO_(2)-H_(2)O(ScCO_(2)-H_(2)O)reactions on fracture geometry and nonlinear flow regimes in deep shale under confining pressures remains inadequately quantified.This study systematically investigates the effects of ScCO_(2)-H_(2)O-shale interactions on fracture morphology and flow properties under confining pressures from 15 MPa to 40 MPa by integrating XRD(X-ray diffraction),micro-CT,3D surface profilometry,and multistage steady-state flow experiments.The results demonstrate that ScCO_(2)-H_(2)O exposure drives pyrite/feldspar dissolution and localized clay precipitation,resulting in fracture branching and macroscopic aperture regularization.Critically,confining pressure dictates the net hydraulic response:under low confining pressure(15-25 MPa),dissolution dominates,enhancing permeability,flow efficiency(Q/VP),and pre-linear flow behavior(n<1).At high confining pressures(30-40 MPa)mechanical compaction and mineral precipitation amplify flow resistance,shifting the flow regime toward quasi-linear behavior,as inertial effects become negligible compared to dominant viscous forces and increased flow resistance.Confining pressure thus critically mediates the dissolution-precipitation balance during ScCO_(2)-H_(2)O treatment,with an optimal window of 15-25 MPa identified for enhancing conductivity while minimizing clogging risk.These findings provide a quantitative framework for predicting stress-dependent flow evolution in chemically altered shale fractures.
文摘Effectively controlling the deformation and temperature of heated structures is crucial for achieving highperformance active cooling through fluid flow.In this study,the topology optimization design of structures considering fluid–structure interactions and heat transfer performance was investigated,and then optimized designs of two-dimensional/three-dimensional cooling impingement systems obtained using the proposed method were obtained.In the optimization model,the objective function was constructed as a weighted combination of the mechanical deformations at specific locations and the average temperature within the designated solid channel structures.Additionally,explicit functional interpolation models were introduced to establish connections between the thermal,fluid,and solid properties,along with the element densities.In the analysis model,the strongly coupled structural mechanical deformation and fluid velocity field were analyzed via a dynamicgrid-based finite element model with a Winslow elliptic smoother to automatically track the fluid–structure interface during the process of optimization.To solve the optimization problems,the globally convergent moving asymptotic optimizer method was used to adjust the design variables on the basis of the sensitivity analysis.A demonstration of the efficacy of the proposed algorithm is provided through the presentation of several optimization examples.Furthermore,two-and three-dimensional cooling impingement systems were designed with the proposed method.
基金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.
