Liquid-solid Coupling Research Articles

Seismic analysis of seabeds with irregular terrain under 2D oblique incident waves by finite/infinite element method

Seismic response of seabeds with irregular terrain (i.e., canyons and hills) under oblique incident P- or SV- waves is analyzed by the finite/infinite element method (FIEM). Initially, based on the exact solution for a corresponding free field, the equivalent seismic forces for the solid-liquid site are calculated. Then, the solid-liquid coupling for the terrain slope boundaries is implemented with key numerical parameters selected for the FIEM. Furthermore, the adequacy of the solid and liquid meshes each was verified by comparing the solutions with existing ones for two extreme cases. Finally, the effect of irregular terrain on the seabed seismic responses was assessed for P- and SV-waves based on the 1999 Chi-Chi Earthquake. Key conclusions include: (1) The difference in the vertical acceleration aymax for SV-waves is the largest for the irregular seabed with one canyon and one hill under super-critical conditions. (2) The canyon shows more prominent amplification for the seabed responses near the seismic source for both P- and·SV-waves.·(3) Deeper liquid drastically amplify the difference in horizontal accelerations on two sides of irregular local terrain for sub-critical P-waves and near-critical SV-waves· Additional observations made for the numerical analysis are given inside the text.

Seismic Response of Foundation Settlement for Liquid Storage Structure in Collapsible Loess Areas

To investigate the impact of foundation settlement in loess areas on the dynamic response of liquid storage structure (LSS) under seismic motion, a finite element analysis model of the liquid–solid coupling of LSS was established using ADINA V9.6 software. By analyzing the dynamic response patterns of LSS under seismic motion with foundation failure, this study examines the effects of foundation failure and the direction of seismic wave incidence on the equivalent stress, maximum shear stress, wall displacement, and liquid sloshing wave height of the structure. The results indicate that among the three foundation failure scenarios, foundation failure at the center of the tank bottom has the least impact on the structural dynamic response. In contrast, foundation failure affecting one-fourth of the tank base has the greatest impact. Furthermore, compared to seismic motion along the X-axis, the dynamic response of the structure is more significantly affected when seismic motion co-occurs along the X-Z-axis.

Overview of Theory, Simulation, and Experiment of the Water Exit Problem

The water exit problem, which is ubiquitous in ocean engineering, is a significant research topics in the interaction between navigators and water. The study of the water exit problem can help to improve the structural design of marine ships and underwater weapons, allowing for better strength and movement status. However, the water exit problem involves complex processes such as three-phase gas–liquid–solid coupling, cavitation, water separation, liquid surface deformation, and fragmentation, making it challenging to study. Following work carried out by many researchers on this issue, we summarize recent developments from three aspects: theoretical research, numerical simulation, and experimental results. In theoretical research, the improved von Karman model and linearized water exit model are introduced. Several classical experimental devices, data acquisition means, and cavitation approaches are introduced in the context of experimental development. Three numerical simulation methods, namely, the BEM (Boundary Element Method), VOF (Volume of Fluid), and FVM (Finite Volume Method) with LES (Large Eddy Simulation) are presented, and the respective limitations and shortcomings of these three aspects are analyzed. Finally, an outlook on future research improvements and developments of the water exit problem is provided.

Dynamic responses of non-isolated and isolated liquid storage tanks under mainshock-aftershock sequences

After the mainshock earthquake, aftershocks often occur frequently, and the destructive effect of aftershocks on structures cannot be ignored. Based on the potential flow theory, the mechanical behavior of liquid is simulated. Considering liquid-solid coupling, material nonlinearity, and liquid sloshing, a three-dimensional numerical calculation model of liquid storage tank is established. No pulse-like and pulse-like mainshock-aftershock sequences are selected, and the seismic responses of non-isolated and isolated liquid storage tanks is compared and studied. The results show that the mainshock-aftershock sequence has a significant impact on the dynamic responses of the liquid storage tank, and the maximum increase rate reaches to 195.4 %. The effect of the mainshock-aftershock sequence on the dynamic responses of non-isolated liquid storage tanks is greater than on the dynamic response of isolated liquid storage tanks. Compared with the pulse-like mainshock-aftershock sequences, the no pulse-like mainshock-aftershock sequences generally have a greater impact on the dynamic responses of the liquid storage tanks. Overall, isolation significantly reduce the impact of mainshock-aftershock sequences on the dynamic responses of liquid storage tanks, and effectively enhance the seismic performance of liquid storage tanks under mainshock-aftershock sequences.

