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.

CAS-ESM2.0 Successfully Reproduces Historical Atmospheric CO_(2) in a Coupled Carbon−Climate Simulation

The atmospheric carbon dioxide(CO_(2))concentration has been increasing rapidly since the Industrial Revolution,which has led to unequivocal global warming and crucial environmental change.It is extremely important to investigate the interactions among atmospheric CO_(2),the physical climate system,and the carbon cycle of the underlying surface for a better understanding of the Earth system.Earth system models are widely used to investigate these interactions via coupled carbon-climate simulations.The Chinese Academy of Sciences Earth System Model version 2(CAS-ESM2.0)has successfully fixed a two-way coupling of atmospheric CO_(2)with the climate and carbon cycle on land and in the ocean.Using CAS-ESM2.0,we conducted a coupled carbon-climate simulation by following the CMIP6 proposal of a historical emissions-driven experiment.This paper examines the modeled CO_(2)by comparison with observed CO_(2)at the sites of Mauna Loa and Barrow,and the Greenhouse Gases Observing Satellite(GOSAT)CO_(2)product.The results showed that CAS-ESM2.0 agrees very well with observations in reproducing the increasing trend of annual CO_(2)during the period 1850-2014,and in capturing the seasonal cycle of CO_(2)at the two baseline sites,as well as over northern high latitudes.These agreements illustrate a good ability of CAS-ESM2.0 in simulating carbon-climate interactions,even though uncertainties remain in the processes involved.This paper reports an important stage of the development of CAS-ESM with the coupling of carbon and climate,which will provide significant scientific support for climate research and China’s goal of carbon neutrality.

Coupled Numerical Simulation of Electromagnetic and Flow Fields in a Magnetohydrodynamic Induction Pump

Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.

Multi-scenario Simulation and Spatial-temporal Analysis of LUCC in China's Coastal Zone Based on Coupled SD-FLUS Model

Increased human activities in China's coastal zone have resulted in the depletion of ecological land resources.Thus,conducting current and future multi-scenario simulation research on land use and land cover change(LUCC)is crucial for guiding the healthy and sustainable development of coastal zones.System dynamic(SD)-future land use simulation(FLUS)model,a coupled simulation model,was developed to analyze land use dynamics in China's coastal zone.This model encompasses five scenarios,namely,SSP1-RCP2.6(A),SSP2-RCP4.5(B),SSP3-RCP4.5(C),SSP4-RCP4.5(D),and SSP5-RCP8.5(E).The SD model simulates land use demand on an annual basis up to the year 2100.Subsequently,the FLUS model determines the spatial distribution of land use for the near term(2035),medium term(2050),and long term(2100).Results reveal a slowing trend in land use changes in China's coastal zone from 2000–2020.Among these changes,the expansion rate of construction land was the highest and exhibited an annual decrease.By 2100,land use predictions exhibit high accuracy,and notable differences are observed in trends across scenarios.In summary,the expansion of production,living,and ecological spaces toward the sea remains prominent.Scenario A emphasizes reduced land resource dependence,benefiting ecological land protection.Scenario B witnesses an intensified expansion of artificial wetlands.Scenario C sees substantial land needs for living and production,while Scenario D shows coastal forest and grassland shrinkage.Lastly,in Scenario E,the conflict between humans and land intensifies.This study presents pertinent recommendations for the future development,utilization,and management of coastal areas in China.The research contributes valuable scientific support for informed,long-term strategic decision making within coastal regions.

Shear behavior of intact granite under thermo-mechanical coupling and three-dimensional morphology of shear-formed fractures

The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.

Multi-dimensional Simulation of Phase Change by a 0D-2D Model Coupling via Stefan Condition

Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.

Spudcan Penetration Simulation Using the Coupled Eulerian-Lagrangian Method with Thermo-Mechanical Coupled Analysis

