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.

Blade Wrap Angle Impact on Centrifugal Pump Performance:Entropy Generation and Fluid-Structure Interaction Analysis

The centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to analyze the influence of the blade wrap angle on flow characteristics and energy distribution of a centrifugal pump evaluated as a low specific speed with a value of 69.This study investigates six impellermodels that possess varying blade wrap angles(95°,105°,115°,125°,135°,and 145°)that were created while maintaining the same volute and other geometrical characteristics.The investigation of energy loss was conducted to evaluate the values of total and entropy generation rates(TEG,EGR).The fluid-structure interaction was considered numerically using the software tools ANSYS Fluent and ANSYSWorkbench.The elastic structural dynamic equation was used to estimate the structural response,while the shear stress transport k–ωturbulence model was utilized for the fluid domain modeling.The findings suggest that the blade wrap angle has a significant influence on the efficiency of the pump.The impeller featuring a blade wrap angle of 145°exhibits higher efficiency,with a notable increase of 3.76%relative to the original model.Variations in the blade wrap angle impact the energy loss,shaft power,and pump head.The model with a 145°angle exhibited a maximum equivalent stress of 14.8MPa and a total deformation of 0.084 mm.The results provide valuable insights into the intricate flow mechanism of the centrifugal pump,particularly when considering various blade wrap angles.

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

An extended multiple-support response spectrum method incorporating fluid-structure interaction for seismic analysis of deep-water bridges

The effects of ground motion spatial variability(GMSV)or fluid-structure interaction(FSI)on the seismic responses of deep-water bridges have been extensively examined.However,there are few studies on the seismic performance of bridges considering GMSV and FSI effects simultaneously.In this study,the original multiple-support response spectrum(MSRS)method is extended to consider FSI effect for seismic analysis of deep-water bridges.The solution of hydrodynamic pressure on a pier is obtained using the radiation wave theory,and the FSI-MSRS formulation is derived according to the random vibration theory.The influence of FSI effect on the related coefficients is analyzed.A five-span steel-concrete continuous beam bridge is adopted to conduct the numerical simulations.Different load conditions are designed to investigate the variation of the bridge responses when considering the GMSV and FSI effects.The results indicate that the incoherence effect and wave passage effect decrease the bridge responses with a maximum percentage of 86%,while the FSI effect increases the responses with a maximum percentage of 26%.The GMSV and FSI effects should be included in the seismic design of deep-water bridges.

Analysis of the Influence of the Blade Deformation on Wind Turbine Output Power in the Framework of a Bidirectional Fluid-Structure Interaction Model

The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.

Analysis and Optimization of Flow-Guided Structure Based on Fluid-Structure Interaction

Gases containing sulfur oxides can cause corrosion and failure of bellows used as furnace blowers in high-temperature environments.In order to mitigate this issue,the behavior of an effective blast furnace blower has been examined in detail.Firstly,the Sereda corrosion model has been introduced to simulate the corrosion rate of the related bellows taking into account the effects of temperature and SO_(2) gas;such results have been compared with effective measurements;then,the average gas velocity in the pipeline and the von Mises stress distribution of the inner draft tube have been analyzed using a Fluid-Structure Interaction model.Finally,the semi-closed internal corrosion environment caused by a 5 mm radial gap between the inner draft tube and the bellows has been considered.The gas flow rate in the residential space has been found to be low(0.5 ms–this value leads to a stable semi-closed internal corrosion environment for exhaust gas exchange);water phase in the exhaust gas is prone to accelerate the corrosion rate.On this basis,a bellows with an optimized inner draft tube has proposed,which includes corrosion-resistant honeycomb buffer rings.

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.

一种评估非接触水下爆炸对海军复合船体影响的设计方法

Despite the non-contact underwater explosion phenomena(UNDEX)have been studied for decades and several numerical methods have been proposed in literature,its effects on military structures,especially composite ones,are even nowadays matter of research.In early design phases,it is not always possible to verify the shock resistance of hull structures modelling the whole phenomenon,in which fluid,gas and solid properties must be properly set in a fully coupled fluid-structure interaction(FSI)numerical model.These ones are extremely complex to set,computationally demanding and certainly not suitable for everyday design practice.In this paper,a simplified finite element(FE)model,easy to use in an early design phase,is proposed.Both,the structure and the fluid are simulated.In this approximation,the fluid behaviour is simplified,using special finite elements,available in a commercial software environment.This choice reduces the computational time and numerical efforts avoiding the problem of combining computational fluid dynamics(CFD)and FE domains and equations in a fully coupled fluid-structure interaction model.A typical parallel body block of a minesweeper is modelled,using two-dimensional multi-layered shell elements to properly account for the composite materials behaviour.For the fluid instead,three dimensional volumetric elements,directly coupled to the structural elements,are placed.In addition,the same calculation is performed,modelling separately fluid in the CFD environment and structures in the finite element one.Thus,realizing a fully coupled fluid-structure interaction model.The results obtained by applying both numerical models are compared with the structural response measured on board of a composite ship during a full-scale shock test.The simplified proposed procedure provides results in satisfactory agreement with experiments,allowing the validation of the model.Approximations are discussed and differences with the real phenomenon and fully coupled CFD+FE method are shown,providing a better understanding of the phenomena.Eventually,the modelling strategy has been considered a valuable and cost-effective tool for the concept and preliminary design of composite structures subject to underwater explosions.

Numerical study of submerged bending vegetation under unidirectional flow