Analysis on transient wave propagation in the soft matter of hydrogels

The precise control and structural design of the soft matter of the hydrogel subjected to dynamic loading require an in-depth understanding of its transient wave characteristics. However, for the complex solid–liquid coupling effects of the soft matter of hydrogels, the traditional single-phase wave model is no longer applicable. In this study, a transient wave propagation model that considers the solid–liquid coupling effects and its computational method is proposed. The wave-governing equations of hydrogels are derived based on the mass conservation and dynamic equilibrium equations, where the motions of solid and liquid phases are independently described. After transforming the governing equations into the equivalent weak form, numerical solutions for transient wave propagation in one- and two-dimensional hydrogels are calculated. The accuracy of the proposed model is verified by a comparison between semi-analytical and commercial finite element method solutions. We observe that the travelling and reflection processes of the two types of compression waves, P1 and P2, with different wave speeds are captured when the dynamic coefficient of permeability kf increases. Furthermore, the influence of solid–liquid coupling effects on the transient responses of hydrogels is discussed. The results show that the hydrogels exhibit the dynamic characteristics of a single-phase medium to a certain extent when kf is sufficiently small.

Experimental investigation of dynamic characteristics of leaching tubing for solution mining of salt cavern carbon and energy storage

Salt caverns are extensively utilized for storing various substances such as fossil energy, hydrogen, compressed air, nuclear waste, and industrial solid waste. In China, when the salt cavern is leached through single-well water solution mining with oil as a cushion, engineering challenges arise with the leaching tubing, leading to issues like damage and instability. These problems significantly hinder the progress of cavern construction and the control of cavern shape. The primary cause of this is the flow-induced vibration instability of leaching tubing within a confined space, which results in severe bending or damage to the tubing. This study presents a model experimental investigation on the dynamic characteristics of leaching tubing using a self-developed liquid-solid coupling physical model experiment apparatus. The experiment utilizes a silicone-rubber pipe (SRP) and a polycarbonate pipe (PCP) to examine the effects of various factors on the dynamic stability of cantilevered pipes conveying fluid. These factors include external space constraint, flexural rigidity, medium outside the pipe, overhanging length, and end conditions. The experiments reveal four dynamic response phenomena: water hammer, static buckling, chaotic motion, and flutter instability. The study further demonstrates that the length of the external space constraint has a direct impact on the flutter critical flow velocity of the cantilevered pipe conveying fluid. Additionally, the flutter critical flow velocity is influenced by the end conditions and different external media.

Design Method of the Stroke Ring Based on Deformation Pre-Compensation

A common consensus is that an optimized curve profile of the stroke ring in the multiple-stroke piston motor can make the output torque more stable. However, the ring generates elastic deformation during operation, which causes the piston component movement trajectory to deviate from the ideal design curve. To address this issue, first, a liquid–solid coupling simulation model was established to obtain the deformation of the ring, and the accuracy of the model was verified through experiments. Second, a stroke ring curve design method based on elastic deformation pre-compensation was proposed. Through this method, a compensated curve can be obtained to make the actual working curve more in line with the ideal curve. Finally, the dynamic characteristics of three different types of multiple-stroke piston motor curves were analyzed—the ideal design curve, the uncompensated working curve, and the compensated working curve. The results showed that the motor torque pulsation rates are 0.821%, 4.723%, and 0.986%, respectively, and the compensated working curve has a relatively reduced pulsation rate of 79.12% compared to the uncompensated working curve, which verifies that this design method can effectively improve motor performance.

Drag on nanoparticles in a liquid: from slip to stick boundary conditions.

