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.

Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

AN ADAPTIVE EFG-FE COUPLING METHOD FOR THE NUMERICAL SIMULATION OF EXTRUSION PROCESSES

An adaptive EFG-FE coupling method is proposed and developed for the numerical simulation of lateral extrusion and forward-backward extrusion. Initially, the simulation has been implemented by using a conventional FE model. During the deforming process, mesh quality is checked at every incremental step. Distorted elements are automatically converted to EFG nodes, whereas, the less distorted elements are reserved. A new algorithm to generate EFG nodes and interface elements is presented. This method is capable of dealing with large deformation and has higher computational efficiency than using an EFG method wholly. Numerical results demonstrate that the adaptive EFG-FE coupling method has reasonable accuracy and is effective for local bulk metal forming such as extrusion processes.

Numerical Simulation of Liquid-Solid Coupling in Multi-Angle Fractures in Pressure-Sensitive Reservoirs

This study clarifies the seepage characteristics of complex fractured pressure-sensitive reservoirs,and addresses a common technological problem,that is the alteration of the permeability degree of the reservoir bed(known to be responsible for changes in the direction and velocity of fluid flows between wells).On the basis of a new pressuresensitive equation that considers the fracture directional pressure-sensitive effect,an oil-gas-water three-phase seepage mathematical model is introduced,which can be applied to pressure-sensitive,full-tensor permeability,ultralow-permeability reservoirs with fracture-induced anisotropy.Accordingly,numerical simulations are conducted to explore the seepage laws for ultralow-permeability reservoirs.The results show that element patterns have the highest recovery percentage under a fracture angle of 45°.Accounting for the pressure-sensitive effect produces a decrease in the recovery percentage.Several patterns are considered:inverted five-seven-and nine-spot patterns and a cross-row well pattern.Finally,two strategies are introduced to counteract the rotation of the direction of the principal permeability due to the fracture directional pressure-sensitive effect.

Numerical Simulation of Coupled Thermo-mechanical Behavior of a Cylinder Billet during Hot-forging Process

A finite volume method is applied to simulate a closed die hot forging process of a cylinder billet. Since variation and distribution of temperature play very important role in hot forging, the code involves a methodology of a coupled system of mechanical and thermal equations. The simulated results are compared with the experimental ones. The distribution of temperature in the billet obtained from the simulation is also discussed.

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.

Numerical simulation of intermediate phase growth in Ti/Al alternate foils

To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.

Numerical Simulation of Projection Welding Processes for Door Hinge of Automobile Based on Coupled Fields Analysis

Projection welding is a variation of electric resistance welding with the dynamic changes of the flow paths for heat and electrical properties with changing temperature caused by the large plastic deformation collapse of projection. As the joint type between the auto door hinge and the inner plate, projection welding may bring welding distortions and would affect the assembly quality of auto body. A comprehensive electric-thermal-mechanical numerical simulation was performed to quantitatively simulate the processes of projection welding by using a coupled finite element method. The mechanism of projection collapse and the formation process of nugget were discussed and good conclusions have been achieved comparing with the test results.

Numerical Simulation of Interaction Between Laminar Flow and Elastic Sheet

A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.

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.

A Fluid-Structure Interaction Simulation of Coal and Gas Outbursts Based on the Interaction between the Gas Pressure and Deformation of a Coal-Rock Mass

Based on the theories of the gas seepage in coal seams and the deformation of the coal-rock medium,the gas seepage field in coal-rock mass is coupled with the deformation field of the coal-rock mass to establish a fluidstructure interaction model for the interaction between coal gas and coal-rock masses.The outburst process in coal-rock masses under the joint action of gas pressure and crustal stress is simulated using the material point method.The simulation results show the changes in gas pressure,velocity distribution,maximum principal stress distribution,and damage distribution during the process of the coal and gas outburst,as well as themovement and accumulation of coal-rock masses after the occurrence of the outburst.It was found that the gas pressure gradient was greatest at theworking face after the occurrence of the outburst,the gas pressures and pressure gradients at each location within the coal seamgradually decreased with time,and the damage distribution was essentially the same as the minimum principal stress distribution.The simulation further revealed that the outburst first occurred in themiddle of the tunnel excavation face and that the speed at which particles of coal mass were ejected was highest at the center and decreased toward the upper and lower sides.The study provides a scientific basis for enhancing our understanding of the mechanism behind coal and gas outbursts,as well as their prevention and control.

Modeling and Simulation of Valve Cycle in Vein Using an Immersed Finite Element Method

A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.

Simulation Monitoring for Rainfall Infiltration in Soil Based on High Density Electrical Method