Submerged vegetation commonly grows and plays a vital role in aquatic ecosystems,but it is also regarded as a barrier to the passing flow.Numerical simulations of flow through and over submerged vegetation were carried out to investigate the effect of vegetation density on flow field.Numerical simulations were computationally set up to replicate flume experiments,in which vegetation was mimicked with flexible plastic strips.The fluid-structure interaction between flow and flexible vegetation was solved by coupling the two modules of the COMSOL packages.Two cases with different vegetation densities were simulated,and the results were successfully validated against the experimental data.The contours of the simulated time-averaged streamwise velocity and Reynolds stress were extracted to highlight the differences in mean and turbulent flow statistics.The turbulence intensity was found to be more sensitive to vegetation density than the time-averaged velocity.The developing length increased with the spacing between plants.The snapshots of the bending vegetation under instantaneous velocity and vorticity revealed that flexible vegetation responded to the effects of eddies in the shear layer by swaying periodically.The first two rows of vegetation suffered stronger approaching flow and were prone to more streamlined postures.In addition,the origin of tip vortices was investigated via the distribution of vorticity.The results reveal the variation of flow properties with bending submerged vegetation and provide useful reference for optimizationofrestorationprojects.

An improved algorithm for fluid-structure interaction of high-speed trains under crosswind

Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.

Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.

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.

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.

OpenIFEM:A High Performance Modular Open-Source Software of the Immersed Finite Element Method for Fluid-Structure Interactions

We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.

Fluid-Structure Interaction Simulation of Aqueous Outflow System in Response to Juxtacanalicular Meshwork Permeability Changes with a Two-Way Coupled Method

Elevated intraocular pressure appears to have a broader impact on increased resistance to aqueous humor outflow through the conventional aqueous outflow system(AOS).However,there is still no consensus about exact location of the increased outflow resistance of aqueous humor,and the mechanism is not perfect.In addition,it is difficult to accurately obtain hydrodynamic parameters of aqueous humor within the trabecular meshwork outflow pathways based on the current technology.In this paper,a two-way fluid-structure interaction simulation was performed to study the pressure difference and velocity in the superficial trabecular meshwork,juxtacanalicular meshwork(JCM)and Schlemm’s canal in response to JCM permeability changes.We obtained the JCM permeability of normal intraocular pressure varied between 1×10?15 m2 and 10×10?15 m2 while permeability of the JCM ranged from 2×10?16 m2 and 3×10?16 m2 under conditions of high intraocular pressure.The study indicated that the fluid dynamics parameters in trabecular meshwork and Schlemm’s canal are most significantly affected by the changes of JCM permeability.Moreover,the study demonstrates that the finite element modeling of AOS provides a practical means for studying the outflow dynamics and the biomechanical environment of the AOS.

Microstructural analysis of a subhorizontal gold-quartz vein deposit at Donalda,Abitibi greenstone belt,Canada:Implications for hydrodynamic regime and fluid-structural relationship

The Donalda gold deposit in the southern part of the Archean Abitibi greenstone belt consists mainly of a subhonzontal gold-quartz vein perpendicular to subvertical shear zones.The 0.3—0.5 m thick vein is characterized by vein-parallel banding structures indicating multiple episodes of fracture opening and mineral precipitation.Measurement of the c-axis of primary growth quartz indicates that quartz preferentially grew perpendicular to the fracture,suggesting open space filling and/or extensional nature of the fracture.Measurement of the orientations of microfractures,veinlets and fluid—inclusion planes（FIPs） crosscutting primary growth quartz indicates that the vein minerals were subject to a vertical maximum principal stress（σ_1）,which is inconsistent with the subhorizontalσ_1 inferred from the regional stress field with N—S shortening.This apparent discrepancy is explained by invoking episodic fluid pressure fluctuation between supralithostatic and hydrostatic regimes accompanied by episodic opening and closing of the sub-horizontal fracture.When fluid pressure was higher than the lithostatic value,the fracture was opened and primary growth minerals were precipitated,whereas when fluid pressure decreased toward the hydrostatic value,the hanging wall of the fracture collapsed,causing collision of protruding primary growth minerals from both sides of the fracture and resulting in formation of vein-parallel deformation bands.The columns where the two facing sides of the fracture collided were subject to higher-than-lithostatic stress due to the bridging effect and reduced support surface area,explaining the development of verticalσ_1.This hypothesis is consistent the fault-valve model,and explains the flipping ofσ_1 without having to change the regional stress field.

Effects of renal artery stenosis on realistic model of abdominalaorta and renal arteries incorporating fluid-structureinteraction and pulsatile non-Newtonian blood flow

The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.

3D numerical simulation on fluid-structure interaction of structure subjected to underwater explosion with cavitation

In the underwater-shock environment, cavitation occurs near the structural surface. The dynamic response of fluid-structure interactions is influenced seriously by the cavitation effects. It is also the difficulty in the field of underwater explosion. With the traditional boundary element method and the finite element method （FEM）, it is difficult to solve the nonlinear problem with cavitation effects subjected to the underwater explosion. To solve this problem, under the consideration of the cavitation effects and fluid compressibility, with fluid viscidity being neglected, a 3D numerical model of transient nonlinear fluid-structure interaction subjected to the underwater explosion is built. The fluid spectral element method （SEM） and the FEM are adopted to solve this model. After comparison with the FEM, it is shown that the SEM is more precise than the FEM, and the SEM results are in good coincidence with benchmark results and experiment results. Based on this, combined with ABAQUS, the transient fluid-structure interaction mechanism of the 3D submerged spherical shell and ship stiffened plates subjected to the underwater explosion is discussed, and the cavitation region and its influence on the structural dynamic responses are presented. The paper aims at providing references for relevant research on transient fluid-structure interaction of ship structures subjected to the underwater explosion.