Stokes' law with stick boundary conditions has been widely accepted for the transport of microscale particles in a liquid. For nanoparticles, however, the hydrodynamic boundary conditions become unclear. In this work, the drag force acting on nanoparticles suspended in a liquid and the hydrodynamic boundary coefficient were calculated by using molecular dynamics simulations. For weak interfacial couplings, slip boundary conditions can be used to describe the particle transport, whereas at strong interfacial couplings, the hydrodynamic boundary coefficient converges to a value greater than the prediction by Stokes' law. In the present paper, we propose a density accumulation length to determine the effective particle size, which makes Stokes' law valid for nanoparticles. For a copper nanoparticle suspended in an argon liquid, the density accumulation length increases to 0.32 nm with increasing solid-liquid coupling strength. Furthermore, it is found that there exists a transition from slip to stick boundary conditions as the solid-liquid intermolecular coupling strength increases. The results presented in this work provide guidance for the prediction and manipulation of the transport properties of nanoparticles in a liquid.

Tribological behaviors of an attapulgite-graphene nanocomposite as an additive for mineral lubricating oil.

In this work, an attapulgite-graphene nanocomposite was prepared. The tribological properties of the prepared attapulgite-graphene nanocomposite as an additive for 200SN mineral lubricating oil were investigated using an SRV-IV tribometer through ball-on-disk contact mode for the first time. The characterization of the prepared nanocomposite indicated that attapulgite nanofibers are enveloped by the graphene nanosheets and present fine combination. The tribological test results show that the friction-reducing and antiwear properties of 200SN were obviously improved by adding the attapulgite-graphene nanocomposite. Through the characterization and analysis of the worn surface and cross-section, it was found that a tribofilm composed of Fe, Fe3O4, FeO, Fe2O3, FeOOH, graphite, graphene, SiO2 and organic compounds was formed on the worn surface. Furthermore, the bonding between the tribofilm and steel matrix is tight. The tribofilm and lubricating oil achieve a solid-liquid coupling lubrication effect, which is responsible for the improvement of the friction-reducing and antiwear properties.

LIQUID-SOLID COUPLING RESPONSE OF SURROUNDING ROCK MASS OF LARGE-DIAMETER RIVER-CROSSING SHIELD TUNNEL

The purpose is to investigate the response of seepage field, displacement field and stress field in the surrounding rock mass during dynamic tunneling in soft soil area. Relied on a large-diameter river-crossing shield tunnel project, considering driving force, shield tail grouting pressure, and the friction resistance between the shield shell and the soil, a three-dimensional fine tunnel model considering the liquid-solid coupling effect in the soil during dynamic tunneling was established by employing the finite difference method. The response characteristics of pore water pressure, displacement and stress in the surrounding rock mass were obtained. The results show that during shield tunneling and shield tail grouting, the pore water pressure in the range of 0.5 times the hole diameter around the tunnel decreases and increases respectively due to the liquid-solid coupling in the surrounding rock mass. When the shield tunneling moves away, the pore water pressure of the soil near the vault decreases, and the pore water pressure near the tunnel arch bottom increases. The impact range of shield tail grouting on the vertical settlement of the upper soil is about 0.5 times the hole diameter. The shield tail grouting can effectively reduce the vertical settlement of the top soil and slow down the vertical uplift of the bottom soil. During shield tunneling the vertical stress distribution of the soil above the vault of the working position and around the excavation surface is funnel-shaped, and the vertical stress around the excavated tunnel decreases.

Research on the Mechanism and Evolution Law of Delayed Water Inrush Caused by Fault Activation with Mining