A novel modeling technique based on the coupled Eulerian-Lagrangian(CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mechanical properties during spudcan penetration in normally consolidated clay soil. In the CEL model, the normal method of assigning an increasing shear strength profile with depth(NA) is defective due to its Eulerian framework. In this paper, a new technique is proposed to update soil material properties by introducing thermo-mechanical coupled analysis(TMCA) to the CEL models. During establishment of the CEL models, the optimal penetration velocity and minimum mesh size are determined through parametric studies. Reasonability and accuracy are then verified through comparison of the preliminary results with the soil flow configuration and penetration resistance(Fv) of a centrifuge test, and the results of the proposed method are compared with those of the remeshing and interpolation technique with small strain(RITSS) method. To achieve a CEL model with satisfactory accuracy, the NA and TMCA methods implemented in the CEL models and the RITSS method are first adopted in weightless soil. Comparison of the findings with those obtained in previous studies shows that the TMCA method can update material properties and predict Fv. The TMCA method is then applied to soils with self-weight and different shear strength profiles. Results show that the proposed method is capable of accurately modeling the large deformation problem of spudcan penetration in non-homogeneous clay.

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.

Dynamic thermo-mechanical coupled simulation of statistically inhomogeneous materials by statistical second-order two-scale method

In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,the random distribution characteristics of particles,including the shape,size,orientation,spatial location,and volume fractions,are all considered.Firstly,the repre.sentation for the microscopic configuration of the statistically inhomogeneous materials is described.Secondly,the SSOTS formulation for the dynamic thermo-mechanical coupled problem is proposed in a constructive way,including the cell problems,effective thermal and mechanical parameters,homogenized problems,and the SSOTS formulas of the temperatures,displacements,heat flux densities and stresses.And then the algorithm procedure corresponding to the SSOTS method is brought forward.The numerical results obtained by using the SSOTS algorithm are compared with those by classical methods.In addition,the thermo-mechanical coupling effect is studied by comparing the results of coupled case with those of uncoupled case.It demonstrates that the coupling effect on the temperatures,heat flux densities,displacements,and stresses is very distinct.The results show that the SSOTS method is valid to predict the dynamic thermo-mechanical coupled performances of statistically inhomogeneous materials.

Influence of water coupling coefficient on the blasting effect of red sandstone specimens

This study investigates the impact of different water coupling coefficients on the blasting effect of red sandstone.The analysis is based on the theories of detonation wave and elastic wave,focusing on the variation in wall pressure of the blasting holes.Using DDNP explosive as the explosive load,blasting tests were conducted on red sandstone specimens with four different water coupling coefficients:1.20,1.33,1.50,and 2.00.The study examines the morphologies of the rock specimens after blasting under these different water coupling coefficients.Additionally,the fractal dimensions of the surface cracks resulting from the blasting were calculated to provide a quantitative evaluation of the extent of rock damage.CT scanning and 3D reconstruction were performed on the post-blasting specimens to visually depict the extent of damage and fractures within the rock.Additionally,the volume fractal dimension and damage degree of the post-blasting specimens are calculated.The findings are then combined with numerical simulation to facilitate auxiliary analysis.The results demonstrate that an increase in the water coupling coefficient leads to a reduction in the peak pressure on the hole wall and the crushing zone,enabling more of the explosion energy to be utilized for crack propagation following the explosion.The specimens exhibited distinct failure patterns,resulting in corresponding changes in fractal dimensions.The simulated pore wall pressure–time curve validated the derived theoretical results,whereas the stress cloud map and explosion energy-time curve demonstrated the buffering effect of the water medium.As the water coupling coefficient increases,the buffering effect of the water medium becomes increasingly prominent.

Metal–Organic Gel Leading to Customized Magnetic‑Coupling Engineering in Carbon Aerogels for Excellent Radar Stealth and Thermal Insulation Performances

Metal–organic gel(MOG)derived composites are promising multi-functional materials due to their alterable composition,identifiable chemical homogeneity,tunable shape,and porous structure.Herein,stable metal–organic hydrogels are prepared by regulating the complexation effect,solution polarity and curing speed.Meanwhile,collagen peptide is used to facilitate the fabrication of a porous aerogel with excellent physical properties as well as the homogeneous dispersion of magnetic particles during calcination.Subsequently,two kinds of heterometallic magnetic coupling systems are obtained through the application of Kirkendall effect.FeCo/nitrogen-doped carbon(NC)aerogel demonstrates an ultra-strong microwave absorption of−85 dB at an ultra-low loading of 5%.After reducing the time taken by atom shifting,a FeCo/Fe3O4/NC aerogel containing virus-shaped particles is obtained,which achieves an ultra-broad absorption of 7.44 GHz at an ultra-thin thickness of 1.59 mm due to the coupling effect offered by dual-soft-magnetic particles.Furthermore,both aerogels show excellent thermal insulation property,and their outstanding radar stealth performances in J-20 aircraft are confirmed by computer simulation technology.The formation mechanism of MOG is also discussed along with the thermal insulation and electromagnetic wave absorption mechanism of the aerogels,which will enable the development and application of novel and lightweight stealth coatings.

A thermo-mechanical damage constitutive model for deep rock considering brittleness-ductility transition characteristics

This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.

The Impact of Deformable Surrounding Rocks on Naturally-Fractured Reservoir Simulation

For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) methods are investigated. The explicit formulation is capable of modelling surrounding media;while the implicit coupling is unable to do so as deformation vector does not appear as a primary variable in the formulation. The governing differential equations and the finite element approximation of the governing equations for each of the methods are presented. Spatial discretization is achieved using the Galerkin method, and temporal discretisation using the finite difference technique. In the explicit model, coupling between flow and deformation is captured through volumetric strain compatibility amongst the phases within the system. In the implicit model, this is achieved by defining the pore space storativity as a function of the formation compressibility and the compressibility of the fluid phases within the pore space. The impact of rock deformability on early, intermediate and late time responses of fractured reservoir is investigated through several numerical examples. Salient features of each formulation are discussed and highlighted. It is shown that the implicit model is unable to capture the constraining effects of a non-yielding, surrounding rock, leading to incorrect projections of reservoir production irrespective of the history matching strategy adopted.

