Simulation of Underground Reservoir Stability of Pumped Storage Power Station Based on Fluid-Structure Coupling

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.

Application of CFD and FEA Coupling to Predict Structural Dynamic Responses of A Trimaran in Uni-and Bi-Directional Waves

To predict the wave loads of a flexible trimaran in different wave fields,a one-way interaction numerical simulation method is proposed by integrating the fluid solver(Star-CCM+)and structural solver(Abaqus).Differing from the existing coupled CFD-FEA method for monohull ships in head waves,the presented method equates the mass and stiffness of the whole ship to the hull shell so that any transverse and longitudinal section stress of the hull in oblique waves can be obtained.Firstly,verification study and sensitivity analysis are carried out by comparing the trimaran motions using different mesh sizes and time step schemes.Discussion on the wave elevation of uni-and bi-directional waves is also carried out.Then a comprehensive analysis on the structural responses of the trimaran in different uni-directional regular wave and bi-directional cross sea conditions is carried out,respectively.Finally,the differences in structural response characteristics of trimaran in different wave fields are studied.The results show that the present method can reduce the computational burden of the two-way fluid-structure interaction simulations.

Fluid-Structure Coupled Analysis of the Transient Thermal Stress in an Exhaust Manifold

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).

VIBRATION ANALYSIS OF TURBINE BASED ON FLUID-STRUCTURE COUPLING

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.

Improved frequency modeling and solution for parallel liquid-filled pipes considering both fluid-structure interaction and structural coupling

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.

Fluid-Structure Interaction Analysis of Flexible Plate with Partitioned Coupling Method

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.

Transverse seismic response of beam aque-duct considering fluid-structure coupling

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.

Methodology for Comparing Coupling Algorithms for Fluid-Structure Interaction Problems

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.

Fluid-structure coupled analysis of underwater cylindrical shells

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.

Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

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.

Fully Coupled Fluid-Structure Interaction Model Based on Distributed Lagrange Multiplier/Fictitious Domain Method

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.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

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.

3D Two-way Coupled TEHD Analysis on the Lubricating Characteristics of Thrust Bearings in Pump-turbine Units by Combining CFD and FEA

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.

A Novel Approach for the Numerical Simulation of Fluid-Structure Interaction Problems in the Presence of Debris

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.

Fluid-structure interaction simulation of three-dimensional flexible hydrofoil in water tunnel

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 calculation of fracturing well based on fluid-structure interaction

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.

Friction coupling vibration characteristics analysis of aviation hydraulic pipelines considering multi factors

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.

Finite element analysis of fluid-structure interaction in buried liquid-conveying pipeline

Long distance buried liquid-conveying pipeline is inevitable to cross faults and under earthquake action,it is necessary to calculate fluid-structure interaction(FSI) in finite element analysis under pipe-soil interaction.Under multi-action of site,fault movement and earthquake,finite element model of buried liquid-conveying pipeline for the calculation of fluid structure interaction was constructed through combinative application of ADINA-parasolid and ADINA-native modeling methods,and the direct computing method of two-way fluid-structure coupling was introduced.The methods of solid and fluid modeling were analyzed,pipe-soil friction was defined in solid model,and special flow assumption and fluid structure interface condition were defined in fluid model.Earthquake load,gravity and displacement of fault movement were applied,also model preferences.Finite element research on the damage of buried liquid-conveying pipeline was carried out through computing fluid-structure coupling.The influences of pipe-soil friction coefficient,fault-pipe angle,and liquid density on axial stress of pipeline were analyzed,and optimum parameters were proposed for the protection of buried liquid-conveying pipeline.

Simulation and Analysis of China Climate Using Two-Way Interactive Atmosphere-Vegetation Model (RIEMS-AVIM)

In this paper, an Atmosphere-Vegetation Interaction Model （AVIM） is coupled to the Regional Integrated Environment Model System （RIEMS）, and a 10-year integration for China is performed using the RIEMS-AVIM. The analysis of the results of the 10-year integration shows that the characters of the spatial distributions of temperature and precipitation over China are well simulated. The patterns of simulated surface sensible and latent heat fluxes match well with the spatial climatological atlas： the values of winter surface sensible and latent heat fluxes are both lower than climatological values over the whole country. Summer surface sensible heat flux is higher than climatological values in western China and lower in eastern China, while summer surface latent heat flux is higher than climatological values in the eastern and lower in the western. Seasonal variations of simulated temperature and precipitation of RIMES-AVIM agree with those of the observed. Simulated temperature is lower than the observed in the Tibetan Plateau and Northwest China for the whole year, slightly lower in the remaining regions in winter, but consistent with the observed in summer. The simulated temperature of RIEMS-AVIM is higher in winter and lower in summer than that of RIEMS, which shows that the simulated temperature of RIEMS-AVIM is closer to the observed value. Simulated precipitation is excessive in the first half of the year, but consistent with the observed in the second half of the year. The simulated summer precipitation of RIEMS-AVIM has significant improvement compared to that of RIEMS, which is less and closer to the observed value. The interannual variations of temperature and precipitation are also fairly well simulated, with temperature simulation being superior to precipitation simulation. The interannual variation of simulated temperature is significantly correlated with the observed in Northeast China, the Transition Region, South China, and the Tibetan Plateau, but the correlation between precipitation simulation and observation is only significant in Northwest China.

MPS-FEM Coupled Method for Study of Wave-Structure Interaction

Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.