The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more ...The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.展开更多
The 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive def...Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive deformation characteristics.This study explores the interaction mechanism between micro-scale variable-sectional shearing flow with hyper-rotation speeds and a three-layer elastic foil assembly through bidirectional aero-elastic coupling in a Multi-layer Thrust Gas Foil Bearing(MTGFB).The bearing capacity of the MTGFB varies non-linearly with the decrease of gas film clearance,while the collaborative deformation of the three-layer elastic foil assembly can deal with different load conditions.As the load capacity increases,the enhanced dynamic pressure effect causes the top foil to evolve from a single arch to multiple arches.The hydrodynamic effects in the gas film evolve to form multiple segmented wedges with different pitch ratios,while the peak pressure of the gas film always occurs near the vaults of the top foil.As the rotational speed frequency approaches the natural frequency,the resonance of the gas film and elastic foil assembly system occurs,and a phase delay occurs between the pressure pulsation and the vibration of foils.The load capacity of the MTGFB also depends on the elastic moduli of the elastic foil assembly.Increasing the elastic modulus decreases the deformation amplitude of the top foil,whereas it increases those of the backboard and middle foil,increasing the load capacity.展开更多
Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these unde...Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these underwater structures. The vibration characteristics of underwater structures are influenced not only by hydrodynamic pressure but also by hydrostatic pressure corresponding to different water depths. In this study, an acoustic finite element method was used to evaluate the underwater structures. Taken the hydrostatic pressure into account in terms of initial stress stiffness, an acoustical fluid-structure coupled analysis of underwater cylindrical shells has been made to study the effect of hydrodynamic pressures on natural frequency and sound radiation. By comparing with the frequencies obtained by the acoustic finite element method and by the added mass method based on the Bessel function, the validity of present analysis was checked. Finally, test samples of the sound radiation of stiffened cylindrical shells were acquired by a harmonic acoustic analysis. The results showed that hydrostatic pressure plays an important role in determining a large submerged body motion, and the characteristics of sound radiation change with water depth. Furthermore, the analysis methods and the results are of significant reference value for studies of other complicated submarine structures.展开更多
We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive ...We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.展开更多
The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computatio...The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.展开更多
The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift...The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute twophase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersedphase flow fields （including the number density of particles） in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.展开更多
The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STA...The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STAR-CCM+software.First,the external characteristic curve of the engine is compared with a one-dimen-sional simulation model,then the parameters of the model are modified until the curve matches the available experimental values.GT-POWER is then used to transfer the inlet boundary data under transient conditions to STAR-CCM+in real-time.The temperature profiles of the inner and outer walls of the exhaust manifold are obtained in this way,together with the thermal stress and thermal deformation of the exhaust manifold itself.Using this information,the original model is improved through the addition of connections.Moreover,the local branch pipes are optimized,leading to significant improvements in terms of thermal stress and thermal deforma-tion of the exhaust manifold(a 7%reduction in the maximum thermal stress).展开更多
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.展开更多
Clearance-fit(side-fit)spline joints are a key component in a permanent magnet synchronous motor(PMSM)in electric submersible pumping wells.The nonuniform spline clearance affects the output performance of a PMSM.A co...Clearance-fit(side-fit)spline joints are a key component in a permanent magnet synchronous motor(PMSM)in electric submersible pumping wells.The nonuniform spline clearance affects the output performance of a PMSM.A concept for energy conservation is optimized in this study to improve the modeling accuracy of electromagnetic torque.However,most existing computation models are one-way model with a magneto-mechanical simulation.In this study,a more accurate two-way coupling method is presented for simulating the electromagnetic and mechanical characteristics of PMSM.Additionally,importance should be attached to this two-way magneto-mechanical coupling methodology in an actual simulation.The coupled power,electrical and magnetic energy,and electromagnetic torque equations are solved iteratively until convergence for PMSMs with a segmented rotor and a non-segmented rotor.The optimal electromagnetic torque is obtained for different rotor configurations with the change of temperatures and rotational speeds.The results show that the output performance and electromagnetic torque of the PMSM are seriously affected by the effects of two-way magneto-mechanical coupling and nonuniform spline clearance.The proposed two-way coupling model gives more reasonable predictions than other one-way models do,because the power transfer between the electrical and magnetic energy can be modeled more accurately.The self-centralizing performance of clearance-fit splines and the sensitivity to the radial clearance magnitude lead to the reduction of the electromagnetic torque for the PMSM.Additionally,the electromagnetic torques decrease with the enhanced rotor temperatures and rotational speeds.The best rotor temperature and rotational speed are chosen through a comparison of the experimental results,and then the optimal electromagnetic torque is provided to ensure the output performance of the PMSM in electric submersible pumping wells.展开更多
The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the...The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the computation by combining computational fluid dynamics （CFD） and FEA which can simulate the TEHD more accurately. In this paper, by using both direct and separate coupled solutions together, steady TEHD lubrication considering the viscosity-temperature effect for a bidirectional thrust bearing in a pump-turbine unit is simulated combining a 3D CFD model for the oil film with a 3D FEA model for the pad and mirror plate. Cyclic symmetry condition is used in the oil film flow as more reasonable boundary conditions which avoids the oil temperature assumption at the leading and trailing edge. Deformations of the pad and mirror plate are predicted and discussed as well as the distributions of oil film thickness, pressure, temperature. The predicted temperature shows good agreement with measurements, while the pressure shows a reasonable distribution comparing with previous studies. Further analysis of the three-coupled-field reveals the reason of the high pressure and high temperature generated in the film. Finally, the influence of rotational speed of the mirror plate on the lubrication characteristics is illustrated which shows the thrust load should be balanced against the oil film temperature and pressure in optimized designs. This research proposes a thrust bearing computation method by combining CFD and FEA which can do the TEHD analysis more accurately.展开更多
Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenome...Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.展开更多
A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in t...A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.展开更多
The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional can...The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.展开更多
Failure pressure is a key parameter in reservoir hydrofracturing operation. Existing analytical methods for calculating the failure pressure are based on the assumption that borehole fluid is under two extreme conditi...Failure pressure is a key parameter in reservoir hydrofracturing operation. Existing analytical methods for calculating the failure pressure are based on the assumption that borehole fluid is under two extreme conditions: non-infiltration or complete infiltration. The assumption is not suitable for the actual infiltration process, and this will cause a great error in practical calculation. It shows that during the injection process, the dynamic variation in effective stress-dependent permeability has an influence on the infiltration, and the influence also brings about calculation errors. Based on the fluid-structure interaction and finite element method (FEM), considering partial infiltration during injection process, a numerical model for calculating rock failure pressure is established. According to the analysis of permeability test results and response-surface method, a new variation rule of rock permeability with the change of effective stress is presented, and the relationships among the permeability, confining pressure and pore pressure are proposed. There are some differences between the dynamic value of permeability-effective-stress coefficient observed herein and the one obtained by the classical theory. Combining with the numerical model and the dynamic permeability, a coupling method for calculating failure pressure is developed. Comparison of field data and calculated values obtained by various methods shows that accurate values can be obtained by the coupling method. The coupling method can be widely applied to the calculation of failure pressure of reservoirs and complex wells to achieve effective fracturing operation.展开更多
As the power transmission system of an aircraft,a hydraulic pipeline system is equivalent to the " blood vessel" of the aircraft. With the development of aircraft hydraulic system to high pressure,high speed...As the power transmission system of an aircraft,a hydraulic pipeline system is equivalent to the " blood vessel" of the aircraft. With the development of aircraft hydraulic system to high pressure,high speed and high power ratio,the fluid-structure interaction vibration mechanism of hydraulic pipeline is more complex and the influence of friction coupling on vibration cannot be ignored. The fluid-structure interaction of hydraulic pipeline will lead to system vibration,lower reliability of system operation and even pipeline rupture. Taking a hydraulic pipeline of C919 aircraft wingtip as the research object,a 14-equation model of fluid-structure interaction vibration considering friction coupling effect is established in this paper. The effects of friction and fluid parameters on the pipeline fluid-structure interaction vibration characteristics are studied and verified by experiments. The research results will provide theoretical guidance for the analysis of the pipeline fluid-structure interaction vibration and have important theoretical significance and great engineering value for promoting the localization process of large aircraft.展开更多
基金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.
基金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.
基金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.
文摘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.
基金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.
文摘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.
基金the financial support from the National Natural Science Foundation of China(No.52206091)the Aeronautical Science Foundation of China(No.201928052008)the Natural Science Foundation of Jiangsu Province,China(No.BK20210303)。
文摘Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive deformation characteristics.This study explores the interaction mechanism between micro-scale variable-sectional shearing flow with hyper-rotation speeds and a three-layer elastic foil assembly through bidirectional aero-elastic coupling in a Multi-layer Thrust Gas Foil Bearing(MTGFB).The bearing capacity of the MTGFB varies non-linearly with the decrease of gas film clearance,while the collaborative deformation of the three-layer elastic foil assembly can deal with different load conditions.As the load capacity increases,the enhanced dynamic pressure effect causes the top foil to evolve from a single arch to multiple arches.The hydrodynamic effects in the gas film evolve to form multiple segmented wedges with different pitch ratios,while the peak pressure of the gas film always occurs near the vaults of the top foil.As the rotational speed frequency approaches the natural frequency,the resonance of the gas film and elastic foil assembly system occurs,and a phase delay occurs between the pressure pulsation and the vibration of foils.The load capacity of the MTGFB also depends on the elastic moduli of the elastic foil assembly.Increasing the elastic modulus decreases the deformation amplitude of the top foil,whereas it increases those of the backboard and middle foil,increasing the load capacity.
基金China National 111 Project Under Grant No. B07019.
文摘Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these underwater structures. The vibration characteristics of underwater structures are influenced not only by hydrodynamic pressure but also by hydrostatic pressure corresponding to different water depths. In this study, an acoustic finite element method was used to evaluate the underwater structures. Taken the hydrostatic pressure into account in terms of initial stress stiffness, an acoustical fluid-structure coupled analysis of underwater cylindrical shells has been made to study the effect of hydrodynamic pressures on natural frequency and sound radiation. By comparing with the frequencies obtained by the acoustic finite element method and by the added mass method based on the Bessel function, the validity of present analysis was checked. Finally, test samples of the sound radiation of stiffened cylindrical shells were acquired by a harmonic acoustic analysis. The results showed that hydrostatic pressure plays an important role in determining a large submerged body motion, and the characteristics of sound radiation change with water depth. Furthermore, the analysis methods and the results are of significant reference value for studies of other complicated submarine structures.
基金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 Shanghai 2021“Science and Technology Innovation Action Plan”:Social Development Science and Technology Research Project(Grant No.21DZ1202703).
文摘The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.
基金The project supported by the National Natural Science Foundation of China(90205024)Russian Foundation for Basic Research(RFBR and(RFBR-NSFC-39004)
文摘The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute twophase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersedphase flow fields （including the number density of particles） in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.
基金supported by the Basic Ability Improvement Project for Young and Middle-Aged Teachers in Guangxi Universities,Project No.2021KY0792.
文摘The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STAR-CCM+software.First,the external characteristic curve of the engine is compared with a one-dimen-sional simulation model,then the parameters of the model are modified until the curve matches the available experimental values.GT-POWER is then used to transfer the inlet boundary data under transient conditions to STAR-CCM+in real-time.The temperature profiles of the inner and outer walls of the exhaust manifold are obtained in this way,together with the thermal stress and thermal deformation of the exhaust manifold itself.Using this information,the original model is improved through the addition of connections.Moreover,the local branch pipes are optimized,leading to significant improvements in terms of thermal stress and thermal deforma-tion of the exhaust manifold(a 7%reduction in the maximum thermal stress).
基金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.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52074161 and 52005281)Taishan Scholar Project of Shandong Province(Grant No.tsqn202211177)+2 种基金Shandong Provincial Plan for Introduction and Cultivation of Young Pioneers in Colleges and Universities(Grant No.2021-Qing Chuang-30613019)Natural Science Foundation of Shandong Province(Grant Nos.ZR2022ME173 and ZR2023QE011)Projects of CNOOC Research Institute Ltd.(Grant Nos.CCL2023RCPS0237RSN and CCL2023RCPS0319RSN)。
文摘Clearance-fit(side-fit)spline joints are a key component in a permanent magnet synchronous motor(PMSM)in electric submersible pumping wells.The nonuniform spline clearance affects the output performance of a PMSM.A concept for energy conservation is optimized in this study to improve the modeling accuracy of electromagnetic torque.However,most existing computation models are one-way model with a magneto-mechanical simulation.In this study,a more accurate two-way coupling method is presented for simulating the electromagnetic and mechanical characteristics of PMSM.Additionally,importance should be attached to this two-way magneto-mechanical coupling methodology in an actual simulation.The coupled power,electrical and magnetic energy,and electromagnetic torque equations are solved iteratively until convergence for PMSMs with a segmented rotor and a non-segmented rotor.The optimal electromagnetic torque is obtained for different rotor configurations with the change of temperatures and rotational speeds.The results show that the output performance and electromagnetic torque of the PMSM are seriously affected by the effects of two-way magneto-mechanical coupling and nonuniform spline clearance.The proposed two-way coupling model gives more reasonable predictions than other one-way models do,because the power transfer between the electrical and magnetic energy can be modeled more accurately.The self-centralizing performance of clearance-fit splines and the sensitivity to the radial clearance magnitude lead to the reduction of the electromagnetic torque for the PMSM.Additionally,the electromagnetic torques decrease with the enhanced rotor temperatures and rotational speeds.The best rotor temperature and rotational speed are chosen through a comparison of the experimental results,and then the optimal electromagnetic torque is provided to ensure the output performance of the PMSM in electric submersible pumping wells.
基金Supported by National Natural Science Foundation of China(Grant No.51439002)Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant Nos.20120002110011,20130002110072)Special Funds for Marine Renewable Energy Projects(Grant no.GHME2012GC02)
文摘The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the computation by combining computational fluid dynamics （CFD） and FEA which can simulate the TEHD more accurately. In this paper, by using both direct and separate coupled solutions together, steady TEHD lubrication considering the viscosity-temperature effect for a bidirectional thrust bearing in a pump-turbine unit is simulated combining a 3D CFD model for the oil film with a 3D FEA model for the pad and mirror plate. Cyclic symmetry condition is used in the oil film flow as more reasonable boundary conditions which avoids the oil temperature assumption at the leading and trailing edge. Deformations of the pad and mirror plate are predicted and discussed as well as the distributions of oil film thickness, pressure, temperature. The predicted temperature shows good agreement with measurements, while the pressure shows a reasonable distribution comparing with previous studies. Further analysis of the three-coupled-field reveals the reason of the high pressure and high temperature generated in the film. Finally, the influence of rotational speed of the mirror plate on the lubrication characteristics is illustrated which shows the thrust load should be balanced against the oil film temperature and pressure in optimized designs. This research proposes a thrust bearing computation method by combining CFD and FEA which can do the TEHD analysis more accurately.
文摘Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.
文摘A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.
基金Project supported by the National Natural Science Foundation of China(No.10832007)the Shanghai Leading Academic Discipline Project(No.B206)
文摘The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.
基金Supported by the National Natural Science Foundation of China (50774064)
文摘Failure pressure is a key parameter in reservoir hydrofracturing operation. Existing analytical methods for calculating the failure pressure are based on the assumption that borehole fluid is under two extreme conditions: non-infiltration or complete infiltration. The assumption is not suitable for the actual infiltration process, and this will cause a great error in practical calculation. It shows that during the injection process, the dynamic variation in effective stress-dependent permeability has an influence on the infiltration, and the influence also brings about calculation errors. Based on the fluid-structure interaction and finite element method (FEM), considering partial infiltration during injection process, a numerical model for calculating rock failure pressure is established. According to the analysis of permeability test results and response-surface method, a new variation rule of rock permeability with the change of effective stress is presented, and the relationships among the permeability, confining pressure and pore pressure are proposed. There are some differences between the dynamic value of permeability-effective-stress coefficient observed herein and the one obtained by the classical theory. Combining with the numerical model and the dynamic permeability, a coupling method for calculating failure pressure is developed. Comparison of field data and calculated values obtained by various methods shows that accurate values can be obtained by the coupling method. The coupling method can be widely applied to the calculation of failure pressure of reservoirs and complex wells to achieve effective fracturing operation.
基金Supported by the National Key Basic Research Program of China(No.2014CB046405)
文摘As the power transmission system of an aircraft,a hydraulic pipeline system is equivalent to the " blood vessel" of the aircraft. With the development of aircraft hydraulic system to high pressure,high speed and high power ratio,the fluid-structure interaction vibration mechanism of hydraulic pipeline is more complex and the influence of friction coupling on vibration cannot be ignored. The fluid-structure interaction of hydraulic pipeline will lead to system vibration,lower reliability of system operation and even pipeline rupture. Taking a hydraulic pipeline of C919 aircraft wingtip as the research object,a 14-equation model of fluid-structure interaction vibration considering friction coupling effect is established in this paper. The effects of friction and fluid parameters on the pipeline fluid-structure interaction vibration characteristics are studied and verified by experiments. The research results will provide theoretical guidance for the analysis of the pipeline fluid-structure interaction vibration and have important theoretical significance and great engineering value for promoting the localization process of large aircraft.
