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.展开更多
This study presents an implicit multiphysics coupling method integrating Computational Fluid Dynamics(CFD),the Multiphase Particle-in-Cell(MPPIC)model,and the Finite Element Method(FEM),implemented with OpenFOAM,Calcu...This study presents an implicit multiphysics coupling method integrating Computational Fluid Dynamics(CFD),the Multiphase Particle-in-Cell(MPPIC)model,and the Finite Element Method(FEM),implemented with OpenFOAM,CalculiX,and preCICE to simulate fluid-particle-structure interactions with large deformations.Mesh motion in the fluid field is handled using the radial basis function(RBF)method.The particle phase is modeled by MPPIC,where fluid-particle interaction is described through momentum exchange,and inter-particle collisions are characterized by collision stress.The structural field is solved by nonlinear FEM to capture large deformations induced by geometric nonlinearity.Coupling among fields is realized through a partitioned,parallel,and non-intrusive iterative strategy,ensuring stable transfer and convergence of interface forces and displacements.Notably,the influence of particles on the structure is not direct but mediated by the fluid,while structural motion directly affects particle dynamics.The results demonstrate that the proposed approach effectively captures multiphysics interaction processes and provides a valuable reference for numerical modeling of coupled fluid-particle-structure systems.展开更多
This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
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%.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between ...The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.展开更多
Spaceborne antennas are essential for remote sensing,deep-space communication,and Earth observation,yet their trajectory planning is complicated by nonlinear base-manipulator coupling and antenna flexibility.To addres...Spaceborne antennas are essential for remote sensing,deep-space communication,and Earth observation,yet their trajectory planning is complicated by nonlinear base-manipulator coupling and antenna flexibility.To address these challenges,this paper proposes a multi-objective trajectory optimization framework.The system dynamics capture both nonlinear rigid-flexible coupling and antenna deformation through a reduced-order formulation.To enhance discretization efficiency,a predictive-terminal hp-adaptive pseudospectral method is employed,assigning collocation density based on task-phase characteristics:finer resolution is applied to dynamic segments requiring higher accuracy,especially near the terminal phase.This enables efficient transcription of the continuous-time problem into a Nonlinear Programming Problem(NLP).The resulting NLP is then solved using a multi-objective optimization strategy based on the nondominated sorting genetic algorithm II,which explores trade-offs among antenna pointing accuracy,energy consumption,and structural vibration.Numerical results demonstrate that the proposed method achieves a reduction of approximately 14.0% in control energy and 41.8%in peak actuation compared to a GPOPS-II baseline,while significantly enhancing vibration suppression.The resulting Pareto front reveals structured trade-offs and clustered solutions,offering robust and diverse options for precision,low-disturbance mission planning.展开更多
Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture...Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.展开更多
Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is sti...Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.展开更多
The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are ca...The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are calculated. Secondly, the influences to runner frequency domain by large flow, small flow and design flow working conditions are compared. Finally the influences to runner modes by centrifugal forces under three rotating speeds of 400 r/rain, 500 r/min and 600 r/rain are compared. The centrifugal force and small flow working condition have greatly influence on the vibration of small runner. With the increase of centrifugal force, the vibration frequency of the runner is sharply increased. Some order frequencies are even close to the runner natural frequency in the air. Because the low frequency vibration will severely damage the stability of the turbine, low frequency vibration of units should be avoided as soon as possible.展开更多
The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more ...The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.展开更多
Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its co...Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct ...Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct transverse seismic response change. The transverse seismic response of a large-scale aqueduct in several work condition are calculated. It shows that the transverse seismic response is greatly influenced by the water mass in the aqueduct, but the shaking water play a TLD role. ff the whole water is appended aqueduct body, it will magnify seismic inertia action. When aqueduct cross-section is selected, the influence of ratio of depth and width to pier seismic response should be considered in order to reduce seismic action.展开更多
A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction pro...A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.展开更多
A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the key...A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.展开更多
The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. ...The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.展开更多
Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenome...Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.展开更多
A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed cou...A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.展开更多
An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted F...An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.展开更多
基金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.
基金supported in part by the Mining Hydraulic Technology and Equipment Engineering Research Center,Liaoning Technical University,Fuxin,China(Grant No.MHTE23-R04)the Fundamental Research Funds for the Central Universities(ID N25BSS068).
文摘This study presents an implicit multiphysics coupling method integrating Computational Fluid Dynamics(CFD),the Multiphase Particle-in-Cell(MPPIC)model,and the Finite Element Method(FEM),implemented with OpenFOAM,CalculiX,and preCICE to simulate fluid-particle-structure interactions with large deformations.Mesh motion in the fluid field is handled using the radial basis function(RBF)method.The particle phase is modeled by MPPIC,where fluid-particle interaction is described through momentum exchange,and inter-particle collisions are characterized by collision stress.The structural field is solved by nonlinear FEM to capture large deformations induced by geometric nonlinearity.Coupling among fields is realized through a partitioned,parallel,and non-intrusive iterative strategy,ensuring stable transfer and convergence of interface forces and displacements.Notably,the influence of particles on the structure is not direct but mediated by the fluid,while structural motion directly affects particle dynamics.The results demonstrate that the proposed approach effectively captures multiphysics interaction processes and provides a valuable reference for numerical modeling of coupled fluid-particle-structure systems.
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
文摘This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
基金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%.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
文摘The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.
基金supported by the National Natural Science Foundation of China(No.62173107).
文摘Spaceborne antennas are essential for remote sensing,deep-space communication,and Earth observation,yet their trajectory planning is complicated by nonlinear base-manipulator coupling and antenna flexibility.To address these challenges,this paper proposes a multi-objective trajectory optimization framework.The system dynamics capture both nonlinear rigid-flexible coupling and antenna deformation through a reduced-order formulation.To enhance discretization efficiency,a predictive-terminal hp-adaptive pseudospectral method is employed,assigning collocation density based on task-phase characteristics:finer resolution is applied to dynamic segments requiring higher accuracy,especially near the terminal phase.This enables efficient transcription of the continuous-time problem into a Nonlinear Programming Problem(NLP).The resulting NLP is then solved using a multi-objective optimization strategy based on the nondominated sorting genetic algorithm II,which explores trade-offs among antenna pointing accuracy,energy consumption,and structural vibration.Numerical results demonstrate that the proposed method achieves a reduction of approximately 14.0% in control energy and 41.8%in peak actuation compared to a GPOPS-II baseline,while significantly enhancing vibration suppression.The resulting Pareto front reveals structured trade-offs and clustered solutions,offering robust and diverse options for precision,low-disturbance mission planning.
基金supported by the Japan Society for the Promotion of Science(JSPS)KAKENHI(Grant Nos.JP23KK0182,JP23K26356,and JP24K00971).
文摘Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金Outstanding Youth Science Fund Subsidization of Sichuan Province, China (No. 05204033).
文摘The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are calculated. Secondly, the influences to runner frequency domain by large flow, small flow and design flow working conditions are compared. Finally the influences to runner modes by centrifugal forces under three rotating speeds of 400 r/rain, 500 r/min and 600 r/rain are compared. The centrifugal force and small flow working condition have greatly influence on the vibration of small runner. With the increase of centrifugal force, the vibration frequency of the runner is sharply increased. Some order frequencies are even close to the runner natural frequency in the air. Because the low frequency vibration will severely damage the stability of the turbine, low frequency vibration of units should be avoided as soon as possible.
基金Project supported by the National Natural Science Foundation of China(No.11972112)the Fundamental Research Funds for the Central Universities of China(Nos.N2103024 and N2103002)the Major Projects of Aero-Engines and Gasturbines(No.J2019-I-0008-0008)。
文摘The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.
基金funded by the BeijingNatural Science Foundation of China(8222003)National Natural Science Foundation of China(41807180).
文摘Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
基金National Natural Science Foundation of China (50279024)Lanzhou Jiaotong University Qinglan Talented Person Fund Project (QL-05-13A).
文摘Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct transverse seismic response change. The transverse seismic response of a large-scale aqueduct in several work condition are calculated. It shows that the transverse seismic response is greatly influenced by the water mass in the aqueduct, but the shaking water play a TLD role. ff the whole water is appended aqueduct body, it will magnify seismic inertia action. When aqueduct cross-section is selected, the influence of ratio of depth and width to pier seismic response should be considered in order to reduce seismic action.
基金open foundation of the Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanicsthe Open Foundation of Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment.
文摘A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.
基金Projects (51175162, 50805045) supported by the National Natural Science Foundation of ChinaProject supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education,China
文摘A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.
基金Sponsored by National Science and Technology Support Plan of China (2009AA04Z143)Science and Technology Support Plan of Hebei Province of China (10212101D)Important Natural Science Foundation of Hebei Province of China (E2006001038)
文摘The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.
文摘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.
基金Projects(51474251,51874351) supported by the National Natural Science Foundation of China
文摘A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.
基金Project (No. 50575143) supported by the National Natural Science Foundation of China
文摘An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.
