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.

A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Coupled multiphysical model for investigation of influence factors in the application of microbially induced calcite precipitation

The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.

An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

Nonlinearly induced entanglement in dissipatively coupled optomechanical system

Nonlinearly induced steady-state photon–phonon entanglement of a dissipative coupled system is studied in the bistable regime. Quantum dynamical characteristics are analysed by solving the mean-field and fluctuation equations of the system. It is shown that dissipative coupling can induce bistable behaviour for the effective dissipation of the system.Under suitable parameters, one of the steady states significantly reduces the dissipative effect of the system. Consequently,a larger steady-state entanglement can be achieved compared to linear dynamics. Furthermore, the experimental feasibility of the parameters is analysed. Our results provide a new perspective for the implementation of steady-state optomechanical entanglement.

Seismic performance evaluation of hybrid coupled shear wall system with shear and flexural fuse-type steel coupling beams

Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.

Effects of asymmetric coupling and boundary on the dynamic behaviors of a random nearest neighbor coupled system

This study investigates the dynamical behaviors of nearest neighbor asymmetric coupled systems in a confined space.First, the study derivative analytical stability and synchronization conditions for the asymmetrically coupled system in an unconfined space, which are then validated through numerical simulations. Simulation results show that asymmetric coupling has a significant impact on synchronization conditions. Moreover, it is observed that irrespective of whether the system is confined, an increase in coupling asymmetry leads to a hastened synchronization pace. Additionally, the study examines the effects of boundaries on the system's collective behaviors via numerical experiments. The presence of boundaries ensures the system's stability and synchronization, and reducing these boundaries can expedite the synchronization process and amplify its effects. Finally, the study reveals that the system's output amplitude exhibits stochastic resonance as the confined boundary size increases.

Customized modulation on plasma uniformity by non-uniform magnetic field in capacitively coupled plasma

A two-dimensional fluid model based on COMSOL Multiphysics is developed to investigate the modulation of static magnetic field on plasma homogeneity in a capacitively coupled plasma(CCP)chamber. To generate a static magnetic field, direct current is applied to a circular coil located at the top of the chamber. By adjusting the magnetic field's configuration, which is done by altering the coil current and position, both the plasma uniformity and density can be significantly modulated. In the absence of the magnetic field, the plasma density exhibits an inhomogeneous distribution characterized by higher values at the plasma edge and lower values at the center. The introduction of a magnetic field generated by coils results in a significant increase in electron density near the coils. Furthermore, an increase in the sets of coils improves the uniformity of the plasma. By flexibly adjusting the positions of the coils and the applied current,a substantial enhancement in overall uniformity can be achieved. These findings demonstrate the feasibility of using this method for achieving uniform plasma densities in industrial applications.

A Coupled Thermomechanical Crack Propagation Behavior of Brittle Materials by Peridynamic Differential Operator

This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic differential operator(PDDO),eliminating the need for calibration procedures.The model employs a multi-rate explicit time integration scheme to handle varying time scales in multi-physics systems.Through simulations conducted on granite and ceramic materials,this model demonstrates its effectiveness.It successfully simulates thermal damage behavior in granite arising from incompatible mineral expansion and accurately calculates thermal crack propagation in ceramic slabs during quenching.To account for material heterogeneity,the model utilizes the Shuffle algorithm andWeibull distribution,yielding results that align with numerical simulations and experimental observations.This coupled thermomechanical model shows great promise for analyzing intricate thermomechanical phenomena in brittle materials.

Coexisting fast–slow dendritic traveling waves in a 3D-array electric field coupled neuronal network