Confined water inrush caused by fault activation is the main form of water disaster in deep mining. With theoretical analysis and similar simulation tests, the mechanism and evolution law of delayed water inrush caused by fault activation are revealed. At the theoretical level, the expansion and extension of the internal microstructure in the fault zone under the action of the mining stress field and seepage field are the essential causes of fault activation. Overlying strata movement and surrounding rock creep failure are the basic reasons for delayed water inrush caused by fault activation, and delayed time caused by surrounding rock creep failure is much longer than that of overlying strata movement. A similar simulation test was carried out with self-development solid–liquid coupling with similar simulation materials; the results show that delayed water inrush caused by fault activation with mining includes three stages. Micro-activation stage: Water inrush weakness point is formed because of the expansion and extension of the micro-fissure and structure at the bottom of the fault zone. Macro-activation stage: With the change in the stress of the waterproof coal pillar and surrounding rock, the micro-fissures and structures in the stress relief area and tension area of the fault zone expand and extend sharply; meanwhile, water intrudes into the interlayer stratification of the floor in the stress relief area, forming a strong laminar flow phenomenon, and cracks in the floor form and expand; finally, water-conducting channels in the fault zone and floor are formed. Water inrush stage: The waterproof coal pillar and water-resisting layer fail and are destroyed, and the first confined water inrush point is located at the junction of the waterproof coal pillar and gob floor.

Failure probability of a liquid storage tank with three-dimensional base-isolation under earthquake action

Liquid storage tanks are widely used in the petrochemical industry, but the fixed-base liquid storage tanks have a high probability of earthquake damage, therefore, the research on isolation and damping control of liquid storage tank has attracted wide attention. Based on the dual control requirements of the seismic response of the liquid storage tank itself and the liquid sloshing wave height, a new type three-dimensional isolation bearing suitable for the liquid storage tank has been obtained. Taking into account the material nonlinearity, the liquid-solid coupling, and the liquid sloshing characteristics, a numerical calculation model of the new type three-dimensional isolated liquid storage tanks is established by ADINA. The seismic vulnerability of the fixed-base and the new type three-dimensional isolated liquid storage tanks are investigated, the failure probability of liquid storage tank under different limit states is discussed, and the annual average failure probability and design life failure probability of the liquid storage tank are finally analyzed. The results show that, the probability of fully damaged wave height and axial stress of fixed-base liquid storage tank is 72.0% and 49.6% respectively, under the 9 degree rare earthquake. Compared to the fixed-base liquid storage tank, the failure probability of the axial compressive stress, hoop stress, and liquid sloshing wave height of the isolated liquid storage tank is significantly reduced, but the displacement of the isolation layer becomes the most easily damaged dynamic response. Besides, the new type three-dimensional isolation can significantly reduce the annual average failure probability and the design life failure probability of the liquid storage tank.

Microchannels Formed Using Metal Microdroplets.

The metal microdroplet deposition manufacturing technique has gained extensive attention due to its potential applications in microstructure fabrication. In order to fabricate components such as microchannel heat sinks and microchannel reactors, this paper investigates the interactions and influences between microdroplets and substrates, as well as between microdroplets themselves. The transient phenomena during the fusion of metal microdroplets in contact with the substrate and the formation of inclined columns, as well as the solid-liquid coupling and morphology formation processes during the collision between microdroplets, are analyzed. The influence of microdroplet spacing on the morphology of microchannels during their formation is specifically studied. A three-dimensional finite element numerical model for the deposition of metal microdroplets forming inclined pillars is established based on the volume of fluid (VOF) method. The model treats the protective gas around the microdroplet as an empty zone and the microdroplet as a single-phase fluid. Simulation analysis is conducted to investigate the forming patterns of unsupported microdroplets at different spacing and their impact on the fusion morphology of microchannel components. Building upon this, a series of validation experiments are conducted using a piezoelectric microdroplet generator to produce uniform aluminum alloy microdroplets with a diameter of approximately 600 μm. A method for fabricating metal microchannel structures is obtained, which is expected to be applied in fields such as scattering structures for high-power electronic devices and microreactors in microchemical fields.

Simulated Prediction of Roof Water Breakout for High-Intensity Mining under Reservoirs in Mining Areas in Western China