Fluid simulation of the effect of a dielectric window with high temperature on plasma parameters in inductively coupled plasma

To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.

Experimental and simulation studies on similitude design method for shock responses of beam-plate coupled structure

The similitude theory helps to understand the physical behaviors of large structures through scaled models. Several papers have studied the similitude of shock issues. However, the dynamic similitude for shock responses of coupled structures is rarely incorporated in open studies. In this paper, scaling laws are derived for the shock responses and spectra of coupled structures. In the presented scaling laws, the geometric distortion and energy loss are considered. The ability of the proposed scaling laws is demonstrated in the simulation and experimental cases. In both cases, the similitude prediction for the prototype's time-domain waveform and spectrum is conducted with the scaled model and scaling laws. The simulation and experimental cases indicate that the predicted shock responses and spectra agree well with those of the prototype, which verifies the proposed scaling laws for predicting shock responses.

Modeling and Simulation of a Transmission Line Response to a 400 kV/400V Capacitor Coupled Substation

The access to electricity in rural areas is extremely limited, but it is crucial for all citizens. The population in rural areas of sub-Saharan African (SSA) countries is generally low, making it economically unfeasible to implement traditional rural electrification (CRE) projects due to the high cost of establishing the necessary distribution infrastructure. To address this cost issue, one alternative technology for rural electrification (URE) that can be explored is the Capacitor Coupled Substation (CCS) technology. CCS is a cost-effective solution for supplying electricity to rural areas. The research is necessitated by the need to offer a cost-effective technology for supplying electricity to sparsely populated communities. This paper examines the impact on the transmission network when a 400 kV/400V CCS is connected to it. The system response when a CCS is connected to the network was modeled using MATLAB/Si-mulink. The results, based on the fixed load of 80 kW, showed negligible interference on the transmission line voltage. However, there was minor impact on the parameters downstream of the tapping point. These findings were further supported by introducing a fault condition to the CCS, which showed that interferences with the CCS could affect the overall stability of the transmission network downstream of the tapping node, similar to the behavior of an unstable load.

Numerical Simulation and Experimental Study of Heat-fluid-solid Coupling of Double Flapper-nozzle Servo Valve

The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo valve would reach 120℃ and the valve core and valve sleeve deform in a short amount of time. So the control precision of servo valve significantly decreases and the clamping stagnation phenomenon of valve core appears. In order to solve the problem of degraded control accuracy and clamping stagnation of servo valve under large temperature difference circumstance, the numerical simulation of heat-fluid-solid coupling by using finite element method is done. The simulation result shows that zero position leakage of servo valve is basically impacted by oil temperature and change of fit clearance. The clamping stagnation is caused by warpage-deformation and fit clearance reduction of the valve core and valve sleeve. The distribution roles of the temperature and thermal-deformation of shell, valve core and valve sleeve and the pressure, velocity and temperature field of flow channel are also analyzed. Zero position leakage and electromagnet＇s current when valve core moves in full-stroke are tested using Electro-hydraulic Servo-valve Characteristic Test-bed of an aerospace sciences and technology corporation. The experimental results show that the change law of experimental current at different oil temperatures is roughly identical to simulation current. The current curve of the electromagnet is smooth when oil temperature is below 80℃, but the amplitude of current significantly increases and the hairy appears when oil temperature is above 80℃. The current becomes smooth again after the warped valve core and valve sleeve are reground. It indicates that clamping stagnation is caused by warpage-deformation and fit clearance reduction of valve core and valve sleeve. This paper simulates and tests the heat-fluid-solid coupling of double flapper-nozzle servo valve, and the obtained results provide the reference value for the design of double flapper-nozzle force feedback servo valve.

Oxidative coupling of methane in a fixed bed reactor over perovskite catalyst:A simulation study using experimental kinetic model

The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separate energy equations for the gas and solid phases coupled with an experimental kinetic model. A lumped kinetic model containing four main species CH4, O2, COx （CO2, CO）, and C2 （C2H4 and C2H6） was used with a plug flow reactor model as well. The results from the model agreed with the experimental data. The model was used to analyze the influence of temperature and feed gas composition on the conversion and selectivity of the reactor performance. The analytical results indicate that the conversion decreases, whereas, C2 selectivity increases by increasing gas hourly space velocity （GHSV） and the methane conversion also decreases by increasing the methane to oxygen ratio.

NUMERICAL SIMULATION OF THE THERMO-MECHANICAL PROCESS FOR BEAM BLANK CONTINUOUS CASTING

The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.