Coexistence of fast and slow traveling waves without synaptic transmission has been found in hhhippocampal tissues,which is closely related to both normal brain activity and abnormal neural activity such as epileptic discharge. However, the propagation mechanism behind this coexistence phenomenon remains unclear. In this paper, a three-dimensional electric field coupled hippocampal neural network is established to investigate generation of coexisting spontaneous fast and slow traveling waves. This model captures two types of dendritic traveling waves propagating in both transverse and longitude directions: the N-methyl-D-aspartate(NMDA)-dependent wave with a speed of about 0.1 m/s and the Ca-dependent wave with a speed of about 0.009 m/s. These traveling waves are synaptic-independent and could be conducted only by the electric fields generated by neighboring neurons, which are basically consistent with the in vitro data measured experiments. It is also found that the slow Ca wave could trigger generation of fast NMDA waves in the propagation path of slow waves whereas fast NMDA waves cannot affect the propagation of slow Ca waves. These results suggest that dendritic Ca waves could acted as the source of the coexistence fast and slow waves. Furthermore, we also confirm the impact of cellular spacing heterogeneity on the onset of coexisting fast and slow waves. The local region with decreasing distances among neighbor neurons is more liable to promote the onset of spontaneous slow waves which, as sources, excite propagation of fast waves. These modeling studies provide possible biophysical mechanisms underlying the neural dynamics of spontaneous traveling waves in brain tissues.

Multi-soliton solutions of coupled Lakshmanan-Porsezian-Daniel equations with variable coefficients under nonzero boundary conditions

This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.

A mechanically coupled MEMS filter with high-Q width extensional mode resonators

This work presents a novel radio frequency(RF)narrowband Si micro-electro-mechanical systems(MEMS)filter based on capacitively transduced slotted width extensional mode(WEM)resonators.The flexibility of the plate leads to multiple modes near the target frequency.The high Q-factor resonators of around 100000 enable narrow bandwidth filters with small size and simplified design.The 1-wavelength and 2-wavelength WEMs were first developed as a pair of coupled modes to form a passband.To reduce bandwidth,two plates are coupled with aλ-length coupling beam.The 79.69 MHz coupled plate filter(CPF)achieved a narrow bandwidth of 8.8 kHz,corresponding to a tiny 0.011%.The CPF exhibits an impressive 34.84 dB stopband rejection and 7.82 dB insertion loss with near-zero passband ripple.In summary,the RF MEMS filter presented in this work shows promising potential for application in RF transceiver front-ends.

Coupled effects of climate change and human activities on vegetation dynamics in the Southwestern Alpine Canyon Region of China

The driving effects of climate change and human activities on vegetation change have always been a focal point of research.However,the coupling mechanisms of these driving factors across different temporal and spatial scales remain controversial.The Southwestern Alpine Canyon Region of China(SACR),as an ecologically fragile area,is highly sensitive to the impacts of climate change and human activities.This study constructed a vegetation cover dataset for the SACR based on the Enhanced Vegetation Index(EVI)from 2000 to 2020.Spatial autocorrelation,Theil-Sen trend,and Mann-Kendall tests were used to analyze the spatiotemporal characteristics of vegetation cover changes.The main drivers of spatial heterogeneity in vegetation cover were identified using the optimal parameter geographic detector,and an improved residual analysis model was employed to quantify the relative contributions of climate change and human activities to interannual vegetation cover changes.The main findings are as follows:Spatially,vegetation cover exceeds 60%in most areas,especially in the southern part of the study area.However,the border area between Linzhi and Changdu exhibits lower vegetation cover.Climate factors are the primary drivers of spatial heterogeneity in vegetation cover,with temperature having the most significant influence,as indicated by its q-value,which far exceeds that of other factors.Additionally,the interaction q-value between the two factors significantly increases,showing a relationship of bivariate enhancement and nonlinear enhancement.In terms of temporal changes,vegetation cover shows an overall improving trend from 2000 to 2020,with significant increases observed in 68.93%of the study area.Among these,human activities are the main factors driving vegetation cover change,with a relative contribution rate of 41.31%,while climate change and residual factors contribute 35.66%and 23.53%,respectively.By thoroughly exploring the coupled mechanisms of vegetation change,this study provides important references for the sustainable management and conservation of the vegetation ecosystem in the SACR.

Abnormal transition of the electron energy distribution with excitation of the second harmonic in low-pressure radio-frequency capacitively coupled plasmas

The self-excited second harmonic in radio-frequency capacitively coupled plasma was significantly enhanced by adjusting the external variable capacitor.At a lower pressure of 3 Pa,the excitation of the second harmonic caused an abnormal transition of the electron energy probability function,resulting in abrupt changes in the electron density and temperature.Such changes in the electron energy probability function as well as the electron density and temperature were not observed at the higher pressure of 16 Pa under similar harmonic changes.The phenomena are related to the influence of the second harmonic on stochastic heating,which is determined by both amplitude and the relative phase of the harmonics.The results suggest that the self-excited high-order harmonics must be considered in practical applications of lowpressure radio-frequency capacitively coupled plasmas.

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.