China is rich in coal resources under water bodies. However, the safety prediction of high-intensity mining under water bodies has long been one of the problems encountered by the coal industry. It is of great significance to realize safe mining under water bodies, improve the recovery rate of coal resources and protect reservoir resources. Therefore, this article takes the No. 5 coal seam and No. 11 mining area of the Wangwa Coal Mine as the research object, and integrates physical simulation, numerical simulation, theoretical analysis, and other methods to study the development height of water-conducting fracture zones in fully mechanized top coal caving mining. Solid–liquid coupling physical simulation tests reveal the failure characteristics of overlying strata in goaf and the seepage law of reservoir water under the influence of mining. By comparing the monitoring data of borehole leakage, the measured data obtained by borehole monitoring with the height data of the water-conducting fracture zone obtained by the traditional empirical formula of three-under standard, the error between the two is as high as −29.39%. In this case, the variance correction coefficient is used to correct the empirical formula, and on this basis, in order to effectively protect the surface water dam and water body, the mining height of the coal seam in the working face with limited height mining is inversely derived. The research results provide a basis for the safety prediction of high-intensity mining under the reservoir dam in the ecologically fragile areas of western China and a scientific guarantee for the formulation of safety measures under such conditions.

Nanoscale Investigation of Bubble Nucleation and Boiling on Random Rough Surfaces.

Surface roughness is one of the significant factors affecting liquid-vapor phase change heat transfer. This paper explores the effect of surface roughness on bubble nucleation and boiling heat transfer, as well as the microscopic mechanism, by constructing random rough surfaces using molecular dynamics (MD) simulation. Bubbles randomly nucleate on a flat surface and tend to nucleate in pits on rough surfaces. The pits on the random rough surface gather more argon atoms than the protrusions, forming low potential energy regions on the surface, thus providing stable nucleation sites for bubbles. As the surface roughness increases, bubble generation, merging, and growth are advanced. In addition, rough surfaces offer a larger effective heat transfer area for the heat transfer process, increase the strength of solid-liquid coupling, and obtain smaller solid-liquid interaction energy. The critical heat flux (CHF) value positively correlates with surface roughness. As the roughness increases, the surface superheat at the onset of CHF decreases accordingly. This paper provides new insights into the mechanism of heat transfer enhancement on rough surfaces and surface design in thermal management.

Analysis of Water Inrush Disaster Mechanism of Inter-Layer Rocks between Close Coal Seams under the Influence of Mining

With the gradual increase in the mining depth of coal resources, the destruction of the rock structure of the inter-layered rock of the near coal seam under the influence of mining has led to the frequent occurrence of water-inrush disasters in mines, which seriously affects the safety of mine production and the safety of the people in the underground. Therefore, it is important to study the mechanism of the water inrush of the rock between the coal seams under the influence of mining to control the occurrence of water inrush disasters and protect the loss of groundwater resources. This paper takes the Hanjiawan coal mine with typical stratigraphic characteristics as the background for research and studies the structural characteristics of interlayer rock breakage and the solid–liquid coupling inrush water disaster mechanism during the mining of 2−2 and 3−1 coals. The study shows that according to the damage degree and destruction depth of the inter-layered rock caused by the mining of the upper and lower coal seams, combined with the slip line theory and the “three bands” collapse theory, the inter-layered rock is classified into a completely fractured inter-layer, a fractured–broken stacked inter-layer, and a fractured–broken–fractured combined inter-layered rock using L≤hm+Hk2′, L>hm+Hk2′, and L≥hm+Hli2′ as the discriminating criteria. Combined with the structural classification of inter-layer rock and the discriminating criteria, we used similar simulation experiments and on-site research to analyze the evolution law and distribution characteristics of four types of inter-layer rock water-inrush fractures in different mines and put forward the classification of inter-layer rock water-inrush channels based on the width, length, and penetration of the fractures. Based on the characteristics of the water-inrush channel of inter-layer rock, we constructed the network-boundary inrush water calculation model of inter-layered rock and network-attach-boundary inrush water calculation model, solved the water movement of the water-inrush channel in the model by transforming the flat flow state, fracture to flow state, and pore-fracture flow state, and finally revealed the mechanism of the disaster by which water-inrush of inter-layer rocked was induced. Finally, we revealed its mechanism of inducing the inter-layer rock inrush water disaster. Our research enriches the theory and research ideas of the water-inrush disaster, provides theoretical support and a basis for the control of water-inrush disasters in similar conditions, and ensures the safe production of mines.

Research on oil boom performance based on Smoothed Particle Hydrodynamics method.

