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.

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.

Towards Production and Energy Coupling System Modeling and Simulation for Energy Optimization in the Process Industry

The production and energy coupling system is used to mainly present energy flow, material flow, information flow, and their coupling interaction. Through the modeling and simulation of this system, the performance of energy flow can be analyzed and optimized in the process industry. In order to study this system, the component based hybrid Petri net methodology (CpnHPN) is proposed, synthesizing a number of extended Petri net methods and using the concept of energy place, material place, and information place. Through the interface place in CpnHPN, the component based encapsulation is established, which enables the production and energy coupling system to be built, analyzed, and optimized on the multi-level framework. Considering the block and brief simulation for hybrid system, the CpnHPN model is simulated with Simulink/Stateflow. To illustrate the use of the proposed methodology, the application of CpnHPN in the energy optimization of chlorine balance system is provided.

Numerical simulation of a sheet metal extrusion process by using thermal-mechanical coupling EAS FEM

The thermal-mechanical coupling finite element method(FEM)was usedto simulate a non-isothermal sheet metal extrusion process. On thebasis of the finite plasticity consistent with multiplicativedecomposition of the deformation gradient, the enhanced as- sumedstrain(EAS)FEM was applied to carry out the numerical simulation. Inorder to make the computation reliable ad avoid hour- glass mode inthe EAS element under large compressive strains, an alterative formof the original enhanced deformation gradient was employed. Inaddition, reduced factors were used in the computation of the elementlocal internal parameters and the enhanced part of elementalstiffness.

The Mesozoic Basin-Mountain Coupling Process of the Southern East China Sea Shelf Basin and its Adjacent Land Area

Basin-mountain coupling is a key issue for basin formation and evolution. The analysis of basin-mountain coupling process, as well as quantitative or semiquantitative restoration of prototype basin and the evolution of continental margin, can be used to interpret the geological process of basin-range conversion and reconstruct early prototype basins, which is a difficult and leadin~ scientific oroblem of basin research.

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.

Formation Mechanism in Alloy Steel Rolling Process Using Thermo-mechanical Coupling Method

Based on the theory of elastic-plastic finite element method, the high-speed hot continuous rolling process of a billet is simulated and analyzed in vertical and horizontal passes. The billet is dragged into the passes by contact friction force between the billet and rollers. The rollers and billet are represented by respectively rigid and deformable bodies, and three-dimensional models are developed for the billet and rollers. The distribution of deformation field, effective strain, rolling force and temperature field are accurately calculated for the whole rolling process （including unstable and stable stages）. In addition, the rolling pressure on the width symmetry center is compared with that in the in-situ experimental measurements. It is revealed that various heat exchange phenomena among the billet, rollers and surroundings can result in unbalanced temperature distribution on the cross section. Rolling force and strain can change significantly when the billet is moved towards or away from the roller gap, and keep almost invariable in the stable stage. It is expected that the simulation results would be useful for practical manufacture and provide the theoretical foundation for improvement of process planning and optimization of process parameters.

Simulation of fluid-structure interaction in a microchannel using the lattice Boltzmann method and size-dependent beam element on a graphics processing unit

Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.

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.

Effects of Tensor Couplings on Nucleonic Direct URCA Processes in Neutron Star Matter

The relativistic neutrino emissivity of the nucleonic direct URCA processes in neutron star matter is investigated within the relativistic Hartree-Fock approximation. We particularly study the influences of the tensor couplings of vector mesons ω and ρ on the nucleonic direct URCA processes. It is found that the inclusion of the tensor couplings of vector mesons w and p can slightly increase the maximum mass of neutron stars. In addition, the results indicate that the tensor couplings of vector mesons ω and ρ lead to obvious enhancement of the total neutrino emissivity for the nucleonic direct URCA processes, which must accelerate the cooling rate of the non- superfluid neutron star matter. However, when considering only the tensor coupling of vector meson ρ, the neutrino emissivity for the nucleonic direct URCA processes slightly declines at low densities and significantly increases at high densities. That is, the tensor coupling of vector meson ρ leads to the slow cooling rate of a low-mass neutron star and rapid cooling rate of a massive neutron star.

Decoupling multimode vibrational relaxations in multi-component gas mixtures: Analysis of sound relaxational absorption spectra

Decoupling the complicated vibrational-vibrational （V-V） coupling of a multimode vibrational relaxation remains a challenge for analyzing the sound relaxational absorption in multi-component gas mixtures. In our previous work [Acta Phys. Sin. 61 174301 （2012）], an analytical model to predict the sound absorption from vibrational relaxation in a gas medium is proposed. In this paper, we develop the model to decouple the V-V coupled energy to each vibrationaltranslational deexcitation path, and analyze how the multimode relaxations form the peaks of sound absorption spectra in gas mixtures. We prove that a multimode relaxation is the sum of its decoupled single-relaxation processes, and only the decoupled process with a significant isochoric-molar-heat can be observed as an absorption peak. The decoupling model clarifies the essential processes behind the peaks in spectra arising from the multimode relaxations in multi-component gas mixtures. The simulation validates the proposed decoupling model.

Coupling catalytic hydrolysis and oxidation for CS_2 removal

CS2 removal was obtained by coupling catalytic hyidation on bi-functional catalyst. On the hydrolysis active sites, CS2 is hydrolyzed to HES, while on the oxidation active sites, HES is oxidized to elemental S or sulfuric acid deposited on the porous support. The above process can be expressed as follows： CS2 H2O→ COS H2O→ H2S O2→ S/SO4^2-. H2S oxidation eliminates its prohibition on C52 hydrolysis so that the rate of coupling removal CS2 is 5 times higher than that of CS2 hydrolysis. The same active energy of hydrolysis and coupling reaction also indicates that HES oxidation does not change the reaction mechanism of CS2 hydrolysis. Temperature has obvious effect on the process while the mole ratio of O2 concentration to CS2 concentration （O/S） does not, especially in excess of 2.5. The formation of sulfuric acid on the catalyst surface poisons hydrolysis active sites and causes the decrease of left OH^-1 concentration on the catalysts surface. Lower temperature is suggested for this bi-functional catalyst owing to the low yield ratio of S/SO4^2-.

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.