Effects of external parameters on plasma characteristics and uniformity in a dual cylindrical inductively coupled plasma

The dual cylindrical inductively coupled plasma source,compared to the conventional structure of inductively coupled plasma source,can significantly improve the uniformity of plasma.It has an enhanced potential for application in processes,such as etching and ashing.A uniform plasma can be obtained by allowing the remote plasma from the upper chamber modulate the main plasma generated in the lower chamber.In this study,a fluid model was employed to investigate a dual cylindrical inductively coupled Ar/O_(2)discharge.The effects of external parameters on electron density,electron temperature,O atomic density,and plasma uniformity in the main chamber were studied,and the reasons were analyzed.The results of this study show that remote power can control the plasma uniformity and increase the plasma density in the main chamber.As the remote power increased,plasma uniformity improved initially and then deteriorated.The main power affected the plasma density at the edge of the main chamber and can modulate the plasma density in the main chamber.The gas pressure affected both the uniformity and density of the plasma.As the gas pressure increased,the plasma uniformity deteriorated,but the free radical density improved.

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

Comparative study of aeroacoustic performance of 1/8 and 1/1 pantographs coupled with cavity

The technology of pantograph sinking in the cavity is generally adopted in the new generation of high-speed trains in China for aerodynamic noise reduction in this region. This study takes a high-speed train with a 4-car formation and a pantograph as the research object and compares the aerodynamic acoustic performance of two scale models, 1/8 and 1/1, using large eddy simulation and Ffowcs Williams–Hawkings integral equation. It is found that there is no direct scale similarity between their aeroacoustic performance. The 1/1 model airflow is separated at the leading edge of the panhead and reattached to the panhead, and its vortex shedding Strouhal number(St) is 0.17. However, the 1/8 model airflow is separated directly at the leading edge of the panhead, and its St is 0.13. The cavity's vortex shedding frequency is in agreement with that calculated by the Rooster empirical formula. The two scale models exhibit some similar characteristics in distribution of sound source energy, but the energy distribution of the 1/8 model is more concentrated in the middle and lower regions. The contribution rates of their middle and lower regions to the radiated noise in the two models are 27.3% and 87.2%, respectively. The peak frequencies of the radiated noise from the 1/1 model are 307 and 571 Hz. The 307 Hz is consistent with the frequency of panhead vortex shedding, and the 571 Hz is more likely to be the result of the superposition of various components. In contrast, the peak frequencies of the radiated noise from the 1/8 scale model are 280 and 1970 Hz. The 280 Hz comes from the shear layer oscillation between the cavity and the bottom frame, and the 1970 Hz is close to the frequency at which the panhead vortex sheds. This shows that the scaled model results need to be corrected before applying to the full-scale model.

Efficient Inverse Analysis for Solving a Coupled Conduction,Convection and Radiation Problem Involving Non⁃gray Participating Media

The presence of non-gray radiative properties in a reheating furnace’s medium that absorbs,emits,and involves non-gray creates more complex radiative heat transfer problems.Furthermore,it adds difficulty to solving the coupled conduction,convection,and radiation problem,leading to suboptimal efficiency that fails to meet real-time control demands.To overcome this difficulty,comparable gray radiative properties of non-gray media are proposed and estimated by solving an inverse problem.However,the required iteration numbers by using a least-squares method are too many and resulted in a very low inverse efficiency.It is necessary to present an efficient method for the equivalence.The Levenberg-Marquardt algorithm is utilized to solve the inverse problem of coupled heat transfer,and the gray-equivalent radiative characteristics are successfully recovered.It is our intention that the issue of low inverse efficiency,which has been observed when the least-squares method is employed,will be resolved.To enhance the performance of the Levenberg-Marquardt algorithm,a modification is implemented for determining the damping factor.Detailed investigations are also conducted to evaluate its accuracy,stability of convergence,efficiency,and robustness of the algorithm.Subsequently,a comparison is made between the results achieved using each method.