To address the issues of fluid-solid coupling, instability in the liquid two-phase flow, poor computational efficiency, treating the free surface as a slip wall, and neglecting the movement of oil booms in simulating oil spill containment, this study adopts the Smoothed Particle Hydrodynamics (SPH) method to establish a numerical model for solid-liquid coupling and liquid two-phase flow, specifically designed for oil boom containment and control. The DualSPHysics solver is employed for numerical simulations, incorporating optimized SPH techniques and eight different skirt configurations of the oil boom into the numerical model of two-phase liquid interaction. By setting relevant parameters in the SPH code to enhance computational efficiency, the variations in centroid, undulation, and stability of undulation velocity for different oil boom shapes are observed. The experimental results demonstrate that the improved oil boom exhibits superior oil containment performance. These findings provide a theoretical basis for the design of oil boom skirt structures.

Implications of liquid impurities filled in breaking cracks on nonlinear acoustic modulation response: Mechanisms, phenomena and potential applications

Supervised strategies for ultrasonic-based cracks detection in existing structures has been one of particular concerns in the engineering community. However, current researches mostly focus on the clean microcracks, while ignoring the fact that such cracks will be contaminated with liquid impurities (lubricants, silicone oil, etc.) due to its exposure to the operational environment. Actually, the presence of liquid may lead to a non-negligible change in the acoustic performance and interfere with damage diagnosis. The central contribution of this work is to reveal the nonlinear acoustic mechanisms and phenomena induced by the gas–liquid–solid coupling interfaces in such liquid-filled cracks and investigate their potential applications for such crack detection. Firstly, the state equation for a liquid-filled cavity in beams is established to determine the nonlinear dependences of stress–strain relationship at crack wall on the viscous and capillary pressure in liquid. Secondly, a modified acoustic modulation method, in which swept signal is used as the pumping wave to sufficiently trigger and enhance the nonlinear behaviours at the gas–liquid–solid interfaces without requiring exact reference information, is utilized for detecting sub-millimeter cracks filled with different liquid. The results suggest that the liquid fillers may cause considerable increase of the damage-related acoustic nonlinearities involving hysteresis nonlinearity, nonlinear elastic behavior and dissipative nonlinearity. Additional analysis provides compelling evidence for the following valuable conclusions: 1) the capillary pressure is a key factor that promotes conducive to the nonlinear elastic behavior that is the primary source of the first-order sum-difference frequency; 2) the viscous pressure contributes to hysteresis as well as the nonlinear dissipation, and the latter promotes the second-order sum frequency. Further, such nonlinear performance especially dissipation offer promising prospects for nonlinear acoustic diagnostic of such cracks. All results obtained in this study will aid in understanding the nonlinear wave processes and further improving the reliability of crack detection in engineering.

Solid-liquid coupling of in-phase and out-of-phase vibrations for two piezoelectric disks

Theoretical analysis is used to solve the solid–liquid coupled vibration characteristics for two piezoelectric disks in a limited space and verified with experimental measurements. Based on the Navier–Stoke equation and equilibrium of motion, the Rayleigh-Ritz and energy methods were developed to analyze the solid–liquid coupled system. Following an entire systematic procedure, two piezoelectric disks immersed in hydrostatic compressible fluids were analyzed using the piezoelectric constitutive equation to perform the in-phase and out-of-phase vibration. The full-field, non-contact, optical technique, called amplitude-fluctuation electronic speckle pattern interferometry (AF-ESPI), was applied to experimentally measure the vibration properties in a solid–liquid coupling system, including the mode shapes and resonant frequencies. The experimental results are presented in good agreement with the theoretical solutions and finite element calculations. Under the solid–liquid coupled system, the distributions of pressure fields and velocity vectors of fluids were also determined in the correspondent vibrating modes. Finally, it is interesting to show the tendency to close the same resonant frequency both on the in-phase and out-of-phase vibrations when the fluid depth increases. The research outcomes can provide a guide in the design of hydrostatic components.

Modeling of solid–liquid coupling and material removal in robotic wet polishing

Modeling of solid–liquid coupling and material removal in robotic wet polishing