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.

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.

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.

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.

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.

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

Stochastic Huge-Resonance Caused by Coupling for a Globally Coupled Linear System

In the paper, we investigate a globally coupled linear system with finite subunits subject to temporal periodic force and with multiplicative dichotomous noise. It is shown that, the global coupling among the subunits can hugely enhance the phenomenon of SR for the amplitude of the average mean field as the functions of the transition rate of the noise and that as the function of the frequency of the signal respectively.

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

Theoretical and Numerical Studies on the Coupling Deformation of Global Lateral Buckling and Walking of Submarine Pipeline

Buckling initiation devices/techniques,including sleepers,distributed buoyancy,snake lay,and residual curvature method(RCM),have recently been widely applied in engineering.These initiated buckles may induce a long pipeline to transform into multiple short pipeline segments,which promote the occurrence of pipeline walking.Thus,a pipeline,which is designed to buckle laterally,may laterally and axially displace over time when subjected to repeated heating and cooling cycles.This study aims to reveal the coupling mechanism of pipeline walking and global lateral buckling.First,an analytic solution is proposed to estimate the walking of pipeline segments between two adjacent buckles.Then,the sensitivity of this method to heating and cooling cycles is analyzed.Results show the applicability of the proposed walking analytical solution of buckling pipelines.Subsequently,an influence analysis of walking on global buckling,including the capacity of buckling initiation,buckling amplitude,buckling mode,and failure assessment of the buckling pipeline,is performed.The results reveal that the effect of walking on the buckling axial force is negligible.However,pipeline walking will aggravate the asymmetry of the pipeline buckling and the failure parameters of the pipeline during the post-buckling.

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.

GLOBAL SOLUTIONS TO A HYPERBOLIC-PARABOLIC COUPLED SYSTEM WITH LARGE INITIAL DATA

This paper is concerned with the existence of global solutions to the Cauchy problem of a hyperbolic-parabolic coupled system with large initial data. To this end, we first construct its local solutions by the standard iteration technique, then we deduce the basic energy estimate by constructing a convex entropy-entropy flux pair to this system. Moreover, the L∞-estimates and H^2-estimates of solutions are obtained through some delicate estimates. In our results, we do not ask the far fields of the initial data to be equal and the initial data can be arbitrarily large which generalize the result obtained in [7].

A Coupled General Circulation Model for the Tropical Pacific Ocean and Global Atmosphere

On the basis of Zeng's theorehcal design, a coupled general circulation model(CGCM) is develO￣ with itscharacteristics different from other CGCMs such as the unified vertical coordinates and subtraction of the standard stratification for both atmosphere and ocean, available energy consideration,and so on.The oceanic comPOnent is a free surface tropical Pacific Ocean GCM betWeen 30W and 30'S with horizontal grid spacing of ic in latitude and 2°in longitude,and with 14 vertical layers.The atmospheric component is a global GCM with low-resolution of 4°in lahtude and 5°in longitude,and tWo layers of equal mass in the verhcal between the surfaCe and 200 hFa.The atmospheric GCM includes comprehensive physical processes.The coupled model is subjected to seasonally-varying cycle.Several coupling experiments,ranging from straight forward coupling without flux correction to one with flux correchon,and to so-called predictor-corrector monthly coupling(PCMC),are conducted tO show the esistence and final controlling of the climate drift in the coupled system.After removing the climate drift with the PCMC SCheme,the coupled model is integrated for more than twenty years.The results show reasonable simulations of the anneal mean and its seasollal cycle of the atmospheric and￣ante circulahon.The model also ProduCeS the coherent intermnual variations of the climate system, manifesting the observed EI Nifio/Southern OSCillation(ENSO).

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.

Application of 1D/3D finite elements coupling for structural nonlinear analysis

An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.

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.