Rainfall infiltration is a porous medium flow problem with variable saturation. Based on the theoretical analysis of the flow field, electrical conductivity of rocks, the electrical field, the paper simulates the coupling relationship between the water saturation in soil and the apparent resistivity distribution. It combines the Richards equation, the Archie formula and the Laplace equation. The experiment simulates the potential field data by the Wenner setting in electrical exploration on a two-layer geologic model with continuous rainfall during 5 days, which shows that the effective saturation in soil is increasing with the rainfall time, while the apparent resistivity is decreasing. This can provide a theoretical basis for the analyzing the rainfall infiltration and porosity of the soil by using high-density electrical method in the future.

THE UPWIND FINITE DIFFERENCE METHOD FOR MOVING BOUNDARY VALUE PROBLEM OF COUPLED SYSTEM

Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.

Geometrically Nonlinear Flutter Analysis Based on CFD/CSD Methods and Wind Tunnel Experimental Verification

This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the aerodynamic simulation was conducted based on computational fluid dynamics,and the structural model was established using the nonlinear finite element model and tangential stiffness matrix.First,the equilibrium position was obtained using the nonlinear static aeroelastic iteration.Second,the structural modal under a steady aerodynamic load was extracted.Finally,the generalized displacement time curve was obtained by coupling the unsteady aerodynamics and linearized structure motion equations.Moreover,if the flutter is not at a critical state,the incoming flow dynamic pressure needs to be changed,and the above steps must be repeated until the vibration amplitude are equal.Furthermore,the high-speed geometrically nonlinear flutter of the wing-body assemblymodel with a high-aspect ratio was investigated,and the correctness of the method was verified using high-speed wind tunnel experiments.The results showed that the geometric nonlinearity of the large deformation of the wing caused in-plane bending to become a key factor in flutter characteristics and significantly decreased the dynamic pressure and frequency of the nonlinear flutter compared to those of the linear flutter.

基于微观流固耦合的超深层致密砂岩气藏应力敏感性分析

超深层致密砂岩气藏具有强应力敏感性,目前常用的研究方法包括压力脉冲实验法和实时在线CT扫描法两种,压力脉冲实验法不能揭示产生应力敏感性的微观机理,而实时在线CT扫描法也无法模拟深部地层高压、高应力的条件。为解决应力敏感性实验研究的不足,基于离散单元法与管道网络模型建立了微观流固耦合算法,编制了模拟器,并对模拟器力学计算和流固耦合模块的正确性进行了验证,分析了应力大小、加载方向和孔隙压力对岩心渗透率的影响,最后从微观上揭示了超深层致密砂岩气藏的应力敏感性机理。研究结果表明:①应力通过增加与之垂直方向上喉道两侧的法向压力,减小喉道的水力半径,进而降低储层的渗透率;②较高的孔隙压力能够阻碍岩石颗粒在应力作用下的移动,从而减缓了孔隙和喉道的变形,使模型保持较高的渗透率;③致密砂岩气藏的渗透率受到应力和地层压力的共同控制,并且具有各向异性,在垂直于最小主应力方向上形成渗透率较大的优势通道;④异常高压阻碍了地应力的压实作用,有利于保护储层孔隙,使地层有较好的储集性能和较高的渗透率。结论认为,根据离散元法结合孔隙网络模型建立的流固耦合方法可为理解超深层致密砂岩应力敏感性提供理论参考,并为超深层致密砂岩气藏的科学高效开发提供指导。

A special numerical method for fluid-structure interaction problems subjected to explosion and impact loading

A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.

The fast method and convergence analysis of the fractional magnetohydrodynamic coupled flow and heat transfer model for the generalized second-grade fluid

In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.

火箭仪器舱结构模型的热振耦合模拟方法

传统的热振耦合方法缺乏对材料热特性变化的考虑,影响仿真分析的精度。为此,文章提出一种热特性参数与振动动态耦合的分析方法,建立动态变化的热振耦合方程,基于闪光瞬态法得到材料的热特性参数并将其应用在结构传热过程的数值仿真计算中。结果表明,是否考虑材料的热特性变化会对结构的传热过程仿真结果产生较大影响,并进一步影响结构的刚度特性仿真结果,从而造成结构模态分析偏差。采用本研究提出的动态变化的热振耦合模拟方法可提高传热仿真分析的精度和可靠度。

某小口径火炮弹丸挤进过程动力学特性分析

为研究某小口径火炮不同坡膛结构下弹丸挤进过程的动力学特性,建立了考虑挤进阻力变化的内弹道方程组并通过vuamp子程序引入到有限元求解中,实现了挤进过程中内弹道计算与弹丸挤进有限元仿真耦合求解。通过对弹丸全膛运动的仿真分析验证了模型的可行性,随后得到并分析不同坡膛锥度、膛线深度及阴线宽度下挤进过程中弹丸的速度、加速度、挤进阻力和膛压曲线及变化影响规律。结果表明,坡膛锥度对挤进过程中弹丸的动力学参数影响不是单调的,而膛线深度及阴线宽度的改变则单调地影响弹丸动力学参数,对优化坡膛结构及弹带设计具有一定的